Differences between revisions 1 and 9 (spanning 8 versions)
Revision 1 as of 2021-06-16 12:15:15
Size: 581
Editor: EricPalmer
Comment:
Revision 9 as of 2021-06-17 12:57:18
Size: 2990
Editor: EricPalmer
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
= Overview =
Line 2: Line 3:
= Global = 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 =
Line 26: Line 44:


------


= 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)

Line 37: Line 118:

== Photometric Data ==

This includes emission angle, incidence angle and phase angle.

Run phasei (located in ORExSPCsupport)
Line 39: Line 127:
= Regional = = Make Pretties =
Line 41: 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 43: 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 45: Line 143:
== Number of Images ==
Line 47: 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 49: 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

Archiving SPC (last edited 2023-05-01 08:35:29 by JohnWeirich)