igapfill-package: Interface for filling missing values of Earth Observation...

In igapfill: Interface for Gap-Filling Earth Observation Datasets

Interface for filling missing values of Earth Observation Datasets

Description

By making extensive use of parallel computing, the functions of this package facilitate the application of the spatio-temporal gap-filling method Gapfill to time series of satellite images (TSSI).

Details

Only GeoTiff files are allowed. In passing, some of the functions of this package can be construed as independent builders of arguments used by Gapfill functions.

Datasets

Spatial subsets of the MOD13Q1 v061 NDVI product, in .tif format, covering Cerro de Garnica National Park (https://simec.conanp.gob.mx/ficha.php?anp=66&reg=11) located at Michoacan, Mexico. These subsets were collected from February 16, 2000 to December 16, 2024. The spatial and temporal resolutions of these images are 250m and 16 days, respectively.

garnica_250m_16_days_NDVI.tif	572 layers of NDVI
garnica_250m_16_days_pixel_reliability.tif	Pixel reliability layers corresponding to
	garnica_250m_16_days_NDVI.tif

Quality assessment summary

The following functions allow to compute the amount of missing values in a TSSI and to determine a sub-set of images to which apply the workflow of this package.

mvSieve	Computes amount of missing values in a TSSI
minmaxBlock	Determines sub-set of images with minimal (or maximal) missing values

Workflow

The following functions allow to define the required directory/folders structure employed by this package. Some of these functions are also useful for better data handling.

create_dirs	Sets up directory tree
dimsReport	Summary of dimensions of images to process
sort_split	Split large TSSI into smaller spatio-temporal chunks
waysToSplit	Briefing of ways to divide rows and columns of images

Interface

These are the functions of this package that allow for filling missing values of spatio-temporal subsets of satellite images using Gapfill.

applyGapfill	Parallel computing-based application of Gapfill
parallel_mosaic	Parallel rasterization and mosaicking (when required) of output of applyGapfill
igapfill	Console-based wrap-up of applyGapfill and parallel_mosaic

Miscellaneous

These functions can be used to obtain some of the arguments required by applyGapfill. In addition to this, these functions can also be employed to define some arguments of Gapfill.

get_3Darray	Assambles 3D array
get_4Darray	Assambles 4D array
get_LAT	Gets latitude information of RasterStack
get_LON	Gets longitude information of RasterStack

Author(s)

Tecuapetla-Gómez, I. [e-mail]

References

Gerber, F., de Jong, R., Schaepman, M.E., Schaepman-Strub, G., Furrer, R. (2018). Predicting missing values in spatio-temporal remote sensing data, IEEE Transactions on Geoscience and Remote Sensing, 1–13.

igapfill documentation built on June 8, 2025, 10:41 a.m.