发明名称 一种基于均值漂移算法和数学形态学的深度图分割方法
摘要 一种基于均值漂移算法和数学形态学的深度图分割方法,涉及深度图分割领域,解决现有分割方法无法自动化程度高的将跨越了整幅图像灰度范围且不平坦的背景分割成一个整体的问题。根据深度图像生成深度图的灰度直方图;利用一维高斯窗函数对灰度直方图进行平滑处理;根据前景阈值T调节步骤一所述的深度图的灰度值,将灰度值高于所述前景阈值T的像素灰度提升10;将获得的深度图利用数学形态学“关”操作进行修正;利用均值漂移算法将获得的深度图进行聚类分割;对分割区域进行灰度值赋值,每个区域所赋灰度值为原深度图中该区域内所有灰度值的平均值;完成对深度图的分割。本发明可广泛应用于对灰度图的前景与背景分割的工作中。
申请公布号 CN103136753B 申请公布日期 2016.02.17
申请号 CN201310056959.2 申请日期 2013.02.25
申请人 哈尔滨工业大学 发明人 关宇东;提纯利;滕艺丹;戴翊轩;李尔佳;杜克;仲小挺;于博良
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于均值漂移算法和数学形态学的深度图分割方法,其特征在于它包括如下步骤:步骤一:根据灰度为f(p,q)的深度图像生成深度图的灰度直方图hist(k),k=0,1…L‑1;L为灰度级,L=256;所述的灰度直方图的横坐标为灰度级,纵坐标为灰度级出现的频率;步骤二:利用一维高斯窗函数对灰度直方图进行平滑处理,并设定前景阈值T;所述前景阈值T为按灰度值从高到低的第一个包含了一定语义范围的极小值,所述一定语义范围是指所含像素数大于整个图像像素数的10%的区域,即hist(k&gt;T)&gt;10%*hist(0…L‑1);所述极小值为hist(k)<sub>min</sub>=min(hist(k‑4),hist(k‑3)...hist(k+4));步骤三:根据前景阈值T调节步骤一所述的深度图的灰度值,将灰度值高于所述前景阈值T的像素的灰度提升10;步骤四:将步骤三获得的深度图利用数学形态学“关”操作进行修正;步骤五:利用均值漂移算法将步骤四获得的深度图进行聚类分割;所述均值漂移算法的公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&equiv;</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>G</mi><mi>H</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>G</mi><mi>H</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000779954670000011.GIF" wi="792" he="292" /></maths>     公式一其中,G<sub>H</sub>(x<sub>i</sub>‑x)=|H|<sup>‑1/2</sup>G(H<sup>‑1/2</sup>(x<sub>i</sub>‑x)),G(x)是一个单位核函数,带宽矩阵H是一个正定的对称d×d矩阵,记为一个对角矩阵<img file="FDA0000779954670000012.GIF" wi="447" he="94" />取正比于单位矩阵H=h<sup>2</sup>I,w(x<sub>i</sub>)≥0是一个赋给采样点x<sub>i</sub>的权重,选用二维的高斯模窗函数;因此,公式一改写为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>M</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>G</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>G</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mi>x</mi></mrow>]]></math><img file="FDA0000779954670000013.GIF" wi="715" he="297" /></maths>    公式二将公式二等式右侧第一项记为m<sub>h</sub>(x),即<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>m</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>G</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>G</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>x</mi></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000779954670000014.GIF" wi="624" he="307" /></maths>      公式三步骤六:给定初始像素x,容许误差ε;计算m<sub>h</sub>(x)并把m<sub>h</sub>(x)赋值给x,当||m<sub>h</sub>(x)‑x||<ε时结束赋值,否则持续将计算得到的m<sub>h</sub>(x)赋值给x;步骤七:步骤六最终获得的x为初始点的收敛点,将相同收敛与同一点的像素归为一类,即划为同一分割区域;步骤八:对分割区域进行灰度值赋值,每个区域所赋灰度值为原深度图中该区域内所有灰度值的平均值;完成对深度图的分割。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号