Книги по разным темам Pages:     | 1 | 2 | Электронный журнал ИССЛЕДОВАНО В РОССИИ 1556 Анализ качества стереореконструкции модельной облачной сцены Козулин Ю. А. (maksar79@mailru.com), Белкин М. С.

ООО Астроинформ СПЕ Введение Одним из видов обработки, интерпретации и анализа данных дистанционного зондирования Земли является стереореконструкция или восстановление трехмерного рельефа исследуемой сцены по совокупности двух или более стереоснимков.

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

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

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

Задача трехмерной реконструкции облачной сцены имеет ряд особенностей, требующих специфических методов решения. Это связано в первую очередь с диффузным, слабоконтрастным характером текстуры облачности, а так же с разрывностью (фрактальностью) рельефа облачной сцены. Кроме того, такого типа сцена достаточно динамична и, строго говоря, требует синхронной стереосъемки. Таким образом, задачи синхронного стереомониторинга облачных структур нуждаются в разработке новых средств 3D реконструкции и формулировке требований к системам стереосъемки на базе оценки их потенциальных возможностей. Данная статья посвящена анализу качества стереореконструкции облачной модельной сцены при различных параметрах геометрии наблюдения.

Электронный журнал ИССЛЕДОВАНО В РОССИИ 1557 Постановка задачи и основные этапы решения Работа имела целью исследование зависимости качества восстановления облачной модельной сцены по стереопарам, полученным при модельной съемке серии изображений для разных положений наблюдателя, определяемых наклонной дальностью и ракурсами съемки.

Решение поставленной задачи состояло из следующих этапов:

Разработан алгоритм стереореконструкции облачных полей.

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

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

Алгоритм обработки реализован в виде программных модулей, написанных на языке IDL (Interactive Data Language). IDL - это среда программирования, ориентированная на обработку массивов, и содержащая богатую библиотеку процедур математического анализа и визуализации данных. Пакет специально предназначен для инженерных и научных расчетов.

Произведен анализ полученных результатов.

Исходные данные В качестве рабочего материала используются данные моделирования облачной синтетической сцены (Synthetic Clouds Scene), предоставленные американскими коллегами.

Эта облачная сцена создана с помощью программного пакета CloudScape [1]. Исходные данные состоят из набора изображений модельной сцены облачного поля и карты высот верхней границы облачности модельного поля. Изображения сцены получены при наблюдении камерами спутников Lead и Trail под предварительно заданными зенитными углами. Наличие карты модельной сцены и изображений модельной облачности, полученных при известных параметрах геометрии съемки, предоставляет уникальную возможность для определения точности и полноты восстановления сцены в зависимости от зенитных углов наблюдения.

Для реконструкции сцены формируются стереопары на основе изображений модельной облачности. Геометрия съемки облачной сцены представлена на рис. 1. Орбита спутников Lead и Trail круговая. Земля в модели облачной сцены предполагается имеющей форму сферы. Высота орбиты спутников - 500.0 км. Модельная облачность представляет собой кучево-дождевые облака - тип Cumulonimbus (Cb) по международной классификации облачности. Облачное поле расположено на высотах от 4.5 до 9.5 км.

Точка пересечения оптических осей двух объективов находится на высоте 7.0 км над поверхностью Земли. Размер исходных изображений стереопар и карты высот верхней границы облачности: 1024 на 1024 пикселя. Поле зрения камер спутников Lead и Trail составляет 2.0 на 2.0 градуса. Пространственное разрешение карты высот - 25 м.

Электронный журнал ИССЛЕДОВАНО В РОССИИ 1558 Под зенитным углом будем понимать угол между перпендикуляром к земной поверхности OgrOorb, проходящим через точку пересечения оптических осей камер, и одной из оптических осей. Восстановление облачности проводится в двух вариантах: для симметричного случая наблюдения (зенитные углы L для Lead и T для Trail лежат по разную сторону и равны) и для несимметричного случая. Изображения сцены получены для следующих зенитных углов наблюдения L и T: 0, 1, 2, 3, 5, 7, 10 и 15 градусов.

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

Введем понятие лугловой базы: под угловой базой будем понимать угол B, равный разности (L - T) с учетом знаков углов. Для несимметричного случая съемки введем угол асимметрии A равный (T + L)/2.

Расчетные варианты и основные этапы обработки Введем единую прямоугольную систему координат с началом, расположенным в точке Ocl пересечения оптических осей объективов камер. Ось Z направлена вверх от поверхности Земли. Ось X направлена в направлении движения спутников (от Trail к Lead) параллельно касательной к точке Oorb (см. рис. 1). Ось Y дополняет систему координат до правой. Данная система координат используется при вычислении поля высот облачности относительно горизонтальной плоскости, проходящей через точку Ocl. Окончательно поле высот пересчитывается относительно горизонтальной плоскости, проходящей через точку Ogr, т.е. относительно поверхности Земли.

