已解决（grads打点信度检验在二维场上如何实现？）

影伊

画纬度-高度场上的t检验区，之前用的是阴影图呈现，现在需要用信度打点表示。搜索了一下，看到他的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=18475介绍了lon-lat图上画打点，方法很好用。楼主说了画二维的需要手动画，我也留言了。这里我想请教一下大家，用grads画二维的打点具体该怎么实现？（小生这厢有礼了~~~，先道个谢）


先交代一下我的数据吧，算出来的t检验值，1979-2012年共34年，通过90%对应的t值为1.717

c.ctl如下dset E:\data\xiebo\djf\Again\grdnew\feb-cu.grd
undef -999
xdef 1 linear -180.0 1.5
ydef 121 linear -90 1.5
zdef 37 levels 1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
tdef 1 linear 00Z01jan1979 1yr
vars 1
u 37 -999 u wind
endvars

zt.ctl如下
dset E:\data\xiebo\djf\Again\grdnew\feb-utjy.grd
undef -999
xdef 1 linear -180.0 1.5
ydef 121 linear -90 1.5
zdef 37 levels 1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
tdef 1 linear 00Z01jan1979 1yr
vars 1
zt 37 -999 fangcha
endvars
***********以上为2个ctl文件（分别是差值和t检验值）*****************

gs文件
'open 路径\zt.ctl'
'set parea 3 8.5 2 6'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set lat 0 80'
'set lev 150 1'
'set zlog on'
'set t 1'
'set ylevs 150 125 100 70 50 30 20 10 7 5 3 2 1'
'set xlevs 0 10 20 30 40 50 60 70 80'
'set gxout contour'
'set clevs -1.717 1.717'
'set ccols 15 0 15'
'set grid horizontal 5 3'
'd zt '
'reinit';

画了水平格点后，可以直观看到等值线里的格点，应该就是所要的打点吧、、、
我又用了一个笨办法，加了下面几句，手动描点，但是肯定是有误差的。那么，准确的方法该如何画二维打点？
'draw mark 3 4.45 2.87 0.05'
'draw mark 3 4.54 2.87 0.05'
'draw mark 3 4.63 2.87 0.05'
'draw mark 3 4.72 2.87 0.05'
'draw mark 3 4.81 2.87 0.05'
'draw mark 3 4.90 2.87 0.05'
*
'draw mark 3 5.35 2.14 0.05'


我就直接在下面补充了（老师教的画打点图）分享一下。

river

楼主没看到那个帖子下面的附件吗，下载下来看一下吧

影伊

river 发表于 2015-11-19 11:50
楼主没看到那个帖子下面的附件吗，下载下来看一下吧

你的意思是说可以实现？

river

影伊 发表于 2015-11-19 13:09
你的意思是说可以实现？

那个贴子里，别人贴出来的图上不就有打点的区域吗，你照着试试不就行了

影伊

river 发表于 2015-11-19 21:37
那个贴子里，别人贴出来的图上不就有打点的区域吗，你照着试试不就行了

谢谢你的提醒，别人的图上是有打点的区域，是在lon-lat平面上实现的，我也已经试过，是可以画的。我说的是纬度-高度剖面上的打点，那个方法并不适用。

xiaocaoqiqiao

谢楼主分享，很实用

river

影伊 发表于 2015-11-22 14:53
谢谢你的提醒，别人的图上是有打点的区域，是在lon-lat平面上实现的，我也已经试过，是可以画的。我说的 ...

GrADS will draw the contents of a shapefile as an overlay on top of an existing plot with the command draw shp. It is necessary to draw a plot first in order to establish the dimensions and scaling of the display. Shapefiles contain 2-dimensional spatial features, so your plot must be varying in the X-Y (lon/lat) domain.
那应该就是说你这样的目前用这种方法还没办法实现。

我建议你使用Bin Guan教授分享drawmark.gs来画。具体用法自己打开gs看。我不上传附件了，Bin Guan教授分享了好多有用的gs，都是他自己编写的。我给你个地址自己去下载吧。http://bguan.bol.ucla.edu/bGASL.html

影伊

river 发表于 2015-11-22 20:24
GrADS will draw the contents of a shapefile as an overlay on top of an existing plot with the comm ...

好的，谢谢

影伊

river 发表于 2015-11-22 20:24
GrADS will draw the contents of a shapefile as an overlay on top of an existing plot with the comm ...

好的，谢谢

samvelly

{:eb502:}