GMT绘图学习记录（2）

今天主要记录一下如何下载卫星影像以及如何以卫星影像为底图用GMT绘图 

1.卫星影像下载

我是在图新地球下载的，具体下载方法参考b站一分钟教你下载高清卫星图！_哔哩哔哩_bilibili

2.坐标系转换

因为下载好的tiff是web墨卡托投影坐标系，必须要将其转换为地理坐标系WGS84下才是经纬度坐标。将下载好的tiff导入arcgis，点击图层-坐标系-地理坐标系-WGS-84

注意：这一步不能用投影栅格做转换，我之前做出来的tif无法导入GMT，不知道什么原因。

3.导出tiff 

特别注意 空间参考要选当前数据框

弹出的选项框，选NO

 4.GMT绘图

gmt begin sat_plot png E1000
    R="117.83/117.88/36.49/36.53"
    gmt basemap -JM8c -R$R -Bxa1mf1m -Bya1mf1m -BWSne
	gmt grdimage ppp.tif 
	# Place modern session commands here
gmt end show

 

完成！！！！