Hi list, I'm trying to convert the attached virtual file into a regular geo-tiff.
Note that the dataset defines both a geo-transform and GCPs but they have different projections. After conversion if I call gdalinfo on the result file I get $ gdalinfo out.tif Driver: GTiff/GeoTIFF Files:out.tif Size is 22000, 21890 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.2572235630016, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (359183.000000000000000,4696272.000000000000000) Pixel Size = (0.500000000000000,-0.500000000000000) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 359183.000, 4696272.000) (359183d 0'0.00"E,4696272d 0'17179869184.00"N) Lower Left ( 359183.000, 4685327.000) (359183d 0'0.00"E,4685327d 0'17179869184.00"N) Upper Right ( 370183.000, 4696272.000) (370183d 0'0.00"E,4696272d 0'17179869184.00"N) Lower Right ( 370183.000, 4685327.000) (370183d 0'0.00"E,4685327d 0'17179869184.00"N) Center ( 364683.000, 4690799.500) (364683d 0'0.00"E,4690799d30'17179869184.00"N) Band 1 Block=22000x1 Type=UInt16, ColorInterp=Gray NoData Value=0 The generated geo-tiff has wrong projection info. Is this a but? Is there some rule to determine which geographic information is authoritative? Best regards -- Antonio Valentino
CSKS1_GEC_B_S2_33_HH_RA_SF_20090507052430_20090507052438_S01_SBI_geo.vrt
Description: Binary data
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev