SAR成像系列：【6】合成孔径雷达（SAR）成像算法-调频变标（Chirp Schaling）算法（附Matlab代码）

HIT夜枭

已于 2022-08-25 16:57:59 修改

阅读量3k

点赞数 2

分类专栏：合成孔径雷达（SAR）成像与目标识别文章标签：算法

于 2022-08-25 16:56:50 首次发布

本文链接：https://blog.csdn.net/zr2006_7/article/details/126481770

版权

合成孔径雷达（SAR）成像与目标识别专栏收录该内容

18 篇文章 327 订阅

订阅专栏

（1）算法介绍

距离多普勒（RD）算法虽然简单，但是有两个缺点：第一距离徙动矫正运算量大，补偿时间长；第二，二次距离压缩时，R变化大，对高阶方位频率无法有效补偿。

调频变标（Chirp Schaling）算法通过对线性调频信号进行频率调制，实现尺度变换和平移，并利用相位相乘代替插值运算完成距离徙动补偿，这样不仅提高了算法的计算效率，而且解决二次距离压缩存在的方位频率、依赖问题。

调频变标（Chirp Schaling）算法的核心是矫正不同距离门上的RCMC差量（补余RCMC），这样所有信号具有相同的RCM；然后通过相位相乘完成RCM的矫正（一致RCMC）。

在实际工程中为了提高运算效率，一般在信号处理系统中首先完成距离压缩操作，这样就损失了距离向的调制信息，需要通过距离延拓重建这些信息，或在RCMC时附加一致方位压缩。

（2）一致RCMC与补余RCMC

在上一篇文章中，给出了距离多普勒域距离等式的精确表达式：

$R_{rd}(R_{0},f_{t})=\frac{R_{0}}{D(f_{t},V_{a})}$

其中， $D(f_{t},V_{a})=\sqrt{1-\frac{c^{2}f_{t}^{2}}{4V_{a}^{2}f_{0}^{2}}}$ .它的二阶泰勒展开近似表达为：

$R_{rd}(R_{0},f_{t})\approx R_{0}+\frac{R_{0}c^{2}f_{t}^{2}}{8V_{a}^{2}f_{0}^{2}}$

它是最短斜距R0的线性函数，也是方位频率ft的二次函数。在成像中，选择场景中心为参考点，如下图。

图中，近点、场景中心点、远点的整体RCM示意图如下图：

定义场景中心的距离徙动为一致RCM，即所有目标的一致RCM都是相同的。整体RCM与一致RCM的差值定义为补余RCM，补余RCM与目标距离相关。

为减少RCMC的平移，定义整体RCM为：目标距离与其参考方位频率处的距离差，表达式为：

$RCM_{total}(R_{0},f_{t})=\frac{R_{0}}{D(f_{t},V_{a})}-\frac{R_{0}}{D(f_{t_{ref}},V_{a})}$

一致RCM为：

$RCM_{bulk}(f_{t})=\frac{R_{ref}}{D(f_{t},V_{ref})}-\frac{R_{ref}}{D(f_{t_{ref}},V_{ref})}$

补余RCM为：

$RCM_{diff}(R_{0},f_{t})=\frac{R_{0}}{D(f_{t},V_{a})}-\frac{R_{0}}{D(f_{t_{ref}},V_{a})}-\frac{R_{ref}}{D(f_{t},V_{ref})}+\frac{R_{ref}}{D(f_{t_{ref}},V_{ref})}$

（3）调频变标（Chirp Schaling）算法

①将回波信号变换到距离多普勒域，并与变标方程相乘，完成补余RCM矫正，变标方程如下：

$s_{cs}(\tau',f_{t})=exp(j\pi K_{m}[\frac{D(f_{t_{ref}},V_{ref})}{D(f_{t},V_{ref})}-1]\tau '^{2})$

②距离向FFT将数据变换到二维频域进行距离匹配滤波、二次距离压缩和一致RCM矫正，上述操作可以合并为一次二维频域相位相乘。二维补偿相位函数为：

