УСК Рыбинск, СК-76-1, WGS-84...

[Провинциал]

Уважаемые Знатоки.
Участвую во внедрении (загрузка/актуализация данных) Портала ИПД Ярославской областиТестовый вариант Не подумайте, что это реклама
Одним из сервисов является загрузка-выгрузка на портал данных из УСК городов и районов области.
По-моему, родной для Портала является WGS-84. Родной для Рыбинска - УСК Рыбинск. Родной для окружающего нас Рыбинского района СК-76-1.
Моя цель - работать на территории Рыбинска с данными трёх СК (УСК Рыбинск, СК-76-1, WGS-84).
Сам я НЕ геодезист. Среди техподдержки Портала профессиональных геодезистов нет.
Делаем, что умеем и как понимаем... Большей частью - "по наитию"...
Сегодня имеем:
1)Используется GeoTools API
2)Описание проекции СК-76-1, откуда-то взятое, с которым объекты ложатся на свои места хорошо:
PROJCS["76_Ярославская_СК_1",
GEOGCS["GCS_Pulkovo_1942",
DATUM["D_Pulkovo_1942",SPHEROID["Krasovsky_1940",6378245.0,298.3]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Gauss_Kruger"],
PARAMETER["False_Easting",1250117.0],
PARAMETER["False_Northing",-6000009.0],
PARAMETER["Central_Meridian",38.55],
PARAMETER["Scale_Factor",1.0],
PARAMETER["Latitude_Of_Origin",0.1333333333333333],
UNIT["Meter",1.0]]
3)Параметры для аффинного перехода от УСК Рыбинск в СК-76-1
дельта X(декарт)=1268308,32
дельта Y(декарт)=422080,63
дельта ФИ=-0,0038827240000000
(мы вычислили на основании координат Росреестровской границы города в СК-76-1 и УСК Рыбинск)
Говорят, что где-то в закромах Росреестра запрятаны официальные КЛЮЧИ ПЕРЕХОДА.
Видимо, это слишком страшная тайна... Даже не пытались её выведать.
С приведёнными выше параметрами при переводах границы туда-обратно координаты вершин совпадают до сантиметра. Нам достаточно.

Правильно ли я понимаю, что для импорта-экспорта из УСК Рыбинск на портал ИПД нужно дополнить описание "Проекции" СК-76-1аффинным преобразованием?
Подскажите, как в таком случае должна выглядеть "проекция" для УСК Рыбинск для работы GeoTools?

PS1.Если всё не так, тогда
Координата 0,0 Рыбинска в WGS-84
38,858046442347586
58,048472604201365
(переведена из УСК Рыбинск в СК-76-1 аффинным переводом, загружена на Портал ИПД через "проекцию СК-76-1")

PS2.Готов выслушивать высказывания в духе "со свиным рылом да в калашный ряд"...
ибо так оно и есть Главное, чтоб результат был

[Донецков]

По МСК-76 зона 1:
PARAMETER["False_Easting",1250117.0],
PARAMETER["False_Northing",-6000009.0],
- должно быть со всеми 0 до сотен метро

по аффинным параметрам, а их GeoTools в принципе поддерживает?

Скорее всего городская СК сделана на базе СК-42 3 град зоны, со смещениями по осям, если их вычислить не надо никаких поворотов и масштабных Кэфф.

[trir]

0. ЧАВО
1. Моя статья
2. OGC 06-103r3
3. Добавление местной координатной системы в GIS

[Провинциал]

..по аффинным параметрам, а их GeoTools в принципе поддерживает?

Это Вы у меня спрашиваете?
из написанного человеческим языком видел в MapInfo нечто подобное
************************
1)в файле параметров указывается тип проекции =1008 (к типу проекции 8 ( Gauss_Kruger) добавляется признак 1000
2)к файлу параметров добавляется запись ....7, A,B,C,D,E,F
где 7 (метры) тип единицы измерения при аффинных преобразованиях
а следующие параметры - коэффициенты в уравнении
X'=A*X+B*Y +C
Y'=D*X+E*Y+F
************************
A=E=cos(фи)=0.999992462
B=D=sin(фи)=-0.003882714
C=1268308.32
F=422080.63
правда, могу перевернуть плюс и минус, потому как даже это уже сверх моего разумения...
Тогда, для нашего слУчая будет 1008...
..7,0.999992462,-0.003882714,1268308.32,0.999992462,-0.003882714,422080.63
************************
Мой опыт юзания разных программ в разных сферах деятельности подсказывает, что подобное должно быть не только в MapInfo

[trir]

http://docs.geotools.org/stable/usergui ... sform.html

[Провинциал]

