主权项 |
一种基于车载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>Δ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><ΔZ′<sub>th</sub>,且ΔZ<sub>2</sub><ΔZ″<sub>th</sub>,判定为地面点,否则为其他类型点;2)若目标点数小于设定的最少点数,或者ΔZ<sub>1</sub><Δ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>Δ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>′</mo></msup></msubsup><mo>=</mo><msub><mi>Z</mi><mi>min</mi></msub><mo>+</mo><msubsup><mi>Δ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>Δ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" />则判定为其他点,否则为建筑物点。 |