Отопление, канализация, водоснабжение

Использование аэрофотосъемки с бла в различных. Топографическая съёмка с использованием бпла

Первая часть статьи «БЕСПИЛОТНЫЕ ЛЕТАТЕЛЬНЫЕ АППАРАТЫ: ПРИМЕНЕНИЕ В ЦЕЛЯХ АЭРОФОТОСЪЕМКИ ДЛЯ КАРТОГРАФИРОВАНИЯ» касалась вопросов общейтеории: были рассмотрены существующие типы БПЛА, приведены пояснения основных терминов, связанных с их использованием, а также дан обзор нескольких моделей БПЛА, успешно применяемых при аэрофотосъемке в картографических целях.

Во второй части статьи будут рассмотрены особенности фотограмметрической обработки беспилотной аэросъемки, даны рекомендации по ее проведению и по установке основного и дополнительного оборудования на борт БПЛА для получения максимальной точности.

А.Ю. Сечин, М.А. Дракин, А.С. Киселева, «Ракурс», Москва, Россия, 2011.

Особенности данных аэросъемки с БПЛА

Принципиально не отличается от съемки с «больших самолетов», но имеет определенные особенности, которые мы далее рассмотрим. Полет БПЛА, как правило, производится с крейсерской скоростью 70-110 км/ч (20-30 м/c) в диапазоне высот 300-1500 м. Для съемки обычно используются неметрические бытовые камеры с размером матрицы 10-20 мегапикселей. Фокусное расстояние камер обычно составляет около 50 мм (в 35 мм эквиваленте), что соответствует размеру пикселя на местности (GSD) от 7 до 35 см.

Часто снимки с БПЛА обрабатываются простыми нестрогими методами (аффинное преобразование снимков на плоскость). В результате, пользователь получает накидные монтажи, которые помимо низкой точности могут содержать разрывы контуров на стыках соседних снимков.

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

Для строгой фотограмметрической обработки данных аэросъемки и получения максимально точных результатов необходимо, чтобы снимки в одном маршруте имели тройное перекрытие, а перекрытие между снимками соседних маршрутов при площадной съемке составляло не менее 20%. На практике, при съемке с БПЛА эти параметры выдерживаются далеко не всегда. Полет БПЛА не устойчив, на него влияют порывы ветра, турбулентность и другие возмущающие факторы. Если съемку с обычных самолетов планируют с перекрытием вдоль маршрута 60%, а между маршрутами 20-30%, то проектировать съемку с БПЛА следует с перекрытием вдоль маршрутов 80%, а между маршрутами – 40%, чтобы, по возможности, исключить разрывы в фототриангуляционном блоке .

На БПЛА, как правило, устанавливаются цифровые камеры Canon. Это связано с легкостью электронного управления камерами этой фирмы. Использование бытовых камер имеет как преимущества (невысокая стоимость, легкость замены при «жесткой посадке»), так и недостатки.

Основным недостатком является то, что бытовые камеры изначально не откалиброваны – неизвестны их точные фокусные расстояния, главная точка, дисторсия. При этом нелинейные искажения оптики (дисторсия), допустимые при бытовой съемке, могут составлять до нескольких десятков пикселей, что на порядок снижает точность результатов обработки. Однако, такие камеры могут быть откалиброваны в лабораторных условиях, что позволяет получать точности обработки, практически такие же, как и для профессиональных малоформатных фотограмметрических камер.

Предпочтительней устанавливать на такие камеры объективы с фиксированным фокусным расстоянием. При съемке следует выставлять фокусировку на бесконечность и отключать функцию «автофокуса».

Второй недостаток используемых на БПЛА камер относится конкретно к камерам Canon– в них, в отличие от профессиональных фотограмметрических камер, используется щелевой затвор, в результате чего экспозиция разных частей изображения производится в разные моменты времени и соответствует разным положениям носителя. Так, если выдержка при съемке составляет 1/250 c, то при скорости БПЛА в 20 м/с смещение камеры при съемке кадра составляет 8 см, что сравнимо с разрешением съемки на малых высотах и вызывает дополнительную систематическую ошибку в снимке. Такие ошибки могут накапливаться в процессе фотограмметрического сгущения (уравнивании) при съемке протяженных территорий. Для того, чтобы уменьшить влияние этого эффекта и для ликвидации «смаза» снимков, следует осуществлять съемку с БПЛА с наименьшими возможными выдержками (не длиннее 1/250 c, максимальная выдержка зависит от высоты). Частично проблему щелевого затвора могли бы решить камеры с центральным затвором, имеющие сравнимое с камерами Canonкачество объектива и матриц. Тем не менее, чтобы избежать «смаза» выдержки все равно следует ограничивать.

Снимки цифровых камер, как любительских, так и профессиональных, имеют прямоугольную форму. «Выгоднее» располагать камеру так, чтобы длинная сторона снимка располагалась поперек полета – это позволяет снимать большую площадь при той же длине маршрута. Съемку следует производить с максимальным качеством – с наименьшим jpegсжатием или в RAW, если последнее возможно.

Современный уровень развития навигационных средств позволяет производить измерения элементов внешнего ориентирования (ЭВО) непосредственно в процессе съемки. Типичные точности таких измерений достигают единиц сантиметров по пространственным координатам X,Yи Zи 0.005 градуса по углам крена, тангажа и рысканья для самых точных систем ApplanixPOSAV, устанавливаемых на «большие самолеты». Часто этого достаточно, чтобы производить обработку без использования опорных точек. В любом случае, наличие таких данных значительно упрощает обработку и позволяет выполнять некоторые этапы обработки полностью в автоматическом режиме. Современные достижения микроэлектроники позволяют собрать механический (точнее MEMS– электронно-механический) гироскоп в корпусе размером в несколько мм, стоимостью от 250 $. Такие гироскопы не дают точность профессиональных, имеют значительный уход (порядка одного градуса за час) при эксплуатации, но существенно упрощают последующую обработку данных. При типовых поставках , Дозор 50 на борт могут быть установлены такие малогабаритные инерциальные системы - IMU(на Дозор-50 ставится IMUразработки ООО «Транзаз Телематика») и высокоточные двухдиапазонные GPS (TOPCONeuro160 на Птеро-E4, встроенный ГЛОНАСС/GPS приемник на Дозор-50). Паспортная точность этих GPS приборов составляет 10 мм + 1,5 мм × B (B – удаление до базовой станции в км) в плане и 20 мм + 1,5 мм × B по высоте. К сожалению, обычно на борт БПЛА устанавливают болеедешевые GPSприемники и не устанавливают IMU датчики. Данные о центрах проекции снимков в телеметрической информации снимаются через протокол NMEAи имеют в таком случае точность до 20-30 м, а углы тангажа, крена и рысканья вычисляются через вектор скорости GPSизмерений. Точность угла рысканья в такой телеметрической информации невысокая и может превышать 10 градусов, а сами значения содержат систематические ошибки, что усложняет последующую обработку данных.

Если при съемке использовался двухдиапазонный GPSприемник в дифференциальном режиме (или PPPобработка данных GPS), то требуется минимальное число опорных точек для получения наиболее точных результатов обработки, обычно достаточно 1-2 точки на 100 снимков, в ряде случаев обработку можно проводить без опорных точек. В случае, когда нет точных центров проекции, требования к планово-высотному обоснованию стандартные: одна планово-высотная точка на 6-10 базисов съемки.

Специфика фотограмметрической обработки данных аэросъемки с БПЛА

