目录
word文档和相关代码可以通过下面的连接下载:
概述
在介绍SAR信号性质与成像处理算法之前,先介绍SAR聚焦质量评价,有了这一套方案才能对后面介绍的成像处理算法进行评价。评价SAR成像处理软件聚焦的重要指标包括分辨率、峰值旁瓣比、积分旁瓣比。
下文已经结合sinc函数,介绍了一维信号的分辨率、峰值旁瓣比、积分旁瓣比的定义和测算方法。
本分重点介绍SAR图像聚焦质量评价插件的接口和使用方法。
插件使用说明
SARImageFocusedQuality1
功能
测试单视复图像聚焦指标,包括点目标峰值,方位向分辨率、峰值旁瓣比和积分旁瓣比,距离向分辨率、峰值旁瓣比和积分旁瓣。
调用方式与接口
调用方式:[Radiopara, Rgpara, Azpara] = SARImageFocusedQuality(sig_in, prf, sampleRate, Vg, sliceSize, interpNum)
输入参数:
sig_in:输入单视复图像,建议以定标器为中心,切出一个切片,提高程序运行效率
prf:脉冲重复频率,单位Hz
sampleRtae:采样率,单位Hz
Vg:波束地速,单位m/s
slizeSize:测指标时,以最大值为中心(定标器),裁切一个尺寸为slizeSize*slizeSize的切片
interpNum:插值倍数,测指标时要对信号进行插值
输出参数:
Radiopara:点目标辐射相关指标
Radiopara.peakVal:峰值点幅度
Radiopara.posa:峰值点方位向像素坐标
Radiopara.posr:峰值点距离向像素坐标
Rgpara:距离向聚焦指标
Rgpara.res:距离向分辨率,单位m
Rgpara.PSLR:距离向峰值旁瓣比,单位dB
Rgpara.ISLR:距离向积分旁瓣比,单位dB
Azpara:方位向聚焦指标
Azpara.res:方位向分辨率,单位m
Azpara.PSLR:方位向峰值旁瓣比,单位dB
Azpara.ISLR:方位向积分旁瓣比,单位dB
工作原理与流程
函数工作的基本流程:
- 找出场景中亮度最大的点,并认为这一个点就是定标器;
- 输出定标器的像素坐标[Radiopara.posa, Radiopara.posr];
- 以定标器为中心,裁切一个尺寸为sliceSize* sliceSize的切片;
- 对切片进行interpNum倍二维插值;
- 插值后再搜索亮度最大值;
- 输出定标的幅度峰值Radiopara.peakVal
- 在插值后的峰值点处,裁切一条距离向信号,测试并输出距离向指标:
- Rgpara.res:距离向分辨率,单位m
- Rgpara.PSLR:距离向峰值旁瓣比,单位dB
- Rgpara.ISLR:距离向积分旁瓣比,单位dB
- 在插值后的峰值点处,裁切一条方位向信号,测试并输出方位向指标
- Azpara.res:方位向分辨率,单位m
- Azpara.PSLR:方位向峰值旁瓣比,单位dB
- Azpara.ISLR:方位向积分旁瓣比,单位dB:
在调用程序时,输入的单视复图像建议是人工截取的切片,且截取后图像中只有一个强目标,否则本插件只会测最强点的指标,忽略其他点。
调用样例
%SAR图像聚焦质量评价
% 使用事先生成的点目标仿真数据
% Young 2022-05-06
close all
clear all
%%
load testData.mat;
%%
sliceSize = 32;
interpNum = 16;
pointSize = 512;
[Radiopara, Rgpara, Azpara] = ...
SARImageFocusedQuality1(sig_image, prf, sampleRate, Vst, sliceSize, interpNum);
fprintf('切片峰值: %f, 坐标:[%d, %d]\n', Radiopara.peakVal, Radiopara.posa, Radiopara.posr);
fprintf('------- 距离向指标 -------\n');
fprintf('峰值: %8.2f ',Rgpara.max);
fprintf(' 分辨率: %4.3f', Rgpara.res);
fprintf(' 峰值旁瓣比:%5.3f', Rgpara.PSLR);
fprintf(' 积分旁瓣比:%5.3f\n', Rgpara.ISLR);
fprintf('------- 方位向指标 -------\n');
fprintf('峰值: %8.2f ',Azpara.max);
fprintf(' 分辨率: %4.3f', Azpara.res);
fprintf(' 峰值旁瓣比:%5.3f', Azpara.PSLR);
fprintf(' 积分旁瓣比:%5.3f\n', Azpara.ISLR);
运行上述程序,插件输出插值后的点目标,如下图所示。
插件输出插值后的一维信号波形和频谱,如下图所示
程序输出的指标如下图所示。
本博客上传的所有原创类资料,仅可用于个人学习、交流等非商业性用途。未经作者同意,不可用于任何商业用途。
如在文献、报告等文字材料中引用本博客发布的原创类资料,请注明出处。