$s_{rc}(f_{\tau},f_{t})=exp(j\frac{\pi D(f_{t},V_{a})}{K_{m}D(f_{t_{ref}},V_{a})}f_{\tau}^{2})exp(j\frac{4\pi}{c}[\frac{1}{D(f_{t},V_{ref})}-\frac{1}{D(f_{t_{ref}},V_{ref})}]R_{ref}f_{\tau})$

③距离向IFFT将数据变回到距离多普勒域，进行方位匹配滤波和附加相位矫正，匹配滤波和相位矫正函数为：

$s_{ac}(f_{t})=exp(j\frac{4\pi R_{0}f_{0}D(f_{t},V_{a})}{c})exp(-j\frac{4\pi K_{m}}{c^{2}}[1-\frac{D(f_{t},V_{ref})}{D(f_{t_{ref}},V_{ref})}][\frac{R_{0}}{D(f_{t},V_{a})}-\frac{R_{ref}}{D(f_{t},V_{a})}]^{2})$

④方位向IFFT将信号变回时域，得到最终SAR图像。

（4）算法仿真

仿真参数如下表所示：

载频	9.875GHz
雷达高度	500Km
斜视角	5°
带宽	30MHz
雷达速度	7100m/s
目标个数	5
脉冲持续时间	20us

仿真结果：

回波信号实部：

距离多普勒域补余RCMC后的信号实部：

变换到二维频域的信号实部：

距离压缩、SRC和一致RCMC操作后的二维频域信号实部：

距离压缩、SRC和一致RCMC操作后的距离多普勒域信号实部：

方位匹配滤波和附加相位矫正后的距离多普勒域信号实部：

成像结果：

Matlab部分代码：

clear;close all;clc;
SNf = -15; % 信噪比

c=3e8;
f0 = 9.875e9; % 雷达工作频率Hz
lamda = c/f0; % 雷达工作波长m
H = 500e3; % 高度
Yc=750e3; % 成像区域中线
R0 = sqrt(Yc^2+H^2); % 中心斜距

theta = asind(H/R0); % 下视角
alpha = 5; %波束斜视角
Rc = R0/cosd(alpha); %景中心到天线距离
Xc = R0*tand(alpha); %景中心的方位坐标
Br=30e6; % 带宽
Vr = 7100; % 雷达有效速度m/s
Tr =20e-6; % 脉冲持续时间s
Kr = Br/Tr; % 线性调频率
Fr = 1.2*Br; % 距离采样频率，1.2为过采样率
Ts = 1/Fr; % 距离采样时间间隔s

