CN1985773A - Celebral operating robot system based on optical tracking and closed-loop control and its realizing method - Google Patents

Celebral operating robot system based on optical tracking and closed-loop control and its realizing method Download PDF

Info

Publication number
CN1985773A
CN1985773A CNA200510122586XA CN200510122586A CN1985773A CN 1985773 A CN1985773 A CN 1985773A CN A200510122586X A CNA200510122586X A CN A200510122586XA CN 200510122586 A CN200510122586 A CN 200510122586A CN 1985773 A CN1985773 A CN 1985773A
Authority
CN
China
Prior art keywords
robot
freedom
degrees
computer
pose
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
CNA200510122586XA
Other languages
Chinese (zh)
Other versions
CN100464720C (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.)
Huazhi Computer Application Technology Co Ltd Tianjin
Tsinghua University
Original Assignee
Huazhi Computer Application Technology Co Ltd Tianjin
Tsinghua University
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 Huazhi Computer Application Technology Co Ltd Tianjin, Tsinghua University filed Critical Huazhi Computer Application Technology Co Ltd Tianjin
Priority to CNB200510122586XA priority Critical patent/CN100464720C/en
Publication of CN1985773A publication Critical patent/CN1985773A/en
Application granted granted Critical
Publication of CN100464720C publication Critical patent/CN100464720C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The cerebral operating robot system based on optical tracking and closed-loop control for receiving medical image information, detecting and determining disease focus position, assisting operation plan and performing operation guidance is one closed-loop pose measuring and real-time feedback controlling robot system comprising one computer, one five-freedom robot, one optical tracking device and one passive marker. The passive marker is mounted in the end of the five-freedom robot, and the five-freedom robot includes one mechanical arm and one mechanical arm controller. The operation plan assisting and operation guiding software includes a digital image inputting and pre-treating module, a disease focus extracting and 3D reconstructing module, an operation planning module, and an operation executing module. The present invention has the advantages of raised precision, low cost, simplified calibration process, etc.

Description

