Size: 677
Comment:
|
Size: 2990
Comment:
|
Deletions are marked like this. | Additions are marked like this. |
Line 1: | Line 1: |
= Overview = Archiving SPC-focus products. There are two major categories (global and regional) with two subdivisions (topography and backplanes -- aka other data). The regional products are derived from bigmaps. === Tools === These tools are stored in GitHub as [archivingTools|https://github.com/StereoPhotoClinometry/archivingTools]. * [[convertLatLon]] - Reads an X,Y,Z (or X,Y,Z,V) list and converts it into lat, lon, distance or lat, lon, value * [[backplanesGMT]] - Reads a bigmap and generates a list of latlon, radius, albedo, slope. Use paste latlon radius > tmp for GMT. * [[backplanesISIS]] - Reads a bigmap and generates 2D matrix for radius, albedo, slope, lat, lon * [[phasei]] - Reads a bigmap and image (the SUMFILE) and generates 2D matrix for incidence, emission, phase angles * coverage_p_low -- Reads all images and maplets. Creates a pgm that has every pixel marked for every image that it has, counter is only 1. * [[map_coverage_p_low]] - Counts the number of images that cover each bigmap pixel. Creates MAPNM-cov.txt and MAPNM-cov.pgm -------- |
|
Line 25: | Line 44: |
Line 30: | Line 51: |
Global ICQ resolution to bin sizes [degrees] || ICQ Q Value || Vertices || Bin Size [deg] || || 512 || 1,579,014 || 0.041 || || 256 || 396,294 || 0.164 || || 128 || 99,846 || 0.649 || || 64 || 25,350 || 2.556 || || 32 || 6,534 || 9.917 || |
|
Line 31: | Line 59: |
{{{ convertLonLat > radius.ll }}} |
|
Line 40: | Line 64: |
{{{ /bin/cp -f SHAPE.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp > radius.ll }}} |
|
Line 43: | Line 74: |
== Sigmas == {{{ convertLatLon tmp 1 > sigma.ll }}} == Number of Images == == Slope == == Albedo == {{{ /bin/cp -f SHAPEA.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp 1 > albedo.ll }}} -------- = Regional = Run backplanesGMT (which is in ORExSPCsupport) (or backplanesISIS) == Topography == |
|
Line 53: | Line 118: |
== Photometric Data == This includes emission angle, incidence angle and phase angle. Run phasei (located in ORExSPCsupport) |
|
Line 55: | Line 127: |
= Regional = | = Make Pretties = |
Line 57: | Line 129: |
== Topography == | {{{ gmt begin GMT_cont gmt set GMT_THEME cookbook gmt grdcontour test.nc gmt end show #gmt grdcontour test.nc -C10 -A50 |
Line 59: | Line 136: |
== Sigmas == | gmt begin GMT_img gmt set GMT_THEME cookbook #gmt makecpt -Crainbow gmt grdimage test.nc -JM6i -B -BWSnE gmt colorbar -DJTC -Bxa -By+lm gmt end show |
Line 61: | Line 143: |
== Number of Images == | |
Line 63: | Line 144: |
== Slope == | gmt begin GMT_img gmt makecpt -Crainbow gmt set GMT_THEME cookbook gmt grdimage test.nc $argR -JM6i -B -BWSnE #gmt colorbar -DJTC -I0.4 -Bxa -By+lm gmt colorbar -DJTC -Bxa -By+lm gmt end show |
Line 65: | Line 152: |
== Albedo == | }}} |
Overview
Archiving SPC-focus products. There are two major categories (global and regional) with two subdivisions (topography and backplanes -- aka other data). The regional products are derived from bigmaps.
Tools
These tools are stored in GitHub as [archivingTools|https://github.com/StereoPhotoClinometry/archivingTools].
convertLatLon - Reads an X,Y,Z (or X,Y,Z,V) list and converts it into lat, lon, distance or lat, lon, value
backplanesGMT - Reads a bigmap and generates a list of latlon, radius, albedo, slope. Use paste latlon radius > tmp for GMT.
backplanesISIS - Reads a bigmap and generates 2D matrix for radius, albedo, slope, lat, lon
phasei - Reads a bigmap and image (the SUMFILE) and generates 2D matrix for incidence, emission, phase angles
- coverage_p_low -- Reads all images and maplets. Creates a pgm that has every pixel marked for every image that it has, counter is only 1.
map_coverage_p_low - Counts the number of images that cover each bigmap pixel. Creates MAPNM-cov.txt and MAPNM-cov.pgm
Grid to GeoTIF
var=radius.ll new=big3 argR="-R280/287/-9/-3" argR="-R0/360/-90/90" argI="-I0.04" date echo "Blockmean" gmt blockmean $var $argR $argI > out echo "sphereinterpolate" cat out | gmt sphinterpolate $argR $argI -Gtest.nc echo "gdal translate" gdal_translate -of GTiff -b 1 -a_srs 60300.prj NETCDF:test.nc $new.tif
Global
Global ICQ resolution to bin sizes [degrees]
ICQ Q Value |
Vertices |
Bin Size [deg] |
512 |
1,579,014 |
0.041 |
256 |
396,294 |
0.164 |
128 |
99,846 |
0.649 |
64 |
25,350 |
2.556 |
32 |
6,534 |
9.917 |
Topography
/bin/cp -f SHAPE.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp > radius.ll
Sigmas
convertLatLon tmp 1 > sigma.ll
Number of Images
Slope
Albedo
/bin/cp -f SHAPEA.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp 1 > albedo.ll
Regional
Run backplanesGMT (which is in ORExSPCsupport) (or backplanesISIS)
Topography
Sigmas
Number of Images
Slope
Albedo
Photometric Data
This includes emission angle, incidence angle and phase angle.
Run phasei (located in ORExSPCsupport)
Make Pretties
gmt begin GMT_cont gmt set GMT_THEME cookbook gmt grdcontour test.nc gmt end show #gmt grdcontour test.nc -C10 -A50 gmt begin GMT_img gmt set GMT_THEME cookbook #gmt makecpt -Crainbow gmt grdimage test.nc -JM6i -B -BWSnE gmt colorbar -DJTC -Bxa -By+lm gmt end show gmt begin GMT_img gmt makecpt -Crainbow gmt set GMT_THEME cookbook gmt grdimage test.nc $argR -JM6i -B -BWSnE #gmt colorbar -DJTC -I0.4 -Bxa -By+lm gmt colorbar -DJTC -Bxa -By+lm gmt end show