<div dir="ltr"><div>Hi Ian,<br><br>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 compute the same distance elements you need for IDW. So Kriging will never become faster than IDW.<br><br>Hence, the only realistic way to use Kriging on large data sets is to base each estimate on local subsets - typically by working in a TIN model, just like in the linear interpolation case. And if your data set is dense, the difference between Kriging and linear interpolation will be miniscule. Kriging will, however, provide a variance estimate for the prediction, which may be realistic assuming you provide a realistic covariance function for the phenomenon under investigation.<br><br>So to make any difference wrt the linear interpolation, you will need to include a larger number of data, e.g. by including neighbouring triangles of the TIN, or using another data structure to subset your observations.<br><br>Some years ago I wrote an informal note about this subject. The title is “Data Assimilation for Updates of Digital Terrain Models” and it is available over at <a href="http://sdfe.dk/media/2916626/kms_technical_report_14.pdf">http://sdfe.dk/media/2916626/kms_technical_report_14.pdf</a><br><br>/Thomas<br><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-06-25 5:19 GMT+02:00 Ian Reese <span dir="ltr"><<a href="mailto:ireese@linz.govt.nz" target="_blank">ireese@linz.govt.nz</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hi Gdal,
<br>
<br>
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 of data we commonly work with. Linear interpolation is a reliable and fast method, however the triangulated results are difficult to work with visually.
<br>
<br>
For some visual comparisons and reference, here is a link to a folder of sample points(shp) for a test area, two examples of kriged data for that test region, and a BASH script to reproduce IDW and Linear interpolations of the same region.<br>
<br>
<a href="https://drive.google.com/open?id=1HKbK4wMmVN6JZSvhxBQDub1WXBLE710H" target="_blank">https://drive.google.com/open?<wbr>id=<wbr>1HKbK4wMmVN6JZSvhxBQDub1WXBLE7<wbr>10H</a><br>
<br>
Here is an example of the Kriging method used by ArcGIS: <br>
<a href="http://desktop.arcgis.com/en/arcmap/latest/extensions/geostatistical-analyst/kriging-in-geostatistical-analyst.htm" target="_blank">http://desktop.arcgis.com/en/<wbr>arcmap/latest/extensions/<wbr>geostatistical-analyst/<wbr>kriging-in-geostatistical-<wbr>analyst.htm</a><br>
<br>
Few graphs explaining best interpolation methods to use with given data.  Gridded vs. scattered.<br>
<a href="https://www.neonscience.org/spatial-interpolation-basics" target="_blank">https://www.neonscience.org/<wbr>spatial-interpolation-basics</a><br>
<br>
If it is possible, we are curious to know what type of funding would be needed.<br>
<br>
Cheers,<br>
<br>
Ian<br>
</div>
<br>
<hr>
<font size="2" face="Verdana" color="Black"><br>
This message contains information, which may be in confidence and may be subject to legal privilege. If you are not the intended recipient, you must not peruse, use, disseminate, distribute or copy this message. If you have received this message in error, please
 notify us immediately (Phone 0800 665 463 or <a href="mailto:info@linz.govt.nz" target="_blank">info@linz.govt.nz</a>) and destroy the original message. LINZ accepts no responsibility for changes to this email, or for any attachments, after its transmission from LINZ. Thank You.<br>
</font>
</div>

<br>______________________________<wbr>_________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/gdal-dev</a><br></blockquote></div><br></div>