<div dir="ltr"><div>Hi Felix,</div><div><br></div><div>I'm really not sure I understand what exactly you are trying to do. <br></div><div>But if you have not tried pyresample yet it might do what you want. See <a href="https://pyresample.readthedocs.io/en/latest/">https://pyresample.readthedocs.io/en/latest/</a></div><div><br></div><div>Christoph<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 3 Dec 2020 at 14:16, Felix <<a href="mailto:felix.divo@sailingteam.tu-darmstadt.de">felix.divo@sailingteam.tu-darmstadt.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello All,<br>
<br>
I'm here for some suggestions on how to sample from a raster band at <br>
_non-regular_ points. Non-regular in this context means that the points <br>
have arbitrary coordinates and do not lie on a regular lattice/mesh like <br>
all pairs of integers would do.<br>
<br>
In my application, there is an irregular set of points for which we <br>
require bathymetric data (i.e. topographical data/elevation). I started <br>
out by manually reading in the ETOPO1 dataset (for now, we might use <br>
Earth2014 later) and then interpolating around each point sequentially. <br>
As we have up to hundred-tousands of points and doing the computation <br>
sequentially in Python/numpy, this is currently running for days (well, <br>
we canceled it at that point).<br>
<br>
Ideally, I'd like to use something like the "gdal_grid" tool, where the <br>
interpolation algorithm can be configured and the heavy-lifting is done <br>
in C rather than in Python. But it only supports creating a new <br>
_regular_ grid. The inverse distance interpolation with cutoff radii <br>
would be ideal.<br>
<br>
Does such a tool/such functionality exist in in GDAL or was my search <br>
rightfully fruitless? I mainly searched in the raster programs, the <br>
(C++) API of GDALRasterBand and of course generally in the internet. <br>
Thank you for your suggestions, in advance!<br>
<br>
Cheers and Stay Safe<br>
Felix Divo<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">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/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="arial, sans-serif"><font color="#444444"><b>Christoph Paulik</b><font> </font><font>|</font><font> </font>Head of Infrastructure<br></font></font></div><span><p class="MsoNormal" style="color:rgb(34,34,34)"><font face="monospace"><font color="#3d85c6"><b>VanderSat.</b></font><font color="#3d85c6"><b> Unmatched water data. Satellite observed.</b><br></font></font></p><p class="MsoNormal" style="color:rgb(34,34,34)"><font face="monospace" color="#3d85c6">Wilhelminastraat 43a, 2011 VK, Haarlem, The Netherlands</font></p></span><p class="MsoNormal"><font face="monospace" color="#3d85c6">M: +31(0)618271928 <br></font><span><font face="monospace"><a href="http://www.vandersat.com/" target="_blank"><font color="#444444">vandersat.com</font></a><font color="#666666"> </font></font><span style="font-family:monospace"><font color="#444444">|</font></span><span style="color:rgb(102,102,102);font-family:monospace"> </span><font face="monospace"><a href="https://www.linkedin.com/company/vandersat" target="_blank"><font color="#444444">Linkedin</font></a><font style="color:rgb(102,102,102)"> </font></font><span style="font-family:monospace"><font color="#444444">|</font></span><span style="font-family:monospace"><font color="#444444"> </font><a href="https://twitter.com/VanderSatNL" target="_blank"><font color="#444444">Twitter</font></a> </span></span></p><p style="color:rgb(0,0,0);font-size:12px;font-family:Calibri,sans-serif;line-height:normal;margin:0in 0in 0.0001pt"><span style="color:rgb(34,34,34);font-family:arial,sans-serif"></span></p><div style="color:rgb(0,0,0);font-size:12px;font-family:Helvetica;line-height:normal;margin:0in 0in 0.0001pt">--------------------------------------------------------------------------------------</div></div></div></div></div>