关于用grads做显著性检验，阴影出现马赛克的问题

肚娘

用的是2.5*2.5的nc格式资料，两个样本场eu1 eu2，通过显著性检验地区用阴影表示，出现了边缘呈马赛克状态，请问这样的马赛克可以改善吗？我看了论坛其他求助，我也添加了smth9做平滑但好像并没有什么改善。。还听说了'set csmooth on'，但是我不知道应该放在哪，
求大神帮忙看看，能否改善马赛克问题？




前面先求了样本平均值eu1，eu2，样本数量为21和18，
接下来求方差sss1和sss2，做显著性检验的gs文件如下：
'define o1=pow((oo1-eu2),2)'
'define p1=pow((pp1-eu2),2)'
'define q1=pow((qq1-eu2),2)'
'define r1=pow((rr1-eu2),2)'
'define sss2=(a1+b1+c1+d1+e1+f1+g1+h1+i1+j1+k1+l1+m1+n1+o1+p1+q1+r1)/18'
'define ooo=sqrt(((20*sss1)+(17*sss2))/37)'
'define ppp=ooo*sqrt((1/21)+(1/18))'
'define ttt=(eu1-eu2)/ppp'

*****显著性检验0.01、0.05、0.1  **** t=2.715、2.026、1.687

'define ttt1=maskout(abs(ttt),abs(ttt)-2.026)'
'define ttt2=smth9(ttt1)'

'set rgb 19 255 194 194'
'set gxout shaded'
'set clevs 0'
'set ccols 19'
'd ttt2'
**************************检验完成************
'set gxout contour'
'd eu1jp-eu2jp'
'set annot 1 3'
'draw xlab (c)'
'printim D:\ziliao\500hpa1\EU33.png white'

肚娘

river 发表于 2017-5-5 20:09
'set csmooth on'放在d 前面就好

我试了一下，感觉图和现在这个没区别..

肚娘

自己顶一下...

river

'set csmooth on'放在d 前面就好

肚娘

river 发表于 2017-5-5 20:09
'set csmooth on'放在d 前面就好

我试了一下，感觉图和现在这个没区别..

高海拔萤火虫

'set gxout shaded’ 后面加'set csmooth on'就可以了

river

肚娘 发表于 2017-5-7 19:03
我试了一下，感觉图和现在这个没区别..

这个应该是maskout函数导致的问题吧。好久不弄也不太知道了，我记得t检验什么的，不是根据样本数什么的最后算出来的显著性检验的值是只有一个那样的，只是分别一正一负。你这个是啥检验啊？

孙亦

'define ttt1=maskout(ttt,abs(ttt)-2.026)'

liy72427

楼主解决了吗？我现在也遇到了这样的问题

xiaopingzjpq

我觉得是数据的分辨率的问题吧，2.5*2.5算出来的可能就是这样的。