GMT加载遥感影像背景（1）

在GMT中，我们可以加载nc或者grd等格网数据，而遥感影像一般是geotiff格式，因此无法在GMT中直接加载遥感影像。本人也一直被这个问题困扰，直到最近看到的一篇微信推送才解决了这个问题。这里我主要介绍两类方法：多波段组合和单波段加载。本节主要介绍多波段的加载。

【1】多波段组合（已经实现）

所需数据：经过前期处理好，经过波段合成的真彩色影像。

下一步需要将其进行坐标转换，即将其转换成wgs84地理坐标。有很多的方法可以实现这一转换，但是我在这里像介绍利用gdal_translate工具直接转换。我在这里使用时gmt5版本，其中没有包含这一函数。而gmt6新版本中包含了这一函数，可以直接调用

坐标转换的命令：

gdalwarp -t_srs EPSG:4326 test.tif test_wgs.tif

test.tif为转换前图像名称，test_wgs.tif是转换后的影像名称，直接修改即可。

另外如果需要进行小区域的裁剪，可以使用一下命令：

gdalwarp -t_srs EPSG:4326 -te 116 40 117 41.2 test_wgs.tif teat_wgs_out.tif

需要修改的是范围大小，输入范围框的左下角和右上角经纬度坐标。

然后利用一下的代码即可完成遥感影像的加载：

set R=116/117/40/41.2

set J=M12c

set PS=map2.ps

gmt gmtset FORMAT_GEO_MAP=ddd:mm:ss

gmt psbasemap -R%R% -J%J% -Bpxa50mf -Bpya50mf -BnWSe -Xc -Yc  -K >%PS% 

gmt pscoast -R%R% -J%J% -O -A500 -W -I1 -K  -Lx1.0i/0.15i+c120/36+w5k >>%ps%

gmt grdimage  teat_wgs_out.tif+b0  teat_wgs_out.tif+b1 teat_wgs_out.tif+b2 -J%J% -R%R% -K -O >> %ps%

gmt psxy ucas.txt -J%J% -R%R% -Sc0.05c -W1,red -O -K>> %ps%

REM # Add jpg

REM gmt psimage app.png -Dx0.3/6.5+w11.4c -O -K >> %ps%

REM # Add small map in Up-Right corner

REM gmt pscoast -Rg -JG120/41/2c  -Dc -A5000 -Gpink -Sblue  -O -X2.5c -Y6.5c -K >>%ps%

REM echo  120:26:30 36:02:57 | gmt psxy -Sa0.1c -Gred -R -J -O  >> %ps%

REM gmt psconvert %ps% -A -P -Tg

创作不易，期待大家的评论收藏分享！！