CN106780509A - Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic - Google Patents
Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic Download PDFInfo
- Publication number
- CN106780509A CN106780509A CN201611088463.3A CN201611088463A CN106780509A CN 106780509 A CN106780509 A CN 106780509A CN 201611088463 A CN201611088463 A CN 201611088463A CN 106780509 A CN106780509 A CN 106780509A
- Authority
- CN
- China
- Prior art keywords
- cloud data
- building
- segmentation
- cloud
- cluster
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/231—Hierarchical techniques, i.e. dividing or merging pattern sets so as to obtain a dendrogram
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention discloses a kind of building object point cloud layer for merging multidimensional characteristic time cluster segmentation method, initial segmentation is carried out to mutual discontinuous cloud data in building cloud data first, then G K clustering algorithms are utilized, ground floor subdivision is carried out to it with reference to the spectral signature of cloud data, the normal direction measure feature and curvature feature for obtaining cloud data carry out second layer subdivision, are repeated twice segmentation and are required until meeting.Present invention building object point cloud layer time cluster segmentation method takes full advantage of the density information of Point Cloud Data from Three Dimension Laser Scanning, can be on the premise of no priori, multiple cloud data block apart from each other and than comparatively dense is carried out into initial segmentation, simultaneously, with reference to the spectral signature and geometric properties of cloud data, to being finely divided through the cloud data block after initial segmentation, until each cloud data block has single geometric properties, untill being modeled using simple Mathematical Modeling, the method can not only extract building cloud data from the environment of surrounding, and building object point cloud can be resolved into different planes, for the reconstruction of building is had laid a good foundation.
Description
Technical field
The present invention relates to a kind of cloud data partitioning algorithm, and in particular to one kind building object point cloud layer time cluster segmentation side
Method.
Background technology
Technology is swept using three-dimensional laser and measurement is scanned to building, can obtain building and retouch cloud data, in order to
Building body characteristicses, i.e. building or its composition structure, the shape at position, position, size, length are extracted based on cloud data
The description of the key geometric sense such as degree, area, volume, it is necessary to carry out the weight of BUILDINGS MODELS to it using building three dimensional point cloud
Build process.Whether correctly the correctness that model is set up directly influences the measurement of building body characteristicses, therefore reasonable utilization is swept
It is extremely important that the cloud data retouched sets up building point cloud model.But building point cloud model is often special containing multiple surface geometries
(such as plane, the face of cylinder, circular conical surface, sphere, free form surface) is levied, if directly carrying out model weight to it using cloud data
Structure, then can cause the difficulty of mathematical notation and the fitting algorithm treatment of surface model to increase, or even cannot use suitable mathematical table
The surface model of this complexity is described up to formula.In order to meet the requirement of Surface Reconstruction, it is necessary to be carried out to three dimensional point cloud
Region segmentation, purpose seeks to the region with similar geometry characteristic in three dimensional point cloud to be split, combines, with after an action of the bowels
Continuous three dimensional point cloud Surface Reconstruction.Therefore, for the architectural characteristic of building cloud data, research is adapted to building
The method of cloud data segmentation is significant.
At present, cloud data partitioning algorithm mainly includes that the partitioning algorithm based on border, the segmentation based on region growing are calculated
Method, cluster segmentation algorithm etc..But existing various partitioning algorithms, the generally a certain feature using cloud data or comprehensive utilization
Multiple features of cloud data, it is intended to which cloud data is disposably separated, this comes for complicated building cloud data
The situation for easily causing segmentation deficiency or over-segmentation is said, the correct segmentation of complex building cloud data is difficult to realize.
The content of the invention
Goal of the invention:The purpose of the present invention is to solve the shortcomings of the prior art, there is provided it is a kind of according to cloud data
Spectral signature and geometric properties carry out the building object point cloud layer time cluster segmentation method of the fusion multidimensional characteristic of accurate segmentation.
Technical scheme:The invention provides a kind of building object point cloud layer for merging multidimensional characteristic time cluster segmentation method, bag
Include following steps:
(1) cloud initial segmentation is put:Using DBSCAN (Density-Based Spatial Clustering of
Applications with Noise) density partitioning algorithm enters to mutual discontinuous cloud data in building cloud data
Row initial segmentation, the cloud data block being separated from each other;This be in order to apart from each other while again than the point cloud number of comparatively dense
According to being clustered, the cloud data of initial segmentation has thus been obtained;
(2) spectral signature segmentation:By the continuous cloud data after initial segmentation, its spectral signature often shows difference
Property, such as fabric structure is different with the color of the environment such as trees, meadow of surrounding, reflected intensity, the same beam of building masonry wall, post
Color, the reflected intensity of son etc. be not also equal.Using G-K (Gustafson-Kessel) clustering algorithm, with reference to cloud data
Spectral signature is further split to it, and building is finely divided between surrounding environment or building different structure, obtains
To the cloud data based on the subdivision of spectral signature multilayer:
(3) geometric properties segmentation:Its normal direction measure feature and curvature feature are calculated according to cloud data, and using G-K clusters
Algorithm, to splitting through spectral signature after cloud data carry out further geometric properties subdivision;
(4) data inspection:To carrying out geometry inspection by the cloud data after subdivision, the cloud data energy after segmentation is judged
It is no to meet the requirement being modeled using simple Mathematical Modeling, if meet requiring, preserve the cloud data after segmentation;It is no
Then, (2nd) and (3) step are repeated, continuation segmentation is carried out to cloud data.
Further, step (2) described spectral signature is extracted using scanner color characteristic and reflected intensity feature.
Further, the extraction of cloud data color characteristic is that the colorful CCD camera carried using scanner shoots tested
The panorama photochrome of object, and the colouring information of testee is obtained, with reference to textures technology, by acquired testee
Color, texture are added in surveyed cloud data, obtain the three-dimensional RGB information of measured object.
Further, the color space of the panorama photochrome for being obtained due to scanning is RGB, is to use for RGB color pattern
Three kinds of primary colours of red, green, blue represent shades of colour, but RGB color can not be combined with color space perceptually well
Get up, therefore need the conversion of RGB to HSV herein:
V=max
Wherein, (R, G, B) is respectively a red, green and blue coordinate for color, and their value is the real number between 0 to 1;
Max is the maximum in R, G and B, and min is the minimum value in R, G and B value;(H, S, V) represents the tone of color, saturation respectively
Degree, lightness.
In general, laser reflection energy is referred to as laser reflection intensity with the ratio of Laser emission energy.Work as laser
When being transmitted into measured target surface with certain wavelength, can be because the degree of roughness on measured target surface be scattered, by a part
Laser light scattering to other directions without being reflected back toward laser, while also can be because of the characteristic on measured target surface (physically or chemically
Characteristic) absorb the factors such as laser energy and cause the laser energy being reflected back in laser to be less than the energy of Laser emission.Cause
This, it is possible to use the equipment that scanner is carried obtains the data reflected intensity feature of point cloud.
Further, G-K clustering algorithms are comprised the following steps:
If being clustered cloud data collection is combined into X={ x1,x2,…,xn, each of which data xkThere is d feature to refer to
Mark, thus its Characteristic index matrix is:
Data set X is divided into c classes (2≤c≤n), if c cluster centre vector is:
If μjk∈ [0,1] represents degree of membership of k-th data for jth class, and meetsThen
Fuzzy partition matrix is:
G-K algorithms clustering criteria obtains minimum value to make following object function:
Wherein:mjRepresent cluster centre, b>1 is Weighted Index, and the overlap between the more big each clusters of b is more;Similarity degree
Flow function is
The distance of k-th data and the cluster centre of jth class is represented, it determines the shape of cluster;Wherein AjIt is one
Individual positive definite matrix, by the cluster covariance matrix F for approximately reflecting each cluster true formjDetermine, work as AjDuring for unit matrix, degree
Flow function uses Euclidean distance;Wherein
ρjIt is a constant for each cluster, in the case where priori is lacked, takes ρj=1 causes each cluster
Capacity it is roughly the same;
Object function Jf(U, M) is minimized can be expressed as about fasciculation problem:
Solved with lagrange Multiplier Methods:
And work asWhen, ujk=1, ulk=0 (l ≠ j)
Iterated to calculate more than, determine corresponding to cloud data cluster number and the cluster of each class cloud data in
The heart.
Further, step (3) carries out the calculating of normal vector using the tangent plane of cloud data, first by partial points cloud number
According to one tangent plane of fitting, arbitrfary point piNormal vector determined by the normal vector estimation of the plane:
Plane fitting:According to cloud data pi(i=0,1 ..., n) in the quadratic sum for arriving the plan range a little most
It is small, determine plane P (ui,vi):
Wherein, n represents the number at cloud data midpoint;
The general type of the plane expression formula is:
Ax+by+cz+d=0
Least square plane fitting object function be:
Ax=0
Wherein:
Using Jacobi method calculating matrix ATThe eigenvalue λ of AiWith corresponding characteristic vector xi(i=1 ..., 4), then definitely
It is worth minimum eigenvalue λiCorresponding characteristic vector xiIt is plane parameter a to be asked, the least square solution of b, c, d;
2. plane normal vector determines:According to the general expression of plane equation, plane normal vector is represented by
To avoid the problem of plane parameter a, b, c, d dependent, unitization treatment is carried out to required plane normal vector, such as following formula institute
Show:
The per unit system arrow of arbitrfary point in cloud data at local surface is:Realize cloud data
The estimation of normal vector.
Further, in the curvature characteristic root strong point cloud of step (3) cloud data each point average curvature and Gaussian curvature meter
Calculating needs to carry out quadric fitting to the vertex neighborhood, after determining the principal curvatures and its principal direction of curved surface S, it is possible to calculate
The curvature characteristic of the data point;
A, Quadratic Surface Fitting:By cloud data pi(i=1,2 ... k) carry out Quadratic Surface Fitting in local neighborhood, intend
Close equation general type be:
S (u, v)=au2+buv+cv2+du+ev
The object function of fitting is:
In formula, u and v is Surface Parameters, and the least square solution of fitting surface can be obtained using singular value decomposition method;
B, Curvature Estimate:According to calculate quadratic surface S parametric equation, using the principal curvatures and principal direction of surface points as
Point piPrincipal curvatures and principal direction, its parameter is calculated as follows:
Su|(0,0)=(1,0,2au+bv+d) |(0,0)=(1,0, d) and Suu|=(0,0,2a)
Sv|(0,0)=(0,1, bu+2cv+e) |(0,0)=(0,1, e) and Svv|(0,0)=(0,0,2c)
Suv|(0,0)=(0,0, b)
Wherein, SuIt is curved surface S to the first derivative of parameter u, SuuIt is second dervative;SvFor curved surface S leads to the single order of parameter v
Number, SvvIt is second dervative;SuvIt is curved surface S to parameter u, the second dervative of v;N is sweared for the per unit system of curved surface;
According to above parameter, can calculate:
E=Su·Su=1+d2And F=Su·Sv=de
G=Sv·Sv=1+e2And
And
Wherein, E, F, G are the first fundamental quantity of curved surface, and L, M, N are the second fundamental quantity of curved surface;
Then the Gaussian curvature and average curvature values at P points are:
Minimum principal curvatures θminWith maximum principal curvatures θmaxComputing formula is respectively:
Beneficial effect:Present invention building object point cloud layer time cluster segmentation method takes full advantage of 3 D laser scanning point cloud number
According to density information, can be on the premise of no priori, by multiple cloud data blocks apart from each other and than comparatively dense
Initial segmentation is carried out, meanwhile, with reference to the spectral signature and geometric properties of cloud data, to through the cloud data block after initial segmentation
It is finely divided, until each cloud data block has single geometric properties, can be built using simple Mathematical Modeling
Untill mould, the method can not only extract building cloud data from the environment of surrounding, and can be by building object point cloud point
Solution, into different planes, is that the reconstruction of building is had laid a good foundation.
Brief description of the drawings
Fig. 1 is the original scan cloud data of building;
Fig. 2 is present invention building object point cloud layer time cluster segmentation flow chart;
Fig. 3 is the cloud data of initial segmentation;
Fig. 4 is segmented for the second layer of cloud data block;
Fig. 5 is segmented for the third layer of cloud data block;
Fig. 6 is split for the geometric properties of data block I;
Fig. 7 is split for the geometric properties of data set II-3- (4);
Fig. 8 is split for the geometric properties of data set II-3- (3).
Specific embodiment
Technical solution of the present invention is described in detail below, but protection scope of the present invention is not limited to the implementation
Example.
Embodiment:The present embodiment uses FARO companies of U.S. FOCUS3DThree-dimensional laser scanner is to building and its surrounding ring
Border is scanned measurement, and the cloud data of acquisition is as shown in Figure 1.First, by building original image understand the building and its
Surrounding environment is divided into 1. top of building, 2. building masonry wall, 3. upstairs ground, 4. passageway barricade, 5. stair, 6. from top to bottom
Seven parts constitute altogether for trees, 7. building outside ground etc., if to be modeled to building, by the same surrounding of building
Environment separate, then building Each part is finely divided, can just access the cloud data of Direct Modeling.Cause
This, is split using the building object point cloud layer time cluster segmentation method for merging multidimensional characteristic to the cloud data, including following
Step, as shown in Figure 2:
1st, initial segmentation
Any one neighborhood of a point K is 20 for 1.97, in space to set the radius Eps in DBSCAN density partitioning algorithms, with
The locus of cloud data is split as input object to original point cloud data, from as shown in Figure 3, the algorithm
Will be apart from each other in original point cloud data while the cloud data such as roof, trees again than comparatively dense is separated, used
Cloud data block I that green (top of building), red (building body and trees) and blueness (isolated trees) are represented respectively,
II, III, due to that cannot be represented with color in figure, therefore are pointed out with arrow.
2nd, the segmentation based on spectral signature
By after initial segmentation, trees, the ground of the main part of building still with periphery etc. mixes, and builds
The structure and its surrounding enviroment of thing are increasingly complex, therefore, segmentation can not use geometric properties.But between building and surrounding environment
Have between spectral signature than larger difference.
(1) ground floor subdivision
Color characteristic (R, G, B) is assigned to cloud data first, and color mode is converted into by (R, G, B) (H, S,
V), characteristic vector (X, Y, Z, H) is collectively constituted by H in locus feature (X, Y, Z) and color characteristic, using G-K cluster sides
Method is further segmented to cloud data block II:
1. it is that 7, Weighted Index b is that 2, stopping criterion for iteration ε is 10 to set initialization clusters number c-6, initial fuzzy division
Matrix U is 210242 rows and 7 column matrix;
2. initial new cluster centre M is calculated0It is the column matrix of 7 row 4:
3. initially fuzzy covariance matrix F and positive definite matrix A is calculated;
4. according to measuring similarity functionFuzzy partition matrix U is updated;
5. calculated through 202 loop iterations, condition ‖ U (I+1)-U (I) ‖ < 1 × 10-6Set up, calculate final cluster
Center M is:
And each puts affiliated cluster number in obtaining cloud data.
As shown in figure 4, cluster identical each point in cloud data is classified as into a class, cloud data block II is divided into
7 parts, but only II-1 and II-2 are relatively complete, i.e., and outside ground and trees can intactly separate from data II;Its
There is not enough and over-segmentation the situation of segmentation in remaining part point cloud data, will split incomplete cloud data and be classified as cloud data block
II-3, in addition it is also necessary to further subdivision.
(2) second layer subdivision
The spectral signature of building main body point cloud II-3 is calculated also according to the above method, thus spectral signature and space are special
Composition characteristic phasor (X, Y, Z, H) is levied, building main body point cloud II-3 is finely divided using G-K clustering algorithms.
As shown in figure 5, building main body has been divided into 6 parts:II-3- (1) is facade wall, II-3- (2)
Upstairs ground, II-3- (3) is stair, and II-3- (4) is upstairs passageway wall, and II-3- (5) is building east side wall, II-3-
(6) it is building west side wall.Wherein facade wall and upstairs ground has single geometric properties.
(3) segmented based on geometric properties point cloud
By after initial segmentation and spectral signature segmentation, cloud data block is for example:I (building roof), II-3- (3) (building
Ladder), and II-3- (4) (passageway wall), still without single geometric properties.Accordingly, it would be desirable to enter to advance using geometric properties
The subdivision of one step:
1. building roof
By after initial segmentation, the normal direction of building roof cloud data block I is calculated according to the method in claim
Amount and curvatureBy the space characteristics and geometric properties composition characteristic phasor of cloud data block IWith reference to
G-K clustering algorithms, are finely divided to cloud data block I.
As shown in fig. 6, building roof divide into four parts:I-1, I-2, I-3 and I-4, wherein, cloud data block I-
1st, I-2 and I-3 are the planes with single geometric properties, their normal vector be respectively (0.212, -0.313,0.926),
(- 0.226,0.314,0.922) and (0.324,0.226,0.919).
2. passageway wall segmentation upstairs
By after second layer subdivision, calculating the normal vector and curvature in passageway barricade cloud data block II-3- (4)
By the space characteristics and geometric properties composition characteristic phasor of cloud data block II-3- (4)With reference to G-K clusters
Algorithm, is finely divided to cloud data block II-3- (4).
As shown in fig. 7, passageway wall II-3- (4) has been divided into three planes with single geometric properties:II-3-
(4)-[1], II-3- (4)-[2] and II-3- (4)-[3].Their normal vector is respectively:(0.573, -0.819,0.018), (-
0.819, -0.573,0.025) and (- 0.572,0.819,0.040).
3. stair segmentation
By after second layer subdivision, calculating the normal vector and curvature of building stair cloud data block II-3- (3)By the space characteristics and geometric properties composition characteristic phasor of cloud data block II-3- (3)With reference to
G-K clustering algorithms, are finely divided to cloud data block II-3- (3).
As shown in figure 8 above, stair are divided into four parts:II-3- (3)-[1] stair right side wall, II-3- (3)-[2]
Stair left side wall, II-3- (3)-[3] stair perpendicular set, II-3- (3)-[4] stair horizontal plane set.Their method
Vector be respectively (- 0.572,0.820, -0.003), (0.572, -0.821,0.002), (- 0.002, -0.002,1.000) and
(0.820,0.572, -0.007).
Analyzed more than, the building object point cloud layer based on multidimensional characteristic proposed by the present invention time cluster segmentation method,
Not only building can be extracted from complicated surrounding environment, detailed segmentation can also be carried out to building cloud data.Finally
The building each several part cloud data with single geometric properties is obtained, is that the reconstruction and feature extraction of building lay good
Basis.
Claims (7)
1. a kind of building object point cloud layer time cluster segmentation method for merging multidimensional characteristic, it is characterised in that:Comprise the following steps:
(1) cloud initial segmentation is put:Using DBSCAN density partitioning algorithm to mutual discrete point cloud in building cloud data
Data carry out initial segmentation, the cloud data block being separated from each other;
(2) spectral signature segmentation:Using G-K clustering algorithms, it is further split with reference to the spectral signature of cloud data,
Building is finely divided between surrounding environment or building different structure, the point cloud based on the subdivision of spectral signature multilayer is obtained
Data:
(3) geometric properties segmentation:Its normal direction measure feature and curvature feature are calculated according to cloud data, and utilize G-K clustering algorithms,
Cloud data after to splitting through spectral signature carries out further geometric properties subdivision;
(4) data inspection:To carrying out geometry inspection by the cloud data after subdivision, judge that can the cloud data after segmentation expire
The requirement that foot is modeled using simple Mathematical Modeling, if meet requiring, preserves the cloud data after segmentation;Otherwise, weight
Multiple (2nd) and (3) step, continuation segmentation is carried out to cloud data.
2. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 1, it is characterised in that:
Color characteristic and reflected intensity feature that step (2) described spectral signature is extracted using scanner.
3. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 2, it is characterised in that:
The extraction of cloud data color characteristic is the colored photograph of panorama that the colorful CCD camera carried using scanner shoots testee
Piece, and the colouring information of testee is obtained, with reference to textures technology, the color of acquired testee, texture are added to
In surveyed cloud data, the three-dimensional RGB information of measured object is obtained.
4. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 3, it is characterised in that:
The color mode of the panorama photochrome that scanner is obtained is HSV patterns by RGB patten transformations:
V=max
Wherein, (R, G, B) is respectively a red, green and blue coordinate for color, and their value is the real number between 0 to 1;max
It is the maximum in R, G and B, min is the minimum value in R, G and B value;(H, S, V) represent respectively the tone of color, saturation degree,
Lightness.
5. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 1, it is characterised in that:
G-K clustering algorithms are comprised the following steps:
If being clustered cloud data collection is combined into X={ x1,x2,…,xn, each of which data xkThere is d characteristic index, thus
Its Characteristic index matrix is:
Data set X is divided into c classes (2≤c≤n), if c cluster centre vector is:
If μjk∈ [0,1] represents degree of membership of k-th data for jth class, and meetsThen obscure
Matrix dividing is:
G-K algorithms clustering criteria obtains minimum value to make following object function:
Wherein:mjRepresent cluster centre, b>1 is Weighted Index, and the overlap between the more big each clusters of b is more;Measuring similarity function
For
The distance of k-th data and the cluster centre of jth class is represented, it determines the shape of cluster;Wherein AjFor one just
Set matrix, by the cluster covariance matrix F for approximately reflecting each cluster true formjDetermine, work as AjDuring for unit matrix, letter is measured
Number uses Euclidean distance;Wherein
ρjIt is a constant for each cluster, in the case where priori is lacked, takes ρj=1 appearance for causing each cluster
Amount is roughly the same;
Object function Jf(U, M) is minimized can be expressed as about fasciculation problem:
Solved with lagrange Multiplier Methods:
And work asWhen, ujk=1, ulk=0 (l ≠ j)
Iterated to calculate more than, determine the cluster centre of the cluster number and each class cloud data corresponding to cloud data.
6. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 1, it is characterised in that:
Step (3) carries out the calculating of normal vector using the tangent plane of cloud data, is cut by local Points cloud Fitting one is micro- first
Plane, arbitrfary point piNormal vector determined by the normal vector estimation of the plane:
1. plane fitting:According to making cloud data pi(i=0,1 ..., n) in the quadratic sum for arriving the plan range a little most
It is small, determine plane P (ui,vi):
Wherein, n represents the number at cloud data midpoint;
The general type of the plane expression formula is:
Ax+by+cz+d=0
Least square plane fitting object function be:
Ax=0
Wherein:
Using Jacobi method calculating matrix ATThe eigenvalue λ of AiWith corresponding characteristic vector xi(i=1 ..., 4), then absolute value is most
Small eigenvalue λiCorresponding characteristic vector xiIt is plane parameter a to be asked, the least square solution of b, c, d;
2. plane normal vector determines:According to the general expression of plane equation, plane normal vector is represented byTo keep away
Exempt from the problem of plane parameter a, b, c, d dependent, unitization treatment is carried out to required plane normal vector, be shown below:
The per unit system arrow of arbitrfary point in cloud data at local surface is:Realize cloud data normal direction
The estimation of amount.
7. the building object point cloud layer time cluster segmentation method of fusion multidimensional characteristic according to claim 1, it is characterised in that:
The average curvature of each point and Gaussian curvature are calculated and needed to the vertex neighborhood in the curvature characteristic root strong point cloud of step (3) cloud data
Quadric fitting is carried out, after determining the principal curvatures and its principal direction of curved surface S, it is possible to which the curvature for calculating the data point is special
Property;
A, Quadratic Surface Fitting:By cloud data pi(i=1,2 ... k) carry out Quadratic Surface Fitting, fitting side in local neighborhood
The general type of journey is:
S (u, v)=au2+buv+cv2+du+ev
The object function of fitting is:
In formula, u and v is Surface Parameters, and the least square solution of fitting surface can be obtained using singular value decomposition method;
B, Curvature Estimate:According to the parametric equation of the quadratic surface S for calculating, using the principal curvatures and principal direction of surface points as point pi
Principal curvatures and principal direction, its parameter is calculated as follows:
Su|(0,0)=(1,0,2au+bv+d) |(0,0)=(1,0, d) and Suu|=(0,0,2a)
Sv|(0,0)=(0,1, bu+2cv+e) |(0,0)=(0,1, e) and Svv|(0,0)=(0,0,2c)
Suv|(0,0)=(0,0, b)
Wherein, SuIt is curved surface S to the first derivative of parameter u, SuuIt is second dervative;SvIt is curved surface S to the first derivative of parameter v,
SvvIt is second dervative;SuvIt is curved surface S to parameter u, the second dervative of v;N is sweared for the per unit system of curved surface;
According to above parameter, can calculate:
E=Su·Su=1+d2And F=Su·Sv=de
G=Sv·Sv=1+e2And
And
Wherein, E, F, G are the first fundamental quantity of curved surface, and L, M, N are the second fundamental quantity of curved surface;
Then the Gaussian curvature and average curvature values at P points are:
Minimum principal curvatures θminWith maximum principal curvatures θmaxComputing formula is respectively:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611088463.3A CN106780509A (en) | 2016-12-01 | 2016-12-01 | Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611088463.3A CN106780509A (en) | 2016-12-01 | 2016-12-01 | Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106780509A true CN106780509A (en) | 2017-05-31 |
Family
ID=58915305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611088463.3A Pending CN106780509A (en) | 2016-12-01 | 2016-12-01 | Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106780509A (en) |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194994A (en) * | 2017-06-16 | 2017-09-22 | 广东工业大学 | A kind of method and device without the demarcation surface points cloud data reconstruction face of cylinder |
CN107492120A (en) * | 2017-07-18 | 2017-12-19 | 北京航空航天大学 | Point cloud registration method |
CN108090960A (en) * | 2017-12-25 | 2018-05-29 | 北京航空航天大学 | A kind of Object reconstruction method based on geometrical constraint |
CN108241871A (en) * | 2017-12-27 | 2018-07-03 | 华北水利水电大学 | Laser point cloud and visual fusion data classification method based on multiple features |
CN108846888A (en) * | 2018-04-23 | 2018-11-20 | 北京建筑大学 | A kind of Gu wooden building components fine size information automation extracting method |
CN109118500A (en) * | 2018-07-16 | 2019-01-01 | 重庆大学产业技术研究院 | A kind of dividing method of the Point Cloud Data from Three Dimension Laser Scanning based on image |
CN109191580A (en) * | 2018-07-26 | 2019-01-11 | 大连交通大学 | A kind of scoliosis orthopedic device digitalized design method |
CN109490072A (en) * | 2018-10-09 | 2019-03-19 | 广东交通职业技术学院 | A kind of civil engineering work detection system and its detection method |
CN110009726A (en) * | 2019-03-08 | 2019-07-12 | 浙江中海达空间信息技术有限公司 | A method of according to the structural relation between plane primitive to data reduction plane |
CN110147775A (en) * | 2019-05-24 | 2019-08-20 | 北京建筑大学 | Utilize refinement method of the space separation method from data reduction indoor navigation element |
CN110222742A (en) * | 2019-05-23 | 2019-09-10 | 星际空间(天津)科技发展有限公司 | Based on point cloud segmentation method, apparatus, storage medium and the equipment for being layered more echoes |
CN110400322A (en) * | 2019-07-30 | 2019-11-01 | 江南大学 | Fruit point cloud segmentation method based on color and three-dimensional geometric information |
CN110910446A (en) * | 2019-11-26 | 2020-03-24 | 北京拓维思科技有限公司 | Method and device for determining building removal area and method and device for determining indoor area of building |
CN111462123A (en) * | 2020-03-30 | 2020-07-28 | 华南理工大学 | Point cloud data segmentation method based on spectral clustering |
CN111692991A (en) * | 2020-06-02 | 2020-09-22 | 哈尔滨工程大学 | Point cloud data acquisition method for measuring batten bonding surface based on white light interference |
WO2020206671A1 (en) * | 2019-04-09 | 2020-10-15 | 北京大学深圳研究生院 | Attribute-based point cloud strip division method |
CN111815776A (en) * | 2020-02-04 | 2020-10-23 | 山东水利技师学院 | Three-dimensional building fine geometric reconstruction method integrating airborne and vehicle-mounted three-dimensional laser point clouds and streetscape images |
CN112070769A (en) * | 2020-09-18 | 2020-12-11 | 福州大学 | Layered point cloud segmentation method based on DBSCAN |
CN112241676A (en) * | 2020-07-07 | 2021-01-19 | 西北农林科技大学 | Method for automatically identifying terrain sundries |
CN112560573A (en) * | 2020-10-29 | 2021-03-26 | 河北省地震局 | Building earthquake damage information detection and extraction method |
CN112686916A (en) * | 2020-12-28 | 2021-04-20 | 淮阴工学院 | Curved surface reconstruction system based on heterogeneous multi-region CT scanning data processing |
CN112785596A (en) * | 2021-02-01 | 2021-05-11 | 中国铁建电气化局集团有限公司 | Dot cloud picture bolt segmentation and height measurement method based on DBSCAN clustering |
CN112907722A (en) * | 2019-11-19 | 2021-06-04 | 广东博智林机器人有限公司 | Building information model generation method, system, device and storage medium |
CN113129393A (en) * | 2020-01-15 | 2021-07-16 | 上海交通大学 | Point cloud data processing method and system |
CN113379748A (en) * | 2020-03-09 | 2021-09-10 | 北京京东乾石科技有限公司 | Point cloud panorama segmentation method and device |
CN113436223A (en) * | 2021-07-14 | 2021-09-24 | 北京市测绘设计研究院 | Point cloud data segmentation method and device, computer equipment and storage medium |
CN113804118A (en) * | 2021-08-16 | 2021-12-17 | 长江水利委员会长江科学院 | Building deformation monitoring method based on three-dimensional laser point cloud geometric features |
CN115346019A (en) * | 2022-09-06 | 2022-11-15 | 南京航空航天大学 | Method, device and system for measuring geometrical parameters of point cloud circular hole |
CN115422387A (en) * | 2022-11-04 | 2022-12-02 | 山东矩阵软件工程股份有限公司 | Point cloud data processing method and system based on multi-dimensional point cloud fusion data |
CN115761023A (en) * | 2022-12-02 | 2023-03-07 | 同济大学 | Three-dimensional point cloud compression system and method based on point cloud matrix singular value characteristics |
CN115797551A (en) * | 2022-11-14 | 2023-03-14 | 国网湖北省电力有限公司超高压公司 | Laser point cloud data automatic modeling method based on two-step unsupervised clustering algorithm |
CN116128886A (en) * | 2023-04-18 | 2023-05-16 | 深圳市其域创新科技有限公司 | Point cloud data segmentation method and device, electronic equipment and storage medium |
CN116310115A (en) * | 2023-03-17 | 2023-06-23 | 合肥泰瑞数创科技有限公司 | Method and system for constructing building three-dimensional model based on laser point cloud |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530899A (en) * | 2013-10-10 | 2014-01-22 | 浙江万里学院 | Geometric featuer-based point cloud simplification method |
CN104573705A (en) * | 2014-10-13 | 2015-04-29 | 北京建筑大学 | Clustering method for building laser scan point cloud data |
CN105844629A (en) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | Automatic segmentation method for point cloud of facade of large scene city building |
-
2016
- 2016-12-01 CN CN201611088463.3A patent/CN106780509A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530899A (en) * | 2013-10-10 | 2014-01-22 | 浙江万里学院 | Geometric featuer-based point cloud simplification method |
CN104573705A (en) * | 2014-10-13 | 2015-04-29 | 北京建筑大学 | Clustering method for building laser scan point cloud data |
CN105844629A (en) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | Automatic segmentation method for point cloud of facade of large scene city building |
Non-Patent Citations (3)
Title |
---|
刘颖: "基于半监督集成支持向量机的土地覆盖遥感分类方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
张强: "基于几何特征的规则激光点云分割方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王春晓: "密度聚类方法在点云数据分割中的应用研究", 《测绘与空间地理信息》 * |
Cited By (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194994B (en) * | 2017-06-16 | 2020-12-15 | 广东工业大学 | Method and device for reconstructing cylindrical surface by using point cloud data without calibration curved surface |
CN107194994A (en) * | 2017-06-16 | 2017-09-22 | 广东工业大学 | A kind of method and device without the demarcation surface points cloud data reconstruction face of cylinder |
CN107492120A (en) * | 2017-07-18 | 2017-12-19 | 北京航空航天大学 | Point cloud registration method |
CN107492120B (en) * | 2017-07-18 | 2020-04-28 | 北京航空航天大学 | Point cloud registration method |
CN108090960A (en) * | 2017-12-25 | 2018-05-29 | 北京航空航天大学 | A kind of Object reconstruction method based on geometrical constraint |
CN108241871A (en) * | 2017-12-27 | 2018-07-03 | 华北水利水电大学 | Laser point cloud and visual fusion data classification method based on multiple features |
CN108846888A (en) * | 2018-04-23 | 2018-11-20 | 北京建筑大学 | A kind of Gu wooden building components fine size information automation extracting method |
CN108846888B (en) * | 2018-04-23 | 2022-03-29 | 北京建筑大学 | Automatic extraction method for fine size information of ancient wood building components |
CN109118500B (en) * | 2018-07-16 | 2022-05-10 | 重庆大学产业技术研究院 | Image-based three-dimensional laser scanning point cloud data segmentation method |
CN109118500A (en) * | 2018-07-16 | 2019-01-01 | 重庆大学产业技术研究院 | A kind of dividing method of the Point Cloud Data from Three Dimension Laser Scanning based on image |
CN109191580A (en) * | 2018-07-26 | 2019-01-11 | 大连交通大学 | A kind of scoliosis orthopedic device digitalized design method |
CN109490072A (en) * | 2018-10-09 | 2019-03-19 | 广东交通职业技术学院 | A kind of civil engineering work detection system and its detection method |
CN110009726A (en) * | 2019-03-08 | 2019-07-12 | 浙江中海达空间信息技术有限公司 | A method of according to the structural relation between plane primitive to data reduction plane |
CN110009726B (en) * | 2019-03-08 | 2022-09-30 | 浙江中海达空间信息技术有限公司 | Method for extracting plane from point cloud according to structural relationship between plane elements |
WO2020206671A1 (en) * | 2019-04-09 | 2020-10-15 | 北京大学深圳研究生院 | Attribute-based point cloud strip division method |
US12039762B2 (en) | 2019-04-09 | 2024-07-16 | Peking University Shenzhen Graduate School | Attribute-based point cloud strip division method |
CN110222742A (en) * | 2019-05-23 | 2019-09-10 | 星际空间(天津)科技发展有限公司 | Based on point cloud segmentation method, apparatus, storage medium and the equipment for being layered more echoes |
CN110222742B (en) * | 2019-05-23 | 2022-12-02 | 星际空间(天津)科技发展有限公司 | Point cloud segmentation method, device, storage medium and equipment based on layered multi-echo |
CN110147775B (en) * | 2019-05-24 | 2021-05-18 | 北京建筑大学 | Method for extracting indoor navigation elements from point cloud by utilizing refined space separation method |
CN110147775A (en) * | 2019-05-24 | 2019-08-20 | 北京建筑大学 | Utilize refinement method of the space separation method from data reduction indoor navigation element |
CN110400322A (en) * | 2019-07-30 | 2019-11-01 | 江南大学 | Fruit point cloud segmentation method based on color and three-dimensional geometric information |
CN110400322B (en) * | 2019-07-30 | 2021-03-16 | 江南大学 | Fruit point cloud segmentation method based on color and three-dimensional geometric information |
CN112907722A (en) * | 2019-11-19 | 2021-06-04 | 广东博智林机器人有限公司 | Building information model generation method, system, device and storage medium |
CN110910446A (en) * | 2019-11-26 | 2020-03-24 | 北京拓维思科技有限公司 | Method and device for determining building removal area and method and device for determining indoor area of building |
CN113129393A (en) * | 2020-01-15 | 2021-07-16 | 上海交通大学 | Point cloud data processing method and system |
CN113129393B (en) * | 2020-01-15 | 2022-12-27 | 上海交通大学 | Point cloud data processing method and system |
CN111815776A (en) * | 2020-02-04 | 2020-10-23 | 山东水利技师学院 | Three-dimensional building fine geometric reconstruction method integrating airborne and vehicle-mounted three-dimensional laser point clouds and streetscape images |
CN113379748B (en) * | 2020-03-09 | 2024-03-01 | 北京京东乾石科技有限公司 | Point cloud panorama segmentation method and device |
CN113379748A (en) * | 2020-03-09 | 2021-09-10 | 北京京东乾石科技有限公司 | Point cloud panorama segmentation method and device |
CN111462123A (en) * | 2020-03-30 | 2020-07-28 | 华南理工大学 | Point cloud data segmentation method based on spectral clustering |
CN111462123B (en) * | 2020-03-30 | 2023-06-20 | 华南理工大学 | Point cloud data segmentation method based on spectral clustering |
CN111692991A (en) * | 2020-06-02 | 2020-09-22 | 哈尔滨工程大学 | Point cloud data acquisition method for measuring batten bonding surface based on white light interference |
CN111692991B (en) * | 2020-06-02 | 2021-09-10 | 哈尔滨工程大学 | Point cloud data acquisition method for measuring batten bonding surface based on white light interference |
CN112241676A (en) * | 2020-07-07 | 2021-01-19 | 西北农林科技大学 | Method for automatically identifying terrain sundries |
CN112070769A (en) * | 2020-09-18 | 2020-12-11 | 福州大学 | Layered point cloud segmentation method based on DBSCAN |
CN112070769B (en) * | 2020-09-18 | 2022-06-03 | 福州大学 | Layered point cloud segmentation method based on DBSCAN |
CN112560573A (en) * | 2020-10-29 | 2021-03-26 | 河北省地震局 | Building earthquake damage information detection and extraction method |
CN112560573B (en) * | 2020-10-29 | 2023-03-28 | 河北省地震局 | Building earthquake damage information detection and extraction method |
CN112686916A (en) * | 2020-12-28 | 2021-04-20 | 淮阴工学院 | Curved surface reconstruction system based on heterogeneous multi-region CT scanning data processing |
CN112686916B (en) * | 2020-12-28 | 2024-04-05 | 淮阴工学院 | Curved surface reconstruction system based on heterogeneous multi-region CT scanning data processing |
CN112785596B (en) * | 2021-02-01 | 2022-06-10 | 中国铁建电气化局集团有限公司 | Dot cloud picture bolt segmentation and height measurement method based on DBSCAN clustering |
CN112785596A (en) * | 2021-02-01 | 2021-05-11 | 中国铁建电气化局集团有限公司 | Dot cloud picture bolt segmentation and height measurement method based on DBSCAN clustering |
CN113436223A (en) * | 2021-07-14 | 2021-09-24 | 北京市测绘设计研究院 | Point cloud data segmentation method and device, computer equipment and storage medium |
CN113804118A (en) * | 2021-08-16 | 2021-12-17 | 长江水利委员会长江科学院 | Building deformation monitoring method based on three-dimensional laser point cloud geometric features |
CN115346019A (en) * | 2022-09-06 | 2022-11-15 | 南京航空航天大学 | Method, device and system for measuring geometrical parameters of point cloud circular hole |
CN115422387A (en) * | 2022-11-04 | 2022-12-02 | 山东矩阵软件工程股份有限公司 | Point cloud data processing method and system based on multi-dimensional point cloud fusion data |
CN115797551A (en) * | 2022-11-14 | 2023-03-14 | 国网湖北省电力有限公司超高压公司 | Laser point cloud data automatic modeling method based on two-step unsupervised clustering algorithm |
CN115797551B (en) * | 2022-11-14 | 2023-11-03 | 国网湖北省电力有限公司超高压公司 | Automatic modeling method for laser point cloud data based on two-step unsupervised clustering algorithm |
CN115761023A (en) * | 2022-12-02 | 2023-03-07 | 同济大学 | Three-dimensional point cloud compression system and method based on point cloud matrix singular value characteristics |
CN116310115B (en) * | 2023-03-17 | 2023-11-24 | 合肥泰瑞数创科技有限公司 | Method and system for constructing building three-dimensional model based on laser point cloud |
CN116310115A (en) * | 2023-03-17 | 2023-06-23 | 合肥泰瑞数创科技有限公司 | Method and system for constructing building three-dimensional model based on laser point cloud |
CN116128886A (en) * | 2023-04-18 | 2023-05-16 | 深圳市其域创新科技有限公司 | Point cloud data segmentation method and device, electronic equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106780509A (en) | Merge the building object point cloud layer time cluster segmentation method of multidimensional characteristic | |
US10755098B2 (en) | Evaluation method of solar energy utilization potential in urban high-density areas based on low-altitude photogrammetry | |
CN110570428B (en) | Method and system for dividing building roof sheet from large-scale image dense matching point cloud | |
CN110111414B (en) | Orthographic image generation method based on three-dimensional laser point cloud | |
Zhou et al. | 2.5 d dual contouring: A robust approach to creating building models from aerial lidar point clouds | |
CN105335966B (en) | Multiscale morphology image division method based on local homogeney index | |
CN105069746B (en) | Video real-time face replacement method and its system based on local affine invariant and color transfer technology | |
Yang et al. | A shape-based segmentation method for mobile laser scanning point clouds | |
Ientilucci et al. | Advances in wide-area hyperspectral image simulation | |
CN109887082A (en) | A kind of interior architecture three-dimensional modeling method and device based on point cloud data | |
CN109446691B (en) | Living standing tree wind resistance performance analysis method based on laser point cloud and aerodynamics | |
CN111047695B (en) | Method for extracting height spatial information and contour line of urban group | |
CN109492852A (en) | A kind of detection method for quality of the water conservancy project structure based on BIM | |
CN106097311A (en) | The building three-dimensional rebuilding method of airborne laser radar data | |
CN106199557A (en) | A kind of airborne laser radar data vegetation extracting method | |
CN106126816B (en) | Repeat the extensive ALS building point cloud modeling method under building automatic sensing | |
CN107292234A (en) | It is a kind of that method of estimation is laid out based on information edge and the indoor scene of multi-modal feature | |
CN107170037A (en) | A kind of real-time three-dimensional point cloud method for reconstructing and system based on multiple-camera | |
CN108154104A (en) | A kind of estimation method of human posture based on depth image super-pixel union feature | |
CN109815847A (en) | A kind of vision SLAM method based on semantic constraint | |
CN115564926A (en) | Three-dimensional patch model construction method based on image building structure learning | |
CN103679210A (en) | Ground object recognition method based on hyperspectral image unmixing | |
CN110363299A (en) | Space reasoning by cases method towards delamination-terrane of appearing | |
CN110207670A (en) | A method of artificial forest forest hat width parameter is obtained using two dimensional image | |
CN105913444A (en) | Livestock body contour reconstruction method and body condition scoring method based on soft laser ranging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |
|
RJ01 | Rejection of invention patent application after publication |