автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Непроекционная реконструкция изображения объектов при монохроматическом зондировании пространства и синтезе апертуры

доктора технических наук
Ющенко, Валерий Павлович
город
Новосибирск
год
2013
специальность ВАК РФ
05.13.18
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Непроекционная реконструкция изображения объектов при монохроматическом зондировании пространства и синтезе апертуры»

Автореферат диссертации по теме "Непроекционная реконструкция изображения объектов при монохроматическом зондировании пространства и синтезе апертуры"

На правах рукописи

Ющенко Валерий Павлович

НЕПРОЕКЦИОННАЯ РЕКОНСТРУКЦИЯ ИЗОБРАЖЕНИЯ ОБЪЕКТОВ ПРИ МОНОХРОМАТИЧЕСКОМ ЗОНДИРОВАНИИ ПРОСТРАНСТВА И СИНТЕЗЕ АПЕРТУРЫ

Специальность 05.13.18 -математическое моделирование, численные методы

и комплексы программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора технических наук

о <зщ 2014

005544280

Новосибирск — 2013

005544280

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет»

Научный консультант: Легкий Владимир Николаевич

доктор технических наук, с.н.с.

Официальные оппоненты: Воскобойников Юрий Евгеньевич

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

Меледии Владимир Гепрневич

доктор технических наук, профессор, «Институт оптико-электронных информационных технологий», (Новосибирск), генеральный директор

Петров Виктор Петрович

доктор технических наук,

Федеральное государственное образовательное бюджетное учреждение высшего профессионального образования «Сибирский государственный университет телекоммуникаций и информатики», профессор

Ведущая организация «Сибирский государственный ордена Трудового

Красного Знамени научно — исследовательский институт метрологии (ФГУП «СНИИМ») (Новосибирск)

Защита состоится " б " марта 2014г. в 10-00 на заседании диссертационного совета Д 003.005.01 при Федеральном государственном бюджетном учреждении науки Институте автоматики и электрометрии СО РАН по адресу: 630090, Новосибирск, проспект Академика Коптюга, 1.

С диссертацией можно ознакомиться в библиотеке ИАиЭ СО РАН.

Автореферат разослан "25" декабря 2013 г.

Ученый секретарь диссертационного совета Д.ф.-м. н.

Ильичев Л.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы диссертации. Поиск новых методов реконструкции изображения объектов, основанных на непроекционных способах сбора данных, расширяет возможности исследования объектов. Важное значение имеет способ сбора данных и связанные с ним методы извлечения информации, результатом обработки которой является реконструкция изображения объектов. Реконструированное изображение является наивысшим достижением информационно-измерительного процесса и даёт исчерпывающую информацию об объекте или о его внутренней структуре. Поэтому поиск новых методов реконструкции изображения - это всегда актуальная задача. В 1979 году Годфри Нью-болд Хаунсфилд и Аллан МакЛеод Кормак получили Нобелевскую премию за проекционные методы реконструкции изображения объектов. Как выяснилось позже, математические основы проекционной томографии были разработаны австрийским математиком Иоганном Радоном в 1917 году, а Годфри Ньюболд Хаунсфилд и Аллан Мак Леод Кормак переоткрыли изобретение И. Радона.

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

Существенный вклад в развитие радиоволновой проекционной томографии сделан в начале восьмидесятых годов американскими учёными Д. Мен-сой, Ш. Халеви и Г. Уэйдом в опубликованной работе «Применение методов когерентной доплеровской томографии для получения изображения на СВЧ». Ими использовался монохроматический зондирующий сигнал. Новизна их метода в том, что вместо просвечивания объектов они работали с отражённой от объекта волной. Они показали, что регистрация проекций по кругу, опоясывающему цель, позволяет восстановить картину цели и её положение в пространстве. Предложенный ими метод относится как к проекционной томографии, так и к круговому апертурному синтезу антенны. К сожалению, проекционные методы реконструкции не применимы в задачах, в которых требуется восстановить пролётные параметры ракеты «Поверхность-Воздух» или «Воздух- Воздух» около цели (рисунки 1 и 2). Это связано с тем, что проекции регистрируются при круговом обходе объекта, а ракеты движутся по прямолинейной траектории. Поэтому есть потребность в разработке непроекционных методов реконструкции изображения объектов, основанных на движении локатора по любой траектории, в том числе по прямолинейной. Из пролётных параметров наибольший интерес представляет промах ракеты, так как он позволяет оценить качество системы наведения и работу системы ближней локации. Под промахом подразумевается кратчайшее расстояние от центра цели до ракеты. Будем считать, что в момент реализации промаха 1 = 0.

Рисунок 1 - Пояснение пролетной си- Рисунок 2 — Пролетная ситуация в туации в системе координат, связанной системе координат ракеты, которую со станцией подсвета цели (СПЦ) требуется восстановить

Пролетная ситуация характеризуется следующими параметрами: а - промах, \!ц— вектор скорости цели, \'р- вектор скорости ракеты, \'от- вектор относительной скорости, а - угол атаки ракеты, В настоящее время промах определяют с помощью триангуляционного теодолита. Точность определения координат с его помощью описывается выражением (8£)0 « ш/,2 / ВЫ, где //-число элементов оптической матрицы, В - база триангуляционной системы, ю - поле зрения. Из этой формулы следует, что погрешность измерения дальности 8Ь возрастает с увеличением расстояния до цели прямо пропорционально квадрату ее удаленности Ь. Пусть погрешность измерения дальности при удаленности цели на 1 км составляла 20 м, тогда при удалении точки встречи ракеты с целью на 10 км от измерителя погрешность измерения дальности увеличится в 100 раз и составит значение 2000 м. При такой погрешности промах определить невозможно. Остаётся единственный путь извлечения информации о промахе - из траекторного доплеровского сигнала, переданного с борта ракеты по телеметрическому каналу.

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

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

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

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

Традиционно синтез апертуры антенны применяется в бортовых РЛС бокового обзора для картографирования земной поверхности. Синтезированная апертура — это условная апертура антенной решетки, которая синтезируется в результате движения по траектории летательного аппарата одним-единственным антенным элементом с изотропной диаграммой направленности. Единственный движущийся элемент, занимая ряд последовательных положений на траектории, ведет себя как полноценная антенная решетка с множеством элементов, расположенных вдоль апертуры (траектории) с некоторым шагом. Как реальная, так и синтезированная антенная решетка имеет узкую диаграмму направленности, угловая ширина которой определяется длиной апертуры антенны: АО = X / с/. Здесь Д0 - ширина диаграммы направленности в радианах, \ - длина волны, с/ - длина апертуры.

Синтез апертуры при малой дальности до объекта имеет много особенностей по сравнению с синтезом апертуры антенны на большом удалении. Эти особенности не позволяют перенести приемы синтеза апертуры, используемые при картографировании земли, на малую дальность. Главное препятствие, мешающее такому переносу, это проблема радиального разрешения.

Радиальное разрешение в системах картографирования достигается за счет импульсной или частотной модуляции зондирующего сигнала. Если при картографировании -земли обходятся радиальным разрешением от 1,5 м до 30 см, то в задачах с реконструкцией изображения цели такого разрешения бывает недостаточно. Для повышения разрешающей способности уменьшают длительность зондирующего импульса. В итоге пространственная протяженность зондирующих импульсов становится сравнимой с длиной волны X. При таком соотношении пространственной протяженности импульса и длины волны синтез апертуры становится невозможным. Более выгодным является путь синтеза апертуры при непрерывном зондирующем сигнале с линейной частотной модуляцией зондирующего сигнала.

Вернёмся к синтезу апертуры при монохроматическом зондировании , где разрешение достигается за счет фокусировки. Максимально достижимое разрешение при этом составляет: ДЯ = X/5. Благодаря фокусировке и высокому разрешению открывается возможность реконструкции изображения цели по

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

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

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

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

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

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

