一种sar图像精确度计算方法
【专利摘要】本发明公开了一种SAR图像精确度计算方法，依据SAR图像中目标的辐射和结构特征构造精确度计算方法，通过计算SAR目标电磁仿真图像的形状精确度、散射统计信息精确度和结构精确度，综合判定SAR目标电磁仿真图像的精确度高低，本发明有效解决了当前无法客观评价SAR目标电磁仿真图像精确度高低的问题。
【专利说明】—种SAR图像精确度计算方法
【技术领域】
[0001]本发明属于SAR信号处理领域，具体涉及一种图像处理方法。
技术背景
[0002]合成孔径雷达SAR(Synthetic Aperture Radar,简称:SAR)因其全天时、全天候和高分辨率的特点，广泛地应用于民用和军用领域。SAR图像反映的是被测区域对电磁波的散射特性分布，目标具有的辐射和结构等目标特征能够使其在雷达自动目标识别中具有极大的应用价值。构建目标在各种不同成像条件下的SAR图像数据库，对于自动目标识别理论和算法性能研究具有重要作用。
[0003]文献“Ross T D, Bradley J J, Hudson L J, et al.SAR ATR: so what’s theproblem ? An MSTAR perspective[C].AeroSense' 99.1nternat1nal Society for Opticsand Photonics, 1999:662-672.”中指出实验获取地面目标SAR图像面临着成本过高、数据量过大、操作复杂等困难。因此，寻求一种新的方便易行、成本较低、可以快速获取SAR图像的方法具有十分重要的意义。
[0004]文献“MishraA M, Bernard M.Database generat1n of bistaticground target signatures[C].Wireless Communicat1ns and AppliedComputat1nal Electromagnetics, 2005.1EEE/ACES Internat1nal Conferenceon.1EEE, 2005:523-528.”中指出，通过电磁计算软件仿真，可以快速得到大量SAR目标电磁仿真图像。该方法方便易行、费用低廉，不仅可以避免真实实验所面临的成本大、操作复杂等问题，也解决了电磁SAR目标电磁仿真图像面临的SAR图像采集困难、时间损耗大的问题。文中直接将SAR目标电磁仿真图像用于自动目标识别，指出SAR目标电磁仿真图像在自动目标识别领域的特征提取和识别算法性能研究中具有重大作用，但并未提出用于评价SAR目标电磁仿真图精确度的方法。因此，寻找一种合理的、客观的评价SAR目标电磁仿真图像精确度高低的方法，具有十分重要的意义。
【发明内容】
[0005]本发明为了解决当前无法客观评价SAR目标电磁仿真图像精确度高低的问题，提出了一种SAR图像精确度计算方法。该方法通过构造形状精确度、散射信息精确度等多种精确度准则，判定通过电磁仿真软件获得的SAR目标电磁仿真图像的精确度高低。
[0006]本发明提供了一种SAR图像精确度计算方法，具体包括以下步骤:
[0007]步骤SI:SAR目标电磁仿真,步骤SI具体包括以下步骤:
[0008]步骤Sll:目标建模；
[0009]步骤S12:设置电磁计算仿真参数，目标模型仿真；
[0010]步骤S13:仿真数据成像；
[0011]步骤S2:SAR目标电磁仿真图像精确度计算,步骤S2具体包含以下分步骤:
[0012]步骤S21:计算形状精确度；
[0013]步骤S22:计算统计信息精确度；
[0014]步骤S23:计算结构精确度；
[0015]步骤S24:计算SAR图像平均精确度；
[0016]进一步地，所述目标建模根据目标三维尺寸，运用三维建模软件实体建模方法建造目标各个部件三维模型，将建立的目标各部件三维模型导入电磁计算软件中，利用电磁计算软件的几何位置调节功能将各部件三维模型组合为整体；
[0017]进一步地，所述形状精确度包括:面积精确度和周长精确度；
[0018]进一步地，所述统计信息精确度包括:图像均值精确度、图像标准差精确度和信息熵精确度；
[0019]更进一步地,所述SAR图像平均精确度Sim计算公式为:
[0020]Sini = (AS + CS+HS+β?Β + c& SSIM)6 ;
[0021]其中，AS表示面积精确度，CS表示周长精确度，HS表示图像信息熵精确度，μ S表示图像均值精确度，σ S表示图像标准差精确度以及表示平均结构精确度。
[0022]本发明公开的一种SAR图像精确度计算方法，通过计算SAR目标电磁仿真图像的形状精确度、统计信息精确度、结构精确度以及SAR图像平均精确度，综合判定通过电磁仿真软件获得的SAR目标电磁仿真图像的精确度高低，该方法有效解决了当前无法客观评价SAR目标电磁仿真图像精确度高低的问题。
[0023]附图表说明
[0024]图1本发明方法的流程图。
[0025]图2本发明实施步骤中Τ72Μ1坦克光学图与三维模型三视图。
[0026]图3本发明实施步骤预处理后的SAR实测图像与SAR目标电磁仿真图像对比图。
具体实施例
[0027]为了方便描述以及理解本发明实施例的内容，首先对以下术语进行解释:
[0028]电大尺寸:当物理尺寸远远大于一个波长的时候，认为该结构是电大尺寸。
[0029]小角度转台成像:转台成像，即假定目标放在一个旋转的平台上，雷达位置保持不动，目标相对于雷达只有旋转运动。当目标相对于雷达旋转角度很小时的成像模式，即称为小角度转台成像。
[0030]本实施例建模目标为MSTAR数据库中的Τ72Μ1主战坦克，使用CATIA软件对Τ72Μ1进行三维建模，将三维模型导入FEKO软件进行电磁计算获得目标RCS数据，再将RCS数据进行SAR成像处理。SAR目标电磁仿真图像与MSTAR SAR实测图像的相似性的验证所有步骤都在MATLAB2010上进行，具体过程如下:
[0031]步骤SI:Τ72Μ1坦克三维建模、模型仿真与仿真数据成像；具体包括以下分步骤:
[0032]步骤Sll:Τ72Μ1坦克三维建模；
[0033]参照Τ72Μ1坦克缩版模型三维尺寸，运用CATIA实体建模方法，通过草图编辑、平面拉伸、旋转、凹槽和倒角等功能，建造Τ72Μ1坦克炮台、车体及车轮三部件三维模型，将CATIA软件导出的Τ72Μ1坦克炮台、车体及车轮三部件三维模型依次导入FEKO软件中，运用FEKO软件中的平移、旋转等位置调节功能，按照真实Τ72Μ1坦克几何位置关系组合炮台、车体及车轮三部件。
[0034]步骤S12 =FEKO软件仿真参数设置，步骤S12具体包括以下分步骤:
[0035]步骤S121:设置频率下限&、频率上限fH、频率(距离向)采样点数Nx、方位向采样点数Ny、带宽B和方位向观测角范围Ω ；
[0036]步骤S122:选择电磁计算方法及与置网格剖分尺寸；
[0037]所述电磁计算方法选择物理光学；
[0038]步骤S123:俯仰角Θ、方位角Ψ与极化方式设置；
[0039]所述俯仰角Θ范围为[0°，90° ]；
[0040]所述方位角f范围为[0°，360° ]；
[0041]所述远场均匀平面电磁波的线性极化方式为，水平极化H和垂直极化V ；
[0042]仿真俯仰角73° (MSTAR数据库目标俯仰角余角)，分辨率为Λ χ = 0.3m、Ay =
0.3m、X波段且中心频率为10GHz，HH极化方式的SAR图像；
[0043]确定FEKO软件仿真参数，即仿真带宽B及频率范围(fu fH)，方位向观测角范围Ω及图像采样点数NxXNy ；
[0044]所述频率范围(fu fH)及频率(距离向)采样点Nx设置:
k C
[0045]由Δ χ = 0.3m，c = 3*108m/s，通过公式可求得 B = 500MHz ；
[0046]由B = 500MHz, fc = 10GHz，通过公式 fc = (fL+fH)/2 和公式 B = (fH-fL)，可求得fL = 9.75GHz, fH = 10.25GHz ；
[0047]由Xniax = 6.95m, Δ χ = 0.3m,通过公式求得最少采样点Nx = 24，本发明取距离向采样点51个；
[0048]因此，设置频率范围(fL, fH)为(9.75GHz, 10.25GHz)，距离向采样点Nx为51个采样点；
[0049]所述方位向观测角范围Ω、方位角f?,方位向采样点Ny设置:
[0050]通过公式Af = I可求得Ω = 2.8652° ；
Y
[0051]由Ymax = 4.5m, Δ y = 0.3m,通过公式ftfF求得最少采样点为Ny = 15，取方位向米样点Ny为51可完全表现方位向为4.5m的目标；
[0052]方位角炉将方位向观测角范围Ω划分为-Ω/2 +炉和+Ω/2 +炉两个区域。炉取值可在[0°，360° ]中任意取值，本发明的f为2°、12°、22°等与MSTAR数据库中真实目标相同方位角的22组值。
[0053]所述网格尺寸及电磁计算方法设置；
[0054]设置电磁计算方法为LEPO法，通过网格尺寸公式l = y可得I = 0.058496m，又
Jh，
Ic
因FEKO软件提供LEPO电磁计算方法网格尺寸范围为可大于~7"设置网格尺寸为0.lm。
J? f
[0055]步骤S13: SAR成像仿真；
[0056]将FEKO仿真的RCS数据导入MATLAB，利用小角度二维FFT算法对仿真数据进行二维傅里叶变换、加窗等成像处理，最终得到目标的SAR目标电磁仿真图像；
[0057]步骤S2 =SAR目标电磁仿真图像精确度计算；
[0058]对方位角f为2°、12°、22°等22组仿真SAR图像和实测SAR图像进行图像分害I]、质心配准和能量归一化预处理；
[0059]采用精确度准则，计算能量归一化后的SAR目标电磁仿真图像形状精确度，散射统计信息精确度和结构精确度；步骤S2具体包括以下分步骤:
[0060]步骤S21:计算形状精确度；
[0061 ] 分别计算能量归一化后SAR目标电磁仿真图像目标区域面积Aa和SAR实测图像目标区域面积Ab，即质心配准后图像中幅度值不为O的像素值总和，通过公式
=求得方位角f为2°、12°、22°等22组每组SAR目标电磁仿真图像面积
-4M
精确度AS ；
[0062]通过边缘检测算子提取能量归一化后SAR目标电磁仿真图像和SAR实测图像边缘，计算SAR目标电磁仿真图像目标边缘长度Ca与真实SAR图像目标边缘长度CB，通过公
j C —C I
式C5 = (l--一^-)χ100%，求得方位角φ为2。>12° >22°等22组每组SAR目标电磁仿真图像周长精确度CS ；
[0063]步骤S22:计算统计信息精确度；
M N^ Α:，I
[0064]通过公式 ΣΣ/Μ公式^分别计算能量归
a =..........—-■■■■—---σ = \——;---
MkNI M XiV
一化后SAR目标电磁仿真图像和SAR实测图像的均值μ、标准差σ，再通过公式?χ 100%、公式於=χ丨00%求得方位角供为2。、12。、22。等22
组每一组SAR目标电磁仿真图像图像均值精确度μ S和图像标准差精确度σ S ；
[0065]对能量归一化后的SAR目标电磁仿真图像和SAR实测图像进行灰度值量化，将灰度值区间转换为[0，255]，通过公式分别计算灰度值量化后的SAR目标电磁仿真图像与SAR实测图像平均信息熵H，通过公式《^ = ^1^x100%求出方位角-
为2°、12°、22°等22组每组SAR目标电磁仿真图像的平均图像信息熵精确度HS ；
[0066]步骤S23:计算结构精确度；
[0067]采用3X3大小的窗口分别从SAR实测图像A和SAR目标电磁仿真图像B的左上角像素点，从左到右、从上往下逐个像素点滑动，通过公式^a3) =爿公式
"Λ +/? +(1、
C7* C7* -?- (J^y1 j
e(A,B) =......................f...............g......;.......-..................公式?(Α,Β) = ^L 和公式 SS頂=I (A, B) c (A, B) s (A, B)，依次
σΛ +σΒ +C2、CFaCJb +€3
计算相同位置处窗口所对应的子块图像的亮度比较I (A，8)、对比度比较。(、B)、结构比较s (A, B)和结构精确度SSIM，其中，C1 = 0.001，C2 = 0.0001，C3 = 0.0002，将滑窗窗口总数
j r
T = 3969带入公Sssim分别求出方位角妒为2°、12°、22°等22组每组
SAR目标电磁仿真图像的平均结构精确度SSIM.
[0068]SSIM 'VSSIM1
1 H ?
[0069]其中，SSMk为第k个窗口的SSM值；
[0070]步骤S24 =SAR图像平均精确度；
[0071]通过公式=+ + + + 痛!M)/6 计算方位角 P 为 2。、12。、
22°等22组每组SAR目标电磁仿真图像形状精确度、散射统计信息精确度和结构精确度的平均精确度Sim。所述方位角f与平均精确度Sim的对应关系见下表:
[0072]
方位(度) 2 12 22 32 42 50 62 72 82 92 100
0,883 0,873 0,867 0.869 0.853 I 0Λ8 0.836 0.831 0.808 0,847 0,825
f均糖度
3 8 7 O I I 2 I 5 3 9 6 8方位(It) 112 122 132 】42 152 162 172 182 192 202 211
0,858 0,822 0.799 0.875 I 0.899 I 0.862 0.896 0.86S 0.844 0,859 0,879
平均鋳确》f
OI 4 I 4 [ 6 j 3 j 40 9 8 4 3
[0073]实验结果表明，本发明所使用的方法可以客观有效地评价SAR目标电磁仿真图像精确度高低。
[0074]综上可以看出，本发明的一种SAR图像精确度计算方法，通过计算SAR目标电磁仿真图像的形状精确度、统计信息精确度、结构精确度以及SAR图像平均精确度，综合判定通过电磁仿真软件获得的SAR目标电磁仿真图像的精确度高低，该方法有效解决了当前无法客观评价SAR目标电磁仿真图像精确度高低的问题。
[0075]本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。
【权利要求】
1.一种SAR图像精确度计算方法，其特征在于，包括以下步骤: 步骤SI =SAR目标电磁仿真，步骤SI具体包括以下分步骤: 步骤Sll:目标建模； 所述目标建模根据目标三维尺寸，运用三维建模软件实体建模方法建造目标各个部件三维模型，将建立的目标各部件三维模型导入电磁计算软件中，利用电磁计算软件的几何位置调节功能将各部件三维模型组合为整体； 步骤S12:设置电磁计算仿真参数，目标模型仿真； 步骤S13:仿真数据成像； 步骤S2 =SAR目标电磁仿真图像精确度计算，步骤S2具体包含以下分步骤: 步骤S21:计算形状精确度； 所述形状精确度包括:面积精确度和周长精确度； 步骤S22:计算统计信息精确度； 所述统计信息精确度包括:图像均值精确度、图像标准差精确度和信息熵精确度； 步骤S23:计算结构精确度； 步骤S24:计算SAR图像平均精确度； 所述SAR图像平均精确度Sim计算公式为:
其中，AS表示面积精确度，CS表示周长精确度，HS表示图像信息熵精确度，μ S表示图像均值精确度，σ S表示图像标准差精确度以及远而表示平均结构精确度。
2.根据权利要求1的一种SAR图像精确度计算方法，其特征在于，所述步骤S12包括以下分步骤: 步骤S121:设置频率下限4、频率上限fH、频率采样点数Nx、方位向采样点数Ny、带宽B和方位向观测角范围Ω ； 步骤S122:选择电磁计算方法及与设置网格剖分尺寸； 所述电磁计算方法选择物理光学； 步骤S123:俯仰角Θ、方位角-与极化方式设置； 所述远场均匀平面电磁波的线性极化方式为，水平极化H和垂直极化V。
3.根据权利要求2的一种SAR图像精确度计算方法，其特征在于，所述俯仰角Θ范围为[0°，90° ]，所述方位角-范围为[0°，360° ]。
4.根据权利要求1的一种SAR图像精确度计算方法，其特征在于，步骤S21所述形状精确度包括:面积精确度AS和周长精确度CS ； 所述面积精确度AS计算公式为:
其中，Aab为预处理后SAR目标电磁仿真图像与SAR实测图像目标相同区域面积，Ab为预处理后SAR实测图像目标区域面积的；
所述周长精确度CS计算公式为
其中，Ca为预处理后SAR目标电磁仿真图像目标边缘长度，Cb为预处理后SAR实测图像目标边缘长度。
5.根据权利要求1的一种SAR图像精确度计算方法，其特征在于，步骤S22所述统计信息精确度包括:图像均值精确度μ S、图像标准差精确度σ S和信息熵精确度HS，分别计算SAR目标电磁仿真图像和真实SAR图像的均值μ、标准差。和平均信息熵H ； 所述SAR图像的均值μ计算公式为:
所述标准差σ计算公式为:
所述平均信息熵H计算公式为
* 其中，M、N为图像矩阵大小，L为灰度值最大值，f(i，j)为图像中(i，j)坐标位置处的灰度值，Pk为灰度值k出现的概率； 所述图像均值精确度μ S计算公式为:
其中，μ A为预处理后SAR目标电磁仿真图像均值，μ B为预处理后SAR实测图像均值； 所述图像标准差精确度σ S计算公式为
； 其中，σ A为预处理后SAR目标电磁仿真图像标准差，O B为预处理后SAR实测图像标准差； 所述图像信息熵精确度Hs计算公式为:
其中，Ha为预处理后SAR目标电磁仿真图像平均信息熵，Hb为预处理后SAR实测图像平均息熵。
6.根据权利要求1的一种SAR图像精确度计算方法，其特征在于，所述步骤S23具体包括以下分步骤: 步骤S231:采用相同大小的窗口分别从SAR实测图像和SAR目标电磁仿真图像的左上角像素点，从左到右、从上往下逐个像素点滑动； 步骤S232:依次计算相同位置处窗口所对应的子块图像的亮度比较、对比度比较、结构比较和结构精确度SSIM ; 所述结构精确度SSM计算公式为:SSM = I (A, B) c (A, B) s (A, B)； 所述亮度比较1(A，B)计算公式为:
所述对比度比较C (A，B)计算公式为
所述结构比较S(A，B)计算公式为:
其中，A为SAR目标电磁仿真图像，B为SAR实测图像，μ Α、μΒ分别为A、B图像中相同滑窗位置的子块图像均值，σ Α、σ B分别为A、B图像中相同滑窗位置的子块图像标准差，σ AB为A、B图像中相同滑窗位置的子块图像的协方差，C1, C2, C3为小于I的正数； 所述μ a公式为
所述μ B公式为
所述σ a计算公式为:
所述。B和计算公式为:σ
所述 σ 计算公式为
步骤S233:求出SAR实测图像和SAR目标电磁仿真图像的整幅图像平均结构精确度SSIMi 所述平均结构精确度?而计算公式为
其中，T为滑窗总数，SSIMk为第k次滑窗的SSM值。
7.根据权利要求1的一种SAR图像精确度计算方法，其特征在于，所述仿真数据的成像处理，采用小角度二维FFT算法，对电磁仿真数据进行二维傅里叶变换、加窗、补零操作。
【文档编号】G06T17/00GK104167019SQ201410384240
【公开日】2014年11月26日 申请日期:2014年8月6日 优先权日:2014年8月6日
【发明者】黄钰林, 刘晓佳, 韩昊, 裴季方, 霍伟博, 武俊杰, 杨建宇 申请人:电子科技大学