dctools.data.coordinates.get_dataset_geometry

dctools.data.coordinates.get_dataset_geometry(ds, coord_sys, max_points=50000)

Robustly extract a geometry from a dataset, avoiding memory errors for huge point clouds.

If the number of points is too large, subsample before computing the geometry.

Parameters:

• ds (xarray.Dataset)

• coord_sys (CoordinateSystem)

• max_points (int)

Return type:

shapely.geometry.base.BaseGeometry