f是Polygon,fblock可以是fblock或者LineString,判断对象是否和f这个Polygon是否有相交等关系:
pythonif f.crosses(fblock) or f.contains(fblock) or f.overlaps(fblock)
g1是GeoDataFrame,qq是GeoDataFrame,取出qq在g1里面的部分,截断的:
pythong1 = GeoDataFrame([{'geometry': fblock}])
res_intersection = geopandas.overlay(g1, qq, how='intersection')
res = list(res_intersection.loc[0, "geometry"].coords)
本文作者:Dong
本文链接:
版权声明:本博客所有文章除特别声明外,均采用 CC BY-NC。本作品采用《知识共享署名-非商业性使用 4.0 国际许可协议》进行许可。您可以在非商业用途下自由转载和修改,但必须注明出处并提供原作者链接。 许可协议。转载请注明出处!