tangxiaodao / nettopologysuite

Automatically exported from code.google.com/p/nettopologysuite
0 stars 0 forks source link

GeometryTransform.TransformGeometry returns always EMPTY #113

Closed GoogleCodeExporter closed 8 years ago

GoogleCodeExporter commented 8 years ago
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

GoogleCodeExporter commented 8 years ago
you can post support questions to related forum, like:
https://groups.google.com/forum/#!forum/nettopologysuite

Original comment by diegogu...@gmail.com on 16 Apr 2012 at 6:20