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

Чернов В.В.

к.т.н., заместитель директора

Глухов О.В.

к.т.н., ведущий инженер отдела контроля качества

А.В. Кононенко

ведущий инженер отдела мониторинга земель

Восточно-Сибирский филиал ФГУП "Госземкадастрсъемка"-ВИСХАГИ

664007, Г. Иркутск ул. Софьи Перовской 30, оф. 301

e-mail:[e-mail]

Аннотация

Разработана специальная технология фотограмметрической съемки и обработки береговых уступов с борта движущегося судна, которая способна обеспечить мониторинг нарушенных земель береговой полосы водоемов метрическими (пространственными) данными с точностью масштаба 1:2000.

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

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

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

Нами была предложена технология стереофотограмметрической съемки участков береговой полосы, подверженной абразионным процессам, с борта специально оборудованного судна [1]. Для ее апробирования была использована следующая аппаратура (см. рис.1):

1.Цифровой фотоаппарат Canon 450D (с матрицей 12 Мп).

2.Комплект спутниковых двухчастотных GPS-приемников геодезического класса Trimble 5700 для получения координат центров фотографирования в процессе движения судна вдоль берега.

Рис.1. Комплект оборудования для стереофотосъемки

3.Портативный компьютер (ноутбук) для управления фотосъемкой с борта судна и сохранения цифровых фотоснимков на «жестком» диске в режиме временной синхронизации с бортовым GPS-приемником.

Для решения поставленной задачи – разработки технологии стереофотограмметрической съемки береговой полосы – было необходимо решить два основных вопроса:

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

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

Первая задача – калибровка цифровой камеры с целью определения элементов внутреннего ориентирования и параметров радиальной дисторсии – решена благодаря программному обеспечению, созданному в Лаборатории компьютерной графики и мультимедиа МГУ (http://graphics.cs.msu.ru). В программе реализована технология калибровки способом многократного фотографирования плоского тестового объекта с большими углами наклонов, с разных ракурсов. Такая методика калибровки была теоретически обоснована в монографии проф. Дубиновского В.Б. [2].

На примере, иллюстрированном рисунками 2–4, калибровка фотокамеры Canon 450D произведена по 27-ми цифровым фотоснимкам плоского тестового объекта в виде шахматки, размером 27х14 черно-белых квадратов со стороной 40 мм. Точность определения фокусного расстояния составила 0.3 pix, координат главной точки – 0.3 pix и 0.2 pix соответственно по осям х и у снимка.

Рис. 2. Автоматическое измерение тестового объекта

при калибровке цифровой фотокамеры

Рис. 3. Плотность измерений при калибровке цифровой фотокамеры

Рис. 4. Результаты калибровки цифровой фотокамеры

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

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

Для ее решения использованы строгие функциональные зависимости аналитической фотограмметрии - уравнения коллинеарности проектирующих лучей [3].

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

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

Для апробирования технологии стереосъемки с борта судна была разработана программа PicCalc, в которой реализована следующая схема обработки цифровых фотоснимков:

1. Производится измерение цифровых фотоснимков

Рис.6. Измерение цифровых снимков в программе PicCalc

2. Определяются приближенные значения угловых элементов внешнего ориентирования , ,  каждого снимка и координат X, Z точек уреза воды. Угловые элементы внешнего ориентирования снимков определяются через вычисленные азимуты при расчете траектории движения судна в программе GrafNav, а приближенные координаты X, Z точек уреза воды вычисляются по формулам прямой фотограмметрической засечки [4].

3. Приближенные значения угловых элементов внешнего ориентирования и координат точек уреза воды, а также известные линейные элементы внешнего ориентирования подставляются в уравнения (1) и вычисляются координаты x и y^ изображений точек на снимке. В силу погрешностей измерений, неточности приближенных величин, принятых в качестве неизвестных, и других причин, вычисленные координаты не будут равны измеренным значениям x и y.

4. Для каждой точки, измеренной на снимке, составляется пара уравнений поправок, принимая разности l_x, l_y между измеренными и вычисленными координатами точек на снимке в качестве их свободных членов

5. Система линейных уравнений (7) решается методом наименьших квадратов

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

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

7. Выполняется оценка точности решения задачи. Для этого из последнего приближения используются весовые коэффициенты Q (диагональные элементы матрицы ) и поправки V.

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

С целью проверки работоспособности алгоритма обработки, а так же для определения оптимальных параметров съемки в марте 2010 года было произведено маркирование экспериментального полигона (на льду замерзшего карьера), определение координат и высот марок с точностью 0.03 м электронным тахеометром, непосредственно стереофотосъемка с точек фотографирования и обработка данных.

В одном проекте PicCalc в обработке одновременно использовалось три снимка (триплет). Предварительно на снимках была выполнена коррекция систематических искажений, вызванная радиальной дисторсией, на основе выполненной ранее калибровки фотокамеры. Измерения снимков выполнены в программе PicCalc

Рис.7. Измерение снимков экспериментального полигона.

Результаты уравнивания триплета ф8-ф2-ф4 по 15-ти опорным точкам и центрам фотографирования представлены в таблице 1. Расстояние между точками фотографирования составляло около 50 м (базисы фотографирования), среднее отстояние базисов от определяемых точек - 190 м, снимок на точке фотографирования ф8 был сделан с разворотом фотокамеры примерно на 15 градусов.

Таблица 1

Вычисленные угловые элементы внешнего ориентирования по опорным точкам и центрам фотографирования

	снимок с т.ф4	снимок с т.ф2	снимок с т.ф8
альфа, градусы	4.90	2.98	-14.82
омега, градусы	-0.95	-0.93	-0.65
каппа, градусы	-1.35	-1.07	-0.96

Результаты уравнивания триплета ф8-ф2-ф4 по 15-ти точкам «уреза воды» представлены в таблице 2, а разница угловых параметров, вычисленных по опорным точкам и только по точкам уреза воды – в таблице 3.

Таблица 2

Вычисленные угловые элементы внешнего ориентирования по точкам «уреза воды» и центрам фотографирования

	снимок с т.ф4	снимок с т.ф2	снимок с т.ф8
альфа, градусы	4.95	2.84	-15.16
омега, градусы	-0.94	-0.94	-0.64
каппа, градусы	-1.32	-1.07	-1.00

Таблица 3

Разница вычисленных угловых элементов внешнего ориентирования по опорным точкам и по точкам «уреза воды»

	снимок с т.ф4	снимок с т.ф2	снимок с т.ф8
альфа, градусы	0.05	-0.14	-0.34
омега, градусы	0.00	0.00	0.00
каппа, градусы	0.03	0.00	-0.04

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

При уравнивании триплета ф8-ф2-ф4 по точкам «уреза воды» и центрам фотографирования были получены следующие точностные показатели:

Таблица 4