MeteoInfoLab脚本：TXT站点数据绘制等值线填色图

静言_GRMC

参考了版主发的脚本示例（http://bbs.06climate.com/forum.p ... 6422&extra=page%3D1），版主的是用micaps资料，我也依葫芦画瓢做了一个，数据存在txt（经、纬度、PM2.5浓度），分享给大家，有代码多余或者不规范的请多多指教~

代码入下：

1. #Get file names
   
2. fn = 'f:/meteoinfoJS/AQI.txt'
   
3. 
   
4. ncol = numasciicol(fn)
   
5. nrow = numasciirow(fn)
   
6. a = asciiread(fn,shape=(nrow,ncol))
   
7. lon = a[:,0]
   
8. lat = a[:,1]
   
9. pm= a[:,2]
   
10. 
    
11. #To grid data
    
12. x = arange(112.9, 114.3, 0.2)
    
13. y = arange(22.5, 24.2, 0.2)
    
14. gtemp,gx,gy = griddata((lon, lat), pm, xi=(x, y), method='idw', radius=0.8)
    
15. #Plot
    
16. axesm()
    
17. bou1_layer = shaperead('f:/MeteoInfo/map/guangzhou_county_new.shp')
    
18. mlayer = shaperead('f:/MeteoInfo/map/guangzhou_county_new.shp')
    
19. geoshow(bou1_layer, edgecolor='lightgray')
    
20. geoshow(mlayer, visible=False)
    
21. 
    
22. levs = [0, 25, 35, 50, 75, 95, 115, 130, 150,200]
    
23. cols = [(255,255,255),(0,255,0),(127,255,0),(255,255,0),(255,215,0),(255,128,0),(255,97,0), \
    
24.     (255,0,0),(176,23,31),(135,38,87),(255,0,255)]
    
25. #layer = contourfm(x, y, gtemp,20)  
    
26. 
    
27. layer = contourfm(x, y, gtemp,levs,colors=cols)
    
28. #slayer = scatterm(lon, lat,pm,colors=['k'], size=10)
    
29. slayer = scatterm(lon, lat,pm,levs,colors=cols, size=8)
    
30. 
    
31. masklayer(mlayer, [layer])
    
32. xlim(112.9, 114.1)
    
33. ylim(22.5, 24)
    
34. title(u'广州市空气质量站点PM2.5浓度',fontname=u'黑体',fontsize=20,bold=False,color='blue')   
    
35. text(113.3, 23.9, u'2017年5月6日13时', fontname=u'黑体', fontsize=16)
    
36. colorbar(layer)
    

复制代码

李小毛123

哈哈，不错不错，加油加油。

MeteoInfo

欢迎这样的经验分享贴！

静言_GRMC

MeteoInfo 发表于 2017-5-16 09:30
欢迎这样的经验分享贴！

谢谢王老师点赞~

︶ㄣ安定■丶

很好很强大

爽歪歪！

请问广州市或者广东省的地图去哪里下载呢？

爽歪歪！

跟大牛要到啦

静言_GRMC

爽歪歪！ 发表于 2017-5-26 14:52
跟大牛要到啦

那就好啊   ~~

静言_GRMC

爽歪歪！ 发表于 2017-5-26 14:52
跟大牛要到啦

那就好啊   ~~

猫猫2003

AQI.txt'能不能把AQI.txt'给我们看下