Обработка аэрофотосъемки с БПЛА в цифровых фотограмметрических системах (ЦФС) в целом аналогична обработке аэрофотосъемки с «больших самолетов». Однако особенности данных с борта БПЛА часто не позволяют использовать автоматические процедуры стандартных пакетов – часть операций (например, расстановку связующих точек) приходится производить в ручном режиме. Ниже мы рассмотрим особенности обработки аэросъемки с БПЛА в ЦФС PHOTOMOD5.2. Именно в этой версии PHOTOMOD введены специальные функции для обработки таких данных, существенно упрощающие и автоматизирующие получение конечной продукции.

Как и при обработке других данных, сначала в ЦФС создается проект, в него вводятся снимки и телеметрическая информация. На основании данных о центрах проекции и углах производится создание накидного монтажа, разбивка по маршрутам. Снимки, попавшие на развороты БПЛА, удаляются в ручном режиме. Неточные угловые элементы внешнего ориентирования приводят к достаточно грубому накидному монтажу (Рис. 1):

Рис. 1. Накидной монтаж по телеметрической информации

Автоматический поиск связующих точек в таких случаях затруднен или требует значительного времени работы компьютера. Для уточнения накидного монтажа в таких случаях в ЦФС PHOTOMOD используется т.н. «автоматический накидной монтаж», который уточняет взаимное расположение снимков (Рис. 2).

Рис. 2. Накидной монтаж после автоматического уточнения

Как мы ранее отмечали, съемка с борта БПЛА производится с увеличенными перекрытиями. Нестабильность полета летательного аппарата иногда может привести к очень большим перекрытиям между соседними снимками, что вызывает сложности в стандартных фотограмметрических пакетах.

Рис. 3. «Перепутывание» снимков при маленьком базисе съемки

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

Рис. 4. Накидной монтаж после первых проходов автоматического измерения связующих точек

На следующих проходах производится доизмерение связующих точек. Несколько проходов необходимы в случае, когда телеметрическая информация не содержит всех углов ориентирования, или углы известны с точностью 10-30 градусов. Если же телеметрическая информация содержит угловые элементы ориентирования с точностью в несколько единиц градуса, то достаточно и одного прохода – надежность автоматических измерений в этом случае повышается. Для борьбы с возможными грубыми ошибками при автоматических измерениях в PHOTOMOD5.2 введено понятие т.н. «доверительной группы связующих точек», когда программа ищет наибольшее число связующих точек для стереопар с наименьшим поперечным параллаксом, остальные связующие точки, не попавшие в группу, считаются ошибочными.

После измерения связующих и опорных точек производится процедура уравнивания. В ЦФС PHOTOMODможно использовать начальное приближение для алгоритма уравнивания как по уточненной схеме блока, так и построенное другими методами. Начиная с версии 5.2 для уравнивания аэросъемки с БПЛА мы рекомендуем использовать новый режим – уравнивание 3D. При уравнивании в PHOTOMODи достаточном числе опорных точек можно использовать самокалибровку. Это дает возможность использования некалиброванных камер. Ожидаемая точность выходных результатов при строгой фотограмметрической обработке составляет приблизительно 1-2 GSDв плане и 2-4 GSDпо высоте. После фотограмметрического уравнивания, результаты которого и определяют точность выходных продуктов, производится построение рельефа (ЦМР) в автоматическом режиме. При необходимости, после уравнивания может быть сделана стереовекторизация – отрисовка в ручном режиме зданий, сооружений, мостов, дамб и других объектов. Построенный рельеф используется для ортотрансформирования снимков. На последнем этапе из ортотрансформированных снимков создается бесшовная мозаика – производится расчет линий порезов, выравнивание яркостей, стыковка контурных объектов. Самокалибровку можно включать и при отсутствии опорных точек, правда, в этом случае можно рассчитать только коэффициенты k 1, k 2 радиальной дисторсии. При использовании камер с щелевым затвором можно дополнительно включить расчет аффинных искажений. В случае стабильности углов ориентирования при съемке такая самокалибровка может повысить точность уравнивания.

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

GPSцентров проекций и дисторсией объектива (до нескольких десятков пикселей). В таких случаях можно применять упрощенную автоматизированную последовательность обработки. Бесшовный накидной монтаж указанной точности при этом получается за счет трансформирования исходных снимков в модуле PHOTOMODGeoMosaic.В этом случае используются простейшие методы трансформирования, не учитывающие рельеф местности, а стыковка контуров осуществляется за счет автоматически рассчитываемых связующих точек вдоль автоматически построенных линий порезов.

Примеры фотограмметрической обработки данных аэросъемки с БПЛА

Рассмотрим несколько примеров . Во всех примерах для обработки использовалась ЦФС PHOTOMOD. Отметим, что различными организациями в компанию «Ракурс» для тестирования было передано более 20 блоков аэросъемки с БПЛА. К сожалению, для многих блоков отсутствовали опорные точки и/или съемка была проведена неоткалиброванными камерами. В таких случаях было невозможно оценить точность конечных результатов обработки.

Первый блок, который мы рассмотрим, был снят с борта БПЛА ZALA421-04ф. Данные для исследований были любезно предоставлены ОАО «Газпром космические системы». Блок состоял из 26 маршрутов. Общее число снимков в блоке составило 595. Использовалась предварительно откалиброванная цифровая камера Canon EOS500D. Высота залета над местностью составила около 500 м, размер пиксела на местности приблизительно равен 8 см. На местности были измерены и промаркированы 25 опорных точек, точность координат опорных точек не превышала 10 см. Общий перепад высот местности протяженностью около 3-х километров достаточно большой ~ 70 метров.

Сначала этот же блок аэросъемки был обработан в автоматическом режиме по упрощенной схеме, без уравнивания и использования опорных точек. Привязка осуществлялась по центрам проекции, трансформирование снимков проводилось сразу в модуле GeoMosaicбез учета рельефа. Последующий контроль полученных «псевдо» ортофотопланов по опорным точкам показал расхождения на опорных точках, превышающие 17 м. Такая невысокая точность ортофотплана обусловлена как большим перепадом высот, так и неточностью измерений центров проекций в полете.

Затем блок был подвергнут строгой фотограмметрической обработке. При уравнивании три из измеренных опорных точек считались контрольными. Среднеквадратическая ошибка уравнивания составила по опорным точкам 15 см, 16 см, 12 см, по контрольным точкам 23 см, 29 см и 57 см. Расхождения на связующих точках составили 8 см, 14 см и 69 см. Общий вид блока представлен на следующем рисунке.

Рис. 5. Общий вид «блока 1»

В процессе уравнивания было обнаружено, что координаты центров проекций из телеметрической информации содержат систематическую ошибку, главная из компонент которой составляет 10,5 метра по высоте Z. Среднеквадратические ошибки на центрах проекции после вычитания систематической ошибки составили 84 см, 239 см и 75 см. Существенно большая ошибка по Y(вдоль полета), скорее всего, связана с неточным определением моментов съемки в телеметрии. Большие ошибки по Zна связующих точках возможно связаны с неточной калибровкой камеры и с накопленной ошибкой при съемке камерой с щелевым затвором. Наибольшие ошибки на связующих точках наблюдаются на краях и в углах снимков.

Рис. 6. Величины ошибок на связующих точках

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

Рис. 7. Фрагмент ортофотоплана без выравнивания яркости

