SAR图像

合成孔径雷达SAR是一种可成像的雷达，它所用的雷达波段大约是300MHz到30GHz。

比如一般用的波段是1~10GHz的合成孔径雷达，大气对这种波段的影响不大。

也就是说如果天上有一个合成孔径雷达卫星，白天黑夜、大气的云雾雨雪等天气变化对雷达看到的结果影响甚微，可忽略不计。

所以合成孔径雷达是一种全天时、全天候的雷达，它所成的图像就是SAR图像了。

SAR图像的场景和照相机拍出来的场景类似，只不过波段不同看到的事物也不一样。

SAR都是斜视的，而光学的可以垂直照射。

SAR卫星方面，我记得最早发射的是加拿大的Radarsat，且有后续计划。

美国有航天飞机上载的SIR-C等合成孔径雷达。

日本现在有ALOS卫星上载的PALSAR合成孔径雷达（1.27GHz）。

德国的有TerraSAR系列。

据我所知，现在分辨率最高的是德国的X波段SAR系统，数据不好弄到；日本的PALSAR的SAR图像可以到官方网站下载到示例数据。

加拿大的Radarsat和美国的SIR-C数据也是可以到网上下载到的。

机载SAR方面，几乎数的上的雷达强国都有自己的系统。

机载SAR图像有日本的，法国的，德国的，美国的，但是在网络上找这种图像要费点功夫，不是很容易！中国虽然也有，但公开的资料较少，公开的图像资料就更少了。

SAR图像处理软件推荐欧空局的一个免费开源关键PolSARpro（我用过的），可以到欧空局网站上下载，里面的pdf有更详细的介绍SAR及其图像处理等内容。

inSAR技术基于Photoshop插件架构的合成孔径雷达(SAR)图像处理与评估系统主要功能.以图像评估插件的开发为例对关键技术进行了分析.结果表明,采用Photoshop插件方式,可以避免复杂的内存管理编程和用户界面设计,充分利用Photoshop的图形处理功能,减少了工作量,并提高系统稳定性和可用性.所以sar是基于photoshop插件的合成孔径雷达SAR(Synthetic Aperture Radar，合成孔径雷达)问：机载合成孔径雷达sar 多久能生一幅图像就是说采图周期怎么算？最佳答案我们按照最简单的条带式来说，每次生成的是一块图像，这块图像在距离向是全部的范围，在方位向则需要根据你系统的运算能力选择合适的长度。

在SAR中没有所谓采图周期这样的概念，采样是一直运行的，数据处理则在后台运行，把要成像的区域向前后各扩展半个孔径，然后做距离压缩方位压缩即可。

当然这是没有考虑距离徙动校正的情况，如果考虑校正，则要选择算法，一般工程上就是RD的改进和CS。

有什么问题可以联系liuchencalay@我国有机载合成孔径雷达吗？有的话性能如何？其他机载对地雷达如何？能被有效干扰吗？请大虾详解。

最佳答案有，可以确认！性能不详，因为我们搜集到的数据都不完全，已经修改过了，我见美国人一架飞机的数据就是让人怀疑，所以我也对这些所有军事技术数据抱有怀疑。

任何机载雷达的能力都是有限的，毕竟它的动力力量都被缩小了，和地面的大型雷达来比，劣势是很明显的，一旦技术不够顶尖，被干扰是肯定的。

机载合成孔径雷达sar 多久能生一幅图像就是说采图周期怎么算我们在SAR 中没有所谓采图周期这样的概念，采样是一直运行的，数据处理则在后台回答者：27253141|四级| 2010-10-16 06:06合成孔径或者逆向合成孔径成像不过是AESA的一种工作模式而已，高性能的AESA都有这种功能，因此，我们的空警2000就有这种功能。

回答者：greyhond|十四级| 2010-10-16 10:21倒也未必是AESA，APG-68(V)9也加上SAR功能了，B-1B的PESA也一直在用SAR功能。

主要难度是数据的实时处理，和高精度的惯性测量单元。

但是做到高精度，再联动IRST，FLIR和GPS确定目标位置，引导对地弹药就是真正困难的地方了。

无可置疑的是AESA就是绘图也快得多什么事合成孔径雷达，百度百科里面看不懂呀浏览次数：111次悬赏分：0 |解决时间：2011-5-29 19:54 |提问者：我很帅很酷的最佳答案合成孔径雷达就是利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达。

合成孔径雷达的特点是分辨率高，能全天候工作，能有效地识别伪装和穿透掩盖物。

合成孔径雷达主要用于航空测量、航空遥感、卫星海洋观测、航天侦察、图像匹配制导等。

它能发现隐蔽和伪装的目标，如识别伪装的导弹地下发射井、识别云雾笼罩地区的地面目标等。

在导弹图像匹配制导中，采用合成孔径雷达摄图，能使导弹击中隐蔽和伪装的目标。

