Тестовые данные для проверки реализации ГОСТ-Р51794

[drug]

Ну насчет разгонять я может и громко сказал, но мое имхо, что уровень военной науки и ВПК недостаточно высок. То, что они порой творят это просто тихий ужас.
Насчет задачи, которую я хочу решить - да, мне нужно всего лишь проверить реализацию ГОСТа. Т.е. я подаю на вход те же данные, что и прислал ErnieBoyd и сравниваю. Уточнять параметры ГОСТа - я очень далек от этого, вы мне очень сильно польстили.
На данный момент получилось, что пересчет от geodetic WGS84 в decart WGS84 совпадает до 8го знака после запятой. А вот при пересчете из географических WGS84 в пространственные ПЗ9002 расхождение до 200 м между моей реализацией ГОСТа и данными ErnieBoyd. В свете чего мы начинаем склоняться к решению воспользоваться proj4, потому что верификация начинает занимать больше времени чем мы можем себе позволить, а совесть мне не позволит просить людей создавать качественные и всеобъемлющие тест-кейсы. А разобравшись с proj4 мы сможем самостоятельно тесты создать для верификации реализации в последующем. В любом случае, большое спасибо за помощь! Очень полезный форум и приятное сообщество.

[Игорь Белов]

пересчет от geodetic WGS84 в decart WGS84 совпадает до 8го знака после запятой … при пересчете из географических WGS84 в пространственные ПЗ9002 расхождение до 200 м

Отлично, широта/долгота WGS 84 замечательно превращаются в XYZ. Дальше надо конвертировать XYZ WGS 84 в XYZ ПЗ-90.02. Открываем Приложение В и читаем:

∆X = -0.36, ∆Y = 0.08, ∆Z = 0.18, ωx = 0, ωy = 0, ωz = 0, ∆S = 0

Очевидно, имеет место параллельный перенос, и координаты XYZ ПЗ-90.02 любых точек должны быть меньше соответствующих XYZ WGS 84 точно на -0.36, 0.08, 0.18.

Либо я что-то неправильно понимаю, либо Ваш модуль трансформирования между СК выдаёт странное.

[drug]

Это вы используете трансформацию декартовых координат - ее я, кстати, не проверял, сейчас сделаю. Ошибка у меня возникла при пересчете по следующей цепочке географические WGS84 -> географические ПЗ90 -> географические на Красовском (СК95) -> декартовы на Красовском (СК95).

Проверил трансформацию декартовых координат по цепочке WGS84->ПЗ90.02->Красовский (в предыдущем посте я соврал, я имел в виду не ПЗ90.02, а Красовский, т.е. СК95). Расхождение по X на ~45м, по Y на ~ 261м, по Z ~161м.

[Игорь Белов]

Так WGS84->ПЗ90.02 даёт правильный результат? Проблема на этапе ПЗ90.02->СК-95? Если так, то что даёт обратное преобразование СК-95->ПЗ90.02?

[drug]

Так WGS84->ПЗ90.02 даёт правильный результат?

Не знаю. Вы же не прислали тестовых данных по этому преобразованию... Как я понял, первый набор точек не подходит, потому что эллипсоид ПЗ90 натянут на оси Красовского, а во втором наборе нет XYZ в ПЗ90.02, только в СК95.

[Игорь Белов]

первый набор точек не подходит, потому что эллипсоид ПЗ90 натянут на оси Красовского

Как раз наоборот, перечитайте. Ну и посчитайте разности XYZ между WGS 84 и ПЗ-90.02 в первом наборе.

во втором наборе нет XYZ в ПЗ90.02

Если нужны, можете посчитать в Excel'е:

координаты XYZ ПЗ-90.02 любых точек должны быть меньше соответствующих XYZ WGS 84 точно на -0.36, 0.08, 0.18

[drug]

Ок, понял, разобрался - конвертация из пространственных WGS84 в ПЗ90.02 работает с точностью до 9го знака после запятой.

[drug]

Да, получается ошибка возникает при переходе от ПЗ90.02 к СК95. Причем ошибка возникает и при обратном переходе на те же значения. Коэффициэнты верные, заданы в метрах и радианах. Формула пересчета простая (перемножение матриц), ошибки там я не нашел. Хм...
Хотя стоп - судя по всему ошибка равна удвоенному значению переносов по осям координат, значит ошибка у меня. Где-то вместо сложения я отнимаю и наоборот.

[drug]

ErnieBoyd, а вы не могли бы привести коэффициенты трансформации из XYZ ПЗ90.02 в СК95? Просто я не могу найти причину - все коэффициенты верные, а перемножение матриц оно и в Африке перемножение матриц. Если поиграться коэффициентами то можно уменьшить разницу, но вроде как ГОСТ это не то место где можно играть...

[Игорь Белов]

Приложение А

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

из       в        ∆X    ∆Y      ∆Z     ωx ωy   ωz    ∆S
СК-42 -> ПЗ-90.02 23.93 -141.03 -79.98 0 -0.35 -0.79 -0.22
СК-95 -> ПЗ-90.02 24.83 -130.97 -81.74 0  0    -0.13 -0.22

[drug]

Нашел ошибку у себя - я эти коэффициенты использовал при переходе от ПЗ-90.02 к СК95, т.е. наоборот. Теперь результаты лучше, но тем не менее - если по Z расхождение начинается с 5го знака после запятой, то по X и Y расхождение 2.5 м. Это может быть следствием разных формул или нужно еще одну ошибку искать у себя?

З.Ы. я правильно предположил, что если к этим коэффициентам поменять знак, то их можно использовать для пересчета от ПЗ-90.02 в СК95?

[Игорь Белов]

я правильно предположил, что если к этим коэффициентам поменять знак, то их можно использовать для пересчета от ПЗ-90.02 в СК95?

Правильно, в ГОСТе это написано русским по белому.