Differences between revisions 2 and 10 (spanning 8 versions)

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

-  ⇤ ← Revision 2 as of 2021-06-16 12:17:14 → 
  Size: 677
  Editor: EricPalmer
  Comment:
+   ← Revision 10 as of 2021-06-17 13:09:59 → ⇥
  Size: 3671
  Editor: EricPalmer
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 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 39:
+Line 63:
+   {{{
/bin/cp -f SHAPE.TXT tmp
vi tmp
1,$s/D/E/g
ZZ
convertLatLon tmp > radius.ll
   }}}
-Line 46:
+Line 79:
+   {{{
convertLatLon tmp 1 > sigma.ll 
   }}}
-Line 47:
+Line 85:
-Line 50:
+Line 90:
-Line 51:
+Line 94:
+   {{{
/bin/cp -f SHAPEA.TXT tmp
vi tmp
1,$s/D/E/g
ZZ
convertLatLon tmp 1 > albedo.ll
   }}}
-Line 55:
+Line 105:
-Line 57:
+Line 108:
+Run backplanesGMT (which is in ORExSPCsupport) (or backplanesISIS)
-Line 58:
+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 61:
+Line 139:
+ * Run backplanesGMT
 * Run std2isis
 * Run isis2ascii
 * unpaste file
 * paste lonlat <new> > sigma.ll
 * Continue with topography script
-Line 63:
+Line 149:
-== Slope ==
-Line 65:
+Line 150:
-== Albedo ==
+ * 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


}}}