博文
【GMT】pstext 通过文件绘制脚本指定区域外的文字
||
最近在用pstext绘制标注文本时遇到了一个问题，如果在脚本中添加标注的范围超过了脚本第一条指定绘图区域的范围，发现标注无法绘制到图形上，如下图所示：
绘图原意是希望能够在图形的上方为图形添加一个题为“Y=sin（x）+zhanjinggang”的标题。因为该数据的覆盖范围是全球覆盖，因此在指定pstext 的第一和第二个参数时，必须在全球的经纬度范围内，否则标签内容超过范围，无法显示。
为了解决这个问题，找了半天解决方案也没有弄成功，后来发现如果利用文本文件指定标签标题内容，然后再用pstext绘制标题内容的话，可以满足绘图要求（如下图黑色标签所示）。
脚本源码：
del .gmt* gmt.*
REM #set the GMT variable for dots per inch (resolution)
gmtset DOTS_PR_INCH 600
REM #set the GMT variable for dots per inch (resolution)
gmtset DOTS_PR_INCH 600
GMTset FRAME_WIDTH 0.03i
gmtset FRAME_PEN 1.50p
gmtset ANNOT_OFFSET_PRIMARY 0.02i
gmtset FRAME_PEN 1.50p
gmtset ANNOT_OFFSET_PRIMARY 0.02i
gmtset HEADER_FONT_SIZE 6
gmtset LABEL_FONT_SIZE 6
gmtset LABEL_FONT_SIZE 6
gmtset ANNOT_FONT_PRIMARY Times-Roman
gmtset ANNOT_FONT_SECONDARY Times-Roman
gmtset HEADER_FONT Times-Roman
gmtset LABEL_FONT Times-Roman
gmtset TICK_LENGTH 0
gmtset ANNOT_FONT_SECONDARY Times-Roman
gmtset HEADER_FONT Times-Roman
gmtset LABEL_FONT Times-Roman
gmtset TICK_LENGTH 0
REM #set the GMT variable for eps paper size Custom_WxH
REM gmtset PAPER_MEDIA Custom_42ix45i
gmtset PS_MEDIA = 10ix12i
REM gmtset PAPER_MEDIA Custom_42ix45i
gmtset PS_MEDIA = 10ix12i
set outps=zhenghui.ps
grdreformat 2008-1.grd temp0.grd=nf -V
rem makecpt -Chaxby -T-6/6/1 -D >t.cpt
rem makecpt -Crainbow -T-6/6/1 -Z >t.cpt
rem makecpt -Chaxby -T-6/6/1 -D >t.cpt
rem makecpt -Crainbow -T-6/6/1 -Z >t.cpt
grdimage temp0.grd -R0/360/-90/90 -E200 -JQ10i -B60g60/50f50WesN -X0.6i -Y2i -K -V -Ct.cpt >%outps%
pscoast -R -JQ -Di -A10000 -W1p,white -V -K -O -P >>%outps%
rem 通过在脚本中添加标注来绘制到图形中，当标注范围超过图形范围是，无法绘制
echo 120 90 30 0 7 LM y=sin(x)+zhanjingang | pstext -G0/0/255 -R -J -K -O >> %outps%
echo 120 90 30 0 7 LM y=sin(x)+zhanjingang | pstext -G0/0/255 -R -J -K -O >> %outps%
rem 通过文件读取来绘制绘图区域外的标注
pstext label2.txt -R -J -K -O -N -Gblack -D0.04c >> %outps%
pstext label2.txt -R -J -K -O -N -Gblack -D0.04c >> %outps%
rem 水平色棒
psscale -D5i/-0.2i/9i/0.08ih -E1 -Ct.cpt -I -B2/:@~m@~gal: -O -V -K >>%outps%
rem 垂直色棒
psscale -D10.5i/2.5i/4.1i/0.08i -E1 -Ct.cpt -I -B2/:@~m@~gal: -O -V >>%outps%
psscale -D10.5i/2.5i/4.1i/0.08i -E1 -Ct.cpt -I -B2/:@~m@~gal: -O -V >>%outps%
echo converting to jpg
ps2raster %outps% -E600 -V -Tj
convert %outps% -density 600 -rotate 90 zhenghui.jpg
pause
https://wap.sciencenet.cn/blog-51026-642603.html
上一篇：【GMT】ps2raster命令设定GMT绘图图片输出格式
扫一扫，分享此博文
全部精选博文导读
相关博文
- • Minerals线下恳谈会：履践致远、与时偕行——对话中国科学院广州地球化学研究所期刊合作学者
- • 聚英才 建高地 | 北京理工大学“特立青年学者”全球招聘开启
- • 700年后日本或濒临灭绝？日本学者推算预测：届时或仅剩1名15岁以下孩子
- • [转载]【同位素视角】非英语母语学者如何区分’e.g.’, ‘i.e.’, ‘namely’与‘such as’等混淆难题
- • 美国佐治亚大学等机构学者：刈割策略对Bulldog 805紫花苜蓿+Tifton 85狗牙根混播草地产量及品质的影响
- • 美国堪萨斯州立大学、密苏里大学等机构学者研究成果：土壤水分管理策略和品种多样性对紫花苜蓿产量、营养品质和农场盈利能力的影