авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

Разработка и исследование алгоритма математической обработки геодезических сетей в разных системах координат

На правах рукописи

Клыпин Игорь Андреевич

РАЗРАБОТКА И ИССЛЕДОВАНИЕ АЛГОРИТМА

МАТЕМАТИЧЕСКОЙ ОБРАБОТКИ

ГЕОДЕЗИЧЕСКИХ СЕТЕЙ

В РАЗНЫХ СИСТЕМАХ КООРДИНАТ

Специальность 25.00.32 — Геодезия

АВТОРЕФЕРАТ

диссертации на соискание учёной степени

кандидата технических наук

Москва — 2011

Работа выполнена на кафедре геодезии Московского государственного университета геодезии и картографии (МИИГАиК)

Научный руководитель: доктор технических наук, профессор Маркузе Юрий Исидорович

Официальные оппоненты: доктор технических наук, профессор Яшкин Станислав Николаевич кандидат технических наук, доцент Калинова Елена Владимировна

Ведущая организация: Московский государственный университет путей сообщения (МИИТ)

Защита диссертации состоится «_» 2011 года в часов на заседании диссертационного совета Д.212.143.03 при Московском государственном университете геодезии и картографии (МИИГАиК) по адресу: 105064, г. Москва, Гороховский пер., 4 (Зал заседаний Учёного совета).

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

Автореферат разослан «_» _ 2011 г.

Учёный секретарь Климков Ю.М.

диссертационного совета

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Построение геодезических сетей основано на линейных, угловых и спутниковых измерениях, но чаще всего на их комплексном использовании, когда традиционные геодезические методы сочетаются с новейшими спутниковыми технологиями. Кроме того, в прикладной геодезии очень часто геодезические сети строятся как самостоятельные, а координаты их пунктов определяются в местной системе координат. Особенно актуально это для больших и быстро развивающихся городов, когда требуется координаты всех имеющихся на обширной территории пунктов привести как можно более в короткие сроки к единой системе координат. При этом часто возникает необходимость объединения геодезических сетей, построенных в разных системах координат, что, в свою очередь, требует развития и совершенствования методов математической обработки геодезических сетей, построенных в разных системах координат.

Цель диссертации — разработка алгоритма математической обработки геодезических сетей, построенных в разных системах координат, на основе современных компьютерных технологий.

Задачи исследования:

1. Обзор ситуации и анализ существующих подходов к математической обработке геодезических сетей, методов и задач уравнительных вычислений;

2. Разработка алгоритма уравнивания и объединения геодезических сетей, построенных в разных системах координат;

3. Моделирование геодезических сетей, построенных в разных системах координат, и исследование на их основе данного алгоритма.

Научная новизна работы заключается в следующем:

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

2. Выполнено исследование по использованию исходной и линеаризованной формулы перехода из одной системы координат в другую. Исследование показало, что линеаризация не оказывает существенного влияния на точность определения координат.

Практическая ценность работы состоит в том, что разработанный алгоритм позволяет выполнять объединение геодезических сетей, построенных в разных системах координат, с контролем грубых ошибок измерений и оценкой точности параметров преобразования.

Апробация работы. Основные выводы и положения диссертационной работы докладывались автором в 2008–2010 г.г. на:

· 63-ей научно-технической конференции студентов, аспирантов и молодых учёных МИИГАиК (2-3 апреля 2008 года);

· 64-ой юбилейной научно-технической конференции студентов, аспирантов и молодых учёных МИИГАиК, посвящённой 230-ой годовщине со дня его основания (7-8 апреля 2009 года);

· 65-ой научно-технической конференции студентов, аспирантов и молодых учёных МИИГАиК, посвящённой 65-летию победы в Великой Отечественной войне (6-7 апреля 2010 года).

Публикации. Содержание и результаты диссертационной работы освещены в 3 статьях, из них 2 — в издании, рекомендованном ВАК по специальности диссертации.

Объём диссертации. Диссертация состоит из введения, трёх глав, заключения, списка литературы и приложений. Общий объём работы составляет 90 страниц машинописного текста, включая 17 таблиц, 6 рисунков и 4 приложения. Список использованных источников включает в себя 75 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

