vector.py 文件源码

def features_on_raster(self, filtercodes=None, raster_extent=None):
        from shapely.geometry import MultiPolygon
        raster_extent = raster_extent if raster_extent else self.demDataExtent

        if filtercodes:
            features = [feat for feat in self.features if self.fieldName_Description in feat['fields'] and feat['fields'][self.fieldName_Description] in filtercodes]
        else:
            features = self.features
        for feat in features:
            if not MultiPolygon(raster_extent).contains(feat['geometry']):
                return False
        return True