polartoolkit.fetch.rema#
- rema(version='1km', region=None, spacing=None, registration=None, **kwargs)[source]#
Load the REMA surface elevation data from Howat et al.[1]. The data are in EPSG3031 and reference to the WGS84 ellipsoid. To convert the data to be geoid-referenced, subtract a geoid model, which you can get from fetch.geoid().
Choose between “1km” or “500m” resolutions with parameter version.
accessed from https://www.pgc.umn.edu/data/rema/
- Parameters:
version (str, optional,) – choose which resolution to fetch, either “1km” or “500m”, by default is “1km”
region (tuple[float, float, float, float], optional) – region to clip the loaded grid to, in format [xmin, xmax, ymin, ymax], by default doesn’t clip
spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3
registration (str, optional,) – choose between ‘g’ (gridline) or ‘p’ (pixel) registration types, by default is the original type of the grid
**kwargs (optional) – additional keyword arguments to pass to the resample_grid function
- Returns:
Returns a loaded, and optional clip/resampled grid of the REMA DEM.
- Return type:
References