Based on optical tracking closed loop control celebral operating robot system and implementation method
Technical field
The invention belongs to the robotics field, be particularly related to a kind ofly, specifically be meant a kind of high precision machines robot system and its implementation that can assist the doctor to carry out meticulous department of cerebral surgery Minimally Invasive Surgery based on optical tracking closed loop control celebral operating robot system and implementation method.
Background technology
The department of cerebral surgery micro-wound operation robot is a kind of robot system that can assist the doctor to implement frameless stereotactic surgery, generally all includes to assist undergo surgery planning and the computer of guiding and the robot that carries out the operation of stereotaxis assisted surgery.Wherein, auxiliary planning and the computer of guiding of undergoing surgery, supporting assisted surgery planning and guiding software are installed, mainly finish the processing and the three-dimensionalreconstruction of patient's brain medical image information, auxiliary doctor carries out virtual operation planning, and sends movement instruction with the control robot motion to robot controller in operation process.This software generally includes digitized video input and functional modules such as pretreatment, focus extraction and three-dimensionalreconstruction, surgery planning and operation enforcement.Digitized video input is the medical image data that receives patient's brain of outside input with the effect of pretreatment module, and converts the form that computer can show to, on this basis image is carried out pretreatment, comprise image reinforcement, remove dry, Histogram statistics etc.; Focus is extracted and the effect of three-dimensionalreconstruction module is to allow the doctor focus be cut apart on image with alternant way, and extracts the head profile, thereby generates the 3 D medical model of focus and head; In the surgery planning module, the doctor can carry out the sign work of target spot and gauge point, and the path planning that can undergo surgery in reality environment; Implement in the module in operation, can send movement instruction with the control robot motion to robot controller according to the good operation pathway of planning in advance in good time, the virtual reality function can provide real-time monitoring for operation simultaneously.
The robot that system adopted is made up of mechanical arm and controller, generally has 5 joint freedom degrees, and all there is driving mechanism in each joint, is used to drive joint motions, can guarantee operating theater instruments with any attitude arrive in the operative space more arbitrarily.Robot controller can be accepted the movement instruction that computer sends, and each joint motions of real time control machine device people arrive assigned address.
The greatest problem that existing celebral operating robot exists is the high precision requirement that positioning accuracy can not satisfy meticulous cerebral surgery operation.Robotic surgical device not only requires the repeatable accuracy height, the absolute precision of apparatus work is required stricter because the surgery planning that carries out in the virtual 3 D medical model space finally will be realized by the motion of robot.Because the factor of mechanism and departure, there are error in the position command of expectation and robot between the actual position that reaches inevitably.This error may reach several millimeters usually to about one centimetre, hinders operating required precision considerably beyond Wicresoft, must manage to be overcome.In commercial Application, the absolute fix error of robot can overcome by the demarcation of variety of way, but robotic surgical device is as a kind of surgical apparatus, and its applicable cases difference is very big, does not generally have strict fixed installation position.Obviously, to each example operation all requires to carry out strictness, loaded down with trivial details field calibration is unaccommodated, moreover the influence of some non-geometric error factors such as mechanism's elastic deformation, joint space and flexibility etc., also be difficult to rely on and demarcate solution.At present, many work all are the methods of considering to improve precision from machinery, rely on mechanism's processing, assembling and high-precision control and the compensation of high precision to guarantee system accuracy, be difficult to reach ideal effect, and increased the design and the manufacturing cost of robot greatly.Existing celebral operating robot is open loop from the overall situation, promptly provide position command to guided robot motion and location process from computer, lack the feedback and the correction of robot end's pose, various error components can't obtain the closed-loop corrected of the overall situation, and absolute pose accuracy is difficult to guarantee.
In addition, from disclosed report (for example Chinese patent publication number CN1243690A), the so-called mapping localization method of determining the spatial position of the relative robot manipulation of focus is too complicated, need use three-dimensional positioning framework, mapping nail or double template, six joint digitized mechanical arms, power control human-computer interaction technology etc., not only time-consuming but also require great effort, also bring misery and pressure to patient; Also have utilization vision measurement technology to carry out process registration (for example Chinese patent publication number CN1554315A), but precision and environmental suitability also need to improve, difficult practical.
Summary of the invention
The purpose of this invention is to provide a kind ofly, can overcome the shortcoming of prior art based on optical tracking closed loop control celebral operating robot system and implementation method.It is that optictracking device is combined with celebral operating robot, utilize the current location and the attitude of optictracking device real-time tracking robot end apparatus, carry out real-time pose closed loop control in the optical measurement space, make the celebral operating robot end accurately to move according to the desirable operation pathway of planning in advance and to locate, realization can be satisfied the high accuracy celebral operating robot system of meticulous cerebral surgery operation requirement.Simultaneously, by means of the optical tracking technology, simplify the process and the method for the mapping transformation of the mapping 3 D medical model space and robot base coordinate space greatly.Use the present invention can make celebral operating robot reach high accuracy, low cost, and save time, convenient, practical, easy to operate.
What the present invention proposed comprises based on optical tracking closed loop control celebral operating robot system:
Robot with five degrees of freedom, computer (selecting model PC P42.0GHz/256M for use), optictracking device (selecting the model NDI POLARIS of company optictracking device for use), optics registration tool (selecting the passive instrument of the model NDI POLARIS of company standard configuration for use), passive marker (selecting the spherical passive marker of the model NDI POLARIS of company standard configuration for use), (that medically uses a kind ofly makes with the tantalum metal material medical science marker, the opaque pellet shapes index point of correlation electrical height, with the organism compatibility, implant into body for a long time), the hardware components that equipment such as moulding pillow and operation table are formed; Assisted surgery planning and guiding software section promptly receive the position of medical image information, mensuration and definite focus, the auxiliary celebral operating robot that undergos surgery and plan and carry out surgical guide; It is characterized in that:
Said computer, robot with five degrees of freedom, optictracking device and passive marker constitute the robot pose measurement and real-time feedback control system of a closed loop, passive marker is installed in the end of robot with five degrees of freedom, is followed the tracks of by optictracking device all the time; Said robot with five degrees of freedom comprises five degree-of-freedom manipulator and mechanical arm controller; Said assisted surgery planning comprises digitized video input and pretreatment module with guiding software, focus extracts and module is implemented in three-dimensionalreconstruction module, surgery planning module and operation.
Use the undergo surgery concrete operations step in stage of the present invention to be:
One. when before undergoing surgery, preparing: on patient's head, paste four medical science markers, brain is carried out medical image scanning, and will scan gained medical image information input computer, the planning of utilization assisted surgery is determined focus with guiding software, the 3 D medical model of reconstruct focus and head, the sign of target spot and medical science sign image position and planning operation pathway undergo surgery.
Two. registering timing signal: allow patient lie on the operation table, head uses moulding pillow and operation table relative fixed, measure the coordinate of four medical science markers on patient's head on the one hand with the optics registration tool, measured value is provided by optictracking device, and send into computer, by the mapping transformation between the COMPUTER CALCULATION 3 D medical model space and the optical measurement space; With the assignment test method of the present invention's design, cooperate on the other hand, measure the mapping transformation between robot with five degrees of freedom pedestal coordinate space and the optical measurement space automatically by robot with five degrees of freedom, optictracking device, passive marker and computer.
Three. when undergoing surgery enforcement: at first, will transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance, and calculate the desirable pose that need reach at current path point robot with five degrees of freedom end by computer; Then, by the desirable pose of computer according to current path point, carry out by of the coordinate transform of optical measurement space, and, obtain the ideal position in each joint of robot with five degrees of freedom by finding the solution the robot with five degrees of freedom inverse kinematics to robot with five degrees of freedom pedestal coordinate space; Then, import the ideal position in each joint to robot with five degrees of freedom, the motion of control robot with five degrees of freedom by computer; Simultaneously, cooperate by optictracking device and the passive marker that is installed in the robot with five degrees of freedom end, the terminal pose of The real time measure robot with five degrees of freedom is also sent into computer, control method with the present invention's design is controlled in real time to the terminal pose of robot with five degrees of freedom, realizes accurate track following and location; At last, robot with five degrees of freedom locking, doctor's operation that under robot with five degrees of freedom auxiliary, undergos surgery.
Five degree-of-freedom manipulator is made up of arm and wrist, has five joints, takes the PPRRR configuration of two Gliding joints and three cradle heads; Arm segment has three joints, is respectively I, II and III joint, and Gliding joint is adopted in the I joint, and the direction of motion is perpendicular to horizontal plane; Gliding joint is adopted in the II joint, and is "T"-shaped vertical with the I joint; Cradle head is adopted in the III joint, and axis is parallel with I joint motions direction; Wrist partly has two joints, is respectively IV and V joint, and cradle head is adopted in the IV joint, and axis is parallel with the III joints axes; Cradle head is adopted in the V joint, and axis is vertical with the IV joints axes; Terminal apparatus is installed in V joint at mechanical arm, and it is parallel with the IV joints axes that axis is installed.All there is independent driving mechanisms in each joint, is made up of motor, decelerator, and arthrodial driving mechanism also includes ball screw.
The mechanical arm controller is made up of programmable logic controller (PLC) PLC and stepper motor driver, adopt three PLC and five stepper motor driver combinations, five joints of corresponding control five degree-of-freedom manipulator, the mechanical arm controller communicates by RS232C Serial Communication Component and outer computer, can accept the joint position order of outer computer input, the motion of control mechanical arm arrives assigned address.
The POLARIS optictracking device of NDI company is a kind of space measurement location instrument of extensive use in the present medical industry.The POLARIS equipment of passive type comprises the position sensor that can launch and receive infrared illumination, the supporting passive instrument that some passive markers is provided or passive marker is installed; POLARIS equipment just can be determined the position and the direction of instrument in real time by the locus of passive marker on the survey tool, and 3D mean square error scope is usually in the 0.35mm scope; By the RS-232/RS-422 communication, the Data Update frequency of Continuous Tracking reaches 60HZ between POLARIS equipment and the computer.
The method of the robot pose measurement of closed loop and real-time feedback control system is: the terminal pose of being measured robot with five degrees of freedom by optictracking device, calculate controlled quentity controlled variable by computer according to pose that measures and the pose that the expectation robot with five degrees of freedom reaches, the control robot with five degrees of freedom further moves to revise deviation.
The passive marker that optictracking device is installed in the robot with five degrees of freedom end by tracking is measured the terminal pose of robot with five degrees of freedom.
Passive marker is installed in the end of robot with five degrees of freedom, and quantity is at least 3.
Passive marker satisfies at the terminal geometry site of installing of robot with five degrees of freedom: the distance between any two passive markers can not be less than 50mm, minimal spatial separation between any two line segments that are made of passive marker line can not be less than 5mm, any two line segments are not parallel, and the angle between any two line segments can not be lower than 0.5 degree.
Implementation method based on the celebral operating robot system of optical tracking closed loop control, it is characterized in that, carry out pose measurement by means of the optical tracking technology, comprising: the closed loop control method of the terminal pose of the automatic assignment test method of the easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings, robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings and robot with five degrees of freedom.
Described three kinds of methods relate to five coordinate systems, as Fig. 2: set up a 3 D medical model coordinate systems { V} in the 3 D medical model space; On real patient's head, set up a patient coordinate system { P}; In the optical measurement space, set up an optical measurement coordinate system { M}; On the pedestal of robot with five degrees of freedom, set up a robot base coordinate sys-tem { R}; Set up an end-of-arm tooling coordinate system { T} at the end of robot with five degrees of freedom; { M} is an absolute reference system to whole system with the optical measurement coordinate system; Wherein, patient coordinate system is based on that the point at four medical science marker places that stick on patient's head describes, and this coordinate system is chosen the some M at any one marker place in four medical science markers 0As coordinate origin, simultaneously with M 0Some M with other three marker places 1, M 2, M 3Line as three coordinate axess to.
Provide the detailed step and the mathematical description of described three kinds of methods below.
One. the spatial mapping transformation of the 3 D medical model space and optical measurement
Carry out pose measurement by means of the optical tracking technology, the easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings is: paste four medical science markers by the doctor in patient head, these four markers are not or not same plane, and any three markers are not on same straight line; Patient head is carried out CT or the scanning of MRI medical image, gained scan-image input computer; Allow patient lie on the operation table, head uses moulding pillow and operation table relative fixed, and four medical science markers using optics registration tool mensuration patient head are at the spatial coordinate of optical measurement, and measured value is imported computer by optictracking device; By the COMPUTER CALCULATION 3 D medical model space and the spatial mapping transformation of optical measurement.
This mapping transformation is made up of two groups of Coordinate Conversion, promptly from the 3 D medical model coordinate systems V} to patient coordinate system the conversion of P} and from patient coordinate system { { conversion of M}, two groups of conversions all are based on four medical science markers pasting on patient's head and determine P} to the optical measurement coordinate system.
Obviously, { definition of P}, four markers are in that { coordinate among the P} is respectively: M according to patient coordinate system 0(0,0,0), M 1(1,0,0), M 2(0,1,0), M 3(0,0,1).
In addition, because four medical science markers can be discerned in the medical scanning image, { coordinate among the V} also can obtain, and supposes to be designated as: M ' in the 3 D medical model coordinate systems for they 0(x V0, y V0, z V0), M ' 1(x V1, y V1, z V1), M ' 2(x V2, y V2, z V2), M ' 3(x V3, y V3, z V3).
Because the 3 D medical model of patient's brain is by patient's brain scans data reconstruction, therefore can think that { { mapping of P} is rigid transformation (comprising translation, rotation and stretching) to the 3 D medical model coordinate systems, can be with a homogeneous transformation matrix for V} and patient coordinate system VT PFinish two position mappings in the coordinate system.By P} to the homogeneous transformation matrix of V} is:
T P V = x v 1 - x v 0 x v 2 - x v 0 x v 3 - z v 0 x v 0 y v 1 - x v 0 y v 2 - y v 0 y v 3 - z v 0 y v 0 z v 1 - x v 0 z v 2 - y v 0 z v 3 - z v 0 z v 0 0 0 0 1 - - - ( 1 )
In like manner, { P} is to the optical measurement coordinate system { transition matrix of M} can to determine patient coordinate system MT P{ position among the M} can be made as: M by o'clock obtaining to four medical science markers respectively with the optics registration tool in the optical measurement coordinate system for four medical science markers 0(x M0, y M0, z M0), M 1(x M1, y M1, z M1), M 2(x M2, y M2, z M2), M 3(x M3, y M3, z M3).Then by P} to the homogeneous transformation matrix of M} is:
T P M = x m 1 - x m 0 x m 2 - x m 0 x m 3 - x m 0 x m 0 y m 1 - y m 0 y m 2 - y m 0 y m 3 - y m 0 y m 0 z m 1 - z m 0 z m 2 - z m 0 z m 3 - z m 0 z m 0 0 0 0 1 - - - ( 2 )
At last, { { homogeneous transformation of M} is V} to the optical measurement coordinate system by the three dimensional virtual models coordinate system MT VBe:
MT VMT P PT VMT P( VT P) -1 (3)
Two. the spatial mapping transformation of robot with five degrees of freedom pedestal coordinate space and optical measurement
Carry out pose measurement by means of the optical tracking technology, the automatic assignment test method of robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings is: at first, four points in the selected robot with five degrees of freedom work space, these four points should be in the measuring range of optictracking device simultaneously, and four points are coplane not, and any three points are conllinear not; By computer according to default program, send the joint position order that arrives above-mentioned four points to robot with five degrees of freedom, control robot with five degrees of freedom end moves to above-mentioned four points, the coordinate of above-mentioned four points in robot with five degrees of freedom pedestal coordinate space under the computer recording successively; Simultaneously, cooperate by optictracking device and the passive marker that is installed in the robot with five degrees of freedom end, the coordinate of its end in the optical measurement space when sequentially determining robot with five degrees of freedom moves to above-mentioned four points, and send into computer; At last, by computer according to above-mentioned four points at robot with five degrees of freedom pedestal coordinate space and the spatial coordinate figure of optical measurement, calculate the spatial mapping transformation of robot with five degrees of freedom pedestal coordinate space and optical measurement, thereby finish automatic mapping.
Suppose that { coordinate among the R} is described four points in robot base coordinate sys-tem RP i(r i, s i, t i), (i=1,2,3,4), simultaneously, the coordinate of described four points in the optical measurement coordinate system is MP i(u i, v i, w i), (i=1,2,3,4), then by robot base coordinate sys-tem R} to the optical measurement coordinate system homogeneous transformation matrix of M} is:
T R M = u 1 u 2 u 3 u 4 v 1 v 2 v 3 v 4 w 1 w 2 w 3 w 4 1 1 1 1 r 1 r 2 r 3 r 4 s 1 s 2 s 3 s 4 t 1 t 2 t 3 t 4 1 1 1 1 - 1 - - - ( 4 )
Three. the closed loop control method of the terminal pose of robot with five degrees of freedom
Carry out pose measurement by means of the optical tracking technology, the closed loop control method of the terminal pose of robot with five degrees of freedom is: with the optical measurement space is the reference space, at first, to transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance by computer, and calculate the desirable pose that need reach at current path point robot with five degrees of freedom end; Then, by the desirable pose of computer according to current path point, carry out by of the coordinate transform of optical measurement space, and, obtain the ideal position in each joint of robot with five degrees of freedom by finding the solution the robot with five degrees of freedom inverse kinematics to robot with five degrees of freedom pedestal coordinate space; Then, import the ideal position in each joint to robot with five degrees of freedom, the motion of control robot with five degrees of freedom by computer; Simultaneously, cooperate by optictracking device and the passive marker that is installed in the robot with five degrees of freedom end, the pose of The real time measure robot with five degrees of freedom end, and send into computer; Then, the desirable pose of the posture information of The real time measure with the current path point of planning in advance compared, obtain both pose deviations by computer; At last, by computer according to the pose deviation according to predetermined control law calculation correction controlled quentity controlled variable, revise the desirable pose of current path point, and begin the motor control of a new round, thereby realize accurate track following and location.
The control law that is adopted during according to pose deviation calculation correction controlled quentity controlled variable by computer is: Position Control adopts proportional-integral-differential to regulate; Attitude control employing ratio is regulated.
The terminal pose of robot with five degrees of freedom is to carry out closed loop control in the optical measurement space.
As Fig. 5, { T} is at optical measurement coordinate system { the pose X=[P among the M} to establish the end-of-arm tooling coordinate system Tφ T] TExpression.P=[xyz wherein] TExpression { the position of T}; φ=[α β γ] TExpression the attitude of T} is made up of one group of Z-Y-X Eulerian angles,
Note X d=[P d Tφ d T] TRepresent desirable pose, P d=[x dy dz d] T, φ d=[α dβ dγ d] TX r=[P r Tφ r T] TThe pose that the expression actual measurement obtains, P r=[x ry rz r] T, φ r=[α rβ rγ r] TΔ X=[Δ P TΔ φ T] TExpression X dAnd X rDeviation, Δ P=[Δ x Δ y Δ z] T, Δ φ=[Δ α Δ β Δ γ] T,
Position deviation can directly be calculated, that is:
ΔP=P d-P r (5)
The calculating of attitude misalignment relates to rotation transformation, uses R d{ T} is with respect to { the desirable rotation transformation matrix of M} in expression; R rThe real transform matrix that expression obtains according to measured value; Δ R represents by R rTo R dConversion, then:
R d = R Z ( α d ) R Y ( β d ) R X ( γ d )
= cos α d cos β d cos α d sin β d sin γ d - sin α d cos γ d cos α d sin β d cos γ d + sin α d sin γ d sin α d cos β d sin α d sin β d sin γ d + cos α d cos γ d sin α d sin β d cos γ d - cos α d sin γ d - sin β d cos β d sin γ d cos β d cos γ d - - - ( 6 )
R r = R Z ( α r ) R Y ( β r ) R X ( γ r )
= cos α r cos β r cos α r sin β r sin γ r - sin α r cos γ r cos α r sin β r cos γ r + sin α r sin γ r sin α r cos β r sin α r sin β r sin γ r + cos α r cos γ r sin α r sin β r cos γ r - cos α r sin γ r - sin β r cos β r sin γ r cos β r cos γ r - - - ( 7 )
Note: ΔR = R d R r - 1 = r 11 r 12 r 13 r 21 r 22 r 23 r 31 r 32 r 33 - - - ( 8 )
Have again:
ΔR = R Z ( Δα ) R Y ( Δβ ) R X ( Δγ )
= cos Δα cos Δβ cos Δα sin Δβ sin Δγ - sin Δα cos Δγ cos Δα sin Δβ cos Δγ + sin Δα sin Δγ sin Δα cos Δβ sin Δα sin Δβ sin Δγ + cos Δα cos Δγ sin Δα sin Δβ cos Δγ - cos Δα sin Δγ - sin Δβ cos Δβ sin Δγ cos Δβ cos Δγ
(9)
Can solve attitude misalignment Δ φ=[Δ α Δ β Δ γ] by (6), (7), (8), the associating of (9) formula T, wherein:
A. when sin Δ β ≠ 0,
Δβ = A tan 2 ( r 31 2 + r 32 2 , r 33 ) ,
Δα=Atan2(r 23/sinΔβ,r 13/sinΔβ),
Δγ=Atan2(r 32/sinΔβ,-r 31/sinΔβ) (10)
B. when sin Δ β=0, if Δ β=0.0, then:
Δα=0.0,
Δγ=Atan2(-r 12,r 11) (11)
If Δ β=180.0 °, then:
Δα=0.0,
Δγ=Atan2(r 12,-r 11) (12)
Position correction and attitude updating problem are discussed respectively below, are adopted the discrete form row to write formula, the subscript k in the bracket represents to control the circulation k step.
Proportional-integral-differential (PID) control is adopted in position correction, proofreaies and correct controlled quentity controlled variable and is designated as U P=[u xu yu z] T, then:
U P ( k ) = K PP ΔP ( k ) + K PI T Σ j = 0 k ΔP ( j ) + K PD ΔP ( k ) - ΔP ( k - 1 ) T - - - ( 13 )
Wherein T is the controlling of sampling cycle of system, and such as pressing the highest measurement data updating rate 60HZ that POLARIS sets, T can be taken as 16.7ms; K PPThe diagonal matrix of forming by three proportional control factors; K PIThe diagonal matrix of forming by three integral control coefficients; K PDThe diagonal matrix of forming by three derivative control coefficients.Following formula is rewritten into incremental form:
U P ( k ) = U P ( k - 1 ) + ( K PP + T K PI + K PD T ) ΔP ( k )
- ( K PP + 2 K PD T ) ΔP ( k - 1 ) + K PD T ΔP ( k - 2 ) - - - ( 14 )
Note A = K PP + TK PI + K PD T , B = K PP + 2 K PD T , C = K PD T , Then further algorithm is write as computer implemented form:
U P ( k ) = AΔP ( k ) + F ( k - 1 ) F ( k ) = U P ( k ) - BΔP ( k ) + CΔP ( k - 1 ) - - - ( 15 )
The desirable F of initial value (k-1)=0, Δ P (k-1)=0, each step of algorithm all will be calculated Δ P (k), U P(k), F (k), wherein F (k) is used for next step and calculates U P(k).
Ratio (P) control is adopted in attitude updating, proofreaies and correct controlled quentity controlled variable and is designated as U φ=[U αu βu γ] T, corresponding spin matrix is Then:
U φ(k)=K φPΔφ(k) (16)
K wherein φ PThe diagonal matrix of forming by three proportional control factors,
Δ R ^ = R Z ( u α ) R Y ( u β ) R X ( u γ ) - - - ( 17 )
Remember that revised pose instruction is X ^ d = [ P ^ d T φ ^ d T ] T , P ^ d = [ x ^ d y ^ d z ^ d ] T , φ ^ d = [ α ^ d β ^ d γ ^ d ] T , Corresponding spin matrix is
Figure A20051012258600148
Then:
P ^ d ( k ) = P d ( k ) + U p ( k ) - - - ( 18 )
R ^ d ( k ) = R z ( α ^ d ( k ) ) R Y ( β ^ d ( k ) ) R X ( γ ^ d ( k ) ) = R d ( k ) Δ R ^ ( k ) - - - ( 18 )
Unite by (6), (17), (19) formula and to find the solution φ ^ d ( k ) = [ α ^ d ( k ) β ^ d ( k ) γ ^ d ( k ) ] T , Its process can be found the solution Δ φ=[Δ α Δ β Δ γ] with reference to above-mentioned the associating by (6), (7), (8), (9) formula T, repeat no more.
At last,
Figure A200510122586001412
Replace X dInstruction is sent as pose, is used for the guided robot motion.
The advantage that the present invention has is: by robot end's pose is carried out real-time closed loop control, making influences robot motion and localized main error component, mechanism and departure as robot body, and the factors such as various measurements and calculations errors in the operation process can both obtain proofreading and correct, thereby really guaranteed the motion and the positioning accuracy of the surgical robot system overall situation, the track following and the positioning accuracy that have solved existing surgical operation robot existence can not satisfy the problem that clinical delicate procedure utilization requires with imitating.On the other hand, adopted after the overall Pose Control, can relax requirement of mechanism design, manufacturing and the control accuracy of robot body etc., thereby reduce the design and the manufacturing cost of robot body, and the operation motility and the operability of having ready conditions and paying close attention to robot more, design the man-machine operating robot train of mechanism of coordinating more.The method that the present invention proposes is suitable for some industrial robot application strict to absolute precision equally.
The present invention can also demarcate robot base coordinate space and the spatial mapping transformation of optical measurement automatically, removed the artificial calibration process of loaded down with trivial details robot base coordinate sys-tem from, and, just can obtain very quickly and easily by the 3 D medical model space and the spatial mapping transformation of optical measurement by pasting the medical science marker in patient head, adopting the optics registration tool to measure the coordinate of medical science marker.
In a word, one aspect of the present invention can obviously improve the track following and the positioning accuracy of robot system, not only solve the critical absolute precision problem that the restriction surgical operation robot is promoted the use of, and helped reducing the design and the manufacturing cost of robot body; Simplify the registration calibration process on the other hand greatly, not only alleviated doctor's workload, psychological burden and the probability that goes wrong, and alleviated patient's misery and pressure.
Description of drawings
Fig. 1 is a structural representation of the present invention.
Fig. 2 is the sketch map of operation principle block diagram of the present invention and the definition of related coordinate system.
Fig. 3 is the FB(flow block) of the assignment test method of the 3 D medical model space of the present invention and the conversion of optical measurement spatial mappings.
Fig. 4 is the FB(flow block) of the automatic assignment test method of robot with five degrees of freedom pedestal coordinate space of the present invention and the conversion of optical measurement spatial mappings.
Fig. 5 is the flow process and the theory diagram of the closed loop control method of the terminal pose of robot with five degrees of freedom of the present invention.
Fig. 6 is a robot with five degrees of freedom mechanical arm structural representation of the present invention.
The specific embodiment
With reference to accompanying drawing the present invention is elaborated:
As shown in the figure, the implication of a label is among equipment composition and the figure, it is by computer 1, robot with five degrees of freedom 2, optictracking device 3, optics registration tool 4, passive marker 5, medical science marker 6, moulding pillow 7, operation table 8 hardware components, form with assisted surgery planning and guiding software 9 parts, and five degree-of-freedom manipulator 11 and mechanical arm controller 10.
Constitute the robot pose measurement and real-time feedback control system of a closed loop by said computer 1, robot with five degrees of freedom 2, optictracking device 3 and passive marker 5, passive marker 5 is installed in the end of robot with five degrees of freedom 2, is followed the tracks of by optictracking device 3 all the time; Said robot with five degrees of freedom 2 comprises five degree-of-freedom manipulator 11 and mechanical arm controller 10; Said assisted surgery planning comprises digitized video input and pretreatment module with guiding software 9, focus extracts and module is implemented in three-dimensionalreconstruction module, surgery planning module and operation.The quantity of passive marker 5 is at least 3.
The concrete operations step is:
One. when before undergoing surgery, preparing: on patient's head, paste four medical science markers 6, brain is carried out medical image scanning, and will scan gained medical image information input computer 1, the planning of utilization assisted surgery is determined focus with guiding software 9, the 3 D medical model of reconstruct focus and head, the sign of target spot and medical science sign image position and planning operation pathway undergo surgery;
Two. registering timing signal: allow patient lie on the operation table 8, head uses moulding pillow 7 and operation table 8 relative fixed, measure the coordinate of four medical science markers 6 on patient's head on the one hand with optics registration tool 4, measured value is provided by optictracking device 3, and send into computer 1, by the mapping transformation between the computer 1 calculating 3 D medical model space and the optical measurement space; With the assignment test method of the present invention's design, cooperate on the other hand, measure the mapping transformation between robot with five degrees of freedom pedestal coordinate space and the optical measurement space automatically by robot with five degrees of freedom 2, optictracking device 3 and computer 1;
Three. when undergoing surgery enforcement: at first, will transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance, and calculate the desirable pose that need reach at current path point robot with five degrees of freedom end by computer 1; Then, by the desirable pose of computer 1 according to current path point, carry out by of the coordinate transform of optical measurement space, and, obtain the ideal position in robot with five degrees of freedom 2 each joints by finding the solution the robot with five degrees of freedom inverse kinematics to robot with five degrees of freedom pedestal coordinate space; Then, by the ideal position of computer 1 to robot with five degrees of freedom 2 each joints of input, 2 motions of control robot with five degrees of freedom; Simultaneously, cooperate by optictracking device 3 and the passive marker 5 that is installed in the robot with five degrees of freedom end, the terminal pose of The real time measure robot with five degrees of freedom is also sent into computer 1, control method with the present invention's design is controlled in real time to the terminal pose of robot with five degrees of freedom, realizes accurate track following and location; At last, robot with five degrees of freedom 2 locking, doctor's operation that under robot with five degrees of freedom 2 auxiliary, undergos surgery.
Five degree-of-freedom manipulator is made up of arm and wrist, has five joints, takes the PPRRR configuration of two Gliding joints and three cradle heads; Arm segment has three joints, is respectively I, II and III joint, and I joint 11-1 adopts Gliding joint, and the direction of motion is perpendicular to horizontal plane; II joint 11-2 adopts Gliding joint, and is "T"-shaped vertical with the I joint; III joint 11-3 adopts cradle head, and axis is parallel with I joint motions direction; Wrist partly has two joints, is respectively IV and V joint, and IV joint 11-4 adopts cradle head, and axis is parallel with the III joints axes; V joint 11-5 adopts cradle head, and axis is vertical with the IV joints axes; Terminal apparatus is installed in V joint at mechanical arm, and it is parallel with the IV joints axes that axis is installed.All there is independent driving mechanisms in each joint, is made up of motor, decelerator, and arthrodial driving mechanism also includes ball screw.
Mechanical arm controller 10 comprises programmed logic controller PLC and stepper motor driver, adopt three PLC and five stepper motor driver combinations, five joints of corresponding control five degree-of-freedom manipulator 11, mechanical arm controller 10 communicates by RS232C Serial Communication Component and computer 1, accept the joint position order of computer 1 input, 11 motions of control five degree-of-freedom manipulator arrive assigned address.
The robot pose measurement of closed loop with the real-time method of feedback control system is: the passive marker 5 that is installed in the robot with five degrees of freedom end by tracking by optictracking device 3 is measured the terminal pose of robot with five degrees of freedom 2, calculate controlled quentity controlled variable by computer 1 according to pose that measures and the pose that expectation robot with five degrees of freedom 2 reaches, control robot with five degrees of freedom 2 further moves to revise deviation; Communicate by the RS232/RS422 serial ports between optictracking device 3 and the computer 1; Communicate by the RS232C serial ports between robot with five degrees of freedom 2 and the computer 1.
Described celebral operating robot system, passive marker 5 satisfies at the robot with five degrees of freedom 2 terminal geometry sites of installing: the distance between any two passive markers 5 can not be less than 50mm, minimal spatial separation between any two line segments that are made of passive marker 5 lines can not be less than 5mm, any two line segments are not parallel, and the angle between any two line segments can not be lower than 0.5 degree.
Implementation method based on the celebral operating robot system of optical tracking closed loop control, carry out pose measurement by means of the optical tracking technology, comprising: the closed loop control method of the automatic assignment test method of the easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings, robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings, the terminal pose of robot with five degrees of freedom.
The easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings is:
A. paste four medical science markers 6 by the doctor in patient head, require these four markers not on same plane, and any three markers are not on same straight line;
B. patient head is carried out CT or the scanning of MRI medical image, gained scan-image input computer 1;
That c. determines four medical science markers is scanned into the coordinate of picture point in the 3 D medical model space;
D. allow patient lie on the operation table 8, head uses moulding pillow 7 and operation table 8 relative fixed, and four medical science markers 6 using optics registration tool 4 mensuration patient head are at the spatial coordinate of optical measurement, and measured value is by optictracking device 3 input computers 1;
E. calculate the 3 D medical model space and the spatial mapping transformation of optical measurement by computer 1.
The automatic assignment test method of robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings is:
A. four points in chosen in advance robot with five degrees of freedom 2 work spaces require these four points should be in the measuring range of optictracking device 3, and four somes coplane not, and any three points are conllinear not;
B. send the joint position order that arrives above-mentioned four points by computer 1 to robot with five degrees of freedom 2, control robot with five degrees of freedom 2 moves to above-mentioned four points successively;
C. note the coordinate of above-mentioned four points in robot with five degrees of freedom pedestal coordinate space by computer 1; Simultaneously, cooperate by optictracking device 3 and the passive marker 5 that is installed in the robot with five degrees of freedom end, the coordinate of its end in optical measurement space when sequentially determining robot with five degrees of freedom 2 moves to above-mentioned four points, and send into computer 1;
D. by computer 1 according to above-mentioned four points at robot with five degrees of freedom pedestal coordinate space and the spatial coordinate figure of optical measurement, calculate robot with five degrees of freedom 2 pedestal coordinate spaces and the spatial mapping transformation of optical measurement, thereby finish automatic mapping.
The closed loop control method of the terminal pose of robot with five degrees of freedom is:
A. will transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance by computer 1, and calculate 1 and go out the desirable pose that need reach at current path point robot with five degrees of freedom end;
B. by the desirable pose of computer 1 according to current path point, carry out by of the coordinate transform of optical measurement space to robot with five degrees of freedom pedestal coordinate space, and, obtain the ideal position in each joint of robot with five degrees of freedom by finding the solution the robot with five degrees of freedom inverse kinematics;
C. by the ideal position of computer 1 to robot with five degrees of freedom 2 each joints of input, control robot with five degrees of freedom 2 moves;
D. cooperate by optictracking device 3 and the passive marker 5 that is installed in the robot with five degrees of freedom end, the pose of The real time measure robot with five degrees of freedom end, and send into computer 1;
E. by computer 1 the desirable pose of the posture information of The real time measure with the current path point of planning in advance compared, obtain both pose deviations;
F. by computer 1 according to the pose deviation according to predetermined control law, Position Control adopts proportional-integral-differential to regulate; Attitude control employing ratio is regulated, and the calculation correction controlled quentity controlled variable is revised the desirable pose of current path point, and begins the motor control of a new round, thereby realizes accurate track following and location.