В работе Д.Менсы подробно рассмотрены вопросы разрешающей способности круговой синтезированной апертуры. О разрешающей способности синтезированной апертуры дают представление функции рассеяния точечных объектов при различных секторах синтеза апертуры. В работе Д.Менсы вместо принятого термина в оптике «функция рассеяния точки» введен термин «передаточная функция точки». Получены функции рассеяния точечных объектов при трехчастотном и пятичастотном зондировании. Однако это все результаты проекционного метода. Потому четвертой задачей является сравнение предлагаемых методов реконструкции с известными проекционными методами.

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

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

Задачи диссертационной работы

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

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

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

4. Сравнение предлагаемых и известных методов реконструкции.

5. Экспериментальная проверка предлагаемых методов и алгоритмов реконструкции изображения объектов.

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

Положения, выносимые на защиту

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

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

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

•рассчитывают опорные траекторные сигналы от точечных отражателей, расположенных на разном удалении от траектории движения локатора;

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

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

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

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

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

•расставляются опорные точки равномерно по всей площади траекторно-го круга в соответствии с декартовой сеткой и шагом не более У5, •определяются полярные координаты опорных точек, •отсчёты траекторных сигналов для опорных точек рассчитываются по формуле радиолокации, в которой расстояние между локатором и опорной точкой определяется в соответствии с геометрическими параметрами кругового перемещения локатора,

•отсчёты сигнала на круговой траектории берутся через равные угловые промежутки в полярной системе координат в соответствии с теоремой Котель-никова,

•начальная точка на круговой траектории для каждого опорного сигнала определяется полярным углом опорной точки,

•расчёт корреляционных функций между объектным и опорным!! траекторными сигналами производится по отсчётам, полученным в полярных координатах,

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

Научная новизна работы

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

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

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

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

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

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

Новизна результатов диссерта11ионной работы подтверждается патентом и авторским свидетельством на изобретение.

Практическая ценность работы

1. Предложенные методы реконструкции изображения цели по траектор-ному сигналу позволяют решить задачи визуализации пролётной ситуации и распознавания объектов.

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

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

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

5. Непроекционные методы реконструкции изображения объектов от-

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

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

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

Методы исследования Теоретические методы основаны:

- на корреляционном анализе натурных и модельных траекторных доплеров-ских сигналов;

- на приемах линейной фильтрации, быстрого преобразования Фурье и статистической радиотехники;

- на использовании математических методов реконструкции изображения по проекциям с применением прямого и обратного преобразования Радона, теоремы о центральном сечении, метода обратного проецирования и метода обращения по Фурье (метод Д.Менсы).

Экспериментальные методы основаны:

- на корреляционной обработке натурных траекторных сигналов и реконструкции по ним изображения объектов с помощью разработанных программ;

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

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

Реализация и внедрение результатов

В ОАО «НИИ Электронных приборов» (г. Новосибирск), в ОАО «Новосибирское производственное объединение «Луч», на ФКП «Новосибирский опытный завод измерительных приборов» результаты диссертации использованы для анализа натурных испытаний систем ближней локации.

Результаты диссертации используются также в научно-методическом комплексе подготовки инженеров по специальности 170100 «Боеприпасы и взрыватели», бакалавров и магистров направлению 220400 «Управление в технических системах» при разработке курсов «Обработка сигналов и изображений», «Математические методы в инженерных расчетах», «Радиоэлектронные информационно-управляющие системы».

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

1.Ющенко, В.П. Метод построения изображения объекта по доплеров-скому сигналу с помощью синтезированной апертуры // Актуальные проблемы электронного приборостроения: труды IV Международной конференции АПЭП-98 / НГТУ.- Новосибирск, 1998,- Т. 10 - С. 47-52.

2. Yushchenko, V. P. A Circular Aperture synthesis for Tomography // On Science and Technology: The 6th Russian-Korean International Symposium KORUS-2002. / At the - Novosibirsk State Technical University Russia, June 24-30, 2002-Vol. 1- P. 374-378.

3. Ющенко, В.П. Особенности в постановке задачи реконструкции при томографии с помощью синтезированной апертуры и ее решение // Информационные системы и технологии: международная науч.-техн. конференция ИСТ' 2003 / НГТУ,- Новосибирск, 22-25 апреля, 2003.- Т. 1- С. 158-162.

4. Ющенко, В.П. Сравнение различных методов и возможностей апер-турного синтеза в ближней локации // Конверсия, оборона, безопасность: всероссийская научно-практическая конференция / ПТУ — Пенза, Россия. 17—19 сентября, 2003.-С. 102-108.

5. Ющенко, В. П. Траекторный контроль с помощью синтезированной апертуры [Текст] / В. П. Ющенко, С. А. Лукьяненко, Н. М. Федуличева // Наука, промышленность, оборона : тр. VII всерос. науч.-техн. конф. - Новосибирск : Изд-во НГТУ, 2006. - С. 105-111.

Публикации. Основные материалы диссертации опубликованы в 24 работах. Из них 12 работ опубликованы в журналах, рекомендованных ВАК для публикации материалов докторских диссертаций, 1 работа - в периодическом издании, 11 статей в сборниках трудов конференций. Имеются 4 отчета по НИР, 1 патент на изобретение, 1 авторское свидетельство СССР на изобретение.

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

Из вышеизложенного следует, что круговой синтез апертуры в спектральной области методом проекций невозможно применить к задачам визуализации пролётной ситуации по следующим причинам:

1. Нельзя работать на малой дальности

2. Известный проекционный метод синтеза апертуры по кругу не применим к прямолинейным траекториям.

Во втором разделе раскрываются п. 1, 2, 3, 4 положений выносимых на защиту, излагаются основы методов синтеза апертуры антенны на малой дальности при монохроматическом зондирующем сигнале. Особенность их в том, что апертура синтезируется в области оригиналов, т. е. в пространственных координатах. Дается решение задачи реконструкции картины объекта с координатной привязкой ее к отрезку траектории, на котором синтезировалась апертура антенны, а именно предлагаются основы, непроекционного метода реконструкции изображения объектов и его приложения для визуализации пролётной ситуации [1, 2, 3,4, 5, 7,13, 14, 16, 18, 20].

Суть предлагаемого метода можно описать следующим образом. Геометрическая схема процесса сбора информации при визуализации пролётной ситуации и реконструкции изображения приведена на рисунке 3.

1-я точка объекта

Сигнал, отраженный от ;-й точки объекта с координатами (хьу,), принимаемый локатором при пролете мимо этой точки по прямолинейной траектории с промахом уh можно записать в виде

¿Ю = + , (1)

(vt + x,-CJ2) +у, X

где под промахом понимается кратчайшее расстояние по нормали от траектории до точки объекта, а(хпу,)— комплексная амплитуда рассеяния в i-й точке, Са -длина синтезированной апертуры, t - время, v — постоянная скорость перемещения приемника с излучателем относительно объекта, \ - длина волны.

Сигнал, отраженный от всего объекта, можно представить в виде

/(') = £/('), (2)

/=1

Введем понятие опорных сигналов. Опорными называются сигналы, отраженные от опорных точек, с координатами (0, /Ду), с единичной амплитудой рассеяния. Опорные точки расставляются в середине синтезированной апертуры

антенны на разном промахе с некоторым шагом Ду. Траекторный сигнал от опорной точки можно описать таким образом:

г.(0 =--^ехр[7—л/<У - С /2)2 + (1А\>)2 ], (3)

Ду — шаг по промаху, /—номер шага по промаху.

Тогда набор одномерных взаимных корреляционных функций между опорными сигналами g|(t) и сигналом ДО может быть представлен так:

Д(т,/Лу) = -\fifte, {(/ + т), 1Ау}Ж. (4)

^ о

