СПОСОБ ВОССТАНОВЛЕНИЯ ИЗОБРАЖЕНИЙ В ДВУХКАНАЛЬНОЙ СКАНИРУЮЩЕЙ СИСТЕМЕ

Изобретение относится к пассивным [1, 2] двухканальным сканирующим системам наблюдения, работающим в оптическом, инфракрасном или миллиметровом диапазонах длин волн.

Система наблюдения состоит из двух измерительных каналов. Каждый канал содержит приемное устройство - приемник, регистрирующий излучаемое объектами поле в ограниченном объеме пространства в соответствии со своей диаграммой направленности (ДН), а также тракт первичной обработки принимаемых сигналов. Приемники одновременно сканируют зону обзора размером M×N элементов дискретизации по азимуту и углу места. Первый приемник движется непрерывно вдоль строки (по азимуту) со съемом данных с малым шагом дискретизации h и осуществляет переход к другой строке с увеличенным шагом k₁⋅h, k₁>1, где (k₁-1) - число пропущенных подряд строк. Второй приемник, наоборот, движется непрерывно вдоль столбца (по углу места) со съемом данных с малым шагом дискретизации h и осуществляет переход к другому столбцу с увеличенным шагом k₂⋅h, k₂>1, где (k₂-1) - число пропущенных подряд столбцов.

По результатам первичной обработки формируются две матрицы наблюдений: Y₁={y₁(i,j)}, , и Y₂={y₂(i,j)}, , . Пропущенные при сканировании строки и столбцы в матрицах Y₁ и Y₂ присутствуют, но не рассматриваются.

Модель наблюдений имеет следующий вид:

, ,

,

, ,

где α(i,j) - весовая функция с областью определения , , представляющая нормированную ДН приемника; X={x(i,j)} - матрица искомого изображения с элементами x(i,j), подлежащими восстановлению на множестве элементов дискретизации , ; P₁(i,j) и p₂(i,j) - шумы аппаратуры в приемных каналах в виде белого шума.

Задача заключается в восстановлении изображения X={x(i,j)} на множестве элементов дискретизации , путем обработки полученных наблюдений Y₁={y₁(i,j)} и Y₂={y₂(i,j)}.

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

Оптимальное решение задачи восстановления изображения X с позиции известных критериев [4] для двумерной модели измерений (1) в пространственной области требует больших предварительных вычислительных затрат и памяти при обращении матриц, что снижает устойчивость решения. При больших значениях М и N реализовать матричные методы не удается.

Оптимальное решение в частотной области [5] требует обработки всех строк и столбцов на множестве , , что достигается интерполяцией пропущенных строк и столбцов матриц Y₁ и Y₂. Несмотря на удобство реализации частотных методов, ошибки интерполяции наблюдений существенно снижают точность восстановления изображений.

На практике целесообразно перейти от (1) к упрощенной модели наблюдений, не учитывающей пропуски строк и столбцов:

,

где α₁(j) и α₂{i) - соответственно центральное горизонтальное и вертикальное сечения ДН, причем α₁(0)=α₂(0); - случайные аддитивные составляющие, порожденные шумами аппаратуры и ошибками аппроксимации при переходе от (1) к (2).

Число оцениваемых параметров в (2) значительно меньше, чем в (1).

В рамках модели (2) ставится одномерная задача независимого восстановления изображений в i-x строках Х₁(i) и j-x столбцах Х₂(j) матрицы X путем обработки i-х строк Y₁(i), матрицы Y₁ и j-x столбцов Y₂(j), матрицы Y₂. Задача решается известными способами.

В качестве прототипа может быть рассмотрен любой способ восстановления одномерного изображения: матричный способ [4] или способ фильтра Винера [5], действие которых применительно к поставленной задаче сводится к следующему:

1. На вход подается i-я строка Y₁(i) первой матрицы наблюдений Y₁ или j-й столбец Y₂(j) второй матрицы наблюдений Y₂ - одномерный массив числовых данных объема М или N.

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

3. В результате действия оператора восстановления на выходе получается i-я строка Х₁(i) или j-й столбец Х₂(j) восстановленного изображения в виде нового массива данных объема М или N.

Совокупность восстановленных строк X₁(i), и столбцов X₂(j), N дает матрицу X_p решетчатой структуры восстановленного изображения, в которой пропущенные строки и столбцы присутствуют, но не рассматриваются.

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

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

Технический результат предлагаемого технического решения достигается тем, что способ восстановления изображений в двухканальной сканирующей системе заключается в том, что при наблюдении зоны обзора с помощью двух приемников, первый из которых дает матрицу наблюдений Y₁ с пропусками строк, а второй - матрицу наблюдений Y₂ с пропусками столбцов, обрабатывают наблюдаемые i-е строки Y₁(i) матрицы Y₁ и наблюдаемые j-е столбцы Y₂(j) матрицы Y₂ оператором восстановления одномерных массивов и в результате этой обработки получают в той же нумерации i-е строки Х₁(i) матрицы Х₁ восстановленного изображения по строкам и j-е столбцы Х₂(j) матрицы Х₂ восстановленного изображения по столбцам, отличающийся тем, что пропущенные строки матрицы Х₁ восстанавливают путем линейной интерполяции соседних не пропущенных строк и получают матрицу , а пропущенные столбцы матрицы Х₂ восстанавливают путем линейной интерполяции соседних не пропущенных столбцов и получают матрицу , затем все соответствующие пары элементов матриц и сравнивают к заданным уровнем γ и выбирают из каждой пары один элемент, наиболее близкий в этому уровню, после чего выбранные элементы помещают в матрицу X, которая представляет восстановленное изображение без пропусков строк и столбцов.

