CN107886565B - Method for ordering rock section optical scanning disordered point cloud - Google Patents
Method for ordering rock section optical scanning disordered point cloud Download PDFInfo
- Publication number
- CN107886565B CN107886565B CN201711153960.1A CN201711153960A CN107886565B CN 107886565 B CN107886565 B CN 107886565B CN 201711153960 A CN201711153960 A CN 201711153960A CN 107886565 B CN107886565 B CN 107886565B
- Authority
- CN
- China
- Prior art keywords
- point
- coordinate
- array
- point cloud
- coordinates
- 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.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000003287 optical effect Effects 0.000 title claims abstract description 16
- 238000005070 sampling Methods 0.000 claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 238000003672 processing method Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Generation (AREA)
Abstract
The invention discloses a rock section optical scanning disordered point cloud ordering processing method. The method comprises the steps of obtaining a rock fracture surface point cloud file, calculating the point cloud length, adjusting deviation, readjusting coordinates, selecting a calculation range, obtaining a z coordinate of a sampling point through linear interpolation and the like. The method can select an ordering processing region, can arrange the ordered point cloud according to the coordinate sequence, and can directly calculate the subsequent fracture surface parameters without readjusting. After the method is adopted to carry out the disordered point cloud ordering treatment, the shape of the rock section is not changed, and the method has stronger practicability in the research of the shape of the rock joint surface and mechanics.
Description
Technical Field
The invention relates to the field of rock fracture mechanics, in particular to a rock section optical scanning disordered point cloud ordering processing method.
Background
The rock is broken under the action of tension, compression and shearing force to form a rock fracture surface. The shape of the fracture surface of the rock has important influence on the strength and deformation characteristics of the joint surface, and is a hot point of research of current scholars. At present, the calculation of the profile morphology characteristics mainly focuses on the three-dimensional average dip angle, roughness parameters, apparent anisotropy, average gradient distortion parameters, JRC and fractal characteristics of a fracture surface. The basic principle is that an optical scanner is used for scanning a rock section to obtain x, y and z coordinates of a series of points on the surface of the fracture surface, namely point cloud, and then calculation is carried out according to the point cloud data of the section. The rough parameters and fractal calculation of the rock fracture surface are carried out based on the ordered point cloud of the rock fracture surface. However, many optical scanners produce point cloud data as a disordered point cloud, which must be converted to an ordered point cloud in the above section parameter calculations.
The existing ordering processing method adopts a local point cloud quadric surface fitting method to calculate the coordinates of a target point, and the processed ordered point cloud is a smooth curved surface and is not suitable for the analysis of the topography characteristics of a rock section.
Disclosure of Invention
The invention aims to provide a rock section disordered point cloud ordering processing method which is simple in algorithm, reliable in result and capable of selecting a processing area, so as to solve the problem that the prior art cannot adopt disordered point cloud to calculate the shape and appearance characteristics of the rock section.
The technical scheme adopted for achieving the purpose of the invention is that the method for ordering the disordered point cloud by the optical scanning of the rock section comprises the following steps:
1) and (3) carrying out optical scanning on the rock fracture surface by adopting an optical scanner to obtain x, y and z coordinates of the rock fracture surface series points in a coordinate system o-xyz. And exporting scanning data to obtain a rock fracture surface point cloud file. Wherein, the x coordinate direction is the rock sample length direction, the y coordinate direction is the rock sample width direction, and the z coordinate is the height of point.
2) And reading the point cloud file, and assigning the data of the point cloud to an array Original. The first row of the array origin is a point cloud x coordinate, the second row is a point cloud y coordinate, and the third row is a point cloud z coordinate.
3) Traversing the first column and the second column of the array origin, searching the maximum value and the minimum value of the x coordinate and the y coordinate, and calculating according to the formula (1) to obtain the length x of the point cloud in the x directionwidthAnd length y in the y directionwidth. Wherein,
in the formula, xmaxIs the maximum value of the x coordinate, xminIs the minimum value of the x coordinate, ymaxIs the maximum value of the y coordinate, yminIs the minimum value of the y coordinate.
4) And adjusting the deviation of the scanning data and the actual sample and readjusting the coordinates to enable the center coordinates of the sample to be zero. And calculating the adjusted x coordinate x according to the formula (2)afteramendAnd the adjusted y-coordinate yafteramend。
In the formula, specimenxlength is the length of the sample in the x direction. specimenylength is the length of the specimen in the y-direction. x is the number ofratioThe length of the sample in the x direction specimenxlength and the length x of the scanning point cloud in the x direction are shown in the specificationwidthThe ratio of (a) to (b). y isratioFor the length y of the sample in the y directionafteramendLength y of scanning point cloud in y directionwidthThe ratio of (a) to (b).
5) The adjusted point coordinates are recorded in the array Amend. Wherein the first column of the array Amend is the adjusted x coordinate xafteramendAnd the second column is the adjusted y coordinate yafteramendAnd the third column is the z coordinate of the point corresponding to the x, y coordinates.
6) And selecting a calculation range. Traversing the array Amend, and writing the ith row of the array Amend into the array Calculatarea if the point represented by the ith row is in the calculation range. The calculation range is a rectangular area which takes the center of the sample as the origin of coordinates and is defined by x and y coordinates. The minimum value of the x coordinate of the rectangular area is mincalvulatex, the maximum value of the x coordinate is maxcalculatatex, the minimum value of the y coordinate is mincalvulatey, and the maximum value of the y coordinate is maxcalculatatey.
7) And performing linear interpolation according to the sampling interval of the ordered point cloud to obtain the z coordinate of the sampling point, and writing the z coordinate into an array Final. Wherein the array Final comprises xsamplingumber rows and ysamplingnumber columns. xsamplingumber is the number of samples of the ordered point cloud in the x-direction, and ysamplingnumber is the number of samples of the ordered point cloud in the y-direction. The values of xsamplingnumber and ysampllingnnumber are calculated using equation (3):
in the formula, xsamplingstep is the sampling interval of the ordered point cloud in the x direction, and ysamplingstep is the sampling interval of the ordered point cloud in the y direction.
Further, the implementation method of the linear interpolation in step 7) is as follows:
1) and (4) calculating the coordinates (samplingx and samplingy) of the sampling point at the ith row and the jth column of the ordered point cloud. Wherein,
2) and (5) calculating the distance D between the projection of each point in the point cloud in the xoy plane and the projection of the sampling point in the xoy plane according to the formula (5). And writes the x, y, z coordinates and Distance D of the point into the 1 st to 4 th columns of the array Distance, respectively. Wherein,
3) traversing the Distance of the array, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point A; and records the x, y, z coordinates of point A and the row number R in the Distance arrayA(ii) a Assigning the array Distance to the array Interdistance, and assigning the array Interdistance to the 4 th column R in the array InterdistanceAThe row has a large value.
4) Traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point B; and records the x, y, z coordinates of point B and the line number R in the Interdistance arrayB(ii) a Assigning to the 4 th column R in the array InterDistanceBThe row has a large value.
5) Checking the x coordinates and the y coordinates of the sampling point, the point A and the point B; if the x coordinate and the y coordinate of the point B are the same as those of the point A, repeating the step 4) until the x coordinate or the y coordinate of the point B is different from that of the point A.
6) Traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point C; and records the x, y, z coordinates of point C and the line number R in the Interdistance arrayC(ii) a Assigning to the 4 th column R in the array InterDistanceCThe row has a large value.
7) And calculating projection areas S1, S2 and S3 of three triangles with the sampling point, the point A and the point B, the sampling point, the point B and the point C as vertexes in the xoy plane, and calculating the projection area S of the triangle with the point A, the point B and the point C as vertexes in the xoy plane.
8) If S is 0 or S1+ S2+ S3 ≠ S, steps 6) and 7) are repeated until S > 0 and S1+ S2+ S3 ═ S.
9) And solving a plane equation passing through the point A, the point B and the point C at the same time.
10) And substituting the x and y coordinates of the sampling points into the plane equation in the step 9), and solving the z coordinate of the sampling points.
11) Writing the z coordinate of the sampling point in the step 10) into the ith row and the jth column of the array Final.
12) Repeating the steps 1) to 11) to calculate the z coordinate of the next sampling point until the last sampling point. The data in the array Final is the z coordinate of the ordered points with the x-direction interval of xsampling and the y-direction interval of ysampling.
The technical effects of the invention are undoubted:
A. the principle is simple, the result is reliable, and disordered point cloud can be converted into ordered point cloud;
B. the sampling interval of the ordered point cloud can be freely defined and does not depend on the existing unordered point cloud;
C. the ordering processing area can be selected freely through the coordinate;
D. the ordered point clouds can be arranged according to the coordinate sequence, and the subsequent fracture surface parameters can be directly calculated without readjustment.
Drawings
FIG. 1 is a flow chart of a point cloud ordering process;
FIG. 2 is a schematic diagram of an optical scanning disordered point cloud image of a rock section;
FIG. 3 is an enlarged view of the partial point cloud of FIG. 2;
FIG. 4 is a diagram showing the relative position relationship of the projection of the sampling point to the point A, the point B and the point C in the xoy plane;
FIG. 5 is a schematic diagram of an image after the ordering process of the disordered point cloud;
fig. 6 is an enlarged view of the partial point cloud of fig. 5.
Detailed Description
The present invention is further illustrated by the following examples, but it should not be construed that the scope of the above-described subject matter is limited to the following examples. Various substitutions and alterations can be made without departing from the technical idea of the invention and the scope of the invention is covered by the present invention according to the common technical knowledge and the conventional means in the field.
The embodiment discloses a rock section optical scanning disordered point cloud ordering processing method, and referring to fig. 1, the method comprises the following steps:
1) and (3) carrying out optical scanning on the rock fracture surface by adopting an optical scanner to obtain x, y and z coordinates of the rock fracture surface series points in a coordinate system o-xyz. And exporting scanning data to obtain a rock fracture surface point cloud file. The schematic diagram of the disordered point cloud image and its partial enlarged view are shown in fig. 2 and 3. Wherein, the x coordinate direction is the rock sample length direction, the y coordinate direction is the rock sample width direction, and the z coordinate is the height of point.
2) And reading the point cloud file, and assigning the data of the point cloud to an array Original. The first row of the array origin is a point cloud x coordinate, the second row is a point cloud y coordinate, and the third row is a point cloud z coordinate.
3) Traversing the first column and the second column of the array origin, searching the maximum value and the minimum value of the x coordinate and the y coordinate, and calculating according to the formula (1) to obtain the length x of the point cloud in the x directionwidthAnd length y in the y directionwidth. Wherein,
in the formula, xmaxIs the maximum value of the x coordinate, xminIs the minimum value of the x coordinate, ymaxIs the maximum value of the y coordinate, yminIs the minimum value of the y coordinate.
4) And adjusting the deviation of the scanning data and the actual sample and readjusting the coordinates to enable the center coordinates of the sample to be zero. And calculating the adjusted x coordinate x according to the formula (2)afteramendAnd the adjusted y-coordinate yafteramend。
In the formula, specimenxlength is the length of the sample in the x direction. specimenylength is the length of the specimen in the y-direction. x is the number ofratioThe length of the sample in the x direction specimenxlength and the length x of the scanning point cloud in the x direction are shown in the specificationwidthThe ratio of (a) to (b). y isratioFor the length y of the sample in the y directionafteramendLength y of scanning point cloud in y directionwidthThe ratio of (a) to (b).
5) The adjusted point coordinates are recorded in the array Amend. Wherein the first column of the array Amend is the adjusted x coordinate xafteramendAnd the second column is the adjusted y coordinate yafteramendAnd the third column is the z coordinate of the point corresponding to the x, y coordinates.
6) And selecting a calculation range, namely selecting an ordering processing area. Traversing the array Amend, and writing the ith row of the array Amend into the array Calculatarea if the point represented by the ith row is in the calculation range. The calculation range is a rectangular area which takes the center of the sample as the origin of coordinates and is defined by x and y coordinates. The minimum value of the x coordinate of the rectangular area is mincalvulatex, the maximum value of the x coordinate is maxcalculatatex, the minimum value of the y coordinate is mincalvulatey, and the maximum value of the y coordinate is maxcalculatatey.
7) And setting a sampling interval, performing linear interpolation according to the ordered point cloud sampling interval to obtain a z coordinate of a sampling point, and writing the z coordinate into an array Final. Wherein the array Final comprises xsamplingumber rows and ysamplingnumber columns. xsamplingumber is the number of samples of the ordered point cloud in the x-direction, and ysamplingnumber is the number of samples of the ordered point cloud in the y-direction. The values of xsamplingnumber and ysampllingnnumber are calculated using equation (3):
in the formula, xsamplingstep is the sampling interval of the ordered point cloud in the x direction, and ysamplingstep is the sampling interval of the ordered point cloud in the y direction.
The implementation method of the linear interpolation in the step 7) comprises the following steps:
1) and (4) calculating the coordinates (samplingx and samplingy) of the sampling point at the ith row and the jth column of the ordered point cloud. Wherein,
2) and (5) calculating the distance D between the projection of each point in the point cloud in the xoy plane and the projection of the sampling point in the xoy plane according to the formula (5). And writes the x, y, z coordinates and Distance D of the point into the 1 st to 4 th columns of the array Distance, respectively. Wherein,
3) traversing the Distance of the array, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point A; and records the x, y, z coordinates of point A and the row number R in the Distance arrayA(ii) a Assigning the array Distance to the array Interdistance, and assigning the array Interdistance to the 4 th column R in the array InterdistanceAThe row has a large value. The step can find the point A with the minimum projection distance with the sampling point in the xoy plane.
4) Traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point B; and records the x, y, z coordinates of point B and the line number R in the Interdistance arrayB(ii) a Assigning to the 4 th column R in the array InterDistanceBThe row has a large value.
5) Checking the x coordinates and the y coordinates of the sampling point, the point A and the point B; if the x coordinate and the y coordinate of the point B are the same as those of the point A, repeating the step 4) until the x coordinate or the y coordinate of the point B is different from that of the point A.
Through the steps 4) and 5), the point B with the second smallest projection distance with the sampling point in the xoy plane can be found, and the point B, the sampling point and the projection of the point A in the xoy plane are not overlapped.
6) Traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point C; and recording the x, y, z coordinates of point C and the position of point C in InterLine number R in Distance arrayC(ii) a Assigning to the 4 th column R in the array InterDistanceCThe row has a large value.
7) Referring to fig. 4, projection areas S1, S2, and S3 in the xoy plane of three triangles having the sampling point, the point a, the point B, the sampling point, the point B, the point C, and the point a as vertexes are calculated, respectively, and a projection area S in the xoy plane of a triangle having the point a, the point B, and the point C as vertexes is calculated.
8) If S is 0 or S1+ S2+ S3 ≠ S, steps 6) and 7) are repeated until S > 0 and S1+ S2+ S3 ═ S.
Through the steps 6), 7) and 8), the point C with the third small projection distance from the sampling point in the xoy plane can be found. Points A, B, C are not collinear, and there is only one plane through point A, B, C. The projections of the sample points in the xoy plane are inside or on the boundaries of a triangle having as vertices the projections of points a, B and C in the xoy plane.
9) And solving a plane equation passing through the point A, the point B and the point C at the same time.
10) And substituting the x and y coordinates of the sampling points into the plane equation in the step 9, and solving the z coordinate of the sampling points.
11) Writing the z coordinate of the sampling point in the step 10) into the ith row and the jth column of the array Final.
12) Repeating the steps 1) to 11) to calculate the z coordinate of the next sampling point until the last sampling point. The data in the array Final is the z coordinate of the ordered points with the x-direction interval of xsampling and the y-direction interval of ysampling.
The schematic diagram of the ordered point cloud and its enlarged partial point cloud obtained through the above steps are shown in fig. 5 and 6. Referring to fig. 2 and 5, the cross-sectional morphology does not change before and after the ordering treatment, the influence on the calculation result is small, and the method has strong practicability in the morphology of the rock joint surface and the mechanical research.
Claims (1)
1. An optical scanning disordered point cloud ordering processing method for a rock section is characterized by comprising the following steps:
1) performing optical scanning on the rock fracture surface by adopting an optical scanner to obtain x, y and z coordinates of the rock fracture surface series points in a coordinate system o-xyz; exporting scanning data to obtain a rock fracture surface point cloud file; the method comprises the following steps of firstly, obtaining a point cloud file, wherein the x coordinate direction is the length direction of a rock sample, the y coordinate direction is the width direction of the rock sample, and the z coordinate is the height of the midpoint of the point cloud file;
2) reading a point cloud file, and assigning the data of the point cloud to an array Original; wherein, the first row of the array origin is a point cloud x coordinate, the second row is a point cloud y coordinate, and the third row is a point cloud z coordinate;
3) traversing the first column and the second column of the array origin, searching the maximum value and the minimum value of the x coordinate and the y coordinate, and calculating according to the formula (1) to obtain the length x of the point cloud in the x directionwidthAnd length y in the y directionwidth;
In the formula, xmaxIs the maximum value of the x coordinate, xminIs the minimum value of the x coordinate, ymaxIs the maximum value of the y coordinate, yminIs the minimum value of the y coordinate;
4) adjusting the deviation of the scanning data and the actual sample and readjusting the coordinates to make the central coordinate of the sample zero; and calculating the adjusted x coordinate x according to the formula (2)afteramendAnd the adjusted y-coordinate yafteramend;
In the formula, specimenxlength is the length of the sample in the x direction; the specimenylength is the length of the sample in the y direction; x is the number ofratioThe length of the sample in the x direction specimenxlength and the length x of the scanning point cloud in the x direction are shown in the specificationwidthThe ratio of (A) to (B); y isratioFor the length y of the sample in the y directionafteramendLength y of scanning point cloud in y directionwidthThe ratio of (A) to (B);
5) recording the adjusted point coordinates into an array Amend; wherein the first column of the array Amend is toneIntegrated x coordinate xafteramendAnd the second column is the adjusted y coordinate yafteramendThe third column is the z coordinate of the point corresponding to the x and y coordinates;
6) selecting a calculation range; traversing the array Amend, and writing the ith row of the array Amend into an array Calculatarea if the point represented by the ith row is within the calculation range; the calculation range is a rectangular area which takes the center of the sample as the origin of coordinates and is defined by x and y coordinates; the minimum value of the x coordinate of the rectangular area is mincalvulatex, the maximum value of the x coordinate is maxcalculatex, the minimum value of the y coordinate is mincalvulatey, and the maximum value of the y coordinate is maxcalculatey;
7) performing linear interpolation according to the sampling interval of the ordered point cloud to obtain the z coordinate of the sampling point, and writing the z coordinate into an array Final; wherein the array Final comprises xsamplingumber rows and ysamplingnumber columns; xsamplingumber is the sampling number of the ordered point cloud in the x direction, and ysamplingnumber is the sampling number of the ordered point cloud in the y direction; the values of xsamplingnumber and ysampllingnnumber are calculated using equation (3):
in the formula, xsamplingstep is the sampling interval of the ordered point cloud in the x direction, and ysamplingstep is the sampling interval of the ordered point cloud in the y direction;
step 7) the method for realizing the linear interpolation comprises the following steps:
7-1) calculating coordinates (samplingx and samplingy) of a sampling point at the jth column of the ith row of the ordered point cloud according to a formula (4); wherein,
7-2) calculating the distance D between the projection of each point in the point cloud in the xoy plane and the projection of the sampling point in the xoy plane according to a formula (5); respectively writing the x, y, z coordinates and the Distance D of the point into the 1 st to 4 th columns of the array Distance; wherein,
7-3) traversing the Distance of the array, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point A; and records the x, y, z coordinates of point A and the row number R in the Distance arrayA(ii) a Assigning the array Distance to the array Interdistance, and assigning the array Interdistance to the 4 th column R in the array InterdistanceALine one large value;
7-4) traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point B; and records the x, y, z coordinates of point B and the line number R in the Interdistance arrayB(ii) a Assigning to the 4 th column R in the array InterDistanceBLine one large value;
7-5) checking the x coordinates and the y coordinates of the sampling point, the point A and the point B; if the x coordinate and the y coordinate of the point B are the same as those of the point A, repeating the step 7-4) until the x coordinate or the y coordinate of the point B is different from that of the point A;
7-6) traversing the array InterDistance, searching the minimum value in the 4 th column of the array, and recording the corresponding point as a point C; and records the x, y, z coordinates of point C and the line number R in the Interdistance arrayC(ii) a Assigning to the 4 th column R in the array InterDistanceCLine one large value;
7-7) respectively calculating projection areas S1, S2 and S3 of three triangles taking the sampling point, the point A and the point B, the sampling point, the point B and the point C as vertexes in the xoy plane, and calculating the projection area S of the triangle taking the point A, the point B and the point C as vertexes in the xoy plane;
7-8) if S ═ 0 or S1+ S2+ S3 ≠ S, repeating steps 7-6) and 7-7) until S > 0 and S1+ S2+ S3 ═ S;
7-9) solving a plane equation passing through the point A, the point B and the point C simultaneously;
7-10) substituting the x and y coordinates of the sampling points into the plane equation in the step 7-9), and solving the z coordinate of the sampling points;
7-11) writing the z coordinate of the sampling point in the step 7-10) into the ith row and the jth column of the array Final;
7-12) repeating the steps 7-1) to 7-11) to calculate the z coordinate of the next sampling point until the last sampling point; the data in the array Final is the z coordinate of the ordered points with the x-direction interval of xsampling and the y-direction interval of ysampling.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711153960.1A CN107886565B (en) | 2017-11-20 | 2017-11-20 | Method for ordering rock section optical scanning disordered point cloud |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711153960.1A CN107886565B (en) | 2017-11-20 | 2017-11-20 | Method for ordering rock section optical scanning disordered point cloud |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107886565A CN107886565A (en) | 2018-04-06 |
CN107886565B true CN107886565B (en) | 2020-12-15 |
Family
ID=61777974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711153960.1A Active CN107886565B (en) | 2017-11-20 | 2017-11-20 | Method for ordering rock section optical scanning disordered point cloud |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107886565B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241059A (en) * | 2018-08-30 | 2019-01-18 | 百度在线网络技术(北京)有限公司 | A kind of building method of point cloud data, device, electronic equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104463894A (en) * | 2014-12-26 | 2015-03-25 | 山东理工大学 | Overall registering method for global optimization of multi-view three-dimensional laser point clouds |
CN105844224A (en) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | Point cloud fast ordering method for on-vehicle LiDAR road points |
CN106896213A (en) * | 2017-02-22 | 2017-06-27 | 中国地质大学(武汉) | A kind of rock mass discontinuity Intelligent Recognition and information extracting method based on cloud data |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9753124B2 (en) * | 2009-07-13 | 2017-09-05 | Celartem, Inc. | LIDAR point cloud compression |
-
2017
- 2017-11-20 CN CN201711153960.1A patent/CN107886565B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104463894A (en) * | 2014-12-26 | 2015-03-25 | 山东理工大学 | Overall registering method for global optimization of multi-view three-dimensional laser point clouds |
CN105844224A (en) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | Point cloud fast ordering method for on-vehicle LiDAR road points |
CN106896213A (en) * | 2017-02-22 | 2017-06-27 | 中国地质大学(武汉) | A kind of rock mass discontinuity Intelligent Recognition and information extracting method based on cloud data |
Non-Patent Citations (2)
Title |
---|
Shear deformation and strength of the interphase between the soil–rock mixture and the benched bedrock slope surface;Duofeng Cen 等;《Acta Geotechnica》;20170430;第2114-2117页 * |
基于空间无序点云的层次构造三角化方法;成宏璟 等;《科学技术与工程》;20100331;第391-413页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107886565A (en) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107876970B (en) | Robot multilayer multi-pass welding seam three-dimensional detection and welding seam inflection point identification method | |
WO2019105044A1 (en) | Method and system for lens distortion correction and feature extraction | |
CN109060836B (en) | Machine vision-based high-pressure oil pipe joint external thread detection method | |
CN102032875B (en) | Image-processing-based cable sheath thickness measuring method | |
CN108898147B (en) | Two-dimensional image edge flattening method and device based on angular point detection | |
CN107767456A (en) | A kind of object dimensional method for reconstructing based on RGB D cameras | |
CN103632366B (en) | A kind of parameter identification method of ellipse target | |
CN110766669B (en) | Pipeline measuring method based on multi-view vision | |
CN105783711B (en) | Three-dimensional scanner correction system and correction method thereof | |
JP6899189B2 (en) | Systems and methods for efficiently scoring probes in images with a vision system | |
CN107657645B (en) | Method for calibrating parabolic catadioptric camera by using properties of conjugate diameters of straight line and circle | |
CN112381847B (en) | Pipeline end space pose measurement method and system | |
JP2013178656A (en) | Image processing device, image processing method, and image processing program | |
CN109986201B (en) | Tracking detection method and device for welding seam, storage medium and laser welding equipment | |
CN110415304B (en) | Vision calibration method and system | |
CN107886565B (en) | Method for ordering rock section optical scanning disordered point cloud | |
CN113379902A (en) | Laser scanning-based weld shape feature extraction and three-dimensional reconstruction method | |
CN115682947A (en) | Welding seam scanning and tracking processing method | |
CN116228652A (en) | Pipeline weld joint surface defect detection method based on three-dimensional point cloud curvature | |
Jin et al. | Automatic multi-stereo-vision reconstruction method of complicated tubes for industrial assembly | |
CN114972457A (en) | Precast beam digital pre-assembly matching method based on design-actual measurement point cloud model | |
Wan et al. | Extraction and segmentation method of laser stripe in linear structured light scanner | |
Tansky et al. | Multi-sensor multi-resolution data fusion modeling | |
CN103018728B (en) | Laser radar real-time imaging and building characteristic extracting method | |
CN117372498A (en) | Multi-pose bolt size measurement method based on three-dimensional point cloud |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |