CN109828271A - Sar成像方法及装置,sar成像系统 - Google Patents
Sar成像方法及装置,sar成像系统 Download PDFInfo
- Publication number
- CN109828271A CN109828271A CN201811648144.2A CN201811648144A CN109828271A CN 109828271 A CN109828271 A CN 109828271A CN 201811648144 A CN201811648144 A CN 201811648144A CN 109828271 A CN109828271 A CN 109828271A
- Authority
- CN
- China
- Prior art keywords
- matched filtering
- filtering coefficient
- indicates
- echo data
- frequency domain
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种SAR成像方法及装置,SAR成像系统。其中,该方法包括:获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。本发明解决了相关技术中SAR过程中的匹配滤波实现可靠性较低的技术问题。
Description
技术领域
本发明涉及雷达成像技术领域,具体而言,涉及一种SAR成像方法及装置,SAR成像系统。
背景技术
SAR(合成孔径雷达)是一种高分辨率成像雷达,可以在能见度极低的气象条件下得到光学照相的高分辨雷达图像。利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成较大的等效天线孔径的雷达,也称综合孔径雷达。合成孔径雷达的特点是分辨率高,能全天候工作,能够有效地识别伪装和穿透掩盖物。而对于雷达,例如,合成孔径雷达的回波数据一般情况下,需要进行滤波后才能应用。现有技术中对回波数据的滤波方式一般是采用在FPGA(现场可编程门阵列)中利用欧拉公式计算匹配滤波函数的函数值SR,并利用计算出的SR与回波数据进行复乘实现匹配滤波。计算SR时会调用CORDIC核并实时计算SR中的sin值和cos值等参数,会使用大量的逻辑资源且延时大,计算SR过程中浪费DSP资源,需要在各次运算(例如乘法、除法、乘方)中对位长进行截位,这会容易引入较大的误差。
针对上述相关技术中SAR成像过程中的匹配滤波实现可靠性较低的问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种SAR成像方法及装置,SAR成像系统,以至少解决相关技术中SAR过程中的匹配滤波实现可靠性较低的技术问题。
根据本发明实施例的一个方面,提供了一种SAR成像方法,包括:获取频域回波数据,其中,所述频域回波数据是对原始回波数据进行预处理后得到的数据,所述原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;将所述匹配滤波系数与所述频域回波数据进行复乘,以对所述频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
可选地,在从FPGA的存储设备中读取匹配滤波系数之前,该SAR成像方法还包括:确定所述匹配滤波系数;其中,确定所述匹配滤波系数包括:利用预定软件确定匹配滤波函数的函数值,其中,所述预定软件包括:MATLAB;基于所述函数值得到所述匹配滤波系数。
可选地,所述匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
可选地,确定所述匹配滤波系数还包括:利用欧拉公式将所述匹配滤波函数转换为第一公式,根据所述第一公式确定所述匹配滤波系数,其中,所述第一公式为:SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;提取所述第一公式中的实部和虚部,将所述实部和所述虚部确定为所述匹配滤波系数。
可选地,在确定所述匹配滤波系数之后,该SAR成像方法还包括:将所述匹配滤波系数转化为所述FPGA的存储设备可识别的预定格式,得到处理后的匹配滤波系数;将所述处理后的匹配滤波系数进行分别存储。
根据本发明实施例的另外一个方面,还提供了一种SAR成像装置,包括:获取单元,用于获取频域回波数据,其中,所述频域回波数据是对原始回波数据进行预处理后得到的数据,所述原始回波数据是SAR接收到的回波数据;读取单元,用于从FPGA的存储设备中读取匹配滤波系数;处理单元,用于将所述匹配滤波系数与所述频域回波数据进行复乘,以对所述频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
可选地,该SAR成像装置还包括:确定单元,用于在从FPGA的存储设备中读取匹配滤波系数之前,确定所述匹配滤波系数;其中,所述确定单元包括:确定模块,用于利用预定软件确定匹配滤波函数的函数值,其中,所述预定软件包括:MATLAB;获取模块,用于基于所述函数值得到所述匹配滤波系数。
可选地,所述匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
可选地,所述确定模块还包括:转换子模块,用于利用欧拉公式将所述匹配滤波函数转换为第一公式,根据所述第一公式确定所述匹配滤波系数,其中,所述第一公式为:SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;确定子模块用于提取所述第一公式中的实部和虚部,将所述实部和所述虚部确定为所述匹配滤波系数。
可选地,该SAR成像装置还包括:转化模块,用于在确定所述匹配滤波系数之后,将所述匹配滤波系数转化为所述现场可编程阵列FPGA的存储设备可识别的预定格式,得到处理后的匹配滤波系数;存储模块,用于将所述处理后的匹配滤波系数进行分别存储。
根据本发明实施例的另外一个方面,还提供了一种SAR成像系统,所述SAR成像系统使用上述中任意一项所述的SAR成像方法。
根据本发明实施例的另外一个方面,还提供了一种存储介质,所述存储介质包括存储的程序,其中,所述程序执行上述中任意一项所述的SAR成像方法。
根据本发明实施例的另外一个方面,还提供了一种处理器,所述处理器用于运行程序,其中,所述程序运行时执行上述中任意一项所述的SAR成像方法。
在本发明实施例中,采用获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。通过本发明实施例提供的SAR成像方法可以实现在进行SAR成像时,利用从预定位置读取匹配滤波系数与频域回波数据进行复乘,以对二维相位图谱进行纠正的目的,而不需要实时计算匹配滤波系数,从而达到了提高SAR雷达成像的精准度的技术效果,进而解决了相关技术中SAR过程中的匹配滤波实现可靠性较低的技术问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的SAR成像方法的流程图;
图2是根据本发明实施例的SAR成像装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例1
根据本发明实施例,提供了一种SAR成像方法的方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的SAR成像方法的流程图,如图1所示,该SAR成像方法包括如下步骤:
步骤S102,获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据。
在步骤S102中,在SAR发射出雷达信号后,会反射回回波信号,对回波信号进行模数转换后得到原始回波数据,在对原始回波数据进行预处理后,将时域回波数据转换为频域回波数据,当频域回波数据到来时会伴随着valid信号,当FPGA(现场可编程阵列)接收到valid信号时,表示即将有待进行滤波的频域回波数据。
步骤S104,从FPGA的存储设备中读取匹配滤波系数。
步骤S106,将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
通过上述步骤,可以在获取频域回波数据之后,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;再将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。相对于相关技术中在需要进行SAR成像过程中的匹配滤波实现时,需要实时计算得到匹配滤波系数,导致的会使用大量的逻辑资源且延时较大,以及在各次运算中对位长进行截位时容易引入较大误差的弊端,通过本发明实施例提供的SAR成像方法可以实现在进行SAR成像时,利用从预定位置读取匹配滤波系数与频域回波数据进行复乘,以对二维相位图谱进行纠正的目的,而不需要实时计算匹配滤波系数,从而达到了提高SAR成像的精准度的技术效果,进而解决了相关技术中SAR过程中的匹配滤波实现可靠性较低的技术问题。
作为一种可选的实施例,在从FPGA的存储设备中读取匹配滤波系数之前,该SAR成像方法还可以包括:确定匹配滤波系数;其中,确定匹配滤波系数包括:利用预定软件确定匹配滤波函数的函数值,其中,预定软件包括:MATLAB;基于函数值得到匹配滤波系数。即,通过本发明实施例节省掉了计算函数值时的实时运算,提前在MATLAB软件或其他的计算工具中完成函数值的计算后,将计算后得到的匹配滤波系数(即,函数值)存储起来,在需要的时候直接从预定存储位置直接读取需要的匹配滤波系数即可。
作为一种可选的实施例,在确定匹配滤波系数之后,由于从MATLAB软件计算出的滤波参数为浮点型数据,而FPGA中进行滤波时所需要的是定点型数据,因此,还需要对得到的匹配滤波系数进行截位,例如,将匹配滤波系数保留到8位,或16位,当然也可以为其他位数。
其中,对得到的匹配滤波系数进行截位的时机可以是在将滤波参数存储到FPGA之前,即,按照一定的顺序将截位后的匹配滤波系数存储在存储设备中;也可以是,在将匹配滤波系数存储到FPGA之后,即,在将匹配滤波系数从存数设备中读取出来后,进行截位。相对于相关技术中使用定点运算过程中需要在各个运算环节进行截位而引入的误差,MATLAB进行计算结果为浮点型,在最终的计算结果处进行截位会使得运算结果精度更高。
其中,匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
优选的,确定匹配滤波系数还可以包括:利用欧拉公式将匹配滤波函数转换为第一公式,根据第一公式确定匹配滤波系数,其中,第一公式为:SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),SR为函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;提取第一公式中的实部和虚部,将实部和虚部确定为匹配滤波系数。例如,在需要处理的频域回波数据的数量比较多的情况下,利用匹配滤波函数SR=exp(j*pi*tp/B*fa2),得到的匹配滤波系数,在进行存储时会浪费大量的逻辑资源且时延较大;此时,可以基于欧拉公式将匹配滤波函数转换为SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),然后获取到第一公式的实部和虚部,仅仅需要存储提取到的第一公式的实部和虚部即可,从而可以减小时延,节约逻辑资源以及DSP(数字信号处理)资源。
具体地,目标距离fa可以通过第二公式确定,其中,第二公式为:fa=linspace(0,Fs,nrn),Fs表示频域回波数据的采样频率,nrn表示距离向的采样点数。
作为一种可选的实施例,在确定匹配滤波系数之后,该SAR成像方法还可以包括:将匹配滤波系数转化为FPGA的存储设备可识别的预定格式,得到处理后的匹配滤波系数;将处理后的匹配滤波系数进行分别存储。
通过本发明实施例提供的SAR成像方法中节省了计算滤波器中匹配滤波函数的函数值的实时运算,提前在MATLAB软件或其他的计算工具中完成函数值的计算后,将计算后的匹配滤波系数存储在FPGA的存储设备中(例如,ROM),在需要对待滤波的频域回波数据进行滤波时,实时从存储设备中读取匹配滤波系数即可。
其中,在FPGA进行计算时,需要对匹配滤波系数进行截位,将经过在MATLAB软件中计算得到的匹配滤波系数由浮点型数据转换为定点型数据,由此可以得到存储设备可以识别的匹配滤波系数,然后按照一定的顺序存储在存储设备中,在需要时实时取出使用即可,从而实现了节省逻辑资源和DSP资源的目的,并且对待进行滤波的频域回波数据进行滤波的时延较小,逻辑控制相对比较简单。
另外,在本发明实施例的另外一个方面,还提供了一种SAR成像系统,该回波数据的滤波系统使用上述中任意一项的SAR成像方法。该SAR成像系统可以应用于SAR中,当SAR设计好后,滤波器的对应的参数就会确定,因此可以利用计算工具(例如,MATLAB软件)将滤波器的函数值计算出来后,存储在FPGA的预定存储位置中,在需要时直接读取即可实现将其与待进行滤波的频域回波数据进行复乘,以对待进行滤波的频域回波数据进行修正。需要说明的是,在将匹配滤波系数与频域回波数据进行复乘时需要注意时序关系。
这种直接将匹配滤波系数直接存储代替实时运算的方式可以是在MATLAB软件中求出匹配滤波函数的表达式的函数值(匹配滤波函数存储在滤波器中),并将计算出的匹配滤波系数的实部和虚部分别取出,并分别存储在两个ROM中。
另外,这种直接存储的方式节省掉了繁杂的数学运算过程,省掉了cordic核的使用,节省了逻辑资源和DSP资源,控制简单,延时小,并且避免了定点运算过程中各个运算环节进行截位而引入的误差,MATLAB进行的计算结果为浮点型,在最终的结果处进行截位运算结果精度更高。控制逻辑更简单。
实施例2
根据本发明实施例还提供了一种SAR成像装置,需要说明的是,本发明实施例的SAR成像装置可以用于执行本发明实施例所提供的SAR成像方法。以下对本发明实施例提供的SAR成像装置进行介绍。
图2是根据本发明实施例的SAR成像装置的示意图,如图2所示,该SAR成像装置包括:获取单元21,读取单元23以及处理单元25。下面对该SAR成像装置进行详细说明。
获取单元21,用于获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据。
读取单元23,用于从FPGA的存储设备中读取匹配滤波系数。
处理单元25,用于将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
在该实施例中,可以利用获取单元21获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;然后利用读取单元23从FPGA的存储设备中读取匹配滤波系数;再利用处理单元25将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。相对于相关技术中在需要进行SAR成像过程中的匹配滤波实现时,需要实时计算得到匹配滤波系数,导致的会使用大量的逻辑资源且延时较大,以及在各次运算中对位长进行截位时容易引入较大误差的弊端,通过本发明实施例提供的SAR成像装置可以实现在进行SAR成像时,利用从预定位置读取匹配滤波系数与频域回波数据进行复乘,以对二维相位图谱进行纠正的目的,而不需要实时计算匹配滤波系数,从而达到了提高SAR成像的精准度的技术效果,进而解决了相关技术中SAR过程中的匹配滤波实现可靠性较低的技术问题。
作为一种可选的实施例,该SAR成像装置还可以包括:确定单元,用于在从FPGA的存储设备中读取匹配滤波系数之前,确定匹配滤波系数;其中,确定单元包括:确定模块,用于利用预定软件确定匹配滤波函数的函数值,其中,预定软件包括:MATLAB;获取模块,用于基于函数值得到匹配滤波系数。
作为一种可选的实施例,匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
作为一种可选的实施例,确定模块还包括:转换子模块,用于利用欧拉公式将匹配滤波函数转换为第一公式,根据第一公式确定匹配滤波系数,其中,第一公式为:SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),SR为函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;确定子模块用于提取第一公式中的实部和虚部,将实部和虚部确定为匹配滤波系数。
作为一种可选的实施例,该SAR成像装置还包括:转化模块,用于在确定匹配滤波系数之后,将匹配滤波系数转化为FPGA的存储设备可识别的预定格式,得到处理后的匹配滤波系数;存储模块,用于将处理后的匹配滤波系数进行分别存储。
上述SAR成像装置包括处理器和存储器,上述获取单元21,读取单元23以及处理单元25等均作为程序单元存储在存储器中,由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
上述处理器中包含内核,由内核去存储器中调取相应的程序单元。内核可以设置一个或以上,通过调整内核参数将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
上述存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM),存储器包括至少一个存储芯片。
根据本发明实施例的另外一个方面,还提供了一种存储介质,存储介质包括存储的程序,其中,程序执行上述中任意一项的SAR成像方法。
根据本发明实施例的另外一个方面,还提供了一种处理器,处理器用于运行程序,其中,程序运行时执行上述中任意一项的SAR成像方法。
在本发明实施例中还提供了一种设备,该设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序,处理器执行程序时实现以下步骤:获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
在本发明实施例中还提供了一种计算机程序产品,当在数据处理设备上执行时,适于执行初始化有如下方法步骤的程序:获取频域回波数据,其中,频域回波数据是对原始回波数据进行预处理后得到的数据,原始回波数据是SAR接收到的回波数据;从FPGA的存储设备中读取匹配滤波系数;将匹配滤波系数与频域回波数据进行复乘,以对频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种SAR成像方法,其特征在于,包括:
获取频域回波数据,其中,所述频域回波数据是对原始回波数据进行预处理后得到的数据,所述原始回波数据是SAR接收到的回波数据;
从FPGA的存储设备中读取匹配滤波系数;
将所述匹配滤波系数与所述频域回波数据进行复乘,以对所述频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
2.根据权利要求1所述的方法,其特征在于,在从FPGA的存储设备中读取匹配滤波系数之前,还包括:确定所述匹配滤波系数;
其中,确定所述匹配滤波系数包括:
利用预定软件确定匹配滤波函数的函数值,其中,所述预定软件包括:MATLAB;
基于所述函数值得到所述匹配滤波系数。
3.根据权利要求2所述的方法,其特征在于,所述匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
4.根据权利要求3所述的方法,其特征在于,确定所述匹配滤波系数还包括:
利用欧拉公式将所述匹配滤波函数转换为第一公式,根据所述第一公式确定所述匹配滤波系数,其中,所述第一公式为:SR=cos(pi*tp/B*fa2)+jsin(pi*tp/B*fa2),SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;
提取所述第一公式中的实部和虚部,将所述实部和所述虚部确定为所述匹配滤波系数。
5.根据权利要求4所述的方法,其特征在于,在确定所述匹配滤波系数之后,还包括:
将所述匹配滤波系数转化为所述FPGA的存储设备可识别的预定格式,得到处理后的匹配滤波系数;
将所述处理后的匹配滤波系数进行分别存储。
6.一种SAR成像装置,其特征在于,包括:
获取单元,用于获取频域回波数据,其中,所述频域回波数据是对原始回波数据进行预处理后得到的数据,所述原始回波数据是SAR接收到的回波数据;
读取单元,用于从FPGA的存储设备中读取匹配滤波系数;
处理单元,用于将所述匹配滤波系数与所述频域回波数据进行复乘,以对所述频域回波数据对应的二维图谱进行修正,并基于修正后的二维图谱进行成像。
7.根据权利要求6所述的装置,其特征在于,还包括:确定单元,用于在从FPGA的存储设备中读取匹配滤波系数之前,确定所述匹配滤波系数;
其中,所述确定单元包括:
确定模块,用于利用预定软件确定匹配滤波函数的函数值,其中,所述预定软件包括:MATLAB;
获取模块,用于基于所述函数值得到所述匹配滤波系数。
8.根据权利要求7所述的装置,其特征在于,所述匹配滤波函数为:SR=exp(j*pi*tp/B*fa2),其中,SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离。
9.根据权利要求8所述的装置,其特征在于,所述确定模块还包括:
转换子模块,用于利用欧拉公式将所述匹配滤波函数转换为第一公式,根据所述第一公式确定所述匹配滤波系数,其中,所述第一公式为:SR=cos(pi*tpB*fa2+jsin(pi*tpB*fa2),SR为所述函数值,j表示虚数单位,pi表示常数π,tp表示回波数据的采样时间,B表示带宽,fa表示目标距离;
确定子模块用于提取所述第一公式中的实部和虚部,将所述实部和所述虚部确定为所述匹配滤波系数。
10.一种SAR成像系统,其特征在于,所述SAR成像系统使用权利要求1至5中任意一项所述的SAR成像方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811648144.2A CN109828271A (zh) | 2018-12-29 | 2018-12-29 | Sar成像方法及装置,sar成像系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811648144.2A CN109828271A (zh) | 2018-12-29 | 2018-12-29 | Sar成像方法及装置,sar成像系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109828271A true CN109828271A (zh) | 2019-05-31 |
Family
ID=66859992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811648144.2A Pending CN109828271A (zh) | 2018-12-29 | 2018-12-29 | Sar成像方法及装置,sar成像系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109828271A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941243A (zh) * | 2014-04-03 | 2014-07-23 | 电子科技大学 | 一种基于sar三维成像的自旋式飞行器测高方法 |
CN106130507A (zh) * | 2016-06-17 | 2016-11-16 | 中国电子科技集团公司第十研究所 | 高速并行位同步滤波器时变系数更新环路 |
CN106772369A (zh) * | 2016-12-16 | 2017-05-31 | 北京华航无线电测量研究所 | 一种基于多视角成像的三维成像方法 |
CN108802726A (zh) * | 2017-12-29 | 2018-11-13 | 西安电子科技大学 | 基于图形处理器gpu的合成孔径雷达成像方法 |
-
2018
- 2018-12-29 CN CN201811648144.2A patent/CN109828271A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941243A (zh) * | 2014-04-03 | 2014-07-23 | 电子科技大学 | 一种基于sar三维成像的自旋式飞行器测高方法 |
CN106130507A (zh) * | 2016-06-17 | 2016-11-16 | 中国电子科技集团公司第十研究所 | 高速并行位同步滤波器时变系数更新环路 |
CN106772369A (zh) * | 2016-12-16 | 2017-05-31 | 北京华航无线电测量研究所 | 一种基于多视角成像的三维成像方法 |
CN108802726A (zh) * | 2017-12-29 | 2018-11-13 | 西安电子科技大学 | 基于图形处理器gpu的合成孔径雷达成像方法 |
Non-Patent Citations (2)
Title |
---|
无: "雷达成像技术实验报告", 《HTTPS://WENKU.BAIDU.COM/VIEW/2A3A6A7D31B765CE0508142A.HTML》 * |
苏斌 等: "基于FPGA的脉冲压缩处理器设计与实现", 《电子测量技术》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1740972B1 (en) | Sub-aperture sidelobe and alias mitigation techniques | |
KR102418200B1 (ko) | 레이더, 특히, 합성 개구 레이더의 수신 데이터의 감소 | |
CN103149561B (zh) | 一种基于场景块稀疏的稀疏微波成像方法 | |
WO1997011387A1 (en) | A sar radar system | |
DE102006006572A1 (de) | Paarweise ZF-Abtastung für Puls-Laufzeit-Füllstandsensoren | |
US20190227164A1 (en) | Automotive Testing Method, System and Computer Program Product | |
Armanious et al. | An adversarial super-resolution remedy for radar design trade-offs | |
CN109444839A (zh) | 目标轮廓的获取方法及装置 | |
CN109683147A (zh) | 杂乱脉冲流信号的实时生成方法、装置及电子设备 | |
CN110208760A (zh) | 一种基于时域上采样的雷达回波仿真方法 | |
CN109978888A (zh) | 一种图像分割方法、装置及计算机可读存储介质 | |
DE102011088346A1 (de) | Vorrichtung zur Erfassung akustischer Signale sowie dazugehöriges Verfahren | |
CN109828271A (zh) | Sar成像方法及装置,sar成像系统 | |
Kelly et al. | A fast decimation-in-image back-projection algorithm for SAR | |
Hussain et al. | Auto-regressive spectral gap filling algorithms for photonics-based highly sparse coherent multi-band radars in complex scenarios | |
CN109425859A (zh) | 一种序贯图像成像方法及装置 | |
US6731233B2 (en) | Method of suppressing jammer signals | |
CN102156872A (zh) | 一种基于多光谱数据的物体识别方法和装置 | |
CN107092012A (zh) | 距离补偿方法和装置 | |
Tello et al. | Use of the multiresolution capability of wavelets for ship detection in SAR imagery | |
Sun et al. | Improved incoherent Back Projection imaging based on self-amplitude weighting and multiplicative Tomography weighting | |
Xuezhi et al. | A new approach to distributed passive radar imaging by 2-D NUFFT | |
Chountasis et al. | The whittaker smoother and the moore-penrose inverse in signal reconstruction | |
Frery et al. | Statistical properties of SAR data and their consequences | |
CN113433527B (zh) | 基于mimo雷达的静杂波抑制方法、装置、存储介质及设备 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190531 |
|
RJ01 | Rejection of invention patent application after publication |