<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi,<div><br></div><div>Assuming the following:</div><div><br></div><div><div>>>> src_ds = gdal.Open('test_3.ter')</div><div>>>> dst_ds = driver.CreateCopy('test_6.hf2', src_ds, 0)</div><div>>>> dst_ds.SetGeoTransform([582495,1,0.5,4512717,0.5,-1])</div><div>0</div><div>>>> srs = osr.SpatialReference()</div><div>>>> srs.SetUTM(18,1)</div><div>0</div><div>>>> srs.SetWellKnownGeogCS('WGS84')</div><div>0</div><div>>>> dst_ds.SetProjection( srs.ExportToWkt())</div><div>0</div><div>>>> dst_ds = None</div><div>>>> src_ds = None</div><div>>>> gdalinfo.main(['foo','test_6.hf2'])</div><div>Driver: HF2/HF2/HFZ heightfield raster</div><div>Files: test_6.hf2</div><div> test_6.hf2.aux.xml</div><div>Size is 4, 4</div><div>Coordinate System is:</div><div>PROJCS["UTM Zone 18, Northern Hemisphere",</div><div> GEOGCS["WGS 84",</div><div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div><div> AUTHORITY["EPSG","7030"]],</div><div> TOWGS84[0,0,0,0,0,0,0],</div><div> AUTHORITY["EPSG","6326"]],</div><div> PRIMEM["Greenwich",0,</div><div> AUTHORITY["EPSG","8901"]],</div><div> UNIT["degree",0.0174532925199433,</div><div> AUTHORITY["EPSG","9108"]],</div><div> AUTHORITY["EPSG","4326"]],</div><div> PROJECTION["Transverse_Mercator"],</div><div> PARAMETER["latitude_of_origin",0],</div><div> PARAMETER["central_meridian",-75],</div><div> PARAMETER["scale_factor",0.9996],</div><div> PARAMETER["false_easting",500000],</div><div> PARAMETER["false_northing",0],</div><div> UNIT["Meter",1]]</div><div>Origin = (0.000000000000000,0.000000000000000)</div><div>Pixel Size = (1.000000000000000,1.000000000000000)</div><div>Metadata:</div><div> VERTICAL_PRECISION=1.000000</div><div>Corner Coordinates:</div><div>Upper Left ( 0.0000000, 0.0000000) ( 79d29'19.48"W, 0d 0' 0.01"N)</div><div>Lower Left ( 0.0000000, 4.0000000) ( 79d29'19.48"W, 0d 0' 0.13"N)</div><div>Upper Right ( 4.0000000, 0.0000000) ( 79d29'19.35"W, 0d 0' 0.01"N)</div><div>Lower Right ( 4.0000000, 4.0000000) ( 79d29'19.35"W, 0d 0' 0.13"N)</div><div>Center ( 2.0000000, 2.0000000) ( 79d29'19.41"W, 0d 0' 0.06"N)</div><div>Band 1 Block=256x1 Type=Float32, ColorInterp=Undefined</div><div> Unit Type: m</div><div>0</div><div>>>> </div></div><div><br></div><div>and assuming I wanted to be nearly in New York rather than La Concordia Peru,</div><div>how to say 'how far north'?</div><div><br></div><div>Thanks in advance,</div><div><br></div><div>Chris</div> </body>
</html>