• Use GMT to grid some point clouds

Example files:
lvis_20221114_191300_00000157.tif
lvis_20221114_191330_00000158.tif
lvis_20221114_191400_00000159.tif
lvis_20221114_191430_00000160.tif

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