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

Пользователь

Сообщений: 3 Регистрация: 09.04.2015

09.04.2015 14:55:04

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

Код MATLAB:

Цитата

a = 6378245.0; % Большая (экваториальная) полуось
b = 6356863.019; % Малая (полярная) полуось
aj = 1/298.3; %сжатие эллипсоида
e = (2*aj - aj^2)^2;

for i=1:n
siB(i) = sin(B(i));
siL(i) = sin(L(i));
coB(i) = cos(B(i));
coL(i) = cos(L(i));

siBd(i) = siB(i).^2;

coz(i) = e*siBd(i);
N(i) = a / sqrt(1 - coz(i));

N = N.';
coz = coz.';

x(i) = (N(i) + H(i))*coB(i)*coL(i);
y(i) = (N(i) + H(i))*coB(i)*siL(i);
z(i) = [(1 - e)*N(i) + H(i)]*(siB(i));

end

Елена Кузнецова

Администратор

Сообщений: 7282 Регистрация: 11.07.2011

09.04.2015 15:25:24

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

Проекция Гаусса-Крюгера (UTM) имеет ограничения области пересчета координат по долготе. Чем дальше от осевого меридиана, тем больше погрешность пересчета. При удалении на несколько десятков градусов
погрешность может быть десятки метров и более.
В пределах зоны погрешность может быть в мм.

Sergey Nifarapovich

Пользователь

Сообщений: 3 Регистрация: 09.04.2015

09.04.2015 16:18:22

Цитата

Елена Кузнецова пишет:

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

Проекция Гаусса-Крюгера (UTM) имеет ограничения области пересчета координат по долготе. Чем дальше от осевого меридиана, тем больше погрешность пересчета. При удалении на несколько десятков градусов

погрешность может быть десятки метров и более.

В пределах зоны погрешность может быть в мм.

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

Елена Кузнецова

Администратор

Сообщений: 7282 Регистрация: 11.07.2011

10.04.2015 10:01:19

Цитата
Sergey Nifarapovich пишет: а можно это как-то учитывать при пересчете? или же может другими формулами пользоваться? мне нужно лишь градусы в метры перевести, для наглядного представления погрешности приёмника

Если Вам нужно отобразить территорию с протяженностью по долготе более 1000 км, то можно применить цилиндрические,
азимутальные или конические проекции.
Можно воспользоваться другими проекциями. Выбор проекции зависит от размеров территории и требований к виду представленных данных.

Sergey Nifarapovich

Пользователь

Сообщений: 3 Регистрация: 09.04.2015

10.04.2015 10:44:26

Цитата

Елена Кузнецова пишет:

Если Вам нужно отобразить территорию с протяженностью по долготе более 1000 км, то можно применить цилиндрические,

азимутальные или конические проекции.

Можно воспользоваться другими проекциями. Выбор проекции зависит от размеров территории и требований к виду представленных данных.

мне нужно просто перевести данные из NMEA протокола (Широта(градусы), Долгота(градусы), Высоты(метры)) в метры

Oleg Belenkov Администратор Сообщений: 4115 Регистрация: 05.06.2008	#6 10.04.2015 12:13:04 Метры в какой системе? Все вышесказанное относится к метрам на плоскости в проекции Гаусса-Крюгера. Она тоже описана в ГОСТ Р 51794-2007. У Вас пересчет в прямоугольные пространственные координаты. Пересчет должен быть взаимнооднозначным с учетом точности формул. О точности в ГОСТ сказано: Для преобразования пространственных прямоугольных координат в геодезические необходимо проведение итераций при вычислении геодезической широты. При преобразованиях координат в качестве допуска прекращения итеративного процесса принимают значение 10-4 угл. с. В этом случае погрешность вычисления геодезической высоты не превышает 0,003 м.