CN104596415A - Method and device for determining lower edge of stack based on laser scanning single lines - Google Patents

Method and device for determining lower edge of stack based on laser scanning single lines Download PDF

Info

Publication number
CN104596415A
CN104596415A CN201410837330.6A CN201410837330A CN104596415A CN 104596415 A CN104596415 A CN 104596415A CN 201410837330 A CN201410837330 A CN 201410837330A CN 104596415 A CN104596415 A CN 104596415A
Authority
CN
China
Prior art keywords
overbar
laser scanning
scanning single
single line
analyzing spot
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
Application number
CN201410837330.6A
Other languages
Chinese (zh)
Other versions
CN104596415B (en
Inventor
谷红伟
许文海
王元臣
华夏
吴厚德
邬启凌
张望
王子乾
李晓伟
李瑛�
王鹤
杨德山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Shenhua Energy Co Ltd
Shenhua Sales Group Co Ltd
Beijing Shenhua Hengyun Energy Technology Co Ltd
Original Assignee
China Shenhua Energy Co Ltd
Shenhua Sales Group Co Ltd
Beijing Shenhua Hengyun Energy Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Shenhua Energy Co Ltd, Shenhua Sales Group Co Ltd, Beijing Shenhua Hengyun Energy Technology Co Ltd filed Critical China Shenhua Energy Co Ltd
Priority to CN201410837330.6A priority Critical patent/CN104596415B/en
Publication of CN104596415A publication Critical patent/CN104596415A/en
Application granted granted Critical
Publication of CN104596415B publication Critical patent/CN104596415B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention discloses a method for determining a lower edge of a stack based on laser scanning single lines. The method comprises the following steps: 1, acquiring a spatial coordinate Pi=[Xi, Yi, Zi] of a scanning point on each laser scanning single line of a target stack; 2, calculating di; 3, comparing the obtained di with a threshold value Ti, wherein Ti=(N-1)/N*h-; when the kth scanning point meets a formula that dk is greater than Tk, the kth scanning point is determined as a boundary point of the scanning single line; a set of the boundary points of all laser scanning single lines of the target stack corresponds to the lower edge of the target stack. The invention also provides a device for determining the lower edge of the stack based on the laser scanning single lines. By using the method and the device, the lower edge of the target stack can be accurately determined.

Description

A kind of method and apparatus determining heap body lower limb based on laser scanning single line
Technical field
The present invention relates to heap swept-volume technical field, be specifically related to a kind of method and apparatus determining to pile body lower limb based on laser scanning single line.
Background technology
Along with speeding up the construction of large-scale coal storage base, the modern management technique that needs to stockpile and store ensures, and be an important step of stockpiling and storing in modern management to making an inventory of coal volume mass, measuring the volume of coal storage base dump heap body and quality rapidly, accurately, is simply the important process that scientific management, assessment of economic benefit and memory space assessment are carried out in coal storage base, has important using value.
Along with the development of laser technology, for dish coal brings a kind of contactless metering system, laser coal checking method conventional at present mainly contains portable disc coal instrument and fixed dish coal instrument.Portable disc coal instrument, by setting relevant sweep parameter and mobile device completes measurement automatically, more accurately obtains coordinate information by GPS location simultaneously, sets up three-dimensional model.Fixed dish coal instrument, can be fixedly mounted on the equipment such as stacker-reclaimer and to carry out making an inventory fast, accurately to fixing coal yard.Although what fixed dish coal instrument decreased portable disc coal instrument manually gets operation ready, there is inevitable scan blind spot, all effectively cannot make an inventory to each dump in whole stockyard.Laser scanner, while scanning dump, inevitably scans the ground data on dump side.When processing these ground datas, some dish coal instrument only will intercept required cloud data in angular range by setting laser scanning angle scope, this mode is applicable to fixed measurement, require that dump Edge divider is easily distinguished, but dump is because of different, therefore sufficiently large in set angle scope, still can introduce ground data; Some dish coal instrument are not then removing ground data, determined the marginal information of dump by the mode of Digital Image Processing after forming three-dimensional laser point cloud, useless point is not removed from laser point cloud itself, ask for heap body volume by the method for image procossing, and ground cloud data inevitably brings error to heap body volume computing in volume computing.
Summary of the invention
In view of this, the present invention proposes a kind of method and apparatus determining to pile body lower limb based on laser scanning single line, by carrying out a series of process to laser scanning single line cloud data, determines the lower limb of target heap body.
The present invention is a kind of determines the method for piling body lower limb based on laser scanning single line, comprising:
Step one: the volume coordinate P of analyzing spot on the laser scanning single line of acquisition target heap body i=[X i, Y i, Z i], wherein P ibe i-th analyzing spot, i=1,2,3 ..., n;
Step 2: calculate according to following formula
h ‾ - = Z ‾ N - - Z i
h ‾ + = Z i - Z ‾ N +
d ‾ i = h ‾ - - h ‾ +
Wherein, with be respectively at P ithe mean value obtained after N number of analyzing spot carries out windowing process is respectively got in front and back;
Step 3: by what obtain with threshold value T icompare, wherein wherein when a kth analyzing spot meets time, this kth analyzing spot is defined as the separation of this scanning single line; Wherein, the lower limb of the set of the separation of all laser scanning single lines of target heap body corresponding target heap body.
Preferably, wherein the span of N is 10 < N < 20.
Beneficial effect of the present invention: embodiment of the present invention single scan line point cloud carries out rim detection, complete local by single line and detect, result is more accurate; The windowing smoothing processing of the embodiment of the present invention, averages to windowing before and after the determination point of single line, thus eliminates the jitter error because hardware error brings, and ensures the accuracy of algorithm process; The threshold adaptive process of the embodiment of the present invention, for added frame size, carries out adaptive analyzing and processing to decision threshold, ensures the accurate reliability of threshold value all the time.Utilize the present invention can obtain reliably, pile the separatrix on body lower limb and ground accurately, remove ground garbage, obtain and pile body information accurately, avoid the terrain surface specifications impact in the past existed in heap body three-dimensional modeling, volume and Mass Calculation.
Accompanying drawing explanation
Fig. 1 is the view of piling body with laser scanner scans of the embodiment of the present invention.
The algorithm principle figure of Fig. 2 embodiment of the present invention.
The algorithm schematic diagram of windowing process is carried out in Fig. 3 embodiment of the present invention.
The one-line scanning point Z coordinate schematic diagram of body without windowing process is piled in Fig. 4 embodiment of the present invention.
The one-line scanning point Z coordinate schematic diagram of body after windowing process is piled in Fig. 5 embodiment of the present invention.
Fig. 6 is the apparatus structure block diagram determining to pile body lower limb based on laser scanning single line of the embodiment of the present invention.
Embodiment
Below in conjunction with accompanying drawing and specific embodiment, technical scheme of the present invention is described in detail.
The present invention is based on laser scanning single line and determine that the principle of piling body lower limb is: use laser scanner scans heap body to obtain piling the laser point cloud of body, due to heap surface profile variations, the vertical of different scanning point highly must present the state of successively decreasing from top to bottom, there is a separation in the position had a common boundary at heap body and ground, more than separation for heap body divides, be above ground portion below separation, if the separation of each sweep trace then can be determined, also just determine the lower limb of whole heap body.Ground data is all rejected the heap body point cloud that can obtain removably face and ring, for piling body three-dimensional modeling, volume and Mass Calculation etc., result accurately and reliably.
Below describe and the present invention is based on the algorithm principle that laser scanning single line determines to pile body lower limb.
First, laser scanner is utilized to obtain the laser scanning single line point cloud coordinate of heap body.With reference to figure 1, body is piled by laser scanner scans, such as hand-held two dimensional laser scanning instrument is advanced one week around dump, the laser point cloud data of dump can be obtained, these cloud datas are set of some sweep traces of dump top-to-bottom (comprising part ground), and the volume coordinate that can obtain any point P on a scanning single line is P i=[X i, Y i, Z i], wherein footmark i represents i-th analyzing spot, i=1,2,3 ..., n.
Utilize the volume coordinate of analyzing spot, mainly Z coordinate calculates d i.With reference to figure 2, for three adjacent analyzing spot P in single scan line i-1, P i, P i+1, Z coordinate is respectively Z i-1, Z i, Z i+1, then have:
h -=Z i-1-Z i
h +=Z i-Z i+1(1)
d i=h --h +
According to each analyzing spot in formula (1) successively traverse scanning single line, the d that each point is corresponding can be obtained ivalue.
Wherein, for ground point, each point height is identical, has h -=Z i-1-Z i=0, h +=Z i-Z i+1=0, the d that therefore ground point is corresponding i=h --h +=0; For the point on heap body, consecutive point height successively decreases, h -value and h +value relatively, therefore piles d corresponding to body point iclose to 0, h in the ideal situation -=h +, then the d of body point is piled i=h --h +=0.
Be appreciated that and calculate the d that on single line, each analyzing spot is corresponding iafterwards, can d corresponding to base area millet cake i=0 and d corresponding to heap body point iclose to 0, determine the point of interface piling body and ground.But, owing to there is equipment and scanning errors, in practical application scene, directly do not adopt d i, but by carrying out windowing process to the analyzing spot on single line, eliminating measured deviation, below describing in detail.
Due to the existence of equipment and scanning errors, the Z coordinate (Z of analyzing spot i-1, Z i, Z i+1) may measured deviation be there is, carry out windowing process for this reason, with reference to figure 3, before and after i-th point, respectively get N number of some windowing, get the mean value of the Z coordinate figure of N number of point, can Z be obtained n-and Z n+, the Z used when replacing formula (1) to calculate respectively i-1and Z i+1, have:
h &OverBar; - = Z &OverBar; N - - Z i
h &OverBar; + = Z i - Z &OverBar; N + - - - ( 2 )
d &OverBar; i = h &OverBar; - - h &OverBar; +
The process that the multiple spot windowing of through type (2) is averaged, can eliminate the height relief deviation put near the frontier point that caused by measuring error, the wherein desirable 10 < N < 20 of frame size.Single scan line Z coordinate schematic diagram before Fig. 4 and Fig. 5 respectively illustrates windowing process and after windowing process, windowing process can improve the accuracy of searching frontier point.
It is corresponding that employing formula (2) calculates each point wherein, in the ideal case, all ground points and heap body point are had point only on ground with heap body intersection, namely separation is corresponding be not 0, and separation is corresponding
Then, in order to determine in single line, which point is the separation on heap body and ground, setting threshold value T i:
T i = ( N - 1 ) / N * h &OverBar; - - - - ( 3 )
Wherein, the frame size N chosen when N is windowing process.For each point, by its correspondence with threshold value T icompare, if wherein a kth point meets then can determine that this kth point is exactly the separation of this scanning single line.All carry out above-mentioned process to each scanning single line in heap body point cloud, can obtain the separation of whole heap body point cloud, the set of these separations is exactly the lower limb of heap body.
Wherein, about threshold value determined by frame size self-adaptation by known its of its set-up mode.Due to heap body point close to 0, frontier point close to and containing forms N, therefore threshold value in formula (3) self-adaptative adjustment can be carried out according to the height of all analyzing spots and selected frame size, ensure the accuracy at every bar sweep trace edge.
According to foregoing description, for a scanning single line, be highly greater than (in Fig. 3, sequence number is less than the point of k) of separation for heap body point, (in Fig. 3, sequence number is greater than the point of k) that be highly less than separation is ground point.During data processing, ground point is rejected, retain heap body point, target that is complete, that do not comprise ground data heap volumetric laser cloud data can be obtained.
Correspondingly, can provide a kind of device determining to pile body lower limb based on laser scanning single line, with reference to figure 6, it comprises: laser scanner, computing module, threshold process module and separation determination module.
Wherein, laser scanner is for obtaining the laser scanning single line point cloud coordinate P of target heap body i=[X i, Y i, Z i]; Computing module is used for calculating according to formula (2) threshold process module is used for according to formula (3) calculated threshold T i; Separation determination module is for by analyzing spot on scanning single line with threshold value T icompare, when a wherein kth point is corresponding this kth point is defined as the separation of this scanning single line.
The embodiment of the present invention processes for single scan line point cloud, and factor affected by environment is little, and accuracy is high, and the three-dimensional modeling and the volume computing that can be follow-up heap body establish good basis.
Above, be described in detail in conjunction with specific embodiments to technical scheme of the present invention, described specific embodiment understands thought of the present invention for helping.The derivation that those skilled in the art make on the basis of the specific embodiment of the invention and modification also belong within scope.

Claims (4)

1. determine a method of piling body lower limb based on laser scanning single line, it is characterized in that, comprising:
Step one: the volume coordinate P of analyzing spot on the laser scanning single line of acquisition target heap body i=[X i, Y i, Z i], wherein P ibe i-th analyzing spot, i=1,2,3 ..., n;
Step 2: calculate according to following formula
h &OverBar; - = Z &OverBar; N - - Z i
h &OverBar; + = Z i - Z &OverBar; N +
d &OverBar; i = h &OverBar; - - h &OverBar; +
Wherein, with be respectively at P ithe mean value of the Z coordinate obtained after N number of analyzing spot carries out windowing process is respectively got in front and back;
Step 3: by what obtain with threshold value T icompare, wherein wherein when a kth analyzing spot meets time, this kth analyzing spot is defined as the separation of this scanning single line;
Wherein, the lower limb of the set of the separation of all laser scanning single lines of target heap body corresponding target heap body.
2. determine the method for piling body lower limb based on laser scanning single line as claimed in claim 1, it is characterized in that, wherein the span of N is 10 < N < 20.
3. determine based on laser scanning single line the device piling body lower limb, it is characterized in that, comprising:
Laser scanner, the volume coordinate P of analyzing spot on the laser scanning single line obtaining target heap body i=[X i, Y i, Z i], wherein P ibe i-th analyzing spot, i=1,2,3 ..., n;
Computing module, for calculating according to following formula
h &OverBar; - = Z &OverBar; N - - Z i
h &OverBar; + = Z i - Z &OverBar; N +
d &OverBar; i = h &OverBar; - - h &OverBar; +
Wherein, with be respectively at P ithe mean value of the Z coordinate obtained after N number of analyzing spot carries out windowing process is respectively got in front and back;
Threshold process module, for calculated threshold
Separation determination module, by what obtain with threshold value T icompare, wherein wherein when a kth analyzing spot meets time, this kth analyzing spot is defined as the separation of this scanning single line;
Heap body lower limb processing module, the separation for all laser scanning single lines according to target heap body forms the lower limb of target heap body.
4. determine based on laser scanning single line the device piling body lower limb as claimed in claim 3, it is characterized in that, wherein the span of N is 10 < N < 20.
CN201410837330.6A 2014-12-29 2014-12-29 A kind of method and apparatus that heap body lower edge is determined based on laser scanning single line Active CN104596415B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410837330.6A CN104596415B (en) 2014-12-29 2014-12-29 A kind of method and apparatus that heap body lower edge is determined based on laser scanning single line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410837330.6A CN104596415B (en) 2014-12-29 2014-12-29 A kind of method and apparatus that heap body lower edge is determined based on laser scanning single line

Publications (2)

Publication Number Publication Date
CN104596415A true CN104596415A (en) 2015-05-06
CN104596415B CN104596415B (en) 2017-06-09

Family

ID=53122359

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410837330.6A Active CN104596415B (en) 2014-12-29 2014-12-29 A kind of method and apparatus that heap body lower edge is determined based on laser scanning single line

Country Status (1)

Country Link
CN (1) CN104596415B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106127770A (en) * 2016-06-27 2016-11-16 感知控股集团有限公司 Contour measuring method and system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101726255A (en) * 2008-10-24 2010-06-09 中国科学院光电研究院 Method for extracting interesting buildings from three-dimensional laser point cloud data
CN101916329A (en) * 2010-07-30 2010-12-15 中国科学院空间科学与应用研究中心 Computational method for modeling volume of material pile
CN103090791A (en) * 2013-01-08 2013-05-08 中联重科股份有限公司 Measurement system, method and device of scattered materials and material piling and taking control system
KR101285131B1 (en) * 2013-04-04 2013-07-11 관동대학교산학협력단 Method of filtering lidar data using fuzzy set theory
CN103268609A (en) * 2013-05-17 2013-08-28 清华大学 Partition method for orderly extracting point cloud on ground
CN103390169A (en) * 2013-07-19 2013-11-13 武汉大学 Sorting method of vehicle-mounted laser scanning point cloud data of urban ground objects
KR101421128B1 (en) * 2014-04-30 2014-07-22 국방과학연구소 Extraction method of building regions using segmented 3d raw datd based on laser radar

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101726255A (en) * 2008-10-24 2010-06-09 中国科学院光电研究院 Method for extracting interesting buildings from three-dimensional laser point cloud data
CN101916329A (en) * 2010-07-30 2010-12-15 中国科学院空间科学与应用研究中心 Computational method for modeling volume of material pile
CN103090791A (en) * 2013-01-08 2013-05-08 中联重科股份有限公司 Measurement system, method and device of scattered materials and material piling and taking control system
KR101285131B1 (en) * 2013-04-04 2013-07-11 관동대학교산학협력단 Method of filtering lidar data using fuzzy set theory
CN103268609A (en) * 2013-05-17 2013-08-28 清华大学 Partition method for orderly extracting point cloud on ground
CN103390169A (en) * 2013-07-19 2013-11-13 武汉大学 Sorting method of vehicle-mounted laser scanning point cloud data of urban ground objects
KR101421128B1 (en) * 2014-04-30 2014-07-22 국방과학연구소 Extraction method of building regions using segmented 3d raw datd based on laser radar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
田茂义等: "一种车载激光扫描数据的滤波处理算法", 《测绘科学技术学报》 *
谭贲等: "车载激光扫描数据的地物分类方法", 《遥感学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106127770A (en) * 2016-06-27 2016-11-16 感知控股集团有限公司 Contour measuring method and system
CN106127770B (en) * 2016-06-27 2019-03-15 感知控股集团有限公司 Contour measuring method and system

Also Published As

Publication number Publication date
CN104596415B (en) 2017-06-09

Similar Documents

Publication Publication Date Title
US20150117753A1 (en) Computing device and method for debugging computerized numerical control machine
CN106595583B (en) RTK measurement receiver inclination measurement method
CN104613930B (en) Method and device for measuring distance as well as mobile terminal
KR101918168B1 (en) Method for performing 3D measurement and Apparatus thereof
CN105866762A (en) Laser-radar automatic calibration method and device thereof
CN110263075B (en) LADCP and USBL combined observation and data processing method
CN105136054A (en) Fine structure deformation monitoring method and system based on ground three-dimensional laser scanning
CN102930540A (en) Method and system for detecting contour of urban building
CN104019745B (en) Based on the free planar dimension measuring method of single visual feel indirect calibration method
CN106646404B (en) A kind of error correcting method and system of two strips front phased-array radar
CN112034431B (en) External parameter calibration method and device for radar and RTK
CN104916216A (en) Map construction method and system thereof
CN110646231B (en) Floor sweeping robot testing method and device
CN104330038A (en) Size measurement method
US20160123722A1 (en) Computing device and method for analyzing thickness
CN105335952A (en) Matching cost calculation method and apparatus, and parallax value calculation method and equipment
CN107367245B (en) Invalid point detection and elimination method in optical three-dimensional profile measurement
CN103065295A (en) Aviation and ground lidar data high-precision automatic registering method based on building angular point self-correction
CN105354575A (en) Image binarization threshold determination method in sea-sky-line detection
CN104596415A (en) Method and device for determining lower edge of stack based on laser scanning single lines
CN113748693B (en) Position and pose correction method and device of roadbed sensor and roadbed sensor
CN113487539A (en) Method, equipment, system and storage medium for analyzing quality of rubber road
CN103714528B (en) Object segmentation device and method
CN107478155A (en) Product inspection method, apparatus and system
CN109000560B (en) Method, device and equipment for detecting package size based on three-dimensional camera

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant