求助：F检验画格点图

月是故乡明

这是我的gs文件，需要通过对两个数据进行F检验，并在每个格点作出标记，满足F检验和不满足F检验分开标记，个人想法是用循环语句在每隔一个经纬度检验一次并作出标记，但看帖子发现还需要将经纬度映射为坐标值来使用draw marktype mark x y 语句，请问怎样实现呢？
'reinit'
'open d:\wb\050300-111.ctl'
'open d:\wb\050300-211.ctl'
'set parea 1 9.4 1 5.8'
'set z 11'
'set grads off'
'set lon 70 140'
'set lat 15 55'
'set grads off'
'set xlint 10'
'set ylint 10'
'set xlopts 1 5 0.225'
'set ylopts 1 5 0.225'
'set clopts 0 0.1 0.2'
'define zave1=ave(ph+phb,t=5,t=372)'
'define zave2=ave(ph.2+phb.2,t=5,t=372)'
'define s1=ave(pow(zave1-ph-phb,2),t=5,t=372)'
'define s2=ave(pow(zave2-ph.2-phb.2,2),t=5,t=372)'
'define F=s1/s2'
*************** a=0.05 level   *************
if(1<F<1.3|0.7629<F<1)
'draw mark 1'
endif
if(F>1.3|0<F<0.7629)
'draw mark 2'
endif
'printim d:\wb\exam.png white'
'print'
'disable print'

river

这个用fortran更简单吧

月是故乡明

river 发表于 2013-10-22 18:16

                               
登录/注册后可看大图

这个用fortran更简单吧

因为要画图，fortran 也只能用在循环语句中，那坐标转换还是实现不了呢，求各位指导

river

月是故乡明 发表于 2013-10-22 20:54

                               
登录/注册后可看大图

因为要画图，fortran 也只能用在循环语句中，那坐标转换还是实现不了呢，求各位指导

你可以在fortran中循环，然后利用fortran自动写出相应的gs，然后画图。

zl.87625

river 发表于 2013-10-22 22:07

                               
登录/注册后可看大图

你可以在fortran中循环，然后利用fortran自动写出相应的gs，然后画图。

您好，请问您那边有木有F检验的fortran程序？急需学习！

river

zl.87625 发表于 2014-1-7 20:46

                               
登录/注册后可看大图

您好，请问您那边有木有F检验的fortran程序？急需学习！

没有，无论什么的检验书上都有公式。带进去计算一下，也就是一个数值而已啊