Inline...
On 08/03/2023 14:33, Even Rouault wrote:
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.
I can't see any problem here:
listgeo HiBar_21022022_A1.tif
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
202245.240786848 421915.720555977 0
ModelPixelScaleTag (1,3):
0.125 0.125 1
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,26): "PCS Name = Israel_TM_Grid"
GeographicTypeGeoKey (Short,1): Code-4141 (Israel 1993)
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
GeogEllipsoidGeoKey (Short,1): Ellipse_GRS_1980
GeogSemiMajorAxisGeoKey (Double,1): 6378137
GeogSemiMinorAxisGeoKey (Double,1): 6356752.31414036
GeogInvFlatteningGeoKey (Double,1): 298.257222101
GeogTOWGS84GeoKey (Double,7): -48 55 52
0 0 0
0
ProjectedCSTypeGeoKey (Short,1): Code-2039 (Israel 1993 / Israeli
TM Grid)
PCSCitationGeoKey (Ascii,15): "Israel_TM_Grid"
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.
PCS = 2039 (Israel 1993 / Israeli TM Grid)
Projection = 18204 (Israeli TM)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 31.734394 ( 31d44' 3.82"N)
ProjNatOriginLongGeoKey: 35.204517 ( 35d12'16.26"E)
ProjScaleAtNatOriginGeoKey: 1.000007
ProjFalseEastingGeoKey: 219529.584000 m
ProjFalseNorthingGeoKey: 626907.390000 m
GCS: 4141/Israel 1993
Datum: 6141/Israel 1993
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
TOWGS84: -48,55,52,0,0,0,0
Projection Linear Units: 9001/metre (1.000000m)
Corner Coordinates:
Upper Left ( 202245.241, 421915.721) ( 35d 1'32.11"E, 29d53'
7.03"N)
Lower Left ( 202245.241, 418790.721) ( 35d 1'32.29"E,
29d51'25.54"N)
Upper Right ( 205370.241, 421915.721) ( 35d 3'28.57"E, 29d53'
7.17"N)
Lower Right ( 205370.241, 418790.721) ( 35d 3'28.72"E,
29d51'25.68"N)
Center ( 203807.741, 420353.221) ( 35d 2'30.42"E, 29d52'16.36"N)
Then, after:
gdal_translate -a_srs EPSG:2039 HiBar_21022022_A1.tif
HiBar_21022022_A1b.tif
The ERROR goes away and I get:
gdalinfo HiBar_21022022_A1b.tif
Driver: GTiff/GeoTIFF
Files: HiBar_21022022_A1b.tif
Size is 25000, 25000
Coordinate System is:
PROJCRS["Israel 1993 / Israeli TM Grid",
BASEGEOGCRS["Israel 1993",
DATUM["Israel 1993",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4141]],
CONVERSION["Israeli TM",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",31.7343936111111,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",35.2045169444444,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1.0000067,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",219529.584,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",626907.39,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["Cadastre, engineering survey, topographic mapping (large
and medium scale)."],
AREA["Israel - onshore; Palestine Territory - onshore."],
BBOX[29.45,34.17,33.28,35.69]],
ID["EPSG",2039]]
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:
Upper Left ( 202245.241, 421915.721) ( 35d 1'32.11"E, 29d53' 7.03"N)
Lower Left ( 202245.241, 418790.721) ( 35d 1'32.29"E, 29d51'25.54"N)
Upper Right ( 205370.241, 421915.721) ( 35d 3'28.57"E, 29d53' 7.17"N)
Lower Right ( 205370.241, 418790.721) ( 35d 3'28.72"E, 29d51'25.68"N)
Center ( 203807.741, 420353.221) ( 35d 2'30.42"E, 29d52'16.36"N)
Band 1 Block=25000x1 Type=Byte, ColorInterp=Red
Band 2 Block=25000x1 Type=Byte, ColorInterp=Green
Band 3 Block=25000x1 Type=Byte, ColorInterp=Blue
Band 4 Block=25000x1 Type=Byte, ColorInterp=Undefined
and the listgeo output is also different:
listgeo HiBar_21022022_A1b.tif
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
202245.240786848 421915.720555977 0
ModelPixelScaleTag (1,3):
0.125 0.125 0
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,30): "Israel 1993 / Israeli TM Grid"
GeogCitationGeoKey (Ascii,12): "Israel 1993"
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
ProjectedCSTypeGeoKey (Short,1): Code-2039 (Israel 1993 / Israeli
TM Grid)
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.
PCS = 2039 (Israel 1993 / Israeli TM Grid)
Projection = 18204 (Israeli TM)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 31.734394 ( 31d44' 3.82"N)
ProjNatOriginLongGeoKey: 35.204517 ( 35d12'16.26"E)
ProjScaleAtNatOriginGeoKey: 1.000007
ProjFalseEastingGeoKey: 219529.584000 m
ProjFalseNorthingGeoKey: 626907.390000 m
GCS: 4141/Israel 1993
Datum: 6141/Israel 1993
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)
Corner Coordinates:
Upper Left ( 202245.241, 421915.721) ( 35d 1'32.11"E, 29d53'
7.03"N)
Lower Left ( 202245.241, 418790.721) ( 35d 1'32.29"E,
29d51'25.54"N)
Upper Right ( 205370.241, 421915.721) ( 35d 3'28.57"E, 29d53'
7.17"N)
Lower Right ( 205370.241, 418790.721) ( 35d 3'28.72"E,
29d51'25.68"N)
Center ( 203807.741, 420353.221) ( 35d 2'30.42"E, 29d52'16.36"N)
So I guess that solves the problem, allowing for reprojection, and
mapping together with other datasets.
Thanks for you help (spot on, as always!)
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.
--
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
+972-523-665918
