OSGeo / libgeotiff

Official repository of the libgeotiff project
178 stars 69 forks source link

Wrong epsg number in pcs #85

Open arturbatory89 opened 1 year ago

arturbatory89 commented 1 year ago

Wrong geotiff keyed information

I am using listgeo and geotifcp to edit keyed information in geotiff raster files. I have a problem with Poland_CS2000 zones 7 and 8, the assigned georeferences from the gtf file show mixed information for zone 7 and 8.

Example:

Geotiff_Information: Version: 1 Key_Revision: 1.0 Tagged_Information: ModelTiepointTag (2,3): 0 0 0
7572000 5544500 0
ModelPixelScaleTag (1,3): 0.05 0.05 0
End_Of_Tags. Keyed_Information: GTModelTypeGeoKey (Short,1): ModelTypeProjected GTRasterTypeGeoKey (Short,1): RasterPixelIsArea GTCitationGeoKey (Ascii,14): "GCS_ETRS_1989" GeographicTypeGeoKey (Short,1): GCS_EUREF89 GeogAngularUnitsGeoKey (Short,1): Angular_Degree GeogEllipsoidGeoKey (Short,1): Ellipse_GRS_1980 GeogSemiMajorAxisGeoKey (Double,1): 6378137
GeogSemiMinorAxisGeoKey (Double,1): 6356752.31
GeogInvFlatteningGeoKey (Double,1): 298.257222
ProjectedCSTypeGeoKey (Short,1): PCS_ETRS89_Poland_CS2000_zone_8 PCSCitationGeoKey (Ascii,28): "ETRS89_Poland_CS2000_zone_7" ProjLinearUnitsGeoKey (Short,1): Linear_Meter ProjNatOriginLongGeoKey (Double,1): 21
ProjNatOriginLatGeoKey (Double,1): 0
ProjFalseEastingGeoKey (Double,1): 7500000
ProjFalseNorthingGeoKey (Double,1): 0
ProjScaleAtNatOriginGeoKey (Double,1): 0.999923
End_Of_Keys. End_Of_Geotiff. PCS = 2178 (ETRS89 / Poland CS2000 zone 7) Projection = 18307 (Poland CS2000 zone 7) Projection Method: CT_TransverseMercator ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N) ProjNatOriginLongGeoKey: 21.000000 ( 21d 0' 0.00"E) ProjScaleAtNatOriginGeoKey: 0.999923 ProjFalseEastingGeoKey: 7500000.000000 m ProjFalseNorthingGeoKey: 0.000000 m GCS: 4258/ETRS89 Datum: 6258/European Terrestrial Reference System 1989 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 ( 7572000.000, 5544500.000) Lower Left ( 7572000.000, 5544000.000) Upper Right ( 7572800.000, 5544500.000) Lower Right ( 7572800.000, 5544000.000) Center ( 7572400.000, 5544250.000)

In addition for zone 8 value for ProjectedCSTypeGeoKey (Short,1): Unknown-2179 I need "PCS_ETRS89_Poland_CS2000_zone_8"

I looked into source files and found epsg_pcs.inc with values:

ValuePair(PCS_ETRS89_Poland_CS2000_zone_5,2176) ValuePair(PCS_ETRS89_Poland_CS2000_zone_6,2177) ValuePair(PCS_ETRS89_Poland_CS2000_zone_7,2177) ValuePair(PCS_ETRS89_Poland_CS2000_zone_8,2178)

Poland_CS2000_zone_7 has epsg: 2178, and Poland_CS2000_zone_8 2179. I've changed it and recompiled but still have this problem.

Could you figure out where is the problem?

Thank you.