CN109597071B - 滑动式弧形阵列微波成像方法 - Google Patents

滑动式弧形阵列微波成像方法 Download PDF

Info

Publication number
CN109597071B
CN109597071B CN201811638042.2A CN201811638042A CN109597071B CN 109597071 B CN109597071 B CN 109597071B CN 201811638042 A CN201811638042 A CN 201811638042A CN 109597071 B CN109597071 B CN 109597071B
Authority
CN
China
Prior art keywords
arc
distance
array
microwave
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811638042.2A
Other languages
English (en)
Other versions
CN109597071A (zh
Inventor
黄平平
杜欣
谭维贤
徐伟
李亚超
张振华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inner Mongolia University of Technology
Original Assignee
Inner Mongolia University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Inner Mongolia University of Technology filed Critical Inner Mongolia University of Technology
Priority to CN201811638042.2A priority Critical patent/CN109597071B/zh
Publication of CN109597071A publication Critical patent/CN109597071A/zh
Application granted granted Critical
Publication of CN109597071B publication Critical patent/CN109597071B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种滑动式弧形阵列微波成像方法,包括以下步骤:S1在空中平台搭载微波成像系统和惯性导航测量系统进行飞行时,通过惯性导航测量系统获得空中平台在空中的当前位置信息;S2根据空中平台的飞行姿态和微波成像系统的参数,设置空中平台的滑动间距;S3根据所述当前位置信息和滑动间距,通过微波开关网络控制弧形接收阵列天线的接收天线阵元与其的发射天线阵元之间的切换,实现对地面场景进行瞬时环状扫描;S4通过微波开关网络控制接收天线阵元接收并存储通过瞬时环状扫描获得的地面场景的散射回波;S5对散射回波进行成像处理以获得空中平台的周围的地面场景的微波图像。

Description

滑动式弧形阵列微波成像方法
技术领域
本发明涉及微波成像技术对地观测技术领域,特别涉及滑动式弧形阵列微波成像方法。
背景技术
专利公开出版物CN104267400A提供了一种用于MIMO-SAR成像的微波信号收发系统、方法及成像系统,其中提出了机载弧形阵列微波成像方法。该方法主要是通过弧形阵列天线收发电磁波,可以实现对平台周围场景的360°范围观测,这相比较于常规线性阵列SAR只能进行单一的前视、侧视、后视成像,大大增加了对场景的观测视角。
但是,现有弧形阵列微波成像方法主要以实现对固定中心的环状区域场景进行观测,未提及当平台进行匀速直线运动时的针对弧形阵列微波的一种连续成像模式及方法。目前,国内外均尚未给出有关于滑动式弧形阵列微波成像模式及方法。因此,尚不能解决当平台匀速直线运动时的360°大视角高分辨率连续观测成像。
发明内容
为了解决现有技术中存在的上述问题和缺陷的至少一个方面,本发明提供了一种滑动式弧形阵列微波成像方法,其能够在空中平台为匀速直线运动时进行360°大视角高分辨率连续观测成像。
所述技术方案如下:
本发明的一个目的是提供了一种滑动式弧形阵列微波成像方法:包括
获取搭载微波成像系统的空中平台飞行的位置信息;
基于空中平台的飞行参数和微波成像系统的参数设置空中平台的滑动间距;
至少根据该滑动间距,对目标场景进行扫描;
收集目标场景响应于所述扫描的散射回波;
基于散射回波,进行成像。
根据本发明的一个方面,具体包括以下步骤:
S1在空中平台搭载微波成像系统和惯性导航测量系统进行飞行时,通过所述惯性导航测量系统获得所述空中平台在空中的当前位置信息;
S2根据所述空中平台的飞行姿态和所述微波成像系统的参数,在所述微波成像系统中设置所述空中平台的滑动间距;
S3根据所述当前位置信息和所述滑动间距,通过所述微波成像系统中的微波开关网络控制所述微波成像系统中的弧形接收阵列天线的接收天线阵元与弧形发射阵列天线的发射天线阵元之间的切换,实现对地面场景进行瞬时环状扫描;
S4通过所述微波开关网络控制所述接收天线阵元接收并存储通过所述瞬时环状扫描获得的所述地面场景的散射回波;
S5对所述散射回波进行成像处理以获得所述空中平台的周围的所述地面场景的微波图像。
进一步地,在步骤S2中,所述飞行姿态包括所述空中平台的飞行高度,所述微波成像系统的参数包括所述弧形阵列的半径、所述弧形接收阵列天线的入射角或所述弧形发射阵列天线的入射角和俯仰向-3dB波束宽度。
具体地,所述成像观测最远地距的表达式为:
Figure BDA0001930517460000021
所述成像观测最近地距的表达式为:
Figure BDA0001930517460000022
其中,rmax为成像观测最远地距,rmin为成像观测最近地距,φin为弧形接收阵列天线或弧形发射阵列天线的入射角,φ-3dB为弧形接收阵列天线或弧形发射阵列天线的俯仰向-3dB波束宽度,H为所述飞行平台的高度,Rarc为所述弧形阵列的半径。
优选地,所述弧形阵列的中心滑动间距的表达式为:
ΔL≤rmax-rmin
其中,ΔL为所述弧形阵列的中心滑动间距,rmax为所述成像观测最远地距,rmin为成像观测最近地距。
具体地,在步骤S5中,所述散射回波的成像处理包括以下步骤:
S51将所述散射回波通过所述微波成像系统中的混频器进行混频处理,获得混频回波信号;
S52将所述混频回波信号沿所述散射回波的距离向进行傅里叶变换,以获得距离频域的回波信号;
S53将所述距离频域的回波信号通过匹配滤波获得滤波后的距离频域信号;
S54对所述滤波后的距离频域信号进行相干叠加,获得所述地面场景的微波图像。
具体地,在步骤S54中,对所述滤波后的距离频域信号进行相干叠加包括以下步骤:
S541通过所述散射回波建立所述地面场景对应的图像空间;
S542将所述图像空间分别沿所述弧形阵列的弧形阵列向和所述距离向进行二维离散化,获得二维图像空间和所述二维图像空间的像素个数;
S543计算所述二维图像空间中的每一个像素位置对应的滤波函数,并通过所述滤波函数与所述滤波后的距离频域信号相干叠加,获得所述二维图像空间中的每一个像素的像素值;
S544根据所述二维图像空间中的每一个像素的像素值获得极坐标格式下的距离-角度图像数据;
S545将所述极坐标格式下的距离-角度图像数据通过插值转换为笛卡尔直角坐标格式下的图像数据。
进一步地,在步骤S3中,所述瞬时环状扫描为所述弧形发射阵列天线辐射电磁波信号对所述地面场景进行360°瞬时环状扫描。
具体地,所述微波成像系统采用线性调频脉冲或调频连续波工作方式。
具体地,在步骤S1中,所述空中平台为匀速直线运动,
所述滑动式阵列微波成像方法还包括以下步骤:
S6将所述微波图像显示在所述微波成像系统中的显示处理模块上。
根据本发明提供的滑动式弧形阵列微波成像方法具有以下优点中的至少一个:
(1)本发明提供的滑动式弧形阵列微波成像方法可以在空中平台匀速直线运动时对平台周围大区域场景进行连续观测成像,而且该方法不仅能够实现弧形阵列天线随平台匀速直线飞行所进行的“滑动式”观测,而且还能保留常规弧形阵列雷达全方位成像的优势,可以全天时全天候,不受迷尘、云、雨、雾的影响进行成像观测;
(2)本发明提供的滑动式弧形阵列微波成像方法填补了关于当平台进行直线运动时呈滑动式弧形阵列微波成像模式的空缺问题,实现了对飞行平台周围大视场区域高分辨率成像数据获取与连续成像处理的能力,为平台的侦察、搜救提供了强有力的支持;
(3)本发明提供的滑动式弧形阵列微波成像方法能够在空中平台为匀速直线运动时进行360°大视角高分辨率连续观测成像。
附图说明
本发明的这些和/或其他方面和优点从下面结合附图对优选实施例的描述中将变得明显和容易理解,其中:
图1是固定中心的单区域弧形阵列微波成像的几何示意图;
图2是根据本发明的一个实施例的滑动式弧形阵列微波成像的几何示意图;
图3是根据本发明的一个实施例的滑动式弧形阵列微波成像方法的流程图;
图4是现有技术中的弧形阵列天线的结构图;
图5是根据图3所示的滑动式弧形阵列微波成像方法形成的微波成像观测区域的示意图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。在说明书中,相同或相似的附图标号指示相同或相似的部件。下述参照附图对本发明实施方式的说明旨在对本发明的总体发明构思进行解释,而不应当理解为对本发明的一种限制。
本发明提供了一种新的滑动式弧形阵列微波成像方法,同固定中心的单区域环扫弧形阵列微波成像(如图1所示)相比,滑动式弧形阵列微波成像通过平台运动实现一种类似于弧形阵列“滑动式”前移的成像模式(如图2所示),该成像模式具体为使空中平台(例如直升飞机)搭载微波成像系统,在一系列特定的位置处对周围场景进行瞬时环状扫描,这种通过平台运动实现一种类似于弧形阵列“滑动式”前移的成像模式,能够大幅度扩大成像区域,同时还保留其常规成像的优势。
参见图3,其示出了根据本发明的一个实施例的滑动式弧形阵列微波成像方法,获取搭载微波成像系统的空中平台飞行的位置信息;基于空中平台的飞行参数和微波成像系统的参数设置空中平台的滑动间距;至少根据该滑动间距,对目标场景进行扫描;收集目标场景响应于所述扫描的散射回波;基于散射回波,进行成像。
所述滑动式弧形阵列微波成像方法,具体可以包括以下步骤:
S1在空中平台(例如直升飞机)搭载微波成像系统和惯性导航测量系统直线匀速飞行,并且将微波成像系统中的弧形阵列天线固定在空中平台的中部(例如腹部)位置时,通过所述惯性导航测量系统获得所述空中平台在空中的当前位置信息。
本发明中所提到的微波成像系统的更多或更详细的内容可以参见专利公开出版物CN104267400A,但是不限于此,本领域技术人员还可以根据需要选择任何可行的或兼容的成像系统。在专利公开出版物中的弧形阵列天线具体如图4所示。在圆弧方向上排列的发射天线阵元和在圆弧方向上排列的接收天线阵元共同组合成弧形阵列天线,每个发射天线阵元或接收天线阵元的通道都由快速微波开关网络控制,可以使弧形阵列天线单元按设定好的工作方式进行分时扫描。
结合滑动式弧形阵列微波成像几何对其进行说明。如图2所示,假设空中平台运动时弧形阵列接收阵列天线或弧形发射阵列天线的位置坐标为(x,y,z),Rarc为弧形阵列的半径,阵列接收阵列天线或弧形发射阵列天线的入射角及俯仰向-3dB波束宽度分别为φin和φ-3dB,目标Pn对应的坐标位置为Pn(xn,yn,zn),空中平台通过装载的弧形阵列天线进行成像,每次成像对应不同的观测场景区域,rmax为观测场景区域的外环半径,rmin为观测场景区域的内环半径,通过多个不同位置的弧形阵列微波成像实现多个观测区域的连续覆盖,从而扩大成像观测范围,这多个不同位置的弧形阵列微波成像飞行过程即可看作为“滑动式”运动。
S2根据所述空中平台的飞行姿态和所述微波成像系统的参数,在所述微波成像系统中设置所述空中平台的滑动间距,具体步骤如下:
S21根据所述空中平台的飞行姿态(例如空中平台的飞行高度)、所述弧形阵列的半径、所述弧形接收阵列天线的入射角或所述弧形发射阵列天线的入射角和俯仰向一3dB波束宽度计算成像观测最远地距(即雷达波束足印的外径)和成像观测最近地距(即雷达波束足印的内径);
所述成像观测最远地距的表达式为:
Figure BDA0001930517460000061
所述成像观测最近地距的表达式为:
Figure BDA0001930517460000062
其中,rmax为成像观测最远地距,rmin为成像观测最近地距,φin为弧形接收阵列天线或弧形发射阵列天线的入射角,φ-3dB为弧形接收阵列天线或弧形发射阵列天线的俯仰向-3dB波束宽度,H为所述飞行平台的高度,Rarc为所述弧形阵列的半径。
S22根据得到的成像观测最远地距和成像观测最近地距,设置相邻每两次进行瞬时环状扫描时形成的弧形阵列的中心滑动间距,这也是空中平台的滑动间距(参见图5)。
所述弧形阵列的中心滑动间距的表达式为:
ΔL≤rmax-rmin
其中,ΔL为所述弧形阵列的中心滑动间距,rmax为所述成像观测最远地距,rmin为成像观测最近地距。
经研究发现,当ΔL≤rmax-rmin时,则表示在空中平台向前直线飞行时,对地面场景区域能够完整地观测;而当ΔL>rmax-rmin时,则表示在空中平台向前直线飞行时,在场景的内圈区域会有观察盲区产生;因此,在空中平台飞行时需要保证弧形阵列采样时中心滑动间距ΔL≤rmax-rmin
S3根据所述位置信息和所述滑动间距,在特定的位置(即空中平台根据滑动间距到达的预设定的位置)处,提供所述微波成像系统中的微波开关网络控制所述微波成像系统中的弧形接收阵列天线的接收天线阵元与弧形发射阵列天线的发射天线阵元之间的切换,通过所述弧形发射阵列天线辐射电磁波信号对所述地面场景进行360°瞬时环状扫描。这种在特定的位置对周围场景进行瞬时环状扫描,通过空中平台运动实现一种类似于弧形阵列“滑动式”前移的成像模式,在该弧形阵列滑动式前移的成像模式中,波束足印在地面呈现为一系列的环形,称这种观测方式为滑动式弧形阵列微波成像。
弧形发射阵列天线发射的电磁波信号可以表示为:
Figure BDA0001930517460000071
其中,Str(t)为电磁波信号,Kr为调频率,fc为雷达工作频率。
S4通过所述微波开关网络控制所述接收天线阵元接收并存储通过所述瞬时环状扫描获得的所述地面场景的散射回波。
S5对所述散射回波进行成像处理获得所述空中平台的周围的所述地面场景的微波图像。具体为:
S51将所述散射回波通过所述微波成像系统中的混频器进行混频处理,获得混频回波信号。
当微波成像系统采用调频连续波(FMCW)工作方式时,弧形阵列天线辐射的微波信号经过观测区域反射,由弧形接收阵列天线进行接收,并经过微波成像系统的混频器,其接收信号可以表达为:
Figure BDA0001930517460000072
当微波成像系统采用线性调频脉冲工作方式时,弧形阵列天线辐射的微波信号经过观测区域反射,由弧形接收阵列天线进行接收,并经过微波成像系统的混频器,其接收信号可以表达为:
Figure BDA0001930517460000073
其中,Sre(t,x,y,z)为弧形接收阵列天线所接收的信号,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,Rn表示位于坐标为(x,y,z)处的弧形阵列等效采样点与Pn(xn,yn,zn)之间的距离,坐标(x,y,z)为空中平台在匀速直线运动时的弧形阵列天线等效采样点的坐标位置,C表示电磁波传播速度,具体地
Figure BDA0001930517460000074
S52若微波成像系统采用线性调频脉冲工作方式,则将所述混频回波信号沿所述散射回波的距离向进行傅里叶变换,并变换到距离频域,以获得距离频域的回波信号SFT_re(Rn,x,y,z);
距离频域的回波信号SFT_re(Rn,x,y,z)的表达式为:
Figure BDA0001930517460000081
其中,SFT_re(Rn,x,y,z)表示距离频域的回波信号,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,FTt表示沿距离向时间变量t进行傅里叶变换,Rn表示位于坐标(x,y,z)处的弧形阵列等效采样点与目标Pn(xn,yn,zn)之间的距离。
S53将经过傅里叶变换后的所述距离频域的回波信号通过匹配滤波获得滤波后的距离频域信号,其所使用的滤波函数H(fr)为:
H(fr)=exp{jπKrt2}
其中,Kr为信号的调频率。
距离频域的回波信号经滤波匹配后得到的滤波后的距离频域信号SFT_rc(Rn,x,y,z)表达式为:
Figure BDA0001930517460000082
其中,SFT_re(Rn,x,y,z)表示距离频域的回波信号,H(fr)为滤波函数,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,fc为雷达工作中心频率,Rn表示位于坐标(x,y,z)处的弧形阵列等效采样点与目标Pn(xn,yn,zn)之间的距离。
S54对所述滤波后的距离频域信号进行相干叠加,获得所述地面场景的微波图像,具体包括以下步骤:
S541通过所述散射回波建立所述地面场景对应的图像空间,即将所观测的地面场景划分为多个格,使之呈网格状,之后将所获得的散射回波信号投影至该对应的网格中,形成了地面场景对应的图像空间。
S542将所述图像空间分别以Δθ、Δr的像素大小沿所述弧形阵列的弧形阵列向和所述距离向进行二维离散化,获得二维图像空间I(nθΔθ,nrΔr),其中nθ=1,…,Nθ,nr=1,…,Nr,Nθ和Nr分别为沿弧形阵列向和距离向离散化后得到的所述二维图像空间的像素个数,Nθ和Nr的表达式为:
Figure BDA0001930517460000091
其中,Nθ为沿弧形阵列向离散化后得到的所述二维图像空间的像素个数,Nr为沿距离向离散化后得到的所述二维图像空间的像素个数,Δθ、Δr为像素大小,rmax和rmin分别为成像观测最远地距和成像观测最近地距。
S543计算所述二维图像空间中的每一个像素位置对应的滤波函数,并通过所述滤波函数与所述滤波后的距离频域信号相干叠加,获得所述二维图像空间中的每一个像素的像素值,具体为:
S5431令nθ=1,nr=1,其中,nθ和nr表示像素计数序号。nθ=1,…,Nθ,nr=1,…,Nr,Nθ和Nr分别为沿弧形阵列向和距离向离散化后得到的所述二维图像空间的像素个数;
S5432计算图像I(nθΔθ,nrΔr)中第(nθ,nr)像素对应的坐标位置(nθΔθ,nrΔr)到弧形阵列天线的等效采样点P(x,y,z)的距离Rm,并以该距离生成匹配的滤波函数HM(θ,r,z;nθΔθ,nrΔr):
Figure BDA0001930517460000092
其中,
Figure BDA0001930517460000093
z0表示地表平面高度。
S5433基于所述匹配滤波函数,求解图像的第(nθ,nr)个像素对应的数值,具体地:
Figure BDA0001930517460000101
其中,i(nθΔθ,nrΔr)为图像,HM(θ,r,z;nθΔθ,nrΔr)为滤波函数,SFT_rc(Rn,x,y,z)为上述滤波后的距离频域信号。
S5434令nr加1,若nr≤Nr,返回至步骤S5432,若nr>Nr,继续执行步骤S5435;
S5435令nθ加1,若nθ≤Nθ,令nr=1,并返回至步骤S5432,若nθ>Nθ,继续执行步骤S544。
S544基于二维图像空间中的每一个像素的像素值,输出滑动式弧形阵列微波成像图像I(nθΔθ,nrΔr),图像I(nθΔθ,nrΔr)为极坐标格式下的距离-角度图像数据;
S545将所述极坐标格式下的距离-角度图像数据转换为笛卡尔直角坐标格式下的x-y图像数据,这里使用二维的sinc插值来实现极坐标的转换的;当极坐标为(ρ,θ)时,对应的直角坐标为(x,y)时,其转换关系可用下式表示:
Figure BDA0001930517460000102
S6将所述微波图像显示在所述微波成像系统中的显示处理模块上,便于观察。
根据本发明的提供的滑动式弧形阵列微波成像方法具有以下优点中的至少一个:
(1)本发明提供的滑动式弧形阵列微波成像方法可以在平台直线运动时对平台周围大区域场景进行连续观测成像,而且该方法不仅能够实现弧形阵列天线随平台匀速直线飞行所进行的“滑动式”观测,而且还能保留常规弧形阵列雷达全方位成像的优势,可以全天时全天候,不受迷尘、云、雨、雾的影响进行成像观测;
(2)本发明提供的滑动式弧形阵列微波成像方法填补了关于当平台进行直线运动时呈滑动式弧形阵列微波成像模式的空缺问题,实现了对飞行平台周围大视场区域高分辨率成像数据获取与连续成像处理的能力,为平台的侦察、搜救提供了强有力的支撑;
(3)本发明提供的滑动式弧形阵列微波成像方法能够在空中平台为匀速直线运动时进行360°大视角高分辨率连续观测成像。
虽然本总体发明构思的一些实施例已被显示和说明,本领域普通技术人员将理解,在不背离本总体发明构思的原则和精神的情况下,可对这些实施例做出改变,本发明的范围以权利要求和它们的等同物限定。

Claims (6)

1.一种滑动式弧形阵列微波成像方法,包括:
获取搭载微波成像系统的空中平台飞行的位置信息;
基于所述空中平台的飞行参数和所述微波成像系统的参数设置所述空中平台的滑动间距;
至少根据所述滑动间距,对目标场景进行扫描;
收集所述目标场景响应于所述扫描的散射回波;
基于所述散射回波,进行成像;
所述滑动式弧形阵列微波成像方法具体包括以下步骤:
S1在所述空中平台搭载微波成像系统和惯性导航测量系统进行飞行时,通过所述惯性导航测量系统获得所述空中平台在空中的当前位置信息;
S2根据所述空中平台的飞行姿态和所述微波成像系统的参数,在所述微波成像系统中设置所述空中平台的滑动间距;其中,
在步骤S2中,设置所述空中平台的滑动间距包括以下步骤:
S21根据所述空中平台的飞行姿态和所述微波成像系统的参数计算成像观测最远地距和成像观测最近地距;其中,
所述成像观测最远地距的表达式为:
Figure QLYQS_1
所述成像观测最近地距的表达式为:
Figure QLYQS_2
其中,rmax为成像观测最远地距,rmin为成像观测最近地距,φin为弧形接收阵列天线或弧形发射阵列天线的入射角,φ-3dB为弧形接收阵列天线或弧形发射阵列天线的俯仰向-3dB波束宽度,H为飞行平台的高度,Rarc为所述弧形阵列的半径;
S22根据所述成像观测最远地距和所述成像观测最近地距,设置相邻的每两次在瞬时环状扫描时形成的弧形阵列的中心滑动间距,获得所述空中平台的滑动间距;
S3根据所述当前位置信息和所述滑动间距,通过所述微波成像系统中的微波开关网络控制所述微波成像系统中的弧形接收阵列天线的接收天线阵元与弧形发射阵列天线的发射天线阵元之间的切换,实现对地面场景进行瞬时环状扫描;
S4通过所述微波开关网络控制所述接收天线阵元接收并存储通过所述瞬时环状扫描获得的所述地面场景的散射回波;
S5对所述散射回波进行成像处理以获得所述空中平台的周围的所述地面场景的微波图像;其中,所述对所述散射回波进行成像处理包括以下步骤:
S51将所述散射回波通过所述微波成像系统中的混频器进行混频处理,获得混频回波信号;其中,
当微波成像系统采用调频连续波工作方式时,弧形阵列天线辐射的微波信号经过观测区域反射,由弧形接收阵列天线进行接收,并经过微波成像系统的混频器,其接收信号可以表达为:
Figure QLYQS_3
当微波成像系统采用线性调频脉冲工作方式时,弧形阵列天线辐射的微波信号经过观测区域反射,由弧形接收阵列天线进行接收,并经过微波成像系统的混频器,其接收信号可以表达为:
Figure QLYQS_4
其中,Sre(t,x,y,z)为弧形接收阵列天线所接收的信号,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,Rn表示位于坐标为(x,y,z)处的弧形阵列等效采样点与Pn(xn,yn,zn)之间的距离,坐标(x,y,z)为空中平台在匀速直线运动时的弧形阵列天线等效采样点的坐标位置,C表示电磁波传播速度,具体公式为:
Figure QLYQS_5
S52将所述混频回波信号沿所述散射回波的距离向进行傅里叶变换,以获得距离频域的回波信号SFT_re(Rn,x,y,z);其中,距离频域的回波信号SFT_re(Rn,x,y,z)的表达式为:
Figure QLYQS_6
其中,SFT_re(Rn,x,y,z)表示距离频域的回波信号,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,FTt表示沿距离向时间变量t进行傅里叶变换,Rn表示位于坐标(x,y,z)处的弧形阵列等效采样点与目标Pn(xn,yn,zn)之间的距离;
S53将所述距离频域的回波信号通过匹配滤波获得滤波后的距离频域信号;其中,所使用的滤波函数H(fr)为:
H(fr)=exp{jπKrt2};
其中,Kr为信号的调频率;
距离频域的回波信号经滤波匹配后得到的滤波后的距离频域信号SFT_rc(Rn,x,y,z)表达式为:
Figure QLYQS_7
其中,SFT_re(Rn,x,y,z)表示距离频域的回波信号,H(fr)为滤波函数,δn(xn,yn,zn)为观测场景中目标Pn(xn,yn,zn)的散射系数,fc为雷达工作中心频率,Rn表示位于坐标(x,y,z)处的弧形阵列等效采样点与目标Pn(xn,yn,zn)之间的距离;
S54对所述滤波后的距离频域信号进行相干叠加,获得所述地面场景的微波图像。
2.根据权利要求1所述的滑动式弧形阵列微波成像方法,其特征在于,
在步骤S2中,所述飞行姿态包括所述空中平台的飞行高度,所述微波成像系统的参数包括所述弧形阵列的半径、所述弧形接收阵列天线的入射角或所述弧形发射阵列天线的入射角,以及俯仰向-3dB波束宽度。
3.根据权利要求1所述的滑动式弧形阵列微波成像方法,其特征在于,
所述弧形阵列的中心滑动间距的表达式为:
ΔL≤rmax-rmin
其中,ΔL为所述弧形阵列的中心滑动间距,rmax为所述成像观测最远地距,rmin为成像观测最近地距。
4.根据权利要求1所述的滑动式弧形阵列微波成像方法,其特征在于,
在步骤S54中,对所述滤波后的距离频域信号进行相干叠加包括以下步骤:
S541通过所述散射回波建立所述地面场景对应的图像空间;
S542将所述图像空间分别沿所述弧形阵列的弧形阵列向和所述距离向进行二维离散化,获得二维图像空间和所述二维图像空间的像素个数;
S543计算所述二维图像空间中的每一个像素位置对应的滤波函数,并通过所述滤波函数与所述滤波后的距离频域信号相干叠加,获得所述二维图像空间中的每一个像素的像素值;
S544根据所述二维图像空间中的每一个像素的像素值获得极坐标系下的距离-角度图像数据;
S545将所述极坐标系下的距离-角度图像数据通过插值转换为笛卡尔直角坐标系下的图像数据。
5.根据权利要求4所述的滑动式弧形阵列微波成像方法,其特征在于,
在步骤S3中,所述瞬时环状扫描为所述弧形发射阵列天线辐射电磁波信号对所述地面场景进行360°瞬时环状扫描。
6.根据权利要求1所述的滑动式弧形阵列微波成像方法,其特征在于,
在步骤S1中,所述空中平台为匀速直线运动,
所述滑动式弧形阵列微波成像方法还包括以下步骤:
S6将所述微波图像显示在所述微波成像系统中的显示处理模块上。
CN201811638042.2A 2018-12-29 2018-12-29 滑动式弧形阵列微波成像方法 Active CN109597071B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811638042.2A CN109597071B (zh) 2018-12-29 2018-12-29 滑动式弧形阵列微波成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811638042.2A CN109597071B (zh) 2018-12-29 2018-12-29 滑动式弧形阵列微波成像方法

Publications (2)

Publication Number Publication Date
CN109597071A CN109597071A (zh) 2019-04-09
CN109597071B true CN109597071B (zh) 2023-04-25

Family

ID=65965452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811638042.2A Active CN109597071B (zh) 2018-12-29 2018-12-29 滑动式弧形阵列微波成像方法

Country Status (1)

Country Link
CN (1) CN109597071B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113933834B (zh) * 2021-10-13 2022-07-29 苏州威陌电子信息科技有限公司 圆柱扫描微波成像方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237888A (zh) * 2014-10-20 2014-12-24 内蒙古工业大学 弧形阵列mimo-sar的成像方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5486832A (en) * 1994-07-01 1996-01-23 Hughes Missile Systems Company RF sensor and radar for automotive speed and collision avoidance applications
JP3224768B2 (ja) * 1997-12-11 2001-11-05 池上通信機株式会社 アレーアンテナを具える探査装置
US7663538B2 (en) * 2004-05-18 2010-02-16 Baolab Microsystems S.L. Electromagnetic signal emitting and/or receiving device and corresponding integrated circuit
EP2320247B1 (en) * 2009-11-04 2017-05-17 Rockwell-Collins France A method and system for detecting ground obstacles from an airborne platform
CN102313885B (zh) * 2010-06-30 2013-06-26 中国科学院电子学研究所 多维度微波成像系统及方法
CN102928841A (zh) * 2012-10-30 2013-02-13 西安电子科技大学 一种基于级数反演的机载圆迹环扫sar成像方法
CN103336278A (zh) * 2013-05-13 2013-10-02 西安电子科技大学 多视角观测下前视三维sar成像方法
CN104122552B (zh) * 2014-07-21 2018-01-26 中国科学院电子学研究所 一种滑动式双站圆周合成孔径雷达成像方法
CN104090277B (zh) * 2014-07-21 2016-09-28 中国科学院电子学研究所 一种滑动式圆周合成孔径雷达成像方法
US10078131B2 (en) * 2015-09-15 2018-09-18 Texas Instruments Incorporated Method and apparatus for FMCW radar processing
EP3349664B1 (en) * 2015-09-17 2020-03-11 Koninklijke Philips N.V. Distinguishing lung sliding from external motion
CN109061635B (zh) * 2018-06-08 2021-07-06 深圳市重投华讯太赫兹科技有限公司 旋转扫描系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237888A (zh) * 2014-10-20 2014-12-24 内蒙古工业大学 弧形阵列mimo-sar的成像方法

Also Published As

Publication number Publication date
CN109597071A (zh) 2019-04-09

Similar Documents

Publication Publication Date Title
CN104269658B (zh) 用于mimo-sar成像的弧形阵列天线
EP2866048B1 (en) Tracking device
US8248298B2 (en) Orthogonal linear transmit receive array radar
CN104280735B (zh) 基于弧形阵列天线的mimo-sar成像方法及装置
US20140125511A1 (en) Apparatus and method for assisting vertical takeoff vehicles
CN104267400B (zh) 用于mimo-sar成像的微波信号收发系统、方法及成像系统
CN109597070B (zh) 一种螺旋扫描式弧形阵列微波成像的方法及装置
CN111257878B (zh) 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法
US20190339385A1 (en) Weather radar apparatus
US20170153325A1 (en) System and method for detecting and visualizing targets by airborne radar
CN104090277B (zh) 一种滑动式圆周合成孔径雷达成像方法
CN104122552B (zh) 一种滑动式双站圆周合成孔径雷达成像方法
CN111090094B (zh) 脉冲多普勒雷达的双波束角度测量方法、系统及存储介质
CN109597076A (zh) 用于地基合成孔径雷达的数据处理方法及装置
CN109597071B (zh) 滑动式弧形阵列微波成像方法
CN110412586B (zh) 一种基于相控发射和聚束成像的水下目标探查方法
CN104237888B (zh) 弧形阵列mimo-sar的成像方法
RU2429990C1 (ru) Многофункциональная радиолокационная станция высокого разрешения с активной фазированной решеткой для пилотируемых и беспилотных летательных аппаратов
US20230072466A1 (en) Radar systems and methods
CN112816957A (zh) 一种基于无人机的高俯仰角度散射测试系统
JP5035782B2 (ja) スプリットビーム方式合成開口レーダ
CN114942440B (zh) 对大幅宽区域快速成像的sar二维波束扫描方法及电子设备
EP4019999B1 (en) Radar based system for processing radar data representing range and radial velocity of a detected object
KR20190116865A (ko) 빔 편이 보상 장치 및 방법
WO2019003194A1 (en) 3D SCANNING RADAR AND METHOD OF DETERMINING ALTITUDE OF PRESENT OBJECT

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant