[gdal-dev] Problem with OSR projection on "+init=IGNF:LAMBE +wktext"

[Jean Michel PIERRET]

Hi,

 

                I'm using OSR (on GDAL 1.6.1) bindings in C# to project
points, and I have problem with the projection "+init=IGNF:LAMBE +wktext" to
"+init=IGNF:LAMB93"

At the first execution, I have one result and on all the others I have
another results. 

 

An exemple :

 

 

                SpatialReference _srin = null;

                SpatialReference _srout = null;

                CoordinateTransformation _ct = null;

 

                double[] x, y;

 

                _srin = new SpatialReference("");

                _srin.ImportFromProj4("+init=IGNF:LAMBE +wktext");

 

                _srout = new SpatialReference("");

                _srout.ImportFromProj4("+init=IGNF:LAMB93");

 

                _ct = new CoordinateTransformation(_srin, _srout);

 

                x = new double[] { 540000 };

                y = new double[] { 1820000 };

 

                _ct.TransformPoints(1, x, y, null);

 

                Console.WriteLine(x[0] + " - " + y[0]); 

                // Return : 586336,591390395 - 6254404,95793043

 

                x = new double[] { 540000 };

                y = new double[] { 1820000 };

 

                _ct.TransformPoints(1, x, y, null);

 

                Console.WriteLine(x[0] + " - " + y[0]);

                //Return : 586281,605619939 - 6254404,35662961

 

                _ct.Dispose();

                _srin.Dispose();

                _srout.Dispose();

 

Thanks for your Help.

 

Jean michel Pierret

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20091103/faf45227/attachment.html