SAR成像(七):RD成像算法
RD算法于1978年处理出第一幅机载SAR数字图像,至今仍在广泛使用,它通过距离和方位上的频域操作,达到了高效的模块化处理要求,同时又具有了一维操作的简便性。该算法根据距离和方位上的大尺度时间差异,在两个一维操作之间使用距离徙动校正(RCMC),对距离和方位进行了近似的分离处理。
由于RCMC是在距离时域-方位频域中实现的,所以也可以进行高效的模块化处理。因为方位频率等同于多普勒频率,所以该处理域又称为“距离多普勒”域。RCMC的“距离多普勒”域实现是RD算法与其他算法的主要区别点,因而称其为距离多普勒算法。
距离相同而方位不同的点目标能量变换到方位频域后,其位置重合,因此频域中的单一目标轨迹校正等效于同一最近斜距处的一组目标轨迹的校正。这是算法的关键,便RCMC能在距离多普勒域高效地实现。
为了提高处理效率,所有的匹配滤波器卷积都通过频域相乘实现,匹配滤波及RCMC都与距离可变参数有关。RD算法区别于其他频域算法的另一王要特点是较易适应距离向参效的变化。所有运算都针对一维数据进行,从而达到了处理的简便和高效。
基本的RD算法流程如下:
在这里插入图片描述
下面给出仿真实验代码
仿真参数:
在这里插入图片描述
close all; clear; clc;
c = 3e8; % 光速
fc = 5e9; %载波频率
B = 200e6; %带宽
lambda = c/fc; %波长
Tp = 1.5e-6; %脉宽
Kr = B/Tp; %调频率
fs = 1.6*B; %采样率
H = 200 ; %飞机高度
Ls = 200 ; %合成孔径长度
v = 100 ; %飞机速度
Lt = Ls/v; %合成孔径时间
%成像区域[X0-Xc,X0+Xc; Y0-Yc,Y0+Yc]
%以合成孔径中心为原点,距离向为x轴,方位向为y轴
Xc = 10000 ;
Yc = 0;
Xo = 100 ;
Yo =100 ;
Rc =sqrt(H2+Xc2); %中心距离
Ka = 2v^2/(Rclambda); %多普勒调频率
Bmax = LtKa; %多普勒最大带宽
fa = ceil(3Bmax); %脉冲重复频率
Rmin = sqrt(H2+(Xc-Xo)2); %观测场景距飞机的最近距离
Rmax = sqrt((Xc+Xo)2+H2+(Yc+Yo+Ls/2)^2); %最远距离
rm = Ls + 2Yo; %雷达走过的总路程长度
tm= 0:1/fa:rm/v-1/fa; %慢时间(合成孔径时间+成像区域时间)
tk= 2Rmin/c-Tp/2:1/fs:2*Rmax/c-1/fs+Tp/2; %快时间(距离门内)
target = [Xc,Yc,0;
Xc+80,Yc+45,0;
Xc-20,Yc-20,0]; %目标坐标
echo = zeros(length™,length(tk),length(target)); %回波
echo_all = zeros(length™,length(tk));
y = -v*(rm/v)/2+vtm; %飞机y轴坐标
for k=1:size(target,1) %目标数
for i=1:length™ %慢时间轴
radar = [0,y(i),H ]; %飞机坐标
Rtm = sqrt(sum((target(k,:)-radar).^2));
echo(i,:,k)=(abs(target(k,2)-y(i))/Xc < 0.01)rectpuls(tk-2Rtm/c,Tp).…
exp(1j2pifc(tk-2Rtm/c)+1jpiKr(tk-2Rtm/c).^2); %回波模型
%注意上边为什么是tk-tr(i)-Tp/2,而不是tk-tr(i)
end
echo_all = echo(:,:,k) + echo_all;
end
r=((tkc/2).2-H2).^(1/2);
%% 距离向匹配滤波
for i = 1:size(echo_all,1)
echo_all(i,:) = echo_all(i,:).(exp(-1j2pifctk)); %回波去载频
end
tt=0:1/fs:Tp-1/fs;
hk=exp(1jpiKrtt.^2); %距离向匹配滤波器
ECHO = zeros(size(echo_all)); %距离向压缩之后
for i=1:length™
ECHO(i,:) = ifft(fft(echo_all(i,:)).*conj(fft(hk,length(tk))));
end
figure;
[R,Y] = meshgrid(r,y);
mesh(R,Y,abs(ECHO));view(0,90);xlim([9900 10100]);
xlabel(‘距离向’);
% ylabel(‘合成孔径 (慢时间), meters’);
title(‘距离压缩后’);
ha = exp(-1jpi2*(v2)/(lambda*Rc)*(0:1/fa:Lt-1/fa).2); %方位向匹配滤波器
for k=1:length(tk)
ECHO(:,k) = fftshift(fft(ECHO(:,k)));
end
%% 距离徙动修正
N=6; %插值核长度为6
ECHO_RCMC=zeros(size(ECHO));
h = waitbar(0,‘Sinc插值中…’); %生成一个进度条
for k = 1:length™
f =linspace(-fa/2,fa/2,length™);
deltaR = (lambdaf(k)/v)^2Rc/8; %注意这个f,到底是哪个f
DU = 2deltaRfs/c;
du = DU-floor(DU);
% kernel_norm = sum(sinc(du-(-N/2:N/2-1)));
for n = N/2+1:length(tk) %快时间
for m = -N/2:N/2-1
if n+floor(DU)+m>length(tk)
ECHO_RCMC(k,n) = ECHO_RCMC(k,n)+ECHO(k,length(tk))*sinc(DU-m);
else
ECHO_RCMC(k,n) = ECHO_RCMC(k,n)+ECHO(k,n+floor(DU)+m)*sinc(du-m);
end
end
end
waitbar(k/length™);
end
close(h); %关闭进度条
figure;
[R,Y] = meshgrid(r,y);
mesh(R,Y,abs(ifft(ECHO_RCMC)));view(0,90);xlim([9900 10100]);
title(‘距离徙动矫正’)
xlabel(‘距离向’);
%% 方位向匹配滤波
for k=1:length(tk)
ECHO(:,k) = abs(ifft(fft(ha,length™)’.*fftshift(ECHO_RCMC(:,k)))); %方位向压缩后
end
figure;
mesh(r,y,ECHO);
xlim([9900 10100]);
xlabel(‘距离向’);
ylabel(‘方位向’);
zlabel(‘幅度’); title(‘最终压缩结果’);
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
点赞
评论
1
分享
收藏
举报
关注
一键三连
点赞Mark关注该博主, 随时了解TA的最新博文
SAR成像仿真-RD算法(Matlab版,含注释)
09-13
合成孔径雷达成像技术是雷达体制中较新的一种,在军事国防领域具有重要的地位,受国内外的广泛的关注。 RD算法是一种基础的SAR成像算法,学生仿真中用的最多,本人将自己做的RD算法整理后分享给大家,希望对
SAR成像RD算法MATLAB仿真
po_huang的博客
9037
SAR成像RD算法MATLAB仿真 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % SAR imaging: RD_algorithm % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%…
SAR成像(七):RD成像算法相关推荐
- 【雷达通信】基于matlab NCP算法SAR回波生成和成像【含Matlab源码 1185期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[雷达通信]基于matlab NCP算法SAR回波生成和成像[含Matlab源码 1185期] 点击上面蓝色字体,直接付费下载,即可. 获 ...
- 单像素成像图像智能处理算法
来源丨 新机器视觉 以前的手机,只在背后有一个安安静静的摄像头,而现在的手机,"二饼" "三筒"甚至"浴霸"层出不穷.越来越强大的手机拍照功 ...
- 三种基于稀疏成像的RCS测量算法思想
1.基于三维稀疏成像的RCS测量方法 算法思路: 1.先设置初始参数: 2.原始回波数据进行脉冲压缩,得到距离向矩阵: 3.对距离向回波数据(脉冲压缩后的数据)进行K倍频域升采样: 3.1 提取距离向 ...
- SAR卫星的不同成像模式
SARclub--小孔经,画山河. SAR卫星的不同成像模式 1.条带模式-strimap 此模式雷达天线可以灵活地调整,改变入射角以获取不同的成像宽 幅.这是基本的SAR成像模式,例如ERS-1RA ...
- 小动物光学成像的发展(荧光成像,生物发光成像)
0 说明 本文内容为Nature上2005的一篇review的简要阅读笔记,属于optical imaging 入门的普及性文章.主要是概念介绍. 1 Basis knowledge 平面成像Phot ...
- java 算法 排序算法_Java七种排序算法以及实现
Java常见七种排序算法以及实现 最近学习一些排序算法,怕自己以后忘记就打算整理起来供自己复习 萌新一枚学习Java没多久,以下仅供参考.如有错误希望大佬指正,欢迎大家在评论区交流探讨. 1.冒泡排序 ...
- 经典算法题每日演练——第七题 KMP算法
原文:经典算法题每日演练--第七题 KMP算法 在大学的时候,应该在数据结构里面都看过kmp算法吧,不知道有多少老师对该算法是一笔带过的,至少我们以前是的, 确实kmp算法还是有点饶人的,如果说红黑树 ...
- 算法设计与分析第七章分支限界算法(完结篇)
算法设计与分析第七章分支限界算法 一.分支界限算法概述 1.分支限界法类似于回溯法,是一种在问题的解空间树上搜索问题解的算法. 分支限界法的求解目标则是找出满足约束条件的一个解,或是在满足约束条件的解 ...
- 第七讲. 经典算法之贪心选择
第七讲. 经典算法之贪心选择 1. 简介 2. 从一个简单例题开始 3. 一个稍难的题目 4. 最重要贪心算法(可作模板) 4.1 最小生成树 4.2 最短路 5. 最后说几句 1. 简介 贪心算法, ...
- 蓝色荧光ps微球/聚苯乙烯蓝色荧光微球用于示踪、体内成像,以及成像仪器和流式细胞仪的校准
蓝色荧光ps微球/聚苯乙烯蓝色荧光微球用于示踪.体内成像,以及成像仪器和流式细胞仪的校准 蓝色荧光微球 聚苯乙烯荧光微球具有荧光强度高,性能稳定,粒径分布窄等特性, 可广泛用于医疗诊断.血流测定.示踪 ...
最新文章
- 新的理念、 新的解决方案、 新的Azure Stack技术预览
- 每天只睡4小时!大佬们都这么拼吗?
- 开发路上踩过的坑要一个个填起来————持续更新······(7月30日)
- mac ox 10.9 安装eclipse cpp launch failed binary not found
- python批量json文件转xml文件脚本(附代码)
- 使用dshow抓取摄像头数据时,回调函数时间为0的问题
- 双向链表删除节点时间复杂度_「十分钟学算法」删除链表的倒数第N个节点
- 斐讯k1路由器刷Breed BootLoader(不死UBoot)
- 关于火狐3,怎么会这样??
- 自我介绍(老师作业,大神自行忽略)
- OTT 端性能优化建设之 Weex 实践之路
- linux gpio-led
- 关于input type=hidden/标签的记录
- uniapp发布苹果IOS测试版 下载无法安装
- 广电为什么禁止投屏_手机投屏到电视机,以前可以投屏,现在不能了,是怎么回事,求解!...
- 版本名称SNAPSHOT、alpha、beta、release、GA含义
- pandas学习之变形
- jQuery boxy弹出层插件中文演示及讲解
- Response响应详解
- 【UOJ311】【UNR #2】积劳成疾
热门文章
- MATLAB 贪婪周边无状态路由协议(GPSR)路由仿真
- android添加快捷按钮,Android中创建快捷方式代码实例
- flume数据采集_大数据采集系统Flume集群部署
- debian linux vnc,Debian 如何配置安装Xfce桌面+VNC远程桌面服务
- php价格结算,PHP根据用户折扣计算商品价格
- linux apache中文名称图片,linux下解决apache不支持中文名的方法详解!
- android stdio 日志追踪,通过检测您的应用生成跟踪日志
- 测试内存兼容软件,Ryzen内存兼容性测试_内存硬盘行情-中关村在线
- 计算机网络信号碰撞问题,计算机网络期末考试试题及问题详解29931(11页)-原创力文档...
- Java中的SPI(Service Provider Interface介绍及示例