<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Gdal,<br>
<br>
I am trying to do a linear interpolation using gdal_grid and I am getting a segmentation fault with output grid cell sizes smaller than 8m.
<br>
<br>
<i>gdal_grid -zfield "depth" -a linear:nodata=-32768 -outsize 526 754 -of GTiff -l inshape inshape.shp outtif.tif --config GDAL_NUM_THREADS ALL_CPUS<br>
Grid data type is "Float64"<br>
Grid size = (526 754).<br>
Corner coordinates = (1671133.229795 5329097.839821)-(1675355.467899 5323053.902837).<br>
Grid cell size = (8.011837 8.005215).<br>
Source point count = 8363817.<br>
Algorithm name: "linear".<br>
Options are "radius=-1.000000:nodata=-32768.000000"<br>
<br>
0...10...20...30...40...50...60...70...80..Segmentation fault (core dumped)<br>
<br>
</i>If I change the number of input points (use a much smaller area) things seem to work out OK.<br>
<br>
<i>Grid cell size = (6.282705 6.842535).<br>
Source point count = 749329<br>
<br>
</i>So, I suspect there is a limitation to the number of points I can use for gdal_grid linear interpolation. Is there a way I can increase the number of input points? Or should I break up the region, interpolate smaller areas, and merge them back together
in the end. I have a large area with a large number of dense points (8 million points at 1m resolution). I am trying to avoid breaking it into as few smaller regions as possible to avoid artifacts at the edges when stitching it back together.
<br>
<br>
I'd really like to use the linear interpolation if possible since it has nice smooth results and is much faster than IDW with smoothing.
<br>
<br>
Cheers,<br>
<br>
Ian<br>
<i></i><i></i></div>
<br>
<hr>
<font face="Verdana" color="Black" size="2"><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 info@linz.govt.nz) 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>
</body>
</html>