CN107784690A - A kind of completion plan automatic drafting method based on three-dimensional point cloud - Google Patents
A kind of completion plan automatic drafting method based on three-dimensional point cloud Download PDFInfo
- Publication number
- CN107784690A CN107784690A CN201711117305.0A CN201711117305A CN107784690A CN 107784690 A CN107784690 A CN 107784690A CN 201711117305 A CN201711117305 A CN 201711117305A CN 107784690 A CN107784690 A CN 107784690A
- Authority
- CN
- China
- Prior art keywords
- fitting
- data
- point cloud
- dimensional point
- cloud
- 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.)
- Granted
Links
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
Abstract
The invention discloses a kind of completion plan automatic drafting method based on three-dimensional point cloud, the height of the architectural plan automatic drafting method energy plan floor needed for provided has both sides advantage by drawing completion plan to three dimensional point cloud processing is automatic：On the one hand avoid carrying out feature extraction to building three-dimensional point cloud, linear fit directly is carried out to two-dimensional points cloud using least square fitting, work difficulty is reduced, improves pattern precision；Another aspect automatic Fitting simultaneously draws plan and substantially increases operating efficiency.
Description
Technical field
It is specially a kind of based on three-dimensional point cloud the present invention relates to Surveying and mapping or engineering survey technical field of data processing
Completion plan automatic drafting method.
Background technology
With the rapid development of society, especially scientific and technical continuous progress, there occurs earth-shaking for the life of people
Change.The continuous quickening of supporter's urbanization, a large amount of populations pour in city so that huge people is born on limited soil
The clothing, food, lodging and transportion -- basic necessities of life pressure of mouth, to social public service, higher want just is proposed especially as the building of underlying hardware thing for this
Ask, therefore, countless overpass, skyscraper are rised sheer from level ground like the mushrooms after rain.And for the building of completion, how to ensure
It is that high standard builds up thing, then what acceptance survey just showed is even more important.
Acceptance survey is an important step of urban architecture engineering, and it is after architectural engineering completion, to building reality
The comprehensive measurement work that body is carried out, including draw building plan and calculate the contents such as construction area.Traditional completion
Measuring method is less efficient, and can not meet the requirement of city acceptance survey gradually.And three-dimensional laser scanning technique has density
High, the features such as precision is high, can three-dimensional coordinate data that is quick, accurate, obtaining to large area, high-resolution measurand surface,
Conventional measurement techniques, which are broken through, for mapping worker provides a kind of brand-new data acquisition means.
Application of the 3 D laser scanning in the acceptance survey of city at present is less, and technology is not mature enough, and survey crew passes through
Three-dimensional point cloud is generally required first to cut three-dimensional model building when drawing completion plan, and then manual drawing floor is put down
The problems such as manual intervention is more, and fit line precision is not high be present in face figure, the method.
The content of the invention
In order to overcome above-mentioned deficiency of the prior art, the invention provides a kind of automaticity height, fitting precision are good
The completion plan automatic drafting method based on three-dimensional point cloud.
The object of the present invention is achieved like this：
A kind of completion plan automatic drafting method based on three-dimensional point cloud, comprises the following steps：
The first step, story height information is obtained in spliced three-dimensional point cloud model, it is determined that needing to draw plan
Story height H；
Second step, the three dimensional point cloud of extraction is projected in XOY plane, obtain the cloud data of two dimension；
3rd step, makes m=1 first, carries out linear least square fitting to two-dimentional scatterplot data, is determined according to error requirement
Error of fitting E and threshold value e, every time fitting are promoted in units of q data point, obtain fit equation f (x), digital simulation
By mistake
Difference
4th step, compare this error of fitting E and threshold value e size, if error of fitting meets condition E<E, then with starting
O'clock it is that section carries out least square fitting to the 2q data point, repeat step three is fitted；Until (n+1) secondary fitting
Error of fitting E>During e, this section of fitting terminates.
5th step, this section is participated in nq data composition data set S of fitting, and the fit equation of minute book segment data is
F (x), data interval are starting point to n-th q data point；
6th step, the 3rd step to the 5th step is repeated, until the fitting of all scatterplot data terminates, include m fit equation and m
Individual fit interval；
7th step, after being fitted end to data, Automatic Drawing is carried out according to fit equation and fit interval.It is final complete
Into the automatic drafting of completion plan.
In the first step, the point cloud of corresponding height is extracted in cloud data, if the height h of scanning element, i.e. scanning element
Z values in coordinate, then h should meet | h-H |<ξ, 0.05m is taken according to a cloud distribution density ξ；
In the second step, cloud data only includes X-coordinate and Y-coordinate.
Built for general rule, this method uses least square linear fit, and wherein fitting function is once linear
Function f (x)=kx+b.
Positive beneficial effect：Architectural plan automatic drafting method provided by the present invention can plan floor needed for basis
Height by drawing completion plan to three dimensional point cloud processing is automatic, have of both advantage：On the one hand avoid to building
Build thing three-dimensional point cloud and carry out feature extraction, linear fit directly is carried out to two-dimensional points cloud using least square fitting, reduces work
Make difficulty, improve pattern precision；Another aspect automatic Fitting simultaneously draws plan and substantially increases operating efficiency.
Brief description of the drawings
Fig. 1 is the flow chart that the present invention is implemented；
Fig. 2 is intercept point cloud top view in Recap of the present invention；
Fig. 3 is intercept point cloud front view in Recap of the present invention；
Fig. 4 is the fit Plane figure of the present invention.
Embodiment
With reference to specific embodiment, the present invention is described further：
As shown in figure 1, a kind of completion plan automatic drafting method based on three-dimensional point cloud, comprises the following steps：
The first step, story height information is obtained in spliced three-dimensional point cloud model, it is determined that needing to draw plan
Story height H；
Second step, the three dimensional point cloud of extraction is projected in XOY plane, obtain the cloud data of two dimension；
3rd step, makes m=1 first, carries out linear least square fitting to two-dimentional scatterplot data, is determined according to error requirement
Error of fitting E and threshold value e, every time fitting are promoted in units of q data point, obtain fit equation f (x), digital simulation
Error
4th step, compare this error of fitting E and threshold value e size, if error of fitting meets condition E<E, then with starting
O'clock it is that section carries out least square fitting to the 2q data point, repeat step three is fitted；Until (n+1) secondary fitting
Error of fitting E>During e, this section of fitting terminates.
5th step, this section is participated in nq data composition data set S of fitting, and the fit equation of minute book segment data is
F (x), data interval are starting point to n-th q data point；
6th step, the 3rd step to the 5th step is repeated, until the fitting of all scatterplot data terminates, include m fit equation and m
Individual fit interval；
7th step, after being fitted end to data, Automatic Drawing is carried out according to fit equation and fit interval.It is final complete
Into the automatic drafting of completion plan.
In the first step, the point cloud of corresponding height is extracted in cloud data, if the height h of scanning element, i.e. scanning element
Z values in coordinate, then h should meet | h-H |<ξ, 0.05m is taken according to a cloud distribution density ξ；
In the second step, cloud data only includes X-coordinate and Y-coordinate.
Built for general rule, this method uses least square linear fit, and wherein fitting function is once linear
Function f (x)=kx+b.
Embodiment
Method flow according to Fig. 1, it is right using floor where drawing 20 meters of certain commodity building threedimensional model as application example
The present invention furthers elucidate：
The first step, will splice, the commercial house three-dimensional point cloud model after denoising saves as pts forms, make H=20.000m, ξ
0.050m is taken, the cloud data that height h in cloud data is 19.950m~20.050m, cloud data are extracted by matlab
The form of expression is as shown in Figures 2 and 3 in Recap softwares.
Second step, read point cloud file and preserve, establish the matrix of one two row and be stored in a cloud first two columns data
In matrix.Now matrix is only comprising a cloud X-coordinate and Y-coordinate.
3rd step, make m=1, q=30, e (u)=0.1, the least square fitting of the 1st section of straight line of progress, i.e., from first
Data point to the 30th data point is fitted, and obtaining fit equation is：f1(x)=k1x+b1, wherein k1=-
2.427272727 b1=467160.245, digital simulation error calculation error of fitting
4th step, compare error of fitting E1With threshold value 0.1, if error of fitting meets condition E1<0.1, then with starting point to
60 data points carry out least square fitting, obtain fit equation：f1 (2)(x)=k1 (2)x+b1 (2), digital simulation error E2.Directly
To the error of fitting E of (n+1) secondary fitting(n+1)>When 0.1, this section of fitting terminates.
5th step, this section is participated in 30n data composition data set S of fitting1, the fit equation of minute book segment data
For f1 (n)(x)=k1 (n)x+b1 (n), data interval is starting point to the 30n data point.
6th step, the 3rd step to the 5th step is repeated, carry out second segment fitting a straight line, obtain fit equation f2(x)=k2x+
b2, and corresponding data section.Until the fitting of all scatterplot data terminates, and comprising m fit equation and m fit interval, every
Line segment equation coefficient is as shown in table 1.
The every line segment equation coefficient of table 1
7th step, after being fitted end to data, Automatic Drawing, such as Fig. 4 are carried out according to fitting function and fit interval
It is shown.It is finally completed the automatic drafting of completion plan.
It is right that architectural plan automatic drafting method energy height of plan floor needed for provided by the present invention passes through
Three dimensional point cloud processing is automatic to draw completion plan, there is both sides advantage：On the one hand avoid to building three-dimensional point cloud
Feature extraction is carried out, linear fit directly is carried out to two-dimensional points cloud using least square fitting, work difficulty is reduced, improves
Pattern precision；Another aspect automatic Fitting simultaneously draws plan and substantially increases operating efficiency.
Above case study on implementation is merely to illustrate the preferred embodiment of the present invention, but the present invention is not limited to above-mentioned embodiment party
Formula, in the field those of ordinary skill possessed knowledge, that is made within the spirit and principles in the present invention is any
Modification, equivalent substitute and improvement etc., it is regarded as the protection domain of the application.
Claims (4)
1. a kind of completion plan automatic drafting method based on three-dimensional point cloud, it is characterised in that comprise the following steps：
The first step, story height information is obtained in spliced three-dimensional point cloud model, it is determined that needing to draw the floor of plan
Height H；
Second step, the three dimensional point cloud of extraction is projected in XOY plane, obtain the cloud data of two dimension；
3rd step, makes m=1 first, and linear least square fitting is carried out to two-dimentional scatterplot data, determines to be fitted according to error requirement
Error E and threshold value e, every time fitting are promoted in units of q data point, obtain fit equation f (x), digital simulation error
4th step, compare this error of fitting E and threshold value e size, if error of fitting meets condition E<E, then with starting point extremely
The 2q data point is that section carries out least square fitting, and repeat step three is fitted；Until the fitting of (n+1) secondary fitting
Error E>During e, this section of fitting terminates；
5th step, this section is participated in nq data composition data set S of fitting, the fit equation of minute book segment data is f
(x), data interval is starting point to n-th q data point；
6th step, the 3rd step to the 5th step is repeated, until the fitting of all scatterplot data terminates, intended comprising m fit equation and m
Close section；
7th step, after being fitted end to data, Automatic Drawing is carried out according to fit equation and fit interval；It is finally completed and completes
The automatic drafting of work plan.
A kind of 2. completion plan automatic drafting method based on three-dimensional point cloud according to claim 1, it is characterised in that：
In the first step, the point cloud of corresponding height is extracted in cloud data, if the height h of scanning element, that is, scans the Z in point coordinates
Value, then h should meet | h-H |<ξ, 0.05m is taken according to a cloud distribution density ξ.
A kind of 3. completion plan automatic drafting method based on three-dimensional point cloud according to claim 1, it is characterised in that：
In the second step, cloud data only includes X-coordinate and Y-coordinate.
A kind of 4. completion plan automatic drafting method based on three-dimensional point cloud according to claim 1, it is characterised in that：
Built for general rule, this method uses least square linear fit, and wherein fitting function is once linear function f (x)
=kx+b.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711117305.0A CN107784690B (en) | 2017-11-13 | 2017-11-13 | Completion plan automatic drawing method based on three-dimensional point cloud |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711117305.0A CN107784690B (en) | 2017-11-13 | 2017-11-13 | Completion plan automatic drawing method based on three-dimensional point cloud |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107784690A true CN107784690A (en) | 2018-03-09 |
| CN107784690B CN107784690B (en) | 2021-04-13 |
Family
ID=61432828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711117305.0A Active CN107784690B (en) | 2017-11-13 | 2017-11-13 | Completion plan automatic drawing method based on three-dimensional point cloud |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107784690B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108846175A (en) * | 2018-05-30 | 2018-11-20 | 链家网（北京）科技有限公司 | A kind of vector house type drawing generating method and device |
| CN110599569A (en) * | 2019-09-16 | 2019-12-20 | 上海市刑事科学技术研究院 | Method for generating two-dimensional plane graph inside building, storage device and terminal |
| CN113610909A (en) * | 2021-09-07 | 2021-11-05 | 中国电建集团昆明勘测设计研究院有限公司 | Point cloud profile generation system and method based on distance search |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1885349A (en) * | 2006-07-05 | 2006-12-27 | 东南大学 | Point cloud hole repairing method for three-dimensional scanning |
| US7386582B2 (en) * | 2003-09-29 | 2008-06-10 | International Business Machines Corporation | Method and structure for producing high performance linear algebra routines using a hybrid full-packed storage format |
| CN104036544A (en) * | 2014-06-25 | 2014-09-10 | 西安煤航信息产业有限公司 | Building roof reconstruction method based on airborne LiDAR data |
| CN104658039A (en) * | 2015-02-12 | 2015-05-27 | 南京市测绘勘察研究院有限公司 | Urban digital map three-dimensional modeling manufacturing method |
-
2017
- 2017-11-13 CN CN201711117305.0A patent/CN107784690B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7386582B2 (en) * | 2003-09-29 | 2008-06-10 | International Business Machines Corporation | Method and structure for producing high performance linear algebra routines using a hybrid full-packed storage format |
| CN1885349A (en) * | 2006-07-05 | 2006-12-27 | 东南大学 | Point cloud hole repairing method for three-dimensional scanning |
| CN104036544A (en) * | 2014-06-25 | 2014-09-10 | 西安煤航信息产业有限公司 | Building roof reconstruction method based on airborne LiDAR data |
| CN104658039A (en) * | 2015-02-12 | 2015-05-27 | 南京市测绘勘察研究院有限公司 | Urban digital map three-dimensional modeling manufacturing method |
Non-Patent Citations (1)
| Title |
|---|
| CHEN, JY 等: "POINT CLOUD MODELING USING ALGEBRAIC TEMPLATE", 《INTERNATIONAL JOURNAL OF INNOVATIVE COMPUTING INFORMATION AND CONTROL》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108846175A (en) * | 2018-05-30 | 2018-11-20 | 链家网（北京）科技有限公司 | A kind of vector house type drawing generating method and device |
| CN110599569A (en) * | 2019-09-16 | 2019-12-20 | 上海市刑事科学技术研究院 | Method for generating two-dimensional plane graph inside building, storage device and terminal |
| CN110599569B (en) * | 2019-09-16 | 2023-09-15 | 上海市刑事科学技术研究院 | Method for generating two-dimensional plan inside building, storage device and terminal |
| CN113610909A (en) * | 2021-09-07 | 2021-11-05 | 中国电建集团昆明勘测设计研究院有限公司 | Point cloud profile generation system and method based on distance search |
| CN113610909B (en) * | 2021-09-07 | 2023-10-20 | 中国电建集团昆明勘测设计研究院有限公司 | Point cloud profile generation system and method based on distance search |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107784690B (en) | 2021-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103020342B (en) | Method for extracting contour and corner of building from ground LiDAR data | |
| Luo et al. | Wind tunnel simulation of the three-dimensional airflow patterns behind cuboid obstacles at different angles of wind incidence, and their significance for the formation of sand shadows | |
| CN110276732B (en) | Mountain area point cloud cavity repairing method considering topographic characteristic line elements | |
| CN106970375A (en) | A kind of method that building information is automatically extracted in airborne laser radar point cloud | |
| CN107784690A (en) | A kind of completion plan automatic drafting method based on three-dimensional point cloud | |
| CN104050474A (en) | Method for automatically extracting island shoreline based on LiDAR data | |
| CN105066912B (en) | The step-length scaling method of rock beam surface scan data in acid etch physical simulation experiment | |
| CN105760581A (en) | Channel drainage basin renovation planning simulating method and system based on OSG | |
| CN108197283A (en) | The method that flood two, three-dimensional dynamic exhibition are carried out using river measured water level | |
| CN103256914B (en) | A kind of method and system calculating silt arrester inundated area based on DEM | |
| KR100657870B1 (en) | Method for sampling the ground height using aviation laser measurement data | |
| CN102353678A (en) | Method for measuring cultural relic diseases | |
| CN108038881A (en) | A kind of computational methods using directional light evaluation rock structural face fluctuating quantity | |
| KR100309529B1 (en) | Realtime shade simulation | |
| CN104299262A (en) | Three-dimensional cloud simulating method based on speed field flow line | |
| CN111369604A (en) | Grid-based topographic feature point extraction method and processing terminal | |
| CN114511995B (en) | Flood classification early warning method based on expressed model | |
| CN116188734A (en) | AR-based house graphic mapping method | |
| CN106846481B (en) | Geological profile generation method | |
| CN115131571A (en) | Building local feature point identification method based on six fields of point cloud pretreatment | |
| CN106846486A (en) | Ancient and well-known trees mapping method and device based on a cloud | |
| CN103218853A (en) | Crop single-root deformable modeling method | |
| KR20070070101A (en) | A contour input device and method thereof | |
| CN105182441A (en) | Three-dimensional projection method for surveyed stratigraphic section drawing | |
| Li et al. | A hierarchical contour method for automatic 3D city reconstruction from LiDAR data |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Nianqing Inventor after: Guo Jiangning Inventor after: Wang Tian Inventor after: Hu Chunxia Inventor before: Guo Jiangning Inventor before: Chen Nianqing Inventor before: Wang Tian Inventor before: Hu Chunxia |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |