[gdal-dev] ERROR 1: No inverse operation

[Even Rouault]

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
>>> gdal-dev at lists.osgeo.org
>>> 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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230308/039f3ebd/attachment-0001.htm>