Возьмем функцию /(г) для одноточечного объекта, т. е. в выражении (2) примем 7=1, а в выражении (1) положим х, =0, V,. Ф 0, я(.тр >', ) = 1 и тогда взаимная корреляционная функция примет вид:

, т

В,(т,1Лу) = Я(т,/Ду) = ~\тё, {(I + т ),/Ду}Л. (5)

^ о

Соотношение (5) представляет собой свертку двух функций времени и является аналогом выражения (б), описывающего реакцию линейного фильтра с импульсной характеристикой g¡(t) при воздействии на него сигнала /(/).

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

5,(т,/Д)0 = —5- | , (6)

С,/2г

где в(и)) = | /(Оехр[-М' + т)К', (7)

-Са/2г С„ /2г

К/(ш) = | gl(t)exv(-jut)dt. (8)

-с.п*

Здесь ш1 = 4тту/Х, 8(ш) - спектр зарегистрированного траекторного сигнала /(!), К,*(ы) - комплексно-сопряженный спектр опорной функции g¡(f) для промаха с номером /, К*(ш) является частотной характеристикой дисперсионного фильтра, перестраиваемого на разные промахи /Дг (дисперсионным считается фильтр, задержка сигнала в котором зависит от частоты).

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

Далее предлагается алгоритм обработки траекторного сигнала и построения изображения. Для построения алгоритма вычисления компоненты /¡(0 и

gí{t) взаимно корреляционной функции (5) представляются в дискретном виде. Шаг дискретизации определяется в соответствии с теоремой Котельникова. Предположим, что исследуемый объект состоит из I точек и траекторный сигнал, отраженный от г'-й точки, в дискретном виде запишется так:

г ,4-к

s,(k) =----7[cos(— J(v(kAt) + x - С/2 f + y;) -

' (v{kAt) + x - Ca / 2) + y. VXV V ' * ; У,)

-} А/) + X - С / 2)2 + у,2)].

Сигнал, отраженный от всех I точек, можно представить суммой

(9)

(10)

Выражения, записанные в дискретном виде, легко программируются. Траекторный сигнал, от трехточечного объекта, рассчитанный по формулам (9) и (10), представлен на рисунках 4, а, а на рисунке 4, б представлен его спектр.

0

а) ^

0,01

б)

о 500 1000 1500 к

Рисунок 4 — Результаты моделирования: а) траекторный сигнал, от трехточечного объекта; б) его спектр

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

g'(k) = ТТ^ла Jcos(^>(fcA0-Cfl/2)2 + ^) +

(v(kAt)-Ca/2) + у, X

+j sin(^(v(kAt) - С / 2f + yf )].

А

В выражении (11) через I обозначается номер опорного сигнала. Всего L опорных сигналов. Каждому I соответствует задаваемое значение промаха у,. Отсчеты промаха у, берутся через равные промежутки Ау.

С помощью прямого преобразования Фурье отсчеты траекторного доп-леровского сигнала (10) превращаются в отсчеты спектра S(n):

S(») = f>(A:)exp(-j^to), TÍ N

где и = 1,2= 1,2,3,...,ТУ.

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

к;(и) = Хй(*)вхр(7^*и), (13)

{-=1 -¡у

где п = 1,2,3,....,Л', Аг = 1,2,3,....,Я.

Далее вычисляется дискретная свертка двух функций, т. е. находятся отсчеты взаимной корреляционной функции.

1 2тт

5ДА:) = —^8(«)К;(«)ехрО—Ап), (14)

реальная (б) где и = 1,2,3,.....к = 1,2,3,.....

Полненные отсчеты взаимно корреляционной функции для одного опорного сигнала g/ соответствуют одной строке разворачиваемого изображения. Число строк изображения равно числу опорных сигналов. Каждый опорный сигнал соответствует конкретному промаху. Значения промаха следуют один за другим от минимального до максимального с некоторым шагом. При этом для компенсации радиального ослабления сигнала в соответствии с формулой локации, амплитуда отсчетов опорного сигнала умножается на величину промаха в четвертой степени. Таким образом, вводится автоматическая регулировка усиления принимаемого сигнала. Далее отсчеты всех взаимно корреляционных функций нормируются по отношению к наибольшему значению из всего набора корреляционных функций. После этого модуль каждого комплексного числа перекодируется в яркостной или цветовой пиксель в соответствии с его величиной по отношению к найденному максимуму. Далее пиксели располагаются по строкам в декартовой системе координат в том же порядке, что и числа корреляционной матрицы. В результате такой развертки воспроизводится картина внутренней структуры объекта в сечении.

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

а)

Рисунок 5 - Диаграмма направленности антенны локатора: идеальная (,а),

реальная (б)

Диаграмма направленности антенны для реальных локаторов приведена на рисунке 5, 6. Результаты модельной реконструкции, которые приводятся ниже, соответствуют идеальной диаграмме направленности (рисунок 5, а). Поверхность, построенную из набора корреляционных функций, называют функцией рассеяния. Этот термин заимствован из оптики. Функция рассеяния точки приведена на рисунке 6, а. На рисунке 6, б дан вид её сверху. Подобные поверхности у Д.Менсы называются передаточными функциями.

Функция рассеяния точки характеризует разрешающую способность метода реконструкции изображения объекта. Как видно из рисунка 6 представленный метод обладает высоким разрешением по координате вдоль траектории. Максимальное разрешение характеризуется «талией восьмерки» (рисунок 6, б).

Р 200

100

0 100 200 к

а) б) а) б)

Рисунок 6 - Результат реконструкции Рисунок 7 - Виды спектра: спектр изображения: функция рассеяния точки функции рассеяния точки (а), его (а), изолинии или вид функции рассея- вид в плане (б)

ния сверху (б)

Следует подчеркнуть отличия в алгоритмах проекционной и непроекционной реконструкции. При проекционном методе регистрируют отсчёты спектра по кольцу. Затем по нему реконструируется функция рассеяния точки. Для этого берут обратное преобразование Фурье от кольцевого спектра. Здесь предлагается обратный алгоритм. Регистрируется траекторный сигнал при движении локатора по прямой траектории. С помощью метода взаимных корреляций производится многократная свёртка зарегистрированного траекторного сигнала с каждым опорным траеьсгорным сигналом и из полученного набора взаимных корреляционных функций строится поверхность в декартовой системе координат, которая даёт функцию рассеяния точки. Прямое двумерное преобразование Фурье от этой функции (15) даёт кольцевой спектр.

В дискретном виде эту математическую операцию можно представить следующим выражением:

НХ,¥) = Л- | | В(х,у)ехр[-ДхЛГ - уГ^Ыу.

(15)

1р)].

(16)

Кольцевой спектр (рисунок 7, а) функции рассеяния (рисунок 6, а) получен в результате двумерного преобразования Фурье (16) от функции рассеяния

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

На рисунке 8 показан результат реконструкции изображения трехточечного объекта. В качестве регистрационных данных взят сигнал, (рисунок 4а).

промах

О

промах

а) б)

Рисунок 8 - Результат реконструкции изображения: функция рассеяния трех точек (а), томограмма трехточечного объекта (б).

На рисунке 9 показан полукольцевой спектр пространственных частот трехточечного объекта.

И*,«»

1 ■ и

А II 1 №

1я 1 8 -.

200

100

\40'

0

а)

100

б)

200

0.328Х

Рисунок 10 - Многоточечный объект в виде тупого угла

Рисунок 9- Виды спектра: спектр пространственных частот трехточечного объекта (а), вид этого спектра сверху (б) На томограмме трехточечного объекта рисунок 8, б хорошо виден интерференционный узор. Интерференция проявляется также в спектральной области рисунок 9, а в виде изрезанноети спектральной функции на полукольце.

Рассмотрим реконструкцию изображения многоточечного объекта, состоящего из ряда точек на ломаной линии в виде тупого угла. Всего тридцать точек, расположенных через 0,328\ (рисунок 10). На рисунке 11 представлено реконструированное изображение этого многоточечного объекта. Она имеет вид двух пересекающихся лучей. Из рисунка 11 можно видеть, что изображение объекта представляет собой диаграмму обратного рассеяния волн, которая образовалась в результате интерференционных явлений. Далее проводится сопоставление результатов модельной реконструкции точечных объектов предлагаемым методом, с результатами реконструкции проекционным методом.

а) б)

Рисунок 11 - Изображение много- Рисунок 12 - Функции рассеяния точки точечного объекта из ломаной ли- полученные: а) методом Д.Менсы;

нии в виде тупого угла б) непроекционным методом

Выражение (17), заимствованное из описания проекционного метода реконструкции по спектру на кольце, позволяет рассчитать и построить функцию рассеяния точки (рисунок 12, а), которая даёт лепестки в виде концентрических окружностей вокруг главного лепестка.

С(|г|,ф) = |4 Т ехр у—И ««(ф - 9)]Ле (17)

И о Х

На рисунке 12, б приведена функция рассеяния точки, полученная методом корреляций при синтезе апертуры по прямолинейной траектории. Из сравнения этих функций видно, что прямолинейный синтез апертуры не дает артефактов в виде концентрических колец, что является его преимуществом.

а) б)

Рисунок 13 - Передаточные функции точки, полученные двумя разными методами при секторе синтеза апертуры 60°: результат Менсы (а)\ результат предлагаемого непроекционного метода (б) Сравним передаточные функции точки, полученные по спектру на кольце (рисунок 13, а) проекционным методом с передаточными функциями, полученными синтезом апертуры при движении локатора по прямолинейной траектории при секторе синтеза апертуры 60° (рисунок 13, б). Результаты проекционной и непроекционной реконструкции совпадают.

Дальнейшее моделирование показало, что при сужении сектора синтеза апертуры происходит ухудшение разрешающей способности по обеим коорди-

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

Далее в разделе 2.5 раскрывается предложенный путь устранения интерференции. Он основан на многоракурсном синтезе апертур антенны [4, 7, 9, 10]. Функция рассеяния точки для идеального процесса реконструкции должна иметь вид 6 -импульса. Такую функцию можно получить, синтезируя апертуру в двух взаимно перпендикулярных направлениях (рисунок 14).

В

Рисунок 14 - Передаточная функция точки при двуракурсном синтезе апертуры

0 к Рисунок 15 - Результаты реконструкции шестнадцатиточечного объекта при одноракурсном синтезе апертуры

Рассмотрим результаты реконструкции изображения шестнадцатиточечного объекта (рисунок 16, а). На рисунке 15 представлен результат восстановления изображения этого объекта по одному траекгорному доплеровскому сигналу. Интерференционная картина на рисунке 15 не позволяет распознать объект.

_0,975 X

Г5ГбТх

• • • « •

. \aS3L

%

Рисунок 17 - Результат реконструкции при четырехра-курсном синтезе апертуры

а) б)

Рисунок 16 — Шестнадцатиточечный объект: а) вид объекта; б) результат реконструкции с помощью двухракурсного синтеза апертуры

На рисунке 16, б представлен результат восстановления этого же объекта, по двум траекторным доплеровским сигналам при ортогональном расположении траекторий синтеза апертур.

Результат реконструкции ещё улучшится, если повернуть координатный базис синтеза апертур в двух взаимно перпендикулярных направлениях на 45" и еще раз сделать синтез апертур в двух взаимно перпендикулярных направлениях, а затем перемножить совпавшие элементы матриц (рисунок 17). Ослабление интерференции объясняется тем, что максимумы функций рассеяния одно-

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

В результате синтеза апертур в двух направлениях получены две матрицы изображения размером 256 х 256 элементов. Если расстояние между элементами в обеих матрицах одинаковое, т. е. шаг по промаху равен шагу вдоль траектории Арг = г>А/, тогда при наложении развернутых на 90° матриц элементы строк совпадают с элементами столбцов и происходит полное совпадение одноименных точек. Совпавшие элементы перемножаются. Другой результат получается, если матрицы развернуть на 45 (рисунок 18).

Чтобы у матриц совпало большее число элементов уменьшают расстояние между элементами какой либо одной матрицы в л/2 раз. Для этого при регистрации шаг по промаху и шаг вдоль траектории уменьшают в -¡2 раз.

£< Ш я @ я ►!< * * ►!< Я

ьь

У___

Рисунок 18 - Совмещение элементов развернутых матриц На участке полного перекрытия матриц совпадение элементов достигает 50 %. При этом перемножаются только совпадающие элементы матриц и отбрасываются не совпадающие элементы. Объект на развернутой матрице изображения будет занимать такое же положение, что и при не развернутой матрице. На первый взгляд разворот матриц с перемножением элементов похож на известный алгоритм обратного проецирования. Однако есть существенное отличие. Главное отличие в том, что при обратном проецировании в выбранной точке объекта суммируются не значения лучевых функций, которые они имеют в точке пересечения лучей, как это делается в предлагаемом методе, а значения интегральных сумм вдоль пересекающихся лучей. Другими словами, в точке пересечения лучей суммируются их проекции. Если представить, что объект состоит из одной точки, то не точка проецируется куда - то по лучам, а наоборот проекции по этим лучам собираются в точке. Отсюда и название «Метод обратного проецирования».

Далее раскрываются п. 1, 2, б положения, выносимого на защиту. Эта часть раздела посвящена опробованию корреляционного метода синтеза апертуры на задачах визуализации пролётной ситуации с использованием натурных сигналов. Визуализация подразумевает восстановление таких параметров пролета, как промах, размеры объекта, положение объекта относительно конечного отрезка траектории. В качестве объекта для исследования взят конус (рисунок 19).

О

-10в

—]

ИГ 1<ИШ

1200

500

1000

1500

Рисунок 19 — Объекта для исследования

Рисунок 20 - Траекторный сигнал, отраженный от объекта

Сигнал от объекта регистрируется доплеровским приемником, работающим с монохроматическим сигналом на длине волны «3 см при перемещении локатора перпендикулярно продольной оси конуса с некоторым промахом.

Траекторные сигналы получены методом проездов по схеме (рисунок 21). Исследуемый объект, в виде конуса

"О-и

Рисунок 21 - Схема проезда локатора мимо исследуемого объекта В первом эксперименте скорость движения локатора была уот =2 м/с. Промах составлял 0,5 м. Отсчеты траекторного сигнала взяты с интервалом 1 мс. Экспериментальный сигнал представлен на рисунке 20. На рисунке 22 представлены результаты реконструкции изображения сечения конуса по этому сигналу, то есть функция рассеяния исследуемого объекта и его реконструированное изображение.

Результаты отображаются в координатах промаха и расстояния вдоль траектории относительного движения локатора и объекта. Поле реконструкции изображения объекта (рисунок 22, б) охватывает пространство размером 4x4 м. Цель отображается в виде сгустка точек повышенной плотности отражения. Как видно из рисунка 22. б область с повышенной плотностью соответствует картине усреднённого поперечного сечения конуса, напоминающего круг. Лепестки диаграммы обратного рассеяния объекта (рисунок 22, б) ориентированны определенным образом. Они веером исходят из центра круга и тем самым помогают локализовать цель. Таким образом, область с повышенной плотностью указывает положение цели в окрестности траектории (горизонтальная ось) и определяет её промах (вертикальная ось).

промах

3.60269-10

-4

-4

i 255 170 85 промах

а)

2.40229-10"

0.20189-10

Р ■ 1г

470

4 м

расстояние по траектории б)

Рисунок 22 - Результат реконструкции: ^функция рассеяния исследуемого объекта; б) реконструированное изображение объекта в сечении

Из реконструированного изображения (рисунок 22, б) видно, что цель находится как раз на расстоянии 0,5м от траектории. Статистика показывает, что погрешность измерения промаха составляет примерно 20% от истинного значения промаха. Интерференция, то есть ДОР одновременно мешает и помогает локализовать цель. Веерный характер расположения лепестков ДОР своим схождением указывает место положения цели. В тоже время лепестки ДОР в принципе препятствуют локализации цели. Оценка погрешности зависит также от того, как на картинке оператор ограничит место локализации цели. Если исходить из практических соображений, то погрешность измерения 20% или 50% лучше, чем полное незнание промаха.

В третьем разделе раскрываются п. 1, 3, 5, 6 положений выносимых на защиту, а именно, теоретические аспекты синтеза апертуры по круговой и квадратной траектории [6, 8]. Формулируется постановка задачи в области оригиналов, вводятся основные понятия и определения, дается математическое описание метода взаимной корреляции при круговом синтезе апертуры. Синтез апертуры по круговой и квадратной траектории интересен не только с общетеоретических позиций. Его применение ориентировано на системы радиовидения, томографию, охрану объектов, досмотр в аэропортах.

Пусть локатор с изотропной диаграммой движется по круговой траектории, облучая монохроматическим сигналом отражающую точку, смещенную относительно центра на расстояние Ьк, и регистрирует траекторный сигнал. Дискретные отсчеты сигнала регистрируются с угловым шагом Дф =(2ю'Ы), где N - число отсчетов, ай9 - радиус круговой траектории, по которой синтезируется апертура. Можно показать, что в этом случае расстояние от реконструируемой точки до локатора будет изменяться по закону

Д(0 = л/й02+6/-2^64со8(/Дф), где Ьк- радиальное удаление к-й точки объекта от центра синтезирования апертуры, I- номер отсчета доплеровского сигнала на круговой траектории от 0 до Лг,

О < г'Дф < 2тт - текущий угол между радиусами Н0 и Ьк при регистрации отсчетов. Эти границы изменения угла соответствуют обходу реконструируемого объекта по круговой траектории в полярных координатах.

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

«(0 = соБ^рО). (17)

Более точная модель, требует учета радиального затухания, которое будет изменяться обратно пропорционально второй степени й(г).

5(0 = (1/й1( 0)СО8ф(/), (18)

... 4тгД(/) х

где ф(г) =-—, X - длина волны.

X

В случае большого количества точек, отраженный от них сигнал $(/), в соответствии с принципом суперпозиции можно рассчитать по формуле:

*(0=2Х(0, (19)

»1=1

где М - общее число реконструируемых точек, т - порядковый номер реконструируемой точки объекта, - траекторный доплеровский сигнал, отраженный от точки с номером т и рассчитанный по формуле (17) или (18).

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

2 N 2К 2т

Вк(ч) = 1>(»>,(»• + *,) = Хв(«)КД»1)ехр(у—1,и) при ¿,=1,2.....2Д (20)

/=1 п=1

N

при и=1,2.....2jV,

2N .it

где S(«) = J^s(i)exp(-j—m) при;'=1,2......27V, (21)

/=1 N

s(i) = 0 при i=N+\,N+2...2N,

2 N

КД«) = Х^(0ехр(-у—ш) при i=N+l,N+2...2N, (22)

í=i N

sk(0 = 0 при /=1,2.....N.

В выражениях (20) (21), (22) приняты следующие обозначения: S(w)- отсчеты спектра сканируемого сигнала (19), Кк(п) - отсчеты частотной характеристики дисперсионного фильтра, настроенного на дальность Ьк, sk(i) - отсчеты опорного сигнала, рассчитанные по формуле (17) или (18) для одной точки, удаленной от центра синтезирования апертуры на Ьк, п - номера отсчетов частоты, i — номера отсчетов сигнала, i¡ — номера отсчетов взаимной корреляционной функции. Таким образом, процесс синтеза апертуры и построения изобра-

жения сводится к многократной свертке траекторного сигнала вида (19) с опорными сигналами вида (17) или (18), рассчитанными при различных расстояниях от центра 6А., и развертке центральных отсчетов взаимной корреляционной функции, перекодированных в тоновые или цветовые пиксели, в декартовых координатах. Каждый центральный отсчет взаимной корреляционной функции соответствует полярным координатам 6аи г[>а. Для перехода от декартовой к полярной сетке координат, используются формулы:

6а =Д5л/г"2 , (23)

=агс*8т> (24)

к

где /=1,2......256 - номер отсчета по оси абсцисс, к-12 ...... 256 - номер отсчета по оси ординат, Д.? - шаг по координатам

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

а) б)

Рисунок 23 - Круговой синтез апертуры: а) функция рассеяния центральной точки; б) функция рассеяния точки, удаленной от центра на 14Х

На рисунке 23, б представлена функция рассеяния точки, удаленной на расстояние 14\от центра синтезирования апертуры. Радиус синтеза апертуры [{„ составлял 25Х. Новым оказался результат, состоящий в том, что для точки, расположенной в центре круговой апертуры и дающей траекторный сигнал с нулевым доплеровским спектром, получается такая же функция рассеяния, как и для точки, удаленной от центра круга, дающей сигнал с широким доплеровским спектром (рисунок 23, а, б). Отсюда следует: 1) нет зависимости разрешающей способности от ширины спектра траекторного сигнала, 2) круговая апертура дает одинаковое разрешение внутри круговой области для всех точек.

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

Рассмотрим точечный объект расположенный в начале координат, функция рассеяния которого приведена на рисунке 23, а. Пространственный спектр этой функции представлен на рисунке 24, а. На рисунке 24, б представлен этот же спектр в плане. На рисунке 24, в представлена фазо-частотная характери-

В(/,/)

стика кольцевого спектра ф(/г--/г) = агс1ё(-~——), где А(/х,/) - действи-

А(/х>/у)

тельная часть спектральной функции, а В(/х, /) - её мнимая часть.

к

|8(Аг,«)| а *

_1. Й§а : и

■п

О п

а) б) в)

Рисунок 24 - Спектральные характеристики: а) кольцевой спектр передаточной функции точки; б) вид кольцевого спектра в плане; в) фазовый спектр

На рисунке 25, а изображен трехточечный объект, реконструированный по траекторному сигналу. Точки удалены от центра синтезирования апертуры на 14Х. Расстояние между точками 2.8Х.

а) б) в)

Рисунок 25 - Круговой синтез апертуры: а) реконструированный по траекторному сигналу трехточечный объект; б) его спектр; в) вид спектра в плане

Как видно из рисунка 25, б спектральная функция на кольце сильно изрезана. Эта изрезанность есть результат интерференционных проявлений.

Четвертый раздел посвящен модельной томографии гомогенных областей с включенной неоднородностью [10, 11]. Точки модели гомогенной области распределены равномерно с шагом Х/7 и имеют одинаковую амплитуду рассеяния за исключением точек включенной неоднородности.

а) б)

Рисунок 26 - Результаты моделирования томографического процесса: а) гомогенная область с пустотой внутри; б) гомогенная область с неоднородностью

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

в пять раз модуль комплексной амплитуды рассеяния точек охватывающей области. Апертуры синтезировались по двум ортогональным траекториям. На рисунке 27 даны характеристики контрастной чувствительности.

1,5

2,5 5 10

Рисунок 27 - Характеристики контрастной чувствительности

Под контрастной чувствительностью подразумевается различимость двух гомогенных областей с незначительным отличием коэффициентов рассеяния волн их точками. Если модули комплексных амплитуд рассеяния точек внутренней области отличаются от модулей комплексных амплитуд рассеяния точек охватывающей области в 1,05, в 1,1, в 1,3, в 1,5, в 2,5, в 5 и в 10 раз, то ряд томограмм с указанными отличиями гомогенных областей будет выглядеть так, как показано на рисунке 27. Заключительная часть четвертого раздела посвящена исследованию устойчивости предложенного метода томографии к шумам.

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

а)

б)

Им&ж&Ж

ШШ

ШИН

б)

Рисунок 29 - Результат реконструкции изображения: а) до НЧ коррекции; б) после НЧ коррекции

Рисунок 28 - Данные эксперимента: а) траекторный сигнал от стержня; б) спектр сигнала

На рисунке 29 представлено реконструированное изображение металлического стержня в сечении, полученное в результате обработки траекторного сигнала, (рисунк 28, а). Сечение стержня отображается в виде сгустка повышенной плотности, прижатого к оси промаха. Изображение смещено влево, так как влево сдвинута реализация траекторного сигнала (рисунок 28, а).

Вторая половина раздела посвящена экспериментальной томографии на электромагнитных волнах при частоте 8-10 ГГц. В разделе выводится мате-

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

0.002 0.001

0

Рисунок 30 - Спектр сигнала от биологического объекта

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

