[gdal-dev] NetCDF: how to read coordinates from "lon" and "lat" variables?

andy aborruso at gmail.com
Tue Dec 3 10:16:37 PST 2024


Hi Michael,
and thank you.

On Tue, 3 Dec 2024 at 19:06, Michael Sumner <mdsumner at gmail.com> wrote:

> You can't use a geotransform unless the data are actually regular. There's
> no saying whether the lon lat arrays describe a regular grid redundantly
> without inspecting them (it's also very common for formats to deliver lon
> lat coords when the data is actually a regular grid in a different crs, but
> it's not easy to automatically determine that situation which is the entire
> entropy-increasing problem with the degenerate coordinate model of
> netcdf-alikes. There's hope though, see xarray PR 9543 for a crack of light
> in the gloom).
>

However, I did not understand one thing. Is my vrt wrong and should I not
use “geotransform,” and is there a more correct way to write it?
Or is this type of netcdf not well readable with gdal (at least until PR
9543 closes)?

Best regards


-- 
___________________

Andrea Borruso
website: https://medium.com/tantotanto
38° 7' 48" N, 13° 21' 9" E, EPSG:4326
___________________

"cercare e saper riconoscere chi e cosa,
 in mezzo all’inferno, non è inferno,
e farlo durare, e dargli spazio"

Italo Calvino
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20241203/4ee3e755/attachment-0001.htm>


More information about the gdal-dev mailing list