CN107127755A - A kind of real-time acquisition device and robot polishing method for planning track of three-dimensional point cloud - Google Patents
A kind of real-time acquisition device and robot polishing method for planning track of three-dimensional point cloud Download PDFInfo
- Publication number
- CN107127755A CN107127755A CN201710334080.8A CN201710334080A CN107127755A CN 107127755 A CN107127755 A CN 107127755A CN 201710334080 A CN201710334080 A CN 201710334080A CN 107127755 A CN107127755 A CN 107127755A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- robot
- point
- 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
- 238000005498 polishing Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000006073 displacement reaction Methods 0.000 claims abstract description 50
- 238000005259 measurement Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 31
- 230000006854 communication Effects 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 11
- 230000003044 adaptive effect Effects 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 238000010845 search algorithm Methods 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000003892 spreading Methods 0.000 claims description 3
- 230000007480 spreading Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000009452 underexpressoin Effects 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1602—Programme controls characterised by the control system, structure, architecture
- B25J9/1607—Calculation of inertia, jacobian matrixes and inverses
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1628—Programme controls characterised by the control loop
- B25J9/163—Programme controls characterised by the control loop learning, adaptive, model based, rule based expert control
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Automation & Control Theory (AREA)
- Manipulator (AREA)
Abstract
The invention discloses a kind of real-time acquisition device of three-dimensional point cloud, including robot, laser displacement sensor, robot real-time control system, the laser displacement sensor is arranged on the end of the robot by fixture, the robot real-time control system connects robot and laser displacement sensor by real-time industrial ethernet bus, for making the reading of laser displacement sensor and the pose of robot obtain synchronization, one-dimensional measurement is expanded into three-dimensional measurement, so as to obtain the three-dimensional point cloud of workpiece by scanning.The invention also discloses a kind of robot polishing method for planning track based on described device.Pose when positional information and normal direction information of the present invention based on workpiece three-dimensional point cloud are polished come planning robot, and corrected polishing track as robot to scan track, precision is high, with real-time and flexibility characteristics, and cost is low, disclosure satisfy that the processing of different workpieces.
Description
Technical field
Acquisition and robot polishing of the present invention applied to three-dimensional point cloud, be both related to based on laser displacement sensor and machine
Real-time Communication for Power obtains the device of three-dimensional point cloud between people, is directed to polishing method for planning track of the robot based on three-dimensional point cloud,
There is provided a set of new three-dimensional point cloud acquisition methods, apply to polishing field for robot and provide new test method and dress
Put.
Background technology
Now, the acquisition equipment of three-dimensional point cloud mainly has three-dimensional laser scanner, and two-dimensional line laser displacement sensor adds auxiliary
Motion is helped, above two method can obtain huge point off density cloud of counting.It can also be obtained in addition by three-coordinates measuring machine
Point cloud, the point cloud positional precision that this method is obtained is high, for sparse less cloud of counting.These equipment prices are expensive so that
Point cloud procurement cost is higher, with the appearance of consumer level products, such as Kinect so that point cloud is developed faster.
But the point cloud that the above method is obtained can be comprising the environmental information residing for target object, it is necessary to be filtered, go to original point cloud
Make an uproar, split and the processing such as registration, could obtain the point cloud of object, post-processing complexity.
Point cloud mainly applies to reverse-engineering, carries out the three-dimensionalreconstruction of existing object, is that it creates digital information storehouse;
Field of machine vision, point cloud can be used for the directions such as target identification, target three-dimensional values and map structuring, and vision is provided for robot
Information, aids in its decision-making and path planning.
Robot technology has been widely used in the fields such as welding, spraying, but the grinding workpieces of complicated shape are main still
By manual polishing and numerical control machine tooling, the former damages workers ' health because polishing noise, dust are big, and to the skill of workman
Art level requirement is higher;The latter's machining accuracy is high, but process equipment is expensive, and flexibility is poor to limit its extensive use.
Industrial robot has multiple degrees of freedom, and flexibility is good, is relatively adapted to apply to polishing field.At present, using robot
Off-line programming technique, successfully applies to complex part polishing by robot technology, but off-line programming technique requires right in advance
Instrument and workpiece coordinate system are demarcated, and the clamping error due to workpiece and processing inconsistency, cause robot polishing one
When criticizing workpiece, there is the inconsistent situation of precision, or even waster occur.
The content of the invention
To solve the problems, such as the acquisition of three-dimensional point cloud and a cloud being applied into robot polishing field, the invention provides one kind
Real-time acquisition device and robot the polishing method for planning track of three-dimensional point cloud based on laser displacement sensor.Design not similar shape
The workpiece of shape is tested, planning robot's scanning track, logical by bus based on bus communication protocol and real time operating system
Letter mode realizes the real-time Communication for Power of robot and laser displacement sensor, workpiece three-dimensional point cloud is obtained, for three-dimensional point cloud is used
In robot, polishing trajectory planning provides external condition.
To realize the above-mentioned purpose of the present invention, one aspect of the present invention provides a kind of real-time acquisition device of three-dimensional point cloud,
Including robot, laser displacement sensor, the robot real-time control system based on real time operating system, the laser displacement is passed
Sensor is arranged on the end of the robot, the artificial six degree of freedom articulated robot of machine, the robot by fixture
Real-time control system connects robot and laser displacement sensor by real-time industrial ethernet bus, for passing laser displacement
The reading of sensor and the pose of robot obtain synchronization, and one-dimensional measurement is expanded into three-dimensional measurement, so as to obtain work by scanning
The three-dimensional point cloud of part.
Further, it is Millisecond between described robot real-time control system and robot and laser displacement sensor
Communication.
Another aspect of the present invention provides a kind of robot polishing method for planning track based on described device, including step
Suddenly:
(1) by the measurement coordinate system of multi-point calibration method Calibration of Laser displacement transducer relative to robot end's coordinate system
Transition matrix;
(2) one-dimensional measurement is expanded to three-dimensional measurement, obtained by robotically-driven laser displacement sensor scanning workpiece, realization
The three-dimensional point cloud of workpiece;
(3) normal estimation is carried out to the three-dimensional point cloud, while the spatial positional information and method that are included based on three-dimensional point cloud
Pose and feed path when being polished to information planning robot, and corrected polishing rail as robot to scan track
Mark.
Further, the step (1) specifically includes:
(11) top in robot end's one demarcation of installation first, control machine people reaches a certain fixing point in space, adopts
With 5 standardizations of robot tool coordinate system, the fixing point is read relative to machine from teaching machine after the completion of tool coordinates system demarcation
The locus of device people's basis coordinates systemBp;
(12) and then in robot end's clamping laser displacement sensor, laser is radiated at the fixation by control machine people
Point, reads pose of the robot current end relative to basis coordinates systemAnd the registration of laser displacement sensorMP, the laser
The registration of displacement transducerMP and the locusBP transformational relation is:
In formulaBP --- coordinates of the spatial point p in robot basis coordinates system;
MCoordinates of p --- the spatial point p in measurement coordinate system;
--- homogeneous transform matrix of robot end's coordinate system relative to basis coordinates system;
--- homogeneous transform matrix of the coordinate system relative to robot end's coordinate system is measured, wherein,
(13) formula (1) is converted toIfThen (1) can be expressed as follows:
(14) multi-point calibration is based on, formula (3) is changed into AX=k form, an over-determined systems are constituted, and using most
A young waiter in a wineshop or an inn multiplies to solve, and form is as follows:
MakeMinimum solution X*As over-determined systems AX=k least square solution;
(15) solve and matrix is obtained after XObtain matrixAfterwards, you can so that the workpiece three-dimensional point cloud of scanning to be turned
Change under robot basis coordinates system.
Further, in step (2), described is real-time by robot the step of one-dimensional measurement is expanded into three-dimensional measurement
Control system, the x of workpiece point cloud, y location is sensed from robot end's pose, the z location of workpiece point cloud by laser displacement
Device is understood, passes through real-time bus mechanics of communication, it is ensured that the synchronization of x, y, z, one-dimensional measurement is expanded into three-dimensional measurement so as to realize.
Further, described step (3) specifically includes step:
(31) spatial positional information of each point based on three-dimensional point cloud and the normal direction information planning machine after optimizing and revising
People's polishing pose;
(32) the processing node of defining point cloud, and processing nodal information is sent to machine successively to scan sequencing
People, robot, which carries out inverse kinematics and moving interpolation, will process the continuous robot polishing rail of node series connection generation one
Mark.
Further, described step (31) specifically includes step:
(301) determine polishing when cutter position, i.e., using each point in three-dimensional point cloud included based on robot base
Mark the spatial positional information of system;Position when make it that cutter reaches polishing by Robotic inverse kinematics control machine people;
(302) posture of cutter during polishing is determined, i.e., each point based on three-dimensional point cloud estimates normal direction information and optimization
The posture of cutter when normal direction information planning robot after adjustment polishes.
Further, described step (302) specifically includes step:
(311) the closest points of k are searched for by improved k Neighborhood-region-search algorithms, i.e., using a point as body-centered, with
Fixed step size builds the square of itself, and k neighborhood search is then carried out in square, if finding k closest points, i.e.,
Search is terminated, otherwise expands the scope of square with fixed step size, k closest points are re-searched for, until in square
Inside search k closest point;
(312) improved PCA algorithms are used, different weighting functions is assigned by the point in the k neighborhoods to certain point, uses
Least square fitting method fits a plane, using the normal direction of institute's fit Plane as the point initial normal direction;It is next based on institute
State initial normal direction and adaptively pinpoint the weighting function of each point in cloud really, then initial normal direction is optimized;
(313) point cloud is redirected, i.e., the normal direction after optimization is adjusted using minimum spanning tree method so that point cloud normal direction
Point to consistent.
Further, described step (312) specifically includes step:
(321) two weighting functions are introduced into PCA algorithms so that approximating method has to the exterior point and noise in a cloud
Certain robustness, least square fitting is represented by:
S.t. | | n | |=1 (5)
ω in formulad(xi) --- Gauss weighting function so that effect of the point near apart from the point to fit Plane is a little big,
Otherwise influence is small, ωd(xi) be expressed as:σdFor apart from bandwidth, its initial value is according to the shape of workpiece
And the distribution situation of scanning element cloud is preset;
ωr(ri) --- the Gauss weighting function relevant with regression criterion, ωr(ri) be expressed as:σrFor
Regression criterion bandwidth, its initial value is preset according to the shape of workpiece and the distribution situation of scanning element cloud;
ri m--- the m times iteration, xiThe regression criterion of point, ri mIt is represented by:ri m=dm+(xi-x)Tnm;
(322) formula (5) construction covariance matrix C is solved, is expressed as follows:
In formula--- the center of neighborhood is represented, is expressed as:
K --- the points in neighborhood are represented,
By solving covariance matrix C, its corresponding characteristic vector of minimum characteristic value is the normal direction n of fit Plane,
I.e. as the initial normal direction of the point, and
(323) a characteristic coefficient F is defined for each point in point cloudcIt is positioned at smooth region, sharp spy to quantify the point
Region or transitional region are levied, characteristic coefficient is represented by:
S=2 in formula, represents when ratio is 0.5, then to judge that the point is located at sharp features region;
δl--- represent the average distance of the point 6 point closest with it;
(324) characteristic coefficient F is determinedcThe first defining point cloud of value in point cloud cluster in a certain vertex neighborhood, when certain point cloud
Points in cluster are more than the business of the points and the quantity of point cloud cluster in neighborhood, then it is assumed that this cloud cluster is principal point cloud cluster;In order to certainly
The quantity of dynamic acquisition point cloud cluster, sets an angle threshold, when 2 points of normal direction angles are more than threshold value, that is, thinks at this 2 points
In different point cloud clusters;A quantity of cloud cluster is obtained, so as to calculate the corresponding characteristic coefficient F of each pointc, a, 2 points of b's
Normal direction angle is represented by:
D (a, b)=cos-1(na,nb); (8)
(325) F is determinedcValue after, it is adaptive calculate each point apart from bandwidth σd, regression criterion bandwidth σr, normal direction
Difference bandwidth σn, it is expressed as follows:
σd=(1+Fc)·δl (9)
σr=rmax/3Fc (10)
σn=| | ni-n||max/3Fc; (11)
(326) the 3rd weighting function is introduced in PCA algorithms, then least square fitting is expressed as:
S.t. | | n | |=1 (12)
ω in formulan(ni) --- the weighting function relevant with normal direction difference, it is expressed as:σnFor normal direction
Difference bandwidth, niFor certain initial normal direction of point, gained is solved by formula (6);
(327) calculated by adaptometer in step (325) each put apart from bandwidth σd, regression criterion bandwidth σr, method
To difference bandwidth σnThe weighting function each put in point cloud is calculated, is replaced in covariance matrix C
Then formula (12) is solved by the covariance matrix C after replacement, obtains the optimization normal direction of each point.
Further, the step (313) specifically includes step:
(331) the maximum point of z coordinate, as root node, adjusts its normal direction and is allowed to and vector (0,0,1) in selected point cloud
Dot product is more than 0, and this will be such that the normal orientation after adjustment points on the outside of workpiece;
(332) weight c=1- is set | ni·nj|, this ensure that the nonnegative weights, and represent two neighbor points when c is smaller
Normal orientation it is more parallel;
(333) traveled through according to weight c size, if father node niWith child node njMeet ni·nj<0, then by njInstead
To;
(334) all nodes are traveled through, normal orientation is carried out to set a threshold value e in Spreading and diffusion, communication process, works as c
During≤e, multiple spot in One Diffusion Process to neighborhood works as c>During e, one point of One Diffusion Process, after the completion of judge whether to meet c≤e, then enter
Row multiple spot spreads, until being diffused into whole point cloud space;
(335) if failing to travel through all points, point cloud is not connected, and partial data forms local isolated island, if do not traveled through
Count number less than point cloud sum an one thousandth, then give up as noise, otherwise search in the data set not traveled through with
Closest point between having traveled through a little, the starting point traveled through as continuation simultaneously goes to step (334).
Compared with prior art, the present invention to scanning element cloud by carrying out normal estimation, the position based on workpiece three-dimensional point cloud
Information and normal direction information carry out pose during planning robot's polishing, and are corrected polishing rail as robot to scan track
Mark, reliable and stable, precision is high, with real-time and flexibility characteristics, and cost simple in construction is low, disclosure satisfy that different workpieces plus
Work.
Brief description of the drawings
Fig. 1 is the point cloud real-time acquisition device schematic diagram of the present invention;
Fig. 2 is caliberating device schematic diagram.
In figure:1- robots, 2- laser displacement sensors;3- robots real-time control system;4- workpiece;5- demarcation top
Point.
Embodiment
To further understand the present invention, the present invention will be further described with reference to the accompanying drawings and examples, but needs
Illustrate, the scope of protection of present invention is not limited to the scope of embodiment statement.
Embodiment
As shown in figure 1, a kind of real-time acquisition device of three-dimensional point cloud, including robot 1, laser displacement sensor 2, it is based on
The robot real-time control system 3 of real time operating system, the laser displacement sensor 2 is arranged on the robot by fixture
1 end, the robot 1 be six degree of freedom articulated robot, the robot real-time control system 3 by real-time industrial with
Too network bus connects robot 1 and laser displacement sensor 2, for making the reading of laser displacement sensor 2 and the position of robot 1
Appearance obtains synchronization, and one-dimensional measurement is expanded into three-dimensional measurement, so as to obtain the three-dimensional point cloud of workpiece 4 by scanning.
Described robot real-time control system between robot 1 and laser displacement sensor 2 for Millisecond with communicating.
During laser scanning, six axle crossmachine people are controlled to move in the way of teaching playback, with laser displacement sensor 2
Workpiece 4 is scanned, the point on workpiece 4 is gathered in real time based on apparatus of the present invention, finally gives the work based on robot basis coordinates system
Part three-dimensional point cloud.Point cloud real-time acquisition device is as shown in Figure 1.
The robot real-time control system 3 of this device is built based on real time operating system and EPA bussing technique,
The x of workpiece point cloud, y location is from robot end's pose, and the z location of workpiece point cloud is led to from laser displacement sensor 2
Cross real-time bus mechanics of communication, it is ensured that the synchronization of x, y, z, one-dimensional measurement is expanded into three-dimensional measurement, and the control so as to realize
The communication cycle of system meets requirement of real-time up to 1ms, can use the three-dimensional point cloud of device scanning generation workpiece.
The shown end clamping laser displacement sensor 2 of robot 1, reads laser displacement by analog input and output module and passes
The analog quantity output signals of sensor are simultaneously connected, then total based on EPA with the pulse signal of six controlled motors of robot
All signal acquisitions are returned to carry out data processing, communication cycle energy by the real time operating system of line technology by bus communication mode
Millisecond is reached, so that the reading synchronization of the pose and laser displacement sensor of robot is realized, by the pose of robot 1 and laser
The distance of displacement transducer 2 is solved by robot transformation equation, obtains the workpiece three-dimensional point cloud under robot basis coordinates system.
The present embodiment utilizes the workpiece of the different surface configuration of SolidWorks Software for Design, for sweep test and polishing
Experiment.This experimental design S-shaped curve surface work pieces, using aluminium alloy, and are processed, for test experiments.
A kind of robot polishing method for planning track based on described device, including step:
(1) by the measurement coordinate system of multi-point calibration method Calibration of Laser displacement transducer 2 relative to the ending coordinates of robot 1
The transition matrix of system, in order to which the point cloud measured under coordinate system is transformed under robot basis coordinates system, must shift to an earlier date Calibration of Laser position
The measurement coordinate system of displacement sensor 2 is relative to the transition matrix of robot end's coordinate system, abbreviation hand and eye calibrating;
(2) one-dimensional measurement is expanded to three-dimensional measurement, obtained by the driving of robot 1 laser displacement sensor 2 scanning workpiece, realization
Take the three-dimensional point cloud of workpiece;
(3) by setting three-dimensional point cloud described in different adaptive weighting function pairs to carry out normal direction to each point in three-dimensional point cloud
Estimation and optimization, while the position when spatial positional information included based on three-dimensional point cloud and the polishing of normal direction information planning robot 1
Appearance and feed path, and corrected continuous polishing track as robot 1 to scan track.
Specifically, the step (1) specifically includes:
(11) demarcation top 5 is installed in the end of robot 1 first, control machine people 1 reaches a certain fixing point in space,
It is relative that from teaching machine the fixing point is read using robot 1 tool coordinates system, 5 standardizations, after the completion of tool coordinates system demarcation
In the locus of the basis coordinates system of robot 1Bp;
(12) and then in the end clamping laser displacement sensor 2 of robot 1, laser is radiated at described by control machine people 1
Fixing point, reads pose of the current end of robot 1 relative to basis coordinates systemAnd the registration of laser displacement sensor 2MP, institute
State the registration of laser displacement sensor 2MP and the locusBP transformational relation is:
In formulaBP --- coordinates of the spatial point p in the basis coordinates system of robot 1;
MCoordinates of p --- the spatial point p in measurement coordinate system;
--- homogeneous transform matrix of the ending coordinates system of robot 1 relative to basis coordinates system;
--- homogeneous transform matrix of the coordinate system relative to ending coordinates system of robot 1 is measured, wherein,
(13) formula (1) is converted toIfThen (1) can be expressed as follows:
(14) in order to reduce calibrated error, based on multi-point calibration, formula (3) is changed into AX=k form, one is constituted and surpasses
Determine equation group, and solved using least square, form is as follows:
MakeMinimum solution X*As over-determined systems AX=k least square solution;
(15) solve and matrix is obtained after XObtain matrixAfterwards, you can so that the workpiece three-dimensional point cloud of scanning to be turned
Change under the basis coordinates system of robot 1.
Specifically, in step (2), described is real by robot 1 the step of one-dimensional measurement is expanded into three-dimensional measurement
When control system, the x of workpiece point cloud, y location passes from the end pose of robot 1, the z location of workpiece point cloud by laser displacement
Sensor 2 is understood, passes through real-time bus mechanics of communication, it is ensured that the synchronization of x, y, z, one-dimensional measurement is expanded into three-dimensional survey so as to realize
Amount.
Can whether the point cloud obtained by said apparatus can really reflect workpiece surface shape facility, as follow-up
The polishing trajectory planning data of robot 1, can be verified by following methods.
When carrying out replication experiment, its three-dimensional CAD model known to the workpiece 4 processed can be discrete by three-dimensional CAD model
Into a cloud as standard point cloud, because standard point cloud is, based on workpiece coordinate system, workpiece can be demarcated in advance relative to machine
The homogeneous transform matrix of the basis coordinates system of people 1, first some points are marked relative to the coordinate of workpiece coordinate system during demarcation on workpiece, then
Control machine people 1 arrives corresponding points, obtains coordinate of this in the basis coordinates system of robot 1.Obtain at same o'clock in two coordinate systems
Under expression, then the transition matrix for calculating two coordinate systems, i.e. workpiece coordinate system and the basis coordinates of robot 1 are connect by svd algorithm
The transition matrix of system.Standard point cloud can be just transformed into after the completion of demarcation under basis coordinates system, then with scanning gained workpiece point cloud
Contrasted, that is, obtain the error existed between scanning element cloud and standard point cloud, so as to verify the feasibility of the device.Due to machine
The presence of the position error of device people 1, the measurement error and calibrated error of laser displacement sensor 2, should allow scanning element cloud and standard
There is error between point cloud, but this error should be controlled within zone of reasonableness.
Specifically, described step (3) specifically includes step:
(31) spatial positional information of each point based on three-dimensional point cloud and by adaptive weighting function optimization adjust after
Normal direction information planning robot 1 polishing pose, the polishing pose planning dependent on three-dimensional point cloud gather accuracy and
The degree of accuracy of point cloud normal estimation, and finally determine the polishing quality of workpiece;
(32) the processing node of defining point cloud, and processing nodal information is sent to robot successively to scan sequencing
1, robot 1 carries out inverse kinematics and moving interpolation and will process one continuous robot 1 of node series connection generation polishing rail
Mark.
Specifically, described step (31) specifically includes step:
(301) position of cutter during polishing is determined, because containing the sky that scanning workpiece surface is each put in three-dimensional point cloud
Between positional information, and it is based on the basis coordinates system of robot 1, so as to utilize each point institute in three-dimensional point cloud to put the positional information of cloud
Comprising the spatial positional information based on the basis coordinates system of robot 1;Knife is caused by resolved motion control robot 1 of robot 1
Tool reaches position during polishing;
(302) posture of cutter during polishing is determined, that is, passes through the pre- of each point of adaptive weighting function pair three-dimensional point cloud
The posture of cutter when normal direction information planning robot 1 after estimating normal direction information and optimizing and revising polishes, therefore accurately point cloud method
There is material impact to the polishing quality of workpiece to estimation.
Specifically, described step (302) specifically includes step:
(311) the closest points of k are searched for by improved k Neighborhood-region-search algorithms, i.e., using a point as body-centered, with
Fixed step size builds the square of itself, and k neighborhood search is then carried out in square, if finding k closest points, i.e.,
Search is terminated, otherwise expands the scope of square with fixed step size, k closest points are re-searched for, until in square
Inside search k closest point, this avoid search for each neighborhood of a point when will travel through a little, improve search
Efficiency;And traditional k Neighborhood-region-search algorithms, be traversal point cloud in institute a little, calculate corresponding Euclidean distance, and carry out ascending order
Arrangement, takes out k closest points and is used as k neighborhoods;This algorithm is required for traversal point cloud when calculating each neighborhood of a point
In each point, take it is longer, do not apply to large-scale point cloud situation;
(312) improved PCA algorithms are used, different weighting functions is assigned by the point in the k neighborhoods to certain point, uses
Least square fitting method fits a plane, using the normal direction of institute's fit Plane as the point initial normal direction;It is next based on institute
State initial normal direction and adaptively pinpoint the weighting function of each point in cloud really, then initial normal direction is optimized;Traditional PCA algorithms
It is in the k neighborhoods of certain point, using least square fitting method, to fit a plane in a cloud, the normal direction of fit Plane is made
For the normal direction of the point, due to inevitably there are noise, exterior point and sharp features in a cloud, so this step is to PCA algorithms
It is improved, different weighting functions is assigned to the point in neighborhood, so as to improve the quality of fit Plane, the final normal direction that improves is estimated
The precision of meter.
(313) point cloud is redirected, i.e., the normal direction after optimization is adjusted using minimum spanning tree method so that point cloud normal direction
Point to unanimously, minimum spanning tree method is:First provide that the normal direction of any is pointed to, diffusion path then generated according to the weight of each point,
It is final to be diffused into whole point cloud space, complete the adjustment of point cloud normal direction.
Specifically, described step (312) specifically includes step:
(321) plane that be fitted PCA algorithms accurately reflects the normal direction at the point, should try one's best in guarantee fitting neighborhood
Point is located on same dough sheet, but conventional measurement means can make space exterior point and noise spot with measured deviation there is scanning
In point cloud, therefore two weighting functions are introduced into PCA algorithms and cause approximating method has to the exterior point and noise in a cloud certain
Robustness, least square fitting is represented by:
S.t. | | n | |=1 (5)
ω in formulad(xi) --- Gauss weighting function so that effect of the point near apart from the point to fit Plane is a little big,
Otherwise influence is small, ωd(xi) be expressed as:σdFor apart from bandwidth, its initial value is according to the shape of workpiece
And the distribution situation of scanning element cloud is preset;
ωr(ri) --- the Gauss weighting function relevant with regression criterion, ωr(ri) be expressed as:σrFor
Regression criterion bandwidth, its initial value is preset according to the shape of workpiece and the distribution situation of scanning element cloud;
ri m--- the m times iteration, xiThe regression criterion of point, ri mIt is represented by:ri m=dm+(xi-x)Tnm;
(322) formula (5) construction covariance matrix C is solved, is expressed as follows:
In formula--- the center of neighborhood is represented, is expressed as:
K --- the points in neighborhood are represented,
By solving covariance matrix C, its corresponding characteristic vector of minimum characteristic value is the normal direction n of fit Plane,
I.e. as the initial normal direction of the point, and
(323) a characteristic coefficient F is defined for each point in point cloudcIt is positioned at smooth region, sharp spy to quantify the point
Region or transitional region are levied, characteristic coefficient is represented by:
S=2 in formula, represents when ratio is 0.5, then to judge that the point is located at sharp features region;
δl--- represent the average distance of the point 6 point closest with it;
(324) characteristic coefficient F is determinedcValue, the point cloud cluster in first defining point cloud in a certain vertex neighborhood, when certain point cloud
Points in cluster are more than the business of the points and the quantity of point cloud cluster in neighborhood, then it is assumed that this cloud cluster is principal point cloud cluster;In order to certainly
The quantity of dynamic acquisition point cloud cluster, sets an angle threshold, when 2 points of normal direction angles are more than threshold value, that is, thinks at this 2 points
In different point cloud clusters;A quantity of cloud cluster is obtained, so as to calculate the corresponding characteristic coefficient F of each pointc, a, 2 points of b's
Normal direction angle is represented by:
D (a, b)=cos-1(na,nb); (8)
(325) F is determinedcValue after, it is adaptive calculate each point apart from bandwidth σd, regression criterion bandwidth σr, normal direction
Difference bandwidth σn, it is expressed as follows:
σd=(1+Fc)·δl (9)
σr=rmax/3Fc (10)
σn=| | ni-n||max/3Fc; (11)
(326) because workpiece inevitably has the sharp features such as edge and turning so that the point in neighborhood is in multiple
On dough sheet, so as to cause larger normal estimation error, therefore the 3rd weighting function is introduced in PCA algorithms, then a most young waiter in a wineshop or an inn
Multiply fitting to be expressed as:
S.t. | | n | |=1 (12)
ω in formulan(ni) --- the weighting function relevant with normal direction difference, it is expressed as:σnFor normal direction
Difference bandwidth, niFor certain initial normal direction of point, gained is solved by formula (6);
(327) calculated by adaptometer in step (325) each put apart from bandwidth σd, regression criterion bandwidth σr, method
To difference bandwidth σnThe weighting function each put in point cloud is calculated, is replaced in covariance matrix C
Then formula (12) is solved by the covariance matrix C after replacement, obtains the optimization normal direction of each point.
When being solved using above-mentioned steps, the weighting function for influenceing each point is apart from bandwidth σ respectivelyd, regression criterion bandwidth σr
With normal direction difference bandwidth σn, above-mentioned steps first according to the shape and the distribution situation of scanning element cloud of workpiece, preset σdAnd σr's
Initial value, initial normal direction is calculated using formula (6), and adaptively pinpoints the bandwidth parameter of each point in cloud really based on initial normal direction, then
Initial normal direction is optimized using formula (7).
Traditional algorithm is the global bandwidth parameter of artificial setting, and it is all that experiment proves that this parameter is not adapted to
Point.The present invention sets different bandwidth parameters to complete normal estimation according to the difference in residing region.By the way of adaptive,
According to the neighborhood characteristics of each point in a cloud, to determine corresponding bandwidth parameter, so as to complete the normal estimation of a cloud.
Specifically, the point cloud normal direction estimated by above-mentioned steps, can both have been pointed in workpiece with ambiguity, i.e. normal direction
Portion, can also point to outer workpiece, so the present invention is adjusted using minimum spanning tree method to normal direction so that point cloud normal direction refers to
To consistent, its principle is by x in the cloud that sets up an officei、xjIt is that both dot products are calculated apart from close 2 points, if ni·nj≈ 1, then
Represent that 2 points of normal orientation is identical, otherwise xiOr xjShould be reverse, therefore, the step (313) specifically includes step:
(331) the maximum point of z coordinate, as root node, adjusts its normal direction and is allowed to and vector (0,0,1) in selected point cloud
Dot product is more than 0, and this will be such that the normal orientation after adjustment points on the outside of workpiece;
(332) weight c=1- is set | ni·nj|, this ensure that the nonnegative weights, and represent two neighbor points when c is smaller
Normal orientation it is more parallel;
(333) traveled through according to weight c size, if father node niWith child node njMeet ni·nj<0, then by njInstead
To;
(334) all nodes are traveled through, normal orientation is carried out to set a threshold value e in Spreading and diffusion, communication process, works as c
During≤e, multiple spot in One Diffusion Process to neighborhood works as c>During e, one point of One Diffusion Process, after the completion of judge whether to meet c≤e, then enter
Row multiple spot spreads, until being diffused into whole point cloud space;
(335) if failing to travel through all points, point cloud is not connected, and partial data forms local isolated island, if do not traveled through
Count number less than point cloud sum an one thousandth, then give up as noise, otherwise search in the data set not traveled through with
Closest point between having traveled through a little, the starting point traveled through as continuation simultaneously goes to step (334).
The present invention obtains its normal direction information by carrying out normal estimation to a cloud, then will believe comprising normal direction information and position
The three-dimensional point of breath is as processing node, then by trajectory planning, and x, y, z, Rx, Ry, Rz information of node are processed in generation, and to sweep
The sequencing retouched will process node and synchronously be sent to robot 1 successively, and robot 1 carries out moving interpolation and inverse kinematics is asked
Solution, generation polishing posture, so as to complete to scan the polishing of workpiece.
Above-mentioned steps are completed, that is, complete the three-dimensional point cloud collection based on apparatus of the present invention, and according only to three-dimensional point cloud
Information generate robot 1 polishing track.
The above embodiment of the present invention is only intended to clearly illustrate example of the present invention, and is not to the present invention
Embodiment restriction.For those of ordinary skill in the field, it can also make on the basis of the above description
Other various forms of changes or variation.There is no necessity and possibility to exhaust all the enbodiments.It is all the present invention
Any modifications, equivalent substitutions and improvements made within spirit and principle etc., should be included in the protection of the claims in the present invention
Within the scope of.
Claims (10)
1. a kind of real-time acquisition device of three-dimensional point cloud, it is characterised in that:Including robot, laser displacement sensor, based on reality
When operating system robot real-time control system, the laser displacement sensor is arranged on the end of the robot by fixture
End, the artificial six degree of freedom articulated robot of described machine;The robot real-time control system passes through real-time industrial ethernet
Bus connects robot and laser displacement sensor, for making the reading of laser displacement sensor and the pose of robot obtain together
Step, three-dimensional measurement is expanded to by one-dimensional measurement, so as to obtain the three-dimensional point cloud of workpiece by scanning.
2. the real-time acquisition device of three-dimensional point cloud according to claim 1, it is characterised in that:Described robot is controlled in real time
System processed between robot and laser displacement sensor for Millisecond with communicating.
3. a kind of robot polishing method for planning track based on the described device of claim 1 or 2, it is characterised in that including step
Suddenly:
(1) measurement coordinate system the turning relative to robot end's coordinate system of multi-point calibration method Calibration of Laser displacement transducer is passed through
Change matrix;
(2) robotically-driven laser displacement sensor scanning workpiece, expands to three-dimensional measurement by one-dimensional measurement, obtains the three of workpiece
Dimension point cloud;
(3) by setting three-dimensional point cloud described in different adaptive weighting function pairs to carry out normal estimation to each point in three-dimensional point cloud
And optimization, while the spatial positional information included based on three-dimensional point cloud and normal direction information planning robot polishing when pose and walk
Cutter track footpath, and corrected continuous polishing track as robot to scan track.
The method for planning track 4. robot according to claim 3 polishes, it is characterised in that the step (1) is specifically wrapped
Include:
(11) top in robot end's one demarcation of installation first, control machine people reaches a certain fixing point in space, using machine
5 standardizations of device people tool coordinates system, the fixing point is read relative to robot after the completion of tool coordinates system demarcation from teaching machine
The locus of basis coordinates systemBp;
(12) and then in robot end's clamping laser displacement sensor, laser is radiated at the fixing point by control machine people,
Read pose of the robot current end relative to basis coordinates systemAnd the registration of laser displacement sensorMP, the laser displacement
The registration of sensorMP and the locusBP transformational relation is:
<mrow>
<mmultiscripts>
<mi>p</mi>
<mi>B</mi>
</mmultiscripts>
<mo>=</mo>
<mmultiscripts>
<mi>T</mi>
<mn>6</mn>
<mn>1</mn>
</mmultiscripts>
<msup>
<mmultiscripts>
<mi>T</mi>
<mi>M</mi>
<mn>6</mn>
</mmultiscripts>
<mi>M</mi>
</msup>
<mi>p</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
In formulaBP --- coordinates of the spatial point p in robot basis coordinates system;
MCoordinates of p --- the spatial point p in measurement coordinate system;
--- homogeneous transform matrix of robot end's coordinate system relative to basis coordinates system;
--- homogeneous transform matrix of the coordinate system relative to robot end's coordinate system is measured, wherein,
<mrow>
<mmultiscripts>
<mi>T</mi>
<mi>M</mi>
<mn>6</mn>
</mmultiscripts>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>13</mn>
</msub>
</mtd>
<mtd>
<mi>x</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>23</mn>
</msub>
</mtd>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>33</mn>
</msub>
</mtd>
<mtd>
<mi>z</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
(13) formula (1) is converted toIfThen (1) can be expressed as follows:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mn>13</mn>
</msub>
<mi>d</mi>
<mo>+</mo>
<mi>x</mi>
<mo>=</mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mn>23</mn>
</msub>
<mi>d</mi>
<mo>+</mo>
<mi>y</mi>
<mo>=</mo>
<mi>b</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mn>33</mn>
</msub>
<mi>d</mi>
<mo>+</mo>
<mi>z</mi>
<mo>=</mo>
<mi>c</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
1
(14) multi-point calibration is based on, formula (3) is changed into AX=k form, an over-determined systems are constituted, and use a most young waiter in a wineshop or an inn
Multiply to solve, form is as follows:
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>|</mo>
<mo>|</mo>
<mi>A</mi>
<mi>X</mi>
<mo>-</mo>
<mi>k</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
MakeMinimum solution X*As over-determined systems AX=k least square solution;
(15) solve and matrix is obtained after XObtain matrixAfterwards, you can so that the workpiece three-dimensional point cloud of scanning to be transformed into
Under robot basis coordinates system.
The method for planning track 5. robot according to claim 3 polishes, it is characterised in that described by one in step (2)
The step of dimension measurement expands to three-dimensional measurement is that y location is by robot by robot real-time control system, the x of workpiece point cloud
End pose understands that the z location of workpiece point cloud passes through real-time bus mechanics of communication from laser displacement sensor, it is ensured that x,
Y, z synchronization, three-dimensional measurement is expanded to so as to realize by one-dimensional measurement.
The method for planning track 6. robot according to claim 3 polishes, it is characterised in that described step (3) is specific
Including step:
(31) spatial positional information of each point based on three-dimensional point cloud and after being optimized and revised by adaptive weighting function
Normal direction information planning robot polishing pose;
(32) the processing node of defining point cloud, and processing nodal information is sent to robot, machine successively to scan sequencing
Device people, which carries out inverse kinematics and moving interpolation, will process the continuous robot polishing track of node series connection generation one.
The method for planning track 7. robot according to claim 6 polishes, it is characterised in that described step (31) is specific
Including step:
(301) determine polishing when cutter position, i.e., using each point in three-dimensional point cloud included based on robot basis coordinates system
Spatial positional information;Position when make it that cutter reaches polishing by Robotic inverse kinematics control machine people;
(302) posture of cutter during polishing is determined, i.e., the estimating to each point of three-dimensional point cloud by adaptive weighting function
The posture of cutter when normal direction information and normal direction information planning robot polishing after optimizing and revising.
The method for planning track 8. robot according to claim 7 polishes, it is characterised in that described step (302) tool
Body includes step:
(311) k closest points are searched for by improved k Neighborhood-region-search algorithms, i.e., using a point as body-centered, with fixation
Step-length builds the square of itself, and k neighborhood search is then carried out in square, if finding k closest points, that is, terminates
Search, otherwise expands the scope of square with fixed step size, k closest points is re-searched for, until being searched in square
Rope is to k closest point;
(312) improved PCA algorithms are used, different weighting functions are assigned by the point in the k neighborhoods to certain point, using minimum
Two, which multiply fitting process, fits a plane, using the normal direction of institute's fit Plane as the point initial normal direction;It is next based on described first
Beginning normal direction adaptively pinpoints the weighting function of each point in cloud really, then initial normal direction is optimized;
(313) point cloud is redirected, i.e., the normal direction after optimization is adjusted using minimum spanning tree method so that point cloud normal direction is pointed to
Unanimously.
The method for planning track 9. robot according to claim 8 polishes, it is characterised in that described step (312) tool
Body includes step:
(321) two weighting functions are introduced into PCA algorithms causes approximating method has to the exterior point and noise in a cloud certain
Robustness, least square fitting is represented by:
<mrow>
<mo>{</mo>
<msup>
<mi>d</mi>
<mi>m</mi>
</msup>
<mo>,</mo>
<msup>
<mi>n</mi>
<mi>m</mi>
</msup>
<mo>}</mo>
<mo>=</mo>
<mi>arg</mi>
<mi> </mi>
<msubsup>
<mi>min&Sigma;r</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
S.t. | | n | |=1 (5)
ω in formulad(xi) --- Gauss weighting function so that effect of the point near apart from the point to fit Plane is a little big, on the contrary
Influence is small, ωd(xi) be expressed as:σdFor apart from bandwidth, its initial value is according to the shape of workpiece and sweeps
The distribution situation of described point cloud is preset;
ωr(ri) --- the Gauss weighting function relevant with regression criterion, ωr(ri) be expressed as:σrFor fitting
Residual error bandwidth, its initial value is preset according to the shape of workpiece and the distribution situation of scanning element cloud;
ri m--- the m times iteration, xiThe regression criterion of point, ri mIt is represented by:ri m=dm+(xi-x)Tnm;(322) formula (5) is solved
Covariance matrix C is constructed, is expressed as follows:
<mrow>
<mi>C</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula--- the center of neighborhood is represented, is expressed as:
K --- the points in neighborhood are represented,
By solving covariance matrix C, its corresponding characteristic vector of minimum characteristic value is the normal direction n of fit Plane, that is, is made
For the initial normal direction of the point, and
(323) a characteristic coefficient F is defined for each point in point cloudcIt is to be located at smooth region, sharp features region to quantify the point
Or transitional region, characteristic coefficient is represented by:
S=2 in formula, represents when ratio is 0.5, then to judge that the point is located at sharp features region;
δl--- represent the average distance of the point 6 point closest with it;
(324) characteristic coefficient F is determinedcValue, the point cloud cluster in first defining point cloud in a certain vertex neighborhood, when in the cloud cluster of certain point
Points be more than neighborhood in points with point cloud cluster quantity business, then it is assumed that this cloud cluster be principal point cloud cluster;In order to automatic
The quantity of point cloud cluster is obtained, an angle threshold is set, when 2 points of normal direction angles are more than threshold value, that is, thinks that be located at this 2 points
In different point cloud clusters;A quantity of cloud cluster is obtained, so as to calculate the corresponding characteristic coefficient F of each pointc, a, the normal direction that 2 points of b
Angle is represented by:
D (a, b)=cos-1(na,nb); (8)
(325) F is determinedcValue after, it is adaptive calculate each point apart from bandwidth σd, regression criterion bandwidth σr, normal direction difference
Bandwidth σn, it is expressed as follows:
σd=(1+Fc)·δl (9)
σr=rmax/3Fc (10)
σn=| | ni-n||max/3Fc; (11)
(326) the 3rd weighting function is introduced in PCA algorithms, then least square fitting is expressed as:
<mrow>
<mo>{</mo>
<msup>
<mi>d</mi>
<mi>k</mi>
</msup>
<mo>,</mo>
<msup>
<mi>n</mi>
<mi>k</mi>
</msup>
<mo>}</mo>
<mo>=</mo>
<mi>arg</mi>
<mi> </mi>
<msubsup>
<mi>min&Sigma;r</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>n</mi>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
S.t. | | n | |=1 (12)
ω in formulan(ni) --- the weighting function relevant with normal direction difference, it is expressed as:σnFor normal direction difference
Bandwidth, niFor certain initial normal direction of point, gained is solved by formula (6);
(327) calculated by adaptometer in step (325) each put apart from bandwidth σd, regression criterion bandwidth σr, normal direction it is poor
Different bandwidth σnThe weighting function each put in point cloud is calculated, is replaced in covariance matrix C
<mrow>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Sigma;&omega;</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Sigma;&omega;</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Then formula (12) is solved by the covariance matrix C after replacement, obtains the optimization normal direction of each point.
The method for planning track 10. robot according to claim 8 polishes, it is characterised in that the step (313) is specific
Including step:
(331) the maximum point of z coordinate, as root node, adjusts its normal direction and is allowed to dot product with vector (0,0,1) in selected point cloud
More than 0, this will be such that the normal orientation after adjustment points on the outside of workpiece;
(332) weight c=1- is set | ni·nj|, this ensure that the nonnegative weights, and when the method for c two neighbor points of smaller expression
It is more parallel to direction;
(333) traveled through according to weight c size, if father node niWith child node njMeet ni·nj<0, then by njReversely;
(334) all nodes are traveled through, normal orientation are carried out to set a threshold value e in Spreading and diffusion, communication process, as c≤e
When, multiple spot in One Diffusion Process to neighborhood works as c>During e, one point of One Diffusion Process, after the completion of judge whether to meet c≤e, then carry out
Multiple spot spreads, until being diffused into whole point cloud space;
(335) if failing to travel through all points, point cloud is not connected, and partial data forms local isolated island, if not traveling through points
Number is then given up less than an one thousandth for point cloud sum as noise, otherwise search in the data set not traveled through with time
Closest point between going through a little, the starting point traveled through as continuation simultaneously goes to step (334).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710334080.8A CN107127755B (en) | 2017-05-12 | 2017-05-12 | Real-time acquisition device of three-dimensional point cloud and robot polishing track planning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710334080.8A CN107127755B (en) | 2017-05-12 | 2017-05-12 | Real-time acquisition device of three-dimensional point cloud and robot polishing track planning method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107127755A true CN107127755A (en) | 2017-09-05 |
CN107127755B CN107127755B (en) | 2023-12-08 |
Family
ID=59731676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710334080.8A Active CN107127755B (en) | 2017-05-12 | 2017-05-12 | Real-time acquisition device of three-dimensional point cloud and robot polishing track planning method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107127755B (en) |
Cited By (49)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107932505A (en) * | 2017-11-15 | 2018-04-20 | 广东工业大学 | Optimal polishing task path planning method and device based on articulated arm robots |
CN108582076A (en) * | 2018-05-10 | 2018-09-28 | 武汉库柏特科技有限公司 | A kind of Robotic Hand-Eye Calibration method and device based on standard ball |
CN108662989A (en) * | 2018-07-02 | 2018-10-16 | 郑州工程技术学院 | A kind of car light profile quality determining method based on 3 D laser scanning |
CN109366492A (en) * | 2018-10-24 | 2019-02-22 | 武汉理工大学 | Casting grinding track online compensation system and method based on robot |
CN109461183A (en) * | 2018-10-23 | 2019-03-12 | 沙洲职业工学院 | A kind of method of space coordinate point and point cloud location point Three Dimensional Contrast |
CN109454642A (en) * | 2018-12-27 | 2019-03-12 | 南京埃克里得视觉技术有限公司 | Robot coating track automatic manufacturing method based on 3D vision |
CN109483369A (en) * | 2018-12-13 | 2019-03-19 | 中国船舶重工集团公司第七六研究所 | A kind of robot polishing system and its control method with 3D vision |
CN109514133A (en) * | 2018-11-08 | 2019-03-26 | 东南大学 | A kind of autonomous teaching method of welding robot 3D curved welding seam based on line-structured light perception |
CN109605390A (en) * | 2018-12-28 | 2019-04-12 | 芜湖哈特机器人产业技术研究院有限公司 | A kind of automobile washing machine people's control system |
CN109605140A (en) * | 2018-12-25 | 2019-04-12 | 珞石(山东)智能科技有限公司 | Based on machine vision and have the function of that the cutter of six shaft mechanical arm of power control puts the first edge on a knife or a pair of scissors method |
CN110000793A (en) * | 2019-04-29 | 2019-07-12 | 武汉库柏特科技有限公司 | A kind of motion planning and robot control method, apparatus, storage medium and robot |
WO2019136716A1 (en) * | 2018-01-12 | 2019-07-18 | 浙江国自机器人技术有限公司 | Cleaning method for self-planning route |
CN110091333A (en) * | 2019-05-17 | 2019-08-06 | 上海交通大学 | The device and method of complex-curved surface weld feature identification and automatic grinding and polishing |
CN110103118A (en) * | 2019-06-18 | 2019-08-09 | 苏州大学 | A kind of paths planning method of milling robot, device, system and storage medium |
CN110111424A (en) * | 2019-05-07 | 2019-08-09 | 易思维(杭州)科技有限公司 | The three-dimensional rebuilding method of arc-shaped object based on line-structured light measurement |
CN110253373A (en) * | 2019-07-15 | 2019-09-20 | 北京石油化工学院 | A kind of system and method for view-based access control model guided robot polishing member welding joints |
CN110281152A (en) * | 2019-06-17 | 2019-09-27 | 华中科技大学 | A kind of robot constant force polishing paths planning method and system based on online examination touching |
CN110355755A (en) * | 2018-12-15 | 2019-10-22 | 深圳铭杰医疗科技有限公司 | Robot hand-eye system calibration method, apparatus, equipment and storage medium |
CN110370287A (en) * | 2019-08-16 | 2019-10-25 | 中铁第一勘察设计院集团有限公司 | Subway column inspection robot path planning's system and method for view-based access control model guidance |
CN110421436A (en) * | 2019-08-29 | 2019-11-08 | 四川智能创新铸造有限公司 | The removal system of robot machining steel-casting increasing meat and riser root |
CN110722554A (en) * | 2019-09-02 | 2020-01-24 | 深圳群宾精密工业有限公司 | Manipulator track editing and correcting method based on laser point cloud data |
CN110744553A (en) * | 2019-12-06 | 2020-02-04 | 大连誉洋工业智能有限公司 | Automatic path planning method for 3D vision robot |
CN110842901A (en) * | 2019-11-26 | 2020-02-28 | 广东技术师范大学 | Robot hand-eye calibration method and device based on novel three-dimensional calibration block |
CN111216124A (en) * | 2019-12-02 | 2020-06-02 | 广东技术师范大学 | Robot vision guiding method and device based on integration of global vision and local vision |
CN111299761A (en) * | 2020-02-28 | 2020-06-19 | 华南理工大学 | Real-time attitude estimation method of welding seam tracking system |
CN111546330A (en) * | 2020-04-15 | 2020-08-18 | 浙江娃哈哈智能机器人有限公司 | Automatic calibration method for coordinate system of chemical part |
CN111558870A (en) * | 2020-04-16 | 2020-08-21 | 华中科技大学 | Robot intelligent polishing system and method for composite material component of airplane body |
CN111843505A (en) * | 2020-07-16 | 2020-10-30 | 武汉数字化设计与制造创新中心有限公司 | In-situ measurement-milling and repairing integrated process method and system for field robot |
CN111898219A (en) * | 2020-07-29 | 2020-11-06 | 华中科技大学 | Area division method and equipment for large-scale complex component robotic surface machining |
CN112223294A (en) * | 2020-10-22 | 2021-01-15 | 湖南大学 | Mechanical arm machining track correction method based on three-dimensional vision |
CN112318226A (en) * | 2020-11-02 | 2021-02-05 | 芜湖哈特机器人产业技术研究院有限公司 | Method for polishing surface of circular workpiece |
CN112446952A (en) * | 2020-11-06 | 2021-03-05 | 杭州易现先进科技有限公司 | Three-dimensional point cloud normal vector generation method and device, electronic equipment and storage medium |
CN113231910A (en) * | 2021-04-29 | 2021-08-10 | 武汉中观自动化科技有限公司 | Method and system for acquiring polishing track of sole edge |
CN113400327A (en) * | 2021-07-07 | 2021-09-17 | 天津大学 | Master-slave teleoperation system and method for grinding and cutting integrated machining of medium-large casting parts |
CN113751890A (en) * | 2020-06-03 | 2021-12-07 | 上海发那科机器人有限公司 | Robot curved surface track cutting method and system based on laser displacement sensor |
CN113894785A (en) * | 2021-10-27 | 2022-01-07 | 华中科技大学无锡研究院 | Control method, device and system for in-situ measurement and processing of blades of water turbine |
CN113910258A (en) * | 2021-10-21 | 2022-01-11 | 上海交通大学 | Double-robot measurement-milling integrated machining system and control method thereof |
CN114049331A (en) * | 2021-11-17 | 2022-02-15 | 长春理工大学 | Method for polishing surface of complex workpiece |
CN114055255A (en) * | 2021-11-18 | 2022-02-18 | 中南大学 | Large-scale complex component surface polishing path planning method based on real-time point cloud |
CN114322847A (en) * | 2022-03-15 | 2022-04-12 | 北京精雕科技集团有限公司 | Vectorization method and device for measured data of unidirectional scanning sensor |
CN114488943A (en) * | 2022-01-07 | 2022-05-13 | 华中科技大学 | Random multi-region efficient polishing path planning method oriented to cooperation condition |
CN115026834A (en) * | 2022-07-02 | 2022-09-09 | 埃夫特智能装备股份有限公司 | Method for realizing program deviation rectifying function based on robot template |
CN115107034A (en) * | 2022-07-18 | 2022-09-27 | 江南大学 | Quantitative iterative learning control method for single mechanical arm |
CN115351781A (en) * | 2022-07-20 | 2022-11-18 | 福州大学 | Industrial robot grinding and polishing path generation method and equipment based on solidworks |
CN115592501A (en) * | 2022-10-11 | 2023-01-13 | 中国第一汽车股份有限公司(Cn) | Top cover brazing self-adaptive polishing method based on 3D line laser vision guidance |
CN115656238A (en) * | 2022-10-17 | 2023-01-31 | 中国科学院高能物理研究所 | Micro-area XRF (X-ray fluorescence) elemental analysis and multi-dimensional imaging method and system |
CN117289651A (en) * | 2023-11-24 | 2023-12-26 | 南通汤姆瑞斯工业智能科技有限公司 | Numerical control machining method and control system for die manufacturing |
CN117422763A (en) * | 2023-12-19 | 2024-01-19 | 商飞智能技术有限公司 | Method and device for positioning polishing area and planning polishing track on surface of die |
CN118066987A (en) * | 2024-04-19 | 2024-05-24 | 中国空气动力研究与发展中心超高速空气动力研究所 | Automatic temperature-sensitive paint film thickness measurement electrical control system and electrical control method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1730248A (en) * | 2005-08-20 | 2006-02-08 | 大连海事大学 | Reverse engineering robot system |
JP2011167815A (en) * | 2010-02-19 | 2011-09-01 | Ihi Corp | Object recognizing robot system |
CN103776378A (en) * | 2014-02-27 | 2014-05-07 | 上海思琢自动化科技有限公司 | Non-contact type flexible on-line dimension measurement system |
US20170060132A1 (en) * | 2015-08-31 | 2017-03-02 | Korea University Research And Business Foundation | Method for detecting floor obstacle using laser range finder |
CN106584273A (en) * | 2017-01-31 | 2017-04-26 | 西北工业大学 | Online visual detecting system for robot polishing |
CN206825428U (en) * | 2017-05-12 | 2018-01-02 | 华南理工大学 | A kind of real-time acquisition device of three-dimensional point cloud |
-
2017
- 2017-05-12 CN CN201710334080.8A patent/CN107127755B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1730248A (en) * | 2005-08-20 | 2006-02-08 | 大连海事大学 | Reverse engineering robot system |
JP2011167815A (en) * | 2010-02-19 | 2011-09-01 | Ihi Corp | Object recognizing robot system |
CN103776378A (en) * | 2014-02-27 | 2014-05-07 | 上海思琢自动化科技有限公司 | Non-contact type flexible on-line dimension measurement system |
US20170060132A1 (en) * | 2015-08-31 | 2017-03-02 | Korea University Research And Business Foundation | Method for detecting floor obstacle using laser range finder |
CN106584273A (en) * | 2017-01-31 | 2017-04-26 | 西北工业大学 | Online visual detecting system for robot polishing |
CN206825428U (en) * | 2017-05-12 | 2018-01-02 | 华南理工大学 | A kind of real-time acquisition device of three-dimensional point cloud |
Cited By (70)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107932505A (en) * | 2017-11-15 | 2018-04-20 | 广东工业大学 | Optimal polishing task path planning method and device based on articulated arm robots |
CN107932505B (en) * | 2017-11-15 | 2021-06-08 | 广东工业大学 | Optimal polishing task path planning method and device based on articulated arm robot |
WO2019136716A1 (en) * | 2018-01-12 | 2019-07-18 | 浙江国自机器人技术有限公司 | Cleaning method for self-planning route |
CN108582076A (en) * | 2018-05-10 | 2018-09-28 | 武汉库柏特科技有限公司 | A kind of Robotic Hand-Eye Calibration method and device based on standard ball |
CN108662989A (en) * | 2018-07-02 | 2018-10-16 | 郑州工程技术学院 | A kind of car light profile quality determining method based on 3 D laser scanning |
CN109461183A (en) * | 2018-10-23 | 2019-03-12 | 沙洲职业工学院 | A kind of method of space coordinate point and point cloud location point Three Dimensional Contrast |
CN109366492A (en) * | 2018-10-24 | 2019-02-22 | 武汉理工大学 | Casting grinding track online compensation system and method based on robot |
CN109514133A (en) * | 2018-11-08 | 2019-03-26 | 东南大学 | A kind of autonomous teaching method of welding robot 3D curved welding seam based on line-structured light perception |
CN109483369B (en) * | 2018-12-13 | 2023-09-29 | 中国船舶集团有限公司第七一六研究所 | Robot polishing system with three-dimensional vision and control method thereof |
CN109483369A (en) * | 2018-12-13 | 2019-03-19 | 中国船舶重工集团公司第七六研究所 | A kind of robot polishing system and its control method with 3D vision |
CN110355755A (en) * | 2018-12-15 | 2019-10-22 | 深圳铭杰医疗科技有限公司 | Robot hand-eye system calibration method, apparatus, equipment and storage medium |
CN109605140A (en) * | 2018-12-25 | 2019-04-12 | 珞石(山东)智能科技有限公司 | Based on machine vision and have the function of that the cutter of six shaft mechanical arm of power control puts the first edge on a knife or a pair of scissors method |
CN109454642A (en) * | 2018-12-27 | 2019-03-12 | 南京埃克里得视觉技术有限公司 | Robot coating track automatic manufacturing method based on 3D vision |
CN109454642B (en) * | 2018-12-27 | 2021-08-17 | 南京埃克里得视觉技术有限公司 | Robot gluing track automatic production method based on three-dimensional vision |
CN109605390A (en) * | 2018-12-28 | 2019-04-12 | 芜湖哈特机器人产业技术研究院有限公司 | A kind of automobile washing machine people's control system |
CN110000793A (en) * | 2019-04-29 | 2019-07-12 | 武汉库柏特科技有限公司 | A kind of motion planning and robot control method, apparatus, storage medium and robot |
CN110111424B (en) * | 2019-05-07 | 2023-06-06 | 易思维(杭州)科技有限公司 | Three-dimensional reconstruction method of arc-shaped object based on line structured light measurement |
CN110111424A (en) * | 2019-05-07 | 2019-08-09 | 易思维(杭州)科技有限公司 | The three-dimensional rebuilding method of arc-shaped object based on line-structured light measurement |
CN110091333A (en) * | 2019-05-17 | 2019-08-06 | 上海交通大学 | The device and method of complex-curved surface weld feature identification and automatic grinding and polishing |
CN110091333B (en) * | 2019-05-17 | 2022-05-06 | 上海交通大学 | Device and method for identifying and automatically grinding and polishing weld joint features on surface of complex curved surface |
CN110281152A (en) * | 2019-06-17 | 2019-09-27 | 华中科技大学 | A kind of robot constant force polishing paths planning method and system based on online examination touching |
CN110103118A (en) * | 2019-06-18 | 2019-08-09 | 苏州大学 | A kind of paths planning method of milling robot, device, system and storage medium |
CN110253373A (en) * | 2019-07-15 | 2019-09-20 | 北京石油化工学院 | A kind of system and method for view-based access control model guided robot polishing member welding joints |
CN110370287B (en) * | 2019-08-16 | 2022-09-06 | 中铁第一勘察设计院集团有限公司 | Subway train inspection robot path planning system and method based on visual guidance |
CN110370287A (en) * | 2019-08-16 | 2019-10-25 | 中铁第一勘察设计院集团有限公司 | Subway column inspection robot path planning's system and method for view-based access control model guidance |
CN110421436A (en) * | 2019-08-29 | 2019-11-08 | 四川智能创新铸造有限公司 | The removal system of robot machining steel-casting increasing meat and riser root |
CN110722554A (en) * | 2019-09-02 | 2020-01-24 | 深圳群宾精密工业有限公司 | Manipulator track editing and correcting method based on laser point cloud data |
CN110842901A (en) * | 2019-11-26 | 2020-02-28 | 广东技术师范大学 | Robot hand-eye calibration method and device based on novel three-dimensional calibration block |
WO2021109575A1 (en) * | 2019-12-02 | 2021-06-10 | 广东技术师范大学 | Global vision and local vision integrated robot vision guidance method and device |
CN111216124A (en) * | 2019-12-02 | 2020-06-02 | 广东技术师范大学 | Robot vision guiding method and device based on integration of global vision and local vision |
CN111216124B (en) * | 2019-12-02 | 2020-11-06 | 广东技术师范大学 | Robot vision guiding method and device based on integration of global vision and local vision |
JP7212236B2 (en) | 2019-12-02 | 2023-01-25 | 広東技術師範大学 | Robot Visual Guidance Method and Apparatus by Integrating Overview Vision and Local Vision |
JP2022516852A (en) * | 2019-12-02 | 2022-03-03 | 広東技術師範大学 | Robot visual guidance method and device by integrating overview vision and local vision |
CN110744553A (en) * | 2019-12-06 | 2020-02-04 | 大连誉洋工业智能有限公司 | Automatic path planning method for 3D vision robot |
CN111299761A (en) * | 2020-02-28 | 2020-06-19 | 华南理工大学 | Real-time attitude estimation method of welding seam tracking system |
CN111546330A (en) * | 2020-04-15 | 2020-08-18 | 浙江娃哈哈智能机器人有限公司 | Automatic calibration method for coordinate system of chemical part |
CN111546330B (en) * | 2020-04-15 | 2022-04-19 | 浙江娃哈哈智能机器人有限公司 | Automatic calibration method for coordinate system of chemical part |
CN111558870B (en) * | 2020-04-16 | 2022-04-15 | 华中科技大学 | Robot intelligent polishing system and method for composite material component of airplane body |
CN111558870A (en) * | 2020-04-16 | 2020-08-21 | 华中科技大学 | Robot intelligent polishing system and method for composite material component of airplane body |
CN113751890B (en) * | 2020-06-03 | 2024-01-23 | 上海发那科机器人有限公司 | Robot curved surface track cutting method and cutting system based on laser displacement sensor |
CN113751890A (en) * | 2020-06-03 | 2021-12-07 | 上海发那科机器人有限公司 | Robot curved surface track cutting method and system based on laser displacement sensor |
CN111843505B (en) * | 2020-07-16 | 2022-04-01 | 武汉数字化设计与制造创新中心有限公司 | In-situ measurement-milling and repairing integrated process method and system for field robot |
CN111843505A (en) * | 2020-07-16 | 2020-10-30 | 武汉数字化设计与制造创新中心有限公司 | In-situ measurement-milling and repairing integrated process method and system for field robot |
CN111898219A (en) * | 2020-07-29 | 2020-11-06 | 华中科技大学 | Area division method and equipment for large-scale complex component robotic surface machining |
CN111898219B (en) * | 2020-07-29 | 2022-04-12 | 华中科技大学 | Area division method and equipment for large-scale complex component robotic surface machining |
CN112223294A (en) * | 2020-10-22 | 2021-01-15 | 湖南大学 | Mechanical arm machining track correction method based on three-dimensional vision |
CN112318226A (en) * | 2020-11-02 | 2021-02-05 | 芜湖哈特机器人产业技术研究院有限公司 | Method for polishing surface of circular workpiece |
CN112446952B (en) * | 2020-11-06 | 2024-01-26 | 杭州易现先进科技有限公司 | Three-dimensional point cloud normal vector generation method and device, electronic equipment and storage medium |
CN112446952A (en) * | 2020-11-06 | 2021-03-05 | 杭州易现先进科技有限公司 | Three-dimensional point cloud normal vector generation method and device, electronic equipment and storage medium |
CN113231910A (en) * | 2021-04-29 | 2021-08-10 | 武汉中观自动化科技有限公司 | Method and system for acquiring polishing track of sole edge |
CN113400327A (en) * | 2021-07-07 | 2021-09-17 | 天津大学 | Master-slave teleoperation system and method for grinding and cutting integrated machining of medium-large casting parts |
CN113910258A (en) * | 2021-10-21 | 2022-01-11 | 上海交通大学 | Double-robot measurement-milling integrated machining system and control method thereof |
CN113894785A (en) * | 2021-10-27 | 2022-01-07 | 华中科技大学无锡研究院 | Control method, device and system for in-situ measurement and processing of blades of water turbine |
CN113894785B (en) * | 2021-10-27 | 2023-06-09 | 华中科技大学无锡研究院 | Control method, device and system for in-situ measurement and processing of turbine blades |
CN114049331A (en) * | 2021-11-17 | 2022-02-15 | 长春理工大学 | Method for polishing surface of complex workpiece |
CN114055255A (en) * | 2021-11-18 | 2022-02-18 | 中南大学 | Large-scale complex component surface polishing path planning method based on real-time point cloud |
CN114488943A (en) * | 2022-01-07 | 2022-05-13 | 华中科技大学 | Random multi-region efficient polishing path planning method oriented to cooperation condition |
CN114488943B (en) * | 2022-01-07 | 2024-01-12 | 华中科技大学 | Random multi-area efficient polishing path planning method oriented to matched working conditions |
CN114322847A (en) * | 2022-03-15 | 2022-04-12 | 北京精雕科技集团有限公司 | Vectorization method and device for measured data of unidirectional scanning sensor |
CN115026834A (en) * | 2022-07-02 | 2022-09-09 | 埃夫特智能装备股份有限公司 | Method for realizing program deviation rectifying function based on robot template |
CN115107034A (en) * | 2022-07-18 | 2022-09-27 | 江南大学 | Quantitative iterative learning control method for single mechanical arm |
CN115351781A (en) * | 2022-07-20 | 2022-11-18 | 福州大学 | Industrial robot grinding and polishing path generation method and equipment based on solidworks |
CN115351781B (en) * | 2022-07-20 | 2024-06-07 | 福州大学 | Method and equipment for generating grinding and polishing path of industrial robot based on solidworks |
CN115592501A (en) * | 2022-10-11 | 2023-01-13 | 中国第一汽车股份有限公司(Cn) | Top cover brazing self-adaptive polishing method based on 3D line laser vision guidance |
CN115656238A (en) * | 2022-10-17 | 2023-01-31 | 中国科学院高能物理研究所 | Micro-area XRF (X-ray fluorescence) elemental analysis and multi-dimensional imaging method and system |
CN117289651B (en) * | 2023-11-24 | 2024-04-16 | 南通汤姆瑞斯工业智能科技有限公司 | Numerical control machining method and control system for die manufacturing |
CN117289651A (en) * | 2023-11-24 | 2023-12-26 | 南通汤姆瑞斯工业智能科技有限公司 | Numerical control machining method and control system for die manufacturing |
CN117422763A (en) * | 2023-12-19 | 2024-01-19 | 商飞智能技术有限公司 | Method and device for positioning polishing area and planning polishing track on surface of die |
CN117422763B (en) * | 2023-12-19 | 2024-05-31 | 商飞智能技术有限公司 | Method and device for positioning polishing area and planning polishing track on surface of die |
CN118066987A (en) * | 2024-04-19 | 2024-05-24 | 中国空气动力研究与发展中心超高速空气动力研究所 | Automatic temperature-sensitive paint film thickness measurement electrical control system and electrical control method |
Also Published As
Publication number | Publication date |
---|---|
CN107127755B (en) | 2023-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107127755A (en) | A kind of real-time acquisition device and robot polishing method for planning track of three-dimensional point cloud | |
CN206825428U (en) | A kind of real-time acquisition device of three-dimensional point cloud | |
CN102494657B (en) | Measuring head radius compensation method for curve surface profile measuring and detecting | |
CN114011608B (en) | Spraying process optimization system based on digital twinning and spraying optimization method thereof | |
CN106248035A (en) | The method and system that a kind of surface profile based on point cloud model accurately detects | |
CN110059879B (en) | Automatic planning method for three-coordinate measurement of vehicle body | |
CN101497279A (en) | Measuring and machining integrated laser three-dimensional marking method and device | |
CN109683552B (en) | Numerical control machining path generation method on complex point cloud model guided by base curve | |
CN104515478A (en) | Automatic three-dimensional measuring method and automatic three-dimensional measuring system for high-precision blade of aviation engine | |
CN108827155A (en) | A kind of robot vision measuring system and method | |
CN102589437A (en) | Calibration method for measuring head center position in light pen-type portable three-coordinate measuring system | |
CN114055255B (en) | Large-scale complex component surface polishing path planning method based on real-time point cloud | |
CN106091931A (en) | A kind of adaptive scanning based on threedimensional model measures system and control method thereof | |
CN103759635A (en) | Scanning measurement robot detection method allowing precision to be irrelevant to robot | |
CN103777570A (en) | Machining error rapid detection and compensation method based on NURBS curved surface | |
CN103434609A (en) | Automatic marking device for ship hull section outer plate | |
CN108994830A (en) | System calibrating method for milling robot off-line programing | |
CN109163675A (en) | A method of angle swing shaft position precision is detected based on laser tracker | |
CN110370316A (en) | It is a kind of based on the robot TCP scaling method vertically reflected | |
CN109883336A (en) | Measuring system and measurement method during a kind of sheet fabrication towards ship surface | |
Vásquez-Gómez et al. | View planning for 3D object reconstruction | |
CN111085902B (en) | Workpiece polishing system for visual online detection and correction | |
CN103236043A (en) | Plant organ point cloud restoration method | |
CN114523475A (en) | Automatic error calibration and compensation device and method for robot assembly system | |
CN104615880A (en) | Rapid ICP (inductively coupled plasma) method for point cloud matching of three-dimensional laser radar |
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 |