Рисунок 31 - Экспериментальная ус- Рисунок 32 -Изображение пальца в тановка для томографирования сечении

Чтобы осуществить корреляционную обработку высокочастотной компоненты сигнала, несущей информацию о внутренней структуре биологического объекта, необходимо сформировать опорные сигналы при длине волны, с которой она распространяется в тканях биологического объекта. Эта волна в 6—7 раз короче волны в воздушной среде и равна примерно 4 мм [12]. Корреляционная обработка сигнала от внутренних тканей опорными сигналами, сформированными при X = 4 мм, позволяет получить разрешение равное 0,8 мм. Такая разрешающая способность вполне приемлема для медицинских исследований.

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

Часть спектра от внутренних тканей

100 200 300 400 500 600

В заключении сформулированы результаты диссертационной работы:

1. Разработаны научные основы непроекционных методов реконструкции изображения объектов по траекторному сигналу при локации объекта монохроматическим сигналом с синтезом апертуры антенны [1, 3, 5, 9, П. 12, 21].

2. Разработан алгоритм синтеза апертуры антенны в пространстве оригиналов на базе применения взаимных корреляционных функций, полученных сверткой объектного траекторного сигнала с набором опорных сигналов от точечных объектов, распределенных на разном удалении от траектории движения, с последующим формированием из набора одномерных корреляционных функций двумерных матриц изображения цели. Показано, что предложенный способ обеспечивает реконструкцию изображения объектов при расположении локатора как в зоне Френеля, так и в зоне Фраунгофера по отношению к апертуре объекта [1,4, 6, 11, 21, 22, 23].

3. Разработанные алгоритмы реконструкции изображения и программное обеспечение реализованы в методике натурных измерений пролетных параметров локатора у цели по реальному траекторному сигналу [1, 3, 4, 7, 8, 10,11].

4. Предложен способ реконструкции изображения объекта при синтезе апертуры по замкнутой траектории (круговой или квадратной) в пространстве оригиналов при монохроматическом зондировании пространства. Разработан алгоритм синтеза апертуры антенны по замкнутой траектории в пространстве оригиналов в полярных координатах с переходом к двумерным матрицам изображения объектов в декартовых координатах. Разработанный алгоритм способа доведен до программной реализации [4, 6,11].

5. С помощью разработанных программно-алгоритмических комплексов исследованы особенности апертурного синтеза на малой дальности от объекта при прямолинейной, круговой или квадратной траекториях носителя. Впервые получены функции рассеяния реконструированных объектов, их пространственные спектры на кольце [1, 2,4, 5, б, 8, 10,11,13, 17, 19, 20].

6. Доказаны следующие свойства изображений точечных объектов:

- отсутствие колец Эйри при синтезе апертуры по прямолинейной траектории [1, 3, 5, 7,8,13];

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

- совпадение передаточных функций точечных объектов, полученных предложенными методами, с передаточными функциями, полученными проекционными методами [5].

7. Предложен метод устранения интерференции при использовании монохроматического зондирующего сигнала с применением многоракурсного синтеза апертуры антенны [2, 5, 7, 8, 10].

8. Разработаны приёмы для обработки наложенных и развернутых на произвольный угол квадратных матриц изображения. Показано, что эта операция неэквивалентна операции обратного проецирования [5].

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

нородностью [8].

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

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

^.Экспериментально показано, что отраженный от биологического объекта сигнал состоит из суммы двух сигналов. Первый из них (с узким спектром в районе нулевой частоты) образован отражениями от границы двух сред (воздух - биологическая ткань). Второй сигнал с широким спектром и малой спектральной плотностью образован отражениями от внутренних тканей [10].

13.Получена картина сечения биологического объекта за счет корреляционной обработки сигнала от внутренних тканей, опорными сигналами, сформированными при длине волны в тканях биологического объекта \=4 мм. Достигнута разрешающая способность 0,8 мм при монохроматическом зондировании на частоте 8,3 ГГц, что приемлемо для использования в медицине [10].

Список опубликованных работ по теме диссертации Публикации в изданиях, входящих в перечень ВАК

1. Ющенко, В. П. Метод построения изображения объекта по доплеровскому сигналу с помощью синтезированной апертуры [Текст] / В. П. Ющенко // Науч. вест. - Новосибирск : Изд-во НГТУ, 1999. - № 3. - С. 36-49.

2. Ющенко, В. П. Интерференцию можно устранить [Текст] / В. П. Ющенко // Науч. вест. - Новосибирск : Изд-во НГТУ, 2000. - № 1 (8). - С. 185-192.

3. Ющенко, В. П. Доплеровское сканирование структуры объекта с помощью синтезированной апертуры [Текст] / В. П. Ющенко // Приборы и техника эксперимента. - 2001. - № 2. - С. 4И5.

4. Ющенко, В. П. Круговой апертурный синтез для целей томографии [Текст] / В. П. Ющенко // Автометрия. - 2002. - Т. 38, № 6. - С. 28-33.

5. Ющенко, В. П. Сопоставление двух методов когерентной томографии [Текст] / В. П. Ющенко // Радиотехника и электроника. - 2004. - № 2. — С. 196-205.

6. Ющенко, В. П. Синтезирование апертуры по квадратной траектории для когерентной томографии [Текст] / В. П. Ющенко // Изв. вузов. Сер. Радиоэлектроника. - 2004. - № 2. - С. 18-27.

7. Ющенко, В. П. Устранение интерференции при томографировании с помощью апертурного синтеза [Текст] / В. П. Ющенко // Науч. вестн. - Новосибирск : Изд-во НГТУ, 2004. - № 1 (16). - С. 181-188.

8. Ющенко, В. П. Восстановление внутренней структуры гомогенных объектов с включенной неоднородностью [Текст] / В. П. Ющенко // Автометрия. -2005.-Т. 41,№ 1.-С. 50-57.

9. Ющенко, В. П. Устойчивость к шумам метода томографирования с помощью синтезированной апертуры [Текст] / В. П. Ющенко // Изв. вузов России. Сер. Радиоэлектроника. - 2005. - № 3. - С. 3-9.

10. Монохромная СВЧ томтографня с синтезом апертуры методом взаимных корреляций [Текст] / В. П. Ющенко, С. А. Литвиненко, В. Б. Ромодин, Л. В. Щебалкова // Автометрия. - 2007. - Т. 43, № 6. - С. 28-38.

11. Yushchenko W. P. Object Image Reconstruction by Using Trajectory Doppler Signal in The Process of Monochromatic Probing and Circular Aperture Synthesis of Antenna. // Electrical and Electronic Engineering. P-ISSN; 2162-9455, e-ISSN: 2162-8459, 2012, 2(2): 1-6, DOI: 10.5923/j.eee. 20120202.01, http://www.sapub.org/journal/archive.aspx?journalid=1008

12. Ющенко, В. П. Обратимость непроекционной реконструкции изображения при монохроматическом зондировании и синтезе апертуры [Текст] / В. П. Ющенко // Сборник научных трудов SWorld. Материалы международной научно-практической конференции «Научные исследования и иж практическое применение. Современное состояние и пути развития 2012» - выпуск 3. Том 11. - Одесса: КУПРИЕНКО. 2012. ЦИТ: 312-119 - С. 10-16.

Публикации в других научных изданиях

13. Ющенко, В. П. Два метода когерентной томографии [Текст] / В. П. Ющенко // Сб. науч. тр. - Новосибирск : Изд-во НГТУ, 2001. - № 3 (25). - С. 2430.

Публикации в материалах конференций

14. Ющенко, В. П. Исследование диаграмм обратного рассеяния волн на малой дальности с помощью синтезированной апертуры [Текст] / В. П. Ющенко // 4-е Уткинские чтения. - Санкт-Петербург. : БГТУ "Военмех", 2009. - Т. 2. - С. 58-62

15. Ющенко, В. П. Помехоустойчивость систем ближней локации с синтезированной апертурой антенны [Текст] / В. П. Ющенко, В.В. Бычков // Проблемы автоматизации и управления в технических системах: тр. международной науч.-техн. конф. - Пенза 2011.

