|
Size: 2990
Comment:
|
Size: 3671
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 86: | Line 86: |
| Line 87: | Line 89: |
| Line 108: | Line 113: |
| * Run backplanesGMT {{{ paste lonlat MAPNM-r > tmp.ll var=tmp.ll argR="-R0/360/-90/90" argI="-I0.04" 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 echo "Done" `date` }}} == Slope == * See topography == Albedo == * See topography |
|
| Line 111: | Line 139: |
| * Run backplanesGMT * Run std2isis * Run isis2ascii * unpaste file * paste lonlat <new> > sigma.ll * Continue with topography script |
|
| Line 113: | Line 149: |
| == Slope == | |
| Line 115: | Line 150: |
| == Albedo == |
* run map_coverage_p_low * unpaste MAPNM-cov.txt * Continue with topography script |
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
- Run backplanesGMT
paste lonlat MAPNM-r > tmp.ll var=tmp.ll argR="-R0/360/-90/90" argI="-I0.04" 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 echo "Done" `date`
Slope
* See topography
Albedo
* See topography
Sigmas
- Run backplanesGMT
- Run std2isis
- Run isis2ascii
- unpaste file
paste lonlat <new> > sigma.ll
- Continue with topography script
Number of Images
- run map_coverage_p_low
- unpaste MAPNM-cov.txt
- Continue with topography script
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