发一个不超出边界的画站点图的grid.grd，不知道有用没

wtf

本人在用站点资料绘图的时候，发现按照GrADS书上的Fortran程序生成的格点（grid.grd），在出图之后有超出中国边界的情况，和书上的例图一样，对这个问题我在本站搜了一下相关的解决办法，用maskout函数或者basemap是可以解决的，但是后来从老师那里拷了一个现成的grid.grd,画出来的图就没有超出边界的情况，用记事本打开和程序生成的对比一上，发现这个里面是有空白行的，应该是把边界外的格点直接删除了？好像没人发过，挂上来各位看看，效果图一张 缺点就是边界的地方会有锯齿 不够美观 grid.grd (11.37 KB, 下载次数: 57)

2012-11-22 12:49 上传

点击文件名下载附件

river

不错哦，不过还是没明白这个grid.grd和平时用的那个有什么不一样

chenggao

在arcgis中校正或平滑应该可以吧

ayzqs

似乎还是超出了一点点?

wtf

river 发表于 2012-11-22 20:54

                               
登录/注册后可看大图

不错哦，不过还是没明白这个grid.grd和平时用的那个有什么不一样

额 那可能是我out了，就是用直接用grads书上那个fortran程序生成的grid的话，画出来的图是会超出边界。。。这个应该是把grid.grd文件里边界外的格点删掉了，这样在插值的时候就不会到边界外面  

wtf

ayzqs 发表于 2012-11-23 09:04

                               
登录/注册后可看大图

似乎还是超出了一点点?

额 这个比较简单粗暴，精度不够  

river

wtf 发表于 2012-11-23 09:45

                               
登录/注册后可看大图

额 那可能是我out了，就是用直接用grads书上那个fortran程序生成的grid的话，画出来的图是会超出边界。。 ...

你应该把生成这个文件的程序也发上来，这样才能看出来区别。而且看你那个插值的图也没有反映出你这个文件和书上有什么不同。我就用书上的画的，也都没有超出边界啊。我觉得主要还是gs文件的编写问题吧。

chengxf

GrADS——没学过，也不懂。但是这种情况在Surfer中也同样存在。所以发表些拙见，供参考。
图中等值线着色后出现某些色块超出边界，某些色块没有达到边界，呈锯齿、马赛克状态。其实这种现象只要在生成.grd文件时注意一下某些参数的设置，就可以解决。GrADS中该设置什么参数，我不知道。在Surfer中，只要修改“网格线几何特征”中的“节点数”，就可以解决。具体如何修改，说出来可能对用GrADS作图的朋友没什么用，但可参考论坛的“Surfer”板块中Surfer11 帮助里的教程(第五课)的内容。参见下图。
     
所以，楼主好心发一个这样的.grd文件，对其他朋友是没有什么用的，只要 Z 值改变了，就得生成新的.grd文件，新的等值线肯定和你这个不一样。还是要从根本上解决。

wtf

river 发表于 2012-11-23 12:09

                               
登录/注册后可看大图

你应该把生成这个文件的程序也发上来，这样才能看出来区别。而且看你那个插值的图也没有反映出你这个文件 ...

哦 这样啊 。。。是gs的问题么 ，这张是我按照书上出的grd画的图，gs我也贴这吧 求大神指教 多谢多谢

reinit
'OPEN e:\helloc\GRID.ctl';'OPEN e:\helloc\rr.ctl';'OPEN e:\helloc\jp.ctl'
*'set parea 1 7 1 7'
'set lon 71 136'
'set lat 15 55'
'set  mpdset cnscs'
'set  map  1  1   2'
'draw  map'
'set rgb 23 85 85 83'
'set rgb 24 125 125 122'
'set rgb 25 165 165 161'
'set rgb 26 205 205 200'
'set rgb 28 246 246 237'
'DEFINE A=OACRES(G,rr.2,1.5)';'DEFINE B=MASKOUT(A,G-0.5)'
'DEFINE C=OACRES(G,jp.3,1.5)';'DEFINE D=MASKOUT(C,G-0.5)'
'DEFINE BB=SMTH9(B)';'DEFINE dd=SMTH9(D)'
'set clab forced'
'set grid off'
'set grads off'
'set xlevs 80 90 100 110 120 130'
'set ylevs 20 30 40 50'
'set xlopts  1 5 0.18'
'set ylopts  1 5 0.18'
'set gxout shaded'
*'set cint 10.'
'set clevs -3 -2 -1 1'
'set rbcols 23 24 25 28 25 '
*'set rbcols 11 0 2'
'd smth9(dd)'
*'cbarn 1 0 4 1'
'set gxout contour'
*'set cint 1.0'
'set csmooth on'
'set cterp on'
'set clopts  1 2 0.1'
'set clskip 2'
*'set cthick 6'
'd smth9(dd)'
'set annot 4 12  0.15'
'set strsiz 0.15 0.20'
'draw string 1.2 6 (a)'
'printim e:\1.png white'
*'enable print e:\helloc\rr.gmf'
*'print'
*'disable print'
;

wtf

chengxf 发表于 2012-11-23 19:24

                               
登录/注册后可看大图

GrADS——没学过，也不懂。但是这种情况在Surfer中也同样存在。所以发表些拙见，供参考。
图中等值线着色后 ...

多谢大神指教