合成孔径雷达还用于深空探测，例如用合成孔径雷达探测月球、金星的地质结构。

合成孔径雷达工作时按一定的重复频率发、收脉冲，真实天线依次占一虚构线阵天线单元位置。

把这些单元天线接收信号的振幅与相对发射信号的相位叠加起来，便合成一个等效合成孔径天线的接收信号。

若直接把各单元信号矢量相加，则得到非聚焦合成孔径天线信号。

在信号相加之前进行相位校正，使各单元信号同相相加，得到聚焦合成孔径天线信号。

地物的反射波由合成线阵天线接收，与发射载波作相干解调，并按不同距离单元记录在照片上，然后用相干光照射照片便聚焦成像。

这一过程与全息照相相似，差别只是合成线阵天线是一维的，合成孔径雷达只在方位上与全息照相相似，故合成孔径雷达又可称为准微波全息设备。

合成孔径雷达：利用遥感平台的移动，将一个小孔径的天线安装在平台侧方，以代替大孔径的天线，提高方位分辨率的雷达。

合成孔径雷达就是利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达。

合成孔径雷达的特点是分辨率高，能全天候工作，能有效地识别伪装和穿透掩盖物。

合成孔径雷达主要用于航空测量、航空遥感、卫星海洋观测、航天侦察、图像匹配制导等。

它能发现隐蔽和伪装的目标，如识别伪装的导弹地下发射井、识别云雾笼罩地区的地面目标等。

合成孔径雷达工作时按一定的重复频率发、收脉冲，真实天线依次占一虚构线阵天线单元位置。

把这些单元天线接收信号的振幅与相对发射信号的相位叠加起来，便合成一个等效合成孔径天线的接收信号。

若直接把各单元信号矢量相加，则得到非聚焦合成孔径天线信号。

在信号相加之前进行相位校正，使各单元信号同相相加，得到聚焦合成孔径天线信号。

地物的反射波由合成线阵天线接收，与发射载波作相干解调，并按不同距离单元记录在照片上，然后用相干光照射照片便聚焦成像。

这一过程与全息照相相似，差别只是合成线阵天线是一维的，合成孔径雷达只在方位上与全息照相相似，故合成孔径雷达又可称为准微波全息设备。

合成孔径雷达干涉测量有什么软件?浏览次数：315次悬赏分：0 |解决时间：2011-4-21 09:52 |提问者：sdzmchy_l最佳答案商业软件主要有加拿大的EARTHVIEW INSAR软件，瑞士的GAMMA等，开源的软件有doris，ROI_PAC等等。

国内一些研究所也有自行开发的软件合成孔径雷达回波数据处理MATLAB程序浏览次数：458次悬赏分：100 |解决时间：2011-4-9 19:02 |提问者：atnbdr各位朋友，谁能帮忙提供以下合成孔径雷达回波信号处理的MATLAB程序。

本人非常需要，最好是CS算法（Chirp-Scaling）。

没有的话RD算法也可以。

我会追加一倍的分数，非常感谢。

最佳答案CS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%产生Stripmap SAR的回波clear allclcthetaT=0;%T平台波束斜视角thetaT=thetaT*pi/180;%radthetaR=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 domainTr=1.33e-6;%LFM信号脉宽1.33us (200m)Br=150e6; %LFM信号带宽150MHzKr=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 domainv=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 %距离向坐标,方位向坐标,sigmaRrc+50,-50,1Rrc+50,50,1];%%产生回波s_ut=zeros(Nr,Na);U=ones(Nr,1)*u;%扩充为矩阵T=t'*ones(1,Na);for i=1:1:Ntarrn=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;%方位向ffts_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;%距离向fftS_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;%距离向ifftS_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;%显示动态范围40dBG=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;%显示动态范围40dBG=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;%显示动态范围40dBG=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;%显示动态范围40dBG=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 allthetaT=0;%T平台波束斜视角thetaT=thetaT*pi/180;%radthetaR=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 domainTr=1.5e-6;%LFM信号脉宽1.5us (200m)Br=150e6; %LFM信号带宽150MHzKr=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 %距离向坐标,方位向坐标,sigmaRrc+50,-50,1Rrc+50,50,1];%%产生回波s_ut=zeros(Nr,Na);U=ones(Nr,1)*u;%扩充为矩阵T=t'*ones(1,Na);for i=1:1:Ntarrn=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)));%距离向FFTsrc_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).')).';%方位向IFFTsubplot(221)G=20*log10(abs(s_ut)+1e-6);gm=max(max(G));gn=gm-40;%显示动态范围40dBG=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;%显示动态范围40dBG=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;%显示动态范围40dBG=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;%显示动态范围40dBG=255/(gm-gn)*(G-gn).*(G>gn);imagesc(x,r-Rc,G),colormap(gray)grid on,axis tight,xlabel('Azimuth')ylabel('Range')title('(d)目标图象')追问怎么直接复制了一个点目标回波程序阿。