make_ice_mask: Downscale an ice mask

In EvolEcolGroup/pastclim: Manipulate Time Series of Palaeoclimate Reconstructions

Downscale an ice mask

Description

Downscaling the ice mask presents some issues. The mask is a binary raster, so any standard downscaling approach will still look very blocky. We can smooth the contour by applying a Gaussian filter. How strong that filter should be is very much a matter of personal opinion, as we do not have any data to compare to. This function attempts to use as sensible default value, but it is worth exploring alternative values to find a good solution.

Usage

make_ice_mask(ice_mask_low_res, land_mask_high_res, d = c(0.5, 3))

Arguments

ice_mask_low_res	a terra::SpatRaster of the low resolution ice mask to downscale (e.g. as obtained with get_ice_mask())
land_mask_high_res	a terra::SpatRaster of the land masks at different times (e.g. as obtained from make_land_mask()). The ice mask will be cropped and matched for the resolution of this land mask.

Value

a terra::SpatRaster of the ice mask (1's), with the rest of the world (sea and land) as NA's

EvolEcolGroup/pastclim documentation built on Nov. 6, 2023, 5:11 a.m.