Patrik, The GDALInfo csharp sample application doesn't exactly follow the C version by producing the same output. However it's pretty straight forward to report the projection similarly using the C# API. I've just added the required changes in the development version:
http://trac.osgeo.org/gdal/changeset/18705 Best regards, Tamas 2010/2/2 Patrik Svensson <pa...@bredband.net> > 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 >
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
