Страница 1 из 1

Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:17
codemonkeydeveloper
Имеются растровые карты в Системе координат пулково 1942 в одной из зон, по этим картам в БД заносились объекты в той же системе координат.
Сейчас производим перевод координат в БД в свою местную систему координат, с отличным осевым меридианом.
У карт в виде картинок возможно перевести только координату верхнего левого угла.
Вопрос лягут ли объекты переведенные в бд на картинки или же будет расхождение?
И еще вопрос сохранят ли площадные объекты при переводе из пулково 1942 в местную систему координат свои геометрические размеры или будет искажение?

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:28
trir
будет расхождение
но это зависит от параметров МСК - есть аффиные преобразования, растояние между севым мередианном
возможно перевести только координату верхнего левого угла.
поворот тоже надо изменить

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:42
gamm
codemonkeydeveloper писал(а):
27 ноя 2017, 09:17
Имеются растровые карты в Системе координат пулково 1942 в одной из зон, по этим картам в БД заносились объекты в той же системе координат.
(1) т.е. на картах есть координаты ...
Сейчас производим перевод координат в БД в свою местную систему координат, с отличным осевым меридианом .
что означает сия фраза? Вы пересчитываете в градусы на Красовском, а потом в свою систему (или просто полиномом пересчитываете по точкам, или еще как-то)?
(2)У карт в виде картинок возможно перевести только координату верхнего левого угла.
это прямо противоречит (1) ...
Вопрос лягут ли объекты переведенные в бд на картинки или же будет расхождение?
лягут, если преобразование обратимо, и либо вы трансформируете растровые карты в новую систему, либо будете при выводе на лету трансформировать координаты в старую.
И еще вопрос сохранят ли площадные объекты при переводе из пулково 1942 в местную систему координат свои геометрические размеры или будет искажение?
что такое "искажение"? любая проекция есть искажение формы на Земле, просто у разных проекций искажения разные. Если у вас разные проекции, то и искажения изменятся, но никто не мешает все внутри делать в старой системе.

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:44
codemonkeydeveloper
но это зависит от параметров МСК - есть аффиные преобразования, растояние между севым мередианном
На сколько я понимаю там только меняется осевой меридиан, там мск для нашего района просто.
Вроде проверил пару точек и на расстоянии около 300 км разница у меня получилась около 5 метров. Это нормально или все тютелька в тютельку совпасть при смене только осевого меридиана?

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:48
codemonkeydeveloper
что такое "искажение"? любая проекция есть искажение формы на Земле, просто у разных проекций искажения разные. Если у вас разные проекции, то и искажения изменятся, но никто не мешает все внутри делать в старой системе.
Ну например есть квадрат 100 на 100 километров в пулково 1942 года в одной из зон, мы переведем координаты вершин квадрата в местную систему координат с отличным осевым меридианом будет ли это квадрат квадратом 100 на 100 км или станет квадратом 101 на 101 км? И не станет ли он ромбом к примеру. Есть например граница земель в виде многоугольника по площади она приличная до 50 км в ширину , лягут ли мои переведенные данные на растровую картинку если я просто её перетащу на новые точки или все таки будет расхождение?

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 09:49
codemonkeydeveloper
что означает сия фраза? Вы пересчитываете в градусы на Красовском, а потом в свою систему (или просто полиномом пересчитываете по точкам, или еще как-то)?
Да в градусы потом в свою систему.

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 10:01
trir
Это нормально
там искажения зависят от расстояния до осевого мередиана
Несмотря на то, что площадь зоны в проекции (зоны Гаусса) будет увеличенной, относительные искажения длин в отдалённых от среднего меридиана точках экватора на границе зоны составит 1/800. Максимальные искажения длин в пределах зоны составляет +0,14 %, а площадей — +0,27 %

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 10:15
codemonkeydeveloper
там искажения зависят от расстояния до осевого мередиана
Ну т.е. если зона это отход от меридиана осевого на 330 км, то на границах зоны погрешность до 462 метров?

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 10:26
trir
да
1.4 м на километр, по этому, когда работают на границе зон - используют промежуточную 3-ёх градусную зону
в МСК применяют аффиные преобразования для компенсации этой погрешности

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 10:44
codemonkeydeveloper
в МСК применяют аффиные преобразования для компенсации этой погрешности
А где про это поподробнее почитать не подскажите?

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 27 ноя 2017, 10:55
trir
в учебнике

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

+proj=tmerc +lat_0=0 +lon_0=63 +k=1 +x_0=500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs

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

+proj=tmerc +lat_0=0 +lon_0=66 +k=1 +x_0=500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +units=m +no_defs
834185.886865;97.8443785721=>500067.933909;97.7093271287
833185.886865;97.8443785721=>499069.310052;97.7101342847
1000=>998.623857
разница - 1.376143

Re: Перевод объектов из пулково 1942 в местную систему координат

Добавлено: 28 ноя 2017, 08:34
Игорь Белов
codemonkeydeveloper писал(а):
27 ноя 2017, 09:17
У карт в виде картинок возможно перевести только координату верхнего левого угла.
Вопрос лягут ли объекты переведенные в бд на картинки или же будет расхождение?
Координат одной точки недостаточно, нужны также разрешения по каждой координате.
Картинки должны развернуться каждая на свой угол. Пожалуй, проще всего их трансформировать в новую координатную систему специальными программами. GDAL подойдёт.
codemonkeydeveloper писал(а):
27 ноя 2017, 09:17
И еще вопрос сохранят ли площадные объекты при переводе из пулково 1942 в местную систему координат свои геометрические размеры или будет искажение?
Если площади объектов посчитаны в проекции, т. е. на плоскости, то непременно изменятся.
Чтобы площади не зависели от проекции, их следует вычислять на эллипсоиде.

Практически достаточно рассчитать площадь в проекции и редуцировать на эллипсоид:

S₂ = S₁ / m² ,

где S₁ — площадь, вычисленная на плоскости; S₂ — площадь на эллипсоиде; m — масштаб изображения

m = 1 + y² / 2 R² + y⁴ / 24 R⁴ + y⁶ / 720 R⁶,

где y — расстояние центроида объекта от осевого меридиана y = Y − Y₀; R — средний радиус кривизны эллипсоида

R = c / (1 + e′² cos² φ₀) .

Здесь φ₀ — широта центроида объекта; c — радиус полярной кривизны; e′ — второй эксцентриситет.