CN113189547A - 一种基于合成带宽频率变标的sar成像的方法及系统 - Google Patents

一种基于合成带宽频率变标的sar成像的方法及系统 Download PDF

Info

Publication number
CN113189547A
CN113189547A CN202110424527.7A CN202110424527A CN113189547A CN 113189547 A CN113189547 A CN 113189547A CN 202110424527 A CN202110424527 A CN 202110424527A CN 113189547 A CN113189547 A CN 113189547A
Authority
CN
China
Prior art keywords
frequency
signal
circuit
distance
echo
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.)
Granted
Application number
CN202110424527.7A
Other languages
English (en)
Other versions
CN113189547B (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.)
Chengdu Huirong Guoke Microsystem Technology Co ltd
Original Assignee
Chengdu Huirong Guoke Microsystem Technology Co ltd
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 Chengdu Huirong Guoke Microsystem Technology Co ltd filed Critical Chengdu Huirong Guoke Microsystem Technology Co ltd
Priority to CN202110424527.7A priority Critical patent/CN113189547B/zh
Publication of CN113189547A publication Critical patent/CN113189547A/zh
Application granted granted Critical
Publication of CN113189547B publication Critical patent/CN113189547B/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
    • 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/28Details of pulse systems
    • G01S7/2813Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
    • 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
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques

Abstract

本发明一种基于合成带宽频率变标的SAR成像的方法及系统,涉及微波成像探测领域。本发明通过基带频域合成技术,向一个重频区域内发射多个窄带线性调频信号,并将重频区域返回的窄带目标回波信号按照基带复信号序列合成宽带回波信号,在不增加固有微波发射端硬件结构的成本基础上,得到宽带回波信号;在不采用运动补偿算法的前提下,通过频率变标波数域成像处理算法来解决了大斜视角下的距离弯曲空变以及距离方位耦合问题,降低了系统的复杂度和响应延迟;本发明还通过频率变标技术,在距离弯曲空变校正上,采用FFT操作和相位函数复数乘法进行弯曲校正,避免了RD算法和PFA算法的插值操作,节省了运算量,能运用于一些实时性要求较高的应用场景。

Description

