车载激光扫描点云中建筑物角点空间位置自动提取方法
【专利摘要】本发明公开了一种车载激光扫描点云中建筑物角点空间位置自动提取方法。首先根据局部建筑物的高度获取局部范围内在竖直方向上的空间分层点云，对空间分层点云进行平面投影，根据空间分层点云的范围与空间分层点云的平面坐标，将离散的激光扫描点云投影到平面；运用平面上单位像素中激光扫描点云的数量定义像素的亮度值，转换成平面投影图像；然后对平面投影图像中亮度较高的区域进行直线检测，提取出具有直线特性的线段，对这些线段分类合并；然后分析提取线段的相对位置关系与相交角度，筛选出相交产生符合建筑物角点特征的线段，并求取线段相交角点在平面投影图像中的坐标，最后根据角点在平面投影图像中的坐标计算角点目标的空间位置信息。
【专利说明】
车载激光扫描点云中建筑物角点空间位置自动提取方法
技术领域
[0001] 本发明属于车载激光扫描点云数据处理技术领域，具体涉及一种车载激光扫描点 云中建筑物角点空间位置自动提取方法。
【背景技术】
[0002] 建筑物角点的测量是地形图、地籍图测量中重要工作之一，目前建筑物角点位置 信息主要通过传统人工测量方法获得，更新周期长、资源损耗大，已经无法满足城市建设和 管理的数据需求。车载移动激光测量系统作为一种先进的测量手段，具有快速、不与测量物 接触、实时、动态、主动、高密度及高精度等特点。基于车载移动激光测量系统的建筑物角点 提取方法主要有投影图像法、特征聚类和扫描线法，然而这些方法没有充分挖掘建筑物点 云数据的形态特征，适应性不高，容易受车载激光扫描点云数据中噪声点的影响。
【发明内容】
[0003] 本发明的目的在于提出一种车载激光扫描点云中建筑物角点空间位置自动提取 方法，以便从海量车载激光扫描点云数据中快速、自动提取建筑物角点的空间位置信息。
[0004] 为了实现上述目的，本发明采用如下技术方案：
[0005] 车载激光扫描点云中建筑物角点空间位置自动提取方法，包括如下步骤：
[0006] a、根据局部建筑物的高度自动获取局部范围内在竖直方向上的空间分层点云数 据，对空间分层点云进行平面投影;根据空间分层点云的范围与空间分层点云的平面坐标 将离散的车载激光扫描点云投影到平面，用单位像素中激光扫描点云的数量定义像素的亮 度值，转换成平面投影图像，并自动获取最优空间分层点云；
[0007] b、最优空间分层点云对生成的平面投影图像进行阈值分割，去掉亮度低的点，即 将最优空间分层点云中平面内投影方向较稀疏的激光扫描点云去掉；
[0008] c、对平面投影图像中亮度较高的区域进行直线检测，提取出有直线特性的线段；
[0009] d、提取出的直线线段分类合并，去掉冗余和不符合建筑物墙壁的线段；
[0010] e、分析这些线段的相对位置关系和相交角度，筛选出相交产生符合建筑物角点特 征的线段，并求这些线段相交产生的角点在平面投影图像中的坐标；
[0011] f、根据角点在平面投影图像中的坐标计算角点目标的空间位置信息。
[0012] 优选地，所述步骤a具体为：
[0013] al、竖直方向空间分层点云自动提取
[0014] 根据实际建筑物空间高度分布对现有的激光扫描点云进行提取：
[0015] 首先对激光扫描点云按照竖直方向整体分层，每5m分一个点云层，获取局部范围 内在竖直方向上的一系列空间分层点云数据；
[0016] 依次对上述空间分层点云进行平面投影、生成基于点数的平面投影图像，即：
[0017] a2、平面投影
[0018] 以激光扫描点云的xy轴对应的平面为投影面，激光扫描点云的z轴方向的负方向 为投影方向，对激光扫描点云进行投影，将建筑物、道路和杆状地物投影到投影坐标系下；
[0019] a3、生成基于点数的平面投影图像
[0020]根据当前激光扫描点云的范围，设任意一个激光扫描点云坐标为(Xp，Yp，Zp)，当 前激光扫描点云范围为{1丨1^^1^此112，1&^，1&^，1&?2}，对应平面投影图像的像素坐 标为(xp，yp)，缩放精度为s;则任意一个激光扫描点云坐标(Xp，Yp，Zp)对应平面投影图像 的像素坐标(xp，yp)分别是:XP = (Xp-MinX) *s，yp = (Yp-MinY) *s;
[0021] 单位像素内点的亮度是落在单位像素内激光扫描点云个数的亮度叠加，设每有一 个激光扫描点云转换成像素点(xp，yp)，则该点亮度值增加30% ;
[0022] 遍历当前所有的激光扫描点云数据，计算其对应的平面投影图像坐标，并进行亮 度叠加，得到基于点数的平面投影图像；
[0023] a4、自动获取最优空间分层点云
[0024]对生成的N个平面投影图像进行粗略的直线检测、块状区域检测和离散点检测，选 取直线分布较多、离散点与块状区域较少的区域作为最优空间分层点云，进行下一步角点 提取。
[0025]优选地，所述步骤d具体为：
[0026] dl、直线线段分类
[0027]首先进行粗分组:遍历所有直线段N1为直线段集合的项数，求其与平 面投影图像坐标X轴正方向的夹角ctk，其中Ctk单位是度％在直线段集合L中任取一条直线段 lx，与直线段集合L中的其他直线段I k的X轴方向的夹角αχ做比较，其中(^单位也是度％ Kk 彡Ni，k乒X，如果有|ax-ak| <3°，且Ix与Ik的端点的最小距离小于l*s，则将Ik与Ix归为一组， 遍历完毕得到一组新的包括I x在内的夹角相差3°内的直线段集合，设新的直线段集合为Lx; 直线段集合L = L-Lx,继续在直线段集合L中取一条新的ly，Ky彡Ni-Len(Lx)，其中Len(L x) 为集合Lx的项数;重复上述过程得到一组新的包括Iy在内的夹角相差3°内且端点相聚较近 的直线段集合L y;按照上述过程对直线段集合L中剩余的直线段做分组，得到一系列的直线 段集合:Ll，L2，L3 · · · LN2 ;
[0028] d2、对上述直线段组集合的每一项Lj进一步细分组，1彡j SN2:计算Lj中直线段与 平面投影图像坐标X轴正方向的夹角的平均值比;在直线段集合L j中任取一条1#，计算该线 段的中点P如遍历L遇内的其他直线段1~计算线段的中点p如基于咖和A构建一条直线 lj U'，计算点Pjx与1扣'的垂距djU;设置垂距阈值dth = 0.25*s;如果有dju<dth，则将Iju与Ijx归 为一组，遍历完毕得到一组包括I jx在内的垂距小于阈值的直线段集合，设该直线段集合为 Ljx;则剩余直线段集合Lj = Lj-Ljx;继续在直线段集合L冲取一条新的ljy，重复上述过程得 到一组新的包括I jy在内的垂距小于垂距阈值的直线段集合Ljy;按照上述过程对Lj中剩余的 直线段做分组，最终将Lj转换为一系列的直线段集合:Ljl，Lj2，Lj3. . .Ljm;
[0029] Lj的每一组按照上述的细分组的方法转换为直线段集合Li，1彡u'彡m'，m'为L JU， 的项数；
[0030] d3、将每一组直线段进行单一化处理，转为一条直线段；
[0031 ]遍历Lw的每一项，即一组直线段:计算Lw中直线段斜率的平均值5^ ;
[0032]如果I Sv I < 1，则基于X轴进行单一化处理:遍历Lv中所有直线段的端点坐标P (x，y)，获得X轴坐标值最小的坐标Pw min和X轴坐标值最大的坐标Pju，max，创建为唯一的直线 段 lju，' ；
[0033] 否则基于Y轴进行单一化处理:遍历LW中所有直线段的端点坐标p(x，y)，获得Y轴 坐标值最小的坐标PjYmin和Y轴坐标值最大的坐标PjYmax，创建为唯一直线段lju〃 ；
[0034] 最终得到单一化处理后的唯一直线段Iw'组成的集合L'。
[0035]优选地，所述步骤e具体为：
[0036] el、在直线段集合1/中任取一条直线段lx'，依次与直线段集合L'中的其他直线段 Γ做比较：如果两直线段对应直线的夹角3>^^，其中，δ为锐角或直角，S th为垂直夹角阈 值，且lx'与Γ的端点的最小距离小于指定距离dth，则计算两直线段对应直线的交点ρ'；遍 历完毕得到一组交点集合P x将lx'从集合L'中移除;在新的集合L'中任取一条直线段1/，重 复上述过程计算得到一组交点集合P y;按照上述过程继续对1/中剩余的直线段进行判断和 求交，得到一系列的交点组的集合Pi，P2，P 3.. .Pn;
[0037] e2、将交点组的集合P1J2J3. ..Pr^组的点Pr组成唯一的点集合P'，1彡r彡η。
[0038] 本发明具有如下优点：
[0039]本发明方法根据激光扫描点云数据，自动获取空间分层数据，并自动选取最优空 间分层数据进行平面投影，通过分析杆建筑物与角点的特征，对基于点数的平面投影图像 进行直线的提取，从而获取角点的空间位置信息，不易受数据噪声点的影响，自动化程度 高，更大程度上充分利用了激光扫描点云的形态特征，达到了较好的提取效果。
【附图说明】
[0040] 图1为本发明中车载激光扫描点云中建筑物角点空间位置自动提取方法的流程示 意图。
【具体实施方式】
[0041] 本发明的基本思想是:通过研究点车载激光扫描点云生成的基于点数的平面投影 图像，运用图像处理的方式，结合建筑物墙壁与角点的几何形态与特征，提出一种车载激光 扫描点云中建筑物角点空间位置自动提取方法。
[0042]结合图1所示，车载激光扫描点云中建筑物角点空间位置自动提取方法包括如下 步骤：
[0043] a、根据局部建筑物的高度自动获取局部范围内在竖直方向上的空间分层点云数 据，对空间分层点云进行平面投影;根据空间分层点云的范围与空间分层点云的平面坐标 将离散的车载激光扫描点云投影到平面，用单位像素中激光扫描点云的数量定义像素的亮 度值，转换成平面投影图像，并自动获取最优空间分层点云；
[0044]该步骤a具体为：
[0045] al、竖直方向空间分层点云自动提取
[0046]根据实际建筑物空间高度分布对现有的激光扫描点云进行提取：
[0047]首先对激光扫描点云按照竖直方向整体分层，每5m分一个点云层，获取局部范围 内在竖直方向上的一系列空间分层点云数据；
[0048]依次对上述空间分层点云进行平面投影、生成基于点数的平面投影图像，即：
[0049] a2、平面投影
[0050] 以激光扫描点云的xy轴对应的平面为投影面，激光扫描点云的Z轴方向的负方向 为投影方向，对激光扫描点云进行投影，将建筑物、道路和杆状地物投影到投影坐标系下；
[0051] a3、生成基于点数的平面投影图像
[0052]根据当前激光扫描点云的范围，设任意一个激光扫描点云坐标为(Xp，Yp，Zp)，当 前激光扫描点云范围为{1丨1^^1^此112，1&^，1&^，1&?2}，对应平面投影图像的像素坐 标为(xp，yp)，缩放精度为s;则任意一个激光扫描点云坐标(Xp，Yp，Zp)对应平面投影图像 的像素坐标(xp，yp)分别是:XP = (Xp-MinX) *s，yp = (Yp-MinY) *s;
[0053]单位像素内点的亮度是落在单位像素内激光扫描点云个数的亮度叠加，设每有一 个激光扫描点云转换成像素点(xp，yp)，则该点亮度值增加30% ;
[0054]遍历当前所有的激光扫描点云数据，计算其对应的平面投影图像坐标，并进行亮 度叠加，得到基于点数的平面投影图像；
[0055] a4、自动获取最优空间分层点云
[0056] 对生成的N个平面投影图像进行粗略的直线检测、块状区域检测和离散点检测，选 取直线分布较多、离散点与块状区域较少的区域作为最优空间分层点云，进行下一步角点 提取。
[0057] b、最优空间分层点云对生成的平面投影图像进行阈值分割，去掉亮度低的点，即 将最优空间分层点云中平面内投影方向较稀疏的激光扫描点云去掉；
[0058] 该步骤b具体为：
[0059]采用自适应阈值的方法进行图像分割，通过计算像素点周围区域的加权平均，然 后减去一个常数来得到自适应阈值；
[0060] 对得到的平面投影图像进行阈值分割得到亮度较大的区域，该区域即对应激光扫 描点云中垂直于Xy平面的单位范围内竖直方向有较多点的数据，符合建筑物墙壁竖直方向 上有较多点的特征。
[0061] c、对平面投影图像中亮度较高的区域进行直线检测，提取出有直线特性的线段；
[0062] 该步骤c具体为：
[0063] cl、图像细化
[0064]对阈值分割后的平面投影图像进行图像细化处理，排除掉同一位置多条直线对角 点提取的影响；
[0065] c2、直线检测
[0066] 对阈值分割后的平面投影图像进行直线检测，这些直线有可能是建筑物墙壁，记 录下这些直线段的端点信息，设直线段的集合为L，直线段分别为h、l 2、l3……InJPLel1, I2,13 . . . Ino
[0067] d、提取出的直线线段分类合并，去掉冗余和不符合建筑物墙壁的线段；
[0068] 该步骤d具体为：
[0069] dl、直线线段分类
[0070]首先进行粗分组:遍历所有直线段N1为直线段集合的项数，求其与平 面投影图像坐标X轴正方向的夹角ak，其中ak单位是度％在直线段集合L中任取一条直线段 lx，与直线段集合L中的其他直线段I k的X轴方向的夹角αχ做比较，其中(^单位也是度％ Kk 彡Ni，k乒X，如果有|ax-ak| <3°，且Ix与Ik的端点的最小距离小于l*s，则将Ik与Ix归为一组， 遍历完毕得到一组新的包括Ix在内的夹角相差3°内的直线段集合，设新的直线段集合为Lx; 直线段集合L = L-Lx,继续在直线段集合L中取一条新的ly，Ky彡Ni-Len(Lx)，其中Len(L x) 为集合Lx的项数;重复上述过程得到一组新的包括Iy在内的夹角相差3°内且端点相聚较近 的直线段集合L y;按照上述过程对直线段集合L中剩余的直线段做分组，得到一系列的直线 段集合：- L: · H、、；
[0071] d2、对上述直线段组集合的每一项Lj进一步细分组，I < jSN2:计算Lj中直线段与 平面投影图像坐标X轴正方向的夹角的平均值比;在直线段集合Lj中任取一条1#，计算该线 段的中点P如遍历L遇内的其他直线段1~计算线段的中点p如基于咖和A构建一条直线 lj U'，计算点Pjx与1扣'的垂距djU;设置垂距阈值dth = 0.25*s;如果有dju<dth，则将Iju与Ijx归 为一组，遍历完毕得到一组包括I jx在内的垂距小于阈值的直线段集合，设该直线段集合为 Ljx;则剩余直线段集合Lj = Lj-Ljx;继续在直线段集合L冲取一条新的ljy，重复上述过程得 到一组新的包括I jy在内的垂距小于垂距阈值的直线段集合Ljy;按照上述过程对Lj中剩余的 直线段做分组，最终将Lj转换为一系列的直线段集合:Ljl，Lj2，Lj3. . .Ljm;
[0072] Lj的每一组按照上述的细分组的方法转换为直线段集合，1彡u'彡m'，m'为LJU， 的项数；
[0073] d3、将每一组直线段进行单一化处理，转为一条直线段；
[0074] 遍历Lw的每一项，即一组直线段:计算Lw中直线段斜率的平均值SJU，；
[0075] 如果I Si I <1，则基于平面投影图像坐标X轴进行单一化处理：
[0076] 遍历Lju1中所有直线段的端点坐标p(x，y)，获得X轴坐标值最小的坐标PjU，min和X轴 坐标值最大的坐标PjYmax，仓Il建为唯一的直线段lju〃 ；
[0077] 否则基于Y轴进行单一化处理:遍历LW中所有直线段的端点坐标p(x，y)，获得Y轴 坐标值最小的坐标PjYmin和Y轴坐标值最大的坐标PjYmax，创建为唯一直线段lju〃 ；
[0078] 最终得到单一化处理后的唯一直线段Iw '组成的集合L'。
[0079] e、分析这些线段的相对位置关系和相交角度，筛选出相交产生符合建筑物角点特 征的线段，并求这些线段相交产生的角点在平面投影图像中的坐标；
[0080] 该步骤e具体为：
[0081] el、在直线段集合1/中任取一条直线段lx'，依次与直线段集合L'中的其他直线段 Γ做比较：如果两直线段对应直线的夹角3>^^，其中，δ为锐角或直角，Sth为垂直夹角阈 值，且l x'与Γ的端点的最小距离小于指定距离dth，则计算两直线段对应直线的交点ρ'；遍 历完毕得到一组交点集合P x将lx'从集合L'中移除;在新的集合L'中任取一条直线段1/，重 复上述过程计算得到一组交点集合P y;按照上述过程继续对L'中剩余的直线段进行判断和 求交，得到一系列的交点组的集合Pi，P2，P 3.. .Pn;
[0082] e2、将交点组的集合P1J2J3. ..Pr^组的点Pr组成唯一的点集合P'，1彡r彡η。
[0083] f、根据角点在平面投影图像中的坐标计算角点目标的空间位置信息。
[0084]该步骤f具体为：
[0085] 建筑物角点坐标设为Pr (xr，yr)，则对应的激光扫描点云坐标(Xr，Yr，Zr)为:Xr = xr/s+MinX，Yr = yr/s+MinY，Z轴方向取对激光扫描点云进行空间分层时Z方向的中间值;求 出的激光扫描点云坐标(Xr，Yr，Zr)即为建筑物角点的空间位置。
[0086]当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应 当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明 显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
【主权项】
1. 车载激光扫描点云中建筑物角点空间位置自动提取方法，其特征在于，包括如下步 骤： a、 根据局部建筑物的高度自动获取局部范围内在竖直方向上的空间分层点云数据，对 空间分层点云进行平面投影;根据空间分层点云的范围与空间分层点云的平面坐标将离散 的车载激光扫描点云投影到平面，用单位像素中激光扫描点云数量定义像素的亮度值，转 换成平面投影图像，并自动获取最优空间分层点云； b、 最优空间分层点云对生成的平面投影图像进行阈值分割，去掉亮度低的点，即将最 优空间分层点云中平面内投影方向较稀疏的激光扫描点云去掉； c、 对平面投影图像中亮度较高的区域进行直线检测，提取出有直线特性的线段； d、 提取出的直线线段分类合并，去掉冗余和不符合建筑物墙壁的线段； e、 分析这些线段的相对位置关系和相交角度，筛选出相交产生符合建筑物角点特征的 线段，并求这些线段相交产生的角点在平面投影图像中的坐标； f、 根据角点在平面投影图像中的坐标计算角点目标的空间位置信息。2. 根据权利要求1所述的车载激光扫描点云中建筑物角点空间位置自动提取方法，其 特征在于，所述步骤a具体为： al、竖直方向空间分层点云自动提取 根据实际建筑物空间高度分布对现有的激光扫描点云进行提取： 首先对激光扫描点云按照竖直方向整体分层，每5m分一个点云层，获取局部范围内在 竖直方向上的一系列空间分层点云数据； 依次对上述空间分层点云进行平面投影、生成基于点数的平面投影图像，即： a2、平面投影 以激光扫描点云的xy轴对应的平面为投影面，激光扫描点云的z轴方向的负方向为投 影方向，对激光扫描点云进行投影，将建筑物、道路和杆状地物投影到投影坐标系下； a3、生成基于点数的平面投影图像 根据当前激光扫描点云的范围，设任意一个激光扫描点云坐标为(Xp，Yp，Zp)，当前激 光扫描点云范围为{1丨1^此1^^112，1^^，1&^，1&?2}，对应平面投影图像的像素坐标为 (xp，yp)，缩放精度为s;则任意一个激光扫描点云坐标(Χρ，Υρ，Ζρ)对应平面投影图像的像 素坐标(xp，yp)分别是：xp = (Xp-MinX)*s，yp = (Yp-MinY)*s; 单位像素内点的亮度是落在单位像素内激光扫描点云个数的亮度叠加，设每有一个激 光扫描点云转换成像素点(xp，yp)，则该点亮度值增加30% ; 遍历当前所有的激光扫描点云数据，计算其对应的平面投影图像坐标，并进行亮度叠 加，得到基于点数的平面投影图像； a4、自动获取最优空间分层点云 对生成的N个平面投影图像进行粗略的直线检测、块状区域检测和离散点检测，选取直 线分布较多、离散点与块状区域较少的区域作为最优空间分层点云，进行下一步角点提取。3. 根据权利要求2所述的车载激光扫描点云中建筑物角点空间位置自动提取方法，其 特征在于，所述步骤d具体为： dl、直线线段分类 首先进行粗分组:遍历所有直线段lk，1 SkSNi，他为直线段集合的项数，求其与平面投 影图像坐标X轴正方向的夹角ak，其中Qk单位是度％在直线段集合L中任取一条直线段l x，与 直线段集合L中的其他直线段lk的X轴方向的夹角αχ做比较，其中αχ单位是度° 乒X，如果有|ax-ak| <3°，且lx与lk的端点的最小距离小于l*s，则将lk与lx归为一组，遍历完 毕得到一组新的包括lx在内的夹角相差3°内的直线段集合，设新的直线段集合为Lx;直线段 集合L = L-Lx，继续在直线段集合L中取一条新的1 y，1彡y彡Ni-Len (Lx)，其中Len (Lx)为集合 Lx的项数;重复上述过程得到一组新的包括ly在内的夹角相差3°内且端点相聚较近的直线 段集合L y;按照上述过程对直线段集合L中剩余的直线段做分组，得到一系列的直线段集 合：Ll，L2，L3 · · · LN2 ; d2、对上述直线段组集合的每一项b进一步细分组，计算b中直线段与平面 投影图像坐标X轴正方向的夹角的平均值比;在直线段集合b中任取一条1#，计算该线段的 中点Pjx;遍历Lj组内的其他直线段1 jU，计算线段的中点pjU;基于pjU和β」构建一条直线1 jU'， 计算点Pjx与ljU'的垂距djU;设置垂距阈值dth = 0.25*s;如果有djU<dth，则将ljU与ljx归为一 组，遍历完毕得到一组包括1#在内的垂距小于阈值的直线段集合，设该直线段集合为L JX; 则剩余直线段集合b = b-Ljx;继续在直线段集合k中取一条新的1 ，重复上述过程得到一 组新的包括1#在内的垂距小于垂距阈值的直线段集合Ljy;按照上述过程对b中剩余的直线 段做分组，最终将Lj转换为一系列的直线段集合:Ljl，Lj2，Lj3. . .Ljm; Lj的每一组按照上述细分组的方法转换为直线段集合LW，m'为LW的项数； d3、将每一组直线段进行单一化处理，转为一条直线段； 遍历Lw的每一项，即一组直线段:计算Lw中直线段斜率的平均值Sw ; 如果|SW|<1，则基于X轴进行单一化处理:遍历LW中所有直线段的端点坐标p(x，y)， 获得X轴坐标值最小的坐标Pju，min和X轴坐标值最大的坐标Pju， max，创建为唯一的直线段 1JU"; 否则基于Y轴进行单一化处理:遍历LW中所有直线段的端点坐标p(x，y)，获得Y轴坐标 值最小的坐标PjYmin和Y轴坐标值最大的坐标PjYmax，创建为唯一直线段lju〃 ； 最终得到单一化处理后的唯一直线段i组成的集合L'。4.根据权利要求3所述的车载激光扫描点云中建筑物角点空间位置自动提取方法，其 特征在于，所述步骤e具体为： el、在直线段集合1/中任取一条直线段lx'，依次与直线段集合L'中的其他直线段Γ做 比较:如果两直线段对应直线的夹角δ>δ*，其中，δ为锐角或直角，Sth为垂直夹角阈值，且 lx'与Γ的端点的最小距离小于指定距离dth，则计算两直线段对应直线的交点p'；遍历完毕 得到一组交点集合P x将lx'从集合L'中移除;在新的集合L'中任取一条直线段1/，重复上述 过程计算得到一组交点集合P y;按照上述过程继续对1/中剩余的直线段进行判断和求交， 得到一系列的交点组的集合Pi，P2，P 3.. .Pn; e2、将交点组的集合PiPlPs. . .?"各组的点Pr组成唯一的点集合P'，1彡r彡η。
【文档编号】G06T7/00GK106056659SQ201610368081
【公开日】2016年10月26日
【申请日】2016年5月27日
【发明人】刘如飞, 田茂义, 卢秀山, 刘甜, 曲杰卿, 侯海龙
【申请人】山东科技大学