Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Ian Reese
Ian From: Even Rouault [even.roua...@spatialys.com] Sent: Tuesday, 26 June 2018 12:21 a.m. To: gdal-dev@lists.osgeo.org Cc: Ian Reese; Jeremy Palmer Subject: Re: [gdal-dev] Kriging Interpolation Ian, Looking at your script, for the IDW, I see you don't use the max_points, nor the radius1 a

Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Ian Reese
Cheer Joakim! That I am not familiar with and will have a look From: Joaquim Manuel Freire Luís [jl...@ualg.pt] Sent: Tuesday, 26 June 2018 10:58 a.m. To: Ian Reese; Even Rouault; gdal-dev@lists.osgeo.org Cc: Jeremy Palmer Subject: RE: [gdal-dev] Kriging

Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Joaquim Manuel Freire Luís
I have plenty of methods for working around speed issues and have |>worked with any number of interpolations methods in GRASS, SAGA and |>Python. You have also the minimum curvature in GMT. Not lightning fast but quicker than krigging and normally better looking, specially on the extrapolation

Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Ian Reese
uault [mailto:even.roua...@spatialys.com] Sent: Tuesday, 26 June 2018 12:21 a.m. To: gdal-dev@lists.osgeo.org Cc: Ian Reese; Jeremy Palmer Subject: Re: [gdal-dev] Kriging Interpolation Ian, Looking at your script, for the IDW, I see you don't use the max_points, nor the radius1 and radius2 to

Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Even Rouault
Ian, Looking at your script, for the IDW, I see you don't use the max_points, nor the radius1 and radius2 to limit the neighbourhoud into which the points are searched. For the linear search, I've profiled the run and found that an enormous amout of time was spent processing the points outside

Re: [gdal-dev] Kriging Interpolation

2018-06-25 Thread Thomas Knudsen
Hi Ian, If you find IDW on N points slow, you will find Kriging on the same number of points unbearably slow, since for Kriging you will have to solve a system of N linear equations (essentially inverting an NxN matrix) to compute the Kriging weights. To set up the matrix you will even need to com

Re: [gdal-dev] Kriging Interpolation

2018-06-24 Thread Chris Marsh
For what is worth, as an end-user of GDAL, I would be very interested in having such capacity in GDAL. Cheers Chris On 24 June 2018 at 21:19, Ian Reese wrote: > Hi Gdal, > > I'd like to know if it is possible to make Kriging Interpolation part of > gdal_grid. We work with a good deal of scatt

[gdal-dev] Kriging Interpolation

2018-06-24 Thread Ian Reese
Hi Gdal, I'd like to know if it is possible to make Kriging Interpolation part of gdal_grid. We work with a good deal of scattered point elevation data and neither IDW or Linear interpolations(smooth or not smooth) produce the desired visual results. Nor are they fast enough for the density