发明名称 基于多光谱图像滤波输出统计分布规律的目标检测方法
摘要 本发明的目的在于提供基于多光谱图像滤波输出统计分布规律的目标检测方法,包括以下步骤:获取被测区域的多光谱图像,并进行辐射校正、分块预处理;采用CEM算法对获取的多个波段图像进行滤波处理;采用EM算法对滤波输出结果的分布参数进行估计;根据参数估计结果,判断目标的有无;确定检测门限,判决目标。本发明不需要目标和背景的空域信息,在目标信号比较微弱时也有较好的检测效果。适合于绿化面积检测、水灾检测等目标数目较多的场合的目标检测。同时,它直接估计出目标成分滤波处理后的统计分布参数,在估计目标的数目、确定检测门限时,具有更高的精度。本发明也可以用于高光谱和超光谱图像的目标检测。
申请公布号 CN102222221A 申请公布日期 2011.10.19
申请号 CN201110183192.0 申请日期 2011.07.01
申请人 哈尔滨工程大学 发明人 崔金辉;李立;张洪伟
分类号 G06K9/00(2006.01)I;G06K9/36(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 代理人
主权项 1.基于多光谱图像滤波输出统计分布规律的目标检测方法,其特征是:(1)获取被测区域的多光谱图像,并进行辐射校正、分块预处理;(2)采用CEM算法对获取的多个波段图像进行滤波处理:待检测的目标光谱信息d,波段数为L,滤波系数为L维列向量、由下式计算得到:<maths num="0001"><![CDATA[<math><mrow><msub><mi>w</mi><mi>cem</mi></msub><mo>=</mo><mfrac><mrow><msup><msub><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>d</mi></mrow><mrow><msup><mi>d</mi><mi>T</mi></msup><msubsup><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>d</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>像元r<sub>i</sub>对应的输出为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><msub><mi>&delta;</mi><mi>CEM</mi></msub><mrow><mo>(</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>w</mi><mi>cem</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><msub><mi>r</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msubsup><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>d</mi></mrow><mrow><msub><mi>d</mi><mi>T</mi></msub><msubsup><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>d</mi></mrow></mfrac><mo>)</mo></mrow><mi>T</mi></msup><msub><mi>r</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>d</mi><mi>T</mi></msup><msubsup><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>r</mi><mi>i</mi></msub></mrow><mrow><msup><mi>d</mi><mi>T</mi></msup><msubsup><mi>R</mi><mrow><mi>L</mi><mo>&times;</mo><mi>L</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>d</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中w<sub>cem</sub>为滤波系数,像元r<sub>i</sub>为L维列向量,<img file="FDA0000073161500000013.GIF" wi="276" he="119" />是L×L维的自相关矩阵,N是图像中的像元总数;(3)采用EM算法对滤波输出结果的分布参数进行估计:步骤(2)CEM滤波处理为线性运算,对图像中成分进行滤波处理,其滤波输出符合一维的高斯分布,整体图像滤波输出的概率密度函数为各个高斯概率密度函数的加权和,即服从高斯混合分布,表示为:<maths num="0003"><![CDATA[<math><mrow><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>a</mi><mi>i</mi></msub><msub><mi>N</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&Sigma;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中μ<sub>i</sub>和∑<sub>i</sub>为高斯分布的均值和方差,M为高斯模型个数,a<sub>i</sub>为某个高斯分布在混合模型中占有的比例,采用EM算法对高斯混合模型进行参数估计,选择高斯模型M,通过EM算法得出各个高斯分布的均值、方差和权值参数的估计结果;(4)根据参数估计结果,判断目标的有无:若对于第j个组分的高斯分布,下式为真时,判断目标存在:h1<a<sub>j</sub><h2 and c1>μ<sub>j</sub>>c2 and μ<sub>j</sub>>μ<sub>b</sub>,其中参数h1和h2是在图像中目标出现比例的上限和下限,c1和c2是目标的幅值范围,μ<sub>b</sub>是背景成分出现的最大均值;(5)确定检测门限,判决目标:通过步骤(4)判决图像内存在目标后,通过a<sub>j</sub>估计出图像内目标的数目,同时根据滤波输出的整体概率密度及第j个组分的高斯分布参数,选择虚警概率或漏警概率,确定检测门限,检测出目标。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室