Способ осуществляют следующим образом:

1. Два приемника одновременно сканируют зону обзора размером M×N элементов дискретизации по азимуту (по j) и углу места (по i).

2. По результатам первичной обработки формируются две матрицы наблюдений: Y₁={y₁(i,j)}, , и Y₂={y₂(i,j)}, , , где k₁, k₂>1. Первый приемник дает матрицу наблюдений Y₁ с пропусками строк, а второй - матрицу наблюдений Y₂ с пропусками столбцов. Пропущенные при сканировании строки или столбцы присутствуют в матрицах Y₁ и Y₂, но не рассматриваются.

3. Наблюдаемые i-e строки Y₁(i) матрицы Y₁ и наблюдаемые j-е столбцы Y₂(j) матрицы Y₂ обрабатывают оператором восстановления одномерных массивов и в результате этой обработки получают в той же нумерации i-e строки X₁(i), матрицы Х₁ восстановленного изображения по строкам и j-е столбцы Х₂(j), матрицы Х₂ восстановленного изображения по столбцам.

4. Пропущенные строки матрицы Х₁ заполняются методом линейной интерполяции: между соседними наблюдаемыми строками матрицы Х₁ помещаются k₁-1 строк по числу k₁-1 недостающих элементов дискретизации угла места, а элементы этих строк находятся линейной интерполяцией соответствующих элементов i-й и (i+1)-й строк матрицы Х₁ по формуле:

Δх₁=(х₁(i+k₁,j)-х₁(i,j))/k₁,

x₁(i+i₁,j)=y(i,j)+Δx₁⋅i₁, , , .

Подобным образом заполняются пропущенные столбцы матрицы Х₂:

Δх₂=(х₂(i+k₂,j)-х₂(i,j))/k₂,

x₂(i,j+j₂)=x₂(i,j)+Δx₂⋅j₂, , , .

5. После интерполяции получается матрица с повышенным разрешением по азимуту и матрица с повышенным разрешением по углу места.

Затем матрицы и объединяются в M×N-матрице оценок искомого изображения X следующим образом.

6. Все пары , соответствующих элементов матриц и , сравниваются с заданным уровнем γ. Из двух элементов и выбирается один элемент, наиболее близкий к этому уровню. Он запоминается как i-й, j-й элемент x(i,j) матрицы X, представляющей восстановленное изображение без пропусков строк и столбцов.

7. Уровень γ назначается человеком-оператором из физических соображений. Уровень может меняться в процессе наблюдения с целью повышения четкости отдельных деталей изображения.

Результаты моделирования. Способ реализовывался с применением операторов восстановления одномерных изображений в пространственной области (матричный метод) и частотной области (метод фильтра Винера). Методы показали близкую точность восстановления изображения вдоль строк и столбцов. Дополнительно небольшим порогом снимались шумовые эффекты на восстановленном изображении. При моделировании ширина ДН составляла (2m+1)×(2n+1)=7×7, отношение сигнал-шум (С-Ш) 30 и 50 при максимальной амплитуде 5, шаг по строкам и столбцам к=1, 2, 3, 4, размер объекта наблюдения L×L=5×5. ДН задавалась экспоненциальной функцией с квадратичным показателем степени. Уровень γ задавался числом 5.

В таблице представлены оценки среднеквадратического отклонения (СКО) ошибки восстановления изображения, полученные сравнением моделируемого и восстановленного изображений на множестве реализаций эксперимента, в зависимости от величины шага k по строкам и столбцам при разном отношении сигнал-шум.

В первой строке значений СКО даны результаты, полученные предлагаемым способом. Вторая строка значений СКО соответствует прототипу, когда пропущенные строки и столбцы предварительно обнулялись, и восстановленное изображение представляло собой решетчатую матрицу. В третьей строке значений СКО представлены результаты, полученные оптимальным матричным способом восстановления, принимаемым за эталон, при совместной обработке двух матриц наблюдения Y₁ и Y₂, подобно предложенному способу, но на основе двумерной модели измерений (1).

Оценка разрешающей способности изображений осуществлялась нахождением минимального расстояния d_min, измеряемого в количестве строк между двумя объектами, при котором они четко различались на восстановленном изображении. Разрешение d_min=2, соответствующее 1/2 ширины ДН, достигалось: для предложенного способа при шаге k=2, для прототипа - при шаге k=1, для эталона - при шаге k=3.

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

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

Литература

1. Шарков Е.А. Радиотепловое дистанционное зондирование Земли: физические основы: в 2 т./ Т. 1. М.: ИКИ РАН, 2014. 544 с.

2. Пассивная радиолокация: методы обнаружения объектов / Под ред. Р.П. Быстрова и А.В. Соколова. М.: Радиотехника. 2008. 320 с.

3. Василенко Г.И., Тараторин A.M. Восстановление изображений. М.: Радио и связь, 1986. 304 с.

4. Патент RU 2 368 917 С1. Способ формирования изображений в многоканальных РТЛС и РЛС / В.К. Клочко. МПК: G01S 13/89. Приоритет 21.12.2007. Опубл.: 27.09.2009. Бюл. №27.

5. Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде MATLAB. М.: Техносфера, 2006. 616 с.