<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Hi all,<br><br>to convert some coordinates from Lambert2 Etendu to WGS84.<br>I wrote the small python script you can see below.<br><br>It works fine on linux (gdal-1.9.2-5) but it doesn't work on win32<br>(GDAL-1.10.1.win32-py2.7.msi and gdal-110-1500-core.msi,downloaded from <br>http://vbkto.dyndns.org/sdk/PackageList.aspx?file=release-1500-gdal-1-10-1-mapserver-6-4-0.zip):<br><br>C:\import>python test_coord.py<br>ERROR 6: No translation for an empty SRS to PROJ.4 format is known.<br>Traceback (most recent call last):<br> File "test_coord.py", line 56, in <module><br> point.Transform(coordTrans)<br> File "z:\Programmi\python27\lib\site-packages\osgeo\ogr.py", line 4615, in Transform<br> return _ogr.Geometry_Transform(self, *args)<br>TypeError: in method 'Geometry_Transform', argument 2 of type 'OSRCoordinateTransformationShadow *'<br><br>env vars are fine:<br><br>C:\import>set gd<br>GDAL_DATA=c:\programmi\gdal\gdal-data<br>GDAL_DRIVER_PATH=c:\programmi\gdal\gdalplugins<br><br>what can I do to let the script working on windows too?<br><br>==========<br><br>import sys, os<br><br>import ogr<br>import osr<br><br>def getCoord(a,b):<br> WGS = osr.SpatialReference()<br> WGS.ImportFromEPSG(4326)<br> IGNF = osr.SpatialReference()<br> IGNF.ImportFromProj4("+init=IGNF:LAMBE +wktext")<br> coordTrans = osr.CoordinateTransformation(IGNF,WGS)<br> point = ogr.Geometry(ogr.wkbPoint)<br> point.AddPoint(a,b)<br> point.Transform(coordTrans)<br> return point.GetX(), point.GetY()<br><br><br>if __name__ == '__main__':<br> print getCoord(602262.89743,2670591.7564)<br> print "expected Lat=51.028160, Lon=2.368641"<br> </div></body>
</html>