发明名称 点云骨架提取方法及装置
摘要 一种点云骨架提取方法,包括:获取输入的点云采样数据;根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:<img file="DDA0000751479860000011.GIF" wi="1247" he="167" /></maths>其中,<img file="DDA0000751479860000012.GIF" wi="1036" he="182" /></maths>其中,i,j为正整数,J为点云采样数据的点集合,q<sub>j</sub>为点集合J中的采样点,I为采样点q<sub>j</sub>的邻域点集合,x<sub>i</sub>为邻域点集合I中的邻域点,<img file="DDA0000751479860000013.GIF" wi="45" he="56" />为邻域点集合I中不属于x<sub>i</sub>的邻域点,R为规整项,γ<sub>i</sub>为权重系数,h为邻域点集合I的邻域半径,σ<sub>i</sub>为分布系数;连接所述骨架分支得到点云骨架。此外,还提供了一种点云骨架提取装置。上述点云骨架提取方法及装置能够提高点云骨架的准确性。
申请公布号 CN103268631B 申请公布日期 2015.09.30
申请号 CN201310196243.2 申请日期 2013.05.23
申请人 中国科学院深圳先进技术研究院 发明人 黄惠;伍世浩;陈宝权;南亮亮
分类号 G06T17/00(2006.01)I 主分类号 G06T17/00(2006.01)I
代理机构 广州华进联合专利商标代理有限公司 44224 代理人 吴平
主权项 一种点云骨架提取方法,包括:获取输入的点云采样数据;根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>X</mi></munder><munder><mo>&Sigma;</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>&Element;</mo><mi>I</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><msub><mi>q</mi><mi>j</mi></msub><mo>&Element;</mo><mi>J</mi></mrow></munder><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>q</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo><mi>&theta;</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>q</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>+</mo><mi>R</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000751479840000011.GIF" wi="1248" he="167" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>&Element;</mo><mi>I</mi></mrow></munder><msub><mi>&gamma;</mi><mi>i</mi></msub><munder><mo>&Sigma;</mo><mrow><msub><mi>x</mi><msup><mi>i</mi><mo>&prime;</mo></msup></msub><mo>&Element;</mo><mi>I</mi><mo>/</mo><mo>{</mo><msub><mi>x</mi><mi>i</mi></msub><mo>}</mo></mrow></munder><mfrac><mrow><mi>&theta;</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>`</mo></mrow></msub><mo>|</mo><mo>|</mo><mo>)</mo></mrow></mrow><mrow><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>`</mo></mrow></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><mn>4</mn><msup><mi>r</mi><mn>2</mn></msup></mrow><msup><mi>h</mi><mn>2</mn></msup></mfrac></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000751479840000012.GIF" wi="1033" he="173" /></maths>其中,i,j为正整数,J为点云采样数据的点集合,q<sub>j</sub>为点集合J中的采样点,I为采样点q<sub>j</sub>的邻域点集合,x<sub>i</sub>为邻域点集合I中的邻域点,x<sub>i'</sub>为邻域点集合I中不属于x<sub>i</sub>的邻域点,R为规整项,γ<sub>i</sub>为权重系数,h为邻域点集合I的邻域半径,σ<sub>i</sub>为分布系数;所述根据迭代公式进行点云收缩的步骤为:设置初始邻域半径<img file="FDA0000751479840000013.GIF" wi="394" he="89" />以h<sub>i</sub>=h<sub>i‑1</sub>+h<sub>0</sub>/2的方式逐步扩大邻域半径,其中,d<sub>bb</sub>为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数;根据扩大的邻域半径以及所述迭代公式进行点云收缩;所述分布系数σ<sub>i</sub>的获取步骤为:根据PCA算法为采样点建立协方差矩阵;根据下述公式得到分布系数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&lambda;</mi><msub><mn>2</mn><mi>i</mi></msub></mrow><mrow><mi>&lambda;</mi><msub><mn>0</mn><mi>i</mi></msub><mo>+</mo><mi>&lambda;</mi><msub><mn>1</mn><mi>i</mi></msub><mo>+</mo><mi>&lambda;</mi><msub><mn>2</mn><mi>i</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000751479840000014.GIF" wi="576" he="142" /></maths>其中,λ0<sub>i</sub>、λ1<sub>i</sub>和λ2<sub>i</sub>为所述协方差矩阵的三个特征值,且λ0<sub>i</sub>≤λ1<sub>i</sub>≤λ2<sub>i</sub>;连接所述骨架分支得到点云骨架。
地址 518055 广东省深圳市南山区西丽大学城学苑大道1068号