Grads显著性区域打点

a1s2d3 · 发表于 2022-5-4 15:25:14

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 a1s2d3 于 2022-5-4 15:28 编辑

用Fortran计算标准化后相关系数、回归系数、t值，再用grads绘制环流场图并打点打点方式主要参考http://t.csdn.cn/pjAdO上用的打点方法需要绘制显著性风场时要先用wnd文件将通过显著性t检验的风场放入一个文件中，未通过显著性检验的风场按缺时处理

huanliuchang.gs (1.55 KB, 下载次数: 21)

伽蓝鸟 · 发表于 2022-5-4 17:07:40

set tile
This command, which is available beginning with GrADS version 2.1, is used to create a new pattern fill. Patterns are defined in the form of "tiles" which are small images that are used repeatedly to fill an area. GrADS supports the automatic generation of tiles for drawing certain simple patterns. More complex patterns can be acheived by providing a custom image file.

a1s2d3 · 发表于 2022-5-4 22:24:45

伽蓝鸟发表于 2022-5-4 17:07
set tile
This command, which is available beginning with GrADS version 2.1, is used to create a new ...

感谢大佬建议，用这个显著性打点确实比循环着一个一个画方便

Bone · 发表于 2022-6-2 22:27:28

本帖最后由 Bone 于 2022-6-2 22:29 编辑

楼主，我按照你的代码改了，我的目标是把大于1的区域全部打点，但是却全部都打上了点，您能帮忙看一看吗？
'reinit'
'open D:/FG/height/hgt/200.ctl'
'set lon 60 180'
'set lat 0 45'
'set lev 200'
'set t 1'
'define hgtclim=ave(hgt,t=1,t=31)'
'set t 31'
'define a=hgt-hgtclim'
'set grads off'
'set grid off'
'set xlint 10'
'set ylint 15'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.18'
'set gxout shaded'
'run define_colors.gs'
'set clevs -6 -4 -3 -2 -1 0 1 2 3 4 6 8 10 12'
'set ccols 47 46 45 43 42 41 21 22 23 24 25 26 27 28 29'
'd a/10'
'set gxout contour'
'set t 31'
'set cint 4'
'set ccolor 1'
'd hgt/10'

*打点
'open D:/FG/height/hgt/200dev.ctl'
'set lon 60 180'
'set lat 0 45'
'set lev 200'
'set t 1'
y2 = 0
while(y2<=45)
x2 = 60
while(x2<=180)
'q w2xy 'x2' 'y2''
abx = subwrd(result,3)
aby = subwrd(result,6)
'set lon 'x2''
'set lat 'y2''
'd dev.2'
aaa = subwrd(result,4)
if(aaa<0)
aaa = -aaa
endif
if(aaa>=1)
'set string 0 c 6'
'set strsiz 0.02'
'draw string 'abx' 'aby' `3D'
endif
x2 = x2+2.5
endwhile
y2 = y2+2.5
endwhile

'cbarn 1.25 0'
'printim D:\FG\height\hgt\result2\21hgt200.png white'
;

a1s2d3 · 发表于 2022-6-4 17:16:10

我感觉没有问题呀，会不会是200dev里面的数据有问题？

		自动登录	找回密码
密码			立即注册

[分享资料] Grads显著性区域打点

登录后查看更多精彩内容~