CN101833095B - 一种基于空域展开的星机联合sar二维频域成像方法 - Google Patents
一种基于空域展开的星机联合sar二维频域成像方法 Download PDFInfo
- Publication number
- CN101833095B CN101833095B CN2010190870579A CN201019087057A CN101833095B CN 101833095 B CN101833095 B CN 101833095B CN 2010190870579 A CN2010190870579 A CN 2010190870579A CN 201019087057 A CN201019087057 A CN 201019087057A CN 101833095 B CN101833095 B CN 101833095B
- Authority
- CN
- China
- Prior art keywords
- synthetic aperture
- aperture radar
- data
- machine combined
- radar
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于空域展开的星机联合双基地合成孔径雷达二维频域成像方法,它是利用距离历史在参考点的空域展开和参考点距离历史的二阶时域近似,将相位表示成一个关于慢时间的二次多项式,并利用驻定相位法得到信号的二维频谱。通过方位和距离压缩得到良好的聚焦效果。最后采用利用散射点位置的线性关系,将距离和方位的坐标从耦合项中剥离出来,得到正确的坐标。从而克服了复杂的二维频谱计算和目标空变补偿的过程,并实现了较大场景下星机联合双基地合成孔径雷达成像。本发明可以应用于合成孔径雷达成像,地球遥感等领域。
Description
技术领域:
本发明属于雷达技术领域,它特别涉及合成孔径雷达(SAR)成像技术中星机联合模式双基地合成孔径雷达成像方法。
背景技术:
双基地合成孔径雷达(Bistatic SAR)是将接收机和发射机分别安装在不同的运动平台上的一种新型合成孔径雷达***。星机联合双基地合成孔径雷达(Spaceborne-Airborne Bistatic SAR)作为双基地合成孔径雷达的一种特殊模式,采用卫星发射、单(多)个飞机接收的工作模式,在继承了双基合成孔径雷达优点的同时,还具有以下独特优点:
1.独特的“远发近收”模式,既充分发挥了卫星站得高、看得远、覆盖面广等优势,又保持了很高的图像信噪比;
2.降低对卫星功率、数据传输容量、处理能力及成本等方面需求;
3.根据客户需求制定观测方案,实施比分布式星载SAR***更灵活的数据采集方式,降低数据获取成本;
4.发挥飞机机动灵活的特点,构建不同于传统条带、聚束及扫描模式的新型工作模式,便于高分辨率和大测绘带SAR***设计实现。
因此星机联合双基地SAR作为一项新概念的空间对地观测或侦察手段,无论在民用还是军事应用领域都有着更为广泛的发展空间,成了最近双基地合成孔径雷达领域的研究热点。
一般地,双基地合成孔径雷达***可以分成两类:移不变双基地合成孔径雷达***和移变双基地合成孔径雷达***。星机联合双基地合成孔径雷达作为典型的移变双基地合成孔径雷达***,发射/接收***的几何关系随时间变化,因此基于线性时不变假设的合成孔径雷达成像方法,如,距离-多普勒方法、波数域方法,不能满足此类双基地合成孔径雷达成像的要求;另外,即使实现了移变双基地合成孔径雷达成像点目标成像,移变双基地合成孔径雷达成像距离向和方位向坐标系的非正交性仍不可避免的导致移变双基地合成孔径雷达图像失真。目前针对移变双基地合成孔径雷达成像方法,特别是星机联合模式下的成像方法比较有限,在有效的时域方面除了后向投影(BP)方法和变尺度逆傅立叶变换方法外很少有其他方法。而且变尺度逆傅立叶变换方法由于没有考虑空域截断误差的 影响,所以并不适应大场景下的星机联合双基地合成孔径雷达成像。见文献“ShiJun;Xiaoling Zhang;Jianyu Yang;“Principle and Methods on Bistatic SARSignal Processing via Time Correlation”,Geoscience and Remote Sensing,IEEE Transactions,Volume 46,pp:3163-3178(2008)”。
频域方法的优势在于运算速度快,但成像精度不高,同时回波二维频谱的公式推导相对复杂,目前以德国学者Otmar Loffled所提出的LBF模型为代表。
Loffled等所提出的LBF模型将相位在发射机和接收机的驻定相位点泰勒展开,将相位表示成两个二次项之和,从而计算得到二维频谱,并引入变尺度傅里叶变换得到良好的聚焦效果。但该方法的频谱推导过程过于冗长和复杂,相位误差较大,且处理过程的有关系数的计算也相当繁琐。见文献“K.Natroshvili,O.Loffeld,H.Nies,A.Medrano Ortiz,and S.Knedlik,“Focusing of general bistaticSAR configuration data with 2-D inverse scaled FFT,”IEEE Trans.Geosci.Remote Sens.,vol.44,no.10,pp.2718-2727,Oct.2006”;以及“R.Wang,O.Loffeld,Q.UI-Ann,H.Nies,A.Medrano-Ortiz,and S.Knedlik,“Analysis andextension of Loffeld’s bistatic formula in spaceborne/airborne configuration,”inProc.EUSAR,Friedrichshafen,Germany,Jun.2008.CD-ROM.”。
刘喆提出的Air-Phase模型将距离历史看成一个整体,在保证相位误差最小的情况下,求得慢时间的展开点,并以此获得了较好的二维频谱解析表达式。但由于时域展开式无法用显式表达散射点之间的空间特性,因此,该方法在成像中需要通过一系列复杂的步骤来完成散射点的空变校正。见文献“Zhe Liu,JianyuYang,Xiaoling Zhang,Yiming Pi.Study on Spaceborne/Airborne HybridBistatic SAR Image Formation in Frequency Domain.IEEE Geoscience andRemote Sensing Letters,2008,5(4):578-582”。
综上所述,对距离历史的时域展开不仅使得频谱的推导变得冗长和复杂,而且无法用显式表达散射点之间的空间特性。而对距离历史采用空域展开则可以避免并克服上述缺点。因此本发明基于距离历史的空域展开,提出了一种新的适用于星机联合SAR的二维频域方法。
发明内容:
为了克服星机联合双基地合成孔径雷达缺乏有效简单的频域成像方法的问题,本发明提供了一种基于空域展开的星机联合双基地合成孔径雷达二维频域成像方法,该方法的特点和思路是:利用距离历史在参考点的空域展开和参考点距离历史的二阶时域近似,将相位表示成一个关于慢时间的二次多项式,并利用驻定相位法得到信号的二维频谱。通过方位和距离压缩得到良好的聚焦效果。最后 采用利用散射点位置的线性关系,将距离和方位的坐标从耦合项中剥离出来,得到正确的坐标。从而克服了复杂的二维频谱计算和目标空变补偿的过程,并实现了较大场景下星机联合双基地合成孔径雷达成像。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、星机联合双基地合成孔径雷达(Spaceborne-Airborne Bistatic SAR)
双基地合成孔径雷达是指雷达发射***和接收***分别安装在不同运动平台上的合成孔径雷达,其中,安装发射***的平台称作发射平台,安装接收***的平台称作接收平台。
星机联合双基地合成孔径雷达是指雷达的发射平台和接收平台分别置于卫星和飞机上的双基地合成孔径雷达的特殊模式。
定义2、星机联合双基地合成孔径雷达***相关参数描述
PRF表示脉冲重复频率;
发射平台雷达的波束指向矢量 发射平台雷达的初始波束指向矢量 接收平台雷达的波束指向矢量 接收平台雷达的初始波束指向矢量 星机联合双基地合成孔径雷达***波束指向矢量 为了记录方便令 表示 在n=0的值。
其他参数:t为快(斜距)时间;n为第n个PRF时刻,也即为慢(方位)时间; 表示载波波数, 表示距离波数,K′=K-K0, 为方位波数;C为光速,f0为载波频率,f为对应于快(斜距)时间的频率,fa为对应于慢(方位)时间的多普勒频率。
定义3、移变模式双基地合成孔径雷达
广义上讲,移变模式双基地合成孔径雷达是指安装发射***和接收***的平台在数据采集过程中相对位置发生变化的合成孔径雷达***。
但实际情况下,发射平台和接收平台的运动轨迹总保持匀速直线运动。因此,本发明中定义“移变模式双基地合成孔径雷达”是指发射***平台和接收***平台的运动轨迹,在数据采集过程中,总保持匀速直线运动的广义移变模式双基地合成孔径雷达。
定义4合成孔径雷达标准距离压缩方法
合成孔径雷达标准距离压缩方法是指利用合成孔径雷达发射参数,主要包括:采用以下公式生成参考信号,并采用匹配滤波技术对合成孔径雷达的距离向信号进行滤波的过程。
其中,f(t)为参考函数,B为雷达发射基带信号的信号带宽,TP为雷达发射信号脉冲宽度,t为快(斜距)时间,取值范围从 到 详见文献“雷达成 像技术”,保铮等编著,电子工业出版社出版。本说明书中将距离压缩后得到的数据称为距离域数据。
定义5合成孔径雷达场景空间
合成孔径雷达场景空间是指现实空间中所有待观测的场景目标点的集合。在不同空间坐标系下有不同的表示,但一旦坐标系确立以后其表示是唯一的。一般情况下为了方便成像取地面坐标系,即距离向-方位向-高度向坐标系。
定义6合成孔径雷达成像空间
合成孔径雷达成像空间是指合成孔径雷达成像方法将场景空间中的散射点投影到的二维平面空间,该空间由合成孔径雷达成像空间中的两个相互正交的坐标基确定,目前典型合成孔径雷达的成像空间包括距离向-方位向投影空间。本发明中用以下数学关系表示成像空间M:
定义7合成孔径雷达成像场景参考点
合成孔径雷达成像场景参考点是指合成孔径雷达成像空间中的某个散射点,作为分析和处理场景中其他散射点的参照。
定义8星机联合双基地合成孔径雷达距离历史的空域展开
星机联合双基地合成孔径雷达距离历史是指任意时刻场景空间中待观测的散射点到卫星平台和飞机平台的距离和。合成孔径雷达距离历史的空域展开是指将距离历史沿场景空间的三个方向做多元泰勒展开(multi-variables Taylor’stheorem)的过程。其中被忽略的二阶及二阶以上的误差称之为空域截断误差。见文献“Shi Jun;Xiaoling Zhang;Jianyu Yang;“Principle and Methods on BistaticSAR Signal Processing via Time Correlation”,Geoscience and Remote Sensing,IEEE Transactions,Volume 46,pp:3163-3178(2008)”
其中,相关符号解释见定义2,需要注意的是这里:
定义9星机联合双基地合成孔径雷达参考点距离历史时域展开
星机联合双基地合成孔径雷达距离历史是指任意时刻场景空间中参考点到卫星平台和飞机平台的距离历史和,记为 星机联合双基地合成孔径雷达距离历史的时域展开是指将距离历史沿慢时间做泰勒展开(Taylor’s theorem)的过程
其中:
A表示初始距离历史;B表示等效速度;D表示等效加速度。其他相关符号见定义2;
定义10星机联合双基地合成孔径雷达距离历史近似表达式
星机联合双基地合成孔径雷达距离历史是指任意时刻场景空间中待观测的散射点到卫星平台和飞机平台的距离历史的近似表示过程。即结合公式(3)和(6)的表达式:
定义11星机联合双基地合成孔径雷达***冲激响应的二维频谱
根据星机联合双基地合成孔径雷达***的特点,定义星机联合双基地合成孔径雷达***冲激响应的二维频谱H(K,Ka)为
其中相关符号见定义2和定义10;
定义12驻定相位原理
驻定相位原理是对具有大时宽带宽积信号的频谱分析方法。具体参见文献“皮亦鸣,杨建宇,合成孔径雷达成像原理,电子科技大学出版社,2007年3月,p22-p24.”
本发明提供的一种基于空域展开的星机联合SAR二维频域成像方法,它包括以下几个步骤:
步骤一、合成孔径雷达成像空间及成像***参数的初始化;
合成孔径雷达成像空间由合成孔径雷达成像空间中的两个相互正交的坐标基确定,定义与发射平台(即卫星)速度方向平行并在地平面内的单位向量作为合成孔径雷达成像空间的第一个坐标基,即 定义在地平面内,并与合成孔径雷达成像空间的第一个坐标基 垂直的单位向量作为合成孔径雷达成像空间的第二个坐标基,记做
初始化成像***参数包括:发射平台速度矢量,即卫星相对地面速度矢量,记做 接收平台速度矢量,即飞机速度矢量记做 发射平台初始位置矢量,记做 接收平台初始位置矢量,记做 雷达发射电磁波的波数,记做K0, 雷达发射基带信号的信号带宽,记做B,雷达发射信号脉冲宽度,记做TP,雷达接收波门持续宽度,记做To,雷达接收***的采样频率,记做fs,雷达***的脉冲重复频率,记做PRF。发射平台初始位置矢量 接收平台初始位置矢量 其中xT yT zT和xR yR zR分别表示发射平台和接收平台的空间三维坐标。采用公式 得到发射平台第n个PRF时刻的位置矢量 采用公式 得到发射平台第n个PRF时刻的位置矢量 n=1,...,Ns表示第n个PRF时刻,Ns为方位向的采样点数。
发射雷达的波束初始指向矢量,记做 接收雷达的波束初始指向矢量 表示雷达的初始波束指向和矢量。发射平台雷达的初始角速度矢量 接收平台雷达的初始角速度矢量 星机联合双基地合成孔径雷达***初始角速度和矢量 并利用公式 分别求得参数:A,B,D。其中:A表示初始距离历史;B表示等效速度;D表示等效加速度; 符号“||”表示求模运算。
根据星机联合双基地合成孔径雷达***方案和星机联合双基地合成孔径雷达观测方案,星机联合双基地合成孔径雷达成像方法需要的初始化成像***参数均为已知。其几何结构如图1所示。
步骤二、星机联合双基地合成孔径雷达原始数据进行距离压缩。
我们将接收机接收到的星机联合双基地合成孔径雷达回波数据记做 采用合成孔径雷达标准距离压缩方法对接收到的星机联合双基地合成孔径雷达回波数据进行压缩,得到距离压缩后的星机联合双基地合成孔径雷达数据,记做 n表示第n个PRF时刻,r为距离自变量, 表示待观测的点目标向量。
步骤三、星机联合双基地合成孔径雷达原始数据进行二维傅里叶变换
步骤四、星机联合双基地合成孔径雷达方位参考函数补偿
步骤五、星机联合双基地合成孔径雷达非线性相位补偿
其中 表示载波波数, 表示距离波数,K′=K-K0, 为方位波数;C为光速,f0为载波频率,f为对应于快(斜距)时间的频率,fa为对应于慢(方位)时间的多普勒频率,A表示初始距离历史;B表示等效速度;D表示等效加速度; 符号“||”表示求模运算。
步骤六、星机联合双基地合成孔径雷达解耦合
由于移变双基地合成孔径雷达成像距离向和方位向坐标系的非正交性将导致移变双基地合成孔径雷达图像失真。因此需要分别将距离相位和方位相位中的耦合项剥离出来,即解耦合。
将步骤五中得到的关于待观测散射点坐标线性相位的数据 沿着行(即距离向)做逆傅里叶变换,得到距离-多普勒域数据 结合步骤一中的相关初始化参数,将数据 中的第ii列的数据分别乘以方位解耦合函数 其中 表示第ii个快时间的取值,ii=1...Nr,Nr为距离采样点;fs为雷达接收***的采样频率。得到方 位解耦合后的数据
将得到的方位解耦合后的数据 沿列做逆傅里叶变换,得到距离频率-方位时间域数据 对数据 中的第jj行数据分别乘以距离解耦合函数 其中 表示第jj个慢时间的取值,jj=1...Ns,Ns为方位采样点,PRF为脉冲重复频率, 为常系数。得到距离解耦后的数据,并对其沿行做逆傅里叶变换,得到最终的成像结果
其中 表示载波波数, 表示距离波数,K′=K-K0, 为方位波数;C为光速,f0为载波频率,f为对应于快(斜距)时间的频率,fa为对应于慢(方位)时间的多普勒频率,A表示初始距离历史;B表示等效速度;D表示等效加速度; 符号“||”表示求模运算。
需要指出的是,由于不同合成孔径雷达发射的信号可能存在互为共轭的关系,本发明步骤中的星机联合双基地合成孔径雷达场景方位向参考函数sref(K′,Ka)的相位中的符号,以及在星机联合双基地合成孔径雷达非线性相位补偿和解耦合中用到的相位Ψ(K)cmp,Ψ*(Ka)cmp以及Θ*(K)cmp的符号应根据实际合成孔径雷达发射的信号的相位函数的符号改变。
本发明的实质和创新点是针对星机联合双基地合成孔径雷达***的空变特性:即其距离徙动曲线随时间和距离同时变化的特点。利用距离历史在参考点的空域展开和参考点距离历史的二阶时域近似,将相位表示成一个关于慢时间的二次多项式,并利用驻定相位法得到信号的二维频谱。通过方位和距离压缩得到良好的聚焦效果。最后采用利用散射点位置的线性关系,将距离和方位的坐标从耦 合项中剥离出来,得到正确的坐标。从而克服了复杂的二维频谱计算和目标空变补偿的过程,并实现了较大场景下星机联合双基地合成孔径雷达成像。本发明可以应用于合成孔径雷达成像,地球遥感等领域。
附图说明:
图1为本发明具体实施方式采用的星机联合双基地合成孔径雷达飞行几何关系图。
其中 和 分别为发射平台和接收平台的在坐标系中的初始位置;其中的xT,yT,zT和xR,yR,zR分别为发射平台和接收平台在坐标系中的三维坐标; 表示场景中的待测点向量; 和 表示场景坐标轴的单位矢量。 和 分别表示发射平台和接收平台雷达波束指向单位向量; 和 分别表示发射机和接收机的速度矢量。
图2是发明具体实施方式采用的星机联合双基地合成孔径雷达***参数表。
图3为本发明所提供方法的流程示意框图。
图4是通过本发明提供的方法得到的9点目标星机联合双基地合成孔径雷达成像结果。
图中9个正方形分布的黑点为布置于地面上的9个正方形分布的散射点的星机联合双基地合成孔径雷达成像结果。这9个点在图中分别用字母A,B,C,D,E,F,G,H,I表示。从图中可以看出,本发明提供的方法可以很好的实现星机联合双基地合成孔径雷达点目标成像处理。
从图中可以看出,本发明提供的方法可以很好的实现星机联合双基地合成孔径雷达面目标成像处理。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLAB7.0上验证正确。具体实施步骤如下:
步骤一、利用计算机仿真产生星机联合双基地合成孔径雷达的仿真数据,仿真所需的***参数如图3所示。并定义合成孔径雷达成像空间及初始化成像系 统参数;
初始化卫星相对地面速度矢量 接收平台速度矢量,即飞机速度矢量 发射平台初始位置矢量 接收平台初始位置矢量 采用公式 得到发射平台第n个PRF时刻的位置矢量 采用公式 得到发射平台的第n个PRF时刻的位置矢量
发射雷达的波束初始指向矢量, 接收雷达的波束初始指向矢量, 星机联合双基地合成孔径雷达***波束初始指向矢量 发射平台雷达的初始角速度矢量 接收平台雷达的初始角速度矢量 星机联合双基地合成孔径雷达***初始角速度矢量 并利用公式 求得参数:A=524110,B=30.6587,D=96.1151,αu=-0.8185,αv=-0.0480, βu=-0.0015,βv=-0.022。
本试验选择的用于初始化星机联合双基地合成孔径雷达成像方法的***参数与表一中提供的参数一致。
步骤二、星机联合双基地合成孔径雷达原始数据进行压缩。
将接收到的合成孔径雷达距离向回波信号,按照脉冲重复周期写成矩阵形式,记做 矩阵行表示PRF时刻,矩阵列表示回波延时时间;采用合成孔径雷达标准距离压缩方法对接收到的星机联合双基地合成孔径雷达原始数据进行压缩,得到距离压缩后的星机联合双基地合成孔径雷达距离域数据,记做
步骤三、星机联合双基地合成孔径雷达原始数据进行二维傅里叶变换
对步骤二得到的星机联合双基地合成孔径雷达距离压缩后的以矩阵形式储存的回波数据 沿着行做Nr=2000点的傅里叶变换,再沿着列做Ns=1200点傅里叶变换,得到星机联合双基地合成孔径雷达回波信号的二维频域数据,记为
步骤四、星机联合双基地合成孔径雷达方位参考函数补偿
利用步骤一中的参数,生成一个Ns×Nr的矩阵,即为方位向参考函数 其中 为Ns点向量, 为Nr点向量。将步骤三得到的星机联合双基地合成孔径雷达回波信号的二维频域数据 与方位向参考函数sref(K′,Ka)相乘,得到方位补偿后的星机联合双基地合成孔径雷达回波信号的二维频域数据
步骤五、星机联合双基地合成孔径雷达非线性相位补偿
将步骤四中得到的方位补偿后的星机联合双基地合成孔径雷达回波信号的二维频域数据 沿着数据的列做Ns=1200点逆傅里叶变换(IFFT),得到星机联合双基地合成孔径雷达回波信号的距离频率-方位时间域数据:
●结合步骤一中的相关初始化参数,并操作以下步骤:
步骤六、星机联合双基地合成孔径雷达解耦合
●结合步骤一中的相关初始化参数,并操作以下步骤
●结合步骤一中的相关初始化参数,并操作以下步骤
通过本发明具体实施方式可以看出,本发明克服星机联合双基地合成孔径雷达缺乏有效,简便的频域成像方法的问题,本发明提供了一种基于空域展开的星机联合双基地合成孔径雷达二维频域成像方法,该方法的特点和思路是:利用距离历史在参考点的空域展开和参考点距离历史的二阶时域近似,将相位表示成一个关于慢时间的二次多项式,并利用驻定相位原理得到信号的二维频谱。通过方位和距离压缩得到良好的聚焦效果。最后采用利用散射点位置的线性关系,将距离和方位的坐标从耦合项中剥离出来,得到正确的坐标。从而克服了复杂的二维频谱计算和目标空变补偿的过程,并实现了较大场景下星机联合双基地合成孔径雷达成像。
Claims (1)
1.一种基于空域展开的星机联合SAR二维频域成像方法,其特征是它包括以下步骤:
步骤一、合成孔径雷达成像空间及成像***参数的初始化;
合成孔径雷达成像空间由合成孔径雷达成像空间中的两个相互正交的坐标基确定,定义与发射平台卫星的速度方向平行并在地平面内的单位向量作为合成孔径雷达成像空间的第一个坐标基,即 定义在地平面内,并与合成孔径雷达成像空间的第一个坐标基 垂直的单位向量作为合成孔径雷达成像空间的第二个坐标基,记做
初始化成像***参数包括:发射平台速度矢量,即卫星相对地面速度矢量,记做 接收平台速度矢量,即飞机速度矢量记做 发射平台初始位置矢量,记做 接收平台初始位置矢量,记做 雷达发射电磁波的波数,记做K0,雷达发射基带信号的信号带宽,记做B,雷达发射信号脉冲宽度,记做TP,雷达接收波门持续宽度,记做To,雷达接收***的采样频率,记做fs,雷达***的脉冲重复频率,记做PRF;发射平台初始位置矢量 接收平台初始位置矢量 其中xT yT zT和xR yR zR分别表示发射平台和接收平台的空间三维坐标;采用公式 得到发射平台第n个PRF时刻的位置矢量 采用公式 得到发射平台第n个PRF时刻的位置矢量 n=1,...,Ns表示第n个PRF时刻,Ns为方位向的采样点数;
发射雷达的波束初始指向矢量,记做 接收雷达的波束初始指向矢量 表示雷达的初始波束指向和矢量;发射平台雷达的初始角速度矢量 接收平台雷达的初始角速度矢量 星机联 合双基地合成孔径雷达***初始角速度和矢量 并利用公式
根据星机联合双基地合成孔径雷达***方案和星机联合双基地合成孔径雷达观测方案,星机联合双基地合成孔径雷达成像方法需要的初始化成像***参数均为已知;
步骤二、星机联合双基地合成孔径雷达原始数据进行距离压缩;
我们将接收机接收到的星机联合双基地合成孔径雷达回波数据记做 采用合成孔径雷达标准距离压缩方法对接收到的星机联合双基地合成孔径雷达回波数据进行压缩,得到距离压缩后的星机联合双基地合成孔径雷达数据,记做 n表示第n个PRF时刻,r为距离自变量, 表示待观测的点目标向量;
步骤三、星机联合双基地合成孔径雷达原始数据进行二维傅里叶变换
步骤四、星机联合双基地合成孔径雷达方位参考函数补偿
步骤五、星机联合双基地合成孔径雷达非线性相位补偿
其中 表示载波波数, 表示距离波数,K′=K-K0, 为方位波数;C为光速,f0为载波频率,f为对应于快时间的频率,fa为对应于慢时间的多普勒频率,A表示初始距离历史;B表示等效速度;D表示等效加速度; 符号“||”表示求模运算;
步骤六、星机联合双基地合成孔径雷达解耦合
由于移变双基地合成孔径雷达成像距离向和方位向坐标系的非正交性将导致移变双基地合成孔径雷达图像失真;因此需要分别将距离相位和方位相位中的 耦合项剥离出来,即解耦合;
将步骤五中得到的关于待观测散射点坐标线性相位的数据 沿着行即距离向做逆傅里叶变换,得到距离-多普勒域数据 结合步骤一中的相关初始化参数,将数据 中的第ii列的数据分别乘以方位解耦合函数 其中 表示第ii个快时间的取值,ii=1...Nr,Nr为距离采样点;fs为雷达接收***的采样频率;得到方位解耦合后的数据
将得到的方位解耦合后的数据 沿列做逆傅里叶变换,得到距离频率-方位时间域数据 对数据 中的第jj行数据分别乘以距离解耦合函数 其中 表示第jj个慢时间的取值,jj=1...Ns,Ns为方位采样点,PRF为脉冲重复频率, 为常系数;得到距离解耦后的数据,并对其沿行做逆傅里叶变换,得到最终的成像结果
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010190870579A CN101833095B (zh) | 2010-04-14 | 2010-04-14 | 一种基于空域展开的星机联合sar二维频域成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010190870579A CN101833095B (zh) | 2010-04-14 | 2010-04-14 | 一种基于空域展开的星机联合sar二维频域成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101833095A CN101833095A (zh) | 2010-09-15 |
CN101833095B true CN101833095B (zh) | 2012-06-27 |
Family
ID=42717234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010190870579A Expired - Fee Related CN101833095B (zh) | 2010-04-14 | 2010-04-14 | 一种基于空域展开的星机联合sar二维频域成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101833095B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012510A (zh) * | 2010-09-25 | 2011-04-13 | 哈尔滨工业大学 | 基于时间——相位导数分布的逆合成孔径雷达成像方法 |
CN102004250B (zh) * | 2010-10-28 | 2012-08-08 | 电子科技大学 | 基于频域展开的星机联合双基地合成孔径雷达成像方法 |
CN102226841B (zh) * | 2011-03-29 | 2013-01-23 | 西安电子科技大学 | 基于高阶多项式距离方程的同步轨道sar成像方法 |
CN102819020B (zh) * | 2012-08-17 | 2013-12-18 | 北京航空航天大学 | 一种俯冲模型方位向非线性变标的合成孔径雷达成像方法 |
CN102854507B (zh) * | 2012-09-12 | 2014-04-09 | 电子科技大学 | 一种基于gpu后向投影双站合成孔径雷达成像方法 |
CN102914775B (zh) * | 2012-10-10 | 2014-04-02 | 北京理工大学 | 一种改进的双基地合成孔径雷达nlcs成像算法 |
CN103226191B (zh) * | 2013-04-17 | 2015-03-18 | 中国科学院电子学研究所 | 基于bp算法的高分辨率p波段sar空变运动误差补偿算法 |
CN103543452B (zh) * | 2013-09-27 | 2015-08-19 | 电子科技大学 | 一种基于多普勒频率展开的双基地合成孔径雷达成像方法 |
CN104199030B (zh) * | 2014-07-22 | 2016-09-28 | 广西大学 | 旋转式合成孔径雷达频域成像方法 |
CN106908789B (zh) * | 2017-02-28 | 2019-08-02 | 西安电子科技大学 | 一种基于空间域支撑区融合的sar成像方法 |
CN107219524B (zh) * | 2017-05-11 | 2020-01-07 | 西北工业大学 | 一种基于全局最小相位近似的sar成像优化方法 |
CN108279406B (zh) * | 2018-04-16 | 2020-06-05 | 北京无线电测量研究所 | 一种雷达的解耦合方法、装置、***及存储介质 |
CN108983238A (zh) * | 2018-07-27 | 2018-12-11 | 山东航天电子技术研究所 | 一种星载Ka波段SAR实时成像*** |
CN110531354B (zh) * | 2019-10-14 | 2021-07-20 | 上海无线电设备研究所 | 一种频控扫描雷达色散信号的二维成像方法 |
CN111220981B (zh) * | 2020-01-20 | 2022-12-02 | 西安电子科技大学 | 基于非正交非线性坐标系输出的中轨星载sar成像方法 |
CN111443339B (zh) * | 2020-05-11 | 2022-04-01 | 中国科学院电子学研究所 | 一种双基sar空变校正成像方法、装置、设备及存储介质 |
CN111965637B (zh) * | 2020-08-20 | 2023-06-30 | 西安电子科技大学 | 基于双/多基雷达弹体解析定位方法 |
CN112083417B (zh) * | 2020-09-18 | 2022-02-01 | 电子科技大学 | 基于波数域拼接的分布式雷达成像拓扑设计方法 |
CN113189588B (zh) * | 2021-04-30 | 2022-05-03 | 电子科技大学 | 集群无人机合成孔径雷达高帧率成像方法 |
CN113391311B (zh) * | 2021-06-21 | 2022-08-09 | 电子科技大学 | 一种分布式雷达广义孔径合成方法 |
CN113506233B (zh) * | 2021-07-08 | 2024-04-19 | 西安电子科技大学 | 基于深度学习的sar自聚焦方法 |
CN114726399B (zh) * | 2022-03-14 | 2023-05-23 | 中国电子科技集团公司第十研究所 | 一种非地面网络的空频二维跳频传输***及传输方法 |
CN114895267B (zh) * | 2022-04-08 | 2024-07-09 | 西安电子科技大学 | 一种基于斜距补偿的星载多模式sar回波快速生成方法 |
CN115902784B (zh) * | 2022-12-12 | 2023-09-08 | 电子科技大学 | 大时宽带宽积雷达匀加速运动超高速目标积累检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101369018A (zh) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | 一种星机联合双基地合成孔径雷达频域成像方法 |
CN101369017A (zh) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | 一种移变模式双基地合成孔径雷达成像方法 |
CN101561503A (zh) * | 2008-04-17 | 2009-10-21 | 电子科技大学 | 一种双站线阵三维成像合成孔径雷达快速成像方法 |
-
2010
- 2010-04-14 CN CN2010190870579A patent/CN101833095B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101369018A (zh) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | 一种星机联合双基地合成孔径雷达频域成像方法 |
CN101369017A (zh) * | 2007-08-17 | 2009-02-18 | 电子科技大学 | 一种移变模式双基地合成孔径雷达成像方法 |
CN101561503A (zh) * | 2008-04-17 | 2009-10-21 | 电子科技大学 | 一种双站线阵三维成像合成孔径雷达快速成像方法 |
Non-Patent Citations (3)
Title |
---|
Kefeng Yang等.A Two-Dimensional Spectrum for General Bistatic SAR Processing.《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》.2010,第7卷(第1期), * |
Shi Jun等.Principle and methods on bistatic SAR signal processing via time correlation.《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》.2008,第46卷(第10期), * |
王卫国等.基于GPS信号的双基地SAR成像方法研究.《装备指挥技术学院学报》.2010,第21卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN101833095A (zh) | 2010-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101833095B (zh) | 一种基于空域展开的星机联合sar二维频域成像方法 | |
US11715228B2 (en) | Imaging systems and related methods including radar imaging with moving arrays or moving targets | |
US11726197B2 (en) | Systems and methods for efficient targeting | |
Ponce et al. | Fully polarimetric high-resolution 3-D imaging with circular SAR at L-band | |
Camp et al. | Wideband radar for ballistic missile defense and range-Doppler imaging of satellites | |
CN102313887B (zh) | 一种星机联合双基地合成孔径雷达成像方法 | |
Chang et al. | SAR image simulation with application to target recognition | |
CN102004250B (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
CN103941243B (zh) | 一种基于sar三维成像的自旋式飞行器测高方法 | |
KR102090318B1 (ko) | 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN103576137B (zh) | 一种基于成像策略的多传感器多目标定位方法 | |
Zhang et al. | Matrix completion for downward-looking 3-D SAR imaging with a random sparse linear array | |
US20130038486A1 (en) | Generating radar cross-section signatures | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN107544068A (zh) | 一种基于频域bp的图像域宽带合成方法 | |
CN106093932A (zh) | 一种扫描波束高分辨率的雷达散射计 | |
CN102854505A (zh) | 一种加权稀疏驱动自聚焦sar成像方法 | |
CN103033811B (zh) | 一种基于类单基等效的双基地合成孔径雷达成像方法 | |
CN101710173B (zh) | 一种星机联合双基地合成孔径雷达时域成像方法 | |
Vu et al. | Fast backprojection algorithm for UWB bistatic SAR | |
Chen et al. | Forward looking imaging of airborne multichannel radar based on modified iaa | |
CN102890270A (zh) | 固定站双基地合成孔径雷达回波模拟方法 | |
Marais et al. | Space based GMTI using scanned pattern interferometric radar (SPIR) | |
Bączyk et al. | Moving target imaging in multistatic passive radar | |
Tan et al. | 3-D range stacking algorithm for forward-looking SAR 3-D imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120627 Termination date: 20150414 |
|
EXPY | Termination of patent right or utility model |