Расчетный вариант №1 (Симметричный случай): Производится обработка стереопар для зенитных углов |T|=|L|: 1, 2, 3, 5, 7, 10 и 15 градусов.

Расчетный вариант №2 (Несимметричный случай): Производится обработка стереопар с угловой базой B, равной 2, 3, 4 и 5 градусов.

Для восстановления поля высот облачности исходные изображения каждой стереопары были предварительно пересчитаны в единую прямоугольную систему Электронный журнал ИССЛЕДОВАНО В РОССИИ 1559 координат с приведением изображений к эпиполярной стереопаре [2]. Рассмотрим сущность приведения к эпиполярной стереопаре.

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

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

Рис. 2. Приведение изображений к эпиполярной стереопаре В геометрии съемки процесс взаимного положения и ориентации камер относят к задаче относительной ориентации, которая позволяет получить преобразование, связывающее два различных положения фиксированной относительно камеры системы координат сцены в различные моменты съемки. Иногда же возникает необходимость связать пространственные координаты, полученные при помощи стереопары, с координатами, измеренными в другой системе. Эту проблему относят к задаче абсолютной ориентации. Рис. 2 иллюстрирует некоторые особенности относительной ориентации в задачах стереореконструкции и понятие об эпиполярных линиях. Если нам известна относительная ориентация камер, то мы можем вычислить положение точки по ее известным координатам на левом и правом изображениях. На рис. 2 точка P является элементом исследуемой сцены, PL и PR - ее отображениями на левом и правом изображениях, SL и SR - центы линз левой и правой камер. (Прямая линия, соединяющая центры объективов, называется стереобазой). Положение правой камеры отображается в точку SL* на левом изображении, а так называемая точка схода для луча от правой камеры к элементу поверхности лежит в точке CL*. Точки SL*, PL и CL* лежат на одной прямой, которая образуется пересечением плоскости изображения с плоскостью, содержащей точки P, SL, SR. Аналогичная прямая получается и для правого изображения. Обе эти линии называют эпиполярными линиями. Объект, изображенный на эпиполярной линии на левом снимке, на правом должен изображаться также на эпиполярной линии, если он вообще присутствует на правом снимке. В процессе стереореконструкции для эффективного восстановления применяют преобразования, приводящие стереопару к Электронный журнал ИССЛЕДОВАНО В РОССИИ 1560 эпиполярной. Данное преобразование модифицирует изображения стереопар таким образом, что для некоторой системы координат эпиполярные линии становятся параллельны одной из ее осей и диспаратность проявляется только вдоль этой оси.

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

Таким образом, в результате приведения обрабатываемых изображений к эпиполярной стереопаре эпиполярные линии трансформированных изображений параллельны оси X (см. рис. 1) и диспаратность проявляется только вдоль этой оси.

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

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

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

Дискретность привязки площадок составляет 1 пиксель. В таблице 1 приведены параметры корреляционного алгоритма для симметричного случая.

Таблица Зенитный угол (|T|=|L|), 1 2 3 5 7 10 Размер области поиска по X, 11 21 31 53 73 103 пикс.

Размер площадки, 15x15 15x15 15x15 15x15 15x15 25x25 25xпикс. X пикс.

Для всех полученных пар сопряженных точек был проведен расчет поля диспаратности. В качестве опорного кадра для всего набора стереопар был взят Lead-кадр.

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

На основе данных геометрии съемки и полей диспаратности по X и Y были определены истинные положения сопряженных точек X и Y и высоты Z фрагментов облачного поля, отображенных в виде сопряженных пар точек в стереопарах. Данные высоты получены относительно точки пересечения линий визирования 7.0 км, поэтому полученное поле Z пересчитывается в высоты относительно поверхности Земли.

Электронный журнал ИССЛЕДОВАНО В РОССИИ 1561 Результаты обработки Симметричный случай Средняя высота и среднеквадратичное отклонение высоты облачной сцены по исходной карте высот составляют 7468 м и 1573 м соответственно. В таблице 2 приведены средние высоты и среднеквадратичные отклонения высоты облачности для восстановленной модельной сцены и для точек карты высот, соответствующих найденным сопряженным точкам.

Таблица Геометрия съемки Зенитный угол (|T| = |L|), 1 2 3 5 7 10 Угловая база B, 2 4 6 10 14 20 Восстановленное поле высот Средняя высота, м 7808 7854 7874 7913 7906 7977 Среднекв. отклонение, м 616 582 555 526 538 466 Карта высот по массиву точек восстановленных фрагментов облачности Средняя высота, м 7802 7809 7829 7875 7867 7940 Среднекв. отклонение, м 648 644 619 586 604 552 На рисунке 3 приведено восстановленное поле высот модельной облачности (зенитный угол наблюдения (|T|, |L|) 3 градуса) в сравнении с реальной картой высот, приведенной на рисунке 4. Градации серых тонов на изображениях соответствуют высоте.

Для каждого из изображений на рисунке приведена шкала соответствия оттенка высоте.

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

Pages:     | 1 | 2 |    Книги по разным темам