Modelling optical propagation for image reconstruction in marine mineral exploration
Modelling optical propagation for image reconstruction in marine mineral exploration
Abstract
This work presents the results of numerical modelling of optical radiation propagation in water taking into account the main physical parameters to evaluate the transmission coefficient of the medium in order to reconstruct underwater images of the seafloor obtained with the use of deep-sea robotic systems within the framework of marine mineral exploration.
It is demonstrated that due to the application of the model based on the numerical solution of the radiation transport equation using the Monte Carlo method and taking into account the main physical parameters of the water medium, the results obtained provide high accuracy in determining the characteristics of optical radiation propagation.
The obtained simulation results for different parameters of the aquatic environment, different distances between the imaging object (source) and the camera lens (receiver) will allow to form a matrix of underwater object colour correction without the necessity of using large computing power and in real time.
1. Введение
Проведение морских геологоразведочных работ по исследованию месторождений железомарганцевых конкреций (ЖМК), кобальтмарганцевых корок (КМК) и глубоководных полиметаллических сульфидов (ГПС) в районах Мирового океана с глубиной до 6000 м осуществляется с помощью специализированного оборудования, осуществляющего гидроакустическую, фото- и видеосъемку морского дна.
Результаты фотосъемки морского дна в дальнейшем обрабатываются и анализируются с целью формирования сведений о рудоносности исследуемого месторождения. Качество съемки оказывает определяющее влияние на эффективность анализа и достоверность выводов. Качество подводной фото- и видеосъемки, в свою очередь, определяется многими факторами, включая особенности освещения, расстояние до объекта съемки, мутность воды и пр., которые приводят к ухудшению контрастности, четкости, яркости фотографий, искажению цвета, появлению дымки , .
В связи с вышеизложенным разработка методики восстановления подводных изображений , с применением современных цифровых интеллектуальных технологий, обеспечивающей минимизацию искажений, вносимых водой и взвешенными в воде частицами, повышение контрастности изображения и восстановление реальных цветов объектов, и ее внедрение в практику геологоразведочных работ повысит их эффективность и достоверность, позволит получить более точные данные для дальнейшей разработки месторождений твердых полезных ископаемых в различных районах Мирового океана.
2. Распространение оптического излучения в воде
Распространение оптического излучения в воде определяется в первую очередь такими физическими процессами, как поглощение и рассеяние оптического излучения.
Поглощение – это физический процесс, при котором фотоны теряют свою энергию, которая преобразуется в другие формы, главным образом в тепловую и химическую (фотосинтез).
Рассеяние вызвано вариациями показателя преломления, которые приводят к изменению направления распространения фотонов.
Явления поглощения и рассеяния оптического излучения приводят к появлению трех нежелательных эффектов:
- потери энергии оптического сигнала из-за наличия поглощения, что приводит к уменьшению дальности связи;
- расширение светового луча вследствие рассеяния в совокупности с ограниченной апертурой линзы приемника приводит к уменьшению количества фотонов, достигающих детектора оптического излучения;
- попадание каждого отдельного фотона на чувствительный элемент приемника из-за наличия рассеяния происходит в разные моменты времени (дисперсия), что может привести к межсимвольным искажениям.
Рисунок 1 - Примеры изображений морского дна, полученные в ходе морских геологоразведочных работ
3. Модель водной среды
Наиболее простым способом методом моделирования распространения оптического излучения под водой является использование закона Бера-Ламберта . Однако в этом случае невозможно учесть явления многолучевого распространения, многократного рассеяния, дисперсии оптического сигнала, которые являются очень важными для точной оценки распространения оптического излучения под водой.
Более точным подходом к моделированию подводного оптического канала распространения оптического излучения является решение уравнения переноса излучения (RTE – Radiative transfer equation) . В связи с тем, что уравнение является интегро-дифференциальным с несколькими независимыми переменными , то аналитическое решение в общем виде не может быть найдено. В связи с вышеизложенным наиболее широкое распространение получили численные методы решения уравнения переноса излучения .
В настоящей работе модель водной среды осуществлено численное решение уравнения переноса излучения с использованием метода Монте Карло , который является вероятностным методом, основанным на трассировке распространения большого количества оптических лучей.
Для точного определения параметров модели на этапе инициализации расчетов осуществляется задание длины волны оптического сигнала, глубины расположения канала связи, данные о концентрации хлорофилла в воде, положение и ориентация источника и приемника оптического излучения. На основе введенного значения глубины расположения канала связи определяется значение коэффициента преломления воды, температура и соленость воды.
Для учета рассеяния используется модель, основанная на функции Хеньи-Гринштейна (Henyey-Greenstein), которая имеет параметр «g» – фактор анизотропии, равный 0,9185.
В связи с тем, что в цветовом пространстве RGB цвета изображения могут быть разложены на три основных составляющих (красная, синяя и зеленая), то моделирование необходимо проводить для трех длин волн, соответствующих красной (690 нм), синей (530 нм) и зеленой (485 нм) областям спектра оптического излучения. Изменение интенсивности каждой отдельной составляющей изображения позволяет обеспечить восстановление естественных цветов подводных объектов за счет компенсации повышенного ослабления красного света при распространении под водой.
Рисунок 2 - Результаты моделирования распространения оптического излучения с длиной волны 690 нм
Примечание: красный свет
Рисунок 3 - Результаты моделирования распространения оптического излучения с длиной волны 530 нм
Примечание: зеленый свет
Рисунок 4 - Результаты моделирования распространения оптического излучения с длиной волны 485 нм
Примечание: синий свет
Рисунок 5 - Гистограмма распределения оптических лучей по времени прихода в точку приема для длины волны 690 нм
Примечание: красный свет
Рисунок 6 - Гистограмма распределения оптических лучей по времени прихода в точку приема для длины волны 530 нм
Примечание: зеленый свет
Рисунок 7 - Гистограмма распределения оптических лучей по времени прихода в точку приема для длины волны 485 нм
Примечание: синий свет
4. Заключение
В настоящей работе представлены результаты численного моделирования распространения оптического излучения через водную среду для оценки коэффициента передачи среды с целью восстановления цветности подводных изображений морского дна, получаемых с применением глубоководных роботизированных комплексов в рамках морских геологоразведочных работ.
В связи с применением модели, основанной на численном решении уравнения переноса излучения с помощью метода Монте Карло и учитывающей основные физические параметры водной среды, получаемые результаты обеспечивают повышенную точность определения характеристик распространения оптического излучения.
Полученные результаты моделирования для различных параметров водной среды, различных расстояний между объектом съемки (источник) и объективом камеры (приемник) позволят сформировать матрицу корректировки цветов подводного объекта без необходимости использования больших вычислительных мощностей и в режиме реального времени.