Differences between revisions 2 and 7 (spanning 5 versions)
Revision 2 as of 2021-06-16 12:17:14
Size: 677
Editor: EricPalmer
Comment:
Revision 7 as of 2021-06-16 15:03:38
Size: 1880
Editor: EricPalmer
Comment:
Deletions are marked like this. Additions are marked like this.
Line 25: Line 25:

Line 30: Line 32:
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 40:
   {{{
   convertLonLat > radius.ll

   }}}
Line 40: Line 45:
   {{{
/bin/cp -f SHAPE.TXT tmp
vi tmp
1,$s/D/E/g
ZZ
convertLatLon tmp > radius.ll
   }}}
Line 43: Line 55:



== 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 ==
Line 53: Line 99:

== Photometric Data ==

This includes emission angle, incidence angle and phase angle.

Run phasei (located in ORExSPCsupport)
Line 55: Line 108:
= Regional = = Make Pretties =
Line 57: Line 110:
== Topography == {{{
gmt begin GMT_cont
gmt set GMT_THEME cookbook
gmt grdcontour test.nc
gmt end show
#gmt grdcontour test.nc -C10 -A50
Line 59: Line 117:
== 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 61: Line 124:
== Number of Images ==
Line 63: Line 125:
== 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 65: Line 133:
== Albedo ==
}}}

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)