甲骨文SDO_CS。变换在从 RD 转换为 WGS84 时会转换坐标

Question

我有一个包含空间数据的Oracle数据库，其中坐标以RD格式存储（srid 28992）。要将它们放在地图上（带有开放街道地图的传单），我需要WGS84（srid 4326）格式的坐标。

我应该能够通过使用SDO_CS来获取该格式的坐标。变换。当我这样做时，所有坐标似乎都向东移动了大约6米。

为了验证表中存储的数据没有问题，我创建了一个查询，该查询采用RD坐标155000463000，并将其转换为WGS84。从该查询中，我得到：521551639130507 538735225721137

如果我对op做同样的变换http://www.gpscoordinaten.nl/converteer-gps-coordinaten.php我得到52.15517，5.38721（这是正确的坐标）

使用的查询是这样的：

SELECT
  t.X,
  t.Y
FROM
  (
    SELECT
      SDO_CS.TRANSFORM( MDSYS.SDO_GEOMETRY( 2001, 28992, MDSYS.SDO_POINT_TYPE(
      155000, 463000, NULL ), NULL, NULL ), 4326 ) AS geometry
    FROM
      dual
  )
  geo,
  TABLE( SDO_UTIL.GETVERTICES( geo.geometry ) ) t;

知道为什么Oracle中的转换会移动坐标吗？

Answer 1

对于那些带着类似问题登陆这里的人来说。

正如公认的答案所指出的，变换参数有各种值。要在Oracle中更改SRID 28992的转换参数，必须执行以下操作：

update mdsys.sdo_datums set
  shift_x = 565.2369,
  shift_y = 50.0087,
  shift_z = 465.658,
  rotate_x = -0.406857330322398,
  rotate_y = 0.350732676542563,
  rotate_z =  -1.8703473836068 ,
  scale_adjust = 4.0812
where datum_id = 6289;

commit;

begin SDO_CS.UPDATE_WKTS_FOR_EPSG_DATUM(6289);
end;

commit;

我花了一些时间才弄清楚Oracle缓存了这些值。为了使这些更改生效，您必须注销并启动一个新会话。

Answer 2

看看这个 http://epsg.io/28992

PROJCS["Amersfoort / RD New",
    GEOGCS["Amersfoort",
        DATUM["Amersfoort",
            SPHEROID["Bessel 1841",6377397.155,299.1528128,
                AUTHORITY["EPSG","7004"]],
            TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725],
            AUTHORITY["EPSG","6289"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4289"]],
    PROJECTION["Oblique_Stereographic"],
    PARAMETER["latitude_of_origin",52.15616055555555],
    PARAMETER["central_meridian",5.38763888888889],
    PARAMETER["scale_factor",0.9999079],
    PARAMETER["false_easting",155000],
    PARAMETER["false_northing",463000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","28992"]]

但是对于Oracle＜code＞来说，从cs_srs中选择wktext3d，其中srid=28992 你得到了什么

PROJCS[ "Amersfoort / RD New", 
    GEOGCS["Amersfoort", 
        DATUM["Amersfoort", 
        SPHEROID[ "Bessel 1841", 6377397.155, 299.1528128, 
            AUTHORITY["EPSG", "7004"]], 
        TOWGS84[593.16, 26.15, 478.54, -1.304398008226, -0.103297414969, -1.144501530423, 4.0775], 
            AUTHORITY["EPSG", "6289"]], 
        PRIMEM["Greenwich", 0.000000, 
            AUTHORITY["EPSG","8901"]], 
        UNIT["degree (supplier to define representation)", 0.0174532925199433, 
            AUTHORITY["EPSG", "9122"]], 
        AXIS["Lat", NORTH], 
        AXIS["Long", EAST], 
        AUTHORITY["EPSG", "4289"]], 
    PROJECTION ["Oblique Stereographic"], 
    PARAMETER ["Latitude_Of_Origin", 52.1561605555555556], 
    PARAMETER ["Central_Meridian", 5.3876388888888889], 
    PARAMETER ["Scale_Factor", 0.9999079], 
    PARAMETER ["False_Easting", 155000.0], 
    PARAMETER ["False_Northing", 463000.0], 
    UNIT["metre", 1.0, 
        AUTHORITY["EPSG", "9001"]], 
    AXIS["X", EAST], 
    AXIS["Y", NORTH], 
    AUTHORITY["EPSG", "28992"]]

您会看到 TOWGS84 参数不同。

也许您应该在https://gis.stackexchange.com/讨论这个问题，或者询问Oracle技术支持，看他们是否在参考资料中放入了错误的数据。

附加组件

在 http://www.epsg-registry.org/ 查看官方 EPSG 数据

当您搜索

• Tpye = 坐标转换 - 单个
• 阿拉 = 荷兰

你得到了一堆转变。我认为你必须做一些更深入的调查。