Hi all, I have problem with geotiff files created with gdal_translate when I try to read them with the C# bindings.
So I tested it with the GDALInfo commands. GDALInfo.exe (The c-api) will output the correct projection information: C:\tmp\release-1500-x64-gdal-mapserver\bin>gdal\apps\GDALInfo.exe w:\Lmv\F \7h0a.tif Driver: GTiff/GeoTIFF Files: w:\Lmv\Foto\7h\7h0a.tif Size is 5500, 5500 Coordinate System is `' GCP Projection = PROJCS["RT90 2.5 gon V", GEOGCS["RT90", DATUM["Rikets_koordinatsystem_1990", SPHEROID["Bessel 1841",6377397.155,299.1528128000009, AUTHORITY["EPSG","7004"]], AUTHORITY["EPSG","6124"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4124"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",15.80827777777778], PARAMETER["scale_factor",1], PARAMETER["false_easting",1500000], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AUTHORITY["EPSG","3021"]] GCP[ 0]: Id=1, Info= (250,250) -> (1550000,6405000,0) GCP[ 1]: Id=2, Info= (5250,250) -> (1555000,6405000,0) GCP[ 2]: Id=3, Info= (250,5250) -> (1550000,6400000,0) GCP[ 3]: Id=4, Info= (5250,5250) -> (1555000,6400000,0) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 5500.0) Upper Right ( 5500.0, 0.0) Lower Right ( 5500.0, 5500.0) Center ( 2750.0, 2750.0) Band 1 Block=5500x1 Type=Byte, ColorInterp=Gray But the c-sharp version has problems with projection string as it only reports this, lacking the projection data: C:\tmp\release-1500-x64-gdal-mapserver\bin>gdal\csharp\GDALInfo.exe w:\Lmv\Foto\ 7h\7h0a.tif Raster dataset parameters: Projection: RasterCount: 1 RasterSize (5500,5500) Using driver GeoTIFF Image Structure Metadata: 0: INTERLEAVE=BAND Corner Coordinates: Upper Left (0, 0) Lower Left (0, 5500) Upper Right (5500, 0) Lower Right (5500, 5500) Center (2750, 2750) GCP Projection: GCP[0]: Id=1, Info= (250,250) -> (1550000,6405000,0) GCP[1]: Id=2, Info= (5250,250) -> (1555000,6405000,0) GCP[2]: Id=3, Info= (250,5250) -> (1550000,6400000,0) GCP[3]: Id=4, Info= (5250,5250) -> (1555000,6400000,0) GCP Equivalent geotransformation parameters: t[0] = 1549750 t[1] = 1 t[2] = 0 t[3] = 6405250 t[4] = 0 t[5] = -1 Band 1 : DataType: Byte ColorInterpretation: Gray Description: Size (5500,5500) BlockSize (5500,1) What can I do to get the correct string to give to SpatialReference? (Tried both GDAL 1.6 and 1.7) Regards Patrik Svensson
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev