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.

• 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

• bustGrid - Takes an ascii grid output and puts each entry on its own line. Needed to take grid and paste it with lat/lon to make grid

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 support/60300.prj  NETCDF:test.nc $new.tif

Global

Global ICQ resolution to bin sizes [degrees]

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 

• Use sigma.ll as the input to GMT

Number of Images

• Run globalStat
• Use global-cov.ll as the input to GMT

  • var=global-cov.ll
    new=tCov
    argR="-R0/360/-90/90"
    argI="-I1"

Best Resolution

• Run globalStat
• Use global-res as the input to GMT

  • var=global-res.ll
    new=tRes
    argR="-R0/360/-90/90"
    argI="-I1"

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

  •   # Make the backplanes
      name=LEADEQ
      echo $name | ~/bin/backplanesGMT

Output:

 Version:   1.50000000    
 Input map name (only 6 char no MAPFILES and .MAP)
 BIGMAP position
     V   120.38036346435547       -510.48150634765625       -54.750602722167969     
     Ux  -2.9186813160777092E-002  0.12163921445608139      -0.99214518070220947     
     Uy  0.97239929437637329       0.23332308232784271        0.0000000000000000     
     Uz  0.23149037361145020      -0.96476125717163086      -0.12509183585643768     
     QSZ          99
 Lon   279.612366       286.917664    
 Lat  -9.55961323      -2.35706925    

We'll need the Lon and Lat range for "argR", and QSZ will give an idea of what you need for "argI".

•   # Convert them into lon/lat/val
    paste $name-lonlat.txt $name-r.txt > tmp.ll
    var=tmp.ll
    argR="-R279.612366/286.917664/-9.55961323/-2.35706925"
    argI="-I199+n"
  
    # Run the rebinning/interpolating routines
    echo "Blockmean"
    gmt blockmean $var $argR $argI  > out
    echo "sphereinterpolate"
    cat out |  gmt sphinterpolate $argR $argI -Gtest.nc

Now you need some math to convert the Lon range number of nodes (i.e. "argI") to deg/px. (286.917664 - 279.612366) / 199 = 0.03671

•   # Convert to geoTiff
    echo "gdal translate"
    degPx=0.0367
    gdal_translate -of ISIS3 -tr $degPx $degPx -r bilinear -b 1 -co TARGET_NAME=Tethys -co DATA_LOCATION=GEOTIFF -a_srs support/60300.prj NETCDF:test.nc $name-r-ISIS.lbl

If you don't want to make it readable by ISIS, and don't care about square pixels, you can run the following

•      gdal_translate -of GTiff -b 1 -a_srs support/60300.prj  NETCDF:test.nc $name-r.tif

prj files can be found here

How to trim size

• Use ArcMap to find the new lat and lon ranges

• Remake geotiff with new lat and lon ranges
• Measure the size of the lat range. Divide this number by the GSD from SPC to get number of pixels.
• Divide number of pixels by lon range to get degPx
• Update config.txt with these new values (i.e. lat and lon ranges, number of pixels across, degPx)
• Make geotiff again

Slope

* See topography

Albedo

* See topography

Sigmas

• Run backplanesGMT (this gets lon/lat only file)
• Take the SIGMAS.pgm and convert it to lon/lat/val files

  •    name=LEADEQ
       std2isis from=SIGMAS.pgm to=tmp.cub
       isis2ascii from=tmp.cub to=tmp1.txt
       sed '1,2d' tmp1.txt > tmp2.txt
        ~/bin/bustGrid tmp2.txt > tmp3.txt
       paste $name-lonlat.txt tmp3.txt > sigma.ll
       var=sigma.ll

• Continue with topography script

Number of Images

• run map_coverage_p_low
• Convert coverage out to lon/lat/value

  •    name=LEADEQ
       bustGrid $name-cov.txt > tmp1.txt
       paste $name-lonlat.txt tmp1.txt > imgNum.ll
       var=imgNum.ll

* 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

Setup

Use Dropbox spcShare directory

• cd /opt/local/spc/
  relink ~/Dropbox/spcShare/Tethys .

Then in the working directory

• relink.sh /opt/local/spc/Tethys/support .
  relink.sh /opt/local/spc/Tethys/mapConfig .