В апреле 2011 кафедрой фотограмметрии Московского государственного университета геодезии и картографии (МИИГАиК) были проведены исследования материалов аэрофотосъёмки, полученных с помощью БПЛА Птеро, с целью оценки качества аэросъёмочных работ и фотограмметрической обработки . Съемка выполнялась с высоты около 900 м над средней плоскостью снимаемой местности с борта БПЛА Птеро цифровой фотокомерой CanonEOS5D. Камера была предварительно откалибрована. Для оценки качества материалов использовался фрагмент блока, состоящий из 2-х маршрутов по 6 снимков в каждом. В качестве опорных использовались 14 точек, плановые координаты XYкоторых были сняты с планов масштаба 1:1000, а высота Zопределялась по материалам воздушного лазерного сканирования, выполненного с точностью около 20-30 см. После фотограмметрического уравнивания среднеквадратические погрешности координат на опорных точках составили по X, Yи Zсоответственно 20 см, 21 см и 50 см. Среднеквадратические погрешности координат связующих точек составили 6 см, 6 см, 15 см. Размер пиксела на местности для этого блока GSDсоставляет около 12 см. Общая схема блока показана на следующем рисунке.

Рис. 8. Схема «блока 2» с опорными и связующими точками

Вопросы метрологического обеспечения

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

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

Отметим здесь интересный проект МИИГАиК – с целью отработки и исследования технологий мониторинга и картографирования местности по материалам беспилотной аэрофотосъемки, начаты работы по созданию специализированного исследовательского полигона. Этот полигон, площадью около 50 кв. км, создается в Заокском районе Тульской области, на базе учебного геополигона МИИГАиК, расположенного в 110 км от Москвы.

Территория полигона представляет собой уникальное многообразие картографических объектов. На этой территории расположены разнообразные населенные пункты: поселок городского типа, деревни, дачные и коттеджные поселки; дорожная сеть в виде железных, шоссейных, проселочных и полевых дорог; линии электропередачи различного напряжения; трубопроводы. На территории полигона имеются лесные массивы, различные гидрографические объекты, многообразные формы рельефа, сельскохозяйственные угодья и производственные объекты.

С целью обеспечения отработки и исследования технологий, основанных на применении БПЛА, на территории полигона начаты работы по созданию высокоточной сети планово-высотных опознаков (в виде естественных контуров местности и маркировочных знаков); ведется топографическая наземная съемка характерных участков местности в масштабе 1: 500 и 1: 2000. На эту же территорию по материалам аэрофотосъемки и космическим снимкам высокого разрешения созданы ортофотопланы и цифровые модели местности. По мере поступления новых съемочных материалов эти работы предполагается выполнять в дежурном режиме.

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

Первые испытания планируется провести в середине июля 2011 г. Планируется провести тестовую аэрофотосъемку территории полигона в различных масштабах с помощью отечественного БПЛА «ПТЕРО» с целью отработки и исследования фотограмметрической технологии создания карт различного масштаба по полученным материалам аэрофотосъемки. Фотограмметрическую обработку полученных снимков предполагается выполнить на цифровой фотограмметрической системе PHOTOMOD. В сентябре предполагается провести испытания и БПЛА «МИИГАиК Х8», разработанного в МИИГАиК.

Созданием полигона и проведением на нем испытаний БПЛА и технологий, основанных на их использовании, МИИГАиК намерен помочь потенциальным пользователям освоить и внедрить новые технологии, а разработчикам летательных аппаратов и съемочных систем адаптировать их к решению актуальных задач производства.

Выводы

Использование БПЛА в качестве аэросъемочной платформы имеет большие перспективы при съемке небольших по протяженности площадных объектов и при съемке линейных объектов. Данные с БПЛА позволяют получать качественные картографические материалы (пространственные данные) при следующих условиях:

    выполнении определенных (вполне посильных) требований к съемочной аппаратуре и процессу съемки (гарантия достаточности перекрытий);

    строгой фотограмметрической обработке. Точность при этом возрастает в десятки раз и может составлять около GSD, как и для обычной аэросъемки и космических снимков.

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

    Использовать на БПЛА калиброванные камеры.

    Производить съемку с выдержкой не длиннее 1/250с.

    Использовать объективы с фиксированным фокусным расстоянием. Если это невозможно, следует фиксировать увеличение (Zoom). Съемка должна производиться с фокусировкой на бесконечность и с отключенным режимом автофокусировки.

    Проектировать съемку с увеличенными перекрытиями (80% вдоль, 40% поперек маршрута).

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

    Желательно использовать двухдиапазонные GPS приемники на борту и дифференциальный режим измерений.

    Желательно использование на борту IMU, пусть и не имеющего высокой точности.

Благодарности

Благодарим компании: «Беспилотные системы ЗАЛА АЭРО», ОАО «Газпром космические системы», «АФМ-Серверс», ООО «Геометр-Центр», НПИ и КЦ «Земинформ», ЗАО «Транзас», ЗАО «Лимб» за помощь в подготовке материала, предоставление данных и полезные обсуждения.

Литература

    Чибуничев А.Г., Михайлов А.П., Говоров А.В. Калибровка цифровых фотокамер: Вторая научно-практическая конференция РОФДЗ. Тезисы докладов. М., 2001 г. с38-39.

    Скубиев С.И., Научно-производственный институт земельно-информационных технологий Государственного университета по землеустройству «Земинформ» (Россия), Использование беспилотных летательных аппаратов для целей картографии. Тезисы XЮбилейной международной научно-технической конференции «От снимка к карте: цифровые фотограмметрические технологии». Гаета, Италия, 2010.

    Результаты полевых исследований БПЛА «Птеро»

Приводим приблизительные расценки на аэросъёмочные работы, осуществляемые при помощи БПЛА.
Расценки могут меняться в зависимости от ПЛОЩАДИ СЪЁМКИ , ТРАНСПОРТНЫХ РАСХОДОВ , ВРЕМЕНИ ГОДА (наличия листвы/снега). Наилучшее соотношение цены/качества/скорости получается в промежутке между таянием снега и появлением листвы на деревьях.
Минимальная площадь съёмки одного объекта 6 км 2 , если есть несколько объектов с расстояниями между ними до 30 км, то мин объём 4 км 2 .

Цена аэрофотосъёмки БПЛА

Цены на аэрофотосъёмку БПЛА приведены на 1 км 2 .
1. Ортофотоплан в WGS84
Масштаб 1:500 (4 см./пикс) – 35 000 рублей 1 км 2 .
Масштаб 1:1000 (7 см./пикс) – 22 000 рублей км 2 .
Масштаб 1:2000 (10 см./пикс) – 17 000 рублей км 2 .
Масштаб 1:5000 (15 см./пикс) – 12 000 рублей км 2 .

2. Создание ортофотоплана с привязкой к местным системам координат:
Масштаб 1:500 + 10 000 рублей км 2 .
Масштаб 1:1000 + 6 000 рублей. км 2 .
Масштаб 1:2000 + 4 500 рублей км 2 .
Масштаб 1:5000 + 3 000 рублей км 2 .

3. Построение ЦМР или горизонталей:
Сечение 0,5 м + 12 000 рублей за км 2 .
Сечение 1 м + 8 000 рублей за км 2 .
Сечение 2 м + 5 000 рублей за км 2 .

