Определить систему координат по 5 реперам

Системы координат, проекции, преобразования, привязка
Ответить
Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Определить систему координат по 5 реперам

Сообщение Константин Силкин » 25 дек 2020, 18:03

Добрый день!

Дано:
  1. Координаты 5 реперных точек, заданные одновременно в двух системах координат (см. списки ниже)
  2. Описание одной из систем координат (см. описание ниже)
Требуется:
  • Определить параметры второй, неизвестной системы координат
Список координат
Известная система координат (СК63-X3 зона 5)
No N E H
1 5010815 5376179 ?
2 5014829 5381954 30.79
3 5015468 5380358 5.91
4 5015527 5380142 17.03
5 5015508 5380039 16.43
Неизвестная система координат
No N E H
1 11874.30 14446.59 ?
2 15852.14 20245.98 30.79
3 16501.17 18653.61 5.91
4 16561.50 18438.14 17.03
5 16543.33 18334.44 16.43

Описание известной системы координат

Код: Выделить всё

PROJCS["1963-X3-5", GEOGCS["GCS_Pulkovo_1942", DATUM["D_Pulkovo_1942", SPHEROID["Krassovsky_1940_1948", 6378245.0, 298.3]], PRIMEM["Greenwich", 0.0], UNIT["Degree", 0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting", 300000], PARAMETER["False_Northing", -9214.69], PARAMETER["Central_Meridian", 35.5], PARAMETER["Scale_Factor", 1], PARAMETER["Latitude_Of_Origin", 0.0], UNIT["Meter", 1.0]]
Прошу помощи

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2227
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1500
Откуда: Казань

Re: Определить систему координат по 5 реперам

Сообщение Игорь Белов » 25 дек 2020, 20:54

Статья по теме: Добавление местной координатной системы в GIS

Координаты в первом списке действительно точностью 1 метр? Тогда упражнения не имеют практического смысла.

Кстати, исходная система координат EPSG:7829 "Pulkovo 1942 / CS63 zone X5"

Код: Выделить всё

PROJCRS["Pulkovo 1942 / CS63 zone X5",
    BASEGEOGCRS["Pulkovo 1942",
        DATUM["Pulkovo 1942",
            ELLIPSOID["Krassowsky 1940",6378245,298.3,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4284]],
    CONVERSION["CS63 zone X5",
        METHOD["Transverse Mercator",
            ID["EPSG",9807]],
        PARAMETER["Latitude of natural origin",0.0833333333333333,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",35.5,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",5300000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["northing (X)",north,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["easting (Y)",east,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["Ukraine - 34°E to 37°E onshore"],
        BBOX[44.33,34,52.25,37.01]],
    ID["EPSG",7829]]
The purpose of computing is insight, not numbers

gamm
Гуру
Сообщения: 4044
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1049
Ваше звание: программист
Откуда: Казань

Re: Определить систему координат по 5 реперам

Сообщение gamm » 25 дек 2020, 20:56

с точностью около 0.3 метра это просто перенос

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Определить систему координат по 5 реперам

Сообщение Константин Силкин » 25 дек 2020, 23:22

Спасибо, Игорь!
Рекомендованную вами статью изучил и там оставил вопросы. Статья недостаточно обобщённая, даётся пример по одной зоне. А что делать для другой не понятно.
Для компактности значения в таблицах округлил. Там было несколько знаков после точки. Реальные значения таковы:

Код: Выделить всё

No	N	E
1	5010815.594	5376179.991
2	5014829.581	5381954.912
3	5015468.749	5380358.446
4	5015527.742	5380142.595
5	5015508.927	5380039.007

То, что зона не та, я догадывался. Перепутали геодезисты, которые скинули данные: СК63-X5-3 и СК63-X3-5
Вы представили её в другом формате, спасибо. Но что мне с этим теперь делать? В смысле что даёт мне нового иное представление тех же параметров?

lam
Гуру
Сообщения: 1018
Зарегистрирован: 01 авг 2012, 13:55
Репутация: 241

Re: Определить систему координат по 5 реперам

Сообщение lam » 26 дек 2020, 14:33

Вычисление параметров проекции - https://racurs.ru/program-products/besp ... lozheniya/

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Определить систему координат по 5 реперам

Сообщение Константин Силкин » 26 дек 2020, 17:38

Спасибо, lam!
Но подскажите ещё как понимать описание:
Программа «Вычисление параметров проекции» позволяет вычислить неизвестные параметры проекции Гаусса-Крюгера. В качестве исходных данных для программы служат значения координат Широта-Долгота точки и плоские координаты X,Y с неизвестными значениями долготы осевого меридиана и параллельных переносов. Вычисления будут корректными только в случае, если исходные данные отнесены к одному и тому же эллипсоиду с одинаковой ориентировкой.
Не понятно где взять значения широты и долготы, если система координат прямоугольная

lam
Гуру
Сообщения: 1018
Зарегистрирован: 01 авг 2012, 13:55
Репутация: 241

Re: Определить систему координат по 5 реперам

Сообщение lam » 26 дек 2020, 19:33

Константин Силкин писал(а):
26 дек 2020, 17:38
Не понятно где взять значения широты и долготы, если система координат прямоугольная
Ну так, переведите метры СК63 в градусы.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Определить систему координат по 5 реперам

Сообщение Константин Силкин » 26 дек 2020, 20:27

Это можно было бы, но подскажите ещё, пожалуйста. Если я в GlobalMapper возьму для этого проекцию Geographic, а датум S-42 Russia, то это будет правильно в смысле использования одного референц-эллипсоида с СК63?

lam
Гуру
Сообщения: 1018
Зарегистрирован: 01 авг 2012, 13:55
Репутация: 241

Re: Определить систему координат по 5 реперам

Сообщение lam » 26 дек 2020, 20:53

Константин Силкин писал(а):
26 дек 2020, 20:27
Geographic, а датум S-42 Russia, то это будет правильно в смысле использования одного референц-эллипсоида с СК63?
Да, правильно.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Определить систему координат по 5 реперам

Сообщение Константин Силкин » 26 дек 2020, 22:27

Спасибо!

lam
Гуру
Сообщения: 1018
Зарегистрирован: 01 авг 2012, 13:55
Репутация: 241

Re: Определить систему координат по 5 реперам

Сообщение lam » 27 дек 2020, 10:34

Константин Силкин писал(а):
26 дек 2020, 22:27
Спасибо!
Напишите свою почту.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Определить систему координат по 5 реперам

Сообщение Константин Силкин » 27 дек 2020, 18:20


stout
Участник
Сообщения: 91
Зарегистрирован: 13 апр 2012, 17:09
Репутация: 52

Re: Определить систему координат по 5 реперам

Сообщение stout » 05 янв 2021, 11:57

C такими параметрами:

CMeridian=36.00005641°
FalseNorthing=-5007804.9408
FalseEasting=-22518.6392

координаты получаются такими
arch.png
arch.png (85.99 КБ) 5691 просмотр

stout
Участник
Сообщения: 91
Зарегистрирован: 13 апр 2012, 17:09
Репутация: 52

Re: Определить систему координат по 5 реперам

Сообщение stout » 05 янв 2021, 12:39

Если изменить центральный меридиан на ровно 36°, то смещения будут

FalseNorthing=-5007804.9691
FalseEasting=-22523.0603

смещение центрального меридиана не должно влиять на северное смещение, но что получилось – то получилось (я кажется посчитал с другим порядком аппроксимации)
а с восточным смещением всё просто
-22518.6392 — -22523.0603=4.4211
(36.00005641°-36°)×π/180×6378245×cos(45°20')= 4.4145

Ответить

Вернуться в «Координаты и привязка»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 5 гостей