CN112581579B - Method for extracting point cloud data of magnetic suspension sliding surface - Google Patents
Method for extracting point cloud data of magnetic suspension sliding surface Download PDFInfo
- Publication number
- CN112581579B CN112581579B CN202011535988.3A CN202011535988A CN112581579B CN 112581579 B CN112581579 B CN 112581579B CN 202011535988 A CN202011535988 A CN 202011535988A CN 112581579 B CN112581579 B CN 112581579B
- Authority
- CN
- China
- Prior art keywords
- cloud data
- point cloud
- sliding surface
- magnetic suspension
- coordinate system
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Control Of Vehicles With Linear Motors And Vehicles That Are Magnetically Levitated (AREA)
Abstract
The invention relates to a method for extracting point cloud data of a magnetic suspension sliding surface, which comprises the following steps: acquiring point cloud data; converting the point cloud data from a space coordinate system to an orbit coordinate system to obtain a fitting range of a sliding surface; constructing a solving model to obtain an initial plane; calculating the distance from all point cloud data in the fitting range to the initial plane, and removing the point cloud data beyond the distance according to a first set threshold value to obtain preprocessed point cloud data; calculating the mean value of the preprocessed point cloud data, calculating the residual error and the standard deviation of the preprocessed point cloud data, and removing the point cloud data of which the residual error is greater than a second set threshold or the standard deviation is greater than a third threshold to obtain single-segment point cloud data; and fitting the single-section point cloud data to obtain integral point cloud data of the magnetic suspension sliding surface. Compared with the prior art, the automatic extraction of the point cloud data of the magnetic suspension sliding surface is realized, the final noise of the point cloud data of the magnetic suspension sliding surface can be effectively reduced, and the defect that the fitting range cannot be matched with point cloud due to the complex environment of the magnetic suspension track is overcome.
Description
Technical Field
The invention relates to the field of magnetic suspension point cloud data extraction, in particular to a magnetic suspension sliding surface point cloud data extraction method.
Background
Magnetic levitation traffic has become an important development direction for rail traffic with the advantages of low noise, fast running speed, low maintenance cost, small running vibration, strong resistance to severe environment and the like. The magnetic suspension sliding surface is a functional surface positioned on two sides of the track, and has an important function for guaranteeing the stable operation of the magnetic suspension train. Due to the high-speed running of the maglev train, direct contact exists in a local area, and a maglev sliding surface can generate certain nonlinear deformation. Therefore, the method for detecting the magnetic suspension sliding surface regularly has extremely important significance for the safe operation of the magnetic suspension train. The magnetic suspension track is measured by adopting a mobile laser scanning mode, a large amount of three-dimensional point cloud information of the magnetic suspension track can be quickly and accurately acquired, and the point cloud information comprises a sliding surface, a guide surface, a stator surface and point cloud information of various structures. The problem of how to extract the point cloud of the sliding surface from the point cloud information of the sliding surface, the guide surface, the stator surface and various structures needs to be solved.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a method for extracting point cloud data of a magnetic suspension sliding surface.
The purpose of the invention can be realized by the following technical scheme:
a method for extracting point cloud data of a magnetic suspension sliding surface comprises the following steps:
s1: acquiring point cloud data;
s2: converting the point cloud data from a space coordinate system to a track coordinate system, and obtaining a fitting range of a sliding surface according to design parameters of the magnetic suspension sliding surface;
s3: constructing a solving model according to the difference between a function value and an actual value of a space plane equation constructed by point cloud data in an orbit coordinate system to obtain an initial plane;
s4: calculating the distance from all point cloud data in the fitting range to the initial plane, and removing unreasonable point cloud data according to a first set threshold value to obtain preprocessed point cloud data;
s5: calculating the mean value of the preprocessed point cloud data, taking the mean value as a standard value, calculating the residual error and the standard deviation of the preprocessed point cloud data, and removing the point cloud data of which the residual error is greater than a second set threshold value or the standard deviation is greater than a third threshold value to obtain single-segment point cloud data;
s6: and fitting the single-section point cloud data to obtain the point cloud data of the magnetic suspension sliding surface.
In S1, the point cloud data is obtained through a laser scanner.
In S2, the track coordinate system comprises an axis I, an axis q and an axis h, wherein the axis I represents mileage, the axis q represents offset from a point to a center line, the axis h represents elevation, and an origin of the track coordinate system is a central point of an initial position.
And S3, normally distributing the elevations of the point cloud data in the fitting range, removing the point cloud data with lower confidence coefficient according to the elevations, and constructing a space plane equation by using the rest point cloud data.
In S3, the spatial plane equation is:
p(l,q,h)=al+bq+ch+m=0
wherein, a, b and c are normal vectors of a plane equation, l, q and h are mileage, offset and elevation respectively, and m is the distance from an origin to a space plane.
In S3, the expression of the solution model is:
wherein e is the difference between the function value and the actual value, and i represents the number of the point cloud data.
In S4, the first threshold is the prior distance t from the point cloud data to the initial plane.
In S5, the second threshold is a variance σ of the residual error.
The third threshold is three times the standard deviation.
And S6, fitting the single-segment point cloud data along the axis I of the track coordinate system.
Compared with the prior art, the invention has the following advantages:
(1) Based on the acquired point cloud data, the automatic extraction of the point cloud data of the magnetic suspension sliding surface can be quickly and efficiently realized; and abnormal point cloud data are removed for multiple times, so that the final noise of the point cloud data of the magnetic suspension sliding surface can be effectively reduced.
(2) The sliding surface is extracted by a sectional fitting mode, and the defect that the fitting range cannot be matched with point clouds due to the complex environment of the magnetic suspension track is overcome.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic diagram of point cloud data acquired by a laser scanner;
FIG. 3 is a schematic diagram of the orbital coordinate system of the present invention;
FIG. 4 is a comparison of the initial plane and a single segment point cloud data plane.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. The present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation manner and a specific operation process are given, but the scope of the present invention is not limited to the following embodiments.
Examples
The present embodiment provides a method for extracting point cloud data of a magnetic levitation sliding surface, as shown in fig. 1, the method includes the following steps:
s1: acquiring point cloud data;
s2: converting the point cloud data from a space coordinate system to a track coordinate system, and obtaining a fitting range of a sliding surface according to design parameters of the magnetic suspension sliding surface;
s3: constructing a solving model according to the difference between a function value and an actual value of a space plane equation constructed by point cloud data in an orbit coordinate system to obtain an initial plane;
s4: calculating the distance from all point cloud data in the fitting range to the initial plane, and removing unreasonable point cloud data according to a first set threshold value to obtain preprocessed point cloud data;
s5: calculating the mean value of the preprocessed point cloud data, taking the mean value as a standard value, calculating the residual error and the standard deviation of the preprocessed point cloud data, and removing the point cloud data of which the residual error is greater than a second set threshold value or the standard deviation is greater than a third threshold value to obtain single-segment point cloud data;
s6: and fitting the single-section point cloud data to obtain the point cloud data of the magnetic suspension sliding surface.
Specifically, the method comprises the following steps:
in S1, a high-precision laser scanner is adopted to dynamically scan the magnetic suspension track to obtain point cloud data.
In S2, in order to determine a fitting range, engineering coordinates (x, y, h) need to be converted into track coordinates (1, q, h), wherein a track coordinate system comprises an axis I, an axis q and an axis h, and the axis I represents mileage and is overlapped with a designed track surface central line; the q axis represents the offset from the point to the central line and is perpendicular to the central line of the rail surface; the h axis represents the elevation, and the origin of the track coordinate system is the central point of the initial position; the established orbital coordinate system is shown in fig. 3.
According to the design parameters of the magnetic levitation sliding surface, obtaining the offset range [ q ] from the boundary of two sides of the sliding surface to the central line under the design scheme 1 ',q 2 ']Given a certain threshold range, [ l ] 1 ',l 2 ']Ensuring the plane characteristic of the sliding surface in a small area range according to [ q ] 1 ',q 2 ']And [ l 1 ',l 2 ']And obtaining a fitting range.
S3, because the point cloud data in the area are basically concentrated near the sliding surface, the elevation of the point cloud data in the fitting range is normally distributed, and the height threshold range [ h ] is determined 1 ',h 2 ']And after removing the point cloud data with lower confidence coefficient according to the elevation, constructing a space plane equation by using the residual point cloud data.
In S3, the spatial plane equation is:
p(l,q,h)=al+bq+ch+m=0
wherein, a, b and c are normal vectors of a plane equation, l, q and h are mileage, offset and elevation respectively, and m is the distance from an origin to a space plane.
In S3, the expression of the solution model is as follows:
wherein e is the difference between the function value and the actual value, i represents the number of point cloud data, and the initial plane equation coefficient is obtained by solving the modelThe initial plane obtained by the method is a sliding plane in an ideal state, and in fact, each coordinate direction has a certain error, so that the initial plane has certain errorsAnd cannot be used as a best-fit plane for the sliding surface.
In S4, the first threshold is the prior distance t from the point cloud data to the initial plane; in S5, the second threshold is a variance σ of the residual error.
The third threshold is three times the standard deviation.
In S6, the plane parameter of the single-segment point cloud data is represented as (a) i ,b i ,c i ) In fig. 4, the left side is an initial plane, and the right side is a plane of single-segment point cloud data; when the magnetic suspension line is designed, the influence of the factors of superelevation, curvature and gradient is caused, the real state of the track plane cannot be simulated by the fitting of the integral point cloud data, the way of piecewise fitting the plane is adopted, the fitting of single-section point cloud data is carried out along the axis I of the track coordinate system, the stepping distance is delta 1 every time, and then the step S3 is carried out until the extraction of the integral sliding surface point cloud data is completed.
The method for extracting the point cloud data of the magnetic suspension sliding surface can effectively eliminate abnormal values, realize accurate extraction of the information of the magnetic suspension sliding surface and provide a basis for detection and evaluation of the magnetic suspension track sliding surface.
Claims (7)
1. A method for extracting point cloud data of a magnetic suspension sliding surface is characterized by comprising the following steps:
s1: acquiring point cloud data;
s2: converting the point cloud data from a space coordinate system to a track coordinate system, and obtaining a fitting range of a sliding surface according to design parameters of the magnetic suspension sliding surface;
s3: constructing a solving model according to the difference between a function value and an actual value of a space plane equation constructed according to the point cloud data in an orbit coordinate system to obtain an initial plane;
s4: calculating the distance from all point cloud data in the fitting range to the initial plane, and removing unreasonable point cloud data according to a first set threshold value to obtain preprocessed point cloud data;
s5: calculating the mean value of the preprocessed point cloud data, taking the mean value as a standard value, calculating the residual error and the standard deviation of the preprocessed point cloud data, and removing the point cloud data of which the residual error is greater than a second set threshold value or the standard deviation is greater than a third threshold value to obtain single-segment point cloud data;
s6: fitting the single-section point cloud data to obtain point cloud data of a magnetic suspension sliding surface;
in S2, the track coordinate system comprises an axis I, an axis q and an axis h, wherein the axis I represents mileage, the axis q represents offset from a point to a center line, the axis h represents elevation, and an origin of the track coordinate system is a central point of an initial position;
in S3, the spatial plane equation is:
p(l,q,h)=al+bq+ch+m=0
wherein a, b and c are normal vectors of a plane equation, l, q and h are mileage, offset and elevation respectively, and m is the distance from an original point to a space plane;
in S3, the expression of the solution model is:
wherein e is the difference between the function value and the actual value, and i represents the number of the point cloud data.
2. The method for extracting the cloud data of the magnetic levitation sliding surface as claimed in claim 1, wherein in S1, the cloud data is obtained by a laser scanner.
3. The method for extracting the point cloud data of the magnetic levitation sliding surface according to claim 1, wherein in S3, the elevation of the point cloud data in the fitting range is normally distributed, and after the point cloud data with lower confidence coefficient is removed according to the elevation, a spatial plane equation is constructed by using the rest point cloud data.
4. The method for extracting point cloud data of a magnetic levitation sliding surface as claimed in claim 1, wherein in S4, the first set threshold is a prior distance t from the point cloud data to an initial plane.
5. The method as claimed in claim 1, wherein in S5, the second threshold is a variance σ of the residual error.
6. The method as claimed in claim 1, wherein the third threshold is three times the standard deviation.
7. The method for extracting point cloud data of a magnetic levitation sliding surface as claimed in claim 1, wherein in S6, fitting of a single segment of point cloud data is performed along an axis l of a track coordinate system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011535988.3A CN112581579B (en) | 2020-12-22 | 2020-12-22 | Method for extracting point cloud data of magnetic suspension sliding surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011535988.3A CN112581579B (en) | 2020-12-22 | 2020-12-22 | Method for extracting point cloud data of magnetic suspension sliding surface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112581579A CN112581579A (en) | 2021-03-30 |
CN112581579B true CN112581579B (en) | 2022-11-18 |
Family
ID=75138877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011535988.3A Active CN112581579B (en) | 2020-12-22 | 2020-12-22 | Method for extracting point cloud data of magnetic suspension sliding surface |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112581579B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594250A (en) * | 2018-05-15 | 2018-09-28 | 北京石油化工学院 | A kind of point cloud data denoising point methods and device |
CN110728693A (en) * | 2019-09-27 | 2020-01-24 | 上海图趣信息科技有限公司 | Method and system for generating three-dimensional point cloud of large-scale driving road surface |
WO2020118970A1 (en) * | 2018-12-10 | 2020-06-18 | 华东理工大学 | On-machine inspection and compensation method employing point clouds and applied to complex surface processing |
EP3677980A1 (en) * | 2018-08-15 | 2020-07-08 | Guangzhou Xaircraft Technology Co., Ltd. | Method, device, and equipment for obstacle or ground recognition and flight control, and storage medium |
CN111461138A (en) * | 2020-03-05 | 2020-07-28 | 云南电网有限责任公司玉溪供电局 | Power line automatic extraction method based on L iDAR point cloud data |
CN111968253A (en) * | 2020-07-09 | 2020-11-20 | 北京工业大学 | Point cloud data-based road surface extraction method and system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109325998B (en) * | 2018-10-08 | 2023-06-30 | 香港理工大学 | Indoor 3D modeling method, system and related device based on point cloud data |
-
2020
- 2020-12-22 CN CN202011535988.3A patent/CN112581579B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594250A (en) * | 2018-05-15 | 2018-09-28 | 北京石油化工学院 | A kind of point cloud data denoising point methods and device |
EP3677980A1 (en) * | 2018-08-15 | 2020-07-08 | Guangzhou Xaircraft Technology Co., Ltd. | Method, device, and equipment for obstacle or ground recognition and flight control, and storage medium |
WO2020118970A1 (en) * | 2018-12-10 | 2020-06-18 | 华东理工大学 | On-machine inspection and compensation method employing point clouds and applied to complex surface processing |
CN110728693A (en) * | 2019-09-27 | 2020-01-24 | 上海图趣信息科技有限公司 | Method and system for generating three-dimensional point cloud of large-scale driving road surface |
CN111461138A (en) * | 2020-03-05 | 2020-07-28 | 云南电网有限责任公司玉溪供电局 | Power line automatic extraction method based on L iDAR point cloud data |
CN111968253A (en) * | 2020-07-09 | 2020-11-20 | 北京工业大学 | Point cloud data-based road surface extraction method and system |
Non-Patent Citations (1)
Title |
---|
《中低速磁浮列车-简支梁***耦合振动试验研究》;耿杰等;《铁道学报》;20180228;第40卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112581579A (en) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113870123B (en) | Automatic detection method for contact net leading height and pulling value based on vehicle-mounted mobile laser point cloud | |
CN113009456B (en) | Vehicle-mounted laser radar data calibration method, device and system | |
CN102561118B (en) | Quick fine tuning method for ballastless track of high-speed rail based on track deviation | |
CN111336990A (en) | Tunnel section convergence rapid analysis method and device | |
CN110490918A (en) | Existing railway track switch frog extracting method based on three-dimensional laser scanning technique | |
CN110807772B (en) | Bounding box-based irrelevant point cloud eliminating method in component size detection | |
CN111047542B (en) | Power line point supplementing method | |
CN108457143B (en) | Track line coordinate measuring system | |
CN102914290A (en) | Metro gauge detecting system and detecting method thereof | |
KR101755687B1 (en) | System of estimating camaera pose and method thereof | |
CN108534765B (en) | Track line coordinate measuring method | |
CN110793501A (en) | Subway tunnel clearance detection method | |
CN114396871A (en) | Prefabricated pier column installation position attitude monitoring method based on three-dimensional laser scanning | |
CN112363178A (en) | Dynamic measurement method for overhead line system tunnel suspension posts | |
CN110490342B (en) | Contact net static geometrical parameter detection method based on Faster R-CNN | |
CN112581579B (en) | Method for extracting point cloud data of magnetic suspension sliding surface | |
CN113291207B (en) | Dynamic measurement method of rigid contact network of subway | |
CN112950493B (en) | Method and device for correcting image distortion of linear array camera of rail train | |
CN112950562A (en) | Fastener detection algorithm based on line structured light | |
CN112158693A (en) | Detection method for elevator guide rail parameters | |
CN116893685A (en) | Unmanned aerial vehicle route planning method and system | |
CN116681912A (en) | Rail gauge detection method and device for railway turnout | |
CN114396892B (en) | Track curvature measuring method for track traffic curve | |
CN113063422B (en) | Clock-error-considered 5G terminal indoor positioning method | |
CN106875378B (en) | A kind of power line foreign matter detecting method |
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 |