一种基于合成带宽频率变标的SAR成像的方法及系统
技术领域
本发明涉及微波成像探测领域,尤其涉及一种基于合成带宽频率变标的SAR成像的方法及系统。
背景技术
随着微波探测技术的不断发展,不同形式的成像技术也不断地涌现出来,其中最为耀眼的当属SAR技术。
1978年,美国宇航局NASA发射了全世界第一颗星载SAR卫星:SEASAT卫星;SEASAT卫星的发射宣告了SAR技术的重要性,并向全世界展示了SAR获得高清晰度地表图像的能力,从此以后,SAR变成为了现代微波探测领域及其重要的一环。
SAR(Synthetic Aperture Radar),即合成孔径雷达,是一种主动式的对地观测系统,可安装在飞机、卫星、宇宙飞船等平台上,并能实现全天时、全天候对地观测的能力、且具备一定的地表穿透能力,SAR之所以这么重要,主要是具备以下几个特点:
1.光学遥感能量来自于太阳,因此不能在夜晚工作,而采用SAR的雷达依靠于主动照射目标区域,信号源来自于自己,故能在黑夜中进行探测,并且以电磁波作为信号源可以几乎无失真地穿透水汽云层,因此在不同天气状况下,SAR依然能够正常运行;
2.SAR根据观察路径不同,可分为距离向和方位向,在距离向上的观察,可利用宽带信号脉冲压缩技术实现高分辨成像,而在方位向时,则可通过载体平台相对目标的运动,积累一定的观测积累角度,等效合成出超大尺寸的观测阵列天线,实现对目标方位的超分辨观测;通过距离向和方位向共同作用,可以实现对目标的二维高分辨观测,获得目标几何形状、运动状态等信息;
3.由于图像的分辨率与波长和雷达作用距离无关,通过选择合适的波长,还可以实现对地表遮蔽物的穿透,对光学区无法成像的目标进行高分辨成像,同时由于物质的光学散射能量与雷达电磁散射能量不同,因此雷达和光学传感器具有互补性,有时甚至具有比光学传感器更强的地表特征区分能力。
正是SAR所具有的这些特点,使得SAR在农业、土壤湿度、林业、地质、水文、洪水和海冰监测、海洋学、舰船和浮油探测、冰雪探测、地面森林覆盖测绘、全球高程测量、地球变化检测(如陆地沉陷、火山活动)等方面,具有重大的应用价值,同时在军事情报侦察、地面动目标检测、武器精确制导等方面也具有重要的应用价值。
中国研究SAR技术起步相对较晚,在1970年左右,国内中科院电子所开始了SAR技术研究,并在随后不久开发了属于我国的SAR卫星,目前我国比较著名的SAR系统有尖兵五号(遥感1号)SAR卫星、尖兵七号(遥感6号)、环境一号(由A、B两颗光学小卫星和一颗S波段的SAR小卫星)组成,机载方面应用于我国自主研发的空警-2000预警机的SAR系统,是我国具有代表性的SAR-GMTI系统。
虽然在军事方面,我国SAR技术取得了长足的进步,但是在民用SAR上,还是处于较为匮乏的状态。
民用SAR具备广泛的应用价值,2008年“5.12”地震时期,我国科学院电子所吴一戎院士便带领机载SAR团队,利用开发的机载SAR设备对地震地区进行了SAR成像探测,并获取了第一手的灾情信息,为后续救灾计划的指定提供了很好的信息支持。但是,SAR技术的民用化普及程度比起军用化SAR技术的发展,是远远不如的。
由于SAR技术的敏感性、技术台阶以及成本问题,现有的民用SAR系统也大多集中在星载SAR方面,如美国商业遥感公司ICEYE公司,国内“智星空间”科技有限公司,其主要产品也是星载SAR图像。但星载SAR由于卫星轨道特性,对特定区域访问时间是有窗口的,需要一定的重访时间,不能做到灵活的观测,且星载SAR经济成本非常的高,其产品价格相比光学图像,昂贵很多,这大大限制了民用SAR的广泛应用,本发明还可兼容脉冲窄带信号体制以及调频连续波体制,能够以相对较低的成本实现对目标的高分辨二维成像。
究其原因,可以将SAR民用化发展限制总结为一下几点:
1.由于SAR与光学遥感在在分辨率上还有一定距离,因此需要通过技术手段来对SAR图像的分辨率进行提升,而提升图像分辨率势必会造成硬件成本上升或者系统复杂度增加;
2.实现大范围斜视角工作需要大量的运算能力,而提高平台稳定性需要精确的运动补偿算法,这势必增加研发成本和硬件成本;
3.满足精度等要求的情况下,势必会增加算法的复杂度,从而导致系统的鲁棒性和实时性下降;
4.民用SAR需要考虑实际成本问题,费用不易高昂。
因此,有必要提供一种新的SAR成像的方法及系统来解决上述之一技术问题。
发明内容
为解决上述之一技术问题,本发明提出一种基于合成带宽频率变标的SAR成像的方法及系统。
本发明提供的一种基于合成带宽频率变标的SAR成像的系统,包括微波发射端、微波接收端和信号处理端。
作为更加具体的解决方案,所述微波发射端包括窄带前端、用于上变频的UDC电路、用于中频放大的中频放大电路、用于调频的混频器电路、用于进行高功率放大的HPA电路和发射天线;所述微波接收端包括接收天线、用于进行高增益低噪声放大的LNA电路、用于选频的混频器电路、用于中频放大的中频放大电路和用于下变频的DDC电路;所述信号处理端能对输入信号通过预设的程序方法进行数字化处理,并将处理后的信号以数字化结果进行展示。
作为更加具体的解决方案,所述窄带前端的输出端接在UDC电路输入端上,UDC电路的输出端与中频放大电路的输入端电性连接,混频器电路的输入端与中频放大电路的输出端电性连接,HPA电路的输入端接在混频器电路的输出端上,HPA电路的输出端通过发射天线进行信号输出。
作为更加具体的解决方案,所述接收天线接在LNA电路的输入端上,LNA电路的输出端与混频器电路的输入端电性连接,混频器的输出端接在中频放大电路的输入端上,中频放大电路的输出端与DDC电路的输入端连接,DDC电路的输出端与信号处理端电性连接。
一种基于合成带宽频率变标的SAR成像的方法,运用于上述的基于合成带宽频率变标的SAR成像的系统中,包括根据频率步进信号的进步频率控制微波发射端向重频区域发射窄带发射信号序列;接收重频区域反射回来的目标回波信号并转换成基带复信号序列;通过对基带复信号序列进行频域子带拼接得到宽带回波信号;通过对宽带回波信号进行上变频频载相位补偿得到将原始回波信号;将原始回波信号作为输入信号输入FSA高分辨成像系统进行成像处理,并得到高分辨率SAR图像。
作为更加具体的解决方案,获得的基带复信号序列需按照序列编号存储在不同分区;将基带复信号序列合成为宽带回波信号时,还需要进行多普勒相位补偿;进行频域子带拼接前,还需对基带复信号序列进行升采样和频谱搬移;频域子带拼接是通过滤波器组进行实现;在进行上变频频载相位补偿前,还需要对宽带回波信号进行傅里叶逆变换来得到时域的宽带回波信号。
作为更加具体的解决方案,通过窄带发射信号序列得到基带复信号序列需要进行如下步骤:
通过微波发射端产生窄带发射信号序列;
通过UDC电路将窄带发射信号序列进行上变频;
中频放大电路将上变频后的信号进行选频和中频信号放大,经过中频放大后,信号接入混频器;
混频器的可调本振频率根据频率步进信号的进步频率进行同步改变,并施加于输入信号上;
进过混频器变频后的信号通过HPA电路进行高功率放大,并通过发射天线将信号向重频区域进行广播;
广播后的窄带发射信号序列会被重频区域内的目标进行反射,通过接收天线目标回波信号;
将目标回波信号转入混频器进行选频,接收端混频器的可调本振同样根据频率步进信号的进步频率进行同步改变;
将选频后的信号通过中频放大电路进行放大;
通过DDC电路对放大后的信号进行下变频,从而得到基带复信号序列。
作为更加具体的解决方案,在通过基带复信号序列得到宽带回波信号时需要进行如下步骤:
将得到的基带复信号序列根据序列编号存储在不同分区;
对不同分区的基带复信号序列进行多普勒相位补偿;
通过傅里叶变换对基带复信号序列进行升采样;
根据分区编号,对基带复信号序列进行频谱搬移;
通过滤波器组对分区编号不同的基带复信号序列进行频域子带拼接,并得到宽带回波信号;
对宽带回波信号进行傅里叶逆变换来得到时域的宽带回波信号。
作为更加具体的解决方案,将原始回波信号作为输入信号输入FSA高分辨成像系统进行成像处理时需要进行如下步骤:
将参考信号与原始回波信号共轭相乘,得到解线频调信号;
对解线频调信号的残余相位进行卷积化;
通过变量替换将解线频调信号卷积由距离时间域转化为距离波数域;
通过方位向傅里叶变换将距离波数域信号变换到二维波数域;
对二维波数域信号解析式进行FSA变标处理并得到变标回波信号;
对变标回波信号进行距离傅里叶逆变换,并将滤波器设置为傅里叶逆变的共轭;
通过在时域匹配进行匹配滤波消除式中的置斜处理项;
将去斜后的信号再进行傅里叶变换得到二维数域信号;
二维数域信号通过乘以参考相位函数来进行逆变标,并得到二维数域的逆变标信号;
逆变标信号通过乘以距离走动校正函数进行距离走动校正;
通过二次压缩参考函数对距离走动校正后的逆变标信号进行距离弯曲校正;
通过对距离弯曲校正后的逆变标信号进行距离向傅里叶变换得到距离向处理回波信号;
通过方位向匹配函数对进过距离向处理后的回波信号进行方位匹配压缩;
对经过匹配压缩后的回波信号进行方位傅里叶逆变换,并得到高分辨率回波信号。
作为更加具体的解决方案,所述窄带发射信号序列经过上变频、中频放大和混频器处理后得到发射窄带线性调频信号;多个发射窄带线性调频信号得到实际的发射信号,所述发射信号为:
Figure BDA0003028799410000051
其中,t为快时间,Tp0为窄带发射信号的脉宽,γ1为窄带发射信号的调频率,fc为系统中心频率,i为频率步进变量,N为频率步进数,Δf为频率步进间隔。
作为更加具体的解决方案,所述基带复信号序列表示为:
Figure BDA0003028799410000052
其中,t为快时间,Tp0为窄带发射信号的脉宽,γ1为窄带发射信号的调频率,fc为系统中心频率,i为频率步进变量,N为频率步进数,Δf为频率步进间隔,
Figure BDA0003028799410000053
σ为目标的散射系数,s(i,t,tm)表示目标回波信号。
作为更加具体的解决方案,对不同分区的基带复信号序列进行多普勒相位补偿,补偿后的基带复信号序列为:
Figure BDA0003028799410000061
其中,σ为目标的散射系数,t为快时间,fc为系统中心频率,
Figure BDA0003028799410000062
σ为目标的散射系数,Tp为合成宽带信号的脉宽,γ为合成宽带信号的调频率。
作为更加具体的解决方案,通过对参考信号与原始回波信号进行共轭相乘,得到解线频调信号,所述解线频调信号为:
Figure BDA0003028799410000063
其中,A为常数项,Rref为参考距离,RΔ=Rt-Rref,Tp为合成宽带信号的脉宽,t为快时间,fc为系统中心频率,γ为合成宽带信号的调频率。
与相关技术相比较,本发明提供的用于一种基于合成带宽频率变标的SAR成像方法及系统具有如下有益效果:
1、本发明通过基带频域合成技术,向一个重频区域内发射多个窄带线性调频信号,并将重频区域返回的窄带目标回波信号按照基带复信号序列合成宽带回波信号,在不增加固有微波发射端硬件结构的成本基础上,得到宽带回波信号,提高了系统距离向分辨率;
2、本发明不采用运动补偿算法,而是通过频率变标波数域成像处理算法来解决了大斜视角下的距离弯曲空变以及距离方位耦合问题,使成像算法能适应较大的斜视角下的应用场景;降低了系统的复杂度和响应延迟,增强了系统的实用性;
3、本发明通过频率变标技术,在距离弯曲空变校正上,采用FFT操作和相位函数复数乘法进行弯曲校正,避免了RD算法和PFA算法的插值操作,大大节省了运算量,能运用于一些实时性要求较高的应用场景;
4、本发明还可兼容脉冲窄带信号体制以及调频连续波体制,能够以相对较低的成本实现对目标的高分辨二维成像。
附图说明
图1为本发明提供的窄带发射信号序列转化为基带复信号序列的较佳实施例系统示意图;
图2为本发明提供的基带复信号序列频域子带拼接的较佳实施例系统示意图;
图3为本发明提供的合成原始回波信号的较佳实施例系统示意图;
图4为本发明提供的FSA高分辨成像系统的较佳实施例系统示意图;
图5为本发明提供的较佳实施例运用场景图。
具体实施方式
下面结合附图和实施方式对本发明作进一步说明。
如图1至图5所示,本实施例提供的一种基于合成带宽频率变标的SAR成像方法及系统,为了控制成本,使SAR成像更好的适应民用需求,本实施在采用低成本窄带前端的前提下,通过基带频域合成技术将多个窄带信号合成为大带宽信号,来实现提高SAR成像分辨率的效果。这个过程可以视为:通过窄带前端发射的窄带发射信号序列来得到基带复信号序列,对基带复信号序列进行相位补偿,通过得到的基带复信号序列进行频域子带拼接并得到宽带回波信号。
通过窄带前端发射的窄带发射信号序列来得到基带复信号序列的过程中,我们需要制定宽带合成步进控制策略,通过宽带合成步进控制策略来控制可调本振频率以达到产生与进步频率关联的基带复信号序列的目的,由于进步策略会直接影响SAR成像的时间与分辨率,故需要在实际的场景需求中进行决定,制定好在制定好进步策略后,我们需要将窄带发射信号序列通过UDC电路进行上变频并将上变频后的信号送往中频放大电路,经中频放大电路处理后的窄带发射信号输入混频器,混频器的本振频率由进控制策略来进行控制,混频器的下一级是HPA电路,HPA电路将信号进行高功率放大以后进行发射;原始的发射窄带线性调频信号形式为:
Figure BDA0003028799410000071
其中,t为快时间,Tp0为窄带发射信号的脉宽,γ1为窄带发射信号的调频率,fc为系统中心频率。
但我们向重频区域发射不止一个个窄带线性调频信号,故实际的发射信号形式为:
Figure BDA0003028799410000072
其中,i为慢时间,也表示频率步进变量,N为频率步进数,每N个发射窄带信号合成一个宽带线性调频信号,Δf为频率步进间隔。
在发射信号在遇到目标以后,会反射回目标回波信号,由实际的发射信号形式,我们可以得知目标回波形式应为:
Figure BDA0003028799410000081
其中,
Figure BDA0003028799410000082
c表示光速,R(tm)表示tm时间回波信号所经历的距离;σ为目标的散射系数。
目标回波信号经由接收端进行接收,接收端到的目标回波信号通过LNA电路进行高增益低噪声放大后接入混频器,由于接收端的混频器本振频率同样根据宽带合成步进控制策略进行变化,故同步改变本振频率和滤波器并对目标回波信号进行滤波,经过滤波后的目标回波信号再接由DDC电路进行下变频从而得到基带复信号序列,基带复信号序列可表示为:
Figure BDA0003028799410000083
需要说明的是,由于发射端和采样端中的中频为均为固定频率,基带与中频之间的系统软硬件可以保持相对固定,因此可以到达节约成本和降低系统复杂度的效果。
作为本实施更进一步的解决方案,在得到基带复信号序列后,需对其进行进行多普勒补偿,以此来消除频率步进导致的多普勒变化,补偿后的基带复信号序列可表示为:
Figure BDA0003028799410000084
经过多普勒补偿后,我们需将基带复信号序列按照信号序列进行排列拼接,得到的宽带回波信号形式为:
Figure BDA0003028799410000085
其中,Tp为合成宽带信号的脉宽,γ为合成宽带信号的调频率。
将宽带回波信号在基带乘以上变频相位,就能得到原始回波信号,从而达到了通过窄带前端来获得宽带的回波信号的目的,原始回波信号的形式表示为:
Figure BDA0003028799410000091
在获得原始回波信号后,我们需要对其进行成像,本实施例采用一种FSA高分辨成像处理算法,该算法不要求信号在距离上是线性调频的,而是直接处理距离解线频调后的SAR数据,故能兼顾容脉冲窄带信号体制以及调频连续波体制两种体制,并能能以相对较低的成本实现对目标的高分辨二维成像。
作为本实施更进一步的解决方案,FSA算法是一种改进的线频调变标算法,故首先进行解线频调,解线频调利用LFM信号的特征,将时间固定、频率和调频率相同的LFM信号作为参考信号,将它和原始回波信号做差频处理的方式,也是一种脉冲压缩方式。
设参考距离为Rref,则参考信号为:
Figure BDA0003028799410000092
差频处理,也就是将回波信号与参考信号的共轭相乘:
sif(t,tm)=sr(t,tm)·s* ref(t,tm),
通过上式可以得到解线频调之后的信号为:
Figure BDA0003028799410000093
其中,RΔ=Rt-Rref,可以看出在一个信号周期里,即距离差RΔ为定值,则差频处理后的信号变成了单频信号,差频输出的信号频率为:
Figure BDA0003028799410000094
在得到解线频调之后的信号后,需要进行残余相位卷积化,以消除残余视频相位(RVP):
解线频调残后的信号为:
Figure BDA0003028799410000101
Figure BDA0003028799410000102
其中,
Figure BDA0003028799410000103
Xa和Xac分别为信号方位位置和孔径中心位置,τ为距离时间,Rref为参考距离,新的距离时间下,是以参考距离为快时间零点。第一项和第三项为有用信息,分别用于距离分辨和方位分辨,第二项为剩余视频相位(RVP),实际处理中需要消除。其中瞬时斜距为:
Figure BDA0003028799410000104
大时间带宽条件下,s(Xa,τ)可以表示为卷积形式:
Figure BDA0003028799410000105
Figure BDA0003028799410000106
作为本实施更进一步的解决方案,由于FSA算法需要在波数域进行,故需要通过变量替换将距离时间域转化为距离波数域:
Figure BDA0003028799410000107
Figure BDA0003028799410000108
KR=KRc+ΔKR
Figure BDA0003028799410000109
则信号可以表示为:
Figure BDA0003028799410000111
Figure BDA0003028799410000112
需要说明的是:此处仅是变量替换,不改变原有信号的实际形式。
将信号变换到二维波数域需要对其进行方位向傅里叶变换,由于卷积符号后一项无方位向时间,因此在对方位向傅里叶变换时,可以作为一个常数,不参与积分变换,因此积分形式如下:
Figure BDA0003028799410000113
则信号二维波数域解析形式为:
Figure BDA0003028799410000114
包络函数为:
Figure BDA0003028799410000115
需要说明的是:此处的
Figure BDA0003028799410000116
是一种近似处理:
Figure BDA0003028799410000117
将近似项带入二维波数域解析式,得到:
Figure BDA0003028799410000121
在转换到二维波数域后,还需进行乘以变标函数来进行变标处理:
变标函数为:
Figure BDA0003028799410000122
变标后的回波信号S1(KX,ΔKR)为:
Figure BDA0003028799410000123
其中,A1=A(KX,αAXΔKR)。
作为本实施更进一步的解决方案,在进行完波数域转换和变标处理后,便可进行残余相位解卷积,而在距离波数域的卷积即是距离域的相乘,故对变标后的信号作距离逆傅里叶变换后,残余视频相位便可通过乘以距离参考函数来去除。
距离向H的形式:
Figure BDA0003028799410000124
对H(ΔKR)进行逆傅里叶变换,推导过程如下:
Figure BDA0003028799410000125
故在时域匹配进行匹配滤波便为上式的共轭形式:
Figure BDA0003028799410000131
对信号继续进行距离向傅里叶变换,得到解卷积后的信号二维波数域形式,即为:
S2(KX,ΔKR)=A1exp(jα(AXRref-RB)ΔKR)
Figure BDA0003028799410000132
Figure BDA0003028799410000133
这时,二维波数域的信号形式即卷积符号前的形式,完成了解卷积。
作为本实施更进一步的解决方案,对比变标前后的信号形式可知变标操作引入了新的关于ΔKR的二次相位误差,故在后续操作前应予以去除,在二维波数域乘以参考函数为:
Figure BDA0003028799410000134
二维波数域乘以上述相位函数,即进行逆变标的过程,逆变标后信号的二维波数域为:
S3(KX,ΔKR)=A1exp(jα(AXRref-RB)ΔKR)
Figure BDA0003028799410000135
Figure BDA0003028799410000136
下面进行校正距离走动:对距离走动进行校正,即在二维波数域乘以线性相位函数,距离走动校正函数为:
HRMC(KX,ΔKR)=exp(-jα(AXRref-Rs)ΔKR)
其中,Rs为天线与参考场景中心的最近距离。
完成距离走动校正后,再进行距离弯曲校正,即二次压缩,二次压缩参考函数为:
Figure BDA0003028799410000137
在完成距离向的全部处理后,需通过距离向傅里叶逆变换变回到距离时域。
由上述步骤可知,在距离向傅里叶变换之前的信号形式为:
S4(KX,ΔKR)=A1exp(jKRcRref)exp(-jα(RB-Rs)ΔKR)·exp(-jAXKRcRB-jKXx0)
对上式进行距离向傅里叶变换,可得到距离向处理后的信号回波为:
Figure BDA0003028799410000141
可见,距离上已经将目标点聚焦到了相应的目标位置,分布在参考距离的两侧,如果目标RB=Rs,那么就正好聚焦在参考位置处,故要将回波函数进行高分辨率成像仅需要进行方位向处理即可。
方位向匹配函数为:
HAS(KX,t)=exp(jAXKRcRB)
方位匹配压缩后的信号为:
Figure BDA0003028799410000142
对回波信号做方位傅里叶逆变换便可得到最后的回波图像为:
Figure BDA0003028799410000143
Figure BDA0003028799410000144
上述回波图像为二维矩阵图像;通过对上述回波图像对应像素点,便可得到高分辨率的SAR二维图像。
需要说明的是:本实施例采用低成本窄带前端,使用基带频域合成技术合成大带宽信号,并使用解线频调技术进行距离向压缩,在不增加带宽成本的前提下,提高系统距离向分辨率;距离向采用解线频调(Dechirp)处理,在足够大的带宽下,可实现厘米级的距离向分辨,同时方位上采用频率变标和波数域相结合的成像处理方式,可以实现对空变的距离徙动的精确校正,实现方位向能量的最大化聚焦,实现对目标的高分辨成像;采用频率变标波数域成像处理算法,有效解决了大斜视角下的距离弯曲空变以及距离方位耦合问题,使得成像算法能适应较大的斜视角下的应用场景;距离向采用解线频调(Dechirp)处理,可以适应非常大的信号带宽,可以以较低的采样率对距离向信号采样,大大降低处理数据量;由于采用频率变标技术,对距离弯曲空变校正上,采用FFT操作和相位函数复数乘法就完成了弯曲校正,避免了RD算法和PFA算法的插值操作,大大节省了运算量,整个成像处理算法数据量降低明显,可以应用到一些对实时性要求较高的应用场景;在系统带宽有效的情况下,结合带宽合成技术,可以在前端性能受限的情况下,实现高分辨成像探测;
本实施例给出的设计实现方案基于现有低成本脉冲体制雷达,不会增加多少成本,复杂的信号处理流程全部数字化实现,技术升级的风险和成本转移到了信号处理上,而实际中信号处理提升的成本是不大的,且数字化实现处理精度更高,相比传统的连续波体制,具有更多的优点。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种基于合成带宽频率变标的SAR成像的系统,其特征在于,包括微波发射端、微波接收端和信号处理端;
所述微波发射端包括窄带前端、用于上变频的UDC电路、用于中频放大的中频放大电路、用于调频的混频器电路、用于进行高功率放大的HPA电路和发射天线;
所述微波接收端包括接收天线、用于进行高增益低噪声放大的LNA电路、用于选频的混频器电路、用于中频放大的中频放大电路和用于下变频的DDC电路;
所述信号处理端能对输入信号通过预设的程序方法进行数字化处理,并将处理后的信号以数字化结果进行展示;
所述窄带前端的输出端接在UDC电路输入端上,UDC电路的输出端与中频放大电路的输入端电性连接,混频器电路的输入端与中频放大电路的输出端电性连接,HPA电路的输入端接在混频器电路的输出端上,HPA电路的输出端通过发射天线进行信号输出;
所述接收天线接在LNA电路的输入端上,LNA电路的输出端与混频器电路的输入端电性连接,混频器的输出端接在中频放大电路的输入端上,中频放大电路的输出端与DDC电路的输入端连接,DDC电路的输出端与信号处理端电性连接。
2.一种基于合成带宽频率变标的SAR成像的方法,运用于权利要求1所述的基于合成带宽频率变标的SAR成像的系统中,其特征在于,包括根据频率步进信号的进步频率控制微波发射端向重频区域发射窄带发射信号序列;接收重频区域反射回来的目标回波信号并转换成基带复信号序列;通过对基带复信号序列进行频域子带拼接得到宽带回波信号;通过对宽带回波信号进行上变频频载相位补偿得到将原始回波信号;将原始回波信号作为输入信号输入FSA高分辨成像系统进行成像处理,并得到高分辨率SAR图像。
3.根据权利要求2所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,获得的基带复信号序列需按照序列编号存储在不同分区;将基带复信号序列合成为宽带回波信号时,还需要进行多普勒相位补偿;进行频域子带拼接前,还需对基带复信号序列进行升采样和频谱搬移;频域子带拼接是通过滤波器组进行实现;在进行上变频频载相位补偿前,还需要对宽带回波信号进行傅里叶逆变换来得到时域的宽带回波信号。
4.根据权利要求2所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,通过窄带发射信号序列得到基带复信号序列需要进行如下步骤:
通过微波发射端产生窄带发射信号序列;
通过UDC电路将窄带发射信号序列进行上变频;
中频放大电路将上变频后的信号进行选频和中频信号放大,经过中频放大后,信号接入混频器;
混频器的可调本振频率根据频率步进信号的进步频率进行同步改变,并施加于输入信号上;
进过混频器变频后的信号通过HPA电路进行高功率放大,并通过发射天线将信号向重频区域进行广播;
广播后的窄带发射信号序列会被重频区域内的目标进行反射,通过接收天线目标回波信号;
将目标回波信号转入混频器进行选频,接收端混频器的可调本振同样根据频率步进信号的进步频率进行同步改变;
将选频后的信号通过中频放大电路进行放大;
通过DDC电路对放大后的信号进行下变频,从而得到基带复信号序列。
5.根据权利要求2所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,在通过基带复信号序列得到宽带回波信号时需要进行如下步骤:
将得到的基带复信号序列根据序列编号存储在不同分区;
对不同分区的基带复信号序列进行多普勒相位补偿;
通过傅里叶变换对基带复信号序列进行升采样;
根据分区编号,对基带复信号序列进行频谱搬移;
通过滤波器组对分区编号不同的基带复信号序列进行频域子带拼接,并得到宽带回波信号;
对宽带回波信号进行傅里叶逆变换来得到时域的宽带回波信号。
6.根据权利要求2所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,将原始回波信号作为输入信号输入FSA高分辨成像系统进行成像处理时需要进行如下步骤:
将参考信号与原始回波信号共轭相乘,得到解线频调信号;
对解线频调信号的残余相位进行卷积化;
通过变量替换将解线频调信号卷积由距离时间域转化为距离波数域;
通过方位向傅里叶变换将距离波数域信号变换到二维波数域;
对二维波数域信号解析式进行FSA变标处理并得到变标回波信号;
对变标回波信号进行距离傅里叶逆变换,并将滤波器设置为傅里叶逆变的共轭;
通过在时域匹配进行匹配滤波消除式中的置斜处理项;
将去斜后的信号再进行傅里叶变换得到二维数域信号;
二维数域信号通过乘以参考相位函数来进行逆变标,并得到二维数域的逆变标信号;
逆变标信号通过乘以距离走动校正函数进行距离走动校正;
通过二次压缩参考函数对距离走动校正后的逆变标信号进行距离弯曲校正;
通过对距离弯曲校正后的逆变标信号进行距离向傅里叶变换得到距离向处理回波信号;
通过方位向匹配函数对进过距离向处理后的回波信号进行方位匹配压缩;
对经过匹配压缩后的回波信号进行方位傅里叶逆变换,并得到高分辨率回波信号。
7.根据权利要求4所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,所述窄带发射信号序列经过上变频、中频放大和混频器处理后得到发射窄带线性调频信号;多个发射窄带线性调频信号得到实际的发射信号,所述发射信号为:
Figure FDA0003028799400000031
其中,t为快时间,Tp0为窄带发射信号的脉宽,γ1为窄带发射信号的调频率,fc为系统中心频率,i为频率步进变量,N为频率步进数,Δf为频率步进间隔。
8.根据权利要求4所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,所述基带复信号序列表示为:
Figure FDA0003028799400000032
其中,t为快时间,Tp0为窄带发射信号的脉宽,γ1为窄带发射信号的调频率,fc为系统中心频率,i为频率步进变量,N为频率步进数,Δf为频率步进间隔,
Figure FDA0003028799400000041
σ为目标的散射系数,s(i,t,tm)表示目标回波信号。
9.根据权利要求5所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,对不同分区的基带复信号序列进行多普勒相位补偿,补偿后的基带复信号序列为:
Figure FDA0003028799400000042
其中,σ为目标的散射系数,t为快时间,fc为系统中心频率,
Figure FDA0003028799400000043
σ为目标的散射系数,Tp为合成宽带信号的脉宽,γ为合成宽带信号的调频率。
10.根据权利要求6所述的一种基于合成带宽频率变标的SAR成像的方法,其特征在于,通过对参考信号与原始回波信号进行共轭相乘,得到解线频调信号,所述解线频调信号为:
Figure FDA0003028799400000044
其中,A为常数项,Rref为参考距离,RΔ=Rt-Rref,Tp为合成宽带信号的脉宽,t为快时间,fc为系统中心频率,γ为合成宽带信号的调频率。
CN202110424527.7A 2021-04-20 2021-04-20 一种基于合成带宽频率变标的sar成像的方法及系统 Active CN113189547B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110424527.7A CN113189547B (zh) 2021-04-20 2021-04-20 一种基于合成带宽频率变标的sar成像的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110424527.7A CN113189547B (zh) 2021-04-20 2021-04-20 一种基于合成带宽频率变标的sar成像的方法及系统