Nk = ceil((2 * 1200/ c + Tr) / Ts); %距离向前后500m
Nf = 2^nextpow2(Nk); % 距离向的采样点个数
tf_ori = [-Nf/2:1:Nf/2-1]*Ts; % 距离向采样时序
tf = [-Nf/2:1:Nf/2-1]*Ts+2*Rc/c; % 实际快时间采样值
fr = ( -Nf/2 : Nf/2-1 )*( Fr/Nf ); % 距离频率
La = 6; % 等效天线尺寸
%前斜视时，采用聚束增加采样个数和方位带宽系数为
Ls = R0*lamda/La*cosd(alpha); % 合成孔径时长度m,Ls=（0.886*R0*lamda）/(La*cos(Theta)）
Ta = Ls/Vr; % 目标照射时间s
% Ta = 0.8; % 目标照射时间s
% Ls =Ta*Vr;
Ka = -2 * Vr^2*(cosd(alpha)).^2 / (lamda * R0); % 方位多普勒调频率Hz
% Ba = 2*Vr*cosd(alpha)/La;%多普勒带宽
Ba=abs(Ka*Ta); % 多普勒频率调制带宽
PRF = ceil(1.2*Ba); % 方位采样率Hz
% PRF = 1000; % 方位采样率Hz
PRT = 1/PRF; % 方位向采样时间间隔s
Ns = 2^nextpow2((200/Vr+Ta)*PRF); % 方位向的采样点个数左右各100m
ts = [-Ns/2 : (Ns/2 - 1)] * PRT; % 方位向采样时序
% ts = [-Ns/2 : (Ns/2 - 1)] * PRT+Xc/Vr; % 实际方位向采样时序
fdc = 2*Vr*sind(alpha)/lamda;%多普勒中心频率式4.33
fdc_ref=fdc;%参考目标的多普勒中心频率
fa=fdc+[-Ns/2:Ns/2-1].'/Ns*PRF; % 方位频率
% fa=[-Ns/2:Ns/2-1].'/Ns*PRF; % 方位频率
% 理论分辨率
rho_r=c/2/Br;
rho_a=Vr*PRT;
% rho_a=La/2;

% 目标参数
X0 = [Xc-100 Xc+100 Xc Xc-100 Xc+100]; % 目标1位置坐标
Z0 = [0 0 0 0 0];
Y0 = [Yc+100 Yc+100 Yc Yc-100 Yc-100];
NT=size(X0,2);

%目标各自的波束中心穿越时刻
eta_c=(X0-sqrt(H^2+Y0.^2)*tand(alpha))/Vr;

%%================================================================
%%生成回波信号
Sb = zeros(Ns,Nf);
sigma = 1; % 回波幅度
for ii=1:NT
R = sqrt((Vr*ts-X0(ii)).^2+Y0(ii).^2+(Z0(ii)-H).^2);
tau = 2*R/c;
Dfast = ones(Ns,1) * tf - tau' * ones(1, Nf);
phase = pi*Kr*Dfast.^2 - (2 * pi *f0 * tau') * ones(1,Nf);
Sb = Sb+sigma * exp(1j*phase) .* (abs(Dfast) <= Tr/2) .* ((abs(ts -eta_c(ii)) <=Ta/2)' * ones(1,Nf));
end
% Sb = awgn(Sb,SNf,0); % 回波加噪

figure
imagesc(real(Sb)),colormap(gray);

%% CS算法
tf_mtx=ones(Ns,1)*tf;
fr_mtx=ones(Ns,1)*fr;
% --------------------------------------------------------------------
% 变换到距离多普勒域，变标处理
% --------------------------------------------------------------------
D_fa=sqrt(1-lamda^2.*(fa).^2/(4*Vr^2));%(7.17) 徙动因子，列向量
D_fa_mtx=D_fa*ones(1,Nf);% 形成矩阵，大小：Na*Nf
D_fdc=sqrt(1-lamda^2*fdc^2/(4*Vr^2));
K_src_mtx=2*Vr^2*f0^3.*D_fa.^3./(c*R0*(fa).^2)*ones(1,Nf);
Km=Kr./(1-Kr./K_src_mtx);%(7.18)矩阵，这是变换到距离多普勒域的距离调频率。
%变标方程 s_sc
H1=exp(j*pi*Km.*(D_fdc./D_fa_mtx-1).*(tf_mtx-2*R0./(c.*D_fa_mtx)).^2);
% 下面将距离多普勒域的信号与变标方程相乘，实现“补余RCMC”
S1=Srd.*H1;

figure;imagesc(abs(S1));
title('距离多普勒域，补余RCMC后，幅度');

具体代码参考：(39条消息) SAR成像中的ChirpScaling算法；适用小斜视角，并对大斜视成像提供参考。-图像处理文档类资源-CSDN文库

HIT夜枭

关注

2
点赞
踩
46

收藏

觉得还不错? 一键收藏
2
评论
SAR成像系列：【6】合成孔径雷达（SAR）成像算法-调频变标（Chirp Schaling）算法（附Matlab代码）

调频变标（Chirp Schaling）算法通过对线性调频信号进行频率调制，实现尺度变换和平移，并利用相位相乘代替插值运算完成距离徙动补偿，这样不仅提高了算法的计算效率，而且解决二次距离压缩存在的方位频率、依赖问题。
复制链接

扫一扫