CN105913437A - 一种道路完好性检测装置 - Google Patents

一种道路完好性检测装置 Download PDF

Info

Publication number
CN105913437A
CN105913437A CN201610230945.1A CN201610230945A CN105913437A CN 105913437 A CN105913437 A CN 105913437A CN 201610230945 A CN201610230945 A CN 201610230945A CN 105913437 A CN105913437 A CN 105913437A
Authority
CN
China
Prior art keywords
image
point
points
module
line segment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610230945.1A
Other languages
English (en)
Inventor
吴本刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201610230945.1A priority Critical patent/CN105913437A/zh
Publication of CN105913437A publication Critical patent/CN105913437A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30248Vehicle exterior or interior
    • G06T2207/30252Vehicle exterior; Vicinity of vehicle
    • G06T2207/30256Lane; Road marking

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明一种道路完好性检测装置,包括图像预处理模块、图像极值点检测模块、图像特征点定位模块、主方向确定模块和特征提取模块,其中所述图像特征点定位模块通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,所述主方向确定模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,并将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段后将多条线段中的最优线段的方向作为特征点的主方向。本发明具有识别检测精度高、速度快的优点。

Description

一种道路完好性检测装置
技术领域
本发明涉及交通领域,具体涉及一种道路完好性检测装置。
背景技术
对于大部分司机来说,道路的完好性意味着行驶的舒适度和安全度,对于无人车来说,道路的完好性检测意味着其能否完成行驶,人类对道路完好性检测的研究从未间断,很多道路完好性检测装置被制造出来。此外,为了对大规模的数据进行处理,需要提高分析处理效率和精度。
发明内容
针对上述问题,本发明提供一种识别检测速度快精度高的道路完好性检测装置。
本发明的目的采用以下技术方案来实现:
提供了一种道路完好性检测装置,用于对道路图像进行识别检测,包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 [ m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) ]
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , σ ) = 1 2 π σ e - x 2 / 2 σ 2 , G ( y , σ ) = 1 2 π σ e - y 2 / 2 σ 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , σ ) + ∂ D ( x , y , σ ) T ∂ x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本发明的有益效果为:
1、设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;
2、提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;
3、设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;
4、设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性。
附图说明
利用附图对本发明作进一步说明,但附图中的实施例不构成对本发明的任何限制,对于本领域的普通技术人员,在不付出创造性劳动的前提下,还可以根据以下附图获得其它的附图。
图1是本发明各模块的连接示意图。
具体实施方式
结合以下实施例对本发明作进一步描述。
实施例1
参见图1,本实施例道路完好性检测装置包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本实施例设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性;本实施例取阈值T1=0.01,T2=10,T3=0.1,道路图像分析检测的精度提高了2%,速度提高了1%。
实施例2
参见图1,本实施例道路完好性检测装置包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本实施例设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性;本实施例取阈值T1=0.02,T2=11,T3=0.08,道路图像分析检测的精度提高了1%,速度提高了1.5%。
实施例3
参见图1,本实施例道路完好性检测装置包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本实施例设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性;本实施例取阈值T1=0.03,T2=12,T3=0.06,道路图像分析检测的精度提高了2.5%,速度提高了3%。
实施例4
参见图1,本实施例道路完好性检测装置包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本实施例设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性;本实施例取阈值T1=0.04,T2=13,T3=0.04,道路图像分析检测的精度提高了1.5%,速度提高了2%。
实施例5
参见图1,本实施例道路完好性检测装置包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定;
优选地,所述道路完好性检测装置还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; m a x , g &OverBar; m a x = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
进一步地,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
本实施例设置的图像预处理模块考虑了视觉习惯以及人眼对不同色彩的感知度同色彩强度的非线性关系,能够更为准确的描述图像;提出了高斯差分尺度空间的简化计算公式,减少了运算量,提高了运算速度,进而提高了图像分析的速度;设置的图像特征点定位模块对极值点进行低对比度点和不稳定边缘点的去除,保证特征点的有效性,其中对图像的灰度值进行增强,能够大大增加图像的稳定性,更为精确的对低对比度点进行去除,进而提高图像分析的准确度;设置主方向确定模块,提出了最优线段的判定公式,以特征点梯度方向直方图中的任意相邻两个峰值连线形成的线段中的最优线段的方向作为特征点的主方向,线段相对于点更加稳定,使得图像对应的特征点的描述符具有可重复性,提高了特征描述符的准确性,进而能够更加快速、准确的对道路图像进行识别检测,具有很高的鲁棒性;本实施例取阈值T1=0.05,T2=14,T3=0.02,道路图像分析检测的精度提高了1.8%,速度提高了1.5%。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案,而非对本发明保护范围的限制,尽管参照较佳实施例对本发明作了详细地说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的实质和范围。

