用OpenTripPlanner生成的空多边形用.is_empty指令生成错误。

Question

我有一个用OpenTripPlanner生成的等时多边形：

{"type":"FeatureCollection","features":[{"type":"Feature","geometry":{"type":"MultiPolygon","coordinates":[]},"properties":{"time":-47},"id":"fid--576b228b_15b66d32d71_-7cbd"}]}

这个多边形被翻译成一个形状匀称的对象，并有以下指令：

isochrone = shapely.geometry.asShape(isochroneJSON['features'][0]['geometry'])

在Spyder中是这样的：

{u'type': u'FeatureCollection', u'features': [{u'geometry': {u'type': u'MultiPolygon', u'coordinates': []}, u'type': u'Feature', u'properties': {u'time': -47}, u'id': u'fid--576b228b_15b66d32d71_-7a54'}]}

在我看来它真的像一个空多边形。我的问题是，我想将它排除在我的其余治疗之外，并检查它是否有效和/或空。和下列指示：

if not isochrone.is_empty:

使用.is_empty shapely指令生成一个错误：

return (self._geom is None) or bool(self.impl['is_empty'](self))
self.__geom__, n = self.factory(self.context)

我完全迷路了，因为the only similar question似乎没有我自己的问题。

Answer 1

空的几何图形有点棘手，在您的特殊情况下(MultiPolygons)，您可以使用shape而不是asShape部分修复它：

import json
from shapely.geometry import MultiPolygon, shape, mapping, asShape

Q = shape({'type': 'MultiPolygon', 'coordinates': []})
print(Q.is_empty)
print(Q.geom_type)
print(json.dumps(mapping(Q)))

这个输出(在空多多边形的情况下，geom_type实际上不等于MultiPolygon)：

True
GeometryCollection
{"type": "MultiPolygon", "coordinates": []}