В главе 1 «Обзор ситуации, методов и задач математической обработки геодезических построений» приводятся сведения о современном состоянии государственных геодезических систем координат в России. Показано, что, несмотря на создание государственной геодезической сети ФАГС и ВГС, которая реализует общеземную систему координат СПК-2005 (пространственную), последняя остаётся фактически недоступной для рядового потребителя.



Показано, что, так как трансформирование координат пунктов спутниковых и наземных сетей с использованием опубликованных на текущий момент параметров преобразования не удовлетворяет практическим требованиям, целесообразно выполнять совместное уравнивание геодезических сетей с определением параметров преобразования.





При этом для преобразования координат из одной системы в другую, используется формула, которую можно представить в виде + (1 + ) = исходной системы связи для идентичного пункта i как (1), где m — разница в линейных масштабах систем координат;

матрица П содержит в себе малые углы поворота и переводит координаты спутниковых построений в систему координат наземной геодезической 1 сети:

= 1 (2).

Для каждой точки i векторы Ti и Si равны = = и, (3) но относятся к системам координат наземной (T) и спутниковой (S) геодезических сетей соответственно.

Вектор a1 постоянен для всех точек и составлен из трёх параметров сдвига:

= (4), представляющих собой, по сути, координаты начала системы координат, в которой выполнены наземные геодезические построения, в системе координат, в которую были переведены координаты пунктов спутниковых построений.

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

Далее формулируются современные задачи уравнительных вычислений:

· уравнивание с учётом ошибок исходных данных при построении геодезических сетей в несколько стадий или уравнивании большой сети постепенным её наращиванием;

· уравнивание с контролем грубых ошибок измерений;

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

· обеспечение устойчивости вычислительного процесса, особенно в случаях ограниченного числа пунктов, когда вместе с вычислением поправок в приближённо известные величины необходимо определять также параметры преобразования;

· и ряд других, не менее важных задач.

Рассмотрены современные методы уравнивания: способ условий с дополнительными неизвестными и рекуррентный способ.

Для способа условий с дополнительными неизвестными (, ) = 0, (5) где y — вектор измерений, x — вектор дополнительных неизвестных, показана возможность контроля грубых ошибок измерений, если вектор () приближённых значений дополнительных неизвестных представить в виде функции общего вида = () () (6) вектора измерений.

При таком подходе математическое ожидание вектора невязок будет равно нулю, а обратная весовая матрица будет определяться как =( ) (7), где,, () () =, =, () (8) =, () (9) а сам вектор невязок определяется согласно =, () (10).

Контроль грубых ошибок по невязкам будет выполняться как =± Q (11) доп.

.

Для рекуррентного уравнивания, разработанного проф. Маркузе Ю.И., приведены основные формулы = (12), = (13), = + = + (14), а также их частный случай, когда учитывается только одно некоррелированное с другими измерение:

= ;

1 =+ ;

=.

(15) Перечислены преимущества данного способа:

· включение и исключение измерений из уже уравненной сети;

· определение обратной весовой матрицы неизвестных не требует составления нормальных уравнений;

· контроль грубых ошибок измерений и исходных данных.

Также рассмотрен метод сопряжённых градиентов при решении систем линейных уравнений вида =, Т (16) где — положительно определённая матрица.

Решение в таком случае находится по следующей последовательности формул:

=, =, = +, Т ) ( ) (,, = = ), (,Т ) (,Т =+, (17) ) (,) (,Т.

= = ( ) (,Т ), При этом значение вектора уточняется на каждом шаге итерации по формуле =. (18) Теоретически точное значение вектора приближённых значений параметров преобразования получается после числа итераций, равного числу t неизвестных.

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

В главе 2 «Разработка алгоритма уравнивания и объединения геодезических сетей» предлагается вариант структуры алгоритма совместной обработки геодезических измерений, вычисления в котором определяются автоматически, исходя из имеющихся данных. Структура алгоритма представлена на рис. 1.

Рис. 1 – Структурная схема программного алгоритма совместного уравнивания геодезических сетей Исходными данными, включаемыми в совместную обработку, будут являться, с одной стороны, геоцентрические координаты, полученные в результате постобработки спутниковых измерений и их точностные характеристики, с другой стороны — плановые координаты xi, yi наземных геодезических построений для этих же пунктов со своими точностными характеристиками, переведённые в прямоугольную пространственную геоцентрическую систему координат по известным формулам перехода, а также нормальные высоты, для перехода от которых к высотам геодезическим, в зависимости от необходимой точности, можно воспользоваться картами аномалий высот, используя соотношение = +, (19) где Hi — геодезическая высота. Или же использовать одну из моделей геоида. Например, EGM2008.

Исходным элементарным блоком данных здесь будут выступать координаты пунктов, поэтому при упрощённой формальной грамматике в соответствии с нотацией Бэкуса–Наура структуру данных любой сети для нашего случая можно упрощённо описать следующей грамматикой:

1) проект::=имя сети{пункт} ;

