Hi, Assuming the following: >>> src_ds = gdal.Open('test_3.ter')>>> dst_ds = >>> driver.CreateCopy('test_6.hf2', src_ds, 0)>>> >>> dst_ds.SetGeoTransform([582495,1,0.5,4512717,0.5,-1])0>>> srs = >>> osr.SpatialReference()>>> srs.SetUTM(18,1)0>>> >>> srs.SetWellKnownGeogCS('WGS84')0>>> dst_ds.SetProjection( >>> srs.ExportToWkt())0>>> dst_ds = None>>> src_ds = None>>> >>> gdalinfo.main(['foo','test_6.hf2'])Driver: HF2/HF2/HFZ heightfield >>> rasterFiles: test_6.hf2 test_6.hf2.aux.xmlSize is 4, 4Coordinate >>> System is:PROJCS["UTM Zone 18, Northern Hemisphere", GEOGCS["WGS 84", >>> DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, >>> AUTHORITY["EPSG","7030"]], TOWGS84[0,0,0,0,0,0,0], >>> AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, >>> AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, >>> AUTHORITY["EPSG","9108"]], AUTHORITY["EPSG","4326"]], >>> PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], >>> PARAMETER["central_meridian",-75], PARAMETER["scale_factor",0.9996], >>> PARAMETER["false_easting",500000], PARAMETER["false_northing",0], >>> UNIT["Meter",1]]Origin = (0.000000000000000,0.000000000000000)Pixel Size = >>> (1.000000000000000,1.000000000000000)Metadata: >>> VERTICAL_PRECISION=1.000000Corner Coordinates:Upper Left ( 0.0000000, >>> 0.0000000) ( 79d29'19.48"W, 0d 0' 0.01"N)Lower Left ( 0.0000000, >>> 4.0000000) ( 79d29'19.48"W, 0d 0' 0.13"N)Upper Right ( 4.0000000, >>> 0.0000000) ( 79d29'19.35"W, 0d 0' 0.01"N)Lower Right ( 4.0000000, >>> 4.0000000) ( 79d29'19.35"W, 0d 0' 0.13"N)Center ( 2.0000000, >>> 2.0000000) ( 79d29'19.41"W, 0d 0' 0.06"N)Band 1 Block=256x1 Type=Float32, >>> ColorInterp=Undefined Unit Type: m0>>> and assuming I wanted to be nearly in New York rather than La Concordia Peru,how to say 'how far north'? Thanks in advance, Chris
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev