CN106091972B - A kind of building change detecting method projecting dot density based on moving window - Google Patents
A kind of building change detecting method projecting dot density based on moving window Download PDFInfo
- Publication number
- CN106091972B CN106091972B CN201610512552.XA CN201610512552A CN106091972B CN 106091972 B CN106091972 B CN 106091972B CN 201610512552 A CN201610512552 A CN 201610512552A CN 106091972 B CN106091972 B CN 106091972B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- point cloud
- brick
- density
- mortar
- Prior art date
Links
- 239000008264 clouds Substances 0.000 claims abstract description 84
- 239000011449 bricks Substances 0.000 claims abstract description 55
- 239000004570 mortar (masonry) Substances 0.000 claims abstract description 32
- 238000004458 analytical methods Methods 0.000 claims abstract description 13
- 230000000875 corresponding Effects 0.000 claims abstract description 7
- 238000000034 methods Methods 0.000 claims abstract description 5
- 239000011159 matrix materials Substances 0.000 claims description 12
- 210000001503 Joints Anatomy 0.000 claims description 8
- 230000001131 transforming Effects 0.000 claims description 8
- 238000006243 chemical reactions Methods 0.000 claims description 3
- 238000000513 principal component analysis Methods 0.000 claims description 3
- 239000010950 nickel Substances 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims description 2
- 239000000284 extracts Substances 0.000 abstract description 2
- 241001269238 Data Species 0.000 abstract 2
- 280000398338 Seismic companies 0.000 description 4
- 238000010586 diagrams Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000004575 stone Substances 0.000 description 3
- 281000044456 Leica companies 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 229910052570 clay Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering processes Methods 0.000 description 1
- 239000004744 fabrics Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000011454 mudbrick Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000002002 slurries Substances 0.000 description 1
- 239000007787 solids Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical means
- G01B11/16—Measuring arrangements characterised by the use of optical means for measuring the deformation in a solid, e.g. optical strain gauge
Abstract
Description
Technical field
The present invention designs a kind of building change detecting method, and in particular to a kind of building variation based on moving window Detection method.
Background technology
In city, building is the main place of mankind's activity, and safe condition is related to mankind's daily life and warp Therefore Ji activity is changed detection to building and has a very important significance, the variation based on building is detected and repaiied It is multiple, it is research emphasis in recent years in the industry, in the hair of fabric structure especially by seismic building deformation and damage It opens up in history, masonry structure is often widely used in the foundation structure of building since ancient times, each with brick, stone, masonry and sun-dried mud brick etc. Kind of block, a kind of assembly made of being built by laying bricks or stones using mortar (mortar, clay slurry etc.) as masonry structure, masonry structure because low cost, Fire resistance, durability and construction are simply widely used, but its masonry strength is relatively low, and shock resistance is poor, therefore, to masonry The research that the building of structure is changed detection has and its important realistic meaning, the existing variation detection based on image Method, according to the level of information processing, can be divided into based on pixel, feature based and based on target variation detection, and for The variation detection of laser radar (Light Detection And Ranging) LiDAR data also gradually starts photogrammetric Studied with computer vision community, and the change detecting method precision based on image is not high enough, target is not clear enough, and easily by The influence of the quality of image, the deformation information of acquisition is also not abundant enough, can not accurately obtain in masonry structure building every piece of brick Deformation information, in addition, the low efficiency for also having seriously affected variation detection of the degree of automation.
Invention content
Goal of the invention:In view of the above-mentioned deficiencies in the prior art, it is an object of the present invention to provide a kind of based on moving window projection The building change detecting method of dot density.
Technical solution:A kind of building change detecting method projecting dot density based on moving window of the present invention, including Following steps:
(1) it uses laser scanner system to carry out the scanning of two phases to same building object, obtains building surface point cloud data, M target is set in variation architecture enclosing, wherein the observation of m >=3, laser scanner system is building surface point Three-dimensional coordinate and laser reflection intensity;
(2) target for utilizing step 1) setting, calculates two phase point cloud coordinate transformation parameter Z, is registrated to a cloud;
(3) choose fixed metope point cloud data, Dimension Reduction Analysis carried out to it using principal component analytical method, obtain feature to Measure (vi, i=1,2,3), the projection for calculating survey station point in plane where fixed metope establishes structure as coordinate origin Coordinate system is converted point cloud data to structure coordinate system by survey station coordinate system;
(4) strength information based on point cloud data classifies to the point cloud for changing metope using K mean cluster method, Isolated brick point cloud and mortar point cloud;
(5) point cloud coordinate projection to Z-direction and Y-direction is defined fixed window by the mortar point cloud for utilizing step (4) to obtain Mouth length LfixWith moving window length Lmove, by moving window, calculate separately along Z-direction and the change of Y-direction point cloud line density Change;
(6) changed according to the line density that step (5) obtains, seek horizontal and vertical cut-off rule between each brick respectively, calculated Four angular coordinates of each brick, establish brick model;
(7) each brick point cloud is obtained according to the brick model that step (6) obtains, calculates each brick center;
(8) the two phases corresponding brick central three-dimensional coordinate obtained according to step (7) obtains deformation information.
Preferably, in (1) two phase scanning process of step, target is fixed.
Preferably, step (2) and three-dimensional coordinate transfer equation described in step (3) are specific as follows:
If matrix A is the point cloud three-dimensional coordinate under A coordinate systems, matrix B is the point cloud three-dimensional coordinate under B coordinate systems, A, B two The three-dimensional coordinate transfer equation of coordinate system is as follows:
(translational movement of Δ x, Δ y and Δ z denotation coordination origins, k are scale factor, and k=0, R are A coordinate systems to B coordinates The spin matrix of system)
Preferably, specific as follows by the calculating of two phase Registration of Measuring Data conversion parameters described in step (2):
Coordinate transformation parameter Z is further write as, Z=[Δ x, Δ y, Δ z, εx, εy, εz, 1] and utilize least square method to sitting It marks conversion parameter Z and carries out parameter Estimation, the valuation that can obtain coordinate transformation parameter Z is:
Z=(ATQ-1A)-1ATQB
In formula, Q is the covariance matrix of m target coordinate measuring error under B coordinate systems, and form is as follows:
Preferably, step (3) is carried converts into structure coordinate system principal component analysis dimensionality reduction and coordinate by survey station coordinate system Origin determine specific method be:
If three-dimensional coordinate { the X of scanning element Xi=(xi,yi, zi) | i=1,2 ..., n }, construct corresponding covariance matrix:
Wherein, For the barycentric coodinates of point set, to Matrix C carry out it is main at Analysis, can acquire three eigenvalue λs1、λ2、λ3It arranges in descending order, obtains λ1≥λ2> λ3> 0, λ3Corresponding feature vector v3, and v3For normal vector, v3For unit vector of the X-axis under survey station coordinate system of structure coordinate system, and the Z axis of structure coordinate system It to be directed toward consistent with survey station coordinate system, Y-axis constitutes right-handed coordinate system perpendicular to determining XOZ planes, calculating survey station point coordinates S (0, 0,0) in the projection S'(x of plane where fixed metopes,ys,zs), as the coordinate origin of structure coordinate system, therefore translate Vector (Δ x, Δ y, Δ z)=(- xs,-ys,-zs), after establishing translation parameters and reference axis rotation parameter, after two phases were registrated Point cloud data rotate to structure coordinate system;
Preferably, step (4) carry the K mean cluster method based on strength information separation metope brick and mortar it is specific Method is:
Using cluster sum of squared errors function E as clustering criteria function, using strength information as categorical attribute, In,xijIt is j-th of sample of the i-th class, miIt is the cluster centre or barycenter of the i-th class, niIt is the i-th class Number of samples, K mean cluster algorithm find k best cluster centres, wherein k=2, by all n samples by iterating This point is assigned to the cluster centre nearest from it so that cluster error sum of squares E is minimum, and process is as follows:
A is randomly assigned k cluster centre mi(i=1,2 ..., k);
B, to each sample xiThe cluster centre nearest from it is found, such is assigned it to;
C recalculates each brand new center:NiIt is the current sample number of the i-th cluster;
D calculates deviation,
E returns to m if E values restraini(i=1,2 ..., k), algorithm terminate, and otherwise return to b;
Preferably, step (5) institute Tilly mortar point cloud calculates point cloud line density variation, tool based on window Mobile Method Body method is as follows:
If the mean breadth of mortar is known quantity Lmortar, define stationary window length LfixWith moving window length Lmove, Three meets following relationship:
Lmortar≈Lwindow+2Lmove
It calculates separately along Z-direction and Y-direction moving window number:
Wherein [] is rounding symbol,
Calculate separately counting out in Z-direction and each window of Y-direction:
nzi(i=1,2 ..., ny),nyi(i=1,2 ..., nz)
Calculate separately the line density of the point along Z-direction and Y-direction:
Density_z=(nz(i-1)+nzi+nz(i+1))/(3Lfix) (i=2,3 ..., (nz-1))
Density_y=(ny(i-1)+nyi+ny(i+1))/(3Lfix) (i=2,3 ..., (ny-1))
For Z-direction and each window of Y-direction, its line density change rate is calculated:
Grad (i, 1)=Density_y (i)-Density_y (i-1)
Grad (i, 2)=Density_y (i+1)-Density_y (i)
Preferably, step (6) is described is changed using line density, calculates horizontal and vertical cut-off rule between each brick, specific Method is as follows:
For the Z-direction or Y-direction of analysis, point cloud line density is apparently higher than non-at the mortar joints being perpendicularly to the direction Seam crossing region, therefore, select in this direction window averag density as threshold value (ntotalFor point Cloud total number), it is residing solid when being more than threshold value when putting cloud density in window, and meeting Grad (i, 1) > 0 and Grad (i, 2) < 0 Determine the range that window is mortar joints, by the point cloud computing mortar joints center line within the scope of mortar joints, then passes through brick Block girth seam center line computation brick four angular coordinates of model, establish brick model.
Advantageous effect:A kind of building change detecting method projecting dot density based on moving window of the present invention, a side On the other hand face high degree of automation takes full advantage of initial data, precision is high, under the premise of accuracy guarantee, it is ensured that brick The deformation information of each brick of stone structure building can be obtained effectively.
Description of the drawings
Fig. 1 is the flow chart of the present invention;
Fig. 2 is the survey station coordinate system schematic diagram of the present invention;
Fig. 3 is the structure coordinate system schematic diagram of the present invention;
Fig. 4 is the brick model schematic of the present invention;
Fig. 5 is a certain metope K mean value classification results that the present invention chooses;
Fig. 6 is that present invention point cloud subpoint line density changes histogram;
Fig. 7 is each brick central point schematic diagram that the present invention extracts;
The brick 3 D deformation information that the positions Fig. 8 present invention obtains.
Specific implementation mode
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.
As shown in Figure 1, a kind of building change detecting method projecting dot density based on moving window, including following step Suddenly:
(1) it uses laser scanner system to carry out the scanning of two phases to same building object, obtains building surface point cloud data, M target is arranged in fixed area around variation building, under normal circumstances m >=4, it is contemplated that the included automatic peace of laser scanner Flat function, and instrument is horizontality when each measurement, therefore, m >=3, the observation of laser scanner system is to build Build the three-dimensional coordinate and laser reflection intensity of object surface point;
(2) target for utilizing step (1) setting, calculates two phase point cloud coordinate transformation parameter Z, is registrated to a cloud;
(3) choose fixed metope point cloud data, Dimension Reduction Analysis carried out to it using principal component analytical method, obtain feature to Measure (vi, i=1,2,3), the projection for calculating survey station point in plane where fixed metope establishes structure as coordinate origin O Coordinate system is converted point cloud data to structure coordinate system by survey station coordinate system;
(4) strength information based on point cloud data classifies to the point cloud for changing metope using K mean cluster method, Isolated brick point cloud and mortar point cloud;
(5) point cloud coordinate projection to Z-direction and Y-direction is defined fixed window by the mortar point cloud for utilizing step (4) to obtain Mouth length LfixWith moving window length Lmove, by moving window, calculate separately along Z-direction and the change of Y-direction point cloud line density Change;
(6) changed according to the line density that step (5) obtains, seek horizontal and vertical cut-off rule between each brick respectively, calculated Four angular coordinates of each brick, establish brick model;
(7) each brick point cloud is obtained according to the brick model that step (6) obtains, calculates each brick center;
(8) the two phases corresponding brick central three-dimensional coordinate obtained according to step (7) obtains deformation information.
By taking " certain field experiment masonry structure building changes in BEFORE AND AFTER EARTHQUAKE to be detected " as an example, the present invention is further elaborated:
{ 1 } building is scanned using Leica C10 laser scanner systems, instrument is set up in position as shown in Figure 5 Device, 4 targets are laid in stability region other than experimental architecture object, are scanned before and after seismic test respectively, and two phases of acquisition build Object surface laser point cloud data is built, observation includes two classes:Three-dimensional coordinate, laser reflection intensity;
{ 2 } as shown in Fig. 2, centered on survey station, vertical direction is Z axis, and the target being arranged using step { 1 } calculates two Phase point cloud coordinate transformation parameter Z is put before putting cloud coordinate registration to seismic test after seismic test in cloud survey station coordinate system, until This, two phase point clouds are under identical coordinate system, in order to which the later stage carries out deformation analysis;
{ 3 } as shown in figure 3, three reference axis are respectively parallel to adjacent three face of the building centered on certain angle point, choosing Fixed metope point cloud data is taken, Dimension Reduction Analysis is carried out using principal component analysis, obtains fixed metope normal vector, i.e. structure coordinate system X-axis unit vector is [0.9348 0.3551 0.0011], and structure coordinate system Z axis unit vector is [0 0 1], therefore Y-axis list Bit vector is [0.3551-0.9348 0], and the subpoint of survey station coordinate system coordinate origin to fixed metope is [8.2766 3.1442 0.0095], according to above-mentioned parameter and transformational relation, point cloud data is converted by survey station coordinate system to structure coordinate system;
{ 4 } as shown in figure 5, according to the strength information of point cloud data, using K mean cluster method respectively to two phase metopes Point cloud is classified, isolated brick point cloud and mortar point cloud;
{ 5 } pass through site inspection and measure, the mean breadth of mortar is 10-12mm between brick, therefore, according to Lmortar≈ Lwindow+2LmoveDefine stationary window length Lfix=0.008m and moving window length Lmove=0.002m is obtained using step { 4 } The mortar point cloud arrived, by moving window, is calculated separately along Z-direction and Y-direction by point cloud coordinate projection to Z-direction and Y-direction Point cloud line density variation, the results are shown in Figure 6;
{ 6 } changed according to the line density that step { 5 } obtains, seek the vertical and horizontal cut-off rule between each brick respectively, and Four angular coordinates for calculating each brick, establish brick model as shown in Figure 4;
{ 7 } it as shown in fig. 7, obtaining each brick point cloud data according to the brick model that step { 6 } obtains, calculates in each brick The heart;
{ 8 } as shown in figure 8, the two phases corresponding brick central three-dimensional coordinate obtained according to step { 7 }, obtains deformation letter Breath.
Claims (7)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512552.XA CN106091972B (en) | 2016-06-30 | 2016-06-30 | A kind of building change detecting method projecting dot density based on moving window |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512552.XA CN106091972B (en) | 2016-06-30 | 2016-06-30 | A kind of building change detecting method projecting dot density based on moving window |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106091972A CN106091972A (en) | 2016-11-09 |
CN106091972B true CN106091972B (en) | 2018-09-21 |
Family
ID=57212867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610512552.XA CN106091972B (en) | 2016-06-30 | 2016-06-30 | A kind of building change detecting method projecting dot density based on moving window |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106091972B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416785B (en) * | 2018-03-26 | 2020-08-11 | 北京进化者机器人科技有限公司 | Topology segmentation method and device for closed space |
CN109035206B (en) * | 2018-06-28 | 2020-10-09 | 中国地震局地震预测研究所 | Building damage state detection method |
CN109684970B (en) * | 2018-12-18 | 2020-08-07 | 暨南大学 | Window length determination method for moving principal component analysis of structural dynamic response |
CN110132233A (en) * | 2019-04-16 | 2019-08-16 | 西安长庆科技工程有限责任公司 | Current relief map drawing practice under a kind of CASS environment based on point cloud data |
CN110060338B (en) * | 2019-04-25 | 2020-11-10 | 重庆大学 | Prefabricated part point cloud identification method based on BIM model |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103307999A (en) * | 2013-06-14 | 2013-09-18 | 河海大学 | Three-dimensional laser scanning control rack and field operation scanning and point cloud registration method for same |
CN103870845A (en) * | 2014-04-08 | 2014-06-18 | 重庆理工大学 | Novel K value optimization method in point cloud clustering denoising process |
CN103940356A (en) * | 2014-02-27 | 2014-07-23 | 山东交通学院 | Building overall-deformation monitoring method based on three-dimensional laser scanning technology |
CN105136054A (en) * | 2015-04-27 | 2015-12-09 | 北京工业大学 | Fine structure deformation monitoring method and system based on ground three-dimensional laser scanning |
CN105423915A (en) * | 2015-11-16 | 2016-03-23 | 天津师范大学 | Accurate positioning method of planar target for ground laser scanning data registration |
-
2016
- 2016-06-30 CN CN201610512552.XA patent/CN106091972B/en active IP Right Grant
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103307999A (en) * | 2013-06-14 | 2013-09-18 | 河海大学 | Three-dimensional laser scanning control rack and field operation scanning and point cloud registration method for same |
CN103940356A (en) * | 2014-02-27 | 2014-07-23 | 山东交通学院 | Building overall-deformation monitoring method based on three-dimensional laser scanning technology |
CN103870845A (en) * | 2014-04-08 | 2014-06-18 | 重庆理工大学 | Novel K value optimization method in point cloud clustering denoising process |
CN105136054A (en) * | 2015-04-27 | 2015-12-09 | 北京工业大学 | Fine structure deformation monitoring method and system based on ground three-dimensional laser scanning |
CN105423915A (en) * | 2015-11-16 | 2016-03-23 | 天津师范大学 | Accurate positioning method of planar target for ground laser scanning data registration |
Also Published As
Publication number | Publication date |
---|---|
CN106091972A (en) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106997049B (en) | Method and device for detecting barrier based on laser point cloud data | |
CN103148804B (en) | Indoor unknown structure identification method based on laser scanning | |
CN104931977B (en) | A kind of obstacle recognition method for intelligent vehicle | |
CN104062973B (en) | A kind of mobile robot based on logos thing identification SLAM method | |
Girardeau-Montaut | CloudCompare | |
CN105180890B (en) | Merge the ORIENTATION OF DISCONTINUITY IN ROCK MASS measuring method of laser point cloud and digitized video | |
CN100561121C (en) | Measurement method of planeness based on Flame Image Process and image recognition | |
Abellán et al. | Terrestrial laser scanning of rock slope instabilities | |
Biosca et al. | Unsupervised robust planar segmentation of terrestrial laser scanner point clouds based on fuzzy clustering methods | |
CN103268729B (en) | Based on mobile robot's tandem type map creating method of composite character | |
González-Jorge et al. | Monitoring biological crusts in civil engineering structures using intensity data from terrestrial laser scanners | |
CN105352504B (en) | The coal mining machine positioning device and method that a kind of inertial navigation merges with laser scanning | |
CN103020342B (en) | Method for extracting contour and corner of building from ground LiDAR data | |
CN101539629B (en) | Remote sensing image change detection method based on multi-feature evidence integration and structure similarity | |
CN103487034B (en) | Method for measuring distance and height by vehicle-mounted monocular camera based on vertical type target | |
CN106053475B (en) | Tunnel defect tunneling boring dynamic device for fast detecting based on active panoramic vision | |
CN103941264B (en) | Positioning method using laser radar in indoor unknown environment | |
CN105136054B (en) | The fine deformation monitoring method of structures and system based on Three Dimensional Ground laser scanning | |
CN105910601B (en) | A kind of indoor ground magnetic positioning method based on Hidden Markov Model | |
CN103940356B (en) | A kind of building bulk deformation monitoring method based on three-dimensional laser scanning technique | |
JP2017072422A (en) | Information processing device, control method, program, and storage medium | |
CN104807460B (en) | Unmanned plane indoor orientation method and system | |
CN102128625B (en) | Initial matching method for use in gravimetric map matching in gravity-aided inertial navigation system | |
CN103488972B (en) | Fingertip Detection based on depth information | |
CN106960468B (en) | A kind of 3 D laser scanning point cloud Accuracy Assessment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C10 | Entry into substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20161109 Assignee: ANHUI JINHE SOFTWARE Co.,Ltd. Assignor: HOHAI University Contract record no.: X2020980010569 Denomination of invention: Building change detection method based on projection point density of moving window Granted publication date: 20180921 License type: Common License Record date: 20210104 |