Publications (2)

Publication Number Publication Date
CN113189547A true CN113189547A (zh) 2021-07-30
CN113189547B CN113189547B (zh) 2023-08-08

Family

ID=76977560

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110424527.7A Active CN113189547B (zh) 2021-04-20 2021-04-20 一种基于合成带宽频率变标的sar成像的方法及系统

Country Status (1)

Country Link
CN (1) CN113189547B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114780911A (zh) * 2022-03-18 2022-07-22 西北工业大学 一种基于深度学习的海洋宽测绘带距离解模糊方法

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05265061A (ja) * 1992-03-19 1993-10-15 Fujitsu Ltd 光増幅器多段中継伝送路
GB9004455D0 (en) * 1990-02-28 1996-04-24 Marconi Gec Ltd Radar system
EP0911659A1 (en) * 1997-10-27 1999-04-28 Lucent Technologies Inc. Optical add/drop multiplexer having complementary stages
US20050007269A1 (en) * 2003-07-09 2005-01-13 Veridian Systems Method and system for providing along-track alignment and formatting of synthetic aperture radar (sar) data, and sar image formation algorithms using such method and system
CN101021561A (zh) * 2007-04-06 2007-08-22 清华大学 一种采用多发多收频分信号的宽带雷达及其成像方法
CN102435991A (zh) * 2011-12-16 2012-05-02 天津职业技术师范大学 天线半圆弧运动提高地面雷达方位分辨率系统
US20120146846A1 (en) * 2009-06-05 2012-06-14 Paul Antonik Method and apparatus for simultaneous synthetic aperture radar and moving target indication
CN102608600A (zh) * 2012-04-09 2012-07-25 西安电子科技大学 基于fpga的步进频率像拼接的实现方法
EP2650695A1 (en) * 2012-08-02 2013-10-16 Institute of Electronics, Chinese Academy of Sciences Imaging method for synthetic aperture radar in high squint mode
CN103472449A (zh) * 2013-09-16 2013-12-25 电子科技大学 一种基于mimo图像域的bp宽带合成方法
CN103716262A (zh) * 2012-10-09 2014-04-09 王晓安 基于时域参数提取的信道估计方法
CN103969644A (zh) * 2013-01-31 2014-08-06 中国人民解放军国防科学技术大学 一种多通道调频连续波sar成像方法
CN107918124A (zh) * 2017-10-26 2018-04-17 西安电子科技大学 带有方位空变校正的机载大斜视高分辨sar成像方法
CN108132467A (zh) * 2017-12-23 2018-06-08 成都汇蓉国科微系统技术有限公司 基于增强型adc的dsp+fpga双基前视sar成像方法及成像装置
KR20180115935A (ko) * 2017-04-14 2018-10-24 연세대학교 산학협력단 단일 경로와 단일 안테나를 이용한 영상 레이다의 고도 측정 시스템 및 방법
CN109613510A (zh) * 2019-01-21 2019-04-12 杭州光预科技有限公司 用于小目标探测或实时跟踪成像的微波光子雷达实现方法及系统
CN110501706A (zh) * 2019-08-20 2019-11-26 中国人民解放军国防科技大学 大角度非均匀转动空间目标isar成像方法
CN111308579A (zh) * 2020-03-20 2020-06-19 南京大桥机器有限公司 一种l波段窄带电子探空仪

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9004455D0 (en) * 1990-02-28 1996-04-24 Marconi Gec Ltd Radar system
JPH05265061A (ja) * 1992-03-19 1993-10-15 Fujitsu Ltd 光増幅器多段中継伝送路
EP0911659A1 (en) * 1997-10-27 1999-04-28 Lucent Technologies Inc. Optical add/drop multiplexer having complementary stages
US20050007269A1 (en) * 2003-07-09 2005-01-13 Veridian Systems Method and system for providing along-track alignment and formatting of synthetic aperture radar (sar) data, and sar image formation algorithms using such method and system
CN101021561A (zh) * 2007-04-06 2007-08-22 清华大学 一种采用多发多收频分信号的宽带雷达及其成像方法
US20120146846A1 (en) * 2009-06-05 2012-06-14 Paul Antonik Method and apparatus for simultaneous synthetic aperture radar and moving target indication
CN102435991A (zh) * 2011-12-16 2012-05-02 天津职业技术师范大学 天线半圆弧运动提高地面雷达方位分辨率系统
CN102608600A (zh) * 2012-04-09 2012-07-25 西安电子科技大学 基于fpga的步进频率像拼接的实现方法
EP2650695A1 (en) * 2012-08-02 2013-10-16 Institute of Electronics, Chinese Academy of Sciences Imaging method for synthetic aperture radar in high squint mode
CN103716262A (zh) * 2012-10-09 2014-04-09 王晓安 基于时域参数提取的信道估计方法
CN103969644A (zh) * 2013-01-31 2014-08-06 中国人民解放军国防科学技术大学 一种多通道调频连续波sar成像方法
CN103472449A (zh) * 2013-09-16 2013-12-25 电子科技大学 一种基于mimo图像域的bp宽带合成方法
KR20180115935A (ko) * 2017-04-14 2018-10-24 연세대학교 산학협력단 단일 경로와 단일 안테나를 이용한 영상 레이다의 고도 측정 시스템 및 방법
CN107918124A (zh) * 2017-10-26 2018-04-17 西安电子科技大学 带有方位空变校正的机载大斜视高分辨sar成像方法
CN108132467A (zh) * 2017-12-23 2018-06-08 成都汇蓉国科微系统技术有限公司 基于增强型adc的dsp+fpga双基前视sar成像方法及成像装置
CN109613510A (zh) * 2019-01-21 2019-04-12 杭州光预科技有限公司 用于小目标探测或实时跟踪成像的微波光子雷达实现方法及系统
CN110501706A (zh) * 2019-08-20 2019-11-26 中国人民解放军国防科技大学 大角度非均匀转动空间目标isar成像方法
CN111308579A (zh) * 2020-03-20 2020-06-19 南京大桥机器有限公司 一种l波段窄带电子探空仪

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
HAIWEN MEI 等: "A New Nonlinear Chirp Scaling Algorithm for Translational Variant Bistatic Forward-Looking SAR with Dive Trajectory", 《2019 6TH ASIA-PACIFIC CONFERENCE ON SYNTHETIC APERTURE RADAR (APSAR)》 *
RULIANG YANG 等: "Spotlight SAR Imaging Algorithms", 《HIGH-RESOLUTION MICROWAVE IMAGING》 *
孟自强 等: "弹载双基前视SAR俯冲段弹道设计方法", 《系统工程与电子技术》 *
张俊: "复杂运动条件下的SAR/ISAR聚焦方法研究", 《中国博士学位论文全文数据库 信息科技辑》 *
徐宗志: "基于实测SAR回波数据的多子带接收合成技术", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
邢孟道 等: "多平台合成孔径雷达成像算法综述", 《雷达学报》 *
陈琦 等: "高分三号卫星聚束SAR误差补偿成像处理技术研究", 《航天器工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114780911A (zh) * 2022-03-18 2022-07-22 西北工业大学 一种基于深度学习的海洋宽测绘带距离解模糊方法
CN114780911B (zh) * 2022-03-18 2024-03-01 西北工业大学 一种基于深度学习的海洋宽测绘带距离解模糊方法

