判断n个兴趣点的经纬度坐标是否在地理经纬度栅格内部python代码
时间: 2024-03-31 08:31:49 浏览: 66
判断指定的经纬度坐标点是否落在指定的多边形区域内
5星 · 资源好评率100%
可以使用Python中的shapely库来实现该功能。具体实现过程如下:
1. 安装shapely库:在命令行中输入`pip install shapely`。
2. 导入必要的库:在Python代码中导入以下库。
```
from shapely.geometry import Point
from shapely.geometry.polygon import Polygon
```
3. 定义栅格边界坐标:将栅格边界的经纬度坐标定义为一个多边形。
```
coords = [(lon1, lat1), (lon2, lat2), (lon3, lat3), ...]
polygon = Polygon(coords)
```
4. 判断点是否在栅格内:将要判断的兴趣点的经纬度坐标定义为一个点对象,然后使用`within`函数判断该点是否在多边形内。
```
point = Point(lon, lat)
if point.within(polygon):
print("该点在栅格内")
else:
print("该点不在栅格内")
```
完整代码示例:
```
from shapely.geometry import Point
from shapely.geometry.polygon import Polygon
# 定义栅格边界坐标
coords = [(116.37363, 39.915606), (116.37363, 39.918311), (116.376434, 39.918311), (116.376434, 39.915606)]
polygon = Polygon(coords)
# 判断点是否在栅格内
point = Point(116.375, 39.916)
if point.within(polygon):
print("该点在栅格内")
else:
print("该点不在栅格内")
```
阅读全文