合成孔径雷达回波数据处理MATLAB程序
作者&投稿：帛露 （若有异议请与网页底部的电邮联系）
以前是有，不知道放哪儿了，那代码还有一点毛病。
弄这个研究，书上都是公式，没人能教你怎么写代码，怎么弄出图像。
跟你说下心得吧，就从怎么读出数据开始一步一步走。自己弄才是硬道理。
学习编程语言，包括你说的API和MFC这些用于做界面的软件，都是为了你以后能更好的从事遥感领域，遥感现在是从研究阶段开始进入实用阶段，这样我们就需要根据成熟的理论，把一些想法变成现实，而遥感领域很多东西都没做出来，或者说做的都不是很成熟，也可以说别人做的一些软件什么的并不能很好的完成你的想法，或者使用上有局限性，这时如果你在做一些大的项目你很可能就需要自己去编程实现，当然也不是你说的用VC编写界面，那麻烦了去了，可以用Qt，API，MFC等等很多在编写界面领域更实用的软件，但是很多甚至绝大部分都已c 为基础，学了c，再学那些可能学的会快一些，C语言仍然是也仍将是全世界最为成功的编程语言之一
CS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%产生Stripmap SAR的回波
clear all
clc
thetaT=0;%T平台波束斜视角
thetaT=thetaT*pi/180;%rad
thetaR=0;%R平台波束斜视角
thetaR=thetaR*pi/180;
c=3e8;%光速
fc=1.5e9;%载频
lambda=c/fc;%波长
%%测绘带区域
X0=200;%方位向[-X0,X0]
Rtc=3000;
Rrc=3000;
Rc=(Rtc+Rrc)/2;
R0=150;%距离向[Rc-R0,Rc+R0]
%%距离向(Range),r/t domain
Tr=1.33e-6;%LFM信号脉宽1.33us (200m)
Br=150e6; %LFM信号带宽 150MHz
Kr=Br/Tr; %调频斜率
Nr=1024;
r=Rc+linspace(-R0,R0,Nr);
t=2*r/c;%t域序列
dt=R0*4/c/Nr;%采样周期
f=linspace(-1/2/dt,1/2/dt,Nr);%f域序列
%%方位向(Azimuth,Cross-Range),x/u domain
v=100;%SAR 平台速度
Lsar=300;%合成孔径长度
Na=512;
x=linspace(-X0,X0,Na);%u域序列
u=x/v;
du=2*X0/v/Na;
fu=linspace(-1/2/du,1/2/du,Na);%fu域序列
ftdc=v*sin(thetaT);
ftdr=-(v*cos(thetaT))^2/lambda/Rtc;
frdc=v*sin(thetaR);
frdr=-(v*cos(thetaR))^2/lambda/Rrc;
fdc=ftdc+frdc;%Doppler调频中心频率
fdr=ftdr+frdr;%Doppler调频斜率
%％目标位置
Ntar=3;%目标个数
Ptar=[Rrc,0,1 %距离向坐标,方位向坐标,sigma
Rrc+50,-50,1
Rrc+50,50,1];
%%产生回波
s_ut=zeros(Nr,Na);
U=ones(Nr,1)*u;%扩充为矩阵
T=t'*ones(1,Na);
for i=1:1:Ntar
rn=Ptar(i,1);xn=Ptar(i,2);sigma=Ptar(i,3);
rtn=rn+Rtc-Rrc;
RT=sqrt(rtn^2+(rtn*tan(thetaT)+xn-v*U).^2);
RR=sqrt(rn^2+(rn*tan(thetaT)+xn-v*U).^2);
R=RT+RR;
DT=T-R/c;
phase=-pi*Kr*DT.^2-2*pi/lambda*R;
s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)<Tr/2).*(abs(v*U-xn)<Lsar/2);
end;
%方位向fft
s_kt=fftshift(fft(fftshift(s_ut).')).';
%CS变换
kc=4*pi/lambda;
kc=kc*ones(1,Na);
kx=fu/v;
p_kx0=-sqrt(kc.^2-kx.^2);%相位项泰勒展开的系数函数
p_kx1=2*kc/c/p_kx0;
p_kx2=-2.*kx.^2/c^2./p_kx0.^3;
C_kx=-(c*p_kx1/2+1);
Ks_r=1-2*Kr*Rc.*p_kx2;
Ks_kx_r=Kr/pi./Ks_r;
r0=Rc;
s2_ut=exp(j*pi*C_kx.*ones(Nr,1)*Ks_kx_r.*(t'*ones(1,Na)-2*r0*(1+C_kx)/c).^2);%设计的线性调频信号
S_cs=s_kt.*s2_ut;
%距离向fft
S_kw=fftshift(fft(fftshift(S_cs)));
%距离向匹配滤波
w=2*pi*f;
rmc_r=exp(j.*w*2*C_kx*r0/c).*exp(j.*w.^2/4/pi/Kr/(1+C_kx));
rmc_r=rmc_r'*ones(1,Na);
S_rmc=S_kw.*rmc_r;
%距离向ifft
S_kt=fftshift(ifft(fftshift(S_rmc)));
d_kxr=4*pi/c^2*Kr*C_kx*(1+C_kx).*(Rc-r0).^2;%CS变换带来的相位误差
S_kt=S_kt.*exp(-j*d_kxr);%消除相位误差
%方位向匹配滤波
FU=ones(Nr,1)*fu;
H_kx=exp(j*pi/fdr*(FU-fdc).^2);%方位向压缩因子
I_ut=S_kt.*H_kx;
I_ut=fftshift(ifft(fftshift(I_ut.'))).';
subplot(221)
G=20*log10(abs(s_ut)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,-G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(a)原始信号')
subplot(222)
G=20*log10(abs(S_rmc)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,-G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(b)距离向匹配滤波后频谱')
subplot(223)
G=20*log10(abs(S_kt)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(c)消除相位误差后频谱')
subplot(224)
G=20*log10(abs(I_ut)+1e-6);
gm=max(max(G));
gn=gm-60;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(d)目标图象')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
RD
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%产生Stripmap SAR的回波
clear all
thetaT=0;%T平台波束斜视角
thetaT=thetaT*pi/180;%rad
thetaR=0;%R平台波束斜视角
thetaR=thetaR*pi/180;
c=3e8;%光速
fc=1.5e9;%载频
lambda=c/fc;%波长
%%测绘带区域
X0=200;%方位向[-X0,X0]
Rtc=3000;
Rrc=3000;
Rc=(Rtc+Rrc)/2;
R0=150;%距离向[Rc-R0,Rc+R0]
%%距离向(Range),r/t domain
Tr=1.5e-6;%LFM信号脉宽 1.5us (200m)
Br=150e6; %LFM信号带宽 150MHz
Kr=Br/Tr; %调频斜率
Nr=512;
r=Rc+linspace(-R0,R0,Nr);
t=2*r/c;%t域序列
dt=R0*4/c/Nr;%采样周期
f=linspace(-1/2/dt,1/2/dt,Nr);%f域序列
%%方位向(Azimuth,Cross-Range),x/u domain
v=100;%SAR 平台速度
Lsar=300;%合成孔径长度
Na=1024;
x=linspace(-X0,X0,Na);%u域序列
u=x/v;
du=2*X0/v/Na;
fu=linspace(-1/2/du,1/2/du,Na);%fu域序列
ftdc=v*sin(thetaT);
ftdr=-(v*cos(thetaT))^2/lambda/Rtc;
frdc=v*sin(thetaR);
frdr=-(v*cos(thetaR))^2/lambda/Rrc;
fdc=ftdc+frdc;%Doppler调频中心频率
fdr=ftdr+frdr;%Doppler调频斜率
%％目标位置
Ntar=3;%目标个数
Ptar=[Rrc,0,1 %距离向坐标,方位向坐标,sigma
Rrc+50,-50,1
Rrc+50,50,1];
%%产生回波
s_ut=zeros(Nr,Na);
U=ones(Nr,1)*u;%扩充为矩阵
T=t'*ones(1,Na);
for i=1:1:Ntar
rn=Ptar(i,1);xn=Ptar(i,2);sigma=Ptar(i,3);
rtn=rn+Rtc-Rrc;
RT=sqrt(rtn^2+(rtn*tan(thetaT)+xn-v*U).^2);
RR=sqrt(rn^2+(rn*tan(thetaT)+xn-v*U).^2);
R=RT+RR;
DT=T-R/c;
phase=pi*Kr*DT.^2-2*pi/lambda*R;
s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)<Tr/2).*(abs(v*U-xn)<Lsar/2);
end;
%%距离压缩
p0_t=exp(j*pi*Kr*(t-2*Rc/c).^2).*(abs(t-2*Rc/c)<Tr/2);%距离向LFM信号
p0_f=fftshift(fft(fftshift(p0_t)));
s_uf=fftshift(fft(fftshift(s_ut)));%距离向FFT
src_uf=s_uf.*(conj(p0_f).'*ones(1,Na));%距离压缩
src_ut=fftshift(ifft(fftshift(src_uf)));%距离压缩后的信号
src_fut=fftshift(fft(fftshift(src_ut).')).';%距离多普勒域
%%二次距离压缩,距离迁移校正原理仿真
src_fuf=fftshift(fft(fftshift(src_uf).')).';%距离压缩后的二维频谱
F=f'*ones(1,Na);%扩充为矩阵
FU=ones(Nr,1)*fu;
p0_2f=exp(j*pi/fc^2/fdr*(FU.*F).^2+j*pi*fdc^2/fc/fdr*F-j*pi/fc/fdr*FU.^2.*F);
s2rc_fuf=src_fuf.*p0_2f;
s2rc_fut=fftshift(ifft(fftshift(s2rc_fuf)));%距离多普勒域
%%方位压缩
p0_2fu=exp(j*pi/fdr*(FU-fdc).^2);%方位向压缩因子
s2rcac_fut=s2rc_fut.*p0_2fu;%方位压缩
s2rcac_fuf=fftshift(fft(fftshift(s2rcac_fut)));%距离方位压缩后的二维频谱
s2rcac_ut=fftshift(ifft(fftshift(s2rcac_fut).')).';%方位向IFFT
subplot(221)
G=20*log10(abs(s_ut)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,-G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(a)原始信号')
subplot(222)
G=20*log10(abs(src_fut)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(fu,r-Rc,-G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(b)距离多普勒域频谱')
subplot(223)
G=20*log10(abs(s2rc_fut)+1e-6);
gm=max(max(G));
gn=gm-40;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(fu,r-Rc,-G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(c)RMC后的RD域频谱')
subplot(224)
G=20*log10(abs(s2rcac_ut)+1e-6);
gm=max(max(G));
gn=gm-60;%显示动态范围40dB
G=255/(gm-gn)*(G-gn).*(G>gn);
imagesc(x,r-Rc,G),colormap(gray)
grid on,axis tight,
xlabel('Azimuth')
ylabel('Range')
title('(d)目标图象')
合成孔径雷达回波数据处理MATLAB程序
v=100;%SAR 平台速度 Lsar=300;%合成孔径长度 Na=512;x=linspace(-X0,X0,Na);%u域序列 u=x\/v;du=2*X0\/v\/Na;fu=linspace(-1\/2\/du,1\/2\/du,Na);%fu域序列 ftdc=v*sin(thetaT);ftdr=-(v*cos(thetaT))^2\/lambda\/Rtc;frdc=v*sin(thetaR);frdr=-(v*cos(thetaR))^2\/lam...
关于合成孔径雷达
详情请查看视频回答
合成孔径雷达采集的是什么数据?
图像数据，回波形成图像
合成孔径雷达干涉(InSAR)应用
合成孔径雷达干涉技术，简称InSAR，是一项利用雷达回波相位信息的精密测量手段。这项技术具有显著优势，能够构建高精度、覆盖大面积的地理数据，并且借助其独特的差分干涉原理，能精确监测地面微小的毫米级位移变化。其应用领域广泛，对于天文地球动力学研究等领域具有重要作用。InSAR的核心原理是将来自不同轨道位...
合成孔径雷达的原理?
合成孔径雷达利用一个小天线沿着长线阵的轨迹等速移动并辐射相参信号, 把在不同位置接收的回波进行相干处理, 从而获得较高分辨率的成像雷达，与其它大多数雷达一样，合成孔径雷达通过发射电磁脉冲和接收目标回波之间的时间差测定距离，其分辨率与脉冲宽度或脉冲持续时间有关，脉宽越窄分辨率越高。合成孔径雷达...
遥感雷达有多厉害?
合成孔径雷达干涉测量技术（INSAR，Interferometric Synthetic Aperture Radar；简称：干涉雷达测量）是以同一地区的两张SAR图像为基本处理数据，通过求取两幅SAR图像的相位差，获取干涉图像，然后经相位解缠，从干涉条纹中获取地形高程数据的空间对地观测新技术。差分干涉雷达测量技术（D-INSAR）是指利用同一地区...
合成孔径雷达深度学习成像研究综述
深度学习驱动的合成孔径雷达成像：一场革命性的探索 SAR成像领域的创新者张群等人引领了一场深度学习的革命，将这一技术应用于各种场景，如静止、运动、三维和ISAR成像。通过深度学习网络的架构，如卷积神经网络（CNN）、循环神经网络（RNN）和生成对抗网络（GAN），SAR图像处理的精度和效率得到了前所未有的...
光学图像与合成孔径雷达图像的区别
成像机制差别：光学影像通常采用中心投影面域成像或推帚式扫描获取数据；而SAR处于信号处理的需要（合成孔径过程，这里就不展开讨论了）不能采用垂直向下的照射方式而只能通过测视主动成像方式发射和接受面域雷达波，并通过信号处理（聚焦、压缩、滤波等）手段后期合成对应于地面目标的复数像元。3、在图像分割上...
合成孔径雷达干涉(InSAR)测量模式
合成孔径雷达干涉（InSAR）是一种精密的地球观测技术，主要依赖于两种不同的测量模式来获取高精度的数据。首先，双天线单航过模式是常见的一种，它通常在航天飞机等飞行平台上实现。这种模式包括一个发射和接收的天线，以及一个专门用于接收信号的辅助天线。发射天线发出的雷达波束经过目标区域后，由接收天线...
合成孔径雷达图像精准处理目录
合成孔径雷达图像精准处理1. SAR图像原理和发展 合成孔径雷达(SAR)的基本原理涉及其成像过程和应用特性。SAR利用卫星在飞行中的运动，通过回波信号分析地面物体信息。其图像特点包括全天候、全天时和高分辨率。当前，星载SAR系统如日本PALSAR、意大利COSMO SAR、德国TerraSAR-X、加拿大RadarSat-2等都有显著发展...
叔王泰乐： CS%产生Stripmap SAR的回波 clear all clc thetaT=0;%T平台波束斜视角 thetaT=thetaT*pi/180;%rad thetaR=0;%R平台波束斜视角 t...
务川仡佬族苗族自治县15086602376： 关于合成孔径雷达的简单问题 - ？
叔王泰乐： 合成孔径雷达就是利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达.合成孔径雷达的特点是分辨率高,能全天候工作,能有效地识别伪装和穿透掩盖物.合成孔径雷达主要用于航空测...
务川仡佬族苗族自治县15086602376： 谁熟悉INSAR?？
叔王泰乐： 合成孔径雷达干涉测量技术(INSAR,Interferometric Synthetic Aperture Radar;简称:干涉雷达测量)是以同一地区的两张SAR图像为基本处理数据,通过求取两幅SAR图像的相位差,获取干涉图像,然后经相位解缠,从干涉条纹中获取地形高程数据的空间对地观测新技术.
务川仡佬族苗族自治县15086602376： 什么是合成孔径雷达? ？
叔王泰乐： 合成孔径雷达与微波成像 微波成像属于遥感技术的一种,利用雷达等对地面目标进行扫描,通过计算回波特性来得到地面目标的物理特性,能够实现这种功能的雷达叫做合成孔径雷达(SAR).合成孔径雷达在军事上的最大作用就是可以全天候监控地面目标,可用于侦查、导航、精确打击等多种任务.目前先进的飞机雷达均有合成孔径能力,很多作战飞机也装备了SAR吊舱,从而具备精确打击能力.近年来还发展出了专门的SAR飞机,例如美国空军的E-8“联合星”空中指挥机,就拥有强大的合成孔径能力,可指挥作战飞机对预定地面目标进行打击.
务川仡佬族苗族自治县15086602376： 求合成孔径雷达成像的原理~~越详细越通俗越好,最好给点儿资料网站 ？
叔王泰乐： 一种高方位分辨率的相干成象雷达.可分为侧视、斜视、多普勒锐化和聚束测绘等工作方式.利用合成的天线技术获取良好的方位分辨率,利用脉冲压缩技术获取良好的距...
务川仡佬族苗族自治县15086602376： 什以是合成孔径雷达?..合成孔径雷达是什么样的,它的工作原理是什 ？
叔王泰乐： 雷达问题我喜欢!我来说两句!记住要选答案啊,呵呵! 合成孔径雷达简称SAR,是新计算机模拟技术支持下诞生的一种新型遥感设备. 为了方便理解,首先简单说说普...
务川仡佬族苗族自治县15086602376： 合成孔径雷达的介绍 - ？
叔王泰乐： 合成孔径雷达( SAR) 是一种高分辨率成像雷达,可以在能见度极低的气象条件下得到类似光学照相的高分辨雷达图像.利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达,也称综合孔径雷达.合成孔径雷达的特点是分辨率高,能全天候工作,能有效地识别伪装和穿透掩盖物.所得到的高方位分辨力相当于一个大孔径天线所能提供的方位分辨力.合成孔径雷达的首次使用是在20世纪50 年代后期,装载在RB-47A和RB-57D 战略侦察飞机上.经过近60 年的发展,合成孔径雷达技术已经比较成熟,各国都建立了自己的合成孔径雷达发展计划,各种新型体制合成孔径雷达应运而生,在民用与军用领域发挥重要作用.
务川仡佬族苗族自治县15086602376： 脉冲多普勒雷达、合成孔径雷达、相控阵雷达三者有什么区别? - ？
叔王泰乐： 我们知道,蜻蜓的每只眼睛由许许多多个小眼组成,每个小眼都能成完整的像,这样就使得蜻蜓所看到的范围要比人眼大得多.与此类似,相控阵雷达的天线阵面也由许多个辐射单元和接收单元(称为阵元)组成,单元数目和雷达的功能有关,...
务川仡佬族苗族自治县15086602376： 合成孔径雷达问题 - ？
叔王泰乐： 主要目的是解除距离徒动所引起的距离向和方位向的耦合,完成方位向聚焦处理,我的理解,希望对你有帮助
务川仡佬族苗族自治县15086602376： 介绍以下脉冲多普勒雷达和合成孔径雷达工作原理,谢谢! ？
叔王泰乐： 利用信号频域特性分辨和检测目标的脉冲雷达.目标和干扰物相对于雷达的径向速度不同,回波信号也有不同的多普勒频率.可用频域过滤的方法选出目标的多普勒频率谱...