4. Рассекречивание материалов ДЗЗ – 40 000 рублей за объект.
5. Отрисовка топографического плана по результатам съёмки БПЛА: оценивается индивидуально – от 5 000 рублей за км 2 .
Таким образом, стоимость 1 км 2 плана масштаба 1:2 000 при выходе на бумагу будет стоить 34 500, т.е. по 345 р/га – такую цену невозможно получить никаким другим способом!

Для больших площадей съёмки в более крупных масштабах (до 1:500) нами разработан способ комбинированных работ, включающих как съемку при помощи БПЛА, так и тахеометрию с проложением магистрального хода и подсъёмкой основных элементов местности.

Качество не раз проверялось инструментально с земли, в том числе силами заказчиков.

Преимущества технологии аэрофотосъёмки беспилотными летательными аппаратами.

Технология аэрофотосъёмки беспилотными летательными аппаратами развивается уже много лет но только в последнее время подошла по своим точностным характеристикам к классическим способам геодезической съемки и на уровне масштабного ряда от 1:500 и мельче сравнялась с ними. На данный момент АФС БПЛА находится на переднем крае развития геодезических технологий, вытесняя в обширном сегменте такие методы как тахеометрия, спутниковое позиционирование в режимах RTK, пилотируемая АФС, воздушное лазерное сканирование, делая их неконкурентоспособными как по стоимости, так и по срокам.
При больших объемах, слабо залесенной и слабо застроенной площади аэрофотосъемка БПЛА делает нерентабельными тахеометрическую и GPS съемку уже при площади в 70 га. В то же время воздушное лазерное сканирование и пилотируемая аэрофотосъёмка могут конкурировать с БПЛА лишь при объемах от 50 кв. км площадных объектов или от 200 пг км линейных. Итак, на данный момент развития технологий АФС БПЛА достаточно вольготно себе чувствует на открытых территориях в диапазоне объема работ от 0,7 до 50 кв. км.

К недостаткам АФС БПЛА стоит отнести ее метеозависимость и сезонозависимость (снежный покров или наличие очень густой растительности чаще всего препятствует получению достаточного для построения инженерно-топографического плана материала). Следует отметить, что эти факторы абсолютно в той же мере препятствует и другим способам съёмки. В масштабном ряду съёмка ограничена диапазоном от 1:5000 до 1:1000 (мельче 1:5000 целесообразнее использовать космические снимки, крупнее 1:1000 необходимы комбинированные с наземными средствами методы).
Растительность также может влиять на результат. В нашей практике мы сталкивались с бамбуковыми зарослями на южных Курильских островах, которые оказались не просвечиваемы практически на 100%, то же касается тростника заболоченных участков Юга России (например, дельта Волги) и, по всей видимости, тропической растительности юга Черноморского побережья. Лесные массивы средней и северной части страны, а также Сибири и Дальнего Востока, как правило, не являются помехой для АФС БПЛА.
Плотная городская застройка может накладывать ограничения на сам процесс полёта, а также скрывать за своими структурами множество элементов, не видимых сверху. Тем не менее, для таких объектов как СНТ, АФС БПЛА становится незаменимым решением ввиду ограничения доступа на каждый участок, значительно ускоряя процесс съёмки.

Применение аэрофотосъёмки БПЛА

Кроме топографической съёмки беспилотные технологии применяются нами и для различных форм мониторинга, объектами которого могут выступать несанкционированные свалки твердых бытовых отходов, линейные объекты - ЛЭП, трубопроводы, транспортная инфраструктура. Также БПЛА решает вопросы определения объемов земляных масс и их динамики при разработках месторождений открытым способом, карьеров. По сравнению с космическим мониторингов БПЛА дает несоизмеримо более актуальную информацию - при заказе космического снимка на конкретную территорию вы можете либо воспользоваться снимком их архива 1-3 месячной, а то и годовой давности или наоборот, ждать долгое время подходящей телеметрии спутника и своей очереди на съёмку вашей территории – данные с беспилотника вы можете получить в работу через несколько часов.
Аграрный сектор в последнее время становится одним из основных потребителей технологии. Агрохолдинги и крупные сельскохозяйственные объединения, обладающие большими площадями незастроенной и открытой местности, заинтересованы не только в производстве инженерно-геодезических изысканий для реконструкции и нового строительства, но и в мониторинге, инвентаризации сельхозугодий, оценки всхожести культур, прогнозе урожайности, мониторинге эрозионных процессов. Нами используется нормализованный вегетационный индекс (NDVI), позволяющий на основе разности ближней инфракрасной и красной областях спектра определять фотосинтезирующую активность биомассы.

Пространственное разрешение аэрофотоснимков с БПЛА

Сравнение снимка Google и аэрофотоснимка с БПЛА

Первая часть статьи «БЕСПИЛОТНЫЕ ЛЕТАТЕЛЬНЫЕ АППАРАТЫ: ПРИМЕНЕНИЕ В ЦЕЛЯХ АЭРОФОТОСЪЕМКИ ДЛЯ КАРТОГРАФИРОВАНИЯ» касалась вопросов общейтеории: были рассмотрены существующие типы БПЛА, приведены пояснения основных терминов, связанных с их использованием, а также дан обзор нескольких моделей БПЛА, успешно применяемых при аэрофотосъемке в картографических целях.

Во второй части статьи будут рассмотрены особенности фотограмметрической обработки беспилотной аэросъемки, даны рекомендации по ее проведению и по установке основного и дополнительного оборудования на борт БПЛА для получения максимальной точности.

А.Ю. Сечин, М.А. Дракин, А.С. Киселева, «Ракурс», Москва, Россия, 2011.

Особенности данных аэросъемки с БПЛА

Аэрофотосъемка с БПЛА принципиально не отличается от съемки с «больших самолетов», но имеет определенные особенности, которые мы далее рассмотрим. Полет БПЛА, как правило, производится с крейсерской скоростью 70-110 км/ч (20-30 м/c) в диапазоне высот 300-1500 м. Для съемки обычно используются неметрические бытовые камеры с размером матрицы 10-20 мегапикселей. Фокусное расстояние камер обычно составляет около 50 мм (в 35 мм эквиваленте), что соответствует размеру пикселя на местности (GSD) от 7 до 35 см.

Часто снимки с БПЛА обрабатываются простыми нестрогими методами (аффинное преобразование снимков на плоскость). В результате, пользователь получает накидные монтажи, которые помимо низкой точности могут содержать разрывы контуров на стыках соседних снимков.

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

Для строгой фотограмметрической обработки данных аэросъемки и получения максимально точных результатов необходимо, чтобы снимки в одном маршруте имели тройное перекрытие, а перекрытие между снимками соседних маршрутов при площадной съемке составляло не менее 20%. На практике, при съемке с БПЛА эти параметры выдерживаются далеко не всегда. Полет БПЛА не устойчив, на него влияют порывы ветра, турбулентность и другие возмущающие факторы. Если съемку с обычных самолетов планируют с перекрытием вдоль маршрута 60%, а между маршрутами 20-30%, то проектировать съемку с БПЛА следует с перекрытием вдоль маршрутов 80%, а между маршрутами – 40%, чтобы, по возможности, исключить разрывы в фототриангуляционном блоке .

На БПЛА, как правило, устанавливаются цифровые камеры Canon. Это связано с легкостью электронного управления камерами этой фирмы. Использование бытовых камер имеет как преимущества (невысокая стоимость, легкость замены при «жесткой посадке»), так и недостатки.