Claims (10)

1, a kind of based on optical tracking closed loop control celebral operating robot system, comprise the hardware components that computer (1), robot with five degrees of freedom (2), optictracking device (3), optics registration tool (4), passive marker (5), medical science marker (6), moulding pillow (7) and operation table equipment such as (8) are formed, with assisted surgery planning and guiding software (9) part, be the position that receives medical image information, mensuration and definite focus, the auxiliary celebral operating robot that undergos surgery and plan and carry out surgical guide; It is characterized in that, constitute the robot pose measurement and real-time feedback control system of a closed loop by said computer (1), robot with five degrees of freedom (2), optictracking device (3) and passive marker (5), passive marker (5) is installed in the end of robot with five degrees of freedom (2), is followed the tracks of by optictracking device (3) all the time; Said robot with five degrees of freedom (2) comprises five degree-of-freedom manipulator (11) and mechanical arm controller (10); Said assisted surgery planning comprises digitized video input and pretreatment module with guiding software (9), focus extracts and module is implemented in three-dimensionalreconstruction module, surgery planning module and operation;
Use the undergo surgery concrete operations step in stage of the present invention to be:
One. when before undergoing surgery, preparing: on patient's head, paste four medical science markers (6), brain is carried out medical image scanning, and will scan gained medical image information input computer (1), the planning of utilization assisted surgery is determined focus with guiding software (9), the 3 D medical model of reconstruct focus and head, the sign of target spot and medical science sign image position and planning operation pathway undergo surgery;
Two. registering timing signal: allow patient lie on the operation table (8), head uses moulding pillow (7) and operation table (8) relative fixed, measure the coordinate of four the medical science markers (6) on patient's head on the one hand with optics registration tool (4), measured value is provided by optictracking device (3), and send into computer (1), by the mapping transformation between computer (1) the calculating 3 D medical model space and the optical measurement space; With the assignment test method of the present invention's design, cooperate on the other hand, measure the mapping transformation between robot with five degrees of freedom pedestal coordinate space and the optical measurement space automatically by robot with five degrees of freedom (2), optictracking device (3) and computer (1);
Three. when undergoing surgery enforcement: at first, will transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance, and calculate the desirable pose that need reach at current path point robot with five degrees of freedom end by computer (1); Then, by the desirable pose of computer (1) according to current path point, carry out by of the coordinate transform of optical measurement space, and, obtain the ideal position in each joint of robot with five degrees of freedom (2) by finding the solution the robot with five degrees of freedom inverse kinematics to robot with five degrees of freedom pedestal coordinate space; Then, import the ideal position in each joint to robot with five degrees of freedom (2), control robot with five degrees of freedom (2) motion by computer (1); Simultaneously, cooperate by optictracking device (3) and the passive marker (5) that is installed in the robot with five degrees of freedom end, the terminal pose of The real time measure robot with five degrees of freedom is also sent into computer (1), control method with the present invention's design is controlled in real time to the terminal pose of robot with five degrees of freedom, realizes accurate track following and location; At last, robot with five degrees of freedom (2) locking, doctor's operation that under robot with five degrees of freedom (2) auxiliary, undergos surgery.
2. celebral operating robot system according to claim 1 is characterized in that, five degree-of-freedom manipulator (11) is made up of arm and wrist, has five joints, takes the PPRRR configuration of two Gliding joints and three cradle heads; Arm segment has three joints, is respectively I, II and III joint, and I joint 11-1 adopts Gliding joint, and the direction of motion is perpendicular to horizontal plane; II joint 11-2 adopts Gliding joint, and is "T"-shaped vertical with the I joint; III joint 11-3 adopts cradle head, and axis is parallel with I joint motions direction; Wrist partly has two joints, is respectively IV and V joint, and IV joint 11-4 adopts cradle head, and axis is parallel with the III joints axes; V joint 11-5 adopts cradle head, and axis is vertical with the IV joints axes; Terminal apparatus is installed in V joint at mechanical arm, and it is parallel with the IV joints axes that axis is installed; All there is independent driving mechanisms in each joint, is made up of motor, decelerator, and arthrodial driving mechanism also includes ball screw.
3. celebral operating robot system according to claim 1, it is characterized in that, mechanical arm controller (10) comprises programmed logic controller (PLC) and stepper motor driver, adopt three PLC and five stepper motor driver combinations, five joints of corresponding control five degree-of-freedom manipulator (11), mechanical arm controller (10) communicates by RS232C Serial Communication Component and computer (1), accept the joint position order of computer (1) input, control five degree-of-freedom manipulator (11) motion arrives assigned address.
4. celebral operating robot system according to claim 1, it is characterized in that, the robot pose measurement of closed loop with the method for real-time feedback control system is: be installed in the terminal pose that the wooden passive marker of holding (5) of robot with five degrees of freedom is measured robot with five degrees of freedom (2) by optictracking device (3) by tracking, calculate controlled quentity controlled variable by computer (1) according to the pose that pose that measures and expectation robot with five degrees of freedom (2) reach, control robot with five degrees of freedom (2) further moves to revise deviation; Communicate by the RS232/RS422 serial ports between optictracking device (3) and the computer (1); Communicate by the RS232C serial ports between robot with five degrees of freedom (2) and the computer (1).
5. celebral operating robot system according to claim 1 is characterized in that, the quantity of passive marker (5) is at least 3.
6. celebral operating robot system according to claim 1, it is characterized in that, passive marker (5) satisfies at the terminal geometry site of installing of robot with five degrees of freedom (2): the distance between any two passive markers (5) can not be less than 50mm, minimal spatial separation between any two line segments that are made of passive marker (5) line can not be less than 5mm, any two line segments are not parallel, and the angle between any two line segments can not be lower than 0.5 degree.
7. implementation method based on the celebral operating robot system of optical tracking closed loop control, it is characterized in that, carry out pose measurement by means of the optical tracking technology, comprising: the closed loop control method of the automatic assignment test method of the easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings, robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings, the terminal pose of robot with five degrees of freedom.
8. according to the said implementation method of claim 7, it is characterized in that the easy assignment test method of the 3 D medical model space and the conversion of optical measurement spatial mappings is:
A. paste four medical science markers (6) by the doctor in patient head, require these four markers not on same plane, and any three markers are not on same straight line;
B. patient head is carried out CT or the scanning of MRI medical image, gained scan-image input computer (1);
That c. determines four medical science markers is scanned into the coordinate of picture point in the 3 D medical model space;
D. allow patient lie on the operation table (8), head uses moulding pillow (7) and operation table (8) relative fixed, four medical science markers (6) of measuring patient head with optics registration tool (4) are at the spatial coordinate of optical measurement, and measured value is by optictracking device (3) input computer (1);
E. calculate the 3 D medical model space and the spatial mapping transformation of optical measurement by computer (1).
9. according to the said method of claim 7, it is characterized in that the automatic assignment test method of robot with five degrees of freedom pedestal coordinate space and the conversion of optical measurement spatial mappings is:
A. four points in chosen in advance robot with five degrees of freedom (2) work space require these four points should be in the measuring range of optictracking device (3), and four somes coplane not, and any three points are conllinear not;
B. send the joint position order that arrives above-mentioned four points by computer (1) to robot with five degrees of freedom (2), control robot with five degrees of freedom (2) moves to above-mentioned four points successively;
C. note the coordinate of above-mentioned four points in robot with five degrees of freedom pedestal coordinate space by computer (1); Simultaneously, cooperate by optictracking device (3) and the passive marker (5) that is installed in the robot with five degrees of freedom end, the coordinate of its end in optical measurement space when sequentially determining robot with five degrees of freedom (2) moves to above-mentioned four points, and send into computer (1);
D. by computer (1) according to above-mentioned four points at robot with five degrees of freedom pedestal coordinate space and the spatial coordinate figure of optical measurement, calculate robot with five degrees of freedom (2) pedestal coordinate space and the spatial mapping transformation of optical measurement, thereby finish automatic mapping.
10. according to the said implementation method of claim 7, it is characterized in that the closed loop control method of robot with five degrees of freedom art end pose is:
A. will transform to the optical measurement space at the operation pathway that the 3 D medical model space is planned in advance by computer (1), and calculating (1) goes out the desirable pose that need reach at current path point robot with five degrees of freedom end;
B. by the desirable pose of computer (1) according to current path point, carry out by of the coordinate transform of optical measurement space to robot with five degrees of freedom pedestal coordinate space, and, obtain the ideal position in each joint of robot with five degrees of freedom by finding the solution the robot with five degrees of freedom inverse kinematics;
C. import the ideal position in each joint by computer (1) to robot with five degrees of freedom (2), control robot with five degrees of freedom (2) motion:
D. cooperate by optictracking device (3) and the passive marker (5) that is installed in the robot with five degrees of freedom end, the pose of The real time measure robot with five degrees of freedom end, and send into computer (1);
E. by computer (1) the desirable pose of the posture information of The real time measure with the current path point of planning in advance compared, obtain both pose deviations;
F. by computer (1) according to the pose deviation according to predetermined control law, Position Control adopts proportional-integral-differential to regulate; Attitude control employing ratio is regulated, and the calculation correction controlled quentity controlled variable is revised the desirable pose of current path point, and begins the motor control of a new round, thereby realizes accurate track following and location.
CNB200510122586XA 2005-12-22 2005-12-22 Celebral operating robot system based on optical tracking and closed-loop control and its realizing method Expired - Fee Related CN100464720C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200510122586XA CN100464720C (en) 2005-12-22 2005-12-22 Celebral operating robot system based on optical tracking and closed-loop control and its realizing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200510122586XA CN100464720C (en) 2005-12-22 2005-12-22 Celebral operating robot system based on optical tracking and closed-loop control and its realizing method

