Estoy tratando de implementar la misma solución que la función SpatialJoin_analysis
de ArcPy en GeoPandas.
Tengo un archivo shapefile que contiene polígonos y multipolígonos. Necesito realizar la siguiente operación espacial de ArcPy utilizando la librería de código abierto GeoPandas
.
arcpy.SpatialJoin_analysis(
target_features=spatial_join_in,
join_features=spatial_join_in,
out_feature_class=spatial_join_out,
join_operation="JOIN_ONE_TO_ONE",
join_type="KEEP_ALL",
field_mapping="#",
match_option="ARE_IDENTICAL_TO"
)
El mismo archivo shapefile se pasa como ambos archivos spatial_join_in.
Esta operación también agrega Join Count
en cada fila, que luego se utiliza para filtrar las áreas de densidad.
¿Cómo podemos lograr una solución similar utilizando GeoPandas?