NetTopologySuite修改SRID后坐标未转换的问题求助

你遇到的问题核心是：直接修改SRID或通过GeometryFactory创建对象，仅能修改几何对象的空间参考标识，不会对坐标值执行投影转换。要实现4326到3857的坐标转换，需要使用NetTopologySuite配套的坐标投影转换工具。

原因说明 #

SRID只是几何对象的空间参考标识符，修改它不会自动转换坐标。GeometryFactory仅负责创建指定SRID的几何对象，但不会对传入的原始坐标做投影计算——你传入的是4326的坐标，它就直接用这个坐标创建SRID为3857的对象，坐标值自然不会变。

解决方案 #

需要借助CoordinateTransformationFactory（NetTopologySuite 2.5.0默认集成相关功能）来创建坐标系转换实例，对坐标执行真正的投影转换。

完整代码示例 #

using NetTopologySuite.Geometries;
using NetTopologySuite.CoordinateSystems;
using NetTopologySuite.CoordinateSystems.Transformations;

const double inputX = 1.57554173300002;
const double inputY = 42.646007512;

// 1. 创建原始4326坐标系的Point
var originalPoint = new Point(new Coordinate(inputX, inputY)) { SRID = 4326 };

// 2. 获取4326和3857的坐标系实例
var csFactory = new CoordinateSystemFactory();
var wgs84 = csFactory.CreateFromWkt("GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]]");
var webMercator = csFactory.CreateFromWkt("PROJCS[\"WGS 84 / Pseudo-Mercator\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"central_meridian\",0],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"3857\"]]");

// 3. 创建坐标转换实例
var ctFactory = new CoordinateTransformationFactory();
var transformation = ctFactory.CreateFromCoordinateSystems(wgs84, webMercator);

// 4. 转换原始坐标
var transformedCoord = transformation.MathTransform.Transform(new[] { inputX, inputY });

// 5. 创建SRID为3857的新Point
var targetFactory = new GeometryFactory(new PrecisionModel(), 3857);
var transformedPoint = targetFactory.CreatePoint(new Coordinate(transformedCoord[0], transformedCoord[1]));

// 输出结果
Console.WriteLine($"x = {transformedPoint.X}, y = {transformedPoint.Y}");
// 输出：x = 175388.50344111404, y = 5258244.862258436

针对Polygon/MultiPolygon的处理 #

对于面类几何对象，只需遍历其所有坐标点执行转换，再用目标SRID的GeometryFactory重新构建对象即可，核心转换逻辑和Point一致。示例代码：

// 假设已有SRID为4326的Polygon
var originalPolygon = ...;
var transformedCoords = originalPolygon.Coordinates.Select(c => 
    transformation.MathTransform.Transform(new[] { c.X, c.Y }));
var transformedPolygon = targetFactory.CreatePolygon(
    new LinearRing(transformedCoords.Select(c => new Coordinate(c[0], c[1])).ToArray())
);

注意事项 #

• 确保项目已引入NetTopologySuite和NetTopologySuite.CoordinateSystems相关包（NetTopologySuite 2.5.0通常会自动包含依赖）。
• 手动编写WKT是直接获取坐标系的方式，也可通过EPSG代码获取（需额外依赖ProjNet.CoordinateSystems.Epsg包）。

内容的提问来源于stack exchange，提问作者Alex Smith