【合成孔径雷达】SAR 距离多普勒算法（RD）+ Chirp Scaling算法（CS）+ InSAR干涉SAR人造场景仿真和干涉处理 + PolSAR极化定标【MATLAB】

原创

已于 2022-05-06 10:51:19 修改 · 置顶 · 7.4k 阅读

168 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #信号处理 #数字信号处理

于 2020-06-24 19:03:13 首次发布

最近发现在优快云等很多平台有抄袭本文内容，或本人GitHub源码的行为。

在此特别声明如下：欢迎转载，欢迎分享，这是对我工作的最大肯定。但是请务必标明出处和来源。严禁任何形式的、无说明的、或抄袭成个人内容的行为，甚至以此谋取虚拟或实际奖赏的行为（比如赚取优快云积分）。

如若发现，一定举报。谢谢！

--- 分割线 ---

距离多普勒算法（RD）是SAR的一种经典算法。关于点目标的仿真和分析，以及针对 Radarsat-1 实际数据的成像及处理，详见github：

详细版本，github地址：SAR-Synthetic-Aperture-Radar/1-SAR成像算法 at master · denkywu/SAR-Synthetic-Aperture-Radar · GitHub
简略版，github地址：

最低0.47元/天解锁文章

79 条评论

wsgxp 2024.03.26
博主您好，我想问一下为什么没有使用距离多普勒算法对NSAR场景进行成像呀

m0_51140283 2024.03.02
作者您好，想问下SAR时域回波数据仿真根据INSAR里面的场景进行修改就行吗，另外能求个联系方式吗有些问题想问[face]emoji:010.png[/face]

莱布尼兹2333 2023.05.29
[face]emoji:010.png[/face]楼主可以来个联系方式吗求教

m0_59443584 2022.11.04
博主您好，我想问问在人造场景原始数据仿真中，你觉得有可能在你写的代码基础上把LFM改成步进频连续波（SFCW）吗？我一开始觉得有可能，因为大框架是一样的，但改了一会儿发现好像不行，LFM和SFCW的差别还是太大了[face]emoji:002.png[/face]想问问博主怎么想的

sparks666_ 2022.08.16
您好，rd算法小斜视角的部分，把距离线数从1024改成256剖面图就画不了了，而且点目标成像偏下不在正中间，这是不是坐标轴设置没有将点坐标设置在零点的问题呢，刚入门请指教[face]emoji:056.png[/face][face]emoji:056.png[/face]

嗯，你是对的 2022.07.12
楼主，请问一下SAR成像原始回波数据请问在哪个文件夹能找到呢。我看了看SAR的那个包里面只有点对象的仿真
- denkywu回复嗯，你是对的 2022.07.13
  在【https://github.com/denkywu/SAR-Synthetic-Aperture-Radar/tree/master/2-InSAR%E5%B9%B2%E6%B6%89SAR-%E4%BA%BA%E9%80%A0%E5%9C%BA%E6%99%AF%E4%BB%BF%E7%9C%9F】这里面，InSAR仿真中有平面场景、圆锥场景的仿真，可以参考这个进行原始回波数据仿真。

主观 2022.04.29
您好，请问在rd算法中，进行rcmc距离徙动校正的时候书上说是用对距离向用插值的方式，但是我看代码插值后的信号的行数和列数与插值前一样，不是应该列数会增多吗，是不是我哪里理解错了？
- 主观回复denkywu 2022.04.29
  感谢！
- denkywu回复主观 2022.04.29
  行列数肯定是保持不变的呀。（我残存的一点印象）rcmc距离徙动校正，校正的意思就是把原来“位置放错的，再重新放到正确的位置”。当然这里不是真正的放错，是因为载荷平台运动，回波的斜距必然是不同的，现在要校正到同一水平下进行后处理。所以需要校正，但是整个网格数（行列数不变）。由于计算过程中可能涉及到非整数，所以用的方法是插值。

dhwl1234 2022.01.17
博主您好，我想问问距离多普勒点目标仿真当中，为什么在进行方位轴FFT时候，需要将数据搬移呀，就是乘以一个相位： s_rc_c = s_rc_c.*exp(-1j*2*pi*fnc.*(ta.'*ones(1,N_rg)));
- dhwl1234回复denkywu 2022.01.17
  嗷嗷感谢博主了
- denkywu回复dhwl1234 2022.01.17
  我怕给你带沟里了。因为真的太多年了。我现在都不确定是不是因为先想这么设置频率轴，才有前者。不过前者和后面的频率轴对应是没问题的。而且这里是关键。当初我在方位向频率轴设置这里花了非常多时间，就因为其对应关系一直做不对，成像一直不对。所以这里确实需要好好理解。你再想想。也可以改代码实践下，看下到底有啥区别或影响
- denkywu回复dhwl1234 2022.01.17
  关于频率轴和前面对应关系那里的理解我认为没有问题。不过逻辑上应该是先因为乘以一个多普勒中心的相位，后面的频率轴才如此设置；（我记得，但真不算很确定）而不是先因为要这样设置频率轴，所以才有前者。我感觉首先这里的逻辑应该是这样。但是，为什么前者需要乘以一个多普勒的相位我是真没啥印象了。如果你说书上没有说的话。你看下我的说明文档里有写吗？
- dhwl1234回复denkywu 2022.01.17
  老哥您好，我看了书上似乎并没有说，但是我思考了一下，认为经过距离压缩之后的信号乘以一个多普勒中心的相位，在进行方位向fft时，意味着信号在频域加上一个多普勒中心频率，这样满足方位向频率轴的设置fnc+fftshift（-NFFT_a/2：NFFT_a/2）*Fa/NFFT_a,这里有一个加上多普勒中心频率的操作，不知道我这样理解会对嘛
- denkywu回复dhwl1234 2022.01.17
  啊，时间久远实在是不记得了。你看下书《合成孔径雷达成像算法与实现》，在rd的方位向部分有没有相关说明。

dhwl1234 2021.12.25
老哥，我想问问在zhibao2这个函数中，求积分旁瓣比时候为什么左旁瓣要除以20呀 ISLR = 10*log10((sum(y(1,P1/20:P1))+sum(y(1,P2:end)))/sum(y(1,P1:P2)));% 积分旁瓣比。
- denkywu回复dhwl1234 2021.12.27
  就是大概限定下旁瓣的范围。最左侧从P1的1/20直到P1处，当然也可以把左侧限制的更广，但是没啥必要。我忘了出处是哪儿了，但是没太多意思，就是限定下旁瓣的边界