出图东部和南部沿海地区没有数据是怎么回事？

LemonTree007

gs 1:将年平均pdsi的区域趋势求出并写入grd文件中

1. 'reinit'
   
2. 'open D:\matlabht\pdsi\pdsi1901-2014.ctl'
   
3. 'set fwrite d:\matlabht\pdsi\liner.grd'
   
4. 'set gxout fwrite'
   
5. 'set x 1 720'
   
6. 'set y 1 360'
   
7. 'set z 1'
   
8. 'set t 1 114'
   
9. 'ltrend pdsi output slope rmse'
   
10. 'd slope'
    
11. 'disable fwrite'
    
12. 'reinit'

复制代码

grd出来后，其ctl写为：
dset d:\matlabht\pdsi\liner.grd
undef -9.99e+8
title  1991-2014 annual ltrend sc_PDSI
xdef 720 linear 0 0.5
ydef 360 linear -89.75 0.5
zdef 1   linear 0 1
tdef 1   linear 00Z01JAN1902 1yr
vars  1
lin=>lin 0 y,x  annual ltrend sc_PDSI**和原始数据维度顺序一致。y,x
endvars

gs2：出图脚本

1. 'reinit'
   
2. 'open d:\matlabht\pdsi\liner.ctl'
   
3. 'set parea 0.9 9.2 0.7 8.0'
   
4. 'set grads off'
   
5. 'set grid off'
   
6. 'set mpdset cnbasemap'
   
7. 'set xlopts 1 5 0.2'
   
8. 'set ylopts 1 5 0.2'
   
9. 'set xlint -10'
   
10. 'set ylint -10'
    
11. 'set gxout shaded'
    
12. 'run define_colors.gs'
    
13. 'set clevs -0.08 -0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08 0.1 0.12'
    
14. 'set ccols 26 25 24 23 22 21 41 42 43 44 45 46'
    
15. 'set lon 70 140'
    
16. 'set lat 10 60'
    
17. 'set lev 0'
    
18. 'cnbasemap lin'
    
19. 'cbarn'
    
20. 'run southsea.gs'
    
21. 'printim d:\matlabht\tp\liner.png white'
    
22. 'set parea off'
    
23. 
    
24. 
    
25. 
    
26. 
    

复制代码

出错图片：东南部沿海以及台湾海南无数据

看云淡风清

楼主，你这个scpdsi0.5在哪下的啊，麻烦告知下谢啦

LemonTree007

看云淡风清 发表于 2018-5-13 12:21
楼主，你这个scpdsi0.5在哪下的啊，麻烦告知下谢啦

这个数据是我导师给我下载的，我刚进去那个网站看了一下，我找到的数据集只有1901-2009年份的，而我导师给我的是1901-2014年的，我也不知道怎么回事。
http://catalogue.ceda.ac.uk/?q=s ... p;sort_by=relevance

看云淡风清

谢谢楼主，我有一个1901-2016年的：https://crudata.uea.ac.uk/cru/data/drought/，不知道这套数据集怎么样？楼主用过吗

LemonTree007

看云淡风清 发表于 2018-5-14 20:39
谢谢楼主，我有一个1901-2016年的：https://crudata.uea.ac.uk/cru/data/drought/，不知道这套数据集怎么样 ...

这个我实在做不出评价，我是新手，第一次做论文，数据还是老师给找的，所以没办法评价数据的优劣。好与不好您用用看嘛

看云淡风清

LemonTree007 发表于 2018-5-16 10:24
这个我实在做不出评价，我是新手，第一次做论文，数据还是老师给找的，所以没办法评价数据的优劣。好与不 ...

嗯嗯，数据只到09年似乎老了些，我还是倾向于用较新的数据，祝你好运啦