主页 / 问题 / 1536130
Accepted
Михаил Ширшов
Asked:2023-08-16 14:03:30 +0000 UTC2023-08-16 14:03:30 +0000 UTC

将 USRN 坐标转换为地理坐标

  • 772

我试图自己弄清楚，但我已经“飘”在边缘了。有来自 USRN 的摘录/CPT 的坐标，其类型为“ x 536359.03 ”和“ y 2205836.12 ”。我找到了一个可以将这些坐标转换为地理坐标的网站。根据该网站，地理坐标“W 56.37917 ”和“ D 43.83337 ”以及坐标结果是正确的。由于我将来会有很多来自 USRN 的此类坐标，因此我想用 Python 来完成。

我在 USRN 中找到了一个坐标系，称为“ MSK-52 zone 2 Nizhny Novgorod Region [MB:6335202] ”（通过反复试验......）及其网站上的“ PROJ.4 ” - “ +proj =tmerc +lat_0= 0 +lon_0=44.55 +k=1 +x_0=2250000 +y_0=-5714743.504 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +单位=m +no_defs“。我发现我上面指出的“ PROJ.4 ”是根据 GOST 51794-2008 （我查看了文档......有一个移液器）并且它被取消了，现在，据我了解 GOST 32453-2013 ，我没有找到PROJ.4它的

我尝试在网站上拿起“ PROJ.4 ” ，在窗口1（坐标系/投影输入）的网站上我输入了上面写的“ PROJ.4 ”，在窗口2中我拿起了“ PROJ.4”。 4 ”对于俄罗斯联邦来说，但我不认为一切都那么简单。

请帮帮我。

这是我的代码：

import pyproj
# Создаем объект проекции для системы координат Гаусса-Крюгера
proj_gk = pyproj.Proj("+proj=tmerc +lat_0=0 +lon_0=45 +k=1 +x_0=5500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0,0,0 +units=m +no_defs")

# МСК-52 Нижегородская область
proj_msk52 = pyproj.Proj("+proj=tmerc +lat_0=0 +lon_0=44.55 +k=1 +x_0=2250000 +y_0=-5714743.504 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +units=m +no_defs")

# Создаем объект проекции для географических координат (широта, долгота)
proj_geo = pyproj.Proj("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")

# Задаем координаты в системе Гаусса-Крюгера
x = 536359.03
y = 2205836.12

# Производим преобразование координат из системы Гаусса-Крюгера в МСК-52 Нижегородская область
x_msk52, y_msk52 = pyproj.transform(proj_gk, proj_msk52, x, y)

# Производим преобразование координат из МСК-52 Нижегородская область в географические координаты
lon, lat = pyproj.transform(proj_msk52, proj_geo, x_msk52, y_msk52)

print("Широта: ", lat)
print("Долгота: ", lon)
python

2 个回答

  1. 找到一个文章网站。发现一条用户评论：

    引用：从“平面坐标系”转换到 WGS84 的任务分为两部分 - 首先，从投影到获得投影的坐标系的转换，然后从该坐标系到 WGS84 的转换（如果原始坐标系不同）。- 蜘蛛埃克布

    写了这段代码：

    import pyproj

x = 536359.03
y = 2205836.12

# МСК-52 зона 2 Нижегородская область
proj_msk_52 = pyproj.Proj('+proj=tmerc +lat_0=0 +lon_0=44.55 +k=1 +x_0=2250000 +y_0=-5714743.504 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +units=m +no_defs')
lon, lat = proj_msk_52(x, y, inverse=True)
# WGS84
proj_wgs84 = pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs')
# из МСК-52 в WGS84
lon, lat = pyproj.transform(proj_msk_52, proj_wgs84, lon, lat)
print(lat)
print(lon)

    它仍然给出左侧坐标“ 47.475971502889465 14.519519023494398 ”，但您需要获得“ W 56.37917 ”和“ D 43.83337 ”。尽管如此，感觉错误是在proj_msk_52中，代码“PROJ.4”是根据 GOST 51794-2008 指示的，但你需要 GOST 32453-2013

    • 2
  2. Best Answer

    基本上......对我来说这个任务（在Python中）是无法解决的。但我找到了一个很好的转换器，你可以将整个KPT驱动到那里并卸载它。这是该网站的链接。您可以关闭此主题

    • 0

相关问题