Преобразование прямоугольных пространственных координат в геодезические(ГОСТ Р 51794-2001)

Здравствуйте, столкнулся с такой проблемой, что после преобразования, лишь Y совпадает по графику с начальными данными. Были ли у кого-то ещё такие проблемы и как вы её решили?(сверху графики входных данных)

Код MATLAB:

 

nifarap, позовите "stout"через личку в эту тему.
Можно ещё кого-нибудь из тем про преобразование координат, чьи сообщения понравятся вам.

 

Телепаты все ушли на фронт.
Поэтому для всех оставшихся опишите задачу целиком: что дано, что хотите получить.
По оси Z прямая линия это результат выбранного масштаба отображения или так оно у вас и дано (что очень странно)?

А из того, что написали, сразу бросается в глаза
e = (2*aj - aj^2)^2;
Должно быть просто
e = (2*aj - aj^2);
С точки зрения числа операций (и точности вычислений) это выражение лучше вычислять как e = aj*(2 - aj);

З.Ы. Квадраты величин наподобие этой записываю как "ee". Попробуйте, вдруг понравится. Для какой нибудь ξ² использую префикс: sqr_xi
З.З.Ы. А зачем вам массив синусов и косинусов? Где-то дальше используете?

 

такое описание, чтобы в дальнейшем не возводить в квадрат, "константа" же,а без квадрата нигде не понадобится.

начальная ось Z , этой входные результаты расчёта высоты, странно, что они постоянны, но чтож поделать.

хочу я перевести из nmea, точнее уже из просто градусов, в метры, для наглядного представления погрешности приёма сигнала,а дан исходный файл с данными nmea. выборка и перевод в чисто градусные сейчас не рассматривается, главное разобраться с переводом в метры.

 

Оффтоп

ГОСТ Р 51794-2001 уже давно не действует ...

 

Но смысл-то квадрата первого эксцентриситета. Сам первый эксцентриситет тоже участвует в формулах, но других.

Стоп-стоп, какой такой высоты? Откуда у вас высота? Z в этих формулах это совсем не высота. (Ну, может быть, только в Питере. Это отдельная смешная история про полковников, генералов и о приёмниках РИРВ для вояк.)

Вам не просто в метры нужно, вам надо в топоцентрические координаты. Когда одна ось на север, другая на восток, а третья — вверх.
см. вложение. Да и эллипсоид Красовского здесь ни при чём.

Пофиг, формулы то с тех пор не изменились, только знаки кое-где подправили.

 

а вот тут я не понял, у меня в ГОСТ нет первого эксцентриситета

не правильно выразился. у меня, получается, строка данных такая(ну Вы её знаете)
$GPGGA,083047.999,5548.20969,N,03732.70195,E,1,03,03.6,-14.5,M,14.5,M,,*71
где:
GGA - NMEA Заголовок
083047.999 - UTC время 08:30:47
5548.20969,N - Широта, 55 градусов 48.20969 минуты северной широты
03732.70195,E - Долгота, 37 градусов 32.70195 минуты восточной долготы
1 - тип решения, StandAlone решение

0 - нет решения​

1 - StandAlone,​

2 - DGPS,​

3 - PPS,​

4 - фиксированный RTK,​

5 - не фиксированный RTK,​

6 - использование данных инерциальных систем,​

7 - ручной режим,​

8 - режим симуляции​

03 - количество используемых спутников
03.6 - геометрический фактор HDOP
-14.5,M - высота над уровенем моря в метрах
14.5,M - высота геоида над эллипсоидом WGS 84

И вот как раз эту высоту геоида я и беру в Z, просто она в начале как H используется,а после этих преобразований в Z.

 

Первый эксцентриситет часто называют просто эксцентриситет.
Особенно тогда, когда в тексте не используется второй эксцентриситет e'
e'² = e²/(1–e²)

Вообще-то высота геоида в данной точке для данной модели геоида величина постоянная. Вам нужна геодезическая высота точки. Находится как сумма ортометрической высоты (высота над уровнем моря) и высоты геоида над эллипсоидом. ( http://www.astronet.ru/db/msg/1190817/node25.html )

 

получается, мне нужно эту формулу для эксцентриситета использовать?

хорошо, спасибо, но всё равно различия графиков входных данных и после преобразования будут разными.

 

нет.
Уже писал:

И не забудьте про высоту.

 

Оффтоп

Это смотря для чего. Для ГСК-2011 и ПЗ-90.11 могёт быть и не совсем так ...

 

Оффтоп

В первом сообщении приведена формула

И это точно не высота.

 

Оффтоп

Видимо имелось равенство высот эллипсоидальной и БС77 в центре круглого зала Пулковской обсерватории для СК-42 ...

 

график latitude схож с Z
но всё равно это не то, может подскажете, как мне эти градусы, градусы, минуты (Широта, Долгота, Высота) всё таки можно перевести в реальные метры?(реальные,в смысле на данное время, чтобы формулы были актуальны З.Ы. отсылка на Родичкин)

 

Оффтоп

Нет. Как я уже сказал, всё было намного смешнее. На индикаторе приёмника выдавались геоцентрические прямоугольные координаты. Во время приёмки(демонстрации) генералы с умным видом кивали головами, слушая объяснения, пока один полковник вдруг не ляпни: "а почему это у вас координата Z несколько тысяч метров? Ведь известно, что Z в Питере — это высота. А таких высот в окрестностях Питера нет." Сам полковник — мужик неплохой, и специалист хороший, но в своей области, в фотограмметрии. Как его занесло экспертом на испытания GNSS-приёмника — один бог ведает. И все дальнейшие робкие попытки разработчика объяснить разницу между различными системами координат наткнулись на железный вердикт генералов, приёмку прекратить, замечание к завтрашнему дню исправить и доложить. Ну и пришлось разработчику firmware выдать на индикатор φ, λ, h. А уж такие слова, как широта, долгота и высота генералы помнили, сверились с картой и остались довольны собой. А разработчика похвалили: "Вот видите, ведь можете, если приказать!"
С тех пор разработчик софта при определённых ситуациях иногда вспоминает эту фразу: "Как известно, зет в Питере — это высота!"

--- Сообщения объединены, , Оригинальное время сообщения: ---

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

 

Это немного вольное трактование протокола.
03.6 - ?геометрический фактор HDOP? - обычно это точность в горизонте (2D)
+ на матлабовских картинках широта и долгота прямо из протокола без разборки на градусы и минуты. Видимо по этой причине скачек по долготе...

 

перевод в чисто градусы присутствует,с этим нет проблем.
а это описание из инструкции к приёмнику брал, чтобы не "наврать"

 

Т.Е. задача стоит визуальзировать nmea лог приемника?

 

да, с переводом градусных данных в метры, для наглядного представления погрешности

 

Лог в студию.