CC BY
108
Наука к Образование
МГТУ им. Н.Э. Баумана
Сетевое научное издание
УДК 621.396.96
Фокусировка многочастотных мультистатических радиоголограмм методом неэквидистантного БПФ
Крайний В. И.1, Семёнов А. Н.1' , ''ьетепоу.ап&Ьтмши
Чапурский В. В.1
1МГТУ им. Н.Э. Баумана, Москва, Россия
Приводится описание нового метода построения радиоизображений на основе мультистатических многочастотных радиоголограмм. Существенным преимуществом метода является возможность использования неэквидистантных раздельных приёмных и передающих элементов антенной решетки. Кратко описаны принципы получения мультистатической голограммы на примере одночастотного случая и восстановления изображения по ней методом обратных проекций. Выводится соотношение связывающее получение радиоизображения к вычислению неквидистантного преобразования Фурье. Приводятся результаты восстановления феноменологической модели 3х мерного крестообразного объекта методом обратных проекций и с помощью неэквидистантного преобразования Фурье.
Ключевые слова: радиоголография, восстановление радиоизображения, мультистатические голограммы, не эквидистантное БПФ, дискретное преобразование Фурье
Введение
К настоящему времени выполнено и опубликовано значительное число научных работ, посвященных общим принципам функционирования и теории построения радиоголографических изображений (РИ) объектов на основе нового метода мультистатических радиоголограмм (МРГ) [1-5]. Метод существенно отличается от известного ранее метода классических радиоголограмм (КРГ) [6,7].
Качественно отличия методов КРГ и МРГ сводятся к следующему. Для метода КРГ используются эквидистантные неподвижные двумерные или сканирующие линейные антенные решетки (АР), состоящие из совмещенных на передачу и прием (приемопередающих) антенных элементов. Для исключения дифракционных максимумов на РИ, восстановленном по методу КРГ, пространственный шаг приемо-передающих элементов в заполненной АР должен быть порядка половины рабочей длины волны. Поэтому для получения требуемого разрешения на сфокусированном РИ, которое зависит от площади регистрации КРГ, может потребоваться большое число приемо-передающих элементов в
Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2015. № 11. С. 292-300.
Б01: 10.7463/1115.0826094
Представлена в редакцию: 16.09.2015 Исправлена: 30.09.2015
© МГТУ им. Н.Э. Баумана
заполненной АР. Число независимых отсчетов дифрагированного поля равно при этом числу приемо-передающих элементов.
В отличие от метода КРГ для получения МРГ возможно использование разреженной АР, состоящей из раздельных и распределенных по АР передающих и приемных элементов. Пространственные отсчеты дифрагированного на объекте электромагнитного поля при методе МРГ получают для всех возможных бистатических пар «передающий элемент - приемный элемент». При этом количество независимых отсчетов поля равно произведению числа передающих элементов на число приемных элементов и может примерно равняться числу отсчетов поля в методе КРГ. Дополнительное улучшение качества РИ для трехмерных объектов в методах КРГ и МРГ достигается за счет применения многочастотного излучения [3,4]. При равном числе независимых отсчетов поля в обоих рассматриваемых методах, одинаковой разрешающей способности и качестве сфокусированных РИ (в части отсутствия дифракционных артефактов) метод МРГ имеет существенные преимущества по общему числу элементов АР. Так, при использовании сканирующей линейной АР выигрыш метода МРГ по числу элементов может быть более чем на порядок [5].
Традиционными алгоритмами фокусировки в методе КРГ являются алгоритмы двумерного быстрого преобразования Фурье (БПФ) [6,7], в то время как для фокусировки МРГ ввиду неэквидистантности по пространству бистатических отсчетов поля до настоящего времени использовался только метод обратных проекций [8]. Метод обратных проекций требует больших затрат машинного времени [9], что в еще большей мере относится и к фокусировке МРГ. По указанной причине в данном докладе в интересах поиска быстрых алгоритмов фокусировки МРГ рассмотрено преобразование характерного для МРГ алгоритма обратных проекций к форме двумерного неэквидистантного БПФ (НБПФ) [10]. Согласно имеющимся в зарубежных изданиях данным вычислительная эффективность НБПФ (Nonuniform FFT - NUFFT) [10-12] близка к эффективности классического БПФ на эквидистантной сетке отсчетов. В данном докладе для конкретных конфигурации и параметров системы регистрации МРГ приведен пример фокусировки многочастотной МРГ и получения РИ на основе алгоритмов обратных проекций и НБПФ для плоского модельного многоточечного объекта и дано сравнение их вычислительной эффективности.
1. Фокусировка МРГ по методу суммы обратных проекций
Основы теории фокусировки РИ методом суммы обратных проекций в задачах МРГ представлены в работе [3]. В данном работе рассматривается случай многочастотной МРГ (каждый передатчик излучает на своей частоте) для которого комплексный корреляционный интеграл ()(гп ), модуль которого \о{г() )| естьРИ объекта, равен:
(1)
Здесь с - скорость распространения электромагнитных волн, a Vn к есть комплексный отсчет МРГ на частоте апк при работе пары «n - ый ПдЭ, к-ый ПрЭ», а |а| означает модуль (длину) вектора a. Положение точки в плоскости наблюдаемого плоского объекта задается вектором Го = 10z + ЛГо, где r0z = ||0,0, z0||T и Дг0 = ||х0,у0, 0||T, а собственно вектор ЛГ0 - есть двумерная координата РИ.
Для совокупности M + 1 точечных отражателей с векторами координат Г коэффициентами отражения am, m = 1, 2, ...M, имеем
м
(m)
и
Kk = Yja>nQxv\ —
m=0
а , n,k ím) í m)
г — г 1 'tn + 1 ) —'rk
С _ _
(2)
Формула (1) решает задачу получения РИ по многочастотной МРГ, являясь
выражением метода суммы обратных проекций для данного случая. Метод не накладывает ограничений на пространственное расположение элементов АС, оно может быть как эквидистантным, так и неэквидистантным.
2. Алгоритм фокусировки МРГ на основе неэквидистантного ДПФ
При решении поставленной задачи создания быстрого алгоритма фокусировки рассмотрим одночастотную МРГ и преобразование формулы (1) к форме двумерного неэквидистантного дискретного преобразования Фурье (НДПФ). Для этого каждое из расстояний |г0 -и |г0 -гЛ| в показателе экспоненты в (1) представим в виде
Г0 rtn \ = Г0 z + ^Г0 rtn\ =
»0 z+Kn
(0)
Г0 Ггк \ = Г0 z + ^Г0 Ггк\ =
Г0 z+Kk
(0)
где вектора
'tn Xtn, ytn ,0
'rk
= || xrn, Угк,0|| Т , Лг0 = ||^ У0, 011
в силу чего
= Лго -»;„ = ||хо -УыМ , А7 = Лго ~ггк = Цхо-угк,
(0)_
(3)
(4)
(5)
Разложения расстояний Френеля [6]:
г +А(0)
Г0z + ^tn
и
г +А(0)
r0z + ^ rk
г +А(0)
Г0 z + ^tn
:(о)
/(2z0),
г +А(0)
0z + Аrk
запишем, используя приближение
(2z0). (6)
\к
Можно показать, что для геометрии наблюдения на рис. 1 оценка погрешности приближения (6) выводится из приближения более высокого порядка, имеющего в данном случае вид:
Ч-.+К
(0)
1
:(о)
2 zn
1-1
:(о)
4 z 2
1-Oz +Кк
(о)
1
Vk
2 zn
1-1
Vk
4 z2
(7)
Как показывают расчеты, для характерных размеров АР порядка 1.0 м и такого же размера МРГ, при удалении г0 = 2,5 м выражения в квадратных скобках (7) отличается от
Т
Т
1,0 всего на 4%. Поэтому применение приближения Френеля (6) влечет за собой столь же небольшую относительную погрешность по фазе комплексного РИ ()(гп ).
Суммируя модули расстояний (7) и учитывая (5), для бистатической суммы расстояний в показателе экспоненты в (1) получим следующее квадратичное приближение:
|г0 -Гг«| + |г0 -Ггк| =
г +А(0)
*0 г + Лгп
+
г +А(0)
Г0 г + Лгк
I * |2 I |2 | |2
о , Аго , К\ + Ы 1
'' 2г0 +J-— +
г
0
2г0 г0
--(Ггп + Ггк, Л1-0). (8)
Для большего удобства и простоты представления результирующего алгоритма обработки представим двойную сумму в (1) по индексам (п, к), изменяющимся в пределах п = 0,1,...,, к = 0,1,...,Ыг в виде однократной суммы по некоторому индексу р, который, следовательно, будет изменяться в пределах р = 0,1,...,Р, где Р = (N +1)(Ыг +1)-1. При этом прямое и обратное преобразования из (п, к) в р и из р в (п, к) имеют соответственно вид:
р = п + к (N +1), (9)
Кр)=р-
N +1
(N +1), к (р ) =
N +1
(10)
где квадратные скобки обозначают целую часть числа. После этого вместо (1) с учетом (8)-(10) можно записать
О(Аг0) = Х Ур ехр
р=0
2 г0 +
|Лг
2
Гп( р) + Ггк (р)
где Ур =¥(п[р),к[р)У Экспонента ехр^'—
2гп
2гп
|Лг
Ггк(р),Лг0)|*, (11)
в (11) выносится из под знака
суммы и с учетом взятия модуля |о(Дг0)| опущена. Введем скорректированную МРГ IVр по формуле
= V ехр <
А
с
гп( р)
гк( р)
2 гп
(12)
Тогда, учитывая, что скалярное произведение в показателе второй экспоненты в (11) имеет вид
(Ггп(р) + Ггк( р), Лг0 ) = ^рх0 +^p>'0, (13)
где
Vр = хгп( р) + хгк (р) и ур = Угп( р) + Угк( р), (14)
окончательно для плоского РИ имеем
|(?(хо>Л))|:
Е ехр | ^ {/-1рхо+уРУо)
(15)
Из (14) и (15) следует, что Q(xñ,yñ) имеет вид двумерного НДПФ [11-12] с неэквидистантным по p распределением двумерных компонент или узлов (,ир,vp). В
свою очередь выходные компоненты координат РИ (х0, у0) дискретизируют с равномерным шагом, исходя из размеров и числа точек РИ по аналогии с тем, как это делается в [12].
3. Пример фокусировки радиоизображения по МРГ с помощью НБПФ
В [12] для быстрого вычисления НДПФ привлекается специальный алгоритм НБПФ, основанный на так называемой быстрой гауссовской привязке неэквидистантных узлов (fast Gaussian gridding). Для этого алогоритма, реализованного программно в среде Matlab [13] были подготовлены исходные данные для вычисления НДПФ в формуле (15) и получения радиоизображения |Ó(x0, j0)|.
Численный машинный эксперимент проводился для следующих исходных данных. Частота излучения f = а>0/(2л) = 10 ГГц ( \ = 3 см). Полоса многочастотного сигнала Af =
6 ГГц и число независимых частот 81. Разреженная АС типа MIMO имеет передающую АР размером 2*2 м с числом элементов 9*9 и межэлементным шагом 25 см, а также расположенную осесимметрично в той же плоскости приемную АР размером 2,3*2,3 м с числом элементов 9*9 и межэлементным шагом 28,8 см. Многоточечный плоский объект расположен параллельно АС на удалении 3 м по нормали к ее центру и имеет форму симметричного прямоугольного креста с одинаковыми размерами 30 см по горизонтали и вертикали, образованного из 12 точечных отражающих элементов с одинаковыми расстояниями между ними по 15 см.
Топографическая диаграмма РИ крестообразного объекта, сфокусированного по методу обратных проекций в соответствие с формулой (1), представлена на рис. 2а,б,в. Аналогичная диаграмма РИ, сфокусированного методом двумерного НБПФ в соответствие с формулой (15), показана на рис. 3а,б,в.
б
Рис. 2 - Радиоизображения крестообразного объекта, полученные методом обратных проекций в плоскости
XY (а), в плоскости XZ (б) и в плоскости YZ (в).
а
в
20 40 60 80 100 120 140 160 160 200 20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200
а б в
Рис. 3 - Радиоизображения крестообразного объекта, полученные методом НБПФ в плоскости ХУ (а), в
плоскости Х2 (б) и в плоскости У2 (в).
Сравнение РИ на рис. 2-3 свидетельствует о несколько лучшем качестве РИ, полученного по методу обратных проекций по точной формуле (1). Было также установлено, что незначительное ухудшение качества РИ при методе НБПФ обусловлено использованием приближения Френеля (15), а не собственно алгоритмом НБПФ. В тоже время оценка времен вычисления показала, что время вычисления РИ на ПЭВМ с процессором 17 по методу обратных проекций составило 4,42с, а по методу НБПФ - всего 0,033 с, что соответствует выигрышу по скорости фокусировки РИ при использовании НБПФ около 150 раз.
Статья выпущена в рамках НИОКТР "Реализация комплексного проекта по созданию высокотехнологичного производства радиолокационного комплекса для системы управления воздушным движением с удаленной диспетчеризацией", выполняемой МГТУ им. Н.Э. Баумана совместно с ОАО "РТИ" в рамках комплексного проекта по созданию высокотехнологичного производства, в целях реализации постановления Правительства Российской Федерации от 9 апреля 2010 г. № 218 «О мерах государственной поддержки развития кооперации российских высших учебных заведений, государственных научных учреждений и организаций, реализующих комплексные проекты по созданию высокотехнологичного производства», при финансовой поддержке по проекту Министерства образования и науки Российской Федерации.
Список литературы
1. Крайний В.И., Семёнов А.Н., Чапурский В.В. Фокусировка одночастотных мультистатических радиоголограмм методом неэквидистантного быстрого преобразования Фурье // VIII Всероссийская научно-техническая конференция «Радиолокация и радиосвязь» (Россия, Москва, 24-26 ноября 2014 г.): тр. Москва, 2014. С. 77-81.
2. Чапурский В.В. Мультистатическая радиоголография // 4-ая Международная конференция «Акустооптические и радиолокационные методы измерений и обработки информации» ARMIMP-2011 (Россия, Суздаль, 20-22 сентябрь 2011 г.): тр. Суздаль, 2011. С. 99-102.
3. Чапурский В.В. Получение радиоголографических изображений объектов на основе разреженных антенных решеток типа MIMO с одночастотным и многочастотным излучением // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2011. № 4. С. 72-91.
4. Чапурский В.В. Синтезирование радиоизображений объектов с помощью линейной антенной решетки типа MIMO // Инженерный журнал: наука и инновации. 2012. № 8. DOI: 10.18698/2308-6033-2012-8-323
5. Разевиг В.В., Бугаев А.С., Чапурский В.В. Сравнительный анализ фокусировки классических и мультистатических радиоголограмм // Радиотехника. 2013. № 8. С. 817.
6. Зверев В.А. Радиооптика. М.: Советское радио, 1975. 304 с.
7. Sheen D.M., McMakin D.L., Hall T.E. Three-Dimentional Millimeter-Wave Imaging for Conceald Weapon Detection // IEEE Transactions on Microwave Theory and Techniques.
2001.Vol. 49, no. 9. P. 1581-1592. DOI: 10.1109/22.942570
8. Курикша А.А. Алгоритм обратной проекции в задачах восстановления пространственного распределения источников волн // Радиотехника и электроника.
2002. Т. 47, № 12. С. 1484-1489.
9. Курикша А.А. Быстрые алгоритмы для оценки пространственного распределения источников сигналов // 3-я Международная конференция «Цифровая обработка сигналов и её применение» (DSPA-2000): докл. Т. 1. Режим доступа: http://www.autex.spb.ru/dspa/dspa2000/part1.htm (дата обращения 01.10.2015) .
10. Dutt A., Rokhlin V. Fast Fourier transform for nonequispaced data // SIAM Journal on Scientific Computing. 1993. Vol. 14, iss. 6. P. 1368-1393. DOI: 10.1137/0914081
11. He Xuezhi, Xu Hao, Liu Changchang, Wang Dongjin, Chen Weidong. A New Approach to Distributed Passive Radar Imaging by 2-D NUFFT // 2010 IEEE 10th International Conference on Signal Processing (ICSP). IEEE Publ., 2010. P. 2067-2070. DOI: 10.1109/IC0SP.2010.5655727
12. Leslie Greengard, June-Yub Lee. Accelerating the Nonuniform Fast Fourier Transform // SIAM Review. 2004. Vol. 46, no. 3. P. 443-453. DOI: 10.1137/S003614450343200X
13. Ferrara M. NUFFT, NFFT, USFFT. 2009 // Matlab Central: website. Режим доступа: http://www.mathworks.com/matlabcentral/fileexchange/25135-nufft--nfft--usfft (дата обращения 30.04.2015).
Science and Education of the Bauman MSTU, 2015, no. 11, pp. 292-300.
DOI: 10.7463/1115.0826094
Received: Revised:
16.09.2015 30.09.2015
Science^Education
of the Bauman MSTU
ISS N 1994-0408 © Bauman Moscow State Technical Unversity
Multi-Static Multi-Frequency Image Reconstruction Based on Non-Uniform FFT
V.I. Krainy1, A.N. Semenov1*, V.V. Chapurskiy1
&ern.enoy.an@bm5tu-nj
:Bauman Moscow State Technical University, Moscow, Russia
Keywords: radio hologram, image restoring, multi-static holography, non-uniform FFT
By now a considerable number of articles on the common principles and theory of radio holographic imaging of objects based on a new method of multi-static radio-holography (MRH) have been are published.
In contrast to the method of image reconstruction using the classical radio holography (CRH) to have the MRH approach it is possible to use a non-equispaced antennas array with independent transmitting and receiving elements. Radio-holography samples are received for each pair of "transmitting element-receiving element". The number of independent samples is equal to the product of the number of transmitter elements and receiver elements and may be approximately equal to the number of samples used in CRH. An additional improvement in the quality of reconstructed image for three-dimensional objects in both approaches is achieved through the use of multi-frequency radiation.
The wide common method of reconstructing images from the classical radio-hologram uses algorithms of two-fold Fast Fourier Transform (FFT), while due to focusing the multi-static radio-hologram a back projection method has been used. The back projection method is time-consuming. However, it can be transformed to the form of two-fold non-uniform FFT (NuFFT). According to data available in the foreign references a computational efficiency of NuFFT is close to the efficiency of the classical FFT on an equidistant grid sampling. In the specific antenna configuration an example of the multi-static multi-frequency focusing and image reconstruction of three dimensional model object with back projection and NuFFT approach are presented.
References
1. Krainii V.I., Semenov A.N., Chapurskii V.V. Focusing single-frequency multistatic radioholograms by nonequidistant fast Fourier transform. 8 Vserossiiskaya nauchno-tekhnicheskaya konferentsiya "Radiolokatsiya i radiosvyaz"': tr. [Proc. of the 8th All-Russian scientific-technical conference "Radiolocation and radio communication"], Russia, Moscow, 24-26 November, 2014, pp. 77-81. (in Russian).
2. Chapurskii V.V. Multistatic radio- holography. 4-aya Mezhdunarodnaya konferentsiya "Akustoopticheskie i radiolokatsionnye metody izmerenii i obrabotki informatsii" ARMIMP-2011: tr. [Proc. of the 4th International Conference "Acoustic-optic and radar measurement methods and information processing" ARMIMP-2011], Russia, Suzdal, 20-22 September, 2011, pp. 99-102. (in Russian).
3. Chapurskii V.V. Receiving the Radio-Holographic Images of Objects on the Basis of Disperse MIMO-Type Antenna Arrays with the Single-Frequency and Multi-Frequency Radiations. Vestnik MGTU im. N.E. Baumana. Ser. Priborostroenie = Herald of the Bauman Moscow State Technical University. Ser. Instrument Engineering, 2011, no. 4, pp. 72-91. (in Russian).
4. Chapurskij V.V. Object radio image synthesizing by means of linear antenna array of the MIMO series. Inzhenernyy zhurnal: nauka i innovatsii = Engineering Journal: Science and Innovation, 2012, no. 8. DOI: 110.18698/2308-6033-2012-8-323 (in Russian).
5. Razevig V.V., Bugaev A.S., Chapursky V.V. The comparative analysis of classical and multistatic microwave holograms focusing. Radiotekhnika = Radioengineering, 2013, no. 8, pp. 8-17. (in Russian).
6. Zverev V.A. Radiooptika [Radiooptics]. Moscow, Sovetskoe radio Publ., 1975. 304 p. (in Russian).
7. Sheen D.M., McMakin D.L., Hall T.E. Three-Dimentional Millimeter-Wave Imaging for Conceald Weapon Detection. IEEE Transactions on Microwave Theory and Techniques,
2001, vol. 49, no. 9, pp. 1581-1592. DOI: 10.1109/22.942570
8. Kuriksha A.A. The Reciprocal Projection Algorithm Used for Reconstructing a Spatial Distribution of Wave Sources. Radiotekhnika i elektronika, 2002, vol. 47, no. 12, pp. 14841489. (English version of journal: Journal of Communications Technology and Electronics,
2002, vol. 47, no. 12, pp. 1361-1363.).
9. Kuriksha A.A. Fast algorithms for estimation of signal's south's distribution. 3-ya
Mezhdunarodnaya konferentsiya "Tsifrovaya obrabotka signalov i ee primenenie" (DSPArd
2000): dokl. T. 1 [Proc. of the 3 International Conference "Digital Signal Processing and its Applications" (DSPA-2000). Vol. 1]. Available at:
http://www.autex.spb.ru/dspa/dspa2000/part1.htm , accessed 01.10.2015. (in Russian).
10. Dutt A., Rokhlin V. Fast Fourier transform for nonequispaced data. SIAM Journal on Scientific Computing, 1993, vol. 14, iss. 6, pp. 1368-1393. DOI: 10.1137/0914081
11. He Xuezhi, Xu Hao, Liu Changchang, Wang Dongjin, Chen Weidong. A New Approach to Distributed Passive Radar Imaging by 2-D NUFFT. 2010 IEEE 10th International Conference on Signal Processing (ICSP). IEEE Publ., 2010, pp. 2067-2070. DOI: 10.1109/ICQSP.2010.5655727
12. Leslie Greengard, June-Yub Lee. Accelerating the Nonuniform Fast Fourier Transform. SIAM Review, 2004, vol. 46, no. 3, pp. 443-453. DOI: 10.1137/S003614450343200X
13. Ferrara M. NUFFT, NFFT, USFFT. 2009. Matlab Central: website. Available at: http://www.mathworks.com/matlabcentral/fileexchange/25135-nufft--nfft--usfft , accessed 30.04.2015.
