Ah, the issue here is that the GeoTIFF keys in the .tif partially define
the projection (or possibly the .tif.aux.xml contain this invalid WKT).
But METHOD["unnamed"] means that the projection method is unknown, so it
is logical that GDAL/PROJ can't compute geographic coordinates from
projected coordinates.
The output of "listgeo HiBar_21022022_A1.tif" might be useful to further
diagnose.
You likely need to gdal_translate -a_srs EPSG:2039 to do something
useful with that file, and report to the data producer that the CRS
encoding is broken
Le 08/03/2023 à 13:18, Micha Silver a écrit :
Hi Even:
Here's the full gdalinfo, with the reported CRS:
gdalinfo HiBar_21022022_A1.tif
Driver: GTiff/GeoTIFF
Files: HiBar_21022022_A1.tif
HiBar_21022022_A1.tif.aux.xml
Size is 25000, 25000
Coordinate System is:
PROJCRS["Israel 1993 / Israeli TM Grid",
BASEGEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]]],
CONVERSION["unnamed",
METHOD["unnamed"]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]
Data axis to CRS axis mapping: 1,2
Origin = (202245.240786847687559,421915.720555976789910)
Pixel Size = (0.125000000000000,-0.125000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
ERROR 1: No inverse operation
Upper Left ( 202245.241, 421915.721)
ERROR 1: No inverse operation
Lower Left ( 202245.241, 418790.721)
ERROR 1: No inverse operation
Upper Right ( 205370.241, 421915.721)
ERROR 1: No inverse operation
Lower Right ( 205370.241, 418790.721)
ERROR 1: No inverse operation
Center ( 203807.741, 420353.221)
and proj version:
proj
Rel. 9.1.1, December 1st, 2022
(a newly installed debian 12)
Thanks,
Micha
On 08/03/2023 13:47, Even Rouault wrote:
Micha,
which CRS is reported and which PROJ version do you use ? It is
presumably one projection method for which PROJ only implements the
forward conversion, that is from geographic to projected coordinates,
but not the reverse way
Even
Le 08/03/2023 à 09:37, Micha Silver a écrit :
I received some Geotiff files, and gdalinfo shows:
Corner Coordinates:
ERROR 1: No inverse operation
Upper Left ( 205370.241, 418790.721)
ERROR 1: No inverse operation
Lower Left ( 205370.241, 415665.721)
ERROR 1: No inverse operation
Upper Right ( 206198.866, 418790.721)
ERROR 1: No inverse operation
Lower Right ( 206198.866, 415665.721)
ERROR 1: No inverse operation
Center ( 205784.553, 417228.221)
What does that mean?
gdalinfo --version
GDAL 3.6.2, released 2023/01/02
Thanks
--
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
+972-523-665918
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
--
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
+972-523-665918
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
