确定网格正方形中的点和多边形答案

【问题标题】：Determining what points and polygons are in a grid square确定网格正方形中的点和多边形
【发布时间】：2020-01-24 06:57:17
【问题描述】：

我正在使用 python 和 geojson 来执行此操作，我想指定一个点，该点将是正方形的中心，假设正方形是 1 英里乘 1 英里，我想列出在中找到的所有点和多边形正方形，包括大于正方形的多边形。

我有多个 geojson 文件，所以需要检查几次，这很好。我一直在玩下面的代码，它检查单元格中心是否靠近正方形的中心，但对于形状奇特的多边形会有问题。我真的很想知道广场上的所有物品/功能。

import json
from shapely.geometry import shape, Point
from shapely.geometry import asShape, mapping

point = Point(14.9783266342289, 16.87265432621112)
max_distance_from_center = 1

with open('cells.geojson') as f:
    js = json.load(f)

for feature in js['features']:
    polygon = asShape(feature['geometry'])
    distance = point.distance(polygon.centroid)
    # print(f'{distance} - {polygon.centroid}')
    if distance < max_distance_from_center:
        print (f'Found cells containing polygon:{feature}')

对于源数据，我使用的是从https://azgaar.github.io/Fantasy-Map-Generator/ 导出的地图，网格应该是 10 英里乘 10 英里。有关如何执行此操作的建议？

更新：

这是一个画得很差的图表。在网格正方形内，我想识别落在正方形边界内的所有标记和多边形，即使它们超出它的一侧。我想列出在网格中存在的所有功能。我用黄色突出显示了这些区域。

Poorly draw image

我查看了交叉点，它可能会这样做。今晚试试。

【问题讨论】：

为清楚起见，您能否更好地定义“发现于”？您是指完全由边界方块包含的特征吗？部分包含？包含中心？
你为什么不用intersects方法而不是计算距离？

标签： python gis geojson shapely

【解决方案1】：

你可以试试这个：

首先，创建网格。

from shapely.geometry import Point
from matplotlib.pyplot as plt

point = Point(0, -10)
square = point.buffer(0.5).envelope

fig, ax = plt.subplots(figsize=(5,5))
gpd.GeoSeries(square).plot(ax=ax)
gpd.GeoSeries(point).plot(ax=ax, color = "black",markersize=30)
plt.grid()
plt.show()

然后，

import geopandas as gpd
# get geodataframe from geojson file
geo_df = gpd.GeoDataFrame.from_file('cells.geojson')
geo_df['grid_yn'] = geo_df['geometry'].apply(lambda x : x.intersects(square))

【讨论】：