Also Published As

Publication number Publication date
CN113189547B (zh) 2023-08-08

Similar Documents

Publication Publication Date Title
EP2998763B1 (en) Phase calibration of a stepped-chirp signal for a synthetic aperture radar
EP3144702B1 (en) Method and device for synthethic aperture radar imaging based on non-linear frequency modulation signal
EP2896971B1 (en) Spaceborne Multi-Channel Synthetic Aperture Radar Imaging Device
EP2650695B1 (en) Imaging method for synthetic aperture radar in high squint mode
CN105785327B (zh) 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
CN107741586B (zh) 基于DBF-TOPS加权的星载Ka InSAR信号处理方法
CN103454632B (zh) 一站固定式调频连续波双基地sar成像方法
Li et al. Modeling and processing of two-dimensional spatial-variant geosynchronous SAR data
CN103336280B (zh) 临近空间慢速平台合成孔径雷达大场景成像方法
CN106932778B (zh) 方位多通道fmcw sar滑动聚束成像方法
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
CN110632594A (zh) 一种长波长星载sar成像方法
CN110187347B (zh) 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法
CN106338731A (zh) 一种多子带并发的mimo‑sar雷达成像方法
CN106950565A (zh) 星载sar成像抖动补偿方法、成像方法
CN113075660B (zh) 基于极化合成孔径雷达sar反演海面风浪参数的方法及装置
CN107102328B (zh) 基于fpga的实时成像信号处理方法及fpga
CN104062657A (zh) 一种合成孔径雷达sar的广义极坐标成像方法
CN108008387B (zh) 一种机载阵列天线下视三维成像方法
CN116299462A (zh) 一种低轨卫星辐射源的被动双基前视sar聚束成像方法
CN103278819B (zh) 基于滑动接收窗的机载高分辨率斜视聚束sar成像方法
CN113189547B (zh) 一种基于合成带宽频率变标的sar成像的方法及系统
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN113702974A (zh) 一种机载/弹载合成孔径雷达图像快速优化方法
CN116794612A (zh) 基于线性变标的合成孔径雷达距离-孔径空变运动补偿方法

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