Publications (2)

Publication Number Publication Date
CN1985773A true CN1985773A (en) 2007-06-27
CN100464720C CN100464720C (en) 2009-03-04

Family

ID=38182807

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200510122586XA Expired - Fee Related CN100464720C (en) 2005-12-22 2005-12-22 Celebral operating robot system based on optical tracking and closed-loop control and its realizing method

Country Status (1)

Country Link
CN (1) CN100464720C (en)

Cited By (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101842755A (en) * 2007-08-28 2010-09-22 萨里大学 Inverse kinematics
CN101999938A (en) * 2010-10-15 2011-04-06 上海交通大学 Manual operating device for bone surgery assisted robot
WO2011069469A1 (en) * 2009-12-11 2011-06-16 Hospital Authority Stereoscopic visualization system for surgery
CN102198660A (en) * 2011-05-04 2011-09-28 上海海事大学 Robotic arm control system and action command control scheme based on brain-computer interface
CN102448398A (en) * 2009-06-01 2012-05-09 皇家飞利浦电子股份有限公司 Distance-based position tracking method and system
CN101870104B (en) * 2009-04-25 2012-09-19 鸿富锦精密工业(深圳)有限公司 Manipulator inverse moving method
CN102848389A (en) * 2012-08-22 2013-01-02 浙江大学 Realization method for mechanical arm calibrating and tracking system based on visual motion capture
CN102922521A (en) * 2012-08-07 2013-02-13 中国科学技术大学 Mechanical arm system based on stereo visual serving and real-time calibrating method thereof
CN103402453A (en) * 2011-03-03 2013-11-20 皇家飞利浦有限公司 System and method for automated initialization and registration of navigation system
CN101711127B (en) * 2007-04-19 2013-12-11 玛口外科股份有限公司 Implant planning using captured joint motion information
CN101612062B (en) * 2008-06-26 2014-03-26 北京石油化工学院 Sensing type six-joint mechanical arm for implementing method for navigating and positioning ultrasonic image
CN103759635A (en) * 2013-12-25 2014-04-30 合肥工业大学 Scanning measurement robot detection method allowing precision to be irrelevant to robot
CN103895023A (en) * 2014-04-04 2014-07-02 中国民航大学 Mechanical arm tail end tracking and measuring system and method based on coding azimuth device
CN104083217A (en) * 2014-07-03 2014-10-08 北京天智航医疗科技股份有限公司 Operation positioning device and method and robot operation system
CN104287833A (en) * 2010-10-21 2015-01-21 伊顿株式会社 Method and device for controlling/compensating movement of surgical robot
CN105025835A (en) * 2013-03-13 2015-11-04 史赛克公司 System for arranging objects in an operating room in preparation for a surgical procedure
CN105318838A (en) * 2015-11-26 2016-02-10 华科精准(北京)医疗科技有限公司 Method and system for single-plane calibration of relation between laser range finder and tail end of mechanical arm
CN105598970A (en) * 2016-02-01 2016-05-25 北京航空航天大学 Full-closed-loop robot system and control method thereof
CN105773629A (en) * 2015-01-14 2016-07-20 库卡罗伯特有限公司 Method for the alignment of a multiaxial manipulator with an input device
CN105852979A (en) * 2016-03-23 2016-08-17 北京柏惠维康科技有限公司 Medical image space localization device and method
CN105997260A (en) * 2016-07-28 2016-10-12 北京柏惠维康科技有限公司 Automatic lifting device, controlling method thereof and operation medical system
CN106462992A (en) * 2014-03-19 2017-02-22 株式会社乐博特思 Apparatus and method for editing and playing back robot motion, and computer-readable recording medium therefor
CN106420054A (en) * 2016-10-19 2017-02-22 胡磊 Anterior cruciate ligament stopping location and ligament tunnel location device combining with preoperative 3D planning information
CN106737683A (en) * 2017-01-11 2017-05-31 吉林省凯迪科技有限公司 The method of correction industrial robot off-line programing error in the field
CN106863332A (en) * 2017-04-27 2017-06-20 广东工业大学 A kind of robot visual orientation method and system
CN106914893A (en) * 2015-12-01 2017-07-04 西门子保健有限责任公司 Method and Medical robot equipment for running Medical robot equipment
CN106999250A (en) * 2014-12-17 2017-08-01 库卡罗伯特有限公司 System for the medical treatment of robot assisted
CN107174343A (en) * 2017-06-26 2017-09-19 南京普爱医疗设备股份有限公司 Three-dimensional orthopedic surgical navigational spatial registration method based on infrared optics alignment system
CN107405161A (en) * 2015-02-25 2017-11-28 奥林巴斯株式会社 Manipulator
CN107690317A (en) * 2015-06-04 2018-02-13 克瑞肖株式会社 The surgical robot system of operation pathway setting device including the operation pathway setting device and the operation pathway establishing method for operating robot
CN107898499A (en) * 2017-12-06 2018-04-13 上海波城医疗科技有限公司 Orthopaedics 3D region alignment system and method
CN107997822A (en) * 2017-12-06 2018-05-08 上海波城医疗科技有限公司 Minimally Invasive Surgery alignment system
CN108042172A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Bone surgery one channel alignment system
CN108042173A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Drilled tunnel alignment system
CN108042218A (en) * 2017-12-05 2018-05-18 北京军秀咨询有限公司 A kind of neurosurgery patient head mark automatic vision positioner and method
CN108042202A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Surgery mechanical arm alignment system
CN108201470A (en) * 2016-12-16 2018-06-26 上海铂联医疗科技有限公司 A kind of autonomous type tooth-implanting robot system and its device and method
CN108938090A (en) * 2018-08-20 2018-12-07 真健康(北京)医疗科技有限公司 A kind of the tumor operation robot system and its object localization method of precise positioning
CN109431604A (en) * 2018-10-09 2019-03-08 北京术锐技术有限公司 A kind of flexible operating robot and its application method for single hole invasive surgery formula
CN109431603A (en) * 2018-10-09 2019-03-08 北京术锐技术有限公司 Flexible robot's surgery systems and its application method for endoscope treatment art formula
CN109591007A (en) * 2017-09-30 2019-04-09 北京柏惠维康科技有限公司 A kind of robot space register method and device
CN109864806A (en) * 2018-12-19 2019-06-11 江苏集萃智能制造技术研究所有限公司 The Needle-driven Robot navigation system of dynamic compensation function based on binocular vision
CN109909646A (en) * 2019-02-02 2019-06-21 宁波吉利汽车研究开发有限公司 Modification method, device and the electronic equipment of robotic gun operating point
CN109938842A (en) * 2019-04-18 2019-06-28 王小丽 Facial surgical placement air navigation aid and device
CN110051436A (en) * 2018-01-18 2019-07-26 上海舍成医疗器械有限公司 Automation collaborative work package and its application in surgical instrument
CN110464465A (en) * 2019-08-14 2019-11-19 重庆文理学院 A kind of total knee replacement operation auxiliary robot
CN110584778A (en) * 2018-06-12 2019-12-20 上海舍成医疗器械有限公司 Method and device for adjusting object posture and application of device in automation equipment
CN111113430A (en) * 2019-12-31 2020-05-08 深圳市优必选科技股份有限公司 Robot and tail end control method and device thereof
CN111134702A (en) * 2019-12-17 2020-05-12 南京佗道医疗科技有限公司 Automatic control system and method for C-arm machine
CN111213098A (en) * 2017-08-17 2020-05-29 快砖知识产权私人有限公司 Communication system for interactive system
CN111249007A (en) * 2018-11-30 2020-06-09 微创(上海)医疗机器人有限公司 Medical robot and mechanical arm for clamping medical instrument
CN111728696A (en) * 2020-06-12 2020-10-02 天津理工大学 Craniotomy robot system and craniotomy method thereof
WO2020228087A1 (en) * 2019-05-15 2020-11-19 北京天智航医疗科技股份有限公司 Tracking device and localization system
CN113297997A (en) * 2021-05-31 2021-08-24 合肥工业大学 6-freedom face tracking method and device of non-contact physiological detection robot
CN113580141A (en) * 2021-08-18 2021-11-02 南京佗道医疗科技有限公司 Pose solving method for 6-axis mechanical arm
CN113693725A (en) * 2021-10-22 2021-11-26 杭州维纳安可医疗科技有限责任公司 Needle insertion path planning method, device, equipment and storage medium
CN113995512A (en) * 2021-02-07 2022-02-01 武汉光脉科技有限公司 Full-automatic tracking minimally invasive laser ablation surgical robot system and treatment method thereof
CN114010961A (en) * 2021-02-07 2022-02-08 武汉光脉科技有限公司 Full-automatic tracking laser pain treatment robot system and treatment method thereof
CN116077190A (en) * 2023-02-27 2023-05-09 上海华易镜升数字科技有限公司 System for improving accuracy of surgical robot system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167295A (en) * 1991-01-28 2000-12-26 Radionics, Inc. Optical and computer graphic stereotactic localizer
US6052611A (en) * 1997-11-28 2000-04-18 Picker International, Inc. Frameless stereotactic tomographic scanner for image guided interventional procedures
EP2070487B1 (en) * 2002-08-13 2014-03-05 NeuroArm Surgical, Ltd. Microsurgical robot system
CN1286438C (en) * 2003-12-26 2006-11-29 北京航空航天大学 Vision registering method for medical robot
CN100496429C (en) * 2005-02-06 2009-06-10 天津市华志计算机应用有限公司 Robot operation locating method of surgical operation navigation system based on optical positioning

Cited By (85)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101711127B (en) * 2007-04-19 2013-12-11 玛口外科股份有限公司 Implant planning using captured joint motion information
US8965582B2 (en) 2007-08-28 2015-02-24 The University Of Surrey Inverse kinematics
CN101842755B (en) * 2007-08-28 2014-12-10 萨里大学 Inverse kinematics
CN101842755A (en) * 2007-08-28 2010-09-22 萨里大学 Inverse kinematics
CN101612062B (en) * 2008-06-26 2014-03-26 北京石油化工学院 Sensing type six-joint mechanical arm for implementing method for navigating and positioning ultrasonic image
CN101870104B (en) * 2009-04-25 2012-09-19 鸿富锦精密工业(深圳)有限公司 Manipulator inverse moving method
CN102448398A (en) * 2009-06-01 2012-05-09 皇家飞利浦电子股份有限公司 Distance-based position tracking method and system
WO2011069469A1 (en) * 2009-12-11 2011-06-16 Hospital Authority Stereoscopic visualization system for surgery
CN101999938A (en) * 2010-10-15 2011-04-06 上海交通大学 Manual operating device for bone surgery assisted robot
CN104287833B (en) * 2010-10-21 2017-04-12 伊顿株式会社 Method and device for controlling/compensating movement of surgical robot
CN104287833A (en) * 2010-10-21 2015-01-21 伊顿株式会社 Method and device for controlling/compensating movement of surgical robot
CN103402453A (en) * 2011-03-03 2013-11-20 皇家飞利浦有限公司 System and method for automated initialization and registration of navigation system
US10290076B2 (en) 2011-03-03 2019-05-14 The United States Of America, As Represented By The Secretary, Department Of Health And Human Services System and method for automated initialization and registration of navigation system
CN102198660A (en) * 2011-05-04 2011-09-28 上海海事大学 Robotic arm control system and action command control scheme based on brain-computer interface
CN102922521A (en) * 2012-08-07 2013-02-13 中国科学技术大学 Mechanical arm system based on stereo visual serving and real-time calibrating method thereof
CN102922521B (en) * 2012-08-07 2015-09-09 中国科学技术大学 A kind of mechanical arm system based on stereoscopic vision servo and real-time calibration method thereof
CN102848389A (en) * 2012-08-22 2013-01-02 浙江大学 Realization method for mechanical arm calibrating and tracking system based on visual motion capture
CN102848389B (en) * 2012-08-22 2015-06-17 浙江大学 Realization method for mechanical arm calibrating and tracking system based on visual motion capture
US10410746B2 (en) 2013-03-13 2019-09-10 Stryker Corporation System and method for arranging objects in an operating room in preparation for surgical procedures
CN105025835A (en) * 2013-03-13 2015-11-04 史赛克公司 System for arranging objects in an operating room in preparation for a surgical procedure
US11183297B2 (en) 2013-03-13 2021-11-23 Stryker Corporation System and method for arranging objects in an operating room in preparation for surgical procedures
CN105025835B (en) * 2013-03-13 2018-03-02 史赛克公司 System for arranging objects in an operating room in preparation for a surgical procedure
CN103759635B (en) * 2013-12-25 2016-10-26 合肥工业大学 The scanning survey robot detection method that a kind of precision is unrelated with robot
CN103759635A (en) * 2013-12-25 2014-04-30 合肥工业大学 Scanning measurement robot detection method allowing precision to be irrelevant to robot
CN106462992B (en) * 2014-03-19 2019-12-10 株式会社乐博特思 Apparatus, method and recording medium for editing and reproducing robot motion
CN106462992A (en) * 2014-03-19 2017-02-22 株式会社乐博特思 Apparatus and method for editing and playing back robot motion, and computer-readable recording medium therefor
CN103895023A (en) * 2014-04-04 2014-07-02 中国民航大学 Mechanical arm tail end tracking and measuring system and method based on coding azimuth device
CN103895023B (en) * 2014-04-04 2015-08-19 中国民航大学 A kind of tracking measurement method of the mechanical arm tail end tracing measurement system based on coding azimuth device
CN104083217A (en) * 2014-07-03 2014-10-08 北京天智航医疗科技股份有限公司 Operation positioning device and method and robot operation system
CN106999250A (en) * 2014-12-17 2017-08-01 库卡罗伯特有限公司 System for the medical treatment of robot assisted
CN105773629A (en) * 2015-01-14 2016-07-20 库卡罗伯特有限公司 Method for the alignment of a multiaxial manipulator with an input device
CN107405161B (en) * 2015-02-25 2020-05-26 奥林巴斯株式会社 Mechanical arm
CN107405161A (en) * 2015-02-25 2017-11-28 奥林巴斯株式会社 Manipulator
CN107690317A (en) * 2015-06-04 2018-02-13 克瑞肖株式会社 The surgical robot system of operation pathway setting device including the operation pathway setting device and the operation pathway establishing method for operating robot
CN107690317B (en) * 2015-06-04 2020-07-10 克瑞肖株式会社 Surgical path setting device, surgical robot system including the same
CN105318838B (en) * 2015-11-26 2017-09-26 华科精准(北京)医疗科技有限公司 The monoplane scaling method of laser range finder and mechanical arm tail end relation
CN105318838A (en) * 2015-11-26 2016-02-10 华科精准(北京)医疗科技有限公司 Method and system for single-plane calibration of relation between laser range finder and tail end of mechanical arm
CN106914893A (en) * 2015-12-01 2017-07-04 西门子保健有限责任公司 Method and Medical robot equipment for running Medical robot equipment
CN106914893B (en) * 2015-12-01 2019-12-10 西门子保健有限责任公司 Method for operating a medical robot system and medical robot system
CN105598970A (en) * 2016-02-01 2016-05-25 北京航空航天大学 Full-closed-loop robot system and control method thereof
CN105852979A (en) * 2016-03-23 2016-08-17 北京柏惠维康科技有限公司 Medical image space localization device and method
CN105997260A (en) * 2016-07-28 2016-10-12 北京柏惠维康科技有限公司 Automatic lifting device, controlling method thereof and operation medical system
CN106420054A (en) * 2016-10-19 2017-02-22 胡磊 Anterior cruciate ligament stopping location and ligament tunnel location device combining with preoperative 3D planning information
CN108201470A (en) * 2016-12-16 2018-06-26 上海铂联医疗科技有限公司 A kind of autonomous type tooth-implanting robot system and its device and method
CN106737683A (en) * 2017-01-11 2017-05-31 吉林省凯迪科技有限公司 The method of correction industrial robot off-line programing error in the field
CN106863332A (en) * 2017-04-27 2017-06-20 广东工业大学 A kind of robot visual orientation method and system
CN107174343A (en) * 2017-06-26 2017-09-19 南京普爱医疗设备股份有限公司 Three-dimensional orthopedic surgical navigational spatial registration method based on infrared optics alignment system
CN111213098A (en) * 2017-08-17 2020-05-29 快砖知识产权私人有限公司 Communication system for interactive system
CN111213098B (en) * 2017-08-17 2024-03-15 快砖知识产权私人有限公司 Communication system for interactive system
CN109591007B (en) * 2017-09-30 2022-02-22 北京柏惠维康科技有限公司 Robot space registration method and device
CN109591007A (en) * 2017-09-30 2019-04-09 北京柏惠维康科技有限公司 A kind of robot space register method and device
CN108042218A (en) * 2017-12-05 2018-05-18 北京军秀咨询有限公司 A kind of neurosurgery patient head mark automatic vision positioner and method
CN108042172A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Bone surgery one channel alignment system
CN107997822A (en) * 2017-12-06 2018-05-08 上海波城医疗科技有限公司 Minimally Invasive Surgery alignment system
CN108042202A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Surgery mechanical arm alignment system
CN108042173A (en) * 2017-12-06 2018-05-18 上海波城医疗科技有限公司 Drilled tunnel alignment system
CN108042202B (en) * 2017-12-06 2019-10-25 上海波城医疗科技有限公司 Surgery mechanical arm positioning system
CN107898499A (en) * 2017-12-06 2018-04-13 上海波城医疗科技有限公司 Orthopaedics 3D region alignment system and method
CN110051436B (en) * 2018-01-18 2020-04-17 上海舍成医疗器械有限公司 Automated cooperative work assembly and application thereof in surgical instrument
CN110051436A (en) * 2018-01-18 2019-07-26 上海舍成医疗器械有限公司 Automation collaborative work package and its application in surgical instrument
CN110584778A (en) * 2018-06-12 2019-12-20 上海舍成医疗器械有限公司 Method and device for adjusting object posture and application of device in automation equipment
CN108938090A (en) * 2018-08-20 2018-12-07 真健康(北京)医疗科技有限公司 A kind of the tumor operation robot system and its object localization method of precise positioning
CN109431603A (en) * 2018-10-09 2019-03-08 北京术锐技术有限公司 Flexible robot's surgery systems and its application method for endoscope treatment art formula
CN109431604A (en) * 2018-10-09 2019-03-08 北京术锐技术有限公司 A kind of flexible operating robot and its application method for single hole invasive surgery formula
CN111249007A (en) * 2018-11-30 2020-06-09 微创(上海)医疗机器人有限公司 Medical robot and mechanical arm for clamping medical instrument
CN111249007B (en) * 2018-11-30 2021-09-07 上海微创医疗机器人(集团)股份有限公司 Medical robot and mechanical arm for clamping medical instrument
CN109864806A (en) * 2018-12-19 2019-06-11 江苏集萃智能制造技术研究所有限公司 The Needle-driven Robot navigation system of dynamic compensation function based on binocular vision
CN109909646A (en) * 2019-02-02 2019-06-21 宁波吉利汽车研究开发有限公司 Modification method, device and the electronic equipment of robotic gun operating point
CN109938842A (en) * 2019-04-18 2019-06-28 王小丽 Facial surgical placement air navigation aid and device
WO2020228087A1 (en) * 2019-05-15 2020-11-19 北京天智航医疗科技股份有限公司 Tracking device and localization system
CN110464465A (en) * 2019-08-14 2019-11-19 重庆文理学院 A kind of total knee replacement operation auxiliary robot
CN111134702A (en) * 2019-12-17 2020-05-12 南京佗道医疗科技有限公司 Automatic control system and method for C-arm machine
CN111134702B (en) * 2019-12-17 2024-04-12 佗道医疗科技有限公司 Automatic control system and method for C-arm machine
US11644840B2 (en) 2019-12-31 2023-05-09 Ubtech Robotics Corp Ltd Method for controlling a robot and its end-portions and device thereof
CN111113430B (en) * 2019-12-31 2021-07-27 深圳市优必选科技股份有限公司 Robot and tail end control method and device thereof
CN111113430A (en) * 2019-12-31 2020-05-08 深圳市优必选科技股份有限公司 Robot and tail end control method and device thereof
CN111728696A (en) * 2020-06-12 2020-10-02 天津理工大学 Craniotomy robot system and craniotomy method thereof
CN113995512A (en) * 2021-02-07 2022-02-01 武汉光脉科技有限公司 Full-automatic tracking minimally invasive laser ablation surgical robot system and treatment method thereof
CN114010961A (en) * 2021-02-07 2022-02-08 武汉光脉科技有限公司 Full-automatic tracking laser pain treatment robot system and treatment method thereof
CN113297997A (en) * 2021-05-31 2021-08-24 合肥工业大学 6-freedom face tracking method and device of non-contact physiological detection robot
CN113580141A (en) * 2021-08-18 2021-11-02 南京佗道医疗科技有限公司 Pose solving method for 6-axis mechanical arm
CN113693725B (en) * 2021-10-22 2022-02-22 杭州维纳安可医疗科技有限责任公司 Needle insertion path planning method, device, equipment and storage medium
CN113693725A (en) * 2021-10-22 2021-11-26 杭州维纳安可医疗科技有限责任公司 Needle insertion path planning method, device, equipment and storage medium
CN116077190A (en) * 2023-02-27 2023-05-09 上海华易镜升数字科技有限公司 System for improving accuracy of surgical robot system
CN116077190B (en) * 2023-02-27 2023-08-08 上海华易镜升数字科技有限公司 System for improving accuracy of surgical robot system

Also Published As

Publication number Publication date
CN100464720C (en) 2009-03-04

Similar Documents

Publication Publication Date Title
CN100464720C (en) Celebral operating robot system based on optical tracking and closed-loop control and its realizing method
CN101980838B (en) Surgery assistance system for guiding a surgical instrument
KR101200961B1 (en) Parallel kinematic machine, calibration method of parallel kinematic machine, and calibration program product
ES2457791T3 (en) Procedure to determine geometric errors in a machine tool or measuring machine
CN108789404A (en) A kind of serial manipulator kinematic calibration method of view-based access control model
CN103759635B (en) The scanning survey robot detection method that a kind of precision is unrelated with robot
US20050273202A1 (en) Method and device for improving the positioning accuracy of a manipulator
CN112833786B (en) Cabin attitude and pose measuring and aligning system, control method and application
CN102052895B (en) Method and device for detecting isocenter of radiological and radiotherapy equipment
CN104408299B (en) Robot location's error compensating method based on distance identification redundancy kinematics parameters
CN105444672B (en) The orthogonal plane scaling method and system of laser range finder and mechanical arm tail end relation
Chen et al. A ballbar test for measurement and identification the comprehensive error of tilt table
CN113331948B (en) Interventional operation robot system, calibration device and calibration method
Wang et al. Identification of strut and assembly errors of a 3-PRS serial–parallel machine tool
CN105318838B (en) The monoplane scaling method of laser range finder and mechanical arm tail end relation
CN110916799A (en) Puncture robot navigation system based on 5G network
CN201858965U (en) Detection device for isocenter points of radiation equipment and radiotherapy equipment
CN113974835B (en) Telecentric motionless point constraint-based surgical robot motion control method
CN105808882A (en) Calibration method and device for movement parameters of reptile-imitated four-foot walking robot
CN109176487A (en) A kind of cooperating joint section scaling method, system, equipment, storage medium
US20190111274A1 (en) Transcranial magnetic stimulation system and positioning assistance method and program
CN113888712A (en) Building structure construction deviation analysis method
CN110279470A (en) Dynamic regulation device, dynamic adjusting system and its application method
JPH11142139A (en) Measurement error correction method for automatic three-dimensional measurement device
CN102538675B (en) Optical detection system using additional light sources

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090304

Termination date: 20141222

EXPY Termination of patent right or utility model