The embodiment of the present invention relates to technical field of image processing, particularly relates to a kind of vehicle-mounted 3D road Real-time Reconstruction method and device.
Summary of the invention
The embodiment of the present invention provides a kind of vehicle-mounted 3D road Real-time Reconstruction method and device, so as driver provide accurately, in real time, road surface ahead information reliably.
An aspect of the embodiment of the present invention is to provide a kind of vehicle-mounted 3D road Real-time Reconstruction method, comprising:
The first depth image M of target road is obtained respectively two neighbouring sample moment
k-1with the second depth image M
k;
By described first depth image M
k-1with described second depth image M
kthe first cloud data set N is obtained respectively through perspective projection transformation
k-1with second point cloud data acquisition N
k;
According to described first cloud data set N
k-1obtain first object cloud data set N'
k-1, and according to described second point cloud data acquisition N
kobtain the second impact point cloud data acquisition N'
k;
Calculate described first object cloud data set N'
k-1in each corresponding PFH feature histogram, and described second impact point cloud data acquisition N'
kin each the 2nd corresponding PFH feature histogram;
If a described PFH feature histogram and described 2nd PFH feature histogram coupling, then corresponding with described PFH feature histogram point is the first matching characteristic point, the point corresponding with described 2nd PFH feature histogram is the second matching characteristic point, and all described first matching characteristic points are formed the first set Q'
k-1, all described second matching characteristic points are formed the second set Q'
k;
According to the first set Q' described in ICP algorithm process
k-1with described second set Q'
kobtain rotation matrix R and translation matrix T, and Q'
k=R*Q'
k-1+ T;
According to described rotation matrix R and described translation matrix T by the first cloud data set N
k-1with second point cloud data acquisition N
kbe normalized to same coordinate system and obtain reconstruction model;
Using the input data of described reconstruction model as algorithm for pattern recognition, according to described algorithm for pattern recognition, pattern-recognition is carried out to described reconstruction model, export information of road surface corresponding to described target road to make described algorithm for pattern recognition.
Another aspect of the embodiment of the present invention is to provide a kind of vehicle-mounted 3D road Real-time Reconstruction device, comprising:
Sampling module, for obtaining the first depth image M of target road respectively two neighbouring sample moment
k-1with the second depth image M
k;
Projection module, for by described first depth image M
k-1with described second depth image M
kthe first cloud data set N is obtained respectively through perspective projection transformation
k-1with second point cloud data acquisition N
k;
Denoising module, for according to described first cloud data set N
k-1obtain first object cloud data set N'
k-1, and according to described second point cloud data acquisition N
kobtain the second impact point cloud data acquisition N'
k;
Histogram calculation module, for calculating described first object cloud data set N'
k-1in each corresponding PFH feature histogram, and described second impact point cloud data acquisition N'
kin each the 2nd corresponding PFH feature histogram;
Matching module, if for a described PFH feature histogram and described 2nd PFH feature histogram coupling, then corresponding with described PFH feature histogram point is the first matching characteristic point, the point corresponding with described 2nd PFH feature histogram is the second matching characteristic point, and all described first matching characteristic points are formed the first set Q'
k-1, all described second matching characteristic points are formed the second set Q'
k;
ICP processing module, for the first set Q' described in foundation ICP algorithm process
k-1with described second set Q'
kobtain rotation matrix R and translation matrix T, and Q'
k=R*Q'
k-1+ T;
Normalizing module, for according to described rotation matrix R and described translation matrix T by the first cloud data set N
k-1with second point cloud data acquisition N
kbe normalized to same coordinate system and obtain reconstruction model; Using the input data of described reconstruction model as algorithm for pattern recognition, according to described algorithm for pattern recognition, pattern-recognition is carried out to described reconstruction model, export information of road surface corresponding to described target road to make described algorithm for pattern recognition.
The vehicle-mounted 3D road Real-time Reconstruction method that the embodiment of the present invention provides and device, by two-dimensional depth image being converted to three-dimensional cloud data set, each corresponding PFH feature histogram in acquisition point cloud data acquisition, the matching characteristic point pair of adjacent two moment point cloud data acquisitions is found by PFH feature histogram coupling, to matching characteristic point to enforcement ICP algorithm, ask for the rotation translation matrix of adjacent two moment point cloud data acquisitions, according to rotating translation matrix, two cloud data set are normalized to the same coordinate system, achieve the fusion of cloud data set, obtain reconstruction model, road information is avoided to be subject to color, texture, the impact of the factors such as intensity of illumination, improve the computational accuracy of road information, for driver provides accurate, in real time, reliable road surface ahead information.
Embodiment
The vehicle-mounted 3D road Real-time Reconstruction method flow diagram that Fig. 1 provides for the embodiment of the present invention; The schematic diagram of the neighborhood point of any point in the impact point cloud data acquisition that Fig. 2 provides for the embodiment of the present invention.The embodiment of the present invention is for technology such as traditional road Identification based on two dimensional image, lane detection, and the reason of sensor imaging mechanism, road information is made to be easy to be subject to the impact of the factors such as color, texture, intensity of illumination, reduce the computational accuracy of road information, can not provide accurately for driver, in real time, road surface ahead information reliably, provide vehicle-mounted 3D road Real-time Reconstruction method, the concrete steps of the method are as follows:
Step S101, obtain the first depth image M of target road respectively two neighbouring sample moment
k-1with the second depth image M
k;
The embodiment of the present invention obtains the depth image of road ahead scene by the TOF imaging device being arranged on front part of vehicle, and TOF imaging device is specifically as follows 3D imaging sensor, 3D imaging sensor is with the depth image of higher employing frequency acquisition road ahead scene, such as, in the K-1 moment, 3D imaging sensor carries out single pass to road ahead scene and obtains the first depth image M
k-1, in the K moment, 3D imaging sensor carries out single pass to road ahead scene and obtains the second depth image M
k.
Step S102, by described first depth image M
k-1with described second depth image M
kthe first cloud data set N is obtained respectively through perspective projection transformation
k-1with second point cloud data acquisition N
k;
Described by described first depth image M
k-1with described second depth image M
kthe first cloud data set N is obtained respectively through perspective projection transformation
k-1with second point cloud data acquisition N
k, comprising: obtain the transformation matrix F of depth image to cloud data set; According to described transformation matrix F and described first depth image M
k-1calculate described first cloud data set N
k-1=F*M
k-1, according to described transformation matrix F and described second depth image M
kcalculate described second point cloud data acquisition N
k=F*M
k.
The parameter of demarcating according to 3D imaging sensor inside and outside location parameter obtain the transformation matrix F of depth image to cloud data set, calculate described first cloud data set N
k-1=F*M
k-1with second point cloud data acquisition N
k=F*M
k.
Step S103, according to described first cloud data set N
k-1obtain first object cloud data set N'
k-1, and according to described second point cloud data acquisition N
kobtain the second impact point cloud data acquisition N'
k;
To described first cloud data set N
k-1carry out noise reduction or simplify pre-service obtaining described first object cloud data set N'
k-1; To described second point cloud data acquisition N
kcarry out noise reduction or simplify pre-service obtaining described second impact point cloud data acquisition N'
k.
Step S104, calculate described first object cloud data set N'
k-1in each corresponding PFH feature histogram, and described second impact point cloud data acquisition N'
kin each the 2nd corresponding PFH feature histogram;
As shown in Figure 2, putting P is first object cloud data set N'
k-1in any point, apart from having H neighborhood point in this P preset range, P
tbe any one point in H the neighborhood point adjacent with a P, what be connected with a P is other neighborhood point, according to the normal n that some P is corresponding
sdeviation between normal corresponding with each the neighborhood point in H neighborhood point respectively obtains a PFH feature histogram corresponding to some P, in like manner, and first object cloud data set N'
k-1in the corresponding PFH feature histogram of each point.Described second impact point cloud data acquisition N' can be calculated according to same method
kin each the 2nd corresponding PFH feature histogram.
If the described PFH feature histogram of step S105 and described 2nd PFH feature histogram coupling, then corresponding with described PFH feature histogram point is the first matching characteristic point, the point corresponding with described 2nd PFH feature histogram is the second matching characteristic point, and all described first matching characteristic points are formed the first set Q'
k-1, all described second matching characteristic points are formed the second set Q'
k;
First object cloud data set N' can be obtained by step S104
k-1in a PFH feature histogram corresponding to each point, and the second impact point cloud data acquisition N'
kin each the 2nd corresponding PFH feature histogram, if first object cloud data set N'
k-1in a PFH feature histogram of a some correspondence and the second impact point cloud data acquisition N'
kin the 2nd PFH feature histogram of a some correspondence identical, then these two points are respectively matching characteristic point, first object cloud data set N'
k-1in matching characteristic point form first set Q'
k-1, the second impact point cloud data acquisition N'
kin matching characteristic point form second set Q'
k.
First set Q' described in step S106, foundation ICP algorithm process
k-1with described second set Q'
kobtain rotation matrix R and translation matrix T, and Q'
k=R*Q'
k-1+ T;
To described first set Q'
k-1with described second set Q'
kcarry out ICP algorithm process, obtain rotation matrix R and translation matrix T, and Q'
k=R*Q'
k-1+ T.
Step S107, according to described rotation matrix R and described translation matrix T by the first cloud data set N
k-1with second point cloud data acquisition N
kbe normalized to same coordinate system and obtain reconstruction model.
Described according to described rotation matrix R and described translation matrix T by the first cloud data set N
k-1with second point cloud data acquisition N
kbe normalized to same coordinate system and obtain reconstruction model, comprising:
According to described rotation matrix R, described translation matrix T and described first cloud data set N
k-1calculate thirdly cloud data acquisition N "
k=R*N
k-1+ T;
According to described thirdly cloud data acquisition N "
kwith described second point cloud data acquisition N
kobtain reconstruction model N "
k+ N
k.
Step S108, using the input data of described reconstruction model as algorithm for pattern recognition, according to described algorithm for pattern recognition, pattern-recognition is carried out to described reconstruction model, export information of road surface corresponding to described target road to make described algorithm for pattern recognition.
The embodiment of the present invention is by being converted to three-dimensional cloud data set by two-dimensional depth image, each corresponding PFH feature histogram in acquisition point cloud data acquisition, the matching characteristic point pair of adjacent two moment point cloud data acquisitions is found by PFH feature histogram coupling, to matching characteristic point to enforcement ICP algorithm, ask for the rotation translation matrix of adjacent two moment point cloud data acquisitions, according to rotating translation matrix, two cloud data set are normalized to the same coordinate system, achieve the fusion of cloud data set, obtain reconstruction model, road information is avoided to be subject to color, texture, the impact of the factors such as intensity of illumination, improve the computational accuracy of road information, for driver provides accurate, in real time, reliable road surface ahead information.
The coordinate schematic diagram of any point in the impact point cloud data acquisition that Fig. 3 provides for the embodiment of the present invention.The schematic diagram of the PFH feature histogram that Fig. 4 provides for the embodiment of the present invention.The structural drawing that principal and subordinate's computing method that Fig. 5 provides for the embodiment of the present invention are suitable for.On the basis of above-described embodiment, described calculating described first object cloud data set N'
k-1in each corresponding PFH feature histogram, comprising:
Define described first object cloud data set N'
k-1the coordinate of middle any point P is (u, v, w), wherein, and u=n
s,
w=u × v, n
sthe normal that a P is corresponding, P
tany one point in H the neighborhood point adjacent with a P, and the normal n that some P is corresponding
swith a P
tcorresponding normal n
tbetween deviation be expressed as
wherein, α=v*n
t,
θ=arctan (w*n
t, u*n
t); The horizontal ordinate dividing two-dimensional coordinate at equal intervals obtains multiple calibration points, and described multiple calibration points is corresponding a different set of respectively
value; Add up for each calibration points in described two-dimensional coordinate the number meeting pre-conditioned impact point in described H neighborhood point and obtain a described PFH feature histogram, it is corresponding that the deviation between the normal that the described pre-conditioned normal corresponding for described impact point is corresponding with some P equals described calibration points
value.
For first object cloud data set N'
k-1in some P, define its coordinate for (u, v, w), wherein, u=n
s,
w=u × v, n
sthe normal that a P is corresponding, P
tany one point in H the neighborhood point adjacent with a P, and the normal n that some P is corresponding
swith a P
tcorresponding normal n
tbetween deviation be expressed as
wherein, α=v*n
t,
θ=arctan (w*n
t, u*n
t).So, the normal n that P is corresponding
sa deviation is had between the normal that any one point in H the neighborhood point be adjacent is corresponding
in like manner, first object cloud data set N'
k-1in H neighborhood point being adjacent of normal corresponding to any one point in any one some correspondence normal between have a deviation
As shown in Figure 4, the horizontal ordinate dividing two-dimensional coordinate at equal intervals obtains multiple calibration points, and preferably, the embodiment of the present invention is divided into 125 calibration points, is specially 0-124, and each calibration points correspondence is a different set of
value; For each calibration points of 125 calibration points, the number of pre-conditioned impact point is met in H neighborhood point in the preset range of statistics point P, this number is as ordinate corresponding to this calibration points, and it is corresponding that the deviation between the concrete pre-conditioned normal corresponding with some P for the normal of impact point equals described calibration points
value, so obtains the PFH feature histogram that some P is corresponding, i.e. described first object cloud data set N'
k-1in each point can portray with a PFH feature histogram, in like manner, described second impact point cloud data acquisition N' can be obtained
kin each the 2nd corresponding PFH feature histogram.
On the basis of the embodiment of the present invention, it is the ability that whole system provides high-speed parallel computing power to calculate with reply complex logic.Owing to needing cloud data amount to be processed huge, as shown in Figure 5, the account form that native system adopts principal and subordinate to combine implements related algorithm, and main computation unit is taken on by universal cpu, for realizing complicated algorithm logic; Taken on by GPU or special parallel computation unit from computing unit, be mainly used in realizing a large amount of double counting; Main computation unit need converge conjunction and is divided into some subsets by point to be processed, is distributed to the sub-computing unit that each is parallel, after being disposed, by main computation unit aggregation process result.
The embodiment of the present invention specifically provides the method how calculating each corresponding PFH feature histogram in impact point cloud data acquisition, in addition, is improve the speed of data processing by principal and subordinate's computing method.
The structural drawing of the vehicle-mounted 3D road Real-time Reconstruction device that Fig. 6 provides for the embodiment of the present invention.The vehicle-mounted 3D road Real-time Reconstruction device that the embodiment of the present invention provides can perform the treatment scheme that vehicle-mounted 3D road Real-time Reconstruction embodiment of the method provides, as shown in Figure 6, vehicle-mounted 3D road Real-time Reconstruction device 50 comprises sampling module 51, projection module 52, denoising module 53, histogram calculation module 54, matching module 55, ICP processing module 56 and normalizing module 57, wherein, sampling module 51 is for obtaining the first depth image M of target road respectively two neighbouring sample moment
k-1with the second depth image M
k; Projection module 52 is for by described first depth image M
k-1with described second depth image M
kthe first cloud data set N is obtained respectively through perspective projection transformation
k-1with second point cloud data acquisition N
k; Denoising module 53 is for according to described first cloud data set N
k-1obtain first object cloud data set N'
k-1, and according to described second point cloud data acquisition N
kobtain the second impact point cloud data acquisition N'
k; Histogram calculation module 54 is for calculating described first object cloud data set N'
k-1in each corresponding PFH feature histogram, and described second impact point cloud data acquisition N'
kin each the 2nd corresponding PFH feature histogram; If matching module 55 is for a described PFH feature histogram and described 2nd PFH feature histogram coupling, then corresponding with described PFH feature histogram point is the first matching characteristic point, the point corresponding with described 2nd PFH feature histogram is the second matching characteristic point, and all described first matching characteristic points are formed the first set Q'
k-1, all described second matching characteristic points are formed the second set Q'
k; ICP processing module 56 is for the first set Q' described in foundation ICP algorithm process
k-1with described second set Q'
kobtain rotation matrix R and translation matrix T, and Q'
k=R*Q'
k-1+ T; Normalizing module 57 for according to described rotation matrix R and described translation matrix T by the first cloud data set N
k-1with second point cloud data acquisition N
kbe normalized to same coordinate system and obtain reconstruction model; Using the input data of described reconstruction model as algorithm for pattern recognition, according to described algorithm for pattern recognition, pattern-recognition is carried out to described reconstruction model, export information of road surface corresponding to described target road to make described algorithm for pattern recognition.
The embodiment of the present invention is by being converted to three-dimensional cloud data set by two-dimensional depth image, each corresponding PFH feature histogram in acquisition point cloud data acquisition, the matching characteristic point pair of adjacent two moment point cloud data acquisitions is found by PFH feature histogram coupling, to matching characteristic point to enforcement ICP algorithm, ask for the rotation translation matrix of adjacent two moment point cloud data acquisitions, according to rotating translation matrix, two cloud data set are normalized to the same coordinate system, achieve the fusion of cloud data set, obtain reconstruction model, road information is avoided to be subject to color, texture, the impact of the factors such as intensity of illumination, improve the computational accuracy of road information, for driver provides accurate, in real time, reliable road surface ahead information.
On the basis of above-described embodiment, projection module 52 is specifically for obtaining the transformation matrix F of depth image to cloud data set; According to described transformation matrix F and described first depth image M
k-1calculate described first cloud data set N
k-1=F*M
k-1, according to described transformation matrix F and described second depth image M
kcalculate described second point cloud data acquisition N
k=F*M
k.
Denoising module 53 is specifically for described first cloud data set N
k-1carry out noise reduction or simplify pre-service obtaining described first object cloud data set N'
k-1; To described second point cloud data acquisition N
kcarry out noise reduction or simplify pre-service obtaining described second impact point cloud data acquisition N'
k.
Histogram calculation module 54 is specifically for defining described first object cloud data set N'
k-1the coordinate of middle any point P is (u, v, w), wherein, and u=n
s,
w=u × v, n
sthe normal that a P is corresponding, P
tany one point in H the neighborhood point adjacent with a P, and the normal n that some P is corresponding
swith a P
tcorresponding normal n
tbetween deviation be expressed as
wherein, α=v*n
t,
θ=arctan (w*n
t, u*n
t); The horizontal ordinate dividing two-dimensional coordinate at equal intervals obtains multiple calibration points, and described multiple calibration points is corresponding a different set of respectively
value; Add up for each calibration points in described two-dimensional coordinate the number meeting pre-conditioned impact point in described H neighborhood point and obtain a described PFH feature histogram, it is corresponding that the deviation between the normal that the described pre-conditioned normal corresponding for described impact point is corresponding with some P equals described calibration points
value.
Normalizing module 57 is specifically for according to described rotation matrix R, described translation matrix T and described first cloud data set N
k-1calculate thirdly cloud data acquisition N "
k=R*N
k-1+ T; According to described thirdly cloud data acquisition N "
kwith described second point cloud data acquisition N
kobtain reconstruction model N "
k+ N
k.
The vehicle-mounted 3D road Real-time Reconstruction device that the embodiment of the present invention provides can specifically for performing the embodiment of the method that above-mentioned Fig. 1 provides, and concrete function repeats no more herein.
The embodiment of the present invention specifically provides the method how calculating each corresponding PFH feature histogram in impact point cloud data acquisition, in addition, is improve the speed of data processing by principal and subordinate's computing method.
In sum, the embodiment of the present invention is by being converted to three-dimensional cloud data set by two-dimensional depth image, each corresponding PFH feature histogram in acquisition point cloud data acquisition, the matching characteristic point pair of adjacent two moment point cloud data acquisitions is found by PFH feature histogram coupling, to matching characteristic point to enforcement ICP algorithm, ask for the rotation translation matrix of adjacent two moment point cloud data acquisitions, according to rotating translation matrix, two cloud data set are normalized to the same coordinate system, achieve the fusion of cloud data set, obtain reconstruction model, road information is avoided to be subject to color, texture, the impact of the factors such as intensity of illumination, improve the computational accuracy of road information, for driver provides accurate, in real time, reliable road surface ahead information, specifically provide the method how calculating each corresponding PFH feature histogram in impact point cloud data acquisition, in addition, improve the speed of data processing by principal and subordinate's computing method.
In several embodiment provided by the present invention, should be understood that, disclosed apparatus and method, can realize by another way.Such as, device embodiment described above is only schematic, such as, the division of described unit, be only a kind of logic function to divide, actual can have other dividing mode when realizing, such as multiple unit or assembly can in conjunction with or another system can be integrated into, or some features can be ignored, or do not perform.Another point, shown or discussed coupling each other or direct-coupling or communication connection can be by some interfaces, and the indirect coupling of device or unit or communication connection can be electrical, machinery or other form.
The described unit illustrated as separating component or can may not be and physically separates, and the parts as unit display can be or may not be physical location, namely can be positioned at a place, or also can be distributed in multiple network element.Some or all of unit wherein can be selected according to the actual needs to realize the object of the present embodiment scheme.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, also can be that the independent physics of unit exists, also can two or more unit in a unit integrated.Above-mentioned integrated unit both can adopt the form of hardware to realize, and the form that hardware also can be adopted to add SFU software functional unit realizes.
The above-mentioned integrated unit realized with the form of SFU software functional unit, can be stored in a computer read/write memory medium.Above-mentioned SFU software functional unit is stored in a storage medium, comprising some instructions in order to make a computer equipment (can be personal computer, server, or the network equipment etc.) or processor (processor) perform the part steps of method described in each embodiment of the present invention.And aforesaid storage medium comprises: USB flash disk, portable hard drive, ROM (read-only memory) (Read-OnlyMemory, ROM), random access memory (RandomAccessMemory, RAM), magnetic disc or CD etc. various can be program code stored medium.
This neighborhood technique personnel can be well understood to, for convenience and simplicity of description, only be illustrated with the division of above-mentioned each functional module, in practical application, can distribute as required and by above-mentioned functions and be completed by different functional modules, inner structure by device is divided into different functional modules, to complete all or part of function described above.The specific works process of the device of foregoing description, with reference to the corresponding process in preceding method embodiment, can not repeat them here.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, the those of ordinary skill of this neighborhood is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.