2) пункт::=[координата X][координата Y][координата Z].

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

Такого подхода и такой совокупности данных вполне достаточно, чтобы программно определить топологию сети, её вид, а также схему вычислений.

Приведённые к линейному виду условные уравнения связи с дополнительными неизвестными при совместном уравнивании геодезических построений будут иметь вид:

+ =0 (20) для каждого идентичного пункта i, с матрицей 1 0 00 =0 1 0 0 01, (21) составленной из координат, которые относятся к спутниковой сети.

Векторы поправок VT и VS здесь равны = = и (22) соответственно.

Вектор содержит в себе поправки к параметрам преобразования и может быть представлен следующим образом:

=.

= (23) Вектор невязок W будет равен = () (24) () для каждого идентичного пункта i, откуда вектор приближённых значений параметров преобразования по формуле () = () (25), где матрицы G1, и содержат в себе уже уточнённые значения координат пунктов.

Матрица G1 здесь представляет собой верхний блок матрицы GS, если её представить в виде = (26).

Далее с целью оптимизации и улучшения обусловленности задачи можно уточнить вектор приближённых значений параметров преобразования, применив итерационный способ нахождения решения по методу сопряжённых градиентов по формулам (17), для чего примем =( ), (27) откуда = (28).

При этом можно говорить о выполнении упрощённого уравнивания, так как попутно уточняются и координаты пунктов, хотя и с нарушением строгости уравнивания.

Тогда для условных уравнений (20) результатам измерений будет соответствовать обратная весовая матрица = (29).

Далее представим вектор приближённых значений параметров преобразования в виде = () () (30) и сформируем матрицу перехода = (31) с целью контроля грубых ошибок измерений, так как в противном случае математическое ожидание вектора невязок оказывается не равно нулю.

В нашем случае матрица перехода будет иметь все нулевые блоки, кроме тех, которые соответствуют пунктам, участвующим в вычислении приближённых параметров преобразования:

0 0 = ( ) (32), После её введения получим обратную весовую матрицу :

= (33), которая в нашем случае содержит блоки, представляющие собой обратные весовые матрицы наземной геодезической сети, спутниковых построений и параметров преобразования:

0, =, (34).

,, Обратная весовая матрица вектора невязок QW будет определяться по формуле =( ) (35), где,, () () =, = () (36) за вычетом условного уравнения, выделенного для определения приближённых значений параметров преобразования.

Далее, используя формулы рекуррентного уравнивания, можно выразить вектор поправок V как = = (37), где — матрица рекуррентного уравнивания, в которой матрица A для нашего случая будет равна =( );

(38) N2 — также матрица рекуррентного уравнивания, которая в нашем случае совпадает с матрицей QW:

=( ) = = = (39).

Ковариационную матрица вектора невязок можно выразить как = (40), откуда допустимую невязку можно оценить по известной формуле =± доп.

. (41) Уравненные значения координат и параметров преобразования теперь несложно получить как = +, (42) причём векторы и состоят из уточнённых координат пунктов сетей, полученных методом сопряжённого градиента.

После этого необходимо лишь преобразовать вектор в вектор по формуле = + (43) или + (1 + ) = (44).

Преобразованные координаты для идентичных пунктов должны совпадать. Также контролем решения задачи будет равенство нулю повторно вычисленного вектора невязок W2.

Обратная весовая матрица уравненных величин будет найдена, исходя из следующей формулы:

