quadrilateral_interpolation

This code provides a fast quadrilateral interpolation to perform satellite image projection from swath geometry to any geographical projection defined by a CRS.

The method and the code implementation are intent to be very fast but may not as accurate as other backward interpolation methods (bilinear, ...). It is suitable for coarse resolution imagery and was initially developed for OLCI. The calculation is based on: Interpolation using an arbitrary quadrilateral

The API follows the conventions of xESMF for the Resampler object and of pyresample for the grid definition. See pyresample for useful tools to create and manipulate grids.

This package was designed to work for satellite images opened with satpy but can be used as a standalone projection tool as long as the basic grid definition information are provided (lat, lon for the source and CRS, pixel size and offset for the destination).

Installation

Install with pip:

pip install quadrilateral_interpolation

Example

An example notebook is provided in the example/ directory. It requires to first download an OLCI image (~700Mb) in the native ESA format.

Contributing

Pull requests are welcome to make this code faster, more accurate, nicer, or more useful.