====== GMT Gridding ======

  * Use GMT to grid some point clouds

<code>
Example files:
lvis_20221114_191300_00000157.tif
lvis_20221114_191330_00000158.tif
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
50      16  123 48      500     232 214 125
500     232 214 125     1200    163 68  0
1200    163 68  0       1800    130 30  30
1800    130 30  30      2800    161 161 161
2800    161 161 161     4000    206 206 206
4000    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>