一种从地面LiDAR数据中提取建筑物轮廓和角点的方法.pdf

(19)中华人民共和国国家知识产权局*CN103020342A*(12)发明专利申请(10)申请公布号CN103020342A(43)申请公布日2013.04.03(21)申请号201210512462.2(22)申请日2012.12.04(71)申请人南京大学地址210093江苏省南京市鼓楼区汉口路22号(72)发明人李满春程亮童礼华陈焱明刘永学王结臣钟礼山张雯陈小雨孙越凡(74)专利代理机构南京同泽专利事务所(特殊普通合伙)32245代理人石敏(51)Int.Cl.G06F17/50(2006.01)权利要求书权利要求书2页2页说明书说明书66页页附图附图66页(54)发明名称一种从地面LiDAR数据中提取建筑物轮廓和角点的方法(57)摘要本发明涉及一种从地面LiDAR数据中提取建筑物轮廓和角点的方法，首先使用分层次的格网密度方法从地面LiDAR数据中提取建筑物轮廓；在此基础上使用轮廓延伸密度方法对提取的建筑物轮廓进行恢复，即得到完整的建筑物轮廓；若需提取角点则将完整的建筑物轮廓投影到三维坐标系的XY平面内寻找二维相交点，如果任两条构成相交点的轮廓的高程差小于1m，则判定两条轮廓在实际的三维空间中相交，两条轮廓的相交点为一个地面角点，并将所述两条轮廓的高程均值作为该地面角点的高程。本发明所用的格网密度方法、格网密度阈值的理论估计确定方法、轮廓密度延伸的方法，保证了从地面LiDAR数据中提取准确的建筑物轮廓线段和高精度的地面角点；并且实现了自动化提取，大大提高了数据处理效率。CN10324ACN103020342A权利要求书1/2页1.一种从地面LiDAR数据中提取建筑物轮廓的方法，包括如下步骤：1a)地面LiDAR点云投影至XY平面——将地面LiDAR点云投影至三维坐标系的XY平面，并保留各个点的X、Y、Z属性；1b)提取粗略轮廓格网——构建1m*1m的粗略格网，计算每个粗略格网中LiDAR投影点的数量，即得到该粗略格网的格网密度，根据建筑物边缘轮廓处的粗略格网密度阈值对所述粗略格网进行筛选，保留格网密度大于所述粗略格网密度阈值的粗略格网，得到粗略轮廓格网；1c)提取精确轮廓格网——在粗略轮廓格网中构建0.2m*0.2m的精细格网，计算精细格网内LiDAR投影点的数量即得到精细格网的格网密度，根据建筑物边缘轮廓处的精细格网密度阈值对所述精细格网进行筛选，保留格网密度大于所述精细格网密度阈值的精细格网，得到精确轮廓格网；1d)格网高差筛选——遍历所有精确轮廓格网，如果精确轮廓格网内的最高LiDAR点和最低LiDAR点的高差大于相应实验区建筑最低高程则保留该精确轮廓格网，否则剔除；1e)获取轮廓线段——对筛选后的精确轮廓格网使用Hough变换得到二维矢量轮廓线段；1f)二维轮廓转化为三维轮廓——寻找步骤1e)中获得的二维矢量轮廓线段周边1m范围内格网，将寻找到的所有格网内LiDAR点最大高程的平均值作为二维矢量轮廓线段的高程，将二维矢量轮廓线段变换为三维建筑物轮廓线段。2.根据权利要求1所述一种从地面LiDAR数据中提取建筑物轮廓的方法，其特征在于：步骤1b)和1c)中的格网密度阈值的确定方法如下：假设O点为仪器中心点，A点为水平垂直于仪器的墙面点，扫描仪对准A点时的角度为0°，B点为格网靠近仪器一侧，C点为格网远离仪器一侧，D点为B点竖直方向上墙面最高点，设OA＝DV，CO＝DM，水平方向格网的边长为DG，建筑高HB，仪器高HL，在A点处水平向相邻两LiDAR点的间距为DR，则格网密度计算方法如下：2a)计算水平方向格网内LiDAR点的列数，记α为扫描仪每次旋转角度的一半，记格网中水平方向上最靠近于B点的角度为β，则线段BC上LiDAR点数为则水平方向格网内LiDAR点的列数为Ncol；2b)计算每一列LiDAR点的数量，第i列LiDAR点的数量为2c)将每列LiDAR点数量相加得到格网处LiDAR点总数网格密度阈值threshod＝rate*N，其中参数rate为描述墙面凹凸、窗户多少的阈值，墙面凹凸越多、窗户2CN103020342A权利要求书2/2页越多，则该参数越小，rate的取值范围为0.2-1。3.根据权利要求3所述一种从地面LiDAR数据中提取建筑物轮廓的方法，其特征在于：步骤1f)完成后，利用轮廓延伸密度方法进行建筑物轮廓的恢复，具体步骤如下：3a)对三维建筑物轮廓线段构建半径为1m的建缓冲区，建缓冲区内LiDAR点数量除以缓冲区体积获得原有轮廓LiDAR点密度；3b)沿轮廓线段方向以单位距离为延伸步长构建半径为1m的缓冲区，缓冲区内LiDAR点数量除以相应缓冲区体积获得待延伸方向的LiDAR点密度，所述单位距离的取值范围为0.1-0.3m；3c)若待延伸方向的LiDAR点密度与原有轮廓LiDAR点密度的差异