Claims (3)

1.一种道路完好性检测装置,用于对道路图像进行识别检测,其特征是,包括:
(1)图像预处理模块,其包括用于将彩色图像转化为灰度图像的图像转换子模块和用于对所述灰度图像进行滤波的图像滤波子模块,所述图像转换子模块的图像灰度转换公式为:
I ( x , y ) = m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) + m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) 2 + 2 &lsqb; m a x ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) - m i n ( R ( x , y ) , G ( x , y ) , B ( x , y ) ) &rsqb;
其中,R(x,y)、G(x,y)、B(x,y)分别代表像素点(x,y)处的红绿蓝强度值,I(x,y)代表像素点(x,y)处的灰度值;
(2)图像极值点检测模块,其通过由高斯差分算子与图像进行卷积建立成的图像的高斯差分尺度空间来检测各极值点的位置,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆大时,所述采样点为极大值点,当采样点相对于与它同尺度的8个相邻点和上下相邻尺度对应的18个点的值皆小时,所述采样点为极小值点,所述高斯差分尺度空间的简化计算公式为:
D(x,y,σ)=(G(x,kσ)-G(x,σ))*I'(x,y)+(G(y,kσ)-G(y,σ))*I'(x,y)
此处
G ( x , &sigma; ) = 1 2 &pi; &sigma; e - x 2 / 2 &sigma; 2 , G ( y , &sigma; ) = 1 2 &pi; &sigma; e - y 2 / 2 &sigma; 2
其中,D(x,y,σ)表示高斯差分尺度空间函数,I'(x,y)为由图像转化子模块输出的图像函数,*表示卷积运算,σ表示尺度空间因子,G(x,σ)、G(y,σ)为定义的尺度可变的高斯函数,k为不变倍增因子;
(3)图像特征点定位模块,其通过剔除所述各极值点中对噪声敏感的低对比度点以及不稳定的边缘点来确定作为特征点的极值点,包括依次连接的用于极值点精确定位的第一定位子模块、用于去除低对比度点的第二定位子模块和用于去除不稳定边缘点的第三定位子模块,其中:
a、所述第一定位子模块通过对所述高斯差分尺度空间函数进行二次泰勒展开并求导得到极值点的精确位置,极值点的尺度空间函数为:
D ( X ^ ) = D ( x , y , &sigma; ) + &part; D ( x , y , &sigma; ) T &part; x X ^
其中,表示极值点的尺度空间函数,D(x,y,σ)T为相对于极值点的偏移量,表示极值点的精确位置;
b、所述第二定位子模块对由图像转换子模快输出的图像依次进行灰度增强、归一化处理后剔除所述低对比度点,增强后的灰度值为:
此处
所述低对比度点的判定公式为:
D ( X ^ ) < T 1 , T 1 &Element; &lsqb; 0.01 , 0.06 &rsqb;
其中,I”(x,y)表示灰度值增强后的图像函数,为包含局部信息的校正系数,M为像素的最大灰度值,所述最大灰度值M=255,mH为图像中灰度值高于128的所有像素的均值,mL是灰度值低于128的所有像素的均值,ψ(x,y)为由图像滤波子模块处理后的图像,T1为设定的阈值;
c、所述第三定位子模块通过计算极值点的位置尺度为2×2的Hessian矩阵H得到该极值点的主曲率,并通过剔除主曲率比值大于设定的阈值T2的极值点来剔除所述不稳定边缘点,其中阈值T2的取值范围为[10,15],所述主曲率比值通过比较矩阵H的特征值之间的比值来确定。
2.根据权利要求1所述的一种道路完好性检测装置,其特征是,还包括:
(1)主方向确定模块,包括依次连接的连线子模块、合并子模块和处理子模块,所述连线在模块用于将关于所述特征点的梯度方向直方图中的任意相邻两个峰值连线以形成多个子线段,所述合并子模块用于将具有相近斜率且相邻的子线段在长度方向上进行合并形成一条线段,所述处理子模块用于将多条线段中的最优线段的方向作为特征点的主方向,所述最优线段的判定公式为:
L Y = L g &OverBar; max , g &OverBar; max = m a x ( g &OverBar; L n ) , g &OverBar; L n = 1 k &Sigma; k = 1 k g k , L n &Element; L &upsi; )
其中,LY表示最优线段,为平均梯度值为的线段,为所述多条线段中第n条线段的平均梯度值,gk为所述第n条线段中的第k条子线段,Lυ为所述多条线段中线段长度大于平均线段长度的线段集合;
(2)特征提取模块,其根据所述主方向来旋转特征点的邻域,并根据旋转后的邻域对所述特征点进行描述,从而生成所述特征点的描述符,与数据库进行对比完成检测。
3.根据权利要求1所述的一种道路完好性检测装置,其特征是,所述具有相近斜率的子线段为斜率差小于预设阈值T3的子线段,所述阈值T3的取值范围为(0,0.1]。
CN201610230945.1A 2016-04-14 2016-04-14 一种道路完好性检测装置 Pending CN105913437A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610230945.1A CN105913437A (zh) 2016-04-14 2016-04-14 一种道路完好性检测装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610230945.1A CN105913437A (zh) 2016-04-14 2016-04-14 一种道路完好性检测装置

