geopandas属性错误:“MultiPolygon”对象没有属性“exterior”

dfuffjeb  于 2023-01-28  发布在  Go
关注(0)|答案(1)|浏览(448)

我有两个GeoDataFrame。一个是爱荷华州州的GeoDataFrame,另一个是北美未来72小时的预尝降雨的GeoDataFrame。我想创建覆盖爱荷华州的降雨预报的GeoDataFrame。但我得到一个错误。

state_rain = gpd.overlay(NA_rain,iowa,how='intersection')      

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-39-ba8264ed63c2> in <module>()
  3 #ws_usa[['WTRSHD_ID','QPF']].groupby('WTRSHD_ID').max().reset_index()
  4 #state_rain = sjoin(usa_r,usa,how='inner',op='intersects')
----> 5 state_rain = gpd.overlay(usa_r,joined_states,how='intersection')
  6 ws_state = gpd.overlay(ws,joined_states,how='intersection')
  7 #print ws_usa.loc[ws_usa.WTRSHD_ID == 'IA-04']['QPF']

C:\Anaconda2\lib\site-packages\geopandas\tools\overlay.pyc in overlay(df1, df2, how, use_sindex)
 95 
 96     # Collect the interior and exterior rings
---> 97     rings1 = _extract_rings(df1)
 98     rings2 = _extract_rings(df2)
 99     mls1 = MultiLineString(rings1)

C:\Anaconda2\lib\site-packages\geopandas\tools\overlay.pyc in     _extract_rings(df)
 50                 # geom from layer is not valid attempting fix by buffer 0"
 51                 geom = geom.buffer(0)
---> 52             rings.append(geom.exterior)
 53             rings.extend(geom.interiors)
 54 

AttributeError: 'MultiPolygon' object has no attribute 'exterior'

我检查了type == '多面体',但GeoDataFrame中均不包含任何类型。

print NA_rain[NA_rain.geometry.type == 'MulitPolygon']
print iowa[iowa.geometry.type == 'MultiPolygon']

Empty GeoDataFrame
Columns: [END_TIME, ID, ISSUE_TIME, PRODUCT, QPF, START_TIME, UNITS, VALID_TIME, geometry]
Index: []
Empty GeoDataFrame
Columns: [sid, AFFGEOID, ALAND, AWATER, GEOID, LSAD, NAME, STATEFP, STATENS, STUSPS, geometry]
Index: []

如果我执行以下操作,交集将工作。

NA_rain.geometry = NA_rain.geometry.map(lambda x: x.convex_hull)

我有两个问题:1.为什么在我的NA_rain GeoDataFrame中没有显示任何多重多边形,2.除了将每个多边形转换为船体(这会破坏多边形的详细轮廓)之外,您建议如何处理多重多边形问题。

w1jd8yoj

w1jd8yoj1#

我同意@jdmcbr。我怀疑NA_rain中至少有一个特性是多重多边形,由于您显示的条件拼写错误(MulitPolygon而不是MultiPolygon),因此未检测到该特性。
如果你的数据框有MultiPolygon,你可以把它们全部转换成Polygon,一个比较麻烦的方法是把list()函数传递给每个MultiPolygon,然后分解成多行:

geom = NA_rain.pop('geometry')
geom = geom.apply(lambda x: list(x) if isinstance(x, MultiPolygon) else x).explode())
NA_rain = NA_rain.join(geom, how='inner')

请注意,第3行中的连接会复制多重多边形的所有多边形的数据框的其他属性,包括要素标识符,您可能希望稍后根据任务进行更改。

相关问题