SAR合成孔径雷达图像点目标仿真报告附matlab代码.docx
- 文档编号:11337512
- 上传时间:2023-02-28
- 格式:DOCX
- 页数:25
- 大小:84.36KB
SAR合成孔径雷达图像点目标仿真报告附matlab代码.docx
《SAR合成孔径雷达图像点目标仿真报告附matlab代码.docx》由会员分享,可在线阅读,更多相关《SAR合成孔径雷达图像点目标仿真报告附matlab代码.docx(25页珍藏版)》请在冰豆网上搜索。
SAR合成孔径雷达图像点目标仿真报告附matlab代码
SAR合成孔径雷达图像点目标仿真报告(附
matlab代码)
SAR图像点目标仿真报告
徐一凡
1SAR原理简介
合成孔径雷达(SyntheticApertureRadar,简称SAR)是一种高分辨率成像雷达技术。
它利用脉冲压缩技术获得高的距离向分辨率,利用合成孔径原理获得高的方位向分辨率,从而获得大面积高分辨率雷达图像。
SAR回波信号经距离向脉冲压缩后,雷达的距离分辨率雷达发射信号带宽决定:
rC,式中r表示雷达的距离分辨率,Br表示雷达发射信号带宽,C表示光速。
同2Br样,SAR回波信号经方位向合成孔径后,雷达的方位分辨率雷达方位向的多谱勒带宽决定:
ava,式中a表示雷达的方位分辨率,Ba表示雷达方位
向多谱勒带宽,va表示BaD,其中D为2方位向SAR平台速度。
在小斜视角的情况下,方位分辨率近似表示为a方位向
合成孔径的长度。
2SAR的几何关系
雷达位置和波束在地面覆盖区域的简单几何模型如图1
所示。
此次仿真考虑的是正侧视的条带式仿真,也就是说倾斜角为零,SAR波束中心和SAR平台运动方向垂直的情况。
图1雷达数据获取的几何关系
建立坐标系XYZ如图2所示,其中X0丫平面为地平面;
SAR平台距地平面高H,以速度V沿X轴正向匀速飞行;P点为SAR平台的位置矢量,设其坐标为(x,y,z);T点为目标
的位置矢量,设其坐标为(xT,yT,zT);几何关系,目标与SAR平台的斜距为:
RPT(xxT)2(yyT)2(zzT)2
(1)图可知:
yO,zH,zTO;令x
其中v为平台速度,s为慢时间变量,vs。
假设xTvs,其中s表示SAR平台的x坐标为xT的时刻;再令r标与SAR的垂直斜距,重写⑴式为:
H2yT2,r表示目
PTR(s;r)r2v2(ssO)2
(2)
R(s;r)就表示任意时刻s时,目标与雷达的斜距。
一般
情况下,vssOr,于是通过
傅里叶技术展开,可将
(2)式可近似写为:
v2R(s;r)rv(ssO)r(ssO)2(3)
2r222可见,斜距是s和r的函数,不同的目标,r也不一样,但当目标距SAR较远时,在观测带内,可近似认为r不变,即rRO。
图2:
空间几何关系(a)正视图(b)侧视图
图2(a)中,Lsar表示合成孔径长度,它和合成孔径时间Tsar的关系是LsarvTsar。
(b)中,为雷达天线半功率点波束角,为波束轴线与Z轴的夹角,即波束视角,Rmin为近距点距离,Rmax为远距点距离,W为测绘带宽度,它们的关系为:
RminHtg
(2)RmaxHtg
(2)(4)WRmaxRminSAR的回波信号模型
SAR在运动中以一定的周期(1/PRF)发射和接收信号,具体过程如图3所示。
发射机
以I的时间发射啁啾脉冲,然后切换天线开关接收回波信号。
图3雷达发射脉冲串的时序
当雷达不处于发射状态时,它接收3反射回波。
发射和
接收回波的时间序列如图4所示。
在机载情况下,每个回波可以在脉冲发射间隔内直接接收到。
但是在星载情况下,于距离过大,某个脉冲的回波要经过6〜10个脉冲间隔才能接
收到。
这里仿真为了方便,默认为机载情况。
图4脉冲雷达的发射与接收周期
假设Tr为chirp信号持续时间,下标r表示距离向;PRF为重复频率,PRT为重复周期。
等于1/PRFo接收序列中,n2*R(s;r)表示发射第i个脉冲时,目标回波相对于发射C序列的延时。
雷达的发射序列
数学表达式为式(5):
s(t)
叩(tn*PRT)
2tp(t)rectejKrtej2fctTr
)表示矩形信号,Kr为距离向的chirp信号调频率,fc为载频。
式中,rect(雷达回波信号发射信号波形,天线
方向图,斜距,目标RCS环境等因素共同决定。
若不考虑环境因素,则单点目标雷达回波信号可写成式
⑹所示:
sr(t)
nwp(tnPRTn)(6)
其中,表示点目标的雷达散射截面,w表示点目标天线
方向图双向幅度加权,n表
示载机发射第n个脉冲时,电磁波再次回到载机时的延
时n2*R(s;r),带入式(6)中得:
Csr(t)nwrect(4tnPRT2R(s;r)/C)Tr(7)
exp[jKr(tnPRT2R(s;r)/C)2]exp[-j
R(s;r)]exp[j2fc(tnPRTn)]式⑺就是单点目标回波信号模型,其中,exp[jKr(tnPRT2R(s;r)/C)2]是
chirp分量,它决定距离向分辨率;exp[-j
4它决定方位向分辨率。
R(s;r)]为多普勒分量。
对于任意一个脉冲,回波信号可表示为式(8)所示:
sr(t,s)A0wr(2R(s;r)/C)wa(ssc)exp{j4f0R(s;r)/C}
exp{jKr(2R(s;r)/C)}2
我们知道,于R(s;r)随慢时间s的变化而变化,所以计算机记录到的回波数据存储
形式如图5所示:
图5目标照射时间内,单个点目标回波能量在信号处理器的二维存储器中的轨迹
4距离徙动及校正
根据图2可知,在倾斜角为零或很小的时候,目标与雷
达的瞬时距离为R(s;r),根据
几何关系可知,R(s;r)r2v2(ss0)2,根据泰勒级数展开
可得:
222v2R(s;r)rv(ss0)r(ss0)2(9)
2r
式(9)可知,不同慢时间对应着不同的R(s;r),并且是
一个双曲线形式或者近似为一
个二次形式。
如图5所示,同一目标的回波存储在计算
机里不在同一直线上,存在距离徙动。
从而定义距离徙动量:
v2R(s,r)(ss0)2
2r为了进行方位向的压缩,方位向的回波数据必须在同
一条直线上,也就是说必须校正距离徙动R(s,r)。
式可知,
不同的最近距离r对应着不同的R(s,r),因此在时域处理距离徙动会非常麻烦。
因此,对方位向进行傅里叶变换,对距离向不进行变换,得到新的域。
于方位向的频率即为多普勒频率,所以这个新的域也称为距离多普勒域。
将斜距R写成多普勒fa的函数,即R(fa,r)。
众所周知,对最近距离为r的点目标P回波多普勒fa是倾斜角的函数,即fa2Vsin,斜距R(fa,r)r/cos,于是
R(fa,r)r/cosr/1sin2r/1(fa2)2V(11)
1r2rfa28V所以距离多普勒域中的我距离徙动为
R(fa,r)=2rfa2,可发现它不随慢时间变换,同一最短距离
r对应着相同大小的距离徙动。
因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正,这样可以节省运算量。
为了对距离徙动进行
校正,需要得到距离徙动单元,即距离徙动体现在存储单元中的移动数值,距离徙动单元可以表示为R(fa,r)/r,这个
值通常为一个分数,于存储单元都是离散的,所以不同通过在存储单元简单的移动得到准确的值。
为了得到准确的徙动校正值,通常需要进行插值运算。
本仿真采用了两种插值方法最近邻点插值和sinc插值,
下面分别进行介绍。
最近邻点插值法的优点是简单而快速,缺点是不够精确。
其中N为整数部分,R(fa,r)/r=Nn,n为
小数部分,整数部分徙动可以直接通过平移消除,对于小数部分则通过四舍五入的方法变为0或者1,这样就可以得到
较为精确的插值。
Sine插值原理如下:
在基带信号下,卷积核是sine函
数
18Vh(x)sinc(x)插值信号为
sin(x)x
g(x)gd(i)sinc(xi)
i
即为所有输入样本的加权平均。
可通过频域来理解,如图6所示,采样信号gd(i)的频
谱Gd(f)等于以采样率重复的信
号频谱。
为了重建信号g(x),只需要一个周期频谱,因此需要理想矩形低通滤波器在频域中提取基带频谱所示。
已知该理想滤波器在时域中是sine函数。
于频域相乘相当于
时域卷积,故插值可以通过与sine核的卷积来实现。
图6理想低通滤波器怎样对采样信号进行插值
5点目标成像matlab仿真
距离多普勒算法距离多普勒算法是在1976年至1978
年为民用星载SAR提出的,它兼顾了成熟、简单、高效和精确等因素,至今仍是使用最广泛的成像算法。
它通过距离和方位上的频域操作,到达了高效的模块化处理要求,同时又具有了一维操作的简便性。
图7示意了RDA的处理流程。
这里主要讨论小倾斜角及短孔径下的基本RDA处理框图。
1.
当数据处在方位时域时,可通过快速卷积进行距离压缩。
也
就是说,距离FFT后随即进行距离向匹配滤波,再利用距离IFFT完成距离压缩。
回波信号为:
s0(t,s)A0wr[t2R(s)/c]wa(ssc)exp{-j4fOR(s)/c}exp{jKr(t-2R(s)/c)}距离向压缩后的信
号为:
2
src(t,s)IFFTt{SO(ft,s)H(ft)}
A0r[t2R(s)/c]wa(ssc)exp{j4f0R(s)/c}ff2H(ft)rect{}exp{j}exp{j2ft0}
|K|TK2.通过方位FFT将数据变换至距离多普勒域,多普勒中心频率估计以及大部分后续操作
都在该域进行。
方位向傅里叶变换后信号为:
S1(t,fs)FFTs{src(t,s)}2Rrd(fs)]Wa(fsfsc)
c4f0R0fs2exp{-j}exp{j}cKaA0pr[t3.在距离多普
勒域进行随距离时间及方位频率变化的RCMC该域中同一距
离上的一组
目标轨迹相互重合。
RCM(将距离徙动曲线拉直到与方位
频率轴平行的方向。
这里可以采用最近邻点插值法或者sinc
插值法,具体插值方法见前面。
假设RCMC插值是精确的,
信号变为:
2R0)Wa(fsfsc)c
4f0R0fs2exp{-j}exp{j}cKaS2(t,fs)A0pr(t
4.通过每一距离门上的频域匹配滤波实现方位压缩。
为进行方位压缩,将RCMC后的
S2(t,fs)乘以频域匹配滤波器Haz(fs)。
fs2Haz(fs)exp{j}
KaS3(t,fs)S2(t,fs)Haz(fs)
A0pr(t2R0/c)Wa(fsfsc)exp{j4f0R0
}c5.最后通过方位IFFT将数据变换回时域,得到压缩后的复图像。
复原后的图像为:
sac(t,s)IFFTs{S3(t,fs)}A0pr(t-2R0/c)pa(s)
exp{-j4f0R0}exp{j2fscs}c
图8距离多普勒算法流程图
ChirpScaling算法距离多普勒算法具有诸多优点,但是距离多普勒算法有两点不足:
首先,当用较长的核函数提高距离徙动校正精度时,运算量较大;其次,二次距离压缩对方位频率的依赖性问题较难解决,从而限制了其对某些大斜视角和长孔径SAR的处理精度。
ChirpScaling算法避免了RCMC中的插值操作,通过对Chirp信号进行频率调制,实现
了对该信号的尺度变换或平移。
图8显示了Chirp
Scaling算法处理流程。
这里主要讨论小倾斜角及短孔径下的基本CSA处理框图。
主要步骤包括四次FFT和三次相位相
乘。
1.通过方位向FFT将数据变换到距离多普勒域。
2.
通过相位相乘实现ChirpScaling操作,使所有目标的距离
徙动轨迹一致化。
这是第一步相位相乘。
用以改变线调频率尺度的ChirpScaling二次相位函数为:
H1(t,fa;Rs)exp[j(fa;RB)a(fa)(t2R(fa;Rs)2)]c
3.通过距离向FFT将数据变到二维频域。
4.通过与参
考函数进行相位相乘,同时完成距离压缩、SRC和一致RCMC这是第二步相位相乘。
用于距离压缩,距离徙动校正的相位函数写为:
H2(fr,fa;Rs)exp[j
1fr2](fa;RB)[1a(fa)]exp[j4Rsa(fa)]fr]c
5.通过距离向IFFT将数据变回到距离多普勒域。
6.
通过与随距离变化的匹配滤波器进行相位相乘,实现方位压缩。
此外,于步骤2中的ChirpScaling操作,相位相乘中还需要附加一项相位校正。
这是第三步相位相乘。
补偿ChirpScaling弓I起的剩余相位函数是:
H2(tr,fa;RB)exp[j
2RBVfaM2fa2]exp[j(fa;RB)]
7.最后通过方位向IFFT将数据变回到二维时域,即SAR
图像域。
图8ChirpScaling算法流程图
简而言之,R-D算法是将徙动曲线逐一校正,CS算法是
以某一徙动曲线为参考,在Doppler域内消除不同距离门的
徙动曲线的差异,令这些曲线成为一组相互“平行”的曲线,
然后在二维频率域内统一的去掉距离徙动。
通俗一点就是,RD算法是将弯曲的信号一根根掰直,而CS算法是先把所有
信号都掰得一样弯,然后再统一掰直。
6仿真结果
使用最近邻点插值的距离多普勒算法仿真结果首先
对5个点目标的回波信号进行了仿真,5个点目标构成了矩
形的4个顶点和中心,其坐标分别如下,格式为:
0
975011009750150100001
0102501100102501
图9的上图是距离向压缩后的图像,从图中可以看到5
条回波信号目标回波信号存在明显的距离徙动,需要进行校正。
图9的下图是通过最近邻点插值法校正后的图像,可以看出图像基本被校正为直线。
距离向压缩,未校正距离徙动的图像-400-200方位向距离向104距离向压缩,校正距离徙动后的图像-400-200方位向距离向104
图9距离向压缩后最近邻点插值的结果
图10为进行方位向压缩后形成的图像,可以明显看出5
个点目标,并且5个点目标构成了矩形的四个顶点及其中心。
方位向0X:
9998
方位向压缩后的图像-400-300-200-100X:
9777丫:
-:
-+04RGB:
,:
9777Y:
:
-+04RGB:
,丫:
:
-+04RGB:
,:
+04丫:
-:
-+04RGB:
,:
+04丫:
:
-+04RGB:
,距离向104
图10通过最近邻点插值生成的点目标图像
使用最近邻点插值的距离多普勒算法仿真结果图11
上图为通过距离压缩后的图像,图11的下图为通过sine插值法校正后的图像。
距离向压缩,未校正距离徙动的图像-400-200方位向距离向104距离向压缩,校正距离徙动后的图像-400-200方位向距离向104
图11距离向压缩后sine插值的结果
图12为进行方位向压缩后形成的图像,可以明显看出5
个点目标,并且5个点目标构成了矩形的四个顶点及其中心。
方位向压缩后的图像-400-300-200-100X:
9777丫:
-:
-+04RGB:
,:
9775Y:
:
-+04RGB:
,方位向0X:
1e+04
Y:
:
-+04RGB:
,:
+04Y:
-:
-+04RGB:
,:
+04Y:
:
-+04RGB:
,距离向104
图12通过sine插值生成的点目标图像
ChirpSealing算法仿真结果同样,在ChirpSealing中,对5个点目标的回波信号进行了仿真,5个点目标构成
了矩形的4个顶点和中心,其坐标分别如下,格式为:
1200011250-501125050
11150-501
1150501图13是仿真的雷达回波信号图
仿真出来的信号50100150200距离向25030035040045050010020XX00400500600方位向7008009001000
图13仿真出来的SAR回波信号
图14是经过第一次相位校正之后,通过距离向压缩后的距离时域-方位时域信号图。
按照理论,该图中所有点的距离徙动都应该一样。
从图中大致可看出,五个点的距离徙动是差不多的。
ChirpScaling后、经过距离向压缩,距离徙动一致
50100150200距离向
25030035040045050010020XX00400500600方位向7008009001000
图。
图14ChirpScaling第2步之后、经过距离向压缩得到的图
图15为步骤5之后,信号距离压缩,距离徙动校正之后的距离多普勒域中的信号
消除距离徙动后的信号50100150200距离向
25030035040045050010020XX00400500600方位向7008009001000
图15距离徙动校正之后的图图16为步骤6之后,消除相位偏移的图。
相位校正后的信号50100150200距离向25030035040045050010020XX00400500600方位向7008009001000
图16消除相位偏移的图
图17为通过ChirpScaling算法生成的点目标图像
生成的点目标50100X:
384丫:
170lndex:
:
,:
642丫:
170lndex:
:
,距离向X:
513Y:
254lndex:
:
,:
386Y:
340lndex:
:
,:
641Y:
340lndex:
:
,方位
向6007008009001000
图17通过ChirpScaling算法生成的点目标图像
几种算法比较讨论了距离多普勒算法和Chirp
Scaling算法,其中距离多普勒算法考虑了最近邻点插值和sine插值两种插值方法。
距离多普勒算法兼顾了成熟、简
单、高效和精确等因素,至今仍被广泛使用,但是距离多普勒算法有两点不足:
首先,当用较长的核函数提高距离徙动校正精度时,运算量较大;其次,二次距离压缩对方位频率的依赖性问题较难解决,从而限制了其对某些大斜视角和长孔径SAR的处理精度。
最近邻点插值的优点是速度快,该
插值的运行时间为秒,缺点是不够精确;sine插值的优点是
精确,该方法的运行时间为秒,缺点是速度慢;ChirpScaling算法避免了插值运算,提高了速度,运行时间为秒,但是其
算法较为复杂。
%%================================================================%文件名:
%%徐一凡
%(功能:
合成孔径雷达距离多普勒算法点目标成像
%%=============================================
clear;clc;closeall;
%%
%常数定义
标区域参数
Xmax=50;
丫0=500;
最短距离%新线参数
D=4;
%方位向天线长
Lsar=lambda*R0/D;%SAR合成孔径长度,《合成孔径雷达成像一一算法与实现》Tsar=Lsar/V;%SARW射时间%%
慢时间域参数
Ka二2*VT/lambda/R0;%多普勒频域调频率
Ba=abs(Ka*Tsar);%多普勒频率调制带宽
PRF=Ba;9脉冲重复频率,PRF其实为多普勒频率的
采样率,又为复频率,所以等于
PRT=1/PRF;%永冲重复时间
ds=PRT;%慢时域的时间步长
Nslow=ceil((Xmax-Xmin+Lsar)/V/ds);%慢时域的采样数,ceil为取整函数,结合的图理解
Nslow=2Anextpow2(Nslow);%nextpow2为最靠近2
的幕次函数,这里为fft变换做准备
sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);
%慢时间域的时间矩阵
PRT=(Xmax-Xmin+Lsar)/V/Nslow;%于Nslow改变
了,所以相应的一些参数也需要更新,周期减小了PRF=1/PRT;ds=PRT;
%淋时间域参数设置
Tr=5e-6;%永冲持续时间5us
Br=30e6;%chirp频率调制带宽为30MHz
Kr=Br/Tr;%chirp调频率
Rmin=sqrt((Yc-丫0)八2+日八2);
Rmax=sqrt((Yc+Y0)A2+HA2+(Lsar/2)A2);
Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt);%快时域的
采样数量
Nfast=2Anextpow2(Nfast);%更新为2的幕次,方
便进行fft变换
tm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast);%快时域的离
散时间矩阵dt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast;%更
新间隔Fsr=1/dt;
%%b辨率参数设置
DY=C/2/Br;%距离向分辨率DX=D/2;%方位向
分辨率%%点目标参数设置
Ntarget=5;%点目标的数量%点目标格式[x,y,反
射系数sigma]
Ptarget=[Xmin,Yc-50*DY,1%点目标位置,这里设
置了5个点目标,构成一个矩形以及矩形的中心
Xmin+50*DX,Yc-50*DY,1
Xmin+25*DX,Yc,1Xmin,Yc+50*DY,1
Xmin+50*DX,Yc+50*DY,1];disp('Parameters:
')%参数显示
domain');disp(Fsr/Br)disp('SamplingNumberin
fast-timedomain');disp(Nfast)disp('SamplingRateinslow-timedomain');disp(PRF/Ba)disp('SamplingNumberinslow-timedomain');disp(Nslow)disp('Range
Resolution');disp(DY)
disp('Cross-rangeResolution');disp(DX)
disp('SARintegrationlength');disp(Lsar)
disp('Positionoftargets');disp(Ptarget)
%%================================================================%生成回波信号
K=Ntarget;%目标数目
N=Nslow;%慢时域的采样数M=Nfast;%快时域
的采样数T=Ptarget;%目标矩阵
Srnm=zeros(N,M);%生成零矩阵存储回波信号for
k=1:
1:
K%总共K个目标
sigma=T(k,3);%得到目标的反射系数
Dslow=sn*V-T(k,1);%方位向距离,投影到方位向
的距离R=sqrt(Dslow.A2+T(k,2)A2+HA2);%实际距
离矩阵
tau=2*R/C;%回波相对于发射波的延时
Dfast=ones(N,1)*tm-tau'*ones(1,M);%(t-tau),
其实就是时间矩阵,ones(N,1)和ones(1,M)都是为了将其扩
展为矩阵
phase=pi*Kr*Dfas
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- SAR 合成孔径雷达 图像 目标 仿真 报告 matlab 代码