Основным недостатком является то, что бытовые камеры изначально не откалиброваны – неизвестны их точные фокусные расстояния, главная точка, дисторсия. При этом нелинейные искажения оптики (дисторсия), допустимые при бытовой съемке, могут составлять до нескольких десятков пикселей, что на порядок снижает точность результатов обработки. Однако, такие камеры могут быть откалиброваны в лабораторных условиях, что позволяет получать точности обработки, практически такие же, как и для профессиональных малоформатных фотограмметрических камер.

Предпочтительней устанавливать на такие камеры объективы с фиксированным фокусным расстоянием. При съемке следует выставлять фокусировку на бесконечность и отключать функцию «автофокуса».

Второй недостаток используемых на БПЛА камер относится конкретно к камерам Canon– в них, в отличие от профессиональных фотограмметрических камер, используется щелевой затвор, в результате чего экспозиция разных частей изображения производится в разные моменты времени и соответствует разным положениям носителя. Так, если выдержка при съемке составляет 1/250 c, то при скорости БПЛА в 20 м/с смещение камеры при съемке кадра составляет 8 см, что сравнимо с разрешением съемки на малых высотах и вызывает дополнительную систематическую ошибку в снимке. Такие ошибки могут накапливаться в процессе фотограмметрического сгущения (уравнивании) при съемке протяженных территорий. Для того, чтобы уменьшить влияние этого эффекта и для ликвидации «смаза» снимков, следует осуществлять съемку с БПЛА с наименьшими возможными выдержками (не длиннее 1/250 c, максимальная выдержка зависит от высоты). Частично проблему щелевого затвора могли бы решить камеры с центральным затвором, имеющие сравнимое с камерами Canonкачество объектива и матриц. Тем не менее, чтобы избежать «смаза» выдержки все равно следует ограничивать.

Снимки цифровых камер, как любительских, так и профессиональных, имеют прямоугольную форму. «Выгоднее» располагать камеру так, чтобы длинная сторона снимка располагалась поперек полета – это позволяет снимать большую площадь при той же длине маршрута. Съемку следует производить с максимальным качеством – с наименьшим jpegсжатием или в RAW, если последнее возможно.

Современный уровень развития навигационных средств позволяет производить измерения элементов внешнего ориентирования (ЭВО) непосредственно в процессе съемки. Типичные точности таких измерений достигают единиц сантиметров по пространственным координатам X,Yи Zи 0.005 градуса по углам крена, тангажа и рысканья для самых точных систем ApplanixPOSAV, устанавливаемых на «большие самолеты». Часто этого достаточно, чтобы производить обработку без использования опорных точек. В любом случае, наличие таких данных значительно упрощает обработку и позволяет выполнять некоторые этапы обработки полностью в автоматическом режиме. Современные достижения микроэлектроники позволяют собрать механический (точнее MEMS– электронно-механический) гироскоп в корпусе размером в несколько мм, стоимостью от 250 $. Такие гироскопы не дают точность профессиональных, имеют значительный уход (порядка одного градуса за час) при эксплуатации, но существенно упрощают последующую обработку данных. При типовых поставках Птеро E4, Дозор 50 на борт могут быть установлены такие малогабаритные инерциальные системы - IMU(на Дозор-50 ставится IMUразработки ООО

«Транзаз Телематика») и высокоточные двухдиапазонные GPS (TOPCONeuro160 на Птеро-E4, встроенный ГЛОНАСС/GPS приемник на Дозор-50). Паспортная точность этих GPS приборов составляет 10 мм + 1,5 мм × B(B– удаление до базовой станции в км) в плане и 20 мм + 1,5 мм × Bпо высоте. К сожалению, обычно на борт БПЛА устанавливают болеедешевые GPSприемники и не устанавливают IMU датчики. Данные о центрах проекции снимков в телеметрической информации снимаются через протокол NMEAи имеют в таком случае точность до 20-30 м, а углы тангажа, крена и рысканья вычисляются через вектор скорости GPSизмерений. Точность угла рысканья в такой телеметрической информации невысокая и может превышать 10 градусов, а сами значения содержат систематические ошибки, что усложняет последующую обработку данных.

Если при съемке использовался двухдиапазонный GPSприемник в дифференциальном режиме (или PPPобработка данных GPS), то требуется минимальное число опорных точек для получения наиболее точных результатов обработки, обычно достаточно 1-2 точки на 100 снимков, в ряде случаев обработку можно проводить без опорных точек. В случае, когда нет точных центров проекции, требования к планово-высотному обоснованию стандартные: одна планово-высотная точка на 6-10 базисов съемки.

Специфика фотограмметрической обработки данных аэросъемки с БПЛА

Обработка аэрофотосъемки с БПЛА в цифровых фотограмметрических системах (ЦФС) в целом аналогична обработке аэрофотосъемки с «больших самолетов». Однако особенности данных с борта БПЛА часто не позволяют использовать автоматические процедуры стандартных пакетов – часть операций (например, расстановку связующих точек) приходится производить в ручном режиме. Ниже мы рассмотрим особенности обработки аэросъемки с БПЛА в ЦФС PHOTOMOD5.2. Именно в этой версии PHOTOMOD введены специальные функции для обработки таких данных, существенно упрощающие и автоматизирующие получение конечной продукции.

Как и при обработке других данных, сначала в ЦФС создается проект, в него вводятся снимки и телеметрическая информация. На основании данных о центрах проекции и углах производится создание накидного монтажа, разбивка по маршрутам. Снимки, попавшие на развороты БПЛА, удаляются в ручном режиме. Неточные угловые элементы внешнего ориентирования приводят к достаточно грубому накидному монтажу (Рис. 1):

Рис. 1. Накидной монтаж по телеметрической информации

Автоматический поиск связующих точек в таких случаях затруднен или требует значительного времени работы компьютера. Для уточнения накидного монтажа в таких случаях в ЦФС PHOTOMOD используется т.н. «автоматический накидной монтаж», который уточняет взаимное расположение снимков (Рис. 2).

Рис. 2. Накидной монтаж после автоматического уточнения

Как мы ранее отмечали, съемка с борта БПЛА производится с увеличенными перекрытиями. Нестабильность полета летательного аппарата иногда может привести к очень большим перекрытиям между соседними снимками, что вызывает сложности в стандартных фотограмметрических пакетах.

Рис. 3. «Перепутывание» снимков при маленьком базисе съемки

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

Рис. 4. Накидной монтаж после первых проходов автоматического измерения связующих точек

На следующих проходах производится доизмерение связующих точек. Несколько проходов необходимы в случае, когда телеметрическая информация не содержит всех углов ориентирования, или углы известны с точностью 10-30 градусов. Если же телеметрическая информация содержит угловые элементы ориентирования с точностью в несколько единиц градуса, то достаточно и одного прохода – надежность автоматических измерений в этом случае повышается. Для борьбы с возможными грубыми ошибками при автоматических измерениях в PHOTOMOD5.2 введено понятие т.н. «доверительной группы связующих точек», когда программа ищет наибольшее число связующих точек для стереопар с наименьшим поперечным параллаксом, остальные связующие точки, не попавшие в группу, считаются ошибочными.

