New gdal_grid interpolation algorithm
|Reported by:||mckelvym||Owned by:||warmerdam|
|Severity:||normal||Keywords:||IDW, Inverse Distance Weighting, Nearest Neighbor, ArcPy, ESRI|
I would like to contribute the attached patch for a new algorithm to be included in gdal_grid. This algorithm option, "invdistnn", is a variation on the existing inverse distance weighting algorithm with the following features:
- Use a quadtree to search for points only in the neighborhood (within radius) of each grid cell. This dramatically improves on the time to produce an interpolated grid since it is not examining the entire point space like "invdist"
- Retains up to "max_points" nearest cells instead of arbitrarily the first "max_points + 1" cells
This implementation will produce interpolated grids that match the output of ESRI's ArcPy? arcpy.gp.Idw_sa method.
The affected files are: