Geotools 变换点失败

Question

我正在使用 java 开发一个应用程序，它允许将 EPSG:4326 点转换为 EPSG:2972。

我的代码有效：

//抓取两个坐标参考系之间的变换

MathTransform mathTransform = CRS.findMathTransform("EPSG:4326","EPSG:2972", true);

//要转换的点

DirectPosition2D srcDirectPosition2D = new DirectPosition2D();

srcDirectPosition2D.setCoordinateReferenceSystem(CoordSysINT);

srcDirectPosition2D.setLocation(4.4665424,-52.4648442);

DirectPosition2D destDirectPosition2D = new DirectPosition2D();

//转换

mathTransform.transform(srcDirectPosition2D, destDirectPosition2D);

//投影点

ProjectedPoint=新点(destDirectPosition2D.getX(), destDirectPosition2D.getY(),null);

结果是：X: 337473.6430296206 Y :493858.9919024287 但它是错误的。
如果我使用一个网站作为http://cs2cs.mygeodata.eu/ 并且我做同样的转换，真实的结果是 337470.842698;493860.962631

我的代码结果对应于 EPSG:4326 和 EPSG:32622 之间的转换，但我不明白为什么？

我的代码有什么问题？

非常感谢

问候

尼古拉斯

Answer 1

EPSG:2972 的 Proj.4 文本是

+proj=utm +zone=22 +ellps=GRS80 +towgs84=2,2,-2,0,0,0,0 +units=m +no_defs

在 Geotool 中，EPSG Database v8.6 将 EPSG:2972 的 WKT 显示为

2972=PROJCS["RGFG95 / UTM zone 22N", GEOGCS["RGFG95", DATUM["Reseau Geodesique Francais Guyane 1995", SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019 "]], TOWGS84[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], AUTHORITY["EPSG","6624"]], PRIMEM["Greenwich", 0.0, AUTHORITY[ "EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["大地经度", EAST], AXIS["大地纬度", NORTH], AUTHORITY["EPSG","4624"]] , PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], PARAMETER["central_meridian", -51.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 0.9996], PARAMETER[" false_easting", 500000.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","2972 "]]

注意突出显示的 toWGS84 转换是不同的。这将解释您所看到的差异。我不知道哪个是正确的，但我感觉 EPSG WKT 是错误的。

这种差异将 EPSG:2972 的行为更改为类似于 EPSG:32622 的行为。更专业的答案是 EPSG:2972 的基准面 (GRS80) 的行为就像 EPSG:32622 的基准面 (WGS84) 的椭球体