После измерения связующих и опорных точек производится процедура уравнивания. В ЦФС PHOTOMODможно использовать начальное приближение для алгоритма уравнивания как по уточненной схеме блока, так и построенное другими методами. Начиная с версии 5.2 для уравнивания аэросъемки с БПЛА мы рекомендуем использовать новый режим – уравнивание 3D. При уравнивании в PHOTOMODи достаточном числе опорных точек можно использовать самокалибровку. Это дает возможность использования некалиброванных камер. Ожидаемая точность выходных результатов при строгой фотограмметрической обработке составляет приблизительно 1-2 GSDв плане и 2-4 GSDпо высоте. После фотограмметрического уравнивания, результаты которого и определяют точность выходных продуктов, производится построение рельефа (ЦМР) в автоматическом режиме. При необходимости, после уравнивания может быть сделана стереовекторизация – отрисовка в ручном режиме зданий, сооружений, мостов, дамб и других объектов. Построенный рельеф используется для ортотрансформирования снимков. На последнем этапе из ортотрансформированных снимков создается бесшовная мозаика – производится расчет линий порезов, выравнивание яркостей, стыковка контурных объектов. Самокалибровку можно включать и при отсутствии опорных точек, правда, в этом случае можно рассчитать только коэффициенты k1, k2 радиальной дисторсии. При использовании камер с щелевым затвором можно дополнительно включить расчет аффинных искажений. В случае стабильности углов ориентирования при съемке такая самокалибровка может повысить точность уравнивания.

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

GPSцентров проекций и дисторсией объектива (до нескольких десятков пикселей). В таких случаях можно применять упрощенную автоматизированную последовательность обработки. Бесшовный накидной монтаж указанной точности при этом получается за счет трансформирования исходных снимков в модуле PHOTOMODGeoMosaic.В этом случае используются простейшие методы трансформирования, не учитывающие рельеф местности, а стыковка контуров осуществляется за счет автоматически рассчитываемых связующих точек вдоль автоматически построенных линий порезов.

Примеры фотограмметрической обработки данных аэросъемки с БПЛА

Рассмотрим несколько примеров обработки аэросъемки с БПЛА. Во всех примерах для обработки использовалась ЦФС PHOTOMOD. Отметим, что различными организациями в компанию «Ракурс» для тестирования было передано более 20 блоков аэросъемки с БПЛА. К сожалению, для многих блоков отсутствовали опорные точки и/или съемка была проведена неоткалиброванными камерами. В таких случаях было невозможно оценить точность конечных результатов обработки.

Первый блок, который мы рассмотрим, был снят с борта БПЛА ZALA421-04ф. Данные для исследований были любезно предоставлены ОАО «Газпром космические системы». Блок состоял из 26 маршрутов. Общее число снимков в блоке составило 595. Использовалась предварительно откалиброванная цифровая камера Canon EOS500D. Высота залета над местностью составила около 500 м, размер пиксела на местности приблизительно равен 8 см. На местности были измерены и промаркированы 25 опорных точек, точность координат опорных точек не превышала 10 см. Общий перепад высот местности протяженностью около 3-х километров достаточно большой ~ 70 метров.

Сначала этот же блок аэросъемки был обработан в автоматическом режиме по упрощенной схеме, без уравнивания и использования опорных точек. Привязка осуществлялась по центрам проекции, трансформирование снимков проводилось сразу в модуле GeoMosaicбез учета рельефа. Последующий контроль полученных «псевдо» ортофотопланов по опорным точкам показал расхождения на опорных точках, превышающие 17 м. Такая невысокая точность ортофотплана обусловлена как большим перепадом высот, так и неточностью измерений центров проекций в полете.

Затем блок был подвергнут строгой фотограмметрической обработке. При уравнивании три из измеренных опорных точек считались контрольными. Среднеквадратическая ошибка уравнивания составила по опорным точкам 15 см, 16 см, 12 см, по контрольным точкам 23 см, 29 см и 57 см. Расхождения на связующих точках составили 8 см, 14 см и 69 см. Общий вид блока представлен на следующем рисунке.

Рис. 5. Общий вид «блока 1»

В процессе уравнивания было обнаружено, что координаты центров проекций из телеметрической информации содержат систематическую ошибку, главная из компонент которой составляет 10,5 метра по высоте Z. Среднеквадратические ошибки на центрах проекции после вычитания систематической ошибки составили 84 см, 239 см и 75 см. Существенно большая ошибка по Y(вдоль полета), скорее всего, связана с неточным определением моментов съемки в телеметрии. Большие ошибки по Zна связующих точках возможно связаны с неточной калибровкой камеры и с накопленной ошибкой при съемке камерой с щелевым затвором. Наибольшие ошибки на связующих точках наблюдаются на краях и в углах снимков.

Рис. 6. Величины ошибок на связующих точках

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

Рис. 7. Фрагмент ортофотоплана без выравнивания яркости

В апреле 2011 кафедрой фотограмметрии Московского государственного университета геодезии и картографии (МИИГАиК) были проведены исследования материалов аэрофотосъёмки, полученных с помощью БПЛА Птеро, с целью оценки качества аэросъёмочных работ и фотограмметрической обработки . Съемка выполнялась с высоты около 900 м над средней плоскостью снимаемой местности с борта БПЛА Птеро цифровой фотокомерой CanonEOS5D. Камера была предварительно откалибрована. Для оценки качества материалов использовался фрагмент блока, состоящий из 2-х маршрутов по 6 снимков в каждом. В качестве опорных использовались 14 точек, плановые координаты XYкоторых были сняты с планов масштаба 1:1000, а высота Zопределялась по материалам воздушного лазерного сканирования, выполненного с точностью около 20-30 см. После фотограмметрического уравнивания среднеквадратические погрешности координат на опорных точках составили по X, Yи Zсоответственно 20 см, 21 см и 50 см. Среднеквадратические погрешности координат связующих точек составили 6 см, 6 см, 15 см. Размер пиксела на местности для этого блока GSDсоставляет около 12 см. Общая схема блока показана на следующем рисунке.

Рис. 8. Схема «блока 2» с опорными и связующими точками

Вопросы метрологического обеспечения

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

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

Отметим здесь интересный проект МИИГАиК – с целью отработки и исследования технологий мониторинга и картографирования местности по материалам беспилотной аэрофотосъемки, начаты работы по созданию специализированного исследовательского полигона. Этот полигон, площадью около 50 кв. км, создается в Заокском районе Тульской области, на базе учебного геополигона МИИГАиК, расположенного в 110 км от Москвы.

Территория полигона представляет собой уникальное многообразие картографических объектов. На этой территории расположены разнообразные населенные пункты: поселок городского типа, деревни, дачные и коттеджные поселки; дорожная сеть в виде железных, шоссейных, проселочных и полевых дорог; линии электропередачи различного напряжения; трубопроводы. На территории полигона имеются лесные массивы, различные гидрографические объекты, многообразные формы рельефа, сельскохозяйственные угодья и производственные объекты.

С целью обеспечения отработки и исследования технологий, основанных на применении БПЛА, на территории полигона начаты работы по созданию высокоточной сети планово-высотных опознаков (в виде естественных контуров местности и маркировочных знаков); ведется топографическая наземная съемка характерных участков местности в масштабе 1: 500 и 1: 2000. На эту же территорию по материалам аэрофотосъемки и космическим снимкам высокого разрешения созданы ортофотопланы и цифровые модели местности. По мере поступления новых съемочных материалов эти работы предполагается выполнять в дежурном режиме.

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

Первые испытания планируется провести в середине июля 2011 г. Планируется провести тестовую аэрофотосъемку территории полигона в различных масштабах с помощью отечественного БПЛА «ПТЕРО» с целью отработки и исследования фотограмметрической технологии создания карт различного масштаба по полученным материалам аэрофотосъемки. Фотограмметрическую обработку полученных снимков предполагается выполнить на цифровой фотограмметрической системе PHOTOMOD. В сентябре предполагается провести испытания БПЛА «Х100» бельгийской фирмы Gatewing и БПЛА «МИИГАиК Х8», разработанного в МИИГАиК.