( ), = (45) или же, согласно теории рекуррентного уравнивания, учитывая (13):

= (46).

Контролем при нахождении обратной весовой матрицы уравненных величин будет соблюдение равенства =( ) (47).

Учитывая диагональные элементы матрицы весов уравненных величин, находятся средние квадратические ошибки отметок каждого пункта, а также средние квадратические ошибки параметров преобразования согласно известной формуле = (48).

Квадратичная форма будет вычисляться по формуле рекуррентного уравнивания для i-ой группы измерений:

= Т (49), где — квадратичная форма, полученная при уравнивании спутниковых построений.

В случае уравнивания плановых сетей уравнения связи для общих пунктов i в линейном виде следует записать в виде:

+ = 0, (50) где = (51) = cos 1, при sin, (52) = (53) 1 m — масштабный фактор, — угол поворота, = 0. (54) В случае уравнивания высотных сетей уравнения связи для + = общих пунктов i в линейном виде следует записать в виде:

(55) = (1 ) для каждого идентичного пункта i, а матрица — (56) =, или (57) если не учитывать масштабный фактор.

Отдельно следует заметить, что в качестве высот H здесь могут выступать либо геодезические, либо нормальные высоты.

По данному алгоритму было выполнено уравнивание моделей нивелирных (рис. 2) и плановых (рис. 3) сетей, построенных в разных системах координат, а также объединение наземной и спутниковой геодезических сетей (рис. 4).

При совместном уравнивании нивелирных сетей значение a = 89,9970, параметра преобразования оказалось равно:

(58) =24.3 мм.

а его средняя квадратическая ошибка составила (59) t3 t r 100,00 100, 102, 0 r 101,03 r 107, r 107, r r5 105, r 107,09 r 100, 0 108, 101, 7 t 100, t1 100,00 102, 96, t 11, 10, 1 6,804 12, t 10, s6 s 13,274 9, s3 s 10,536 16, s 12,286 s 15, s 10, s 14, t 10, Рис. 2 – Схема нивелирных сетей T2 T T2 P P T T3 T4 P3 P T P4 P T T Рис. 3 – Схема плановых сетей преобразования получились следующими: a = +49,735 м, a = +99,932 м, При совместном уравнивании плановых сетей значения параметров = 179°59 54,3, = 1,234 = 0,012 м, = 0,013 м, = 1,4, = 1,134 10.

. Их средние квадратические ошибки составили значения параметров преобразования оказались равны a = +26,684 м, При объединении наземной и спутниковой геодезических сетей a = 132,595 м, a = 76,228 м, = 0,389, = 0,426, = 0,874, = 1,234 10, а их средние квадратические ошибки — = 0,009 м, = 0,011 м, = 0,012 м, = 0,005, = 0,004, = 0,008 и = 5,821 10.

Рис. 4 – Схема сети Глава 3 «Программная реализация алгоритма и его анализ»

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

При этом рассмотренный в главе 2 алгоритм был реализован как программный модуль “CGM” (от англ. "conjugate gradient method" — метод сопряжённых градиентов), написанный лично автором на языке C++.

Для геодезических сетей, представленных на рис. 2-4, было произведено сравнение уравненных результатов, полученных в программе “CGM” и в программах, разработанных проф. Маркузе Ю.И., по уравниванию нивелирных (“COMBINE1”) и плановых (“COMBINE4”) сетей и объединению наземных и спутниковых сетей (“GPS-3D”).

Как видно из таблиц 1-3, результаты уравнивания в каждом случае практически совпадают.

