<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Hi,</p>
<p>I guess this discussion, and past similar ones, calls for an
enhancement. A new API function, like CPLErr
GDALRasterInterpolateAtPoint(GDALRasterBandH, double dfPixel,
double dfLocation, GDALRIOResampleAlg eInterpolation, double
*pdfValue), that could be used by a new mode of gdallocationinfo.
The GDALRPCGetDEMHeight() function in alg/gdal_rpc.cpp is a
plausible candidate implementation for bilinear and bicubic (we
could potentially restrict to that at the moment).<br>
</p>
<p>Even<br>
</p>
<div class="moz-cite-prefix">Le 24/04/2024 à 10:33, Javier Jimenez
Shaw via gdal-dev a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CADRrdKvRc1eGhqxOhZAEPkeaVu7YKvdGk9J-XR1Y_QmSEHKbYA@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr">
<div>Hi</div>
<div><br>
</div>
<div>I would like to read in QGIS or GDAL an interpolated value
in a DSM (well, actually it is a geoid model, but it is the
same behaviour). See that I do not want the pixel value, but
the linear interpolation among the neighbour pixels, assuming
that the pixel value is in the center of the pixel.<br>
</div>
<div>For instance, this file</div>
<div><a
href="https://www.isgeoid.polimi.it/Geoid/Asia/Japan/japan2000_g.html"
moz-do-not-send="true" class="moz-txt-link-freetext">https://www.isgeoid.polimi.it/Geoid/Asia/Japan/japan2000_g.html</a></div>
<div><br>
</div>
<div>Is there any way to get it (without implementing the
interpolation myself)?</div>
<div>If I understood correctly gdallocationinfo is not
interpolating, just giving the pixel value.<br>
</div>
<div><br>
</div>
<div>Thanks<br>
</div>
<div><br>
</div>
<div>
<div>
<div dir="ltr" class="gmail_signature"
data-smartmail="gmail_signature">
<div dir="ltr">
<div>.___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ...
.... ._ .__</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>