|
GMT 画彩色图时裁剪有效部分
grdimagecut.txtmakecpt
cut.txt
此文件是对整个区域的判断，使用的方法是，对有效范围标记为NAN值，对要裁剪的无效部分标记为1.0。 格式为(longitude, latitude, mark)
makecpt
gmt makecpt -Cjet -T2.0/5.0/1.0 -I -Z > cut.cpt
其中，jet对应的color种类，2.0/5.0/1.0对应的是显示的颜色的上下界面和分段尺度（不重要），此项目的是使得cut.cpt的背景色（B）为white，且保证无效部分（1.0）不在cpt范围内。
grdimage
gmt xyz2grd cut.txt -Gcut.grd -R$R -I0.1
gmt grdimage cut.grd -J$J -R$R -Ccut.cpt -Q -Z -K -O >> $ps
其中，xyz2grd是将xyz格式的文件转换为grd文件，grdimage是在已经存在彩图的基础上进行裁剪，方法就是利用cut.grd的内容，跳过有效范围（-Q），将无效范围（1.0） 画成白色（cpt的背景色）。
|