CN110095775B - 基于混合坐标系的颠簸平台sar快速时域成像方法 - Google Patents
基于混合坐标系的颠簸平台sar快速时域成像方法 Download PDFInfo
- Publication number
- CN110095775B CN110095775B CN201910354530.9A CN201910354530A CN110095775B CN 110095775 B CN110095775 B CN 110095775B CN 201910354530 A CN201910354530 A CN 201910354530A CN 110095775 B CN110095775 B CN 110095775B
- Authority
- CN
- China
- Prior art keywords
- aperture
- sub
- coordinate system
- distance
- image
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 81
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 230000004927 fusion Effects 0.000 claims description 27
- 238000005070 sampling Methods 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 8
- 238000007906 compression Methods 0.000 claims description 8
- 230000001427 coherent effect Effects 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 32
- 238000010586 diagram Methods 0.000 description 12
- 238000004088 simulation Methods 0.000 description 12
- 238000000034 method Methods 0.000 description 10
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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/9004—SAR image acquisition techniques
- G01S13/9017—SAR image acquisition techniques with time 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于混合坐标系的颠簸平台SAR快速时域成像方法,该方法的具体步骤为:对基带回波信号进行距离向匹配滤波处理,得到距离压缩后的信号;划分子孔径,构建混合坐标系成像网格;将每个子孔径数据在混合坐标系成像网格上进行后向投影;将混合坐标系下的子孔径图像通过混合坐标系下的快速后向投影或者快速分解后向投影算法,得到全空间分辨率的SAR图像。本发明通过将快速时域算法中的极坐标子孔径成像网格替换为混合坐标网格,优化了时域成像算法中的二维插值,提高了时域算法的执行效率,改善时域算法的处理速度和成像质量。
Description
技术领域
本发明属于SAR成像技术领域,特别涉及基于混合坐标系的颠簸平台 SAR快速时域成像方法,可用于无人机等颠簸平台的SAR成像处理。
背景技术
合成孔径雷达(SAR)作为一种主动式微波有源传感器,可以全天时、 全天候、远距离获取观测场景的二维图像,已广泛应用于各种雷达载体。合 成孔径雷达成像算法主要分为两类:频域算法和时域算法。由于时域算法对 成像条件没有限制,且成像过程中不需要脉冲积累,近年来逐渐成为研究热 点。然而时域成像算法的运算复杂度相对较高,导致其在应用中存在一定的 局限性。
为提高时域算法的处理速度,现有思路是利用划分子孔径的思想对原始 后向投影(BP)算法进行改进,并提出了快速后向投影(FBP)算法和快速 分解后向投影(FFBP)算法。FBP算法先将每幅子孔径图像重建在以子孔径 中心为原点的局部极坐标系;再通过距离维和角度维的二维插值将极坐标子 孔径图像变换到笛卡尔坐标系;最后在笛卡尔坐标系将全部子孔径图像进行 相干相加,得到全空间分辨率的SAR图像。FFBP算法先将整个合成孔径按 照一定的分解系数划分成若干较短的子孔径,并将子孔径对应的距离压缩数 据投影到以其孔径中心为原点的局部极坐标网格;再通过距离维和角度维的 二维插值将先前阶段的子孔径图像不断融合,不断提高子孔径图像的角分辨 率;最后通过二维插值将极坐标图像变换到笛卡尔坐标系,得到全空间分辨 率SAR图像。
然而,FBP算法和FFBP算法极坐标子孔径图像递归融合及极坐标到笛卡 尔坐标转换过程中的两维插值在实现过程中占据较大比重,不可避免地增加 计算量,同时引入了插值误差,最终影响图像质量,造成旁瓣抬高和积分旁 瓣比性能变差,在实际应用中存在一定的局限性。
发明内容
为了解决上述问题,本发明的目的是提出一种基于混合坐标系的颠簸平 台SAR快速时域成像方法。本发明通过将快速时域算法中的极坐标子孔径 成像网格替换为混合坐标网格,优化了时域成像算法中的二维插值,提高了 时域算法的执行效率,改善时域算法的处理速度和成像质量。
为了达到上述目的,本发明采用以下技术方案予以解决。
基于混合坐标系的颠簸平台SAR快速时域成像方法,包括以下步骤:
步骤1,采用颠簸平台SAR发射并接收回波信号;对回波信号进行解调, 得到解调后的基带回波信号s0(τ,tm),对解调后的基带回波信号s0(τ,tm)进行距 离向匹配滤波处理,得到距离压缩后的信号s(τ,tm);其中,τ为距离快时间, tm为方位慢时间。
步骤2,设定距离压缩后的全孔径数据s(τ,tm)长度为L,将距离压缩后的 全孔径数据s(τ,tm)均匀划分成K个子孔径数据,其中,每个子孔径的长度 l=L/K;并构建混合坐标系成像网格。
步骤3,将每个子孔径数据在混合坐标系成像网格上进行后向投影,得到 每个子孔径数据对应的混合坐标系下的子孔径图像Ik(y,θk);其中,y表示混 合坐标系成像网格中像素点到理想航迹的最近距离,θk表示垂直飞行平台速 度方向与像素点到第k个子孔径中心连线的夹角。
步骤4,将混合坐标系下的子孔径图像Ik(y,θk)通过混合坐标系下的快速 后向投影或者快速分解后向投影,得到全空间分辨率的SAR图像I(x,y);其 中,x表示像素点到孔径中心距离在速度方向的投影距离。
与现有技术相比,本发明的有益效果为:
(1)本发明通过引入混合坐标系,改善了颠簸平台SAR快速时域成像算 法的处理效率和成像效果。
(2)本发明通过提出混合坐标系,将快速时域算法的极坐标成像网格替 换为混合坐标网格,将子孔径图像形成过程的距离插值与极坐标到笛卡尔坐 标转换(包括极坐标到极坐标的转换)时的距离插值合二为一,使得之后的 坐标转换和图像融合只需要进行角度维插值,提高了时域算法的执行效率, 减少了插值误差,改善了成像质量。
附图说明
下面结合附图和具体实施例对本发明做进一步详细说明。
图1为传统的极坐标系和传统的笛卡尔坐标系示意图;其中,图1(a) 为传统的极坐标系示意图,图1(b)为传统的笛卡尔坐标系示意图;
图2为本发明中的混合坐标系示意图;
图3为本发明的混合坐标系下颠簸平台SAR成像几何模型示意图;
图4为本发明基于混合坐标系的颠簸平台SAR快速时域成像方法的流程 图;
图5为本发明实施例中混合坐标系下图像融合前和融合后的示意图;其 中,图5(a)为本发明实施例中混合坐标系下图像融合前的示意图,图5(b) 为本发明实施例中混合坐标系下图像融合后的示意图;
图6(a)为本发明的混合坐标系下快速时域成像仿真实验中像素点距离 向剖面图,图6(b)本发明的为混合坐标系下快速时域成像仿真实验中像素 点方位向剖面图。
具体实施方式
下面结合附图对本发明的实施例及效果作进一步详细描述。
参照图3所示,为本发明的混合坐标系下颠簸平台SAR成像几何 模型示意图;无人机载合成孔径雷达工作于条带模式,载机以水平速度 v沿x方向(雷达的方位向)飞行,P点为场景中的某一像素点。载机 飞行高度为H,景中心斜距为Rs。
参照图4,为本发明基于混合坐标系的颠簸平台SAR快速时域成像方法 的流程图。该方法的具体实现步骤如下:
步骤1,采用颠簸平台SAR发射并接收回波信号;对回波信号进行解调, 得到解调后的基带回波信号s0(τ,tm),对解调后的基带回波信号s0(τ,tm)进行距 离向匹配滤波处理,得到距离压缩后的信号s(τ,tm)。
具体地,利用无人机载合成孔径雷达发射线性调频信号,利用合成孔径 雷达接收对应的回波信号;对回波信号进行解调,得到解调后的基带回波信 号s0(τ,tm):
其中,τ为距离快时间,tm为方位慢时间;A0为像素点的后向散射系数, ωr(·)和ωa(·)分别为基带回波信号的距离窗函数和方位窗函数,R(tm)为tm时刻 雷达与像素点的瞬时斜距,c表示光速,t0表示波束中心穿越时刻,fc为信号 载频,Kr为距离向调频率。
将基带回波信号s0(τ,tm)进行快速傅立叶变换并乘以距离向匹配滤波函 数,消除距离二次项,再进行逆傅立叶变换,得到距离压缩后的信号s(τ,tm)
步骤2,设定距离压缩后的全孔径数据s(τ,tm)长度为L,将距离压缩后的 全孔径数据s(τ,tm)均匀划分成K个子孔径数据,其中,每个子孔径的长度 l=L/K;构建混合坐标系成像网格。
具体地说,假设距离压缩后的全孔径数据s(τ,tm)长度为L,设定子孔径长 度l,则全孔径数据可以划分为K=L/l个子孔径数据。
传统快速时域成像算法的成像网格是基于极坐标系或者笛卡尔坐标系建 立的,其坐标系的建立具体如下:
其一,参考图1(a),建立极坐标系,其坐标原点位于孔径中心点,场 景中的任一点表示为(r,θ),其中,r为偏离孔径中心的距离,θ为垂直飞行平 台速度方向与像素点到孔径中心连线的夹角,即角度维。
其二,参考图1(b),建立笛卡尔坐标系,其坐标原点位于孔径中心点, 以平台运动方向为x轴方向,以垂直于雷达轨迹方向为y轴方向,场景中任一 点表示为(x,y),其中,x为像素点到孔径中心距离在速度方向的投影距离,即 方位维,y为像素点到孔径中心距离在垂直雷达轨迹方向的投影距离,即像素 点到理想航迹的最近距离,即距离维。笛卡尔坐标与极坐标的转换关系如下:
不同于传统成像网格,本发明的混合坐标系成像网格建立的具体步骤为:
参考图2,建立混合坐标系,混合坐标系以孔径中心为原点,以垂直于雷 达轨迹方向为y轴方向,在距离维和角度维进行等间隔划分,建立起混合坐 标系成像网格(y,θ);y和θ满足关系:
y=r·cosθ;
其中,y表示偏离孔径的距离,即像素点到理想航迹的最近距离,即距离 维,θ表示垂直飞行平台速度方向与像素点到孔径中心连线的夹角,即角度维, r为像素点到孔径中心的距离。
步骤3,将每个子孔径数据在混合坐标系成像网格上进行后向投影,得到 每个子孔径数据对应的混合坐标系下的子孔径图像Ik(y,θk)。
其中,Ik(y,θk)表示第k个混合坐标系下的子孔径图像,k表示第k个子孔 径,θk表示垂直飞行平台速度方向与像素点到第k个子孔径中心连线的夹角; R(y,θk;tm)表示tm时刻雷达天线相位中心到像素点(y,θk)的瞬时斜距,其表达式 为:
所述混合坐标系下快速时域成像算法的子孔径图像的Nyquist采样需求 为:
其中,λ为载波波长,Δy表示混合坐标系成像网格中两个相邻像素点在 距离向的最短距离,Δθk表示混合坐标系成像网格中两个相邻像素点在角度维 的最小间隔角度。
对图像Ik(y,θk)作二维傅立叶变换,并根据驻定相位原理求得混合坐标系 下快速时域成像算法的子孔径图像的Nyquist采样需求,只有满足Nyquist采 样需求时,混合坐标系下的子孔径图像才能较为完美地反映实际场景。
由以上采样条件可知,角度维采样与子孔径长度l有关,而距离采样只取 决于带宽,与子孔径大小无关。
对于固定的成像角,角度维网格划分得越细密,像素间隔越小,角度维 采样率越高。考虑到混合坐标系下的子孔径图像具有距离全分辨和角度维低 分辨的特点,可以通过减小角度维的采样点数来提高运算效率。
步骤4,将混合坐标系下的子孔径图像Ik(y,θk)通过混合坐标系下的快速 后向投影或快速分解后向投影,得到全空间分辨率的SAR图像I(x,y)。
一种实施方式下,所述快速后向投影的具体步骤为:
首先,采用角度维插值将每个混合坐标系下的子孔径图像Ik(y,θk)转换成 对应的笛卡尔坐标系下的子孔径图像Ik(x,y):
其中,R(x,y;tm)表示tm时刻雷达天线相位中心到笛卡尔坐标系成像网格中 像素点(x,y)的瞬时斜距。
由于每个混合坐标系下的子孔径图像均在混合坐标系下的成像网格上进 行后向投影成像,其距离维坐标y表示像素点到第k个子孔径的最近距离(像 素点到理想航迹的最近距离),即与笛卡尔坐标系中距离维坐标y含义相同, 因此在进行坐标转换时不需要进行距离维插值,只进行角度维插值即可。
其次,对K个笛卡尔坐标系下的子孔径图像Ik(x,y)进行相干相加,得到 一幅全空间分辨率的SAR图像I(x,y):
其中,x表示像素点到孔径中心距离在速度方向的投影距离。
另一实施方式下,所述快速分解后向投影的具体步骤为:
首先,对K个混合坐标系下的子孔径图像Ik(y,θk)进行递归融合,得到一 幅混合坐标下的全空间分辨率的SAR图像I(G)(y,θ);
其次,将混合坐标下的全空间分辨率的SAR图像I(G)(y,θ),通过角度维 插值转换成笛卡尔坐标下的全空间分辨率的SAR图像,即为全空间分辨率的 SAR图像I(x,y)。其表达式为:
上述步骤中,所述递归融合为基二融合。
所述基二融合的具体步骤为:
假设将K个混合坐标系下的子孔径图像的基二融合划分为G个阶段,基 二融合的初始阶段划分得到的子孔径数目为K(1)=2G-1,初始阶段的第k个子孔 径的雷达天线相位中心为对应的混合坐标系成像网格的 角度维坐标划分为
其中,G为大于1的整数;K(1)表示初始阶段的子孔径数,初始阶段每个 子孔径长度为l(1)=L/K(1)。
初始阶段的图像融和的具体步骤如下:
参考图5(b),为混合坐标系下融合后图像示意图,从图中可以看出, 融合后的子孔径长度等于融合前两个子孔径长度之和,融合后图像的角度维 分辨率为融合前的2倍。由于融合前的两个子孔径图像均在混合坐标系成像 网格中进行成像,其距离维的划分相同,即均为y,表示像素点到第k个子孔 径的最近距离,因此在进行图像融合时不需要进行距离维插值,只进行角度 维插值即可。
对初始阶段的第2k-1个和第2k个子孔径的混合坐标系下的子孔径图像进 行角度维插值融合,得到第2阶段第k个子孔径的混合坐标系下的子孔径图 像:
将初始阶段的K(1)个子孔径图像分成K(1)/2对,对每对子孔径图像分别进 行角度维插值融合,则得到K(1)/2个子孔径图像,即为第2阶段的待融合图像。
对第2阶段的待融合图像进行融合得到K(1)/4个子孔径图像,作为第3阶 段的待融合图像,依次类推,最终得到一幅融合后的图像即为混合坐标下的 全空间分辨率的SAR图像I(G)(y,θ)。
验证过程:
为了说明本发明的快速时域成像方法采用混合坐标系的优势,此处对采 用混合坐标系所带来的插值运算效率的提升进行分析。
以FBP算法中极坐标系成像网格得到的子孔径成像结果(简称极坐标子 孔径图像)到笛卡尔坐标子孔径图像的坐标转换的运算量为例。
假设孔径长度为Npulse,笛卡尔坐标系下重建的成像网格点数为N×N,每 个子孔径长度为m,则总共需划分Npulse/m个子孔径。设极坐标子孔径图像的 像素点数为N×m,将该子孔径图像投影到笛卡尔坐标网格,其距离维和角度 维插值均采用升采样后的线性插值,插值倍数为Num,计算得到距离维插值 的浮点运算为(一次实数加法或一次实数乘法作为一次浮点运算):
Operation_r=m·(5·N·log2N+5·Num·N·log2(Num·N)+6·N);
角度维插值的浮点运算为:
Operation_a=N·(5·m·log2m+5·Num·m·log2(Num·m)+6·N);
总运算量为Operation_r+Operation_a。采用混合成像坐标系后,子孔径图像坐标转换(即从混合坐标系下的子孔径图像转换为笛卡尔坐标子孔径图像)只 包含角度维插值,故此时运算量为Operation_a。由此可得混合坐标系下的运算 效率为
而针对FFBP算法,由于每一级处理在采用混合坐标系后均可以节省距离 维插值运算,效率提升会更加明显。
仿真实验
下面通过像素点仿真成像实验进一步说明本发明的正确性和有效性。
1)像素点仿真成像仿真条件
表1雷达参数
雷达系统仿真参数如表1所示,参照图2,为混合坐标系下SAR成像几 何模型示意图;其中,X轴为颠簸SAR平台(载机)运动方向,Y轴垂直于 雷达平台运动方向,Z轴为平台高度方向,SAR平台高度为H。载机以速度v 沿X轴匀速直线飞行,设定Rs为场景中心点对应的距离。
2)仿真内容
仿真实验:为验证本发明的有效性,此处通过快速时域成像法在混合坐 标系下的仿真实验,对其运算效率的改进作分析。
3)仿真结果分析
参照图6a,为本发明的混合坐标系下快速时域成像仿真实验中像素点距 离向剖面图;参照图6b,为本发明的混合坐标系下快速时域成像仿真实验中 像素点方位向剖面图。设孔径长度为4096,分别将其划分成长度为32、64、 128、256的子孔径(SA),并建立4096*4096的成像网格。其中,倒三角的 虚线表示极坐标(PC)下的成像结果,菱形的点线表示混合坐标系(HC)下 的成像结果,二者几乎完全相同,且混合坐标系下的成像结果峰值旁瓣比 (PSLR)优于极坐标系成像结果。
表2处理时间对比
PC/s | HC/s | 运算效率 | |
SA=32 | 18.36 | 10.76 | 41.4% |
SA=64 | 12.98 | 7.44 | 42.7% |
SA=128 | 9.25 | 4.92 | 46.8% |
SA=256 | 9.85 | 5.28 | 46.4% |
通过计算图像坐标转换时的插值运算时间,当子孔径长度为64时,混合 坐标系下的运算时间为7.44s,而极坐标系下的运算时间为12.98s,其运算效 率为42%。其他子孔径长度的插值运算时间对比如表2所示,由表2可知, 混合坐标系极大地提高了处理时间,减少了算法的复杂度,提高了运算效率。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部 分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一 计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例 的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种 可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并 不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围 内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此, 本发明的保护范围应以所述权利要求的保护范围为准。
Claims (5)
1.基于混合坐标系的颠簸平台SAR快速时域成像方法,其特征在于,包括以下步骤:
步骤1,采用颠簸平台SAR发射并接收回波信号;对回波信号进行解调,得到解调后的基带回波信号s0(τ,tm),对解调后的基带回波信号s0(τ,tm)进行距离向匹配滤波处理,得到距离压缩后的信号s(τ,tm);其中,τ为距离快时间,tm为方位慢时间;
步骤2,设定距离压缩后的全孔径数据s(τ,tm)长度为L,将距离压缩后的全孔径数据s(τ,tm)均匀划分成K个子孔径数据,其中,每个子孔径的长度l=L/K;并构建混合坐标系成像网格;
所述构建混合坐标系成像网格的具体步骤为:
建立混合坐标系,混合坐标系以孔径中心为原点,以垂直于雷达轨迹方向为y轴方向,在距离维和角度维进行等间隔划分,建立起混合坐标系成像网格(y,θ);其中,y和θ满足关系:
y=r·cosθ;
其中,y表示像素点到理想航迹的最近距离,即距离维,θ表示垂直飞行平台速度方向与像素点到孔径中心连线的夹角,即角度维,r为像素点到孔径中心的距离;
步骤3,将每个子孔径数据在混合坐标系成像网格上进行后向投影,得到每个子孔径数据对应的混合坐标系下的子孔径图像Ik(y,θk);
其中,y表示混合坐标系成像网格中像素点到理想航迹的最近距离,θk表示垂直飞行平台速度方向与像素点到第k个子孔径中心连线的夹角;
所述混合坐标系下的子孔径图像Ik(y,θk)的表达式为:
其中,Ik(y,θk)表示第k个混合坐标系下的子孔径图像,k表示第k个子孔径,θk表示垂直飞行平台速度方向与像素点到第k个子孔径中心连线的夹角;v表示载机沿x方向飞行的水平速度;R(y,θk;tm)表示tm时刻雷达天线相位中心到像素点(y,θk)的瞬时斜距,其表达式为:
步骤4,将混合坐标系下的子孔径图像Ik(y,θk)通过混合坐标系下的快速后向投影或者快速分解后向投影,得到全空间分辨率的SAR图像I(x,y);
其中,x表示像素点到孔径中心距离在速度方向的投影距离;
所述快速后向投影的具体步骤为:
首先,采用角度维插值将每个混合坐标系下的子孔径图像Ik(y,θk)转换成对应的笛卡尔坐标系下的子孔径图像Ik(x,y):
其中,R(x,y;tm)表示tm时刻雷达天线相位中心到笛卡尔坐标系成像网格中像素点(x,y)的瞬时斜距;
其次,对K个笛卡尔坐标系下的子孔径图像Ik(x,y)进行相干相加,得到一幅全空间分辨率的SAR图像I(x,y):
其中,x表示像素点到孔径中心距离在速度方向的投影距离;
所述快速分解后向投影的具体步骤为:
首先,对K个混合坐标系下的子孔径图像Ik(y,θk)进行递归融合,得到一幅混合坐标下的全空间分辨率的SAR图像I(G)(y,θ);
其次,将混合坐标下的全空间分辨率的SAR图像I(G)(y,θ),通过角度维插值转换成笛卡尔坐标下的全空间分辨率的SAR图像,即为全空间分辨率的SAR图像I(x,y)为:
5.根据权利要求1所述的基于混合坐标系的颠簸平台SAR快速时域成像方法,其特征在于,所述递归融合为基二融合;所述基二融合的具体步骤为:
假设将K个混合坐标系下的子孔径图像的基二融合划分为G个阶段,则初始阶段划分得到的子孔径数目为K(1)=2G-1,初始阶段的第k个子孔径的雷达天线相位中心为对应的混合坐标系成像网格的角度维坐标划分为其中,G为大于1的整数;K(1)表示初始阶段的子孔径数,初始阶段每个子孔径长度为l(1)=L/K(1);为初始阶段的第k个子孔径的雷达天线相位中心;
则初始阶段的图像融和的具体步骤为:
其次,对初始阶段的第2k-1个和第2k个子孔径的混合坐标系下的子孔径图像进行角度维插值融合,得到第2阶段第k个子孔径的混合坐标系下的子孔径图像:
最后,将初始阶段的K(1)个子孔径图像分成K(1)/2对,对每对子孔径图像分别进行角度维插值融合,则得到K(1)/2个子孔径图像,即为第2阶段的待融合图像;
依次类推,对第2阶段的待融合图像进行第2到G阶段的图像融合,得到一幅融合后的图像即为混合坐标下的全空间分辨率的SAR图像I(G)(y,θ)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354530.9A CN110095775B (zh) | 2019-04-29 | 2019-04-29 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354530.9A CN110095775B (zh) | 2019-04-29 | 2019-04-29 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110095775A CN110095775A (zh) | 2019-08-06 |
CN110095775B true CN110095775B (zh) | 2023-03-14 |
Family
ID=67446324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910354530.9A Active CN110095775B (zh) | 2019-04-29 | 2019-04-29 | 基于混合坐标系的颠簸平台sar快速时域成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110095775B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456354A (zh) * | 2019-09-17 | 2019-11-15 | 上海无线电设备研究所 | 一种太赫兹圆迹sar快速后向投影成像方法 |
CN111537999B (zh) * | 2020-03-04 | 2023-06-30 | 云南电网有限责任公司电力科学研究院 | 一种稳健高效的分解投影自动聚焦方法 |
CN111736151B (zh) * | 2020-06-16 | 2022-03-04 | 西安电子科技大学 | 一种高效全局直角坐标投影融合的改进ffbp成像方法 |
CN113189588B (zh) * | 2021-04-30 | 2022-05-03 | 电子科技大学 | 集群无人机合成孔径雷达高帧率成像方法 |
CN113552564B (zh) * | 2021-06-23 | 2022-10-25 | 南昌大学 | Sar对复杂地形场景时域快速成像方法、系统、终端及应用 |
CN114578355B (zh) * | 2022-03-03 | 2022-10-21 | 西安电子科技大学 | 高超声速飞行器合成孔径雷达快速时域成像方法 |
CN116087954B (zh) * | 2023-04-06 | 2023-06-30 | 中国科学院空天信息创新研究院 | 一种实孔径系统近距离偏焦合成孔径成像方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6772181B1 (en) * | 1999-10-29 | 2004-08-03 | Pentomics, Inc. | Apparatus and method for trigonometric interpolation |
CN101915920A (zh) * | 2010-07-02 | 2010-12-15 | 北京航空航天大学 | 一种地球同步轨道合成孔径雷达卫星的高分辨率成像方法 |
CN102288963A (zh) * | 2011-07-21 | 2011-12-21 | 西安电子科技大学 | 基于子孔径参数估计的双基地isar图像融合方法 |
CN104049241A (zh) * | 2014-05-29 | 2014-09-17 | 电子科技大学 | 目标位置坐标未知的双基地合成孔径雷达的空间同步方法 |
CN104833974A (zh) * | 2015-05-08 | 2015-08-12 | 西安电子科技大学 | 基于图像谱压缩的sar成像快速后向投影方法 |
CN105842694A (zh) * | 2016-03-23 | 2016-08-10 | 中国电子科技集团公司第三十八研究所 | 一种基于ffbp sar成像的自聚焦方法 |
CN106154269A (zh) * | 2016-09-23 | 2016-11-23 | 南京特艺科技有限公司 | 一种适用于圆周合成孔径雷达快速时域成像方法 |
CN106802416A (zh) * | 2017-02-21 | 2017-06-06 | 电子科技大学 | 一种快速因式分解后向投影sar自聚焦方法 |
CN108205135A (zh) * | 2018-01-22 | 2018-06-26 | 西安电子科技大学 | 基于无插值融合快速后向投影的雷达视频成像方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2120063A1 (en) * | 2008-05-15 | 2009-11-18 | The European Community, represented by the European Commission | Radar-imaging of a scene in the far-field of a one-or two-dimensional radar array |
CN104007440B (zh) * | 2014-06-03 | 2016-05-18 | 西安电子科技大学 | 一种加速分解后向投影聚束合成孔径雷达成像方法 |
-
2019
- 2019-04-29 CN CN201910354530.9A patent/CN110095775B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6772181B1 (en) * | 1999-10-29 | 2004-08-03 | Pentomics, Inc. | Apparatus and method for trigonometric interpolation |
CN101915920A (zh) * | 2010-07-02 | 2010-12-15 | 北京航空航天大学 | 一种地球同步轨道合成孔径雷达卫星的高分辨率成像方法 |
CN102288963A (zh) * | 2011-07-21 | 2011-12-21 | 西安电子科技大学 | 基于子孔径参数估计的双基地isar图像融合方法 |
CN104049241A (zh) * | 2014-05-29 | 2014-09-17 | 电子科技大学 | 目标位置坐标未知的双基地合成孔径雷达的空间同步方法 |
CN104833974A (zh) * | 2015-05-08 | 2015-08-12 | 西安电子科技大学 | 基于图像谱压缩的sar成像快速后向投影方法 |
CN105842694A (zh) * | 2016-03-23 | 2016-08-10 | 中国电子科技集团公司第三十八研究所 | 一种基于ffbp sar成像的自聚焦方法 |
CN106154269A (zh) * | 2016-09-23 | 2016-11-23 | 南京特艺科技有限公司 | 一种适用于圆周合成孔径雷达快速时域成像方法 |
CN106802416A (zh) * | 2017-02-21 | 2017-06-06 | 电子科技大学 | 一种快速因式分解后向投影sar自聚焦方法 |
CN108205135A (zh) * | 2018-01-22 | 2018-06-26 | 西安电子科技大学 | 基于无插值融合快速后向投影的雷达视频成像方法 |
Non-Patent Citations (3)
Title |
---|
Integrating Autofocus Techniques With Fast Factorized Back-Projection for High-Resolution Spotlight SAR Imaging;Lei Zhang;《IEEE Geoscience and Remote Sensing Letters》;20130719;第10卷(第6期);全文 * |
一种改进的快速分解后向投影SAR成像算法;李浩林等;《电子与信息学报》;20130630(第06期);全文 * |
圆周合成孔径雷达的快速时域成像算法;陈乐平等;《国防科技大学学报》;20180430;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110095775A (zh) | 2019-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110095775B (zh) | 基于混合坐标系的颠簸平台sar快速时域成像方法 | |
CN108107431B (zh) | 一种圆柱扫描sar三维成像快速实现方法 | |
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN104391297B (zh) | 一种划分子孔径pfa雷达成像方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN111856461B (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
CN104698459B (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
CN111190181B (zh) | 颠簸平台无人机载sar实时成像处理方法 | |
CN111736151B (zh) | 一种高效全局直角坐标投影融合的改进ffbp成像方法 | |
CN103293528B (zh) | 扫描雷达超分辨成像方法 | |
CN109188432A (zh) | 一种平行双基聚束sar快速bp成像方法 | |
CN110208796B (zh) | 基于奇异值逆滤波的扫描雷达超分辨成像方法 | |
CN105093225A (zh) | 基于双重稀疏约束的逆合成孔径雷达自聚焦成像方法 | |
CN108957452A (zh) | 一种合成孔径雷达自适应ffbp成像方法 | |
CN102798861B (zh) | 一种基于最优图像空间双基地合成孔径雷达成像方法 | |
CN106154269A (zh) | 一种适用于圆周合成孔径雷达快速时域成像方法 | |
CN111220979B (zh) | 一种用于曲线合成孔径雷达的成像方法 | |
CN114325704B (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN116679265A (zh) | 一种适用于地形起伏场景下的sar时域快速回波模拟方法 | |
CN105044720A (zh) | 一种基于直角坐标系的后向投影成像方法 | |
CN109188436B (zh) | 适用于任意平台轨迹的高效双基sar回波生成方法 | |
CN102928839A (zh) | 多通道波束指向sar的全孔径成像方法 | |
CN112285707B (zh) | 基于gps导航信号的无源多基地高分辨成像方法 | |
CN116359921A (zh) | 基于加速轨迹双基前视合成孔径雷达的快速时域成像方法 | |
CN105068071A (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 |