1. I'm trying to convert Lambert 93 shapes to WGS 84
2. I use Shape files come from French IGN site :
http://professionnels.ign.fr/DISPLAY/000/528/175/5281750/GEOFLADept_FR_Corse_AV_
L93.zip
3. I build an ICoordinateTransformation
4. I call GeometryTransform.TransformGeometry to the IGeometry
4. The geometry I get is "GEOMETRY EMPTY"
Here is the code
string v_sourceCoordSys = "PROJCS[\"RGF93 /
Lambert-93\",GEOGCS[\"RGF93\",DATUM[\"Reseau_Geodesique_Francais_1993\",SPHEROID
[\"GRS
1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0,0,0,0,0,0,0
],AUTHORITY[\"EPSG\",\"6171\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901
\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY
[\"EPSG\",\"4171\"]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"sta
ndard_parallel_1\",49],PARAMETER[\"standard_parallel_2\",44],PARAMETER[\"latitud
e_of_origin\",46.5],PARAMETER[\"central_meridian\",3],PARAMETER[\"false_easting\
",700000],PARAMETER[\"false_northing\",6600000],UNIT[\"metre\",1,AUTHORITY[\"EPS
G\",\"9001\"]],AUTHORITY[\"EPSG\",\"2154\"]]";
string targetCoordSys = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS
84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"63
26\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01
745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]]";
ICoordinateSystem csSource =
ProjNet.Converters.WellKnownText.CoordinateSystemWktReader.Parse(v_sourceCoordSy
s) as ICoordinateSystem;
ICoordinateSystem csTarget =
ProjNet.Converters.WellKnownText.CoordinateSystemWktReader.Parse(targetCoordSys)
as ICoordinateSystem;
ICoordinateTransformation trans = new
CoordinateTransformationFactory().CreateFromCoordinateSystems(csSource,
csTarget);
using (ShapefileDataReader shapeDataReader = new
ShapefileDataReader(_shapeFile, GeometryFactory.Default))
{
while (shapeDataReader.Read())
{
IGeometry geomOut = GeometryTransform.TransformGeometry(GeometryFactory.Default, shapeDataReader.Geometry, trans.MathTransform);
// returns GEOMETRY EMPTY
}
}
I have the latest NTS V1.X (today april 14th)
Thank you
Original issue reported on code.google.com by xavier.f...@gmail.com on 14 Apr 2012 at 9:50
Original issue reported on code.google.com by
xavier.f...@gmail.com
on 14 Apr 2012 at 9:50