利用圆图形maskout数据

MeteoInfo

对于格点数据，如果希望获取某个站点周边固定距离范围内的数据，可以以站点坐标为中心生成一个圆图形（Circle），然后利用geolib的maskout函数将圆图形外的格点数据maskout为NaN值，圆图形内部的格点数据保持原值，这样利用maskout后的格点数据进行平均就可以得到圆内部数据的平均值了（NaN数据不参与计算）。示例脚本如下：

1. fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')
   
2. f = addfile(fn)
   
3. ps = f['PS'][0,'10:60','60:140']
   
4. 
   
5. # maskout
   
6. circle = Circle([110,40], 10, facecolor=None)
   
7. ps = geolib.maskout(ps, [circle])
   
8. ps_mean = ps.mean()
   
9. print(ps_mean)
   
10. 
    
11. # plot test
    
12. ax = axesm()
    
13. geoshow('country', edgecolor=(0,0,255))
    
14. ps[ps==nan] = 0   # for plotting
    
15. layer = scatter(ps, cmap='WhBlGrYeRe', edgecolor='gray')
    
16. ax.add_patch(circle)
    
17. title('Maskout with patches')
    
18. yticks([10,30,50])
    
19. colorbar(layer, shrink=0.8)

fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')<br /> f = addfile(fn)<br /> ps = f['PS'][0,'10:60','60:140']<br /> <br /> # maskout<br /> circle = Circle([110,40], 10, facecolor=None)<br /> ps = geolib.maskout(ps, [circle])<br /> ps_mean = ps.mean()<br /> print(ps_mean)<br /> <br /> # plot test<br /> ax = axesm()<br /> geoshow('country', edgecolor=(0,0,255))<br /> ps[ps==nan] = 0   # for plotting<br /> layer = scatter(ps, cmap='WhBlGrYeRe', edgecolor='gray')<br /> ax.add_patch(circle)<br /> title('Maskout with patches')<br /> yticks([10,30,50])<br /> colorbar(layer, shrink=0.8)

tulalang

王老师好，类似功能可以用纯java代码实现吗

MeteoInfo

tulalang 发表于 2023-3-30 09:37
王老师好，类似功能可以用纯java代码实现吗

可以，MeteoInfo是开源的，你如果有兴趣可以看相关源代码。

tulalang

MeteoInfo 发表于 2023-3-30 10:47
可以，MeteoInfo是开源的，你如果有兴趣可以看相关源代码。

好的，谢谢王老师~