[gdal-dev] transform point problem using gdal 1.5 with osr
prairie Last
prairie99 at gmail.com
Wed Dec 17 02:23:10 EST 2008
I am trying to tranform one coordinate system to another (wgs84 to utm).
Found one example on the internet (adapted a little). It got an error when I
try to run the program. Since I am new to GDAL, I cannot figure out how to
correct this. Please see the following program. I also tried to add one more
value to "res.ct.TransformPoint(lon, lat)", but it did not solve the
problem. Any help or hint appreciated. I also attach the error to the end.
from osgeo import osr
import math
import pyproj
def UTMZone(lon, lat):
return int((lon + 180) / 6) + 1
if __name__ == '__main__':
lat = 41.395304
lon = 73.310067
wgs84 = osr.SpatialReference()
utm = osr.SpatialReference()
wgs84.ImportFromProj4("+proj=latlong +datum=WGS84")
utm.ImportFromProj4("+proj=utm +zone="+str(UTMZone(lon,lat))+"
+datum=WGS84")
ct = osr.CoordinateTransformation(wgs84,utm)
res= ct.TransformPoint(lon, lat)
print res
Error Message:
Traceback (most recent call last):
File "C:\Python25\trans.py", line 28, in <module>
res= ct.TransformPoint(lon, lat)
File "C:\Python25\lib\site-packages\osgeo\osr.py", line 602, in
TransformPoint
return _osr.CoordinateTransformation_TransformPoint(*args)
NotImplementedError: Wrong number of arguments for overloaded function
'CoordinateTransformation_TransformPoint'.
Possible C/C++ prototypes are:
TransformPoint(double [3])
TransformPoint(double [3],double,double,double)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20081217/cf5859be/attachment.html
More information about the gdal-dev
mailing list