发明名称 基于K-means聚类和偏振信息提取的伪装目标识别方法
摘要 本发明公开了一种基于K-means聚类和偏振信息提取的伪装目标识别方法,首先获得目标偏振参量图像数据,然后获得偏振度图像、偏振角图像和椭率角图像,对于得到的多维图像I(P,θ,ε)用K-means进行聚类,计算新的聚类中心C<sub>mi</sub>,重复计算新的聚类中心,直至聚类中心C<sub>mi</sub>不再变化,迭代终止,得到最终的聚类融合图像。本发明在主观视觉和客观评价上有了很大的提高,聚类融合的结果由于融入了偏振度,偏振角和椭率角信息,使得图像中的偏振信息更完整。
申请公布号 CN105787499A 申请公布日期 2016.07.20
申请号 CN201410830874.X 申请日期 2014.12.26
申请人 南京理工大学 发明人 顾国华;赵蓉;杨蔚;周骁骏;张骏
分类号 G06K9/62(2006.01)I;G06T7/40(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 唐代盛
主权项 一种基于K‑means聚类和偏振信息提取的伪装目标识别方法,其特征在于步骤如下:步骤1,获得目标偏振参量图像数据:利用偏振成像系统对目标进行偏振成像,得到偏振成像系统中偏振片的透光轴与水平方向的夹角分别为0°、45°、90°、135°的线偏振分量光强图像I<sub>0</sub>、I<sub>45</sub>、I<sub>90</sub>和I<sub>135</sub>,加上1/4波片后,分别测量左旋和右旋圆偏振分量光强图像I<sub>L</sub>和I<sub>R</sub>,进而得到斯托克斯参量图像数据I、Q、U、V;步骤2,获得偏振度图像、偏振角图像和椭率角图像:由步骤1中的斯托克斯参量图像数据I、Q、U、V计算得到偏振度图像I(P)、偏振角图像I(θ)和椭率角图像I(ε);步骤3,对于步骤2中得到的多维图像I(P,θ,ε)用K‑means进行聚类,随机选取初始聚类中心C<sub>1i</sub>(P<sub>0</sub>,θ<sub>0</sub>,ε<sub>0</sub>),i=1,2...k,其中k为类别数,下标1表示第一次迭代,P<sub>0</sub>,θ<sub>0</sub>,ε<sub>0</sub>分别表示选取的偏振度图像、偏振角图像和椭率图像中某一位置像素的像素值;求解每个样本数据与初始聚类中心C<sub>1i</sub>的欧几里得距离<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>P</mi><mi>j</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>j</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>&epsiv;</mi><mi>j</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0000645247920000011.GIF" wi="856" he="79" /></maths>i=1,2...k,j=1,2...N,N为每幅图像中的总像素数,其中C<sub>1i</sub>(1)、C<sub>1i</sub>(2)和C<sub>1i</sub>(3)分别表示初始聚类中心C<sub>1i</sub>的第一、二、三个元素,P<sub>j</sub>、θ<sub>j</sub>和ε<sub>j</sub>分别表示偏振度图像、偏振角图像和椭率图像中像素按每行叠加后第j个像素的像素值;将每个数据对象分配到离该数据对象最近的初始聚类中心C<sub>1i</sub>(P<sub>0</sub>,θ<sub>0</sub>,ε<sub>0</sub>),i=1,2...k;步骤4,计算新的聚类中心<img file="FDA0000645247920000012.GIF" wi="755" he="120" />其中m为迭代次数,r<sub>i</sub>表示离C<sub>(m‑1)i</sub>,i=1,2...k最近的像素点个数,C<sub>(m‑1)i</sub>,i=1,2...k为第(m‑1)次迭代的第i类中心;并将每个数据对象分配到离它最近的聚类中心C<sub>mi</sub>,i=1,2...k;计算各个聚类中心之间的欧几里得距离d(C<sub>i</sub>,C<sub>j</sub>),i,j=1,2...k,i≠j,C<sub>i</sub>,C<sub>j</sub>分别表示第i类和第j类聚类中心;步骤5:重复步骤4,直至聚类中心C<sub>mi</sub>不再变化,迭代终止,得到最终的聚类融合图像。
地址 210094 江苏省南京市孝陵卫200号