Созданием полигона и проведением на нем испытаний БПЛА и технологий, основанных на их использовании, МИИГАиК намерен помочь потенциальным пользователям освоить и внедрить новые технологии, а разработчикам летательных аппаратов и съемочных систем адаптировать их к решению актуальных задач производства.

Использование БПЛА в качестве аэросъемочной платформы имеет большие перспективы при съемке небольших по протяженности площадных объектов и при съемке линейных объектов. Данные с БПЛА позволяют получать качественные картографические материалы (пространственные данные) при следующих условиях:

· выполнении определенных (вполне посильных) требований к съемочной аппаратуре и процессу съемки (гарантия достаточности перекрытий);

· строгой фотограмметрической обработке. Точность при этом возрастает в десятки раз и может составлять около GSD, как и для обычной аэросъемки и космических снимков.

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

· Использовать на БПЛА калиброванные камеры.

· Производить съемку с выдержкой не длиннее 1/250с.

· Использовать объективы с фиксированным фокусным расстоянием. Если это невозможно, следует фиксировать увеличение (Zoom). Съемка должна производиться с фокусировкой на бесконечность и с отключенным режимом автофокусировки.

· Проектировать съемку с увеличенными перекрытиями (80% вдоль, 40% поперек маршрута).

· Желательно использовать камеры с центральным затвором.

· Желательно использовать двухдиапазонные GPS приемники на борту и дифференциальный режим измерений.

· Желательно использование на борту IMU, пусть и не имеющего высокой точности.

Благодарности

Благодарим компании: «Беспилотные системы ЗАЛА АЭРО», ОАО «Газпром космические системы», «АФМ-Серверс», ООО «Геометр-Центр», НПИ и КЦ «Земинформ», ЗАО «Транзас», ЗАО «Лимб» за помощь в подготовке материала, предоставление данных и полезные обсуждения.

Литература

1. Чибуничев А.Г., Михайлов А.П., Говоров А.В. Калибровка цифровых фотокамер: Вторая научно-практическая конференция РОФДЗ. Тезисы докладов. М., 2001 г. с38-39.

2. Скубиев С.И., Научно-производственный институт земельно-информационных технологий Государственного университета по землеустройству «Земинформ» (Россия), Использование беспилотных летательных аппаратов для целей картографии. Тезисы XЮбилейной международной научно-технической конференции «От снимка к карте: цифровые фотограмметрические технологии». Гаета, Италия, 2010.

3. Результаты полевых исследований БПЛА «Птеро»

В одном мы уверены точно: высокая цена не всегда означает высокое качество.

Мы окунемся в индустрию и узнаем, как дроны покажут себя при съемке.

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


Описание

Топографическая съемка является неотъемлемой частью всех проектов в области землеустройства.

В этом примере мы рассмотрим участок земли, на котором должен был быть построен новый поселок. До начала работ было необходимо провести точную топографическую съемку по нескольким причинам:


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

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

Геодезия 101

Для традиционной топографической съемки требуется сбор координат точек в заранее определенной сетке. В этом случае использовалась сетка размером 150х150 см:

Измерения производились каждые 150 сантиметров, на каждом перекрестке:


Всего на площади съемки 34,5 Га было собрано 1632 координаты.

Без дрона, снимающего со скоростью 20 точек/час (1 точка, каждые 3 минуты), сбор данных занял бы приблизительно 82 часа.


82 часа традиционной съемки означают, что инженер вынужден ждать как минимум неделю, чтобы приступить к обработке данных. Далее понадобится еще 3-4 дня, прежде чем работа будет сделана.

Проведя ту же съемку с использованием БПЛА, полевая команда смогла предоставить разработчику более быстрый вариант обзора.

Прежде всего, не нужно было собирать 1600 точек по всей площади. Вместо этого потребовалась съемка всего 10 наземных меток, расположенных в зоне обзора:

Для более крупных проектов Наземные опорные точки (GCP) лучше расставить по сетке.

10 наземных меток или 1632 точки:

10 опорных меток могут быть сделаны за 1-2 часа.

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

Завершив сбор GCP, были собраны точки традиционным методом в участках со стоячей водой - комбинация двух методов, описанных выше.

Конечные собранные точки:


В итоге мы получили 117 точек (10 GCP + 107 на участках со стоячей водой).

Время на съемку:

Теоретически: 10 наземных меток + сбор точек = 1-2 часа

Фактически: 117 точек (10 GCP + 107 на участках со стоячей водой) при скорости сбора 20 точек / час = 5,85 часа

Традиционный метод: 1,632 точки при скорости сбора 20 точек / час = 81,6 часа



В течение часа были завершены все действия с БПЛА, включая сборку, предполетные проверки, запуск, посадку, разборку и первоначальную сшивку карты.


Таким образом мы получили:

БПЛА (1 час) + сбор точек (5,8 часа) =

Общее время полевых работ: 6,8 часа

Сравнение:

34,5 Га/ полевые работы с использованием БПЛА = 6,8 часа

34,5 Га/ полевые работы по традиционному методу = 81,6 часа

Общая экономия: 74,8 часа

Анализ данных

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

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


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

После 2 дней дополнительной обработки, результирующие топографические контуры стали точными в пределах 9 сантиметров, как по горизонтали (X, Y), так и по вертикали (Z):





Общие сроки выполнения проекта:

Метод с использованием БПЛА::

Полевые работы (6.8 часов) + обработка данных (24 часа) =

30,8 часов (около 4 дней)

Обычный метод:

Полевые работы (81,6 часов) + Обработка данных (24 часа)=

105,6 часов (около 13 дней)



Используя технологию с использованием беспилотника, инженер получил окончательный топографический обзор примерно за 75 часов

По полученным данным выяснилось, что:

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

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

3. Инженер узнал о недорогой и рентабельной съемке с БПЛА и планирует снова использовать этот метод для проведения окончательного «встроенного» топографического исследования в ближайшие недели.

Здесь Вы можете больше и лучшие модели беспилотников.

Данные съёмки с БПЛА, показанные на этой странице, предоставлены . Технология обработки материалов съемки в ПО Agisoft PhotoScan предоставлена ООО "Плаз" .

Применение беспилотных летательных аппаратов (БПЛА) позволяет существенно снизить затраты на производство аэрофотосъемочных работ. С точки зрения традиционной фотограмметрии качество подобной съемки вероятнее всего будет оценено, как неприемлемое, поскольку на БПЛА, как правило, устанавливаются камеры бытового сегмента, не используется гиростабилизирующая аппаратура, при съемке нередки отклонения оптических осей от вертикали в несколько градусов, что значительно усложняет процесс первичной обработки снимков. Однако для современного фотограмметрического программного обеспечения эти недостатки не представляют значительных проблем. Более того, развитие цифровых методов фотограмметрической обработки уже привело к появлению программ и программных комплексов, способных обрабатывать даже такие "некачественные" данные аэрофотосъемки в высокоавтоматизированном режиме, при минимальном участии оператора.

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

  • БПЛА для выполнения аэрофотосъемки;
  • ПО Agisoft PhotoScan в качестве инструмента обработки материалов съемки;
  • инструментарий ГИС Панорама для векторизации ортофотопланов и получения топокарт.

