Differences

This shows you the differences between two versions of the page.

--- code:work:lvisf:2023:gmt_gridding:start [2023/03/16 17:19] – created david
+++ code:work:lvisf:2023:gmt_gridding:start [2023/03/16 17:22] (current) – david
@@ Line 2: / Line 2: @@
   * Use GMT to grid some point clouds
 <code>
@@ Line 11: / Line 9: @@
 lvis_20221114_191400_00000159.tif
 lvis_20221114_191430_00000160.tif
+</code>
+<code>
+# mby.cpt by M. Burak YIKILMAZ
+#
+# (original title gmt_color.cpt)
+#
+# COLOR_MODEL = RGB
+-150    24  140 205     -100    135 206 250
+-100    135 206 250     -90     176 226 255
+-90     176 226 255     -50     0   97  71
+-50       0  97   71     50     16  123 48
+      16  123 48      500     232 214 125
+     232 214 125     1200    163 68  0
+    163 68  0       1800    130 30  30
+    130 30  30      2800    161 161 161
+    161 161 161     4000    206 206 206
+    206 206 206     5000    255 255 255
+B       -
+F       -
+N       128     128     128
+</code>
+<code>
+export GMT_CPT_ELEVATION="/home/lvis/repos/lvisf/rtp/scripts/v2/gedi.cpt"
+export GMT_CPT_HEIGHT="/home/lvis/repos/lvisf/rtp/scripts/v2/ht.cpt"
+export GMT_GRID_NODATA="-ENaN"
+export GMT_GRID_SIZE="-I10e/10e"
+export GMT_GRID_SEARCH="-S100e/100e"
+rtp_ql2filter lvis_20221114_200600_00000263.ql2 lvis_20221114_200600_00000263.geoloc.h5 lvis_20221114_200600_00000263.b3d
+REGION=$(gmt info lvis_20221114_200600_00000263.b3d -I0.01 -bi3d)
+gmt nearneighbor lvis_20221114_200600_00000263.b3d -N4+m2 -bi3d -Glvis_20221114_200600_00000263.grd $GMT_GRID_SIZE $GMT_GRID_SEARCH $GMT_GRID_NODATA -V $REGION
+gmt grdgradient lvis_20221114_200600_00000263.grd -A-45 -Glvis_20221114_200600_00000263_I.grd
+gmt grdimage -C$GMT_CPT_ELEVATION -Q -JQ267.66/29.34/5i lvis_20221114_200600_00000263.grd -Ilvis_20221114_200600_00000263_I.grd -Alvis_20221114_200600_00000263.tif
 </code>