Для trir
Мне стало как-то неловко и неуютно
Вы слишком высокого мнения о моём трудолюбии и любознательности
Коллеги, в моей работе появился маленький эпизод: данные в УСК Рыбинск разместить/забрать на ресурсе, работающем в WGS-84 и использующим geo tools
Мне это нужно сделать с минимальными затратами.
Мне кажется, что дополнить работающую "проекцию" СК-76-1 указанием выполнить смещение и поворот - это простейший путь.
Прошу подсказать, что надо дописать проекцию СК-76-1 чтобы она сработала для координат в УСК Рыбинск?

[Провинциал]

Для trir
предыдущий текст был написан ДО Вашего сообщения про аффинные преобразования.
Спасибо.
Теперь дальше.
Нужно проводить последовательно 2 преобразования (проекция, аффинное) по двум наборам параметров?
Или существует правило записи их их в один набор параметров (проекция+аффинное)

[trir]

В WKT это не предусмотренно

Нужно проводить последовательно 2 преобразования (проекция, аффинное) по двум наборам параметров?

да

[Провинциал]

Trir (с большой буквы, заметьте), спасибо. Пока пауза. Будем трясти...

[Донецков]

подобное должно быть не только в MapInfo

- увы, она одна из немногих , которая это использует... , а 90% другого ПО, особенно если оно использует GDAL этого не понимает...

[stout]

Правильно ли я понимаю, что для импорта-экспорта из УСК Рыбинск на портал ИПД нужно дополнить описание "Проекции" СК-76-1аффинным преобразованием?
Подскажите, как в таком случае должна выглядеть "проекция" для УСК Рыбинск для работы GeoTools?

Совсем не обязательно. У вас не просто аффинное преобразование. Оно ортогональное. Что намного грамотнее, чем подгонка аффинным преобразованием, которое в общем случае не является конформным.
С большой долей вероятности эти два последовательных преобразования возможно описать просто проекцией типа Гаусса-Крюгера.
Так как вы говорите, что возможно ошибаетесь в знаках коэффициентов ортогонального преобразования, то приведите координаты 10 точек (можно и больше) в СК-76-1 и в УСК Рыбинск. Понятное дело, что точки совсем не обязательно должны принадлежать реальным объектам. Часть точек желательно разместить на границе области, а ещё лучше — за границей области. Чем больше расстояние между точками, тем точнее результат. При переводе сохраните все значащие цифры, какие возможно. Это также сказывается на точности определения параметров проекции.
Посмотрю, можно ли описать параметрами только проекции.

По МСК-76 зона 1:
PARAMETER["False_Easting",1250117.0],
PARAMETER["False_Northing",-6000009.0],
- должно быть со всеми 0 до сотен метро

Почти на 100% уверен, что это результат несоответствия используемого "D_Pulkovo_1942" реальным значениям для этого места.

[Донецков]

результат несоответствия используемого "D_Pulkovo_1942" реальным значениям для этого места

- тут нужно смотреть, что на что натягивать, неправильную подложку на правильный вектор (если есть возможность корректировать привязку подложки), либо правильный вектор, на неправильную подложку (что реализовано в приведенном примере)

[Провинциал]

Доброго всем дня.
Спасибо за участие в делах наших...

1) stout попросил 10 (или больше) точек - пожалуйста:
1295 точек границы города в УСК Рыбинск и в СК-76-1
Граница между Городом и Районом запрошена через Технокад-муниципалитет у Росреестра
для города Рыбинск (получаем в УСК) и для Рыбинского района (получаем в СК-76-1)
https://yadi.sk/d/MtKI9FHMkfLgK
Обратите внимание: X и Y - декартовые, разделитель ";", дес.разд. "."
Ну, уж количество знаков после дес.раздел. - сколько есть...
На основании этих данных я и вычислил параметры для аффинного перехода от УСК Рыбинск в СК-76-1
дельта X(декарт)=1268308,32
дельта Y(декарт)=422080,63
дельта ФИ=-0,0038827240000000
Уж как уж умею, с помощью Ёкселя...
И с помощью же Ёкселя с этими смещением и поворотом точки пересчитываются туда-обратно...
(Ну, могут появиться изменения в "первых сантиметрах", как я нахватался во время "совещаний" с местными геодезистами старой закваски).
Я догадываюсь, что в каком-то из этих наборов и у Росреестра координаты должны быть не круглыми в сантиметрах, но, ничего другого у меня нет.

2) trir, дал ссылочку на примеры аффинных преобразований
http://docs.geotools.org/stable/usergui ... sform.html
я там не увидел поворота осей.
Можно примерчик смещение+поворот
(от масштабирований, зеркал и обрезаний Бог миловал...)

3) для Донецков Мне нужны точные ВЕКТОРА. Вопрос совмещения фона (подстилки) с векторами это уже другая песня.
Причём, УСК Рыбинск с СК-76-1 должны совмещаться ТАК КАК у Росреестра. Даже если их вариант не самый точный...

[trir]

[Донецков]

Не пытались узнать у Ваших земляков см. ссылка
Если они ставили базу ГНСС в городе, то явно вычисляли параметры на городскую систему