发明名称 基于变化检测技术的SAR图像舰船目标鉴别方法
摘要 本发明提供一种基于变化检测技术的SAR图像舰船目标鉴别方法。技术方案是:首先利用ROI切片灰度数据构建变化检测量,生成标准灰度差异图像;然后,通过对差异图像进行阈值分割,获取潜在目标像素数量和真实目标像素数量;最后,根据真实目标像素数量与潜在目标像素数量之比,利用目标的像素聚集度对ROI切片进行鉴别。本发明能够有效提高SAR图像中舰船目标和海面杂波虚警的区分能力和鉴别方法的适应能力。
申请公布号 CN104240257A 申请公布日期 2014.12.24
申请号 CN201410519668.7 申请日期 2014.09.30
申请人 中国人民解放军国防科学技术大学 发明人 匡纲要;熊博莅;赵凌君;陆军;张小强;吴健
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 一种基于变化检测技术的SAR图像舰船目标鉴别方法,所述SAR是指合成孔径雷达,其特征在于,包括下述步骤:第一步,构建变化检测量:记SAR图像ROI切片为I,切片为正方形,尺寸为N×N,ROI切片中像素位置坐标为(i,j)的灰度表示为I(i,j),其中i,j=1,2,3...N,ROI是指感兴趣区域,利用下式计算ROI切片四角位置的四个正方形区域的灰度均值μ<sub>1</sub>,μ<sub>2</sub>,μ<sub>3</sub>,μ<sub>4</sub>,其中M=[N/4],[·]表示四舍五入取整操作:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&mu;</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>M</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>M</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mn>3</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>M</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>M</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000581975600000011.GIF" wi="803" he="710" /></maths>    (公式一)则背景杂波均值μ为:μ=(μ<sub>1</sub>+μ<sub>2</sub>+μ<sub>3</sub>+μ<sub>4</sub>)/4    (公式二)利用下式计算似然比变化检测量η(i,j):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&eta;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>&mu;</mi><mrow><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow><mi>&mu;</mi></mfrac></mrow>]]></math><img file="FDA0000581975600000012.GIF" wi="707" he="165" /></maths>    (公式三)其中η(i,j),i,j=1,2,3...N;利用线性变换方法,将似然比变化检测量η(i,j)的取值范围线性变换至[0,255],得到标准灰度差异图像D(i,j),i,j=1,2,…N;第二步,自适应阈值分割:对标准灰度差异图像D(i,j),令p<sub>k</sub>为D(i,j)=k的概率,即灰度级k的出现概率,k=0,1,…255;对任意阈值T∈[1,254],利用下式计算熵H(T):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mi>ln</mi><msub><mi>P</mi><mi>T</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>T</mi></msub><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>H</mi><mi>T</mi></msub><msub><mi>P</mi><mi>T</mi></msub></mfrac><mo>+</mo><mfrac><mrow><mi>H</mi><mo>-</mo><msub><mi>H</mi><mi>T</mi></msub></mrow><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>T</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000581975600000013.GIF" wi="866" he="164" /></maths>    (公式四)其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>T</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>T</mi></munderover><msub><mi>p</mi><mi>k</mi></msub><mo>,</mo><mi>H</mi><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mn>255</mn></munderover><msub><mi>p</mi><mi>k</mi></msub><mo>&CenterDot;</mo><mi>ln</mi><msub><mi>p</mi><mi>k</mi></msub><mo>,</mo><msub><mi>H</mi><mi>T</mi></msub><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>T</mi></munderover><msub><mi>p</mi><mi>k</mi></msub><mi>ln</mi><msub><mi>p</mi><mi>k</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000581975600000021.GIF" wi="1257" he="165" /></maths>令阈值<img file="FDA0000581975600000022.GIF" wi="503" he="120" />用阈值T<sub>0</sub>对标准灰度差异图像D(i,j)进行分割得到二值图像B:即灰度值大于阈值T<sub>0</sub>的像素赋值为1,否则为0;在二值图像B中设赋值为1的像素为潜在目标像素,统计潜在目标像素数量,记为N<sub>1</sub>;第三步,目标像素聚集度分析:对二值图像B进行区域生长,计算二值图像B中心区域连通像素的个数,记为N<sub>2</sub>,计算目标像素聚集度<img file="FDA0000581975600000023.GIF" wi="219" he="159" />根据实际场景设定判决阈值t,若ρ>t,则认为当前的ROI切片中存在舰船目标,否则不存在舰船目标。
地址 410073 湖南省长沙市开福区德雅路109号