Javier,
The error message comes from the XYZ driver which is confused and
believes it might be for it, but this is definitely not a .xyz file
There are quite a lot of variations of ASCII grid formats. This one
reminds me a bit of the IGN France ones that are handled by the
https://github.com/OSGeo/gdal-extra-drivers/blob/master/src/ignfheightasciigrid/ignfheightasciigrid.cpp
driver (note that it has been retired as out of tree), but I can see
that this .gri file has min_lat, max_lat first whereas the French ones
have min_lon, max_lon :-)
I've tried with
https://www.isgeoid.polimi.it/Geoid/America/Argentina/public/GEOIDEAR16_20160419.isg,
linked from
https://www.isgeoid.polimi.it/Geoid/America/Argentina/argentina2016_g.html
. Hopefully it has the same content as the .gri file. I had to push in
https://github.com/OSGeo/gdal/pull/7730 a fix to open it as the rounding
of coordinates & resolution is quite rough in those .isg files.
with the fix:
$ gdalinfo GEOIDEAR16_20160419.isg -mm
Driver: ISG/International Service for the Geoid
Files: GEOIDEAR16_20160419.isg
Size is 1440, 2220
Coordinate System is:
GEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["geodetic latitude (Lat)",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["geodetic longitude (Lon)",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (-76.009853500000006,-20.008733500000002)
Pixel Size = (0.016667000000000,-0.016667000000000)
Metadata:
MODEL_NAME=GEOIDEAR16
MODEL_TYPE=gravimetric
Corner Coordinates:
Upper Left ( -76.0098535, -20.0087335) ( 76d 0'35.47"W, 20d 0'31.44"S)
Lower Left ( -76.0098535, -57.0094735) ( 76d 0'35.47"W, 57d 0'34.10"S)
Upper Right ( -52.0093735, -20.0087335) ( 52d 0'33.74"W, 20d 0'31.44"S)
Lower Right ( -52.0093735, -57.0094735) ( 52d 0'33.74"W, 57d 0'34.10"S)
Center ( -64.0096135, -38.5091035) ( 64d 0'34.61"W, 38d30'32.77"S)
Band 1 Block=1440x1 Type=Float32, ColorInterp=Undefined
Computed Min/Max=-6.007,45.203
NoData Value=-9999
Unit Type: meters
Note that CRS is hard coded to WGS 84 by the driver, and should be
altered to something more relevant.
Even
Le 09/05/2023 à 17:21, Javier Jimenez Shaw a écrit :
Hi
I want to read the Argentinian Geoid model with GDAL.
It is an ascii raster floating-point file with just one band.
https://dnsg.ign.gob.ar/apps/geoidear/documentos/GEOIDE-Ar16.gri (22 MB)
The content is starts with:
-57.0007400000 -20.0166670000 -76.0000000000 -52.0161870000
0.01666700000000001 0.01666700000000000
13.278 13.310 13.321 13.322 13.357 13.383 13.429 13.442
13.477 13.504 13.516 13.549 13.574 13.549 13.570 13.602
13.623 13.584 13.574 13.661 13.766 13.802 13.823 13.821
13.835 13.881 13.907 13.944 14.005 14.060 14.085 14.070
14.098 14.153 14.221 14.259 14.277 14.396 14.446 14.451
14.390 14.471 14.572 14.645 14.681 14.674 14.702 14.745
14.796 14.843 14.872 14.872 14.952 15.065 15.110 15.053
15.056 15.050 15.095 15.126 15.125 15.154 15.222 15.261
15.244 15.369 15.480 15.519 15.464 15.462 15.505 15.477
and finishes with
19.069 19.065 19.057 19.047 19.039 19.031 19.021 19.012
19.004 18.997 18.990 18.987 18.987 18.988 18.991 18.998
19.001 19.009 19.027 19.042 19.067 19.141 19.209 19.208
19.179 19.156 19.147 19.150 19.152 19.156 19.164 19.155
Total number of lines: 399602
It seems to have a width of 1440 and height of 2220.
gdalinfo says
$ gdalinfo GEOIDE-Ar16.gri
ERROR 1: Maximum number of characters allowed reached.
ERROR 1: Couldn't determine X spacing
gdalinfo failed - unable to open 'GEOIDE-Ar16.gri'.
I do not know which exact format it is. I can deduce the meaning of
the values (specially knowing that it covers Argentina ;), but I would
like to read it directly if possible.
Is that a know format for GDAL?
Thanks,
Javier.
.___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ... .... ._ .__
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