Publications (1)

Publication Number Publication Date
CN105913437A true CN105913437A (zh) 2016-08-31

Family

ID=56746908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610230945.1A Pending CN105913437A (zh) 2016-04-14 2016-04-14 一种道路完好性检测装置

Country Status (1)

Country Link
CN (1) CN105913437A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615607A (zh) * 2018-11-09 2019-04-12 福建和盛高科技产业有限公司 一种基于单张图像自定义特征的噪声检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101470896A (zh) * 2007-12-24 2009-07-01 南京理工大学 基于视频分析的机动目标飞行模式预测方法
CN103020945A (zh) * 2011-09-21 2013-04-03 中国科学院电子学研究所 一种多源传感器的遥感图像配准方法
CN104978709A (zh) * 2015-06-24 2015-10-14 北京邮电大学 描述符生成方法及装置
CN105046681A (zh) * 2015-05-14 2015-11-11 江南大学 一种基于SoC的图像显著性区域检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101470896A (zh) * 2007-12-24 2009-07-01 南京理工大学 基于视频分析的机动目标飞行模式预测方法
CN103020945A (zh) * 2011-09-21 2013-04-03 中国科学院电子学研究所 一种多源传感器的遥感图像配准方法
CN105046681A (zh) * 2015-05-14 2015-11-11 江南大学 一种基于SoC的图像显著性区域检测方法
CN104978709A (zh) * 2015-06-24 2015-10-14 北京邮电大学 描述符生成方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴京辉: "视频监控目标的跟踪与识别研究", 《中国博士学位论文全文数据库 信息科技辑》 *
张建兴: "基于注意力的目标识别算法及在移动机器人的应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615607A (zh) * 2018-11-09 2019-04-12 福建和盛高科技产业有限公司 一种基于单张图像自定义特征的噪声检测方法

Similar Documents

Publication Publication Date Title
CN103605977B (zh) 一种车道线的提取方法及装置
CN107330376B (zh) 一种车道线识别方法及系统
CN105261020B (zh) 一种快速车道线检测方法
EP3036730B1 (en) Traffic light detection
CN105678285B (zh) 一种自适应的道路鸟瞰图变换方法和道路车道检测方法
CN107392139B (zh) 一种基于霍夫变换的车道线检测方法及终端设备
WO2022126377A1 (zh) 检测车道线的方法、装置、终端设备及可读存储介质
CN104899554A (zh) 一种基于单目视觉的车辆测距方法
CN104036246B (zh) 一种基于多特征融合及多态均值的车道线定位方法
Kortli et al. A novel illumination-invariant lane detection system
CN105844337A (zh) 一种智能垃圾分类装置
CN108052904A (zh) 车道线的获取方法及装置
CN103927509A (zh) 眼睛定位方法及装置
CN104112141A (zh) 一种基于公路监控设备的货车安全带系挂状态的检测方法
CN104915642A (zh) 前方车辆测距方法及装置
CN113988112A (zh) 车道线的检测方法、装置、设备及存储介质
Chen et al. A novel lane departure warning system for improving road safety
CN110733416A (zh) 一种基于逆透视变换的车道偏离预警方法
CN111652033A (zh) 基于OpenCV的车道线检测方法
CN105928099A (zh) 一种智能空气净化机
CN105844651A (zh) 一种图像分析装置
CN105844260A (zh) 一种多功能智能清洁机器人装置
CN108090425B (zh) 一种车道线检测方法、装置及终端
CN105335704A (zh) 一种基于双线性插值的车道线识别方法与装置
CN105913437A (zh) 一种道路完好性检测装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160831