MeteoInfoLab脚本示例：Point in polygon or not

MeteoInfo

脚本示例在地图上添加一个多边形和一个点，并判断点是否在多边形内。

1. axesm(tickfontsize=12)
   
2. lworld = shaperead('D:/Temp/map/country1.shp')
   
3. geoshow(lworld, edgecolor='k', facecolor='g')
   
4. lat = array([30, 0, 18, 48, 30])
   
5. lon = array([60, 70, 130, 120, 60])
   
6. poly = geoshow(lat, lon, displaytype='polygon', color=[150,230,230,230], edgecolor='r', size=2)
   
7. plon = 80
   
8. plat = 30
   
9. geoshow(plat, plon, size=14, color='r', marker='S')
   
10. isin = inpolygon(plon, plat, poly.getShape())
    
11. print isin
    
12. xlim(-180, 180)
    
13. ylim(-90, 90)
    
14. xticks(arange(-180, 181, 30))
    
15. yticks(arange(-90, 91, 30))
    
16. title('Point in polygon: ' + str(isin))

axesm(tickfontsize=12)<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld, edgecolor='k', facecolor='g')<br /> lat = array([30, 0, 18, 48, 30])<br /> lon = array([60, 70, 130, 120, 60])<br /> poly = geoshow(lat, lon, displaytype='polygon', color=[150,230,230,230], edgecolor='r', size=2)<br /> plon = 80<br /> plat = 30<br /> geoshow(plat, plon, size=14, color='r', marker='S')<br /> isin = inpolygon(plon, plat, poly.getShape())<br /> print isin<br /> xlim(-180, 180)<br /> ylim(-90, 90)<br /> xticks(arange(-180, 181, 30))<br /> yticks(arange(-90, 91, 30))<br /> title('Point in polygon: ' + str(isin))

mezhe

学习学习！！

limeiman

简直不要太好啊，感谢