2

我有两个 GeoDataFrame。一个是爱荷华州，另一个是预测未来 72 小时北美的降雨。我想创建一个覆盖爱荷华州的降雨预报的 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 == 'MultiPolygon'，但 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.为什么没有任何 MultiPolygons 出现在我的 NA_rain GeoDataFrame 中，以及 2.除了将每个 Polygon 变成一个 convex_hull 之外，这会破坏 Polygon 的详细轮廓，您建议如何处理 MultiPolygon 问题。

4

1 回答 1

0

我同意@jdmcbr。我怀疑其中的至少一个功能是 MultiPolygon，由于您显示的条件拼写错误（而不是） NA_rain，因此未检测到该功能。MulitPolygonMultiPolygon

如果您的数据框有 MultiPolygons，您可以将它们全部转换为多边形。一个肮脏的方法是将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 行中的连接复制了 MultiPolygon 的所有多边形的数据框的其他属性，包括特征标识符，您可能希望稍后更改，具体取决于您的任务。

于 2021-04-07T02:48:55.640 回答