Таблица 1 – Сравнение результатов уравнивания, полученных в программах “CGM” и “GPS-3D” № № по X, м по Y, м по Z, м по X, м по Y, м по Z, м п/п п/п –0,001 –0,001 –0, 1 +0,003 +0,003 12 +0, –0, 2 +0,002 +0,000 +0,002 13 +0,003 +0, –0,001 –0, 3 +0,000 14 +0,001 +0,000 +0, –0,000 –0,002 –0, 4 +0,001 15 +0,000 +0, –0,000 –0,003 –0,000 –0,002 –0, 5 +0,001 –0,002 –0,001 –0,000 –0, 6 +0,001 +0,001 –0, 7 +0,000 +0,002 +0,000 18 +0,001 +0, –0,003 –0, 8 +0,002 19 +0,001 +0,001 +0, –0,000 –0,002 –0, 9 +0,000 +0,003 20 +0, –0,002 –0,000 –0,002 –0, 10 +0,001 21 +0, –0, 11 +0,001 +0, Таблица 2 – Сравнение Таблица 3 – Сравнение результатов уравнивания, результатов уравнивания, полученных в программах полученных в программах “CGM” и “COMBINE1” “CGM” и “COMBINE4” по H, м по H, м по x, м по y, м r1 +0,000 s1 +0,001 T1 +0,002 +0, –0,002 –0, T –0, r2 +0,000 s T3 +0,000 +0, –0, r3 s3 +0, T4 +0,002 +0, –0, r4 s4 +0, –0,000 –0, T –0, r5 +0,000 s –0,001 –0, T –0,001 –0, r6 s –0,001 –0, Р –0, r7 +0,002 s7 Р2 +0,001 +0, –0, r8 s8 +0,001 Р3 +0,000 +0, 1 +0,001 3 +0,000 –0,000 –0, Р –0,000 Р5 +0,002 +0, –0,001 –0, Р Также было проведено исследование расхождений между координатами пунктов, вычисленными по формулам (43) и (44).

= ;

Показано, что эти расхождения определяются как = ;

=.

(60) среднем 3,87 10 м (по X), 3,18 10 м (по Y) и 0,47 10 м (по Для сети, представленной на рис. 3, эти расхождения составили в Z). Для большинства случаев эти расхождения не только не превышают средние квадратические ошибки определения координат пунктов, но и оказываются существенно меньше ошибок округления. Следовательно, можно полагать, что результаты, получаемые по формулам (43) и (44) совпадают.

Также на примере сети, представленной на рис. 4, было произведено сравнение результатов уравнивания в пространстве (без учёта координаты Z) и на плоскости (аналогично уравниванию плановых сетей). Максимальное расхождение в полученных координатах составило 0.002 м (для X) и 0.001 м (для Y), что для большинства случаев оказывается несущественным. Таким образом, когда не требуется знать высоты пунктов, с точки зрения экономии процессорного времени компьютера и сокращения расчётов при объединении наземных и спутниковых геодезических сетей целесообразно отдавать предпочтение уравниванию на плоскости.

Заключение.

Основные результаты теоретических и экспериментальных исследований, выполненных в диссертационной работе, заключаются в следующем:

1. Предложена структура алгоритма совместной математической обработки геодезических сетей, построенных в разных системах координат, ход вычислений в которой может определяться автоматически, исходя из имеющихся данных.

2. Разработан алгоритм совместного уравнивания геодезических сетей, построенных в разных системах координат, на основе модифицированного способа условий с дополнительными неизвестными с применением формул рекуррентного уравнивания и метода сопряжённых градиентов, позволяющий выполнять контроль грубых ошибок и уточнять приближённые значения параметров преобразования и координаты общих пунктов сетей на стадии вычислений. Проведённые исследования показали эффективность использования данного алгоритма.

3. Показано, что использование исходной и линеаризованной формулы перехода из одной системы координат в другую приводит к совпадающим результатам и не оказывает существенного влияния на точность определения координат пунктов.

4. Показано, что в случае, когда нет необходимости знать высоты пунктов, при объединении наземных и спутниковых геодезических сетей целесообразно выполнять уравнивание на плоскости.

Список публикаций по теме диссертации:

1. Клыпин И.А. К вопросу об уравнивании геодезических сетей, построенных в разных системах координат // Известия вузов.

Геодезия и аэрофотосъёмка. – 2010. – №6. – С. 11–13.

2. Клыпин И.А. Объединение наземных геодезических сетей и спутниковых построений // Известия вузов. Геодезия и аэрофотосъёмка. – 2011. – №5. – С. 30-31.

3. Клыпин И.А., Коныжева М.В. Объединение нивелирных сетей, созданных в разных системах координат // Сборник статей по итогам научно-технических конференций: приложение к журналу «Известия вузов. Геодезия и аэрофотосъёмка». – Выпуск 3. – 2010. – С. 67–

 

Похожие работы:





 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.