MeteoInfoLab脚本示例：区域填色图

MeteoInfo

读取美国区域的温度格点数据和美国州的行政区域shape图层，给图层增加温度字段'temp'，利用维度数组的maskout函数获取每个州的平均温度，并放入温度字段中，创建一个温度数据的图例，利用geoshow函数显示区域填色图。

1. #Read temporature data from a surfer grid data
   
2. f = addfile_surfer('D:/Temp/ascii/usgrid.dat')
   
3. tdata = f['var'][:,:]
   
4. 
   
5. #Read US states layer from shape file
   
6. us = shaperead('D:/Temp/map/states.shp')
   
7. 
   
8. #Add temp field
   
9. us.addfield('temp', 'float')
   
10. 
    
11. #Average temporature for each state and add to the temp field
    
12. for i in range(us.shapenum()):
    
13.     rpoly = us.shapes()
    
14.     mdata = tdata.maskout(rpoly)
    
15.     tave = mdata.ave()
    
16.     us.setcellvalue('temp', i, tave)
    
17. 
    
18. #Plot
    
19. axesm()
    
20. geoshow('country', facecolor='lightgray', edge=False)
    
21. levs = arange(270, 300, 2)
    
22. cols = makecolors(len(levs)+1)
    
23. ls = makesymbolspec('polygon', field='temp', levels=levs, colors=cols,
    
24.     edge=True, edgecolor='gray')
    
25. geoshow(us, symbolspec=ls)
    
26. xlim(-128, -65)
    
27. ylim(24, 50)
    
28. title('Average temporature of states')
    
29. colorbar(us)

#Read temporature data from a surfer grid data<br /> f = addfile_surfer('D:/Temp/ascii/usgrid.dat')<br /> tdata = f['var'][:,:]<br /> <br /> #Read US states layer from shape file<br /> us = shaperead('D:/Temp/map/states.shp')<br /> <br /> #Add temp field<br /> us.addfield('temp', 'float')<br /> <br /> #Average temporature for each state and add to the temp field<br /> for i in range(us.shapenum()):<br />     rpoly = us.shapes()<i><br />     mdata = tdata.maskout(rpoly)<br />     tave = mdata.ave()<br />     us.setcellvalue('temp', i, tave)<br /> <br /> #Plot<br /> axesm()<br /> geoshow('country', facecolor='lightgray', edge=False)<br /> levs = arange(270, 300, 2)<br /> cols = makecolors(len(levs)+1)<br /> ls = makesymbolspec('polygon', field='temp', levels=levs, colors=cols, <br />     edge=True, edgecolor='gray')<br /> geoshow(us, symbolspec=ls)<br /> xlim(-128, -65)<br /> ylim(24, 50)<br /> title('Average temporature of states')<br /> colorbar(us)

qxtlyf

王老师，可以用站点资料吗?即每个县一个数据带表平均值

CHch

新手小白求问大神，for i in range(China.shapenum()):添加这句话之后报错>>> run script...
  File "E:\data\paifangyuanqingdan\py\11.py", line 10
    for i in range(China.shapenum()):
                                    ^
SyntaxError: mismatched input '<EOF>' expecting INDENT
这是什么问题呢

吴小小

qxtlyf 发表于 2018-3-10 08:15
王老师，可以用站点资料吗?即每个县一个数据带表平均值

我也想知道