@cpp.js/package-gdal 1.0.0-beta.20

package/dist/prebuilt/Emscripten-x86_64/include/gdal_alg.h

@@ -0,0 +1,730 @@
+/******************************************************************************
+ * $Id$
+ *
+ * Project:  GDAL Image Processing Algorithms
+ * Purpose:  Prototypes, and definitions for various GDAL based algorithms.
+ * Author:   Frank Warmerdam, warmerdam@pobox.com
+ *
+ ******************************************************************************
+ * Copyright (c) 2001, Frank Warmerdam
+ * Copyright (c) 2008-2012, Even Rouault <even dot rouault at spatialys.com>
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included
+ * in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
+ * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
+ * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ ****************************************************************************/
+
+#ifndef GDAL_ALG_H_INCLUDED
+#define GDAL_ALG_H_INCLUDED
+
+/**
+ * \file gdal_alg.h
+ *
+ * Public (C callable) GDAL algorithm entry points, and definitions.
+ */
+
+#ifndef DOXYGEN_SKIP
+#include "gdal.h"
+#include "cpl_minixml.h"
+#include "ogr_api.h"
+#include <stdint.h>
+#endif
+
+CPL_C_START
+
+int CPL_DLL CPL_STDCALL GDALComputeMedianCutPCT(
+    GDALRasterBandH hRed, GDALRasterBandH hGreen, GDALRasterBandH hBlue,
+    int (*pfnIncludePixel)(int, int, void *), int nColors,
+    GDALColorTableH hColorTable, GDALProgressFunc pfnProgress,
+    void *pProgressArg);
+
+int CPL_DLL CPL_STDCALL GDALDitherRGB2PCT(
+    GDALRasterBandH hRed, GDALRasterBandH hGreen, GDALRasterBandH hBlue,
+    GDALRasterBandH hTarget, GDALColorTableH hColorTable,
+    GDALProgressFunc pfnProgress, void *pProgressArg);
+
+int CPL_DLL CPL_STDCALL GDALChecksumImage(GDALRasterBandH hBand, int nXOff,
+                                          int nYOff, int nXSize, int nYSize);
+
+CPLErr CPL_DLL CPL_STDCALL GDALComputeProximity(GDALRasterBandH hSrcBand,
+                                                GDALRasterBandH hProximityBand,
+                                                char **papszOptions,
+                                                GDALProgressFunc pfnProgress,
+                                                void *pProgressArg);
+
+CPLErr CPL_DLL CPL_STDCALL GDALFillNodata(
+    GDALRasterBandH hTargetBand, GDALRasterBandH hMaskBand,
+    double dfMaxSearchDist, int bDeprecatedOption, int nSmoothingIterations,
+    char **papszOptions, GDALProgressFunc pfnProgress, void *pProgressArg);
+
+CPLErr CPL_DLL CPL_STDCALL GDALPolygonize(GDALRasterBandH hSrcBand,
+                                          GDALRasterBandH hMaskBand,
+                                          OGRLayerH hOutLayer, int iPixValField,
+                                          char **papszOptions,
+                                          GDALProgressFunc pfnProgress,
+                                          void *pProgressArg);
+
+CPLErr CPL_DLL CPL_STDCALL
+GDALFPolygonize(GDALRasterBandH hSrcBand, GDALRasterBandH hMaskBand,
+                OGRLayerH hOutLayer, int iPixValField, char **papszOptions,
+                GDALProgressFunc pfnProgress, void *pProgressArg);
+
+CPLErr CPL_DLL CPL_STDCALL GDALSieveFilter(
+    GDALRasterBandH hSrcBand, GDALRasterBandH hMaskBand,
+    GDALRasterBandH hDstBand, int nSizeThreshold, int nConnectedness,
+    char **papszOptions, GDALProgressFunc pfnProgress, void *pProgressArg);
+
+/*
+ * Warp Related.
+ */
+
+typedef int (*GDALTransformerFunc)(void *pTransformerArg, int bDstToSrc,
+                                   int nPointCount, double *x, double *y,
+                                   double *z, int *panSuccess);
+
+/*! @cond Doxygen_Suppress */
+#define GDAL_GTI2_SIGNATURE "GTI2"
+
+typedef struct
+{
+    GByte abySignature[4];
+    const char *pszClassName;
+    GDALTransformerFunc pfnTransform;
+    void (*pfnCleanup)(void *pTransformerArg);
+    CPLXMLNode *(*pfnSerialize)(void *pTransformerArg);
+    void *(*pfnCreateSimilar)(void *pTransformerArg, double dfSrcRatioX,
+                              double dfSrcRatioY);
+} GDALTransformerInfo;
+
+/*! @endcond */
+
+/*! @cond Doxygen_Suppress */
+void CPL_DLL GDALDestroyTransformer(void *pTransformerArg);
+int CPL_DLL GDALUseTransformer(void *pTransformerArg, int bDstToSrc,
+                               int nPointCount, double *x, double *y, double *z,
+                               int *panSuccess);
+void *GDALCreateSimilarTransformer(void *psTransformerArg, double dfSrcRatioX,
+                                   double dfSrcRatioY);
+/*! @endcond */
+
+/* High level transformer for going from image coordinates on one file
+   to image coordinates on another, potentially doing reprojection,
+   utilizing GCPs or using the geotransform. */
+
+void CPL_DLL *
+GDALCreateGenImgProjTransformer(GDALDatasetH hSrcDS, const char *pszSrcWKT,
+                                GDALDatasetH hDstDS, const char *pszDstWKT,
+                                int bGCPUseOK, double dfGCPErrorThreshold,
+                                int nOrder);
+void CPL_DLL *GDALCreateGenImgProjTransformer2(GDALDatasetH hSrcDS,
+                                               GDALDatasetH hDstDS,
+                                               char **papszOptions);
+void CPL_DLL *GDALCreateGenImgProjTransformer3(
+    const char *pszSrcWKT, const double *padfSrcGeoTransform,
+    const char *pszDstWKT, const double *padfDstGeoTransform);
+
+void CPL_DLL *GDALCreateGenImgProjTransformer4(
+    OGRSpatialReferenceH hSrcSRS, const double *padfSrcGeoTransform,
+    OGRSpatialReferenceH hDstSRS, const double *padfDstGeoTransform,
+    const char *const *papszOptions);
+
+void CPL_DLL GDALSetGenImgProjTransformerDstGeoTransform(void *,
+                                                         const double *);
+void CPL_DLL GDALDestroyGenImgProjTransformer(void *);
+int CPL_DLL GDALGenImgProjTransform(void *pTransformArg, int bDstToSrc,
+                                    int nPointCount, double *x, double *y,
+                                    double *z, int *panSuccess);
+
+void GDALSetTransformerDstGeoTransform(void *, const double *);
+void GDALGetTransformerDstGeoTransform(void *, double *);
+
+/* Geo to geo reprojection transformer. */
+void CPL_DLL *GDALCreateReprojectionTransformer(const char *pszSrcWKT,
+                                                const char *pszDstWKT);
+void CPL_DLL *
+GDALCreateReprojectionTransformerEx(OGRSpatialReferenceH hSrcSRS,
+                                    OGRSpatialReferenceH hDstSRS,
+                                    const char *const *papszOptions);
+void CPL_DLL GDALDestroyReprojectionTransformer(void *);
+int CPL_DLL GDALReprojectionTransform(void *pTransformArg, int bDstToSrc,
+                                      int nPointCount, double *x, double *y,
+                                      double *z, int *panSuccess);
+
+/* GCP based transformer ... forward is to georef coordinates */
+void CPL_DLL *GDALCreateGCPTransformer(int nGCPCount,
+                                       const GDAL_GCP *pasGCPList,
+                                       int nReqOrder, int bReversed);
+
+/* GCP based transformer with refinement of the GCPs ... forward is to georef
+ * coordinates */
+void CPL_DLL *GDALCreateGCPRefineTransformer(int nGCPCount,
+                                             const GDAL_GCP *pasGCPList,
+                                             int nReqOrder, int bReversed,
+                                             double tolerance, int minimumGcps);
+
+void CPL_DLL GDALDestroyGCPTransformer(void *pTransformArg);
+int CPL_DLL GDALGCPTransform(void *pTransformArg, int bDstToSrc,
+                             int nPointCount, double *x, double *y, double *z,
+                             int *panSuccess);
+
+/* Thin Plate Spine transformer ... forward is to georef coordinates */
+
+void CPL_DLL *GDALCreateTPSTransformer(int nGCPCount,
+                                       const GDAL_GCP *pasGCPList,
+                                       int bReversed);
+void CPL_DLL GDALDestroyTPSTransformer(void *pTransformArg);
+int CPL_DLL GDALTPSTransform(void *pTransformArg, int bDstToSrc,
+                             int nPointCount, double *x, double *y, double *z,
+                             int *panSuccess);
+
+/*! @cond Doxygen_Suppress */
+#ifdef GDAL_COMPILATION
+#define RPCInfoV1ToMD RPCInfoToMD
+#else
+#define RPCInfoToMD RPCInfoV2ToMD
+#endif
+char CPL_DLL **RPCInfoV1ToMD(GDALRPCInfoV1 *psRPCInfo);
+char CPL_DLL **RPCInfoV2ToMD(GDALRPCInfoV2 *psRPCInfo);
+/*! @endcond */
+
+/* RPC based transformer ... src is pixel/line/elev, dst is long/lat/elev */
+
+/*! @cond Doxygen_Suppress */
+#ifdef GDAL_COMPILATION
+#define GDALCreateRPCTransformerV1 GDALCreateRPCTransformer
+#else
+#define GDALCreateRPCTransformer GDALCreateRPCTransformerV2
+#endif
+
+void CPL_DLL *GDALCreateRPCTransformerV1(GDALRPCInfoV1 *psRPC, int bReversed,
+                                         double dfPixErrThreshold,
+                                         char **papszOptions);
+/*! @endcond */
+
+void CPL_DLL *GDALCreateRPCTransformerV2(const GDALRPCInfoV2 *psRPC,
+                                         int bReversed,
+                                         double dfPixErrThreshold,
+                                         char **papszOptions);
+
+void CPL_DLL GDALDestroyRPCTransformer(void *pTransformArg);
+int CPL_DLL GDALRPCTransform(void *pTransformArg, int bDstToSrc,
+                             int nPointCount, double *x, double *y, double *z,
+                             int *panSuccess);
+
+/* Geolocation transformer */
+
+void CPL_DLL *GDALCreateGeoLocTransformer(GDALDatasetH hBaseDS,
+                                          char **papszGeolocationInfo,
+                                          int bReversed);
+void CPL_DLL GDALDestroyGeoLocTransformer(void *pTransformArg);
+int CPL_DLL GDALGeoLocTransform(void *pTransformArg, int bDstToSrc,
+                                int nPointCount, double *x, double *y,
+                                double *z, int *panSuccess);
+
+/* Approximate transformer */
+void CPL_DLL *GDALCreateApproxTransformer(GDALTransformerFunc pfnRawTransformer,
+                                          void *pRawTransformerArg,
+                                          double dfMaxError);
+void CPL_DLL GDALApproxTransformerOwnsSubtransformer(void *pCBData,
+                                                     int bOwnFlag);
+void CPL_DLL GDALDestroyApproxTransformer(void *pApproxArg);
+int CPL_DLL GDALApproxTransform(void *pTransformArg, int bDstToSrc,
+                                int nPointCount, double *x, double *y,
+                                double *z, int *panSuccess);
+
+int CPL_DLL CPL_STDCALL GDALSimpleImageWarp(
+    GDALDatasetH hSrcDS, GDALDatasetH hDstDS, int nBandCount, int *panBandList,
+    GDALTransformerFunc pfnTransform, void *pTransformArg,
+    GDALProgressFunc pfnProgress, void *pProgressArg, char **papszWarpOptions);
+
+CPLErr CPL_DLL CPL_STDCALL
+GDALSuggestedWarpOutput(GDALDatasetH hSrcDS, GDALTransformerFunc pfnTransformer,
+                        void *pTransformArg, double *padfGeoTransformOut,
+                        int *pnPixels, int *pnLines);
+
+/** Flag for GDALSuggestedWarpOutput2() to ask to round-up output size */
+#define GDAL_SWO_ROUND_UP_SIZE 0x1
+/** Flag for GDALSuggestedWarpOutput2() to ask to force square pixels  */
+#define GDAL_SWO_FORCE_SQUARE_PIXEL 0x2
+
+CPLErr CPL_DLL CPL_STDCALL GDALSuggestedWarpOutput2(
+    GDALDatasetH hSrcDS, GDALTransformerFunc pfnTransformer,
+    void *pTransformArg, double *padfGeoTransformOut, int *pnPixels,
+    int *pnLines, double *padfExtent, int nOptions);
+
+/*! @cond Doxygen_Suppress */
+CPLXMLNode CPL_DLL *GDALSerializeTransformer(GDALTransformerFunc pfnFunc,
+                                             void *pTransformArg);
+CPLErr CPL_DLL GDALDeserializeTransformer(CPLXMLNode *psTree,
+                                          GDALTransformerFunc *ppfnFunc,
+                                          void **ppTransformArg);
+/*! @endcond */
+
+CPLErr CPL_DLL GDALTransformGeolocations(
+    GDALRasterBandH hXBand, GDALRasterBandH hYBand, GDALRasterBandH hZBand,
+    GDALTransformerFunc pfnTransformer, void *pTransformArg,
+    GDALProgressFunc pfnProgress, void *pProgressArg, char **papszOptions);
+
+/* -------------------------------------------------------------------- */
+/*      Contour Line Generation                                         */
+/* -------------------------------------------------------------------- */
+
+/** Contour writer callback type */
+typedef CPLErr (*GDALContourWriter)(double dfLevel, int nPoints, double *padfX,
+                                    double *padfY, void *);
+
+/** Contour generator opaque type */
+typedef void *GDALContourGeneratorH;
+
+GDALContourGeneratorH CPL_DLL
+GDAL_CG_Create(int nWidth, int nHeight, int bNoDataSet, double dfNoDataValue,
+               double dfContourInterval, double dfContourBase,
+               GDALContourWriter pfnWriter, void *pCBData);
+CPLErr CPL_DLL GDAL_CG_FeedLine(GDALContourGeneratorH hCG,
+                                double *padfScanline);
+void CPL_DLL GDAL_CG_Destroy(GDALContourGeneratorH hCG);
+
+/*! @cond Doxygen_Suppress */
+typedef struct
+{
+    void *hLayer;
+
+    double adfGeoTransform[6];
+
+    int nElevField;
+    int nElevFieldMin;
+    int nElevFieldMax;
+    int nIDField;
+    int nNextID;
+} OGRContourWriterInfo;
+
+CPLErr CPL_DLL OGRContourWriter(double, int, double *, double *, void *pInfo);
+/*! @endcond */
+
+CPLErr CPL_DLL GDALContourGenerate(
+    GDALRasterBandH hBand, double dfContourInterval, double dfContourBase,
+    int nFixedLevelCount, double *padfFixedLevels, int bUseNoData,
+    double dfNoDataValue, void *hLayer, int iIDField, int iElevField,
+    GDALProgressFunc pfnProgress, void *pProgressArg);
+
+CPLErr CPL_DLL GDALContourGenerateEx(GDALRasterBandH hBand, void *hLayer,
+                                     CSLConstList options,
+                                     GDALProgressFunc pfnProgress,
+                                     void *pProgressArg);
+
+/* -------------------------------------------------------------------- */
+/*      Viewshed Generation                                             */
+/* -------------------------------------------------------------------- */
+
+/** Viewshed Modes */
+typedef enum
+{
+    GVM_Diagonal = 1,
+    GVM_Edge = 2,
+    GVM_Max = 3,
+    GVM_Min = 4
+} GDALViewshedMode;
+
+/** Viewshed output types */
+typedef enum
+{
+    GVOT_NORMAL = 1,
+    GVOT_MIN_TARGET_HEIGHT_FROM_DEM = 2,
+    GVOT_MIN_TARGET_HEIGHT_FROM_GROUND = 3
+} GDALViewshedOutputType;
+
+GDALDatasetH CPL_DLL GDALViewshedGenerate(
+    GDALRasterBandH hBand, const char *pszDriverName,
+    const char *pszTargetRasterName, CSLConstList papszCreationOptions,
+    double dfObserverX, double dfObserverY, double dfObserverHeight,
+    double dfTargetHeight, double dfVisibleVal, double dfInvisibleVal,
+    double dfOutOfRangeVal, double dfNoDataVal, double dfCurvCoeff,
+    GDALViewshedMode eMode, double dfMaxDistance, GDALProgressFunc pfnProgress,
+    void *pProgressArg, GDALViewshedOutputType heightMode,
+    CSLConstList papszExtraOptions);
+
+bool CPL_DLL GDALIsLineOfSightVisible(
+    const GDALRasterBandH, const int xA, const int yA, const double zA,
+    const int xB, const int yB, const double zB, int *pnxTerrainIntersection,
+    int *pnyTerrainIntersection, CSLConstList papszOptions);
+
+/************************************************************************/
+/*      Rasterizer API - geometries burned into GDAL raster.            */
+/************************************************************************/
+
+CPLErr CPL_DLL GDALRasterizeGeometries(
+    GDALDatasetH hDS, int nBandCount, const int *panBandList, int nGeomCount,
+    const OGRGeometryH *pahGeometries, GDALTransformerFunc pfnTransformer,
+    void *pTransformArg, const double *padfGeomBurnValues,
+    CSLConstList papszOptions, GDALProgressFunc pfnProgress,
+    void *pProgressArg);
+
+CPLErr CPL_DLL GDALRasterizeGeometriesInt64(
+    GDALDatasetH hDS, int nBandCount, const int *panBandList, int nGeomCount,
+    const OGRGeometryH *pahGeometries, GDALTransformerFunc pfnTransformer,
+    void *pTransformArg, const int64_t *panGeomBurnValues,
+    CSLConstList papszOptions, GDALProgressFunc pfnProgress,
+    void *pProgressArg);
+
+CPLErr CPL_DLL GDALRasterizeLayers(
+    GDALDatasetH hDS, int nBandCount, int *panBandList, int nLayerCount,
+    OGRLayerH *pahLayers, GDALTransformerFunc pfnTransformer,
+    void *pTransformArg, double *padfLayerBurnValues, char **papszOptions,
+    GDALProgressFunc pfnProgress, void *pProgressArg);
+
+CPLErr CPL_DLL GDALRasterizeLayersBuf(
+    void *pData, int nBufXSize, int nBufYSize, GDALDataType eBufType,
+    int nPixelSpace, int nLineSpace, int nLayerCount, OGRLayerH *pahLayers,
+    const char *pszDstProjection, double *padfDstGeoTransform,
+    GDALTransformerFunc pfnTransformer, void *pTransformArg, double dfBurnValue,
+    char **papszOptions, GDALProgressFunc pfnProgress, void *pProgressArg);
+
+/************************************************************************/
+/*  Gridding interface.                                                 */
+/************************************************************************/
+
+/** Gridding Algorithms */
+typedef enum
+{
+    /*! Inverse distance to a power */ GGA_InverseDistanceToAPower = 1,
+    /*! Moving Average */ GGA_MovingAverage = 2,
+    /*! Nearest Neighbor */ GGA_NearestNeighbor = 3,
+    /*! Minimum Value (Data Metric) */ GGA_MetricMinimum = 4,
+    /*! Maximum Value (Data Metric) */ GGA_MetricMaximum = 5,
+    /*! Data Range (Data Metric) */ GGA_MetricRange = 6,
+    /*! Number of Points (Data Metric) */ GGA_MetricCount = 7,
+    /*! Average Distance (Data Metric) */ GGA_MetricAverageDistance = 8,
+    /*! Average Distance Between Data Points (Data Metric) */
+    GGA_MetricAverageDistancePts = 9,
+    /*! Linear interpolation (from Delaunay triangulation. Since GDAL 2.1 */
+    GGA_Linear = 10,
+    /*! Inverse distance to a power with nearest neighbor search for max points
+     */
+    GGA_InverseDistanceToAPowerNearestNeighbor = 11
+} GDALGridAlgorithm;
+
+/** Inverse distance to a power method control options */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridInverseDistanceToAPowerOptions) */
+    size_t nSizeOfStructure;
+    /*! Weighting power. */
+    double dfPower;
+    /*! Smoothing parameter. */
+    double dfSmoothing;
+    /*! Reserved for future use. */
+    double dfAnisotropyRatio;
+    /*! Reserved for future use. */
+    double dfAnisotropyAngle;
+    /*! The first radius (X axis if rotation angle is 0) of search ellipse. */
+    double dfRadius1;
+    /*! The second radius (Y axis if rotation angle is 0) of search ellipse. */
+    double dfRadius2;
+    /*! Angle of ellipse rotation in degrees.
+     *
+     * Ellipse rotated counter clockwise.
+     */
+    double dfAngle;
+    /*! Maximum number of data points to use.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPoints;
+    /*! Minimum number of data points to use.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPoints;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+} GDALGridInverseDistanceToAPowerOptions;
+
+/** Inverse distance to a power, with nearest neighbour search, control options
+ */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridInverseDistanceToAPowerNearestNeighborOptions) */
+    size_t nSizeOfStructure;
+    /*! Weighting power. */
+    double dfPower;
+    /*! The radius of search circle. */
+    double dfRadius;
+    /*! Smoothing parameter. */
+    double dfSmoothing;
+
+    /*! Maximum number of data points to use.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPoints;
+    /*! Minimum number of data points to use.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPoints;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+    /*! Maximum number of data points to use for each of the 4 quadrants.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPointsPerQuadrant;
+    /*! Minimum number of data points to use for each of the 4 quadrants.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPointsPerQuadrant;
+} GDALGridInverseDistanceToAPowerNearestNeighborOptions;
+
+/** Moving average method control options */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridMovingAverageOptions) */
+    size_t nSizeOfStructure;
+    /*! The first radius (X axis if rotation angle is 0) of search ellipse. */
+    double dfRadius1;
+    /*! The second radius (Y axis if rotation angle is 0) of search ellipse. */
+    double dfRadius2;
+    /*! Angle of ellipse rotation in degrees.
+     *
+     * Ellipse rotated counter clockwise.
+     */
+    double dfAngle;
+    /*! Maximum number of data points to use.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPoints;
+    /*! Minimum number of data points to average.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPoints;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+    /*! Maximum number of data points to use for each of the 4 quadrants.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPointsPerQuadrant;
+    /*! Minimum number of data points to use for each of the 4 quadrants.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPointsPerQuadrant;
+} GDALGridMovingAverageOptions;
+
+/** Nearest neighbor method control options */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridNearestNeighborOptions) */
+    size_t nSizeOfStructure;
+    /*! The first radius (X axis if rotation angle is 0) of search ellipse. */
+    double dfRadius1;
+    /*! The second radius (Y axis if rotation angle is 0) of search ellipse. */
+    double dfRadius2;
+    /*! Angle of ellipse rotation in degrees.
+     *
+     * Ellipse rotated counter clockwise.
+     */
+    double dfAngle;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+} GDALGridNearestNeighborOptions;
+
+/** Data metrics method control options */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridDataMetricsOptions) */
+    size_t nSizeOfStructure;
+    /*! The first radius (X axis if rotation angle is 0) of search ellipse. */
+    double dfRadius1;
+    /*! The second radius (Y axis if rotation angle is 0) of search ellipse. */
+    double dfRadius2;
+    /*! Angle of ellipse rotation in degrees.
+     *
+     * Ellipse rotated counter clockwise.
+     */
+    double dfAngle;
+    /*! Minimum number of data points to average.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPoints;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+    /*! Maximum number of data points to use for each of the 4 quadrants.
+     *
+     * Do not search for more points than this number.
+     */
+    GUInt32 nMaxPointsPerQuadrant;
+    /*! Minimum number of data points to use for each of the 4 quadrants.
+     *
+     * If less amount of points found the grid node considered empty and will
+     * be filled with NODATA marker.
+     */
+    GUInt32 nMinPointsPerQuadrant;
+} GDALGridDataMetricsOptions;
+
+/** Linear method control options */
+typedef struct
+{
+    /*! Added in GDAL 3.6 to detect potential ABI issues. Should be set to
+     * sizeof(GDALGridLinearOptions) */
+    size_t nSizeOfStructure;
+    /*! In case the point to be interpolated does not fit into a triangle of
+     * the Delaunay triangulation, use that maximum distance to search a nearest
+     * neighbour, or use nodata otherwise. If set to -1, the search distance is
+     * infinite. If set to 0, nodata value will be always used.
+     */
+    double dfRadius;
+    /*! No data marker to fill empty points. */
+    double dfNoDataValue;
+} GDALGridLinearOptions;
+
+CPLErr CPL_DLL GDALGridCreate(GDALGridAlgorithm, const void *, GUInt32,
+                              const double *, const double *, const double *,
+                              double, double, double, double, GUInt32, GUInt32,
+                              GDALDataType, void *, GDALProgressFunc, void *);
+
+/** Grid context opaque type */
+typedef struct GDALGridContext GDALGridContext;
+
+GDALGridContext CPL_DLL *
+GDALGridContextCreate(GDALGridAlgorithm eAlgorithm, const void *poOptions,
+                      GUInt32 nPoints, const double *padfX, const double *padfY,
+                      const double *padfZ, int bCallerWillKeepPointArraysAlive);
+
+void CPL_DLL GDALGridContextFree(GDALGridContext *psContext);
+
+CPLErr CPL_DLL GDALGridContextProcess(GDALGridContext *psContext, double dfXMin,
+                                      double dfXMax, double dfYMin,
+                                      double dfYMax, GUInt32 nXSize,
+                                      GUInt32 nYSize, GDALDataType eType,
+                                      void *pData, GDALProgressFunc pfnProgress,
+                                      void *pProgressArg);
+
+GDAL_GCP CPL_DLL *GDALComputeMatchingPoints(GDALDatasetH hFirstImage,
+                                            GDALDatasetH hSecondImage,
+                                            char **papszOptions,
+                                            int *pnGCPCount);
+
+/************************************************************************/
+/*  Delaunay triangulation interface.                                   */
+/************************************************************************/
+
+/** Triangle fact */
+typedef struct
+{
+    int anVertexIdx[3];   /**< index to the padfX/padfY arrays */
+    int anNeighborIdx[3]; /**< index to GDALDelaunayTriangulation.pasFacets, or
+                             -1 */
+    /* anNeighborIdx[k] is the triangle to the opposite side */
+    /* of the opposite segment of anVertexIdx[k] */
+} GDALTriFacet;
+
+/** Triangle barycentric coefficients.
+ *
+ * Conversion from cartesian (x,y) to barycentric (l1,l2,l3) with :
+ *  l1 = dfMul1X * (x - dfCxtX) + dfMul1Y * (y - dfCstY)
+ *  l2 = dfMul2X * (x - dfCxtX) + dfMul2Y * (y - dfCstY)
+ *  l3 = 1 - l1 - l2
+ */
+typedef struct
+{
+    double dfMul1X; /**< dfMul1X */
+    double dfMul1Y; /**< dfMul1Y */
+    double dfMul2X; /**< dfMul2X */
+    double dfMul2Y; /**< dfMul2Y */
+    double dfCstX;  /**< dfCstX */
+    double dfCstY;  /**< dfCstY */
+} GDALTriBarycentricCoefficients;
+
+/** Triangulation structure */
+typedef struct
+{
+    int nFacets;             /**< number of facets */
+    GDALTriFacet *pasFacets; /**< array of nFacets facets */
+    GDALTriBarycentricCoefficients
+        *pasFacetCoefficients; /**< arra of nFacets barycentric coefficients */
+} GDALTriangulation;
+
+int CPL_DLL GDALHasTriangulation(void);
+
+GDALTriangulation CPL_DLL *GDALTriangulationCreateDelaunay(int nPoints,
+                                                           const double *padfX,
+                                                           const double *padfY);
+int CPL_DLL GDALTriangulationComputeBarycentricCoefficients(
+    GDALTriangulation *psDT, const double *padfX, const double *padfY);
+int CPL_DLL GDALTriangulationComputeBarycentricCoordinates(
+    const GDALTriangulation *psDT, int nFacetIdx, double dfX, double dfY,
+    double *pdfL1, double *pdfL2, double *pdfL3);
+int CPL_DLL GDALTriangulationFindFacetBruteForce(const GDALTriangulation *psDT,
+                                                 double dfX, double dfY,
+                                                 int *panOutputFacetIdx);
+int CPL_DLL GDALTriangulationFindFacetDirected(const GDALTriangulation *psDT,
+                                               int nFacetIdx, double dfX,
+                                               double dfY,
+                                               int *panOutputFacetIdx);
+void CPL_DLL GDALTriangulationFree(GDALTriangulation *psDT);
+
+/*! @cond Doxygen_Suppress */
+#ifndef CPL_WARN_DEPRECATED_GDALOpenVerticalShiftGrid
+#define CPL_WARN_DEPRECATED_GDALOpenVerticalShiftGrid CPL_WARN_DEPRECATED
+#endif
+/*! @endcond */
+
+GDALDatasetH CPL_DLL GDALOpenVerticalShiftGrid(const char *pszProj4Geoidgrids,
+                                               int *pbError)
+    /*! @cond Doxygen_Suppress */
+    CPL_WARN_DEPRECATED_GDALOpenVerticalShiftGrid(
+        "GDALOpenVerticalShiftGrid() will be removed in GDAL 4.0")
+    /*! @endcond */
+    ;
+
+/*! @cond Doxygen_Suppress */
+#ifndef CPL_WARN_DEPRECATED_GDALApplyVerticalShiftGrid
+#define CPL_WARN_DEPRECATED_GDALApplyVerticalShiftGrid CPL_WARN_DEPRECATED
+#endif
+/*! @endcond */
+
+GDALDatasetH CPL_DLL GDALApplyVerticalShiftGrid(GDALDatasetH hSrcDataset,
+                                                GDALDatasetH hGridDataset,
+                                                int bInverse,
+                                                double dfSrcUnitToMeter,
+                                                double dfDstUnitToMeter,
+                                                const char *const *papszOptions)
+    /*! @cond Doxygen_Suppress */
+    CPL_WARN_DEPRECATED_GDALApplyVerticalShiftGrid(
+        "GDALApplyVerticalShiftGrid() will be removed in GDAL 4.0")
+    /*! @endcond */
+    ;
+
+CPL_C_END
+
+#endif /* ndef GDAL_ALG_H_INCLUDED */