使用idw插值画图出现问题

半颗麦芽糖

王老师您好，我今天使用自己处理好的站点数据画图的时候发现图中表现出来的值皆为正值，但是数据本身以负值为主，附上数据和脚本，请王老师指导，谢谢！ mod1.txt (2.15 KB, 下载次数: 4)

2017-4-19 18:26 上传

点击文件名下载附件

1. 
   
2. <P> fn = 'F:/tuer/test/mod1.txt'
   
3. ncol = numasciicol(fn)
   
4. nrow = numasciirow(fn)
   
5. a = asciiread(fn,shape=(nrow,ncol))
   
6. lat = a[:,0]
   
7. lon = a[:,1]
   
8. co = a[:,2]
   
9. x = arange(97, 111, 0.5)
   
10. y = arange(21, 36, 0.5)
    
11. gco,gx,gy = griddata((lon, lat), co, xi=(x, y), method='idw', pointnum=4)
    
12. axesm()
    
13. lcq = shaperead('F:/11/MeteoInfo/map/chongqing/Chongqing_province.shp')
    
14. lsc = shaperead('F:/11/MeteoInfo/map/sichuan/Sichuan_province.shp')
    
15. lyn = shaperead('F:/11/MeteoInfo/map/yunan/Yunnan_province.shp')
    
16. lgz = shaperead('F:/11/MeteoInfo/map/guizhou/Guizhou_province.shp')
    
17. geoshow(lcq,edgecolor='k')
    
18. geoshow(lsc,edgecolor='k')
    
19. geoshow(lyn,edgecolor='k')
    
20. geoshow(lgz,edgecolor='k')
    
21. layer = contourfm(x, y, gco)
    
22. colorbar(layer)
    
23. xlim(97,110)
    
24. ylim(21,35)</P>

复制代码

MeteoInfo

刚测试了一下，插值后还是已负值为主，不过负值都比较小，正值比较大。

半颗麦芽糖

MeteoInfo 发表于 2017-4-19 18:49
刚测试了一下，插值后还是已负值为主，不过负值都比较小，正值比较大。

感谢王老师的回复，刚才看了一下，确实如您所说，修改了色标间隔以后就可以了，再次感谢！