• Archiving SPC

Clear message

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 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 coverage_p_low
• Use coverage_p.ll as the input to GMT

  • var=radius.ll
    new=tCov
    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

  •   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

• bustGrid MAPNM-cov.txt > new.txt

• paste latlon.txt new.txt > new.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