CN104865571A - 多通道多子带滑动聚束模式sar成像方法 - Google Patents
多通道多子带滑动聚束模式sar成像方法 Download PDFInfo
- Publication number
- CN104865571A CN104865571A CN201510155399.5A CN201510155399A CN104865571A CN 104865571 A CN104865571 A CN 104865571A CN 201510155399 A CN201510155399 A CN 201510155399A CN 104865571 A CN104865571 A CN 104865571A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- doppler
- frequency
- mfrac
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 26
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 14
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 14
- 230000006835 compression Effects 0.000 claims abstract description 11
- 238000007906 compression Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000002592 echocardiography Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 10
- 238000013507 mapping Methods 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9052—Spotlight mode
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9047—Doppler beam sharpening mode
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多通道多子带滑动聚束模式SAR成像方法,主要解决步进频率多通道合成孔径雷达在滑动聚束模式下的方位重建和频带合成难度大的问题。其实现过程是:1.对多通道天线进行子孔径划分,并对每个子孔径接收的回波依次进行等效相位中心的误差补偿和距离向脉冲压缩;2.对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱移到多普勒基带;3.计算出空域滤波的权重系数;4.利用空域滤波权重系数和基带信号进行多普勒解模糊;5.对多普勒解模糊后信号进行频带合成,完成二维频谱的重建并成像。本发明避免了高分辨率宽测绘数据的巨大计算量,提高了宽测绘带上的目标高分辨率,可用于机载或星载平台SAR多通道多子带成像。
Description
技术领域
本发明属于雷达信号处理领域,涉及步进频合成孔径雷达系统中的频带合成技术,特别涉及多通道多子带滑动聚束SAR成像,可用于机载或星载平台SAR多通道多子带成像处理。
背景技术
合成孔径雷达在对地观测方面的应用越来越广泛,为了缩短对于地面的再访时间,要求观测带的宽度要大。高分辨也是对地进行精确观测的要求,为了获得距离高分辨率的SAR图像,在距离向通过发射大带宽的信号,并对回波进行匹配滤波而获得高分辨距离像。实现大带宽信号的一种比较经济的方式是通过发射步进频率SF的波形信号,利用步进信号的频带合成技术,得到大带宽信号。这样可以降低对雷达接收机带宽的要求,同时降低AD采样频率,便于系统实现。步进信号的频带合成技术在SAR上面都得到了广泛应用。
传统的SAR工作模式如滑动聚束是通过波束控制,其波束指向随慢时间发生变化,从而提高方位分辨率。但是,波束方位变化导致其方位信号带宽急剧增大,通常比采样频率PRF要大很多,因此方位信号通常是混叠的。这使得采用常规成像方法,距离多普勒算法,调频变标算法,距离徙动算法等成像方法不能得到无混叠的图像。一种解混叠的方法就是在系统设计时要提高采样频率PRF,然而大的采样频率PRF将导致方位数据量极大地增加,而且也容易带来距离模糊、数据下传受限等问题。
为了获得宽测绘带,目前国内外开展多发多收MIMO-SAR技术研究,MIMO-SAR很难同时达到二维高分辨率。为了同时达到SAR成像结果的高分辨率和宽测绘带两个要求,需要将滑动聚束模式与MIMO-SAR技术相结合,然而由于天线波束的旋转将导致方位频带的增加,很难实现方位重建和频带合成,从而不能有效实现高分辨率宽测绘带的SAR成像。
发明内容
本发明的目的在于针对上述MIMO-SAR在滑动聚束模式下的方位重建和频带合成难度大的不足,提出一种多通道多子带滑动聚束模式SAR成像方法,获得高分辨率宽测绘带的SAR成像结果。
实现本发明目的的技术关键是:结合MIMO-SAR系统的多通道和SF波形信号多子带的优点,通过方位多通道同时发射一组步进频率SF波形信号,参考图1,成像步骤包括如下:
(1)对步进频率多通道合成孔径雷达SF-MIMO-SAR的天线进行子孔径划分,并对子孔径接收的回波依次进行多通道天线等效相位中心的误差补偿和距离向脉冲压缩;
(2)对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱搬移到多普勒基带;
(3)对以上得到的多普勒基带进行方位快速傅里叶FFT处理,即将多普勒基带的第i子孔径信号转化到多普勒频域,并为了减少空域滤波器中的权重系数的计算量,继续对以上方位FFT处理后的结果按如下补偿函数H3,n(fa,m)进行补偿:
式中,X′nm是第n和m两个通道的等效位置,v是雷达运动速度,fa是方位多普勒频率,fc是距离频率,ti为第i子孔径的方位中心时间;
(4)对补偿后的基带信号数据进行多普勒解模糊操作,计算空域滤波的权重系数矩阵其中:wn(m)是第n个发射通道,第m个接收通道的回波的权重系数;
(5)利用权重系数与补偿后的基带信号相乘,即进行空域滤波处理,得到多普勒解模糊后的基带信号多普勒谱S2,n(t,fd):
式中:Rref是雷达作用的参考距离,wr(·)为距离向的窗函数,wazi(·)为方位向的窗函数,fdc是多普勒中心频率,fc是距离频率,ti为第i子孔径的方位中心时间,fd为方位多普勒频率。
(6)对多普勒解模糊后相邻基带信号进行频带合成,依次在距离频谱上进行小数部分的频移,并将小数部分的频移后相邻基带信号转换到距离频域,进行相干叠加,完成二维频谱重建;
(7)对重建后的二维频谱进行成像。
本发明与现有技术相比具有如下优点:
本发明由于结合MIMO-SAR和步进频率信号各自的优点,利用划分子孔径进行多普勒解模糊处理,利用有效的二维重建算法,避免了高分辨率宽测绘带大数据的巨大计算量,从而实现了高分辨率宽测绘带的成像目标,同时由于本发明采样一种改进的频域带宽合成方法,避免了距离向的上采样,可有效地进行子带信号合成。
附图说明
图1是本发明的实现流程图;
图2是本发明的多通道滑动聚束成像收发阵列的等效相位中心示意图;
图3是本发明的子孔径划分后时频图;
图4是本发明第一子带空域滤波前后的点目标二维频谱图;
图5是本发明的整个频带的信号二维频谱图和点目标距离图像;
图6是本发明在场景坐标(-1km,616km),(0km,617km)和(1km,618km)三点的等高线图;
图7是本发明机载滑动聚束SAR数据的结果图。
具体实施方式
参照图1,本发明对SAR数据处理成像步骤如下:
步骤1,对步进频率多通道合成孔径雷达SF-MIMO-SAR的天线进行子孔径划分;并对每个子孔径接收的回波依次进行多通道天线等效相位中心的误差补偿和距离向脉冲压缩;
1a)参照图2,设发射阵列的阵元数和接收阵列的阵元数均为N,根据收发双方的不同位置Xn和Xm,计算出两者的等效相位中心位置:X′nm=(Xn+Xm)/2;
其中,Xn为第n个发射阵元的位置,n=1,2,...,N;Xm为第m个接收阵元的位置,m=1,2,...,N,N是大于零的整数;
1b)对收发双方的不同位置Xn和Xm两者的等效相位中心X′nm进行相位误差补偿,补偿函数H1,n(m)是:
式中,m为第m个接收阵元的编号,n为第n个接收阵元的编号,Rn为第n个发射阵元与地面目标的距离,n=1,2,...,N;Rm为第m个发射阵元与地面目标的距离,m=1,2,...,N,收发双方两者的等效相位中心与地面目标的距离Rnm,λ是雷达载频的波长;
1c)对相位误差补偿后等效相位中心进行距离压缩处理,压缩后的目标信号s1,n(t,ta,m)公式如下:
其中,为第n个发射阵元和第m个接收阵元的等效相位中心X′nm与某个场景目标X的斜距,wr(·)为距离向的窗函数,wazi(·)为方位向的窗函数,t为快时间,ta慢时间。
步骤2,对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱搬移到多普勒基带。
2a)参考图3,计算每个子孔径的带宽:
Bsub≈Bins+krotTsub+fdc(fc+NΔfc/2,ti)-fdc(fc-NΔfc/2,ti)
式中,Bins为子孔径瞬时带宽,krot=2v2fc/cRrot为多普勒中心波束的变化率,Tsub为子孔径录取数据时间,fdc多普勒中心频率,ti为第i子孔径的方位中心时间,fc雷达载波频率,N为天线阵元数,Δfc是相邻基带信号频谱差值;
2b)设定频移函数为:H2,n(ta,m)=exp(-j2πtafdc(fc,ti)),式中,fdc是多普勒中心频率,fc是距离频率,ta是方位时间,ti为第i子孔径的方位中心时间,m是第m个接收通道的编号;
2c)将子孔径带宽Bsub与多普勒频移函数H2,n(ta,m)进行相乘处理,完成多普勒频移操作,即使多普勒谱搬移到多普勒基带。
步骤3,对步骤二中得到的多普勒基带进行方位快速傅里叶FFT处理,即将多普勒基带的第i子孔径信号转化到多普勒频域,并为了减少空域滤波器中的权重系数的计算量,继续对以上方位FFT处理后的结果按如下补偿函数H3,n(fa,m)进行补偿:
式中,X′nm是第n和m两个通道的等效位置,v是雷达运动速度,fa是方位多普勒频率,fc是距离频率,ti为第i子孔径的方位中心时间。
步骤4,对补偿后的基带信号数据进行多普勒解模糊操作。
4a)计算空域滤波的权重系数矩阵:其中:wn(m)是第n个发射通道,第m个接收通道的回波的权重系数;
4b)利用权重系数与补偿后的基带信号相乘,即进行空域滤波处理,得到多普勒解模糊后的基带信号多普勒谱S2,n(t,fd):
式中:Rref是雷达作用的参考距离,wr(·)为距离向的窗函数,wazi(·)为方位向的窗函数,fdc是多普勒中心频率,fc是距离频率,ti为第i子孔径的方位中心时间,fd为方位多普勒频率,v是载机的飞行速度,c为光速,X为目标的横向向位置。
步骤5,完成二维频谱重建。
5a)对多普勒解模糊后相邻基带信号进行频带合成,在距离频谱上进行小数部分的频移,将多普勒解模糊后的结果乘以如下频移函数H4,n(t),完成小数部分的频移,其中频移函数公式如下:
H4,n(t)=exp(j2πt(Δfc,n-KnΔfr))
式中,Δfr是一个频点单元之间的频率间隔,Δfc,n是相邻基带信号频谱差值,Kn为Δfc,n/Δfr的整数部分,此处频移的数据远小于一个频点单元的值。
5b)将小数部分的频移后相邻基带信号转换到距离频域,进行相干叠加,完成二维频谱重建。
步骤6,对重建后的二维频谱进行成像。
对重建后的二维频谱利用常规的合成孔径雷达SAR成像算法,例如:距离多普勒成像算法和距离徙动成像算法,对二维频谱进行成像处理,实现成像区域的目标完全聚焦,完成高分辨率宽测绘带的成像要求。
以下通过仿真和实测数据来进一步验证本发明的有效性。
一、点目标仿真数据
在2千米×2千米的仿真场景中,在方位向和距离向放置3×3矩阵的9个点目标,仿真的基本参数如表一所示:
表一
二.仿真内容和结果分析
仿真1,用本发明的方法对子带信号进行多普勒解模糊处理,结果如图4,其中:图4(a)是第一子带信号在空域滤波前的二维频域图像,横轴表示方位多普勒,单位是Hz,纵轴表示距离频率,单位是GHz;图4(b)是第一子带信号在空域滤波后的二维频域图像,横轴表示方位多普勒,单位是Hz,纵轴表示距离频率,单位是GHz。
从图4的两幅图中可以看出本发明成功地解决了多普勒模糊的问题,完成多普勒频谱重建。
仿真2,用本发明的方法对两个子带信号在距离频率上进行频带合成,结果如图5,其中:图5(a)是整个频带信号的二维频域图,横轴表示方位多普勒,单位是Hz,纵轴表示距离频率,单位是GHz;图5(b)是点目标在频带合成之前后的距离分辨率切面图图像,横轴表示分辨率的相关距离,单位是米,纵轴表示幅度,单位是分贝。
从图5可以明显看出本发明频带合成后的距离分辨率有显著提高。
仿真3,用本发明的方法对点目标进行成像,结果如图6,其中:图6(a)是场景坐标(-1km,616km)的等高线图,横轴表示方位向,单位是方位分辨单元,纵轴表示距离向,单位是距离分辨单元;图6(b)是场景坐标(0km,617km)的等高线图,横轴表示方位向,单位是方位分辨单元,纵轴表示距离向,单位是距离分辨单元;图6(c)是场景坐标(1km,618km)的等高线图,横轴表示方位向,单位是方位分辨单元,纵轴表示距离向,单位是距离分辨单元。
从图6可以明显看出点目标在距离向和方位向聚焦很好,该仿真验证了本发明的有效性。
三、实测数据
由于国内暂时还没有多通道多子带滑动聚束模式的SAR实测数据,所以本实验利用现有的单通道机载滑动聚束SAR数据来构造数据。操作步骤如下:
首先,将距离时域数据转换到距离频域上,在距离向上将整个频带数据划分为两个子带数据,分别标为第一子带和第二子带;
然后,将第二子带的数据转换到多普勒域,乘上一个时延函数去构成一个不同于第一子带等效相位中心的数据,从而获得两个子带信号的滑动聚束模式的SAR数据。系统参数如表二所示:
表二
四、实测数据处理结果分析
用本发明的方法对实测数据进行处理,结果如图7,其中,图7(a)为用本发明的多普勒解模糊后第一子带数据聚焦的图像,横轴表示方位向,纵轴表示距离向;图7(b)为频带合成后的数据聚焦图像,横轴表示方位向,纵轴表示距离向;图7(c)为频带合成前后图像中一个点目标的距离剖面对比图,横轴表示幅度,单位是分贝;纵轴表示距离分辨率,单元是米。
从图7中可以看出用本发明的方法处理实测数据,可以获得更高的分辨率,因此,实测数据验证了本发明的有效性。
Claims (4)
1.一种多通道多子带滑动聚束模式SAR成像,包括如下步骤:
(1)对步进频率多通道合成孔径雷达SF-MIMO-SAR的天线进行子孔径划分,并对子孔径接收的回波依次进行多通道天线等效相位中心的误差补偿和距离向脉冲压缩;
(2)对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱搬移到多普勒基带;
(3)对以上得到的多普勒基带进行方位快速傅里叶FFT处理,即将多普勒基带的第i子孔径信号转化到多普勒频域,并为了减少空域滤波器中的权重系数的计算量,继续对以上方位FFT处理后的结果按如下补偿函数H3,n(fa,m)进行补偿:
式中,X′nm是第n和m两个通道的等效位置,υ是雷达运动速度,fa是方位多普勒频率,fc是距离频率,ti为第i子孔径的方位中心时间;
(4)对补偿后的基带信号数据进行多普勒解模糊操作,计算空域滤波的权重系数矩阵其中:wn(m)是第n个发射通道,第m个接收通道的回波的权重系数;
(5)利用权重系数与补偿后的基带信号相乘,即进行空域滤波处理,得到多普勒解模糊后的基带信号多普勒谱S2,n(t,fd):
式中:Rref是雷达作用的参考距离,wr(·)为距离向的窗函数,wazi(·)为方位向的窗函数,fdc是多普勒中心频率,fc是距离频率,ti为第i子孔径的方位中心时间,fd为方位多普勒频率;
(6)对多普勒解模糊后相邻基带信号进行频带合成,依次在距离频谱上进行小数部分的频移,并将小数部分的频移后相邻基带信号转换到距离频域,进行相干叠加,完成二维频谱重建;
(7)对重建后的二维频谱进行成像。
2.根据权利要求1中所述的多通道多子带滑动聚束模式SAR成像,其中所述步骤(1)中对子孔径接收的回波进行多通道天线等效相位中心的误差补偿和距离压缩,按如下步骤进行:
(1a)设发射阵列的阵元数和接收阵列的阵元数均为N,根据收发双方的不同位置Xn和Xm,计算出两者的等效相位中心位置:X′nm=(Xn+Xm)/2;
其中,Xn为第n个发射阵元的位置,n=1,2,...,N;Xm为第m个接收阵元的位置,m=1,2,...,N,N是大于零的整数;
(1b)对收发双方的不同位置Xn和Xm两者的等效相位中心X′nm进行相位误差补偿,补偿函数H1,n(m)是:
式中,m为第m个接收阵元的编号,n为第n个接收阵元的编号,Rn为第n个发射阵元与地面目标的距离,n=1,2,...,N;Rm为第m个发射阵元与地面目标的距离,m=1,2,...,N,收发双方两者的等效相位中心与地面目标的距离Rnm,λ是雷达载频的波长;
(1c)对相位误差补偿后等效相位中心进行距离压缩处理,压缩后的目标信号s1,n(t,ta,m)公式如下:
其中, 为第n个发射阵元和第m个接收阵元的等效相位中心X′nm与某个场景目标X的斜距,wr(·)为距离向的窗函数,wazi(·)为方位向的窗函数,t为快时间,ta慢时间。
3.根据权利要求1中所述的多通道多子带滑动聚束模式SAR成像,其中步骤(2)所述的对距离向脉冲压缩后的信号进行多普勒频移操作,按如下步骤进行:
(2a)计算子孔径带宽:
Bsub≈Bins+krotTsub+fdc(fc+NΔfc/2,ti)-fdc(fc-NΔfc/2,ti)
式中,Bins为子孔径瞬时带宽,krot=2υ2fc/cRrot为多普勒中心波束的变化率,Tsub为子孔径录取数据时间,fdc多普勒中心频率,Rrot为旋转中心距离,ti为第i子孔径的方位中心时间,υ雷达运动速度,fc雷达载波频率,N为天线阵元数,Δfc是相邻基带信号频谱差值;
(2b)设定频移函数为:H2,n(ta,m)=exp(-j2πtafdc(fc,ti))
式中,fdc是多普勒中心频率,fc是距离频率,ta是方位时间,ti为第i子孔径的方位中心时间,m是第m个接收通道的编号;
(2c)将子孔径带宽Bsub与多普勒频移函数H2,n(ta,m)进行相乘处理,完成多普勒频移操作,即使多普勒谱搬移到多普勒基带。
4.根据权利要求1中所述的多通道多子带滑动聚束模式SAR成像,其中所述步骤(6)中在距离频谱上进行小数部分的频移,是将多普勒解模糊后的结果乘以如下频移函数H4,n(t),完成小数部分的频移,其中频移函数公式如下:
H4,n(t)=exp(j2πt(Δfc,n-KnΔfr))
式中,Δfr是一个频点单元之间的频率间隔,Δfc,n是相邻基带信号频谱差值,Kn为Δfc,n/Δfr的整数部分,此处频移的数据远小于一个频点单元的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510155399.5A CN104865571A (zh) | 2015-04-02 | 2015-04-02 | 多通道多子带滑动聚束模式sar成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510155399.5A CN104865571A (zh) | 2015-04-02 | 2015-04-02 | 多通道多子带滑动聚束模式sar成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104865571A true CN104865571A (zh) | 2015-08-26 |
Family
ID=53911532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510155399.5A Pending CN104865571A (zh) | 2015-04-02 | 2015-04-02 | 多通道多子带滑动聚束模式sar成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104865571A (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105242252A (zh) * | 2015-09-30 | 2016-01-13 | 西安电子科技大学 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN106443672A (zh) * | 2016-08-30 | 2017-02-22 | 西安电子科技大学 | 一种方位多通道sar信号自适应重构方法 |
CN106597442A (zh) * | 2016-12-21 | 2017-04-26 | 中国航空工业集团公司雷华电子技术研究所 | 一种方位多通道脉内聚束sar成像方法 |
CN106842196A (zh) * | 2015-12-03 | 2017-06-13 | 中国航空工业集团公司雷华电子技术研究所 | 一种用于sar子孔径成像的方位模糊抑制方法 |
CN106908789A (zh) * | 2017-02-28 | 2017-06-30 | 西安电子科技大学 | 一种基于空间域支撑区融合的sar成像方法 |
CN107229036A (zh) * | 2017-05-27 | 2017-10-03 | 西安电子科技大学 | 基于信号处理的多通道阵列雷达幅相误差在线检测方法 |
CN107462888A (zh) * | 2017-07-19 | 2017-12-12 | 中国人民解放军国防科学技术大学 | 滑动聚束sar的步进扫描影响参数计算方法 |
CN107741586A (zh) * | 2017-09-29 | 2018-02-27 | 王辉 | 基于DBF‑TOPS加权的星载Ka InSAR信号处理方法 |
CN107942327A (zh) * | 2017-10-20 | 2018-04-20 | 西安电子科技大学 | 基于脉冲相位编码的单通道hrws‑sar成像方法 |
CN108469608A (zh) * | 2018-03-21 | 2018-08-31 | 电子科技大学 | 一种运动平台雷达多普勒质心精确估计方法 |
CN108646242A (zh) * | 2018-05-10 | 2018-10-12 | 北京航空航天大学 | 一种针对复杂目标的多子波段雷达数据融合成像方法 |
CN109960843A (zh) * | 2017-12-26 | 2019-07-02 | 中国船舶重工集团公司七五〇试验场 | 一种基于正交原理的多普勒频移数值仿真方法 |
CN109975804A (zh) * | 2019-03-04 | 2019-07-05 | 广东工业大学 | 一种多平台星座sar融合相参成像方法 |
CN110045374A (zh) * | 2019-05-28 | 2019-07-23 | 辽宁工程技术大学 | 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法 |
RU199139U1 (ru) * | 2020-01-27 | 2020-08-19 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны им. Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Радиоприемное устройство импульсно-доплеровской РЛС с многоканальной весовой обработкой |
CN111929682A (zh) * | 2020-07-22 | 2020-11-13 | 西安瑞得空间信息技术有限公司 | 基于距离向分块的机载dbf-tops sar成像方法 |
CN111965644A (zh) * | 2020-07-15 | 2020-11-20 | 中国科学院国家空间科学中心 | 一种行波天线合成孔径雷达成像方法 |
CN112180368A (zh) * | 2020-09-10 | 2021-01-05 | 中国科学院空天信息创新研究院 | 一种数据处理方法、设备、系统及存储介质 |
CN115128608A (zh) * | 2022-08-31 | 2022-09-30 | 南京隼眼电子科技有限公司 | 车载多普勒分集mimo合成孔径雷达时域成像方法 |
CN115407279A (zh) * | 2022-09-01 | 2022-11-29 | 中国人民解放军海军工程大学 | 一种脉冲压缩雷达距离多普勒耦合误差补偿方法 |
CN115436942A (zh) * | 2022-09-21 | 2022-12-06 | 南京航空航天大学 | 一种提高mimo-sar俯仰向有效等效相位中心数目的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776755A (zh) * | 2010-02-05 | 2010-07-14 | 西安电子科技大学 | 高速平台超高分辨率sar数据预处理方法 |
CN101907704A (zh) * | 2010-06-11 | 2010-12-08 | 西安电子科技大学 | 多模式合成孔径雷达仿真成像评估方法 |
CN103954938A (zh) * | 2014-01-26 | 2014-07-30 | 西安电子科技大学 | 一种sar回波信号的多子带接收合成方法 |
-
2015
- 2015-04-02 CN CN201510155399.5A patent/CN104865571A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776755A (zh) * | 2010-02-05 | 2010-07-14 | 西安电子科技大学 | 高速平台超高分辨率sar数据预处理方法 |
CN101907704A (zh) * | 2010-06-11 | 2010-12-08 | 西安电子科技大学 | 多模式合成孔径雷达仿真成像评估方法 |
CN103954938A (zh) * | 2014-01-26 | 2014-07-30 | 西安电子科技大学 | 一种sar回波信号的多子带接收合成方法 |
Non-Patent Citations (2)
Title |
---|
JIAJIA ZHANG等: ""An Efficient Signal Reconstruction Algorithm for Stepped Frequency MIMO-SAR in the Spotlight and Sliding Spotlight Modes"", 《INTERNATIONAL JOURNAL OF ANTENNAS AND PROPAGATION》 * |
NICOLAS GEBERT等: ""Multichannel Azimuth Processing in ScanSAR and TOPS Mode Operation"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105242252B (zh) * | 2015-09-30 | 2017-09-12 | 西安电子科技大学 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
CN105242252A (zh) * | 2015-09-30 | 2016-01-13 | 西安电子科技大学 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
CN106842196B (zh) * | 2015-12-03 | 2019-12-17 | 中国航空工业集团公司雷华电子技术研究所 | 一种用于sar子孔径成像的方位模糊抑制方法 |
CN106842196A (zh) * | 2015-12-03 | 2017-06-13 | 中国航空工业集团公司雷华电子技术研究所 | 一种用于sar子孔径成像的方位模糊抑制方法 |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN106443672A (zh) * | 2016-08-30 | 2017-02-22 | 西安电子科技大学 | 一种方位多通道sar信号自适应重构方法 |
CN106443672B (zh) * | 2016-08-30 | 2019-03-29 | 西安电子科技大学 | 一种方位多通道sar信号自适应重构方法 |
CN106597442A (zh) * | 2016-12-21 | 2017-04-26 | 中国航空工业集团公司雷华电子技术研究所 | 一种方位多通道脉内聚束sar成像方法 |
CN106597442B (zh) * | 2016-12-21 | 2018-11-16 | 中国航空工业集团公司雷华电子技术研究所 | 一种方位多通道脉内聚束sar成像方法 |
CN106908789A (zh) * | 2017-02-28 | 2017-06-30 | 西安电子科技大学 | 一种基于空间域支撑区融合的sar成像方法 |
CN106908789B (zh) * | 2017-02-28 | 2019-08-02 | 西安电子科技大学 | 一种基于空间域支撑区融合的sar成像方法 |
CN107229036A (zh) * | 2017-05-27 | 2017-10-03 | 西安电子科技大学 | 基于信号处理的多通道阵列雷达幅相误差在线检测方法 |
CN107229036B (zh) * | 2017-05-27 | 2020-07-03 | 西安电子科技大学 | 基于信号处理的多通道阵列雷达幅相误差在线检测方法 |
CN107462888B (zh) * | 2017-07-19 | 2019-06-25 | 中国人民解放军国防科学技术大学 | 滑动聚束sar的步进扫描影响参数计算方法 |
CN107462888A (zh) * | 2017-07-19 | 2017-12-12 | 中国人民解放军国防科学技术大学 | 滑动聚束sar的步进扫描影响参数计算方法 |
CN107741586A (zh) * | 2017-09-29 | 2018-02-27 | 王辉 | 基于DBF‑TOPS加权的星载Ka InSAR信号处理方法 |
CN107942327A (zh) * | 2017-10-20 | 2018-04-20 | 西安电子科技大学 | 基于脉冲相位编码的单通道hrws‑sar成像方法 |
CN109960843A (zh) * | 2017-12-26 | 2019-07-02 | 中国船舶重工集团公司七五〇试验场 | 一种基于正交原理的多普勒频移数值仿真方法 |
CN109960843B (zh) * | 2017-12-26 | 2023-04-18 | 中国船舶重工集团公司七五〇试验场 | 一种基于正交原理的多普勒频移数值仿真方法 |
CN108469608A (zh) * | 2018-03-21 | 2018-08-31 | 电子科技大学 | 一种运动平台雷达多普勒质心精确估计方法 |
CN108646242A (zh) * | 2018-05-10 | 2018-10-12 | 北京航空航天大学 | 一种针对复杂目标的多子波段雷达数据融合成像方法 |
CN108646242B (zh) * | 2018-05-10 | 2022-02-08 | 北京航空航天大学 | 一种针对复杂目标的多子波段雷达数据融合成像方法 |
CN109975804A (zh) * | 2019-03-04 | 2019-07-05 | 广东工业大学 | 一种多平台星座sar融合相参成像方法 |
CN110045374A (zh) * | 2019-05-28 | 2019-07-23 | 辽宁工程技术大学 | 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法 |
CN110045374B (zh) * | 2019-05-28 | 2023-07-25 | 辽宁工程技术大学 | 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法 |
RU199139U1 (ru) * | 2020-01-27 | 2020-08-19 | Федеральное государственное казённое военное образовательное учреждение высшего образования "Военная академия воздушно-космической обороны им. Маршала Советского Союза Г.К. Жукова" Министерства обороны Российской Федерации | Радиоприемное устройство импульсно-доплеровской РЛС с многоканальной весовой обработкой |
CN111965644A (zh) * | 2020-07-15 | 2020-11-20 | 中国科学院国家空间科学中心 | 一种行波天线合成孔径雷达成像方法 |
CN111929682A (zh) * | 2020-07-22 | 2020-11-13 | 西安瑞得空间信息技术有限公司 | 基于距离向分块的机载dbf-tops sar成像方法 |
CN111929682B (zh) * | 2020-07-22 | 2022-07-19 | 西安瑞得空间信息技术有限公司 | 基于距离向分块的机载dbf-tops sar成像方法 |
CN112180368B (zh) * | 2020-09-10 | 2022-07-15 | 中国科学院空天信息创新研究院 | 多通道滑动聚束sar的数据处理方法、设备、系统及介质 |
CN112180368A (zh) * | 2020-09-10 | 2021-01-05 | 中国科学院空天信息创新研究院 | 一种数据处理方法、设备、系统及存储介质 |
CN115128608B (zh) * | 2022-08-31 | 2022-11-18 | 南京隼眼电子科技有限公司 | 车载多普勒分集mimo合成孔径雷达时域成像方法 |
CN115128608A (zh) * | 2022-08-31 | 2022-09-30 | 南京隼眼电子科技有限公司 | 车载多普勒分集mimo合成孔径雷达时域成像方法 |
CN115407279A (zh) * | 2022-09-01 | 2022-11-29 | 中国人民解放军海军工程大学 | 一种脉冲压缩雷达距离多普勒耦合误差补偿方法 |
CN115407279B (zh) * | 2022-09-01 | 2024-04-19 | 中国人民解放军海军工程大学 | 一种脉冲压缩雷达距离多普勒耦合误差补偿方法 |
CN115436942A (zh) * | 2022-09-21 | 2022-12-06 | 南京航空航天大学 | 一种提高mimo-sar俯仰向有效等效相位中心数目的方法 |
CN115436942B (zh) * | 2022-09-21 | 2023-12-01 | 南京航空航天大学 | 一种提高mimo-sar俯仰向有效等效相位中心数目的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104865571A (zh) | 多通道多子带滑动聚束模式sar成像方法 | |
CN103728618B (zh) | 一种高分辨率、宽测绘带的星载sar体制实现方法 | |
CN107741586B (zh) | 基于DBF-TOPS加权的星载Ka InSAR信号处理方法 | |
CN102680974B (zh) | 一种星载滑动聚束合成孔径雷达的信号处理方法 | |
CN104714231B (zh) | 一种基于完全互补序列和相位补偿的mimo sar成像方法 | |
CN104166141A (zh) | 基于子带合成的多输入多输出合成孔径雷达系统设计方法 | |
CN103235305B (zh) | 星载超高分辨率滑动聚束sar成像方法 | |
CN103885061B (zh) | 基于改进相位梯度自聚焦的双基sar运动补偿方法 | |
CN106932778B (zh) | 方位多通道fmcw sar滑动聚束成像方法 | |
CN102226841A (zh) | 基于高阶多项式距离方程的同步轨道sar成像方法 | |
CN104678393B (zh) | 斜视滑动聚束sar的子孔径波数域成像方法 | |
CN106597442A (zh) | 一种方位多通道脉内聚束sar成像方法 | |
CN109613507B (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
CN102608603A (zh) | 一种基于完全互补序列的多通道合成孔径雷达成像方法 | |
CN104898107A (zh) | 一种多发多收合成孔径激光雷达信号处理方法 | |
CN104133215A (zh) | 基于距离徙动微调和子带分割的同步轨道雷达成像方法 | |
CN108535726A (zh) | 基于幂权Fourier变换的ISAR成像方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN105974409B (zh) | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 | |
CN105005045A (zh) | 基于信号预处理的高速目标isar步进频信号合成方法 | |
Liu et al. | Novel range profile synthesis algorithm for linearly stepped-frequency modulated inversed synthetic aperture radar imaging of remote manoeuvring target | |
Yang et al. | A subaperture imaging algorithm to highly squinted TOPS SAR based on SPECAN and deramping | |
RU2309425C2 (ru) | Способ калибровки радиопеленгатора-дальномера | |
CN109633641B (zh) | 一种太赫兹频段的旋翼叶片逆合成孔径雷达成像算法 | |
Wang | Bandwidth synthesis for stepped chirp signal: A multichannel sampling prospective |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150826 |
|
WD01 | Invention patent application deemed withdrawn after publication |