def start_stop_on_raster(self, raster_extent=None):
from shapely.geometry import Point, MultiPolygon
raster_extent = raster_extent if raster_extent else self.demDataExtent
if self.features and raster_extent:
for feat in self.features:
if not (MultiPolygon(raster_extent).contains(Point(feat['geometry'].coords[0])) and
MultiPolygon(raster_extent).contains(Point(feat['geometry'].coords[-1]))):
return False
return True
评论列表
文章目录