发明名称 基于车载LiDAR点云数据的分类方法
摘要 本发明公开了一种基于车载LiDAR点云数据的分类方法,包括步骤:获取街景的车载LiDAR点云数据,并对其进行预处理,以去除点云数据中冗余数据和噪声,在预处理过程中计算点云的法向量、曲率和密度;对经过预处理的点云数据进行分类,根据点云的法向量、曲率和密度提取电力线、建筑物立面点以及树木点。本发明的分类方法简单,计算量小,且快捷有效,为街景面片的提取做好了铺垫,使得面片提取更加精准。
申请公布号 CN104463872A 申请公布日期 2015.03.25
申请号 CN201410752056.2 申请日期 2014.12.10
申请人 武汉大学 发明人 姚剑;陈梦怡;李礼;鲁小虎
分类号 G06T7/00(2006.01)I;G06T5/00(2006.01)I;G06K9/62(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 湖北武汉永嘉专利代理有限公司 42102 代理人 许美红;唐万荣
主权项 一种基于车载LiDAR点云数据的分类方法,其特征在于,包括步骤:S1、获取街景的车载LiDAR点云数据,并对其进行预处理,以去除点云数据中冗余数据和噪声,在预处理过程中计算点云的法向量、曲率和密度;S2、对经过预处理的点云数据进行分类,根据得到的点云的法向量、曲率和密度提取电力线、建筑物立面点以及树木点,其中提取建筑物立面点包括以下步骤:S21、找到点云数据在XOY平面上的最小矩形外包框(BoundingBox);S22、对全部点云数据在XOY平面上格网化,格网宽高求解公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><mi>ceil</mi><mrow><mo>(</mo><mfrac><msub><mi>box</mi><mi>x</mi></msub><msub><mi>GSD</mi><mi>xy</mi></msub></mfrac><mo>+</mo><mn>2.01</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000629367210000011.GIF" wi="537" he="126" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mi>ceil</mi><mrow><mo>(</mo><mfrac><msub><mi>box</mi><mi>y</mi></msub><msub><mi>GSD</mi><mi>xy</mi></msub></mfrac><mo>+</mo><mn>2.01</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000629367210000012.GIF" wi="527" he="128" /></maths>式中,box<sub>x</sub>,box<sub>y</sub>是外包框的x,y方向长度,GSD<sub>xy</sub>为格网间距;S23、将全部点云数据与格网建立对应的索引关系,分离地面点和非地面点,包括以下两种情况:1)找到最小高程值Z<sub>min</sub>和最大高程值Z<sub>max</sub>,在Z方向上进行修正,得到修正后的最小高程值<img file="FDA0000629367210000013.GIF" wi="104" he="76" />和最大高程值<img file="FDA0000629367210000014.GIF" wi="133" he="74" />设定经验阈值ΔZ′<sub>th</sub>,ΔZ″<sub>th</sub>,ΔZ<sub>1</sub>=Z<sub>max</sub>‑Z<sub>min</sub>,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&Delta;Z</mi><mn>2</mn></msub><mo>=</mo><msub><mi>Z</mi><mi>min</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>min</mi><mo>^</mo></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000629367210000017.GIF" wi="428" he="93" /></maths>若ΔZ<sub>1</sub>&lt;ΔZ′<sub>th</sub>,且ΔZ<sub>2</sub>&lt;ΔZ″<sub>th</sub>,判定为地面点,否则为其他类型点;2)若目标点数小于设定的最少点数,或者ΔZ<sub>1</sub>&lt;ΔZ″<sub>th</sub>,<img file="FDA0000629367210000015.GIF" wi="145" he="90" />为设定的地面点高程最大值阈值,<img file="FDA0000629367210000016.GIF" wi="281" he="93" />两种基准的地面点最大高程值,计算公式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>Z</mi><mi>max</mi><mi>flr</mi></msubsup><mo>=</mo><msubsup><mi>z</mi><mi>min</mi><mo>^</mo></msubsup><mo>+</mo><msubsup><mi>&Delta;Z</mi><mi>max</mi><mi>flr</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA00006293672100000211.GIF" wi="570" he="119" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>Z</mi><mi>max</mi><msup><mi>flr</mi><mo>&prime;</mo></msup></msubsup><mo>=</mo><msub><mi>Z</mi><mi>min</mi></msub><mo>+</mo><msubsup><mi>&Delta;Z</mi><mi>max</mi><mi>flr</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000629367210000021.GIF" wi="497" he="112" /></maths>若,某点高程值<img file="FDA0000629367210000022.GIF" wi="258" he="84" />则为地面点;若<img file="FDA0000629367210000023.GIF" wi="234" he="90" />利用点的曲率和法向量,若<img file="FDA0000629367210000024.GIF" wi="236" he="92" />且<img file="FDA0000629367210000025.GIF" wi="285" he="91" /><img file="FDA0000629367210000026.GIF" wi="99" he="90" />和<img file="FDA0000629367210000027.GIF" wi="135" he="92" />为设定的曲率阈值和法向量变化阈值,R<sub>i</sub>为某点的曲率,n<sub>i</sub>为对应的法向量;否则点被判定为建筑物点;S24、步骤S23中将点分类为非地面点和地面点,非地面点里也分离出了一部分建筑物立面点,现对非地面点分析找到所有的建筑物立面点:将非地面点云数据格网化,重新建立与格网对应的索引关系,由于建筑物立面面积较大,可划分为若干小区域立方体(blob),设置blob的大小阈值为S<sub>th</sub>,建筑物立面的高程差值为<img file="FDA0000629367210000028.GIF" wi="134" he="77" /><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&Delta;Z</mi><mi>max</mi></msub><mo>=</mo><msub><mi>Z</mi><mi>max</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>min</mi><mo>^</mo></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000629367210000029.GIF" wi="489" he="78" /></maths>若<img file="FDA00006293672100000210.GIF" wi="323" he="71" />则判定为其他点,否则为建筑物点。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学