如何用srai将Shapely LineString转为H3六边形?报错求助
解决Shapely LineString转H3序列时的AttributeError问题
核心问题:函数误用
你大概率搞反了srai库的函数用途:h3_to_geoseries是将H3索引序列转换为对应多边形的GeoSeries,而要把Shapely LineString转换成覆盖它的H3六边形索引,应该用geoseries_to_h3函数。误用函数会导致内部处理逻辑错误,触发MultiPolygon无exterior属性的报错。
正确实现代码
import geopandas as gpd from shapely.geometry import LineString import srai # 示例LineString对象 line = LineString([(0, 0), (1, 1)]) # 转为GeoSeries格式 geo_series = gpd.GeoSeries([line]) # 指定H3分辨率(可根据需求调整,范围0-15),转换为H3索引序列 h3_indices = srai.geoseries_to_h3(geo_series, resolution=10)
若确实需要处理MultiPolygon的情况
如果你的流程中生成了MultiPolygon对象,可先拆分它为单个Polygon再处理:
from shapely.geometry import MultiPolygon import geopandas as gpd import srai # 你的MultiPolygon对象 multi_poly = MultiPolygon([...]) # 转为GeoSeries并拆分为单个Polygon split_geo_series = gpd.GeoSeries([multi_poly]).explode(index_parts=True) # 转换为H3索引序列 h3_indices = srai.geoseries_to_h3(split_geo_series, resolution=10)
内容的提问来源于stack exchange,提问作者Simon
相关产品推荐
相关产品推荐