16. Ющенко, В. П. Рациональное построение доплеровских систем ближней локации угломерного типа с синтезом апертуры [Текст] / В. П. Ющенко, В.В. Бычков // Проблемы автоматизации и управления в технических системах: тр. международной науч.-техн. конф. - Пенза 2011.

17. Ющенко, В. П. Интроскопия с помощью синтеза апертуры для поиска предметов скрытых под одеждой [Текст] / Н. М. Федуличева, С. А. Лукья-ненко, В. П. Ющенко // Наука, промышленность, оборона : тр. VII всерос. науч.-техн. конф. - Новосибирск : Изд-во НГТУ, 2006. - С. 87-93.

18. Ющенко, В. П. Траекторньш контроль с помощью синтезированной апертуры [Текст] / В. П. Ющенко, С. А. Лукьяненко, Н. М. Федуличева // Наука, промышленность, оборона : тр. УП всерос. науч.-техн. конф. — Новосибирск : Изд-во НГТУ, 2006. - С. 105-111.

19. Ющенко, В. П. Применение синтеза апертуры в ближней локации [Текст] / В. П. Ющенко // Вторые Рдултовские чтения : науч.-техн. конф. - СПб. : Изд-во БГТУ, 2008. - С. 52-58.

20. Ющенко, В. П. Прямолинейный апертурный синтез для определения пролетной ситуации [Текст] / В. П. Ющенко // Актуальные проблемы развития

автономных информационных и управляющих систем : сб. тр. V всерос. науч.-техн. конф. - СПб.: Изд-во Военмех, 2003. - Т. 2. - 5 с.

21. Ющенко, В. П. Сравнение различных методов и возможностей апертурно-го синтеза в ближней локации [Текст] / В. П. Ющенко // Конверсия, оборона, безопасность : всерос. науч.-практ. конф. - Пенза : Изд-во ИГУ, 2003.-С. 102-108

22. Ющенко, В. П. Метод построения изображения объекта по доплеровскому сигналу с помощью синтезированной апертуры [Текст] / В. П. Ющенко // АПЭП-98: Актуальные проблемы электронного приборостроения : тр. IV междунар. конф. - Новосибирск, 1998. - Т. 10. - С. 47-52

23. Ющенко, В. П. Особенности в постановке задачи реконструкции при томографии с помощью синтезированной апертуры и ее решение [Текст] / В. П. Ющенко // ИСГ2003: Информационные системы и технологии : междунар. науч.-техн. конф. - Новосибирск : Изд-во НГТУ, 2003. - Т. 1. - С. 158-162.

24. Yushchenko, V. P. A Circular Aperture synthesis for Tomography // OnScience and Technology: The 6th Russian-Korean International Symposium KORUS22 2002. / At the - Novosibirsk State Technical University Russia, June 24-30, 2002-Vol. l.-P. 374-378.

Патенты и авторские свидетельства

25. Пат. 2066060 Российская Федерация, МПК6 G01S 13/90. Способ картографирования с помощью синтезированной апертуры [Текст] / Ющенко В. П ; заявитель и патентообладатель Ющенко Валерий Павлович. - № 93008773/09 ; заявл. 16.02.1993 ; опубл. 27.08.1996, Бюл. № 24 - 10 с.

26. А. с. 325757 СССР, МКИ3 6 G 01 S 13/04. Способ картографирования с помощью синтезированной апертуры [Текст] / В. П. Ющенко. — № 4513396 ; заявл. 24.04.89 ; опубл. 5.05.91. - 10 с.

отчет о НИР

27.Исследование возможностей апертурного синтеза для целей томографии [Текст] : отчет о НИР / Новосиб. электротехн. ин-т ; рук. В. П. Ющенко ; ис-полн. : А. И. Харламов, М. И. Гольберг, П. А. Дмитриенко. - Новосибирск, 1992. - 58 с. - № гос. регистрации 01920005752

28.Исследование помехоустойчивости автономных систем управления: отчет о НИР Шифр «Нормаль» / НГТУ; рук. В.Н.Легкий; исполнители: Шумейко В.А., Ющенко В.П., Борисова И.В. - Новосибирск, 2012г. — 188с.

29. Исследование по воздействию излучения на объекты: отчет о НИР / ОАО «НИИЭП»; рук. В.Н.Легкий; исполнители: Колышев Д.В., Ющенко В.П., Борисова И.В. - Новосибирск, 2012г. - 196с.

30. Исследование алгоритмов и систем признаков детального распознавания воздушных целей: отчет о НИР / НП СибНЦ «Безопасность»; рук. В.Н.Легкий; исполнители: Ющенко В.П., Борисова И.В., Санков О.В. - Новосибирск, 2012г. - 196с.

Отпечатано в типографии Новосибирского Государственного технического университета 630073, г.Новосибирск, пр. К. Маркса, 20, Тел./факс (383) 346-08-57 Формат 60 х 84/16. Объем 2. п.л. Тираж 120 экз. Заказ 1001. Подписано в печать 04.12.2013 г.

Текст работы Ющенко, Валерий Павлович, диссертация по теме Математическое моделирование, численные методы и комплексы программ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ

УНИВЕРСИТЕТ»

НЕПРОЕКЦИОННАЯ РЕКОНСТРУКЦИЯ ИЗОБРАЖЕНИЯ ОБЪЕКТОВ ПРИ МОНОХРОМАТИЧЕСКОМ ЗОНДИРОВАНИИ ПРОСТРАНСТВА И

СИНТЕЗЕ АПЕРТУРЫ

Специальность 05.13.18 - математическое моделирование, численные методы

и комплексы программ

Диссертация на соискание ученой степени доктора технических наук

05201450664

На правах рукописи

Ющенко Валерий Павлович

Научный консультант

доктор технических наук Легкий В.Н.

Новосибирск - 2013

ВВЕДЕНИЕ.................................................................................................................7

1 АНАЛИЗ ИЗВЕСТНЫХ МЕТОДОВ СИНТЕЗА АПЕРТУРЫ НА МАЛОЙ ДАЛЬНОСТИ КАК ФАКТОР ПОИСКА ПУТЕЙ ИХ РАЗВИТИЯ..............30

1.1 Апертурный синтез в науке и технике..............................................................30

1.2 Историческая справка об использовании апертурного синтеза на малой дальности....................................................................................................................31

1.3 Анализ состояния вопросов и тенденции развития синтеза апертуры на малой дальности........................................................................................................33

1.4 Апертурный синтез в спектральной области (метод Д.Менсы)..........................38

1.5 Проблемы применения методов И.Радона и Д.Менсы для реконструкции изображения с помощью синтеза апертуры на малой дальности........................45

1.6 Обзор методов построения изображений в сечении, постановка задачи исследования.............................................................................................................48

1.6.1 Алгоритм восстановления, основанный на преобразовании Радона..........50

1.6.2 Алгоритм обратного проецирования.............................................................53

1.6.3 Реализация оператора обратного проецирования.........................................56

1.6.4 Реконструкция с помощью операции свертки и обратного.........................58

проецирования...........................................................................................................58

Выводы.......................................................................................................................60

2 РЕКОНСТРУКЦИЯ РАДИОЛОКАЦИОННОГО ИЗОБРАЖЕНИЯ ОБЪЕКТА ПРИ СИНТЕЗЕ АПЕРТУРЫ ДЛЯ ПРЯМОЛИНЕЙНОЙ ТРАЕКТОРИИ..........................................................................................................61

2.1 Теоретическое обоснование метода..................................................................61

2.2 Алгоритм обработки сигнала и построения изображения..............................68

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

2.4 Применение теоремы о центральном сечении к результатам..........................90

томографии методом прямолинейного синтеза апертуры.....................................90

2.5 Проблемы интерференции при когерентной томографии..............................96

2.6 Операции с матрицами изображения после синтеза апертуры для устранения интерференции...................................................................................102

2.7 Отличия предлагаемого метода от известного метода Радона, метода обратного проецирования и метода обращения по Фурье..................................109