Аэрофотосъемка с использованием БПЛА

В техническом плане процесс аэрофотосъемки с использованием БПЛА состоит из трех этапов: подготовительного, собственно съемки, и постобработки полученных данных.

Подготовительный этап
На данном этапе производится:

  • изучение имеющихся материалов; формирование или сбор требований к материалам, которые нужно получить по результатам съемки – тип и масштаб карты, границы объекта съемки; приведение их в технические требования к съемочным материалам: разрешение, координаты контура участка съемки, перекрытие снимков, точность определения координат центров фотографирования, требования к наземной опорной сети (при комбинированной съемке, например, когда привязка фотоплана производится по точкам наземной опорной сети, требования к точности определения КЦФ вообще не предъявляются);
  • формирование полетного задания для БПЛА. Выполняется программой – планировщиком полета, входящей в состав комплекса. Оператор должен выбрать используемый комплекс БПЛА (в случае, если программа позволяет работать с несколькими конфигурациями БПЛА и фотоаппаратуры), задать на карте контур участка съемки и примерное положение стартовой площадки, установить требуемое разрешение и перекрытие, после чего программа рассчитывает план полета и проверяет его выполнимость.

Выполнение аэрофотосъемки
По прибытии на стартовую площадку производится:

  • уточнение положения стартовой площадки, задание точки возвращения и ввод данных о скорости и направлении ветра на рабочей высоте, если таковые известны;
  • автоматическое уточнение плана полета и повторная проверка его выполнимости;
  • старт БПЛА с пускового устройства;
  • выполнение съемки в автоматическом режиме;
  • посадка.

Выполнение съемки местности с использованием БПЛА

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

Постобработка данных
Заключается в:

  • снятии данных (фотоснимки и журнал полета) с бортовых носителей информации;
  • визуальной оценке качества фотографий и отбраковке "технических" кадров, если такие записаны. Под техническими кадрами понимаются снимки, сделанные вне пределов участка съемки - при подлете к участку, на дугах разворота и т.п.;
  • генерация файла привязки центров фотографирования. В ходе полета аппаратура управления ведет запись различных параметров, среди которых – координаты, скорость и параметры ориентирования летательного аппарата. После окончания съемки из файла журнала полета необходимо выбрать координаты, соответствующие моментам фотографирования, и приписать их конкретным снимкам. Такая обработка, как правило, выполняется в той же программе – планировщике полетного задания.

В соответствии с требованиями отраслевых инструкций , для получения топокарт масштаба 1:2000 необходима фотооснова, имеющая разрешение 15 см/пикс и имеющая погрешность определения координат в каждой точке не выше 60 см. Такое разрешение легко обеспечивается при съемке с БПЛА с использованием компактных фотоаппаратов. Например, съемка камерами типа Canon S-95 или Sony NEX-5 (с объективом SEL30M35) с высоты порядка 200-300 м дает снимки, имеющие разрешение 5 см/пикс.

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

Обработка аэрофотоснимков в ПО Agisoft PhotoScan

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

Построение и привязка модели местности в программе состоит из трех основных этапов:

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

Затем полученная модель используется для генерации ортофотопланов и матриц высот.

С точки зрения оператора процесс работы с программой выглядит следующим образом:

  • Выбор системы координат и загрузка данных привязки центров фотографирования
  • Формирование точечной модели поверхности Земли
  • При наличии наземной опорной сети – установка отметок опорных точек на фотоснимках и загрузка координат точек опорной сети
  • Оптимизация модели (уравнивание параметров привязки)
  • Генерация полигональной модели поверхности Земли
  • Экспорт данных – ортофотоплан, матрица высот
  • Приведенные скриншоты окна программы наглядно иллюстрируют процесс обработки материалов аэрофотосъемки на примере съемки полигона "Заокский", материалы которой предоставлены ОАО "Газпром космические системы" . Обработка данных материалов на ПК, оснащенном 4-хядерным процессором Intel Core i7 2600K и имеющем 16 Гб оперативной памяти, заняла порядка трех-четырех часов – от загрузки фотографий до экспорта ортофотоплана и цифровой модели местности в формате GeoTiff. Из этого времени около одного часа ушло на дешифирование и маркирование опорных точек – ручной труд оператора, а остальное время заняло выполнение расчетов.

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

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

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

    1) Фотографии загружены. В свойствах проекта видно, что проект состоит из блоков (chunks) – обрабатываемых независимо частей проекта со своими фотографиями, моделью, СК, параметрами калибровки оптики и т.п. В данном проекте - один блок, состоящий из 415 фотографий. Метки NA (not aligned) рядом с фотографиями показывают, что положение этих снимков в пространстве модели еще не известно.

    2) Выбор системы координат

    4) Метки в форме синих шариков отображают взаимное расположение точек съемки (КЦФ), после уравнивания они будут заменены метками другого вида, соответствующим положению плоскостей кадров

    5) После выполнения первого этапа обработки – первичного уравнивания и построения точечной модели, формируется облако точек, описывающее модель, и набор параметров взаимного ориентирования фотоснимков. Положение выбранного снимка отображается в области просмотра модели. Снимки, которые не удалось уравнять, по-прежнему отображаются сферами/шариками, и в списке фотографий отмечены меткой NA (not aligned). В данном проекте таких нет

    6) Установка маркеров (меток опорных точек). Если известно положение маркеров на снимках (в системе координат снимка), можно просто импортировать эти данные в PhotoScan. Если маркеры еще не дешифрированы, придется задавать их расположение прямо в программе. Для каждого маркера достаточно отметить их положение на одном-двух снимках, и PhotoScan автоматически определяет их положение на других снимках, выделяя снимки, на которых присутствует выбранный маркер, специальными метками. На каждом снимке можно подтвердить итли уточнить автоматически выбранное положение маркера

    7) Маркеры расставлены. Можно выполнять построение модели местности

    8) Модель готова. Ее можно экспортировать как матрицу высот (цифровую модель местности) и сформировать на основе этой модели ортофотоплан местности.

    9) В завершение можно построить текстуру модели и рассматривать ее прямо в программе.

    10) Внутреннее представление модели поверхности Землки в PhotoScan - сеть триангуляции Делоне, TIN модель

    11) Ортофотоплан всего участка съемки.

    12) Цифровая модель местности всего участка съемки

    Получение карт на основе ортофотопланов в ГИС Панорама

    Комплекс автоматизированного дешифрирования и векторизации по данным ДЗЗ, разработанный на базе ГИС "Панорама" , предназначен для автоматической векторизации линейных и площадных объектов по цветным растровым изображениям земной поверхности.

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

    • предварительная обработка растра;
    • классификация;
    • обработка растра классификации;
    • преобразование растра в вектор;
    • векторная обработка.

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

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

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

    К загруженному ортофотоплану в формате GeoTiff применяется технология классификации и распознования

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

    На следующем этапе производиться фильтрация лишней информации, ее сглаживание и перевод с в линейный и площадной вид

    Для этого используются морфологические операции - изменение бинарного состояния пикселя на основе анализа состояния его соседей. К таким операциям относятся:

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

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

    Результат автоматизированного дешифрирования и векторизации ортофотопланов можно посмотреть и отредактировать в ГИС "Панорама"

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

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

    Литература
    1. Инструкция по фотограмметрическим работам при создании цифровых топографических карт и планов ГКИНП (ГНТА)-02-036-02. Москва, ЦНИИГАиК, 2002