Source code for safe.storage.clipping
"""Raster clipping by polygons
"""
from safe.common.polygon import clip_grid_by_polygons
#from safe.common.polygon import clip_lines_by_polygon
# FIXME (Ole): Order should be reversed and this should move into
# interpolation module.
# Then retire this one
# I THINK WE CAN RETIRE THIS NOW (3/9/12)
[docs]def clip_raster_by_polygons(R, P):
"""Separate raster grid points by polygons
Args:
* R: Raster layer
* P: Polygon layer
Returns:
* L: List of point vectors and their associated grid values.
One item for each polygon
"""
res = clip_grid_by_polygons(R.get_data(),
R.get_geotransform(),
P.get_geometry(as_geometry_objects=True))
# Return
return res