2.8 Сравнение точечной и непрерывной моделей объекта.................................116

2.9 Проверка корреляционного метода синтеза апертуры антенны в задачах траекторного контроля для реальных сигналов...................................................127

Выводы.....................................................................................................................134

3 СИНТЕЗ АПЕРТУРЫ АНТЕННЫ ДЛЯ ЗАМКНУТОЙ ТРАЕКТОРИИ ДВИЖЕНИЯ ЛОКАТОРА...................................................................................136

3.1 Постановка задачи кругового синтеза апертуры в области..........................136

пространственных функций...................................................................................136

3.2 Сугь кругового синтеза апертуры в области пространственных функций... 137

3.2 Математическое и модельное описание кругового синтеза.........................138

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

3.4 Синтезирование апертуры для квадратной траектории................................159

3.5 Сопоставление алгоритма реконструкции на основе взаимных корреляционных функций с алгоритмами реконструкции изображения по проекциям................................................................................................................171

Выводы.....................................................................................................................180

4 МОДЕЛЬНОЕ ТОМОГРАФИРОВАНИЕ ГОМОГЕННЫХ ОБЪЕКТОВ С ВКЛЮЧЕННОЙ НЕОДНОРОДНОСТЬЮ......................................................183

4.1 Существо вопроса и постановка задачи.........................................................183

4.2 Выбор метода томографирования...................................................................183

4.3 Модель гомогенной области............................................................................184

4.4 Результаты модельного томографирования гомогенной области с включенной неоднородностью..............................................................................187

4.5 Устойчивость к шумам непроекционного метода реконструкции изображения объектов............................................................................................189

4.6 Результат томографии точечных объектов, наблюдаемых в адитивном гауссовом шуме.......................................................................................................192

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

4.8 Количественная мера оценки качества изображения в шумах...................198

4.9 Характеристика контрастного разрешения для метода................................200

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

Выводы.....................................................................................................................201

5 ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ..........................................202

5.1 Цели и задачи эксперимента, возможные пути получения экспериментальных результатов...........................................................................202

5.2 Сущность экспериментов Д.Менсы................................................................205

5.3 Выбор типа волн и среды томографирования................................................209

5.4 Эксперименты по реконструкции акустического изображения объекта в задачах траекторного контроля.............................................................................210

5.4.1 Выбор способов проведения экспериментов и построение экспериментальной установки...............................................................................210

5.4.2 Устройство экспериментальной установки.................................................216

5.4.3 Экспериментальные траекторные сигналы и их спектры, сопоставление с модельными траекторными сигналами, результаты реконструкции................223

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

5.5 Эксперименты по томографированию в диапазоне частот 8-10 ГГц электромагнитного излучения...............................................................................243

5.5.1 Синтез апертуры с учетом границы двух сред............................................247

5.5.2 Результаты экспериментов по томографированию....................................253

5.5.3 Объяснение ячеистой структуры изображения..........................................266

5.5.4 Анализ результатов эксперимента...............................................................267

Выводы.....................................................................................................................269

ЗАКЛЮЧЕНИЕ......................................................................................................271

Основные результаты диссертационной работы.................................................271

СПИСОК ЛИТЕРАТУРЫ.......................................................................................274

ПРИЛОЖЕНИЯ.......................................................................................................288

ПРИЛОЖЕНИЕ А...................................................................................................288

Список сокращений................................................................................................288

ПРИЛОЖЕНИЕ Б....................................................................................................289

Символы операций..................................................................................................289

ПРИЛОЖЕНИЕ В...................................................................................................290

Программа и результаты расчета двумерной взаимной корреляционной функции для точечного объекта............................................................................290

ПРИЛОЖЕНИЕ Г....................................................................................................291

Программа и результаты расчета передаточной функции точки.......................291

ПРИЛОЖЕНИЕ Д...................................................................................................292

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

ПРИЛОЖЕНИЕ Е....................................................................................................293

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

ПРИЛОЖЕНИЕ Ж..................................................................................................294

Сведения о практическом использовании результатов диссертации................294

ВВЕДЕНИЕ

Актуальность темы диссертации. Поиск новых методов реконструкции изображения объектов, основанных на непроекционных способах сбора данных, расширяет возможности исследования объектов. Имеет большое значение способ сбора данных и связанные с ним методы извлечения информации, результатом обработки которой является реконструкция изображения объектов. Реконструированное изображение является наивысшим достижением информационно-измерительного процесса и даёт исчерпывающую информацию об объекте или о его внутренней структуре. Поэтому поиск новых методов реконструкции изображения - это всегда актуальная задача. В 1979 году Годфри Ныо-болд Хаунсфилд и Аллан МакЛеод Кормак получили Нобелевскую премию за проекционные методы реконструкции изображения объектов. Как выяснилось позже, математические основы проекционной томографии были разработаны австрийским математиком Иоганном Радоном в 1917 году, а Годфри Ныоболд Хаунсфилд и Аллан МакЛеод Кормак переоткрыли изобретение И. Радона.

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

пациента в пределах от 0° до 360°. Далее, используя преобразование Радона, восстанавливают картину сечения пациента по проекциям.

Предпринимаются попытки вместо рентгеновского излучения использовать электромагнитные волны. В Томском университете [96, 97] реализован радиоволновой томограф с монохроматическим излучением в трёхсантиметро-вом диапазоне волн по выше описанной методики просвечивания исследуемого объекта. При этом использовалась двойная фокусировка, чтобы получить тонкий проекционный луч, так как толщина луча определяет разрешающую способность проекционного метода реконструкции изображения. К сожалению, не представлены результаты работы с биологическими объектами.

Существенный вклад в развитие радиоволновой проекционной томографии сделан в начале восьмидесятых годов американскими учёными Д. Мен-сой, Ш. Халеви и Г. Уэйдом в опубликованной работе «Применение методов когерентной доплеровской томографии для получения изображения на СВЧ». Ими предлагалось использовать монохроматический зондирующий сигнал для реконструкции изображения объектов проекционными методами с помощью кругового апертурного синтеза. Новизна их работы состояла в том, что вместо просвечивания объектов они стали работать с отражённой от объекта волной. Проекционный метод с синтезом апертуры по кругу был опробован ими в диапазоне СВЧ на волне 3 см. Получены передаточные функции одноточечных объектов (функции рассеяния одноточечных объектов), при различных секторах синтеза апертуры.

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

По материалам открытой печати зарождение апертурного синтеза на малой дальности можно отнести к 1970 годам с появлением работ по голографии и томографии. В это время вышли работы, указывавшие на тесную связь голографии и апертурного синтеза [18, 19, 50, 60, 112].

Появились также статьи, указывающие на голографическую трактовку эффекта Доплера [99, 112]. В настоящее время доплеровские измерительные системы ближней локации работают по-прежнему на упрощенных принципах синтеза апертуры, не используя многих информационных возможностей, которые может предоставить полноценный синтез апертуры антенны. Объяснить это можно, если исходить из целесообразности. Нужно ли расширять количество информации об объекте, если его все равно придется разрушать? Совершенно другой оттенок вопрос целесообразности приобретает в задачах визуализации пролётной ситуации на этапе испытаний и отработки всей системы и отдельных узлов зенитного комплекса. В этом случае информация о промахе и моменте взаимодействия может оказаться для разработчиков особенно ценной, позволяющей корректировать его в направлении, повышения эффективности.

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

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

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

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

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

/ког^\2/Д\ = с/Д/, где /ког - пространственная длина когерентности, \ - длина волны, f— частота зондирующего сигнала, Д/- абсолютная нестабильность частоты, с - скорость света.

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

на частоте 10 ГГц будет равна А/ = 10~б/ = 104Гц. Тогда длина пространственной когерентности составит. 1К0Г = с/А/ = 3-108/104 = 30км. Так как взаимодействие локатора с целыо происходит на участке траектории длиной 1 км, который значительно меньше длины пространственной когерентности, то можно считать, что условие когерентности обеспечивается. Учитывая связь когерентности и �