信度打点图画不出来，按照之前的精华帖修改的，不知道自己哪里出问题了，求指导，，

LemonTree007

1. 'reinit'
   
2. 'open D:\matlabht\pdsi\pdsi1950-2014.ctl'
   
3. 'enable print section.met'
   
4. 'set grads off'
   
5. 'set grid off'
   
6. 'set lon 70 140'
   
7. 'set lat 10 60'
   
8. 'set lev 0'
   
9. 'set t 1'
   
10. 'set gxout shaded'
    
11. 'set xlint -10'
    
12. 'set ylint -10'
    
13. 'set xlopts 1 4 0.2'
    
14. 'set ylopts 1 4 0.2'
    
15. 'set mpdset cnbasemap'
    
16. 'run define_colors.gs'
    
17. 'set clevs -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5'
    
18. 'set ccols 27 26 25 24 23 22 42 43 44 45 46 47'
    
19. 'cnbasemap ave(pdsi,t=1,t=65)'
    
20. 'cbarn'
    
21. 'run southsea.gs'
    
22. 'draw string 10.1 0.7 sc_PDSI'
    
23. <font color="#ff0000">**红色代码是关于t检验和信度打点的，请帮忙看看是哪里出错了，出来的只有平均分布图，上面没有点，可以确定这段代码是有错误的--
    
24. 'set lon 70 140'
    
25. 'set lat 10 60'
    
26. 'set lev 0'
    
27. 'set t 1'
    
28. 'define ppave = ave(pdsi, t=1, t=65)'  **样本均值
    
29. 'ss = 0'
    
30. while ('t'<=65)
    
31. 'define ss = pow(pdsi(t)-ppave,2)+ss'
    
32. 'set  t  't+1
    
33. endwhile
    
34. 'ss = ss/65'  **样本方差
    
35. 'define test=ppave/sqrt(ss)*sqrt(65)'  **t统计量
    
36. 'set  lat  10 60'
    
37. 'set  lon  70  140'
    
38. 'set  lev 0'
    
39. 'set t 1'
    
40. 'myp=maskout(test,abs(test)-1.67)'
    
41. ****0.1显著性水平
    
42. 'set gxout shp'
    
43. 'set shp -pt shppt'
    
44. 'd myp'
    
45. '!gxps -b 0.10 -c -i section.met -o ttest.ps'</font>
    
46. 'printim D:\matlabht\test\pdsi.png white'
    
47. pull dummy
    
48. 'c'

复制代码

LemonTree007

是因为缺少什么文件吗？错误显示 Invalid command found 11000

ynccjy

我也是按帖子上的方法画不了图，还报错，是不是grads版本问题，或者是没有创建shapefile文件的原因呀？