SAR成像（RMA算法）

RMA 算法叫做距离徙动算法，是一种高分辨率的频域成像算法。 它的核心是通过 STOLT 插值实现无近似的距离方位的解耦合操作。

算法流程图：

 主要的关键在于匹配滤波器的设计和STLOT插值.

二维频域数据经过匹配滤波之后，就相当于得到了粗滤的相位补偿，可以将参考中心附近的点目标相位被很好的聚焦，而远离中心点，聚焦效果越差。所以我们要SOTLOT插值对所有场景中的回波相位进行精准聚焦。SOTLOT插值主要通过对距离频域的映射完成距离徙动的校正以及残余方位压缩。

%创建匹配滤波器
for ii = 1:size(S,2) %遍历每个时间步行以查找 phi_if
    for jj = 1:size(S,1) %逐步遍历当前时间步长行中的每个方位维
        phi_mf(jj,ii) = Rs*sqrt((Kr(ii))^2 - (Kx(jj))^2);
        Krr(jj,ii) = Kr(ii); %生成 2d Kr 用于绘图
        Kxx(jj,ii) = Kx(jj); %生成 2d Kx 用于绘图
    end
end
smf = exp(1i*phi_mf); 

%stolt插值
Ky_even = linspace(kstart, kstop, 1024); %为真实数据的插值创建均匀间隔的 Ky

count = 0;
for ii = 1:zpad
    count = count + 1;
    Ky(count,:) = sqrt(Kr.^2 - Kx(ii)^2);
    S_st(count,:) = (interp1(Ky(count,:), S_mf(ii,:), Ky_even));
end

 

具体的公式推导有很多论文，这里就不在推导，下面拿MIT的课程实验代码学习：

引言：

学生将以 3 人一组的方式构建雷达，并在雷达套件实现时进行一系列实验。 您将把雷达套件带到现场并执行其他实验，例如测量经过的汽车的速度或绘制移动目标的范围。 最终的 SAR 成像竞赛将测试您从校园周围选择的目标场景形成 SAR 图像的能力； 最详细和最具创意的图像获胜。

 

雷达特征：

• 相干 FMCW 架构
• S波段
• 使用咖啡罐作为发射和接收天线
• 6个微型电路元件
• 1 个四运算放大器作为视频增益级和抗混叠滤波器
• 带触发输出的模拟斜坡发生器
• 连接到计算机的音频输入以进行视频数字化和触发
• 记录您的实验的 .wav 文件
• MATLAB® 脚本读取 .wav 文件并整理触发脉冲和脉冲组以处理 3 种操作模式
• 3 种操作模式包括： 多普勒与时间、距离与时间、合成孔径雷达成像

效果如下所示：

该代码的流程：方位向FFT-距离向FFT-匹配滤波-STOLT插值-2D-IFFT

 方位向FFT:

 方位距离2D-FFT:

匹配滤波之后：

 STOLT插值后2D-IFFT:

 这是其它数据的成像效果图：

（代码和数据我经过修改后上传,包括代码数据实验报告）下载链接：

https://download.csdn.net/download/weixin_44670677/87952554

PS:引用链接：
Projects | Build a Small Radar System Capable of Sensing Range, Doppler, and Synthetic Aperture Radar Imaging | Supplemental Resources | MIT OpenCourseWare