SAR点目标分析

点目标分析

单一孤立散射体（如角反射器）称为点目标。通过测量点目标响应估计SAR图像质量。点目标测量给出的重要质量参数包括（1）冲激响应宽度（IRW）（2）峰值旁瓣比（PSLR）（3）积分旁瓣比（ISLR）

频域补零插值

点目标峰值在处理后的数据中仅占一两个像素，所以必须对峰值邻域内的目标响应进行插值，以便进行参数测量：

1. 一维 ，频谱中心是高频，fftshift函数是将零频点移至中间。插值是通过增加信号的高频分量使线条更加平滑。ifft(fft(sig)(1,m/2),zeros(1,n),fft(sig)(m/2+1,end))
2. 二维 ，频谱中心是高频，四个角是低频，补零将原来的频谱拆分四块放在四个角，其余都是零。对于点目标，低频能量更高。imagesc(data)会发现点目标黄色，没有目标处蓝色

Data=fft2(data);
Dataq(1:round(m/2),1:round(n/2))=Data(1:round(m/2),1:round(n/2));
Dataq(m*freq-m+round(m/2)+1:m*freq,1:round(n/2))=Data(round(m/2)+1:m,1:round(n/2));
Dataq(1:round(m/2),n*freq-n+round(n/2)+1:n*freq)=Data(1:round(m/2),round(n/2)+1:n);
Dataq(m*freq-m+round(m/2)+1:m*freq,n*freq-n+round(n/2)+1:n*freq)=Data(round(m/2)+1:m,round(n/2)+1:n);
dataq=ifft2(Dataq);

IRW

扩展后的数据可以经过峰值点截取水平和垂直方向上的幅相剖面。IRW是冲激响应宽度，指冲激响应的3dB宽度，在SAR处理中又叫做图像分辨率。
IRW的单位是像素，也可以是米。
fft(sig)后需要确定主瓣，然后确定3dB带宽。

寻找主瓣代码思路

dif = diff(sig(1:xm)) < 0;
ind = find(dif == 1,1,'last'); %找到最后一个单调递减
x1 = ind+1;
dif = diff(sig(xm+1:end)) > 0;
ind = find(dif == 1,1,'first');%找到第一个单调递增
x2 = ind + xm;

寻找3dB

20log₁₀sig+3转为dB量纲，后在峰值点前后找到3dB对应的点。

PSLR

最大旁瓣与主瓣高度比，dB表示。SAR系统中的PSLR必须小于-13dB（sinc函数对应的PSLR）
，以使弱目标不会被邻近的强目标掩盖。一般在-20dB左右。

[peak1,peak2] = findpeaks(sig,'SortStr','descend');

ISLR

ISLR为积分旁瓣比。
$$ISLR=10{\log }_{10}\frac{totalpwr-mainpwr}{mainpwr}$$
totalpwr和mainpwr分别代表功率总和和主瓣功率。