根据dataFrame的结果，有效地更新地理框架中值的方法。Within方法

Question

我有两个大的GeodataFrame：

一个来自一个ShapeFile，每个多边形具有称为“ Asapp”的浮子值。

其次是带有3x3米的渔网网格的质心和圆柱“ Asapp”零。

我需要的是填充第二个基数的“ asapp”，其中质心位于第一个基数内。

下面的代码执行此操作，但嘲笑的速度低15个多边形（最小的Shapefile之一具有更多的THA 20000多边形）。

# fishnet_grid is a dict created by GDAL with a raster with 3m pixel size
cells_in_wsg = np.array([(self.__convert_geom_sirgas(geom, ogr_transform), int(fid), 0.0) for fid, geom in fishnet_grid.items()])

# transforming the grid raster (which are square polygons) in a GeoDataframe of point using the centroids of the cells
fishnet_base = gpd.GeoDataFrame({'geometry': cells_in_wsg[..., 0], 'id': cells_in_wsg[..., 1], 'asapp': cells_in_wsg[..., 2]})
fishnet = gpd.GeoDataFrame({'geometry': fishnet_base.centroid, 'id': fishnet_base['id'], 'asapp': fishnet_base['asapp']})

# as_applied_data is the polygons GeoDataFrame
# the code below takes a lot of time to complete
for as_applied in as_applied_data.iterrows():
    fishnet.loc[fishnet.within(as_applied[1]['geometry']), ['asapp']] += [as_applied[1]['asapp']]

还有另一种方法可以更好地做到这一点吗？

tys！

Answer 1

我解决了问题。

我阅读了有关使用geopandas.overlay（ htttps：/htttps：/ /geopandas.org/en/stable/docs/user_guide/set_operations.html ）与许多多边形一起工作，但问题是它仅适用于多边形，我有多边形和点。

因此，我的解决方案是从点创建非常小的多边形（2厘米侧正方形），然后使用覆盖层。

最终代码：

# fishnet is now a GeoDataFrame of little squares
fishnet = gpd.GeoDataFrame({'geometry': cells_in_wsg[..., 0], 'id': cells_in_wsg[..., 1]})

#intersection has only the little squares that intersects with all as_applied_data polygons and the value in those polygons
intersection = gpd.overlay(fishnet, as_applied_data, how='intersection')

# now this is as easy as to calculate the mean and put it back in the fishnet using the merge
values = fishnet.merge(intersection.groupby(['id'], as_index=False).mean())
#and values has the the little squares, the geom_id and the mean values of the intersections!

效果很好！