gdalwarp is your friend ! See http://gdal.org/gdalwarp.html
Or if I see well from the relationship between the pixel and georeferenced coordinatesif your image doesn't really need to be warped, you could also try gdal_translate -a_srs srs_def -a_ullr ulx uly lrx lry. Le Tuesday 16 September 2008 21:20:58 John Mitchell, vous avez écrit : > Hi, > > For a raster file header how do I transform the GCP Projection and GCP into > a format where instead the Coordinate System contains corner coordinates? > > *Below is an example of a portion of a header file containing GCP > Projection and GCP and below the dotted line is how I would like it > transformed.* > > Driver: GTiff/GeoTIFF > Size is 11601, 9027 > Coordinate System is `' > GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS > 84",6378137,298.2572235630016,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6 >326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EP >SG","4326"]] GCP[ 0]: Id=1, Info= > (264,8127) -> (-95.8337916667,28.533825,0) > GCP[ 1]: Id=2, Info= > (264,275) -> (-95.8337916667,29.1006305556,0) > GCP[ 2]: Id=3, Info= > (11245,275) -> (-94.9336861111,29.1006305556,0) > GCP[ 3]: Id=4, Info= > (11245,8127) -> (-94.9336861111,28.533825,0) > GCP[ 4]: Id=5, Info= > (10436,8127) -> (-94.9999972222,28.533825,0) > > * > ........................................................................... >... * > > Driver: GTiff/GeoTIFF > Size is 1800, 900 > 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 = (-99.000000000000000,27.000000000000000) > Pixel Size = (0.001111111111111,-0.001666666666667) > Metadata: > AREA_OR_POINT=Area > TIFFTAG_SOFTWARE=Adobe Photoshop CS Windows > TIFFTAG_DATETIME=2005:12:26 22:09:51 > TIFFTAG_XRESOLUTION=72 > TIFFTAG_YRESOLUTION=72 > TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) > Corner Coordinates: > Upper Left ( -99.0000000, 27.0000000) ( 99d 0'0.00"W, 27d 0'0.00"N) > Lower Left ( -99.0000000, 25.5000000) ( 99d 0'0.00"W, 25d30'0.00"N) > Upper Right ( -97.0000000, 27.0000000) ( 97d 0'0.00"W, 27d 0'0.00"N) > Lower Right ( -97.0000000, 25.5000000) ( 97d 0'0.00"W, 25d30'0.00"N) > Center ( -98.0000000, 26.2500000) ( 98d 0'0.00"W, 26d15'0.00"N) > Band 1 Block=1800x1 Type=Byte, ColorInterp=Red > Band 2 Block=1800x1 Type=Byte, ColorInterp=Green > Band 3 Block=1800x1 Type=Byte, ColorInterp=Blue _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
