numpy简单应用---等值线(多边形)平滑

原创

MiaoGIS

修改于 2020-11-04 18:01:32

1.3K0

修改于 2020-11-04 18:01:32

文章被收录于专栏：Python in AI-IOT Python in AI-IOTPython in AI-IOT

风向条件下的气体扩散模型计算得到的数据，相当于等值线需要在地图上展示。如下图，可以看到，离中心最近的多边形毛糙锯齿严重。

使用滑动窗口平均方法可以平滑多边形，而且滑动窗口大小根据每一个多边形的点数量动态计算，效果如下。

对比如下图

代码如下：

import numpy as np
import matplotlib.pyplot as plt

for i,vvL in enumerate(v):
	w=i+3
	ids=list(range(len(vvL)))
	vvA=np.array(vvL)
	vvA=np.vstack((vvA[-w:],vvA,vvA[:w]))
	vvA=list(map(lambda x:vvA[x:x+1+w+w].mean(axis=0).tolist(),ids))
	vvA=np.array(vvA)
	vvA=np.vstack((vvA,vvA[0]))
	plt.plot(vvA[:,0],vvA[:,1])

原创声明：本文系作者授权腾讯云开发者社区发表，未经许可，不得转载。

如有侵权，请联系 cloudcommunity@tencent.com 删除。

地图