WO2023029784A1 - Surgery robot tracking and moving method and system - Google Patents

Surgery robot tracking and moving method and system Download PDF

Info

Publication number
WO2023029784A1
WO2023029784A1 PCT/CN2022/106427 CN2022106427W WO2023029784A1 WO 2023029784 A1 WO2023029784 A1 WO 2023029784A1 CN 2022106427 W CN2022106427 W CN 2022106427W WO 2023029784 A1 WO2023029784 A1 WO 2023029784A1
Authority
WO
WIPO (PCT)
Prior art keywords
coordinate system
actuator
model
spatial position
target area
Prior art date
Application number
PCT/CN2022/106427
Other languages
French (fr)
Chinese (zh)
Inventor
张逸凌
刘星宇
Original Assignee
北京长木谷医疗科技有限公司
张逸凌
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 北京长木谷医疗科技有限公司, 张逸凌 filed Critical 北京长木谷医疗科技有限公司
Publication of WO2023029784A1 publication Critical patent/WO2023029784A1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/14Surgical saws ; Accessories therefor
    • A61B17/15Guides therefor
    • A61B17/154Guides therefor for preparing bone for knee prosthesis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/30Surgical robots
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/70Manipulators specially adapted for use in surgery
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/46Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor
    • A61F2/4603Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor for insertion or extraction of endoprosthetic joints or of accessories thereof
    • A61F2/461Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor for insertion or extraction of endoprosthetic joints or of accessories thereof of knees
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/246Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101Computer-aided simulation of surgical operations
    • A61B2034/102Modelling of surgical devices, implants or prosthesis
    • A61B2034/104Modelling the effect of the tool, e.g. the effect of an implanted prosthesis or for predicting the effect of ablation or burring
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101Computer-aided simulation of surgical operations
    • A61B2034/105Modelling of the patient, e.g. for ligaments or bones
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/107Visualisation of planned trajectories or target regions
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/108Computer aided selection or customisation of medical implants or cutting guides
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A61B2034/2046Tracking techniques
    • A61B2034/2055Optical tracking systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A61B2034/2046Tracking techniques
    • A61B2034/2065Tracking using image or pattern recognition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A61B2034/2068Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis using pointers, e.g. pointers having reference marks for determining coordinates of body points
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/46Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor
    • A61F2002/4632Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor using computer-controlled surgery, e.g. robotic surgery
    • A61F2002/4633Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor using computer-controlled surgery, e.g. robotic surgery for selection of endoprosthetic joints or for pre-operative planning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30008Bone

Definitions

  • the present application relates to the field of computer technology, in particular to a method and system for tracking and moving a surgical robot.
  • Traditional TKA surgery uses osteotomy plates for modular osteotomy, mainly referring to the patient’s preoperative radiographic X-ray films, measuring the bony landmarks during the operation, and manually placing the osteotomy plates for operation.
  • the traditional method mainly Relying on the skill and experience of the surgeon to judge the balance of osteotomy, prosthesis position and soft tissue, resulting in low surgical accuracy.
  • the main purpose of the present application is to provide a method and system for tracking and moving a surgical robot that can improve surgical accuracy.
  • a method for tracking and moving a surgical robot is provided.
  • the method for tracking and moving the surgical robot according to the present application includes:
  • the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time determine the real-time spatial position of the actuator at the end of the manipulator and the bone in the world coordinate system, and convert them according to the registration result Go to the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
  • the mechanical arm is controlled to limit the movement of the actuator within the current target area.
  • the method also includes:
  • the adjustment path is displayed in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
  • the steps of controlling the mechanical arm to limit the movement of the actuator within the current target area include:
  • the manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions
  • the direction in which the actuator cuts into the current target area is recorded as the depth direction
  • the direction within the current target area and perpendicular to the cutting direction is recorded as the horizontal direction
  • the direction perpendicular to the current target area is recorded as the vertical direction
  • the offset includes the deviation in the depth direction.
  • the preset stiffness value of the virtual spring in the depth direction and the preset stiffness value of the virtual spring in the lateral direction are all in the range of 0N/m to 500N/m;
  • the preset stiffness value of the virtual spring in the vertical direction ranges from 4000N/m to 5000N/m;
  • the preset stiffness value of the virtual spring in the direction of rotation with the vertical axis as the axis ranges from 0Nm/rad to 20Nm/rad;
  • the preset stiffness value of the virtual spring in the direction of rotation with the depth direction as the axis and the stiffness value of the virtual spring in the direction of rotation with the horizontal axis as the axis both range from 200Nm/rad to 300Nm/rad.
  • the step of registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system includes:
  • the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone;
  • the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system determine the spatial position of the marking point set in the world coordinate system when the surgical probe performs a marking operation on the solid bone;
  • the spatial position of the line point set in the world coordinate system is finely registered with the 3D model, and the registration result is obtained.
  • a system for tracking and moving a surgical robot is provided.
  • the system for tracking and moving the surgical robot according to the present application includes:
  • the registration module is configured to register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain registration result;
  • the tracking module is configured to determine the real-time spatial position of the actuator at the end of the mechanical arm and the bone in the world coordinate system according to the position of the tracker on the end of the mechanical arm and the tracker on the bone acquired by the tracking camera in real time, according to The registration result is transformed into the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
  • the motion control module is configured to control the mechanical arm according to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, so as to limit the movement of the actuator to the current target area.
  • the system also includes:
  • the plane alignment module is configured to determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model when the mechanical arm moves to the bone before the actuator runs;
  • the display module is configured to display the adjustment path in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is the same as the current target area. noodle.
  • the motion control module is further configured to:
  • the manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions
  • a computer device includes a memory and a processor, the memory stores a computer program that can run on the processor, and the processor implements the steps in each of the above method embodiments when executing the computer program.
  • a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the steps in the foregoing method embodiments are implemented.
  • the above-mentioned method, system, computer equipment and storage medium for tracking and moving the surgical robot perform registration on the three-dimensional model coordinate system in which the pre-acquired three-dimensional model of the skeleton is located and the world coordinate system in which the skeleton of the entity is located, so as to realize the registration of the world coordinate system.
  • the real-time spatial position of the actuator and bone in the world coordinate system can be tracked in real time, so that the spatial position of the actuator and bone in the 3D model coordinate system can be obtained according to the registration result.
  • determine the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm to control the mechanical arm and limit the movement of the actuator to the current target area. It can not only prevent damage to ligaments, blood vessels, nerves and other tissues, but also effectively avoid excessive osteotomy, improving the accuracy and safety of surgery.
  • Fig. 1 is an application environment diagram of a method for surgical robot tracking and movement in an embodiment
  • FIG. 2 is a schematic flow chart of a method for tracking and moving a surgical robot in an embodiment
  • Fig. 3 is a structural block diagram of a system for surgical robot tracking and movement in an embodiment
  • Figure 4 is an internal block diagram of a computer device in one embodiment.
  • the method for tracking and moving a surgical robot can be applied to a surgical robot system as shown in FIG. 1 , and the surgical robot system can be applied to joint surgery, for example, knee joint replacement surgery.
  • the surgical robot system includes a host computer main control system 102, a manipulator system 104 and an optical navigation and positioning system 106, wherein the host computer main control system 102 communicates with the manipulator system 104 and the optical navigation and positioning system 106 through the network, and the manipulator system 104 communicates with the optical navigation positioning system 106 through the network.
  • the host computer main control system 102 mainly includes a host computer and a display screen.
  • the upper computer is configured to perform various calculations and processing on images, and the upper computer also stores a library of prostheses, which is configured to plan and select models before surgery.
  • the manipulator system includes a manipulator control unit and a manipulator.
  • the manipulator control unit is configured to receive an osteotomy start signal sent by a host computer control system to control the movement of the manipulator.
  • the optical navigation and positioning system 106 is provided with a tracking camera, a tracer and a display screen; a tracer is installed on the end of the mechanical arm and bones (such as the femur and tibia of the knee joint), and a plurality of optical balls are arranged on the tracer , the upper computer determines the spatial position of the knee joint and the end of the mechanical arm in the world coordinate system by tracking the position of the optical ball through the tracking camera (binocular infrared camera).
  • the display screen is configured to display a three-dimensional model of the patient's skeleton, such as the knee joint.
  • the host computer can register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located with the world coordinate system where the skeleton of the entity is located through the optical navigation and positioning system 106, so as to register the world coordinate system to the 3D model Coordinate system to get the registration result.
  • the host computer also uses the optical navigation and positioning system 106 to determine the world coordinates of the actuators and bones at the end of the manipulator in the manipulator system 104 according to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time.
  • the real-time spatial position in the 3D model coordinate system it is transformed into the 3D model coordinate system according to the registration result, and the spatial position of the actuator and bone in the 3D model coordinate system is obtained; according to the spatial position of the current target area, the actuator at the end of the manipulator
  • the spatial position of the robot is controlled to limit the movement of the actuator to the current target area.
  • a method for tracking and moving a surgical robot including the following steps:
  • Step 202 registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain a registration result.
  • the registration process of registering the world coordinate system to the 3D model coordinate system can be divided into two stages: coarse registration stage and fine registration stage.
  • coarse registration stage the coarse registration matrix is obtained by roughly registering the spatial positions of the tracked intraoperative marker points in the world coordinate system and the acquired preoperative planning points in the 3D model coordinate system.
  • Fine registration performs fine registration with the 3D model on the spatial position of the tracked line point set in the world coordinate system according to the coarse registration matrix, and then obtains the registration result.
  • the three-dimensional skeleton model can be a skeleton model of a knee joint, a skeleton model of a hip joint, or other types of skeleton models.
  • the intraoperative marking points are multiple points marked on the bone by the doctor during the operation with a surgical probe.
  • the preoperative planning points are points planned in advance in the three-dimensional model for registration.
  • the line-marking point set is determined by the doctor using the surgical probe to carry out the line-line operation on the bone during the operation.
  • a plurality of optical balls are installed on the surgical probe, and the host computer determines the spatial positions of the intraoperative marking points and the marking point set in the world coordinate system according to the position of the optical balls tracked by the tracking camera.
  • step 202 registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system includes : According to the spatial position of the surgical probe in the world coordinate system acquired by the tracking camera, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone; obtain the 3D model of the bone The spatial position of the preoperative planning point on the bone in the 3D model coordinate system; the spatial position of the preoperative planning point in the 3D model coordinate system is roughly registered with the spatial position of the intraoperative marker point in the world coordinate system to obtain a rough Registration matrix; according to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the scribing point set in the world coordinate system when the surgical probe is scribing on the solid bone; The quasi-matrix precisely registers the spatial position
  • the host computer tracks the spatial position of the surgical probe in the world coordinate system through the tracking camera in the optical navigation positioning system, so that the surgical probe can be determined according to the spatial position of the surgical probe in the world coordinate system obtained by the tracking camera.
  • the spatial position of the intraoperative marker point in the world coordinate system when performing the point collection operation on the solid bone.
  • the spatial position of the preoperative planning point in the 3D model coordinate system is compared with the intraoperative marker point in the world coordinate system
  • Coarse registration is carried out at the spatial position of , and the coarse registration matrix is obtained.
  • the doctor’s intraoperative operation can be aligned with the 3D model, and the initial coordinate system between the world coordinate system and the 3D model coordinate system can be obtained. Transformation relationship, that is, the coarse registration matrix.
  • the second stage of fine registration can be performed. In the fine registration stage, no preoperative planning is required, and calibration equipment such as surgical probes are used to draw lines on the solid bone surface during the operation.
  • the surgical probe When the surgical probe is performing a marking operation on the solid bone, according to the spatial position of the tracer on the surgical probe in the world coordinate system during the marking operation process acquired by the tracking camera, determine the bone of the solid knee joint The spatial position of the dashed point set on the world coordinates.
  • roughly registering the spatial position of the preoperative planning point in the three-dimensional model coordinate system with the spatial position of the intraoperative marker point in the world coordinate system includes: by presetting the three-dimensional space
  • the point cloud search method triangulates the preoperative planning points according to the spatial position of the preoperative planning points in the 3D model coordinate system, and triangulates the intraoperative marking points according to the spatial position of the intraoperative marking points in the world coordinate system
  • the actual operation triangle sequence corresponding to the intraoperative marker points and the planning triangle sequence corresponding to the preoperative planning points are obtained;
  • the preoperative planning points are compared in the 3D model coordinate system according to the planning triangle sequence Correct the spatial position of the corrected preoperative planning point to obtain the corrected preoperative planning point; register the intraoperative marker point corresponding to the actual operation triangle sequence with the corrected preoperative planning point.
  • Triangulation processing refers to forming a triangle from every three points.
  • the principle of triangle composition is that the perimeter is the largest, and the points in the triangles can overlap, so as to obtain the practical triangle sequence and preoperative planning points corresponding to the marked points in the operation.
  • the corresponding planning triangle sequence is
  • the preoperative planning point is triangulated according to the spatial position of the preoperative planning point in the 3D model coordinate system, and the surgical operation is performed according to the spatial position of the intraoperative marker point in the world coordinate system.
  • Triangulation processing is performed on the marked points in the operation to obtain the practical triangle sequence corresponding to the intraoperative marked point and the planning triangle sequence corresponding to the preoperative planning point, including: according to the spatial position of the preoperative planning point in the three-dimensional model coordinate system, the preoperative planning point
  • the first three points of the intraoperative marker point form a triangle
  • the first three points of the intraoperative marker point form a triangle according to the spatial position of the intraoperative marker point in the world coordinate system
  • Starting from the fourth point select two points from the previous points
  • Points form a triangle with the current point to obtain the actual operation triangle sequence corresponding to the marked point in the operation and the planning triangle sequence corresponding to the preoperative planning point
  • the triangle formation sequence of the actual operation triangle sequence and the planning triangle sequence are the same.
  • the triangulation method of intraoperative marker points and preoperative planning points is the same.
  • the first three points automatically form a triangle
  • starting from the fourth point can be Select two points from the points to form a triangle with the current point.
  • the selection principle is that the perimeter of the triangle formed after selection is the largest. According to this principle, several triangle sequences are obtained.
  • the way of generating the triangular sequence of marked points during operation is the same as the way of planning points before operation.
  • correcting the spatial position of the preoperative planning point in the three-dimensional model coordinate system according to the planning triangle sequence includes: through the preset three-dimensional space point cloud search method , determine the second neighborhood space point set on the 3D model according to the spatial position of the preoperative planning point in the 3D model coordinate system; filter out the second target point set from the second neighborhood space point set; The spatial position of the previous planning point under the coordinates of the three-dimensional model is corrected to the position of the second target point set.
  • a second neighborhood space point set on the three-dimensional model of the preoperative planning point in the coordinate system of the three-dimensional model is determined through a preset three-dimensional space point cloud search method.
  • the second neighborhood space point set includes a large number of points.
  • the planning triangle sequence includes multiple triangles, and each triangle includes three triangle points.
  • the target point corresponding to each triangle point of the current triangle is screened in the second neighborhood space point set to obtain the first A set of target points.
  • the default screening strategy is that the triangle formed by the screened three target points is congruent with the triangle in the practical triangle sequence. Since the error of congruent triangles is extremely small, the spatial positions of the three triangular points of the current triangle under the coordinates of the 3D model can be corrected to the positions of the corresponding target points, and the correction process can be repeated to achieve continuous alignment of a large number of triangles in the planned triangle sequence. The spatial position of the pre-planning point under the coordinates of the three-dimensional model is corrected, and then the corrected pre-operative planning point closest to the intraoperative marker point is obtained.
  • the intraoperative marker points corresponding to the actual operation triangle sequence are registered with the corrected preoperative planning points through the registration algorithm to obtain the registration results.
  • the registration algorithm may be ICP (Iterative Closest Point, iterative closest point algorithm).
  • the preoperative planning points can become transparent.
  • the preoperatively planned three-dimensional model may include a three-dimensional femoral model and a three-dimensional tibial model, and the three-dimensional femoral model may be as shown in FIG. 3 , and the points in the figure are femoral marker points.
  • the three-dimensional tibial model can be shown in Figure 4, and the points in the figure are tibial marker points.
  • the registration points become transparent.
  • the tibial marker points in the intraoperative marker points are registered with the tibial planning points. After the registration is completed, the registration points become transparent.
  • the preoperative planning points are corrected according to the planning triangle sequence, and the corrected preoperative planning points are obtained. Since the triangle is unique and sufficient The stability of the registration is improved, and the preoperative planning points are corrected in advance, which effectively improves the accuracy of registration.
  • the second stage of fine registration can be performed.
  • the scribing area where the scribing operation can be performed is the key bone area on the bone surface, that is, the area containing key bone points.
  • the position of the tracer on the surgical probe is tracked by the tracking camera in the optical navigation positioning system, and the tracer on the surgical probe is in the world coordinate system during the marking process obtained according to the tracking camera Determine the spatial position of the line point set on the skeleton of the entity in world coordinates to obtain the line point set.
  • the surgical probe can be used to sample at the frequency S, and the point collection operation can be performed online, and the entire line segment can be subdivided into several point sets, so as to obtain the line point set.
  • the neighborhood space point set of the dashed point set on the 3D model can be determined first, and then the dashed point set can be aligned according to the neighborhood space point set and the spatial position of the dashed point set in the world coordinate system.
  • the spatial position in the three-dimensional model coordinate system is corrected, and then the corrected line point set is registered with the space position of the line point set in the world coordinate system.
  • finely registering the spatial position of the lined point set in the world coordinate system with the 3D model according to the coarse registration matrix includes: The spatial position in the world coordinate system is reflected back to the 3D model coordinate system, and the position of the dashed point set in the 3D model coordinate system is obtained; according to the position of the dashed point set in the 3D model coordinate system, the neighborhood space is performed on the 3D model Search to obtain the first neighborhood space point set; according to the space position of the first neighborhood space point set and the dashed point set in the world coordinate system, the spatial position of the dashed point set in the 3D model coordinate system is corrected to obtain The corrected set of dashed points; register the corrected set of dashed points with the spatial position of the set of dashed points in the world coordinate system.
  • the coarse registration matrix represents the conversion relationship between the world coordinate system and the 3D model coordinate system obtained by coarse registration.
  • the spatial position of the dashed point set in the world coordinate system can be reflected back to the 3D model coordinate system, so as to obtain the position of the dashed point set in the 3D model coordinate system.
  • the neighborhood space search can be performed on the 3D model according to the position of the dashed point set in the 3D model coordinate system to obtain the first neighborhood space point set.
  • the first neighborhood space point set is a neighborhood space point set corresponding to the dashed line point set in the three-dimensional model coordinate system.
  • correcting the spatial position of the dashed point set in the three-dimensional model coordinate system according to the spatial position of the first neighborhood spatial point set and the dashed point set in the world coordinate system includes: according to the dashed point set At the spatial position in the world coordinate system, triangular pairing is performed on the points in the dashed point set to obtain a paired triangle sequence; according to the first neighborhood space point set and the paired triangle sequence, the spatial position of the dashed point set in the 3D model coordinate system is obtained. Make corrections.
  • the dashed point set is composed of points on multiple line segments, for example, may include points in three line segments.
  • the points in the set of dashed points are paired in triangles, and a point is selected in each line segment, and every three points form a triangle.
  • the principle of composition is that the perimeter of the triangle is the largest.
  • a sequence of paired triangles is obtained.
  • the sequence of paired triangles includes a plurality of triangles.
  • the spatial position of the preoperative planning points under the coordinates of the 3D model is corrected through the second neighborhood space point set, and the line point set is adjusted in 3D according to the first neighborhood space point set and the paired triangle sequence.
  • the spatial position in the model coordinate system is corrected.
  • correcting the spatial position of the dashed point set in the three-dimensional model coordinate system according to the first neighborhood space point set and the paired triangle sequence includes: filtering out the first target point set from the first neighborhood space point set;
  • the paired triangle sequence corrects the spatial position of the dashed point set in the three-dimensional model coordinate system to the position of the first target point set.
  • the first neighborhood space point set includes a large number of points.
  • the paired triangle sequence includes multiple triangles, and each triangle includes three triangle points.
  • the target point corresponding to each triangle point of the current triangle can be screened in the second neighborhood space point set according to the paired triangle sequence to obtain the first A set of target points.
  • the default screening strategy is that the triangle formed by the screened three target points is congruent with the triangle in the paired triangle sequence.
  • the spatial positions of the three triangle points of the current triangle under the coordinates of the three-dimensional model can be respectively corrected to the positions corresponding to the target points in the first target point set, and the correction process can be repeated to achieve A large number of triangles continuously correct the spatial position of the dashed point set in the 3D model coordinates, making the spatial position of the dashed point set reflected in the 3D model coordinate system more accurate.
  • the corrected set of dashed points and the spatial position of the set of dashed points in the world coordinate system are registered through a registration algorithm to obtain a registration result.
  • the registration algorithm may be ICP (Iterative Closest Point, iterative closest point algorithm).
  • the registration result can be the transformation relationship between the final world coordinate system and the three-dimensional coordinates, and the accuracy of the intraoperative operation can be improved through the registration result.
  • the spatial position of the scribed point set on the skeleton of the entity in the world coordinate system is obtained through the scribe operation, so that the spatial position of the scribed point set in the world coordinate system is calculated according to the coarse registration matrix Performing fine registration with the 3D model, compared with the traditional point-taking registration algorithm, greatly improves the registration efficiency, and the registration accuracy is also greatly improved.
  • Step 204 according to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time, determine the real-time spatial position of the actuator at the end of the manipulator and the bone in the world coordinate system, and according to the registration result Transform it into the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system.
  • the optical navigation and positioning system tracks the tracker on the end of the manipulator and the tracker on the bone in real time through the tracking camera, and the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time , determine the real-time spatial position of the actuator and bones at the end of the manipulator in the world coordinate system.
  • the actuator at the end of a robotic arm could be a saw blade.
  • the real-time spatial position of the actuator at the end of the manipulator in the world coordinate system and the real-time spatial position of the bone in the world coordinate system are respectively transformed into the three-dimensional model coordinate system, and then the actuator at the end of the manipulator is obtained.
  • the spatial position in the 3D model coordinate system, and the spatial position of the bone in the 3D model coordinate system are respectively transformed into the three-dimensional model coordinate system, and then the actuator at the end of the manipulator is obtained.
  • Step 206 according to the spatial position of the current target area and the spatial position of the actuator at the end of the robotic arm, control the robotic arm to limit the movement of the actuator within the current target area.
  • the offset of the actuator relative to the current target area is determined according to the current spatial position of the actuator and the spatial position of the current target area of the bone; according to the offset The amount is used to control the robotic arm to limit the movement of the actuator to the current target area.
  • the actuator at the end of the manipulator can be determined
  • the current spatial position of the actuator and the current spatial position of each target area so that the spatial position of the actuator and the spatial position of the current target area can be determined in real time, and then based on the spatial position of the actuator, the spatial position of the current target area can be used to determine the relative position of the actuator.
  • the offset of the current target area can be used to determine the relative position of the actuator.
  • the target area may be a pre-planned osteotomy plane of the knee joint.
  • the osteotomy planes may include femoral osteotomy planes and tibial osteotomy planes, and the number of tibial osteotomy planes may be one plane area, namely the tibial osteotomy planes.
  • the number can include 5 plane areas, including the frontal femoral osteotomy plane, femoral anterior oblique osteotomy plane, femoral posterior condyle osteotomy plane, femoral posterior oblique osteotomy plane, and distal femoral osteotomy plane .
  • the pre-planned operation sequence of multiple target areas is displayed in the three-dimensional model, and the current target area is a target area selected from the multiple target areas in response to the operator.
  • the adjustment path of the actuator is determined according to the current spatial position of the actuator and the current target area in the three-dimensional solid model; it is displayed in the three-dimensional model The path is adjusted so that the operator operates the mechanical arm according to the adjusted path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
  • the displacement amount so as to determine the adjustment path of the actuator according to the displacement amount.
  • the adjustment path is displayed in the three-dimensional model to guide the doctor to support the robotic arm, so that the robotic arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
  • the step of controlling the mechanical arm to limit the movement of the actuator within the current target area includes:
  • the stiffness-damping model of the virtual spring is also called Cartesian Impedance Control Mode (CICM).
  • CICM Cartesian Impedance Control Mode
  • the behavior of the robot is compliance-sensitive and reacts to external influences such as obstacles or process forces. Applying an external force can cause the robot to deviate from the planned orbital path.
  • a relatively large stiffness value is set, and the stiffness value is greater than a predetermined threshold, so as to limit the movement of the actuator in a direction perpendicular to the current target area , so as to effectively prevent the actuator from deviating from the current target area.
  • the control robot After the actuator is aligned with the current target area, run the actuator. At this time, the control robot enters the state of the virtual spring damping model. In this state, the entire mechanical arm can be regarded as an approximate virtual spring, in any direction When a force is applied, the virtual spring obeys Hooke's law. In the direction perpendicular to the current target area, if the stiffness in this direction is large, the actuator will have a small deviation in this direction, so that the actuator can be stably limited to the current target area, and avoid The actuator moves in a direction perpendicular to the current target area, so as to minimize the actuator beyond the target area and reduce accidental injuries to patients.
  • the direction in which the actuator cuts into the current target area is recorded as the depth direction
  • the direction within the current target area and perpendicular to the cutting direction is recorded as the horizontal direction
  • the direction perpendicular to the current target area is recorded as the vertical direction
  • the displacement direction includes the offset value in the depth direction, the offset value in the horizontal direction, and the offset value in the vertical direction; the offset value for rotating around the depth direction, the offset value for rotating around the horizontal direction, and the offset value for rotating around the vertical direction.
  • the preset stiffness values of the virtual springs in the depth direction and the preset stiffness values of the virtual springs in the transverse direction both range from 0 N/m to 500 N/m. According to Hooke's law, when the force is constant, the smaller the stiffness, the larger the spring deformation. Therefore, setting the stiffness in the depth direction as small as possible can help the displacement of the actuator in this direction. In the transverse direction, the stiffness of the setting is also relatively small, which also helps the actuator to move in this direction for cutting.
  • the preset stiffness value of the virtual spring in the vertical direction ranges from 4000N/m to 5000N/m. According to Hooke's law, when the force is constant, the greater the stiffness, the smaller the spring deformation. Therefore, setting the stiffness in the Z direction as large as possible can help to avoid the displacement of the actuator in the Z direction, because if the displacement occurs in the Z direction, the saw blade will be directly separated from the current target area, which is likely to cause harm to the patient , which is not allowed.
  • the preset stiffness value of the virtual spring in the rotation direction with the vertical direction as the axis ranges from 0Nm/rad to 20Nm/rad, so that the saw blade can rotate in the current target area with the vertical direction Z as the axis.
  • the preset stiffness value of the virtual spring in the rotation direction with the depth direction as the axis and the stiffness value of the virtual spring in the rotation direction with the horizontal axis as the axis are both 200Nm/rad to 300Nm/rad, limiting the saw blade to the depth direction
  • the displacement of rotating on the axis and rotating on the horizontal axis further prevents the saw blade from breaking away from the current target area and ensures the safety of osteotomy.
  • the value range may also be other range values.
  • the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located is registered with the world coordinate system where the skeleton of the entity is located, so as to realize the registration of the world coordinate system to the 3D model coordinate system, which can be performed during osteotomy.
  • the real-time spatial positions of the actuators and bones in the world coordinate system are tracked in real time, so that the spatial positions of the actuators and bones in the 3D model coordinate system are obtained according to the registration results, and then the spatial position of the current target area, mechanical
  • the spatial position of the actuator at the end of the arm to control the robotic arm confines the motion of the actuator to the current target region. It can not only prevent damage to ligaments, blood vessels, nerves and other tissues, but also effectively avoid excessive osteotomy, improving the accuracy and safety of surgery.
  • the above method also includes the step of preoperative planning, which step includes:
  • the medical image of the bone is obtained, the medical image is segmented and three-dimensionally reconstructed to obtain the three-dimensional model of the bone; based on the three-dimensional model, the key parameters of the bone are determined; the type and model of the three-dimensional prosthesis model are determined based on the key parameters of the bone; the selected The three-dimensional prosthesis model is implanted into the three-dimensional model; the placement position and placement angle of the three-dimensional prosthesis model are adjusted based on the key bone parameters and the type and model of the three-dimensional prosthesis model.
  • the scanned image can be segmented through the neural network model, and can be segmented into regions of different granularities, such as the femoral region and the tibial region, or It can also be divided into femoral region, tibial region, fibula region and patella region as required; and then 3D reconstruction can be performed on the images of each segmented region to obtain a 3D model of each bone region.
  • the scanned image after acquiring the bone CT or nuclear magnetic image data of the target user, the scanned image can be segmented through the neural network model, and can be segmented into regions of different granularities as required, for example, in the acquired CT image of the knee joint Finally, it can be divided into femoral region and tibial region, or can be divided into femoral region, tibial region, fibula region and patella region according to needs; and then 3D reconstruction can be performed on the images of each segmented region to obtain a 3D model of each bone region.
  • Key bone parameters can include key anatomical points of bones, key axes of bones and bone size parameters.
  • Key anatomical points of bones can be identified based on deep learning algorithms, such as neural network models, and the identified key anatomical points of bones can be marked on the 3D model .
  • Bone size can include left and right femur diameter, femur anteroposterior diameter, tibial left and right diameter and tibial anteroposterior diameter.
  • the line connecting the medial and lateral borders of the tibia is determined, and the anteroposterior diameter of the tibia is determined according to the line connecting the anterior and posterior borders of the tibia.
  • the key axis of the bone is determined based on the key anatomical points of the bone, and the key angle of the bone is determined based on the key axis of the bone. However, based on the key axis of the bone and the key angle of the bone, it is helpful to determine the type and model of the three-dimensional prosthesis model.
  • the three-dimensional prosthesis model of the knee joint generally includes a three-dimensional femoral prosthesis model, a three-dimensional tibial prosthesis model, and a spacer model connecting the three-dimensional tibial prosthesis model and the three-dimensional femoral prosthesis model.
  • the three-dimensional prosthesis model can be a prosthesis model for total knee joint replacement currently available on the market.
  • the types of three-dimensional femoral prosthesis models include ATTUNE-PS, ATTUNE-CR, SIGMA-PS150, etc.
  • the models of ATTUNE-PS include 1, 2, 3, 3N, 4, 4N, 5, 5N, 6, 6N.
  • the selected three-dimensional skeletal prosthesis model is implanted into the three-dimensional bone model, and the placement position and placement angle of the three-dimensional prosthesis model are adjusted based on the key parameters of the skeleton and the type and model of the three-dimensional prosthesis model.
  • the three-dimensional visual display of the matching adjustment process and matching effect of the bone and the prosthesis is realized.
  • tibial varus angle femoral valgus angle
  • left and right tibial diameter tibial anteroposterior diameter
  • the three-dimensional model includes a three-dimensional femoral model
  • the three-dimensional prosthesis model includes a three-dimensional femoral prosthesis model
  • key bone parameters include femoral key parameters
  • femoral key parameters include femoral mechanical axis, femoral Condyle line, posterior condyle connection line, left and right femoral diameter and femoral anteroposterior diameter
  • the steps of adjusting the placement position and placement angle of the three-dimensional prosthesis model based on the key bone parameters and the type and model of the three-dimensional prosthesis model include: based on the left and right femoral diameter and the femoral Adjust the placement position of the 3D femoral prosthesis model; adjust the varus or valgus angle of the 3D femoral prosthesis model so that the cross section of the 3D femoral prosthesis model is perpendicular to the mechanical axis of the femur; adjust the varus angle of the 3D femoral
  • the installation position of the femoral prosthesis model satisfies that the femoral prosthesis model can cover the left and right sides of the femur and the front and back of the femur, the installation position is appropriate.
  • the femoral valgus angle and femoral varus angle can be determined according to the relative angle between the central axis of the femoral prosthesis model in the upper and lower direction of the coronal plane and the femoral force line, and according to the transverse axis of the femoral prosthesis model
  • the external rotation angle and internal rotation angle are determined by the relative angle to the condylar line;
  • the femoral flexion angle is determined by the angle between the femoral mechanical axis and the central axis of the femoral prosthesis model in the sagittal front-posterior direction.
  • the installation angle of the three-dimensional femoral prosthesis model is appropriate. For example, when the varus/valgus angle is adjusted to 0°, and the PCA is generally adjusted to 3°, then it is determined as the placement position of the femoral prosthesis model Adjust the installation angle to a suitable position.
  • the three-dimensional model also includes a three-dimensional tibial model
  • the three-dimensional femoral prosthesis model also includes a three-dimensional tibial prosthesis model
  • the bone key parameters also include tibial key parameters
  • the tibial key parameters include tibial mechanical axis, Left and right tibial diameter and tibial anteroposterior diameter
  • the steps of adjusting the placement position and placement angle of the 3D prosthesis model based on the key bone parameters and the type and model of the 3D prosthesis model include: adjusting the 3D tibial prosthesis based on the left and right tibial diameter and the tibial anteroposterior diameter Placement of the model; adjust the varus or valgus angle of the three-dimensional tibial prosthesis so that the mechanical axis of the tibia is perpendicular to the cross-section of the three-dimensional tibial prosthesis.
  • the posterior inclination angle of the tibial prosthesis can also be determined according to the design principles of the tibial prosthesis, and the adjustment of the flexion angle of the tibial prosthesis can be based on the patient's physiological The characteristics are determined, adjust to 0° or other, avoid notch (gap), Over.
  • the method further includes: performing a simulated osteotomy based on the matching relationship between the 3D prosthesis model and the 3D prosthesis model, Obtain a three-dimensional bone postoperative simulation model; perform motion simulation on the three-dimensional femoral postoperative simulation model including extension and flexion; determine the extension gap in the extension state, and determine the flexion gap in the flexion state; compare the extension gap and flexion Gap, to verify the matching of the three-dimensional prosthesis model.
  • the femoral osteotomy thickness is determined according to the design principles of the femoral prosthesis, and different femoral prostheses may correspond to different osteotomy thicknesses;
  • the osteotomy plane can be determined.
  • the simulated osteotomy is performed based on the matching relationship between the three-dimensional prosthesis model and the three-dimensional model, and a three-dimensional skeleton postoperative simulation model is obtained.
  • the motion simulation can be performed, and the extension gap and flexion gap can be determined through the extension position simulation map and the flexion position simulation map. Based on the extension gap and the flexion gap, it is determined whether the three-dimensional prosthesis model fits the osteotomized three-dimensional model. By simulating the installation effect of the prosthesis, it can be observed from different angles whether the size and position of the prosthesis are appropriate, whether there is collision or dislocation of the prosthesis, and then it is possible to accurately determine whether the prosthesis fits the bone. The user can determine whether the 3D prosthesis model can be adjusted through the final simulation image.
  • the preoperative planning method further includes: determining the three-dimensional coordinates of the center point of the femoral medullary cavity based on the three-dimensional femoral model; creating an intramedullary positioning analog rod by a circular fitting method; using the intramedullary positioning analog rod Determine the opening point of the femur.
  • the position of the needle entry point of the intramedullary positioning analog rod in the femur can also be determined during knee arthroplasty, wherein the apex of the intercondylar notch can be used as the needle entry point of the intramedullary positioning analog rod, The position of the point can be used as the opening point of the femur.
  • the intramedullary locating analog rod and the opening point of the femur are visualized on the three-dimensional bone model to guide the doctor to open the pulp.
  • a system for tracking and moving a surgical robot including: a registration module 302, a tracking module 304 and a motion control module 306; wherein:
  • the registration module 302 is configured to register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain the registration quasi-result.
  • the tracking module 304 is configured to determine the real-time spatial position of the actuator at the end of the mechanical arm and the bone in the world coordinate system according to the positions of the tracker on the end of the mechanical arm and the tracker on the bone acquired by the tracking camera in real time, Transform it into the 3D model coordinate system according to the registration result, and obtain the spatial positions of the actuator and bones in the 3D model coordinate system.
  • the motion control module 306 is configured to control the mechanical arm according to the spatial position of the current target area and the spatial position of the actuator at the end of the robotic arm, so as to limit the movement of the actuator within the current target area.
  • the above system also includes:
  • the plane alignment module is configured to determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model when the mechanical arm moves to the bone before the actuator runs;
  • the display module is configured to display the adjustment path in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is the same as the current target area. noodle.
  • the motion control module 306 is further configured to:
  • the manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions
  • the direction in which the actuator cuts into the current target area is marked as the depth direction
  • the direction within the current target area and perpendicular to the cutting direction is marked as the horizontal direction
  • the direction perpendicular to the current target area is marked as the vertical direction
  • the amount includes the offset value in the depth direction, the offset value in the horizontal direction, the offset value in the vertical direction, the offset value rotated around the depth direction, the offset value rotated around the horizontal direction, the offset value rotated around the vertical direction
  • the depth direction The value ranges of the preset stiffness values of the virtual springs above and the preset stiffness values of the virtual springs in the lateral direction are both 0N/m ⁇ 500N/m; the value ranges of the preset stiffness values of the virtual springs in the vertical direction are 4000N/m ⁇ 5000N/m; the value range of the preset stiffness value of the virtual spring in the rotation direction with the vertical axis as the axis is 0Nm/rad ⁇ 20Nm/rad; the preset stiffness value of the
  • registration module 304 is further configured to:
  • the spatial position of the surgical probe in the world coordinate system obtained by the tracking camera, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone; obtain the bone in the 3D model of the bone
  • the spatial position of the preoperative planning point on the 3D model coordinate system the spatial position of the preoperative planning point in the 3D model coordinate system is roughly registered with the spatial position of the intraoperative marker point in the world coordinate system, and the rough registration is obtained.
  • Quasi-matrix according to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the marking point set in the world coordinate system when the surgical probe is scribing on the solid bone; according to the coarse registration
  • the matrix precisely registers the spatial position of the line point set in the world coordinate system with the 3D model to obtain the registration result.
  • Each module in the above-mentioned system for tracking and moving the surgical robot can be fully or partially realized by software, hardware and a combination thereof.
  • the above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, and can also be stored in the memory of the computer device in the form of software, so that the processor can invoke and execute the corresponding operations of the above-mentioned modules.
  • a computer device is provided.
  • the computer device may be a server, and its internal structure may be as shown in FIG. 4 .
  • the computer device includes a processor, memory, network interface and database connected by a system bus.
  • the processor of the computer device is configured to provide computing and control capabilities.
  • the memory of the computer equipment includes a non-volatile storage medium and an internal memory.
  • the non-volatile storage medium stores an operating system, computer programs and databases.
  • the internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium.
  • the database of the computer device is configured to store data of a method of tracking and movement of a surgical robot.
  • the network interface of the computer device is configured to communicate with an external terminal via a network connection.
  • FIG. 4 is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation to the computer equipment on which the solution of the application is applied.
  • the specific computer equipment can be More or fewer components than shown in the figures may be included, or some components may be combined, or have a different arrangement of components.
  • a computer device including a memory and a processor, the memory stores a computer program, and the processor implements the steps in the foregoing embodiments when executing the computer program.
  • a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the steps in the foregoing embodiments are implemented.
  • Nonvolatile memory can include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory.
  • Volatile memory can include random access memory (RAM) or external cache memory.
  • RAM is available in many forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Chain Synchlink DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
  • SRAM Static RAM
  • DRAM Dynamic RAM
  • SDRAM Synchronous DRAM
  • DDRSDRAM Double Data Rate SDRAM
  • ESDRAM Enhanced SDRAM
  • SLDRAM Synchronous Chain Synchlink DRAM
  • Rambus direct RAM
  • DRAM direct memory bus dynamic RAM
  • RDRAM memory bus dynamic RAM

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Molecular Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Robotics (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Transplantation (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Quality & Reliability (AREA)
  • Cardiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Vascular Medicine (AREA)
  • Dentistry (AREA)
  • Multimedia (AREA)
  • Manipulator (AREA)

Abstract

A surgery robot tracking and moving method and system. The method comprises: registering a three-dimensional model coordinate system where a pre-acquired three-dimensional skeleton model is located with a world coordinate system where a solid skeleton is located, so as to register the world coordinate system to the three-dimensional model coordinate system, and thereby obtaining a registration result (S202); determining a real-time spatial location of an actuator at a tail end of a robotic arm and of the skeleton in the world coordinate system according to the locations of a tracer at the tail end of the robotic arm and a tracer on the skeleton that are acquired by a tracking camera in real time, and converting the real-time spatial locations in the three-dimensional model coordinate system according to the registration result, so as to obtain spatial locations of the actuator and the skeleton in the three-dimensional model coordinate system (S204); and controlling the robotic arm according to a spatial location of the current target area and the spatial location of the actuator at the tail end of the robotic arm, so as to limit the movement of the actuator to be within the current target area (S206). By means of the method, the accuracy of surgery can be improved.

Description

手术机器人追踪和移动的方法及系统Method and system for tracking and moving surgical robot
相关申请的交叉引用Cross References to Related Applications
本申请要求于2021年09月03日提交的申请号为202111035623.9,名称为“手术机器人追踪和移动的方法及系统”的中国专利申请的优先权,其通过引用方式全部并入本文。This application claims the priority of the Chinese patent application with application number 202111035623.9, entitled "Method and System for Tracking and Moving Surgical Robot", filed on September 03, 2021, which is incorporated herein by reference in its entirety.
技术领域technical field
本申请涉及计算机技术领域,尤其是涉及一种手术机器人追踪和移动方法及系统。The present application relates to the field of computer technology, in particular to a method and system for tracking and moving a surgical robot.
背景技术Background technique
传统TKA手术(全膝关节置换术)采用截骨板进行模块截骨,主要参考患者术前的影像学X线片,术中骨性标志进行测量后手动放置截骨板进行操作,传统方式主要依靠手术医生的技术与经验来判断截骨、假体位置和软组织的平衡情况,导致手术精准性较低。Traditional TKA surgery (total knee arthroplasty) uses osteotomy plates for modular osteotomy, mainly referring to the patient’s preoperative radiographic X-ray films, measuring the bony landmarks during the operation, and manually placing the osteotomy plates for operation. The traditional method mainly Relying on the skill and experience of the surgeon to judge the balance of osteotomy, prosthesis position and soft tissue, resulting in low surgical accuracy.
发明内容Contents of the invention
本申请的主要目的在于提供一种能够提高手术精准性的手术机器人追踪和移动的方法及系统。The main purpose of the present application is to provide a method and system for tracking and moving a surgical robot that can improve surgical accuracy.
为了实现上述目的,根据本申请的一个方面,提供了一种手术机器人追踪和移动的方法。In order to achieve the above purpose, according to one aspect of the present application, a method for tracking and moving a surgical robot is provided.
根据本申请的手术机器人追踪和移动的方法包括:The method for tracking and moving the surgical robot according to the present application includes:
对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果;Register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located with the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain a registration result;
根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置;According to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time, determine the real-time spatial position of the actuator at the end of the manipulator and the bone in the world coordinate system, and convert them according to the registration result Go to the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。According to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, the mechanical arm is controlled to limit the movement of the actuator within the current target area.
所述方法还包括:The method also includes:
在执行器运行前,当机械臂运动至骨骼处时,根据执行器、当前目标区域在三维实 体模型中的当前空间位置确定执行器的调节路径;Before the actuator runs, when the mechanical arm moves to the bone, determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model;
在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。The adjustment path is displayed in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内的步骤包括:According to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, the steps of controlling the mechanical arm to limit the movement of the actuator within the current target area include:
在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行执行器;After the plane of the actuator is coplanar with the current target area, when the robot arm is operated, run the actuator;
启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。Start the Cartesian damping control mode modeled on virtual springs and dampers. The manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions The offset Δx outputs the feedback force F opposite to the operated direction, F=Δx*C, so as to limit the movement of the actuator within the current target area.
执行器向当前目标区域切入的方向记为进深方向,在当前目标区域内并与切入方向垂直的方向记为横向,与当前目标区域垂直的方向记为垂直方向;偏移量包括进深方向的偏移值、横向的偏移值、垂直方向的偏移值、绕进深方向旋转的偏移值、绕横向旋转的偏移值、绕垂直方向旋转的偏移值;The direction in which the actuator cuts into the current target area is recorded as the depth direction, the direction within the current target area and perpendicular to the cutting direction is recorded as the horizontal direction, and the direction perpendicular to the current target area is recorded as the vertical direction; the offset includes the deviation in the depth direction. Offset value, horizontal offset value, vertical offset value, offset value rotated around the depth direction, offset value rotated around the horizontal direction, offset value rotated around the vertical direction;
进深方向上的虚拟弹簧的预设刚度值、横向上的虚拟弹簧的预设刚度值的取值范围均为0N/m~500N/m;The preset stiffness value of the virtual spring in the depth direction and the preset stiffness value of the virtual spring in the lateral direction are all in the range of 0N/m to 500N/m;
垂直方向上的虚拟弹簧的预设刚度值的取值范围为4000N/m~5000N/m;The preset stiffness value of the virtual spring in the vertical direction ranges from 4000N/m to 5000N/m;
以垂直方向为轴旋转方向上的虚拟弹簧的预设刚度值的取值范围为0Nm/rad~20Nm/rad;The preset stiffness value of the virtual spring in the direction of rotation with the vertical axis as the axis ranges from 0Nm/rad to 20Nm/rad;
以进深方向为轴旋转方向上的虚拟弹簧的预设刚度值、以横向为轴旋转方向上的虚拟弹簧的刚度值的取值范围均为200Nm/rad~300Nm/rad。The preset stiffness value of the virtual spring in the direction of rotation with the depth direction as the axis and the stiffness value of the virtual spring in the direction of rotation with the horizontal axis as the axis both range from 200Nm/rad to 300Nm/rad.
对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系的步骤包括:The step of registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system includes:
根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行采点操作时的术中标记点在世界坐标系的空间位置;According to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone;
获取骨骼的三维模型中骨骼上的术前规划点在三维模型坐标系的空间位置;Obtain the spatial position of the preoperative planning point on the bone in the three-dimensional model of the bone in the three-dimensional model coordinate system;
将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准,得到粗配准矩阵;Roughly register the spatial position of the preoperative planning point in the 3D model coordinate system with the spatial position of the intraoperative marker point in the world coordinate system to obtain a coarse registration matrix;
根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行划线操作时的划线点集在世界坐标系的空间位置;According to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the marking point set in the world coordinate system when the surgical probe performs a marking operation on the solid bone;
根据粗配准矩阵将划线点集在世界坐标系下的空间位置与三维模型进行精配准, 得到配准结果。According to the coarse registration matrix, the spatial position of the line point set in the world coordinate system is finely registered with the 3D model, and the registration result is obtained.
为了实现上述目的,根据本申请的另一方面,提供了一种手术机器人追踪和移动的系统。In order to achieve the above purpose, according to another aspect of the present application, a system for tracking and moving a surgical robot is provided.
根据本申请的手术机器人追踪和移动的系统包括:The system for tracking and moving the surgical robot according to the present application includes:
配准模块,被配置为对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果;The registration module is configured to register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain registration result;
追踪模块,被配置为根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置;The tracking module is configured to determine the real-time spatial position of the actuator at the end of the mechanical arm and the bone in the world coordinate system according to the position of the tracker on the end of the mechanical arm and the tracker on the bone acquired by the tracking camera in real time, according to The registration result is transformed into the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
运动控制模块,被配置为根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。The motion control module is configured to control the mechanical arm according to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, so as to limit the movement of the actuator to the current target area.
所述系统还包括:The system also includes:
平面对齐模块,被配置为在执行器运行前,当机械臂运动至骨骼处时,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定执行器的调节路径;The plane alignment module is configured to determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model when the mechanical arm moves to the bone before the actuator runs;
显示模块,被配置为在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。The display module is configured to display the adjustment path in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is the same as the current target area. noodle.
所述运动控制模块还被配置为:The motion control module is further configured to:
在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行所述执行器;After the plane of the actuator is coplanar with the current target area, when the mechanical arm is operated, the actuator is operated;
启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。Start the Cartesian damping control mode modeled on virtual springs and dampers. The manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions The offset Δx outputs the feedback force F opposite to the operated direction, F=Δx*C, so as to limit the movement of the actuator within the current target area.
一种计算机设备,包括存储器和处理器,所述存储器存储有可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述各个方法实施例中的步骤。A computer device includes a memory and a processor, the memory stores a computer program that can run on the processor, and the processor implements the steps in each of the above method embodiments when executing the computer program.
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述各个方法实施例中的步骤。A computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps in the foregoing method embodiments are implemented.
上述手术机器人追踪和移动的方法、系统、计算机设备和存储介质,对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,实现将世界坐标系配准到三维模型坐标系,可以在截骨过程中,实时追踪执行器、骨骼在世界坐 标系下的实时空间位置,从而根据配准结果获得执行器、骨骼在三维模型坐标系中的空间位置,进而确定当前目标区域的空间位置、机械臂末端的执行器的空间位置,以控制机械臂,将执行器的运动限定在当前目标区域内。既可以防止对韧带、血管、神经等组织的损伤,也可以有效避免截骨量过多,提高了手术的精准度和安全性。The above-mentioned method, system, computer equipment and storage medium for tracking and moving the surgical robot perform registration on the three-dimensional model coordinate system in which the pre-acquired three-dimensional model of the skeleton is located and the world coordinate system in which the skeleton of the entity is located, so as to realize the registration of the world coordinate system. Accurate to the 3D model coordinate system, during the osteotomy process, the real-time spatial position of the actuator and bone in the world coordinate system can be tracked in real time, so that the spatial position of the actuator and bone in the 3D model coordinate system can be obtained according to the registration result. Then determine the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm to control the mechanical arm and limit the movement of the actuator to the current target area. It can not only prevent damage to ligaments, blood vessels, nerves and other tissues, but also effectively avoid excessive osteotomy, improving the accuracy and safety of surgery.
附图说明Description of drawings
构成本申请的一部分的附图用来提供对本申请的进一步理解,使得本申请的其它特征、目的和优点变得更明显。本申请的示意性实施例附图及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:The accompanying drawings, which constitute a part of this application, are included to provide a further understanding of the application and make other features, objects and advantages of the application apparent. The drawings and descriptions of the schematic embodiments of the application are used to explain the application, and do not constitute an improper limitation to the application. In the attached picture:
图1是一个实施例中手术机器人追踪和移动的方法的应用环境图;Fig. 1 is an application environment diagram of a method for surgical robot tracking and movement in an embodiment;
图2为一个实施例中手术机器人追踪和移动的方法的流程示意图;FIG. 2 is a schematic flow chart of a method for tracking and moving a surgical robot in an embodiment;
图3为一个实施例中手术机器人追踪和移动的系统的结构框图;Fig. 3 is a structural block diagram of a system for surgical robot tracking and movement in an embodiment;
图4为一个实施例中计算机设备的内部结构图。Figure 4 is an internal block diagram of a computer device in one embodiment.
具体实施方式Detailed ways
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiment of the application will be clearly and completely described below in conjunction with the accompanying drawings in the embodiment of the application. Obviously, the described embodiment is only It is an embodiment of a part of the application, but not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by persons of ordinary skill in the art without creative efforts shall fall within the scope of protection of this application.
术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。The terms "comprising" and "having", as well as any variations thereof, are intended to cover a non-exclusive inclusion, for example, a process, method, system, product or device comprising a series of steps or elements need not be limited to those steps explicitly listed or units, but may include other steps or units not explicitly listed or inherent to the process, method, product or apparatus.
在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。In the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and embodiments.
本申请提供的手术机器人追踪和移动的方法,可以应用于如图1所示的手术机器人系统中,手术机器人系统可应用于关节手术中,例如,膝关节置换手术。手术机器人系统包括上位机主控系统102、机械臂系统104和光学导航定位系统106,其中,上位机主控系统102分别与机械臂系统104和光学导航定位系统106通过网络进行通信,机械臂系统104与光学导航定位系统106通过网络进行通信。The method for tracking and moving a surgical robot provided in this application can be applied to a surgical robot system as shown in FIG. 1 , and the surgical robot system can be applied to joint surgery, for example, knee joint replacement surgery. The surgical robot system includes a host computer main control system 102, a manipulator system 104 and an optical navigation and positioning system 106, wherein the host computer main control system 102 communicates with the manipulator system 104 and the optical navigation and positioning system 106 through the network, and the manipulator system 104 communicates with the optical navigation positioning system 106 through the network.
上位机主控系统102主要包括上位机和显示屏幕。其中,上位机被配置为对图像进行各种运算处理,上位机中还存储有假体库,被配置为手术前进行规划选型。机械臂系统包括机械臂控制单元和机械臂,机械臂控制单元被配置为接收上位机主控系统发送的截骨启动信号,控制机械臂运动。光学导航定位系统106设置有跟踪相机、示踪器和显示屏;机械臂末端、骨骼(如膝关节的股骨、胫骨)上均安装有示踪器,示踪器上设置有多个光学小球,上位机通过跟踪相机(双目红外相机)追踪光学小球的位置确定膝关节、机械臂末端在世界坐标系下的空间位置。显示屏被配置为显示患者骨骼(如膝关节)的三维模型。The host computer main control system 102 mainly includes a host computer and a display screen. Among them, the upper computer is configured to perform various calculations and processing on images, and the upper computer also stores a library of prostheses, which is configured to plan and select models before surgery. The manipulator system includes a manipulator control unit and a manipulator. The manipulator control unit is configured to receive an osteotomy start signal sent by a host computer control system to control the movement of the manipulator. The optical navigation and positioning system 106 is provided with a tracking camera, a tracer and a display screen; a tracer is installed on the end of the mechanical arm and bones (such as the femur and tibia of the knee joint), and a plurality of optical balls are arranged on the tracer , the upper computer determines the spatial position of the knee joint and the end of the mechanical arm in the world coordinate system by tracking the position of the optical ball through the tracking camera (binocular infrared camera). The display screen is configured to display a three-dimensional model of the patient's skeleton, such as the knee joint.
可选的,上位机通过光学导航定位系统106可以对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果。上位机还通过光学导航定位系统106根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂系统104中机械臂末端的执行器、骨骼在世界坐标系中的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置;根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。Optionally, the host computer can register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located with the world coordinate system where the skeleton of the entity is located through the optical navigation and positioning system 106, so as to register the world coordinate system to the 3D model Coordinate system to get the registration result. The host computer also uses the optical navigation and positioning system 106 to determine the world coordinates of the actuators and bones at the end of the manipulator in the manipulator system 104 according to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time. According to the real-time spatial position in the 3D model coordinate system, it is transformed into the 3D model coordinate system according to the registration result, and the spatial position of the actuator and bone in the 3D model coordinate system is obtained; according to the spatial position of the current target area, the actuator at the end of the manipulator The spatial position of the robot is controlled to limit the movement of the actuator to the current target area.
在一个实施例中,如图2所示,提供了一种手术机器人追踪和移动的方法,包括如下步骤:In one embodiment, as shown in Figure 2, a method for tracking and moving a surgical robot is provided, including the following steps:
步骤202,对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果。 Step 202, registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain a registration result.
将世界坐标系配准到三维模型坐标系的配准过程可以分为两个阶段:粗配准阶段和精配准阶段。粗配准阶段通过将追踪到的术中标记点在世界坐标系的空间位置与获取的术前规划点在三维模型坐标系的空间位置进行粗配准,得到粗配准矩阵。精配准根据粗配准矩阵将追踪到的划线点集在世界坐标系下的空间位置与三维模型进行精配准,进而得到配准结果。The registration process of registering the world coordinate system to the 3D model coordinate system can be divided into two stages: coarse registration stage and fine registration stage. In the coarse registration stage, the coarse registration matrix is obtained by roughly registering the spatial positions of the tracked intraoperative marker points in the world coordinate system and the acquired preoperative planning points in the 3D model coordinate system. Fine registration performs fine registration with the 3D model on the spatial position of the tracked line point set in the world coordinate system according to the coarse registration matrix, and then obtains the registration result.
骨骼三维模型可以膝关节的骨骼模型、髋关节的骨骼模型或其他类的骨骼模型。术中标记点为医生在术中利用手术探针在骨骼上标记的多个点。术前规划点为预先在三维模型中规划的用于配准的点。划线点集为医生在术中利用手术探针在骨骼上进行划线操作确定的。手术探针上安装有多个光学小球,上位机根据跟踪相机追踪到的光学小球的位置,确定术中标记点、划线点集分别在世界坐标系下的空间位置。The three-dimensional skeleton model can be a skeleton model of a knee joint, a skeleton model of a hip joint, or other types of skeleton models. The intraoperative marking points are multiple points marked on the bone by the doctor during the operation with a surgical probe. The preoperative planning points are points planned in advance in the three-dimensional model for registration. The line-marking point set is determined by the doctor using the surgical probe to carry out the line-line operation on the bone during the operation. A plurality of optical balls are installed on the surgical probe, and the host computer determines the spatial positions of the intraoperative marking points and the marking point set in the world coordinate system according to the position of the optical balls tracked by the tracking camera.
在其中一个实施例中,步骤202,对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系包括:根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的 骨骼上进行采点操作时的术中标记点在世界坐标系的空间位置;获取骨骼的三维模型中骨骼上的术前规划点在三维模型坐标系的空间位置;将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准,得到粗配准矩阵;根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行划线操作时的划线点集在世界坐标系的空间位置;根据粗配准矩阵将划线点集在世界坐标系下的空间位置与三维模型进行精配准,得到配准结果。In one of the embodiments, step 202, registering the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system includes : According to the spatial position of the surgical probe in the world coordinate system acquired by the tracking camera, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone; obtain the 3D model of the bone The spatial position of the preoperative planning point on the bone in the 3D model coordinate system; the spatial position of the preoperative planning point in the 3D model coordinate system is roughly registered with the spatial position of the intraoperative marker point in the world coordinate system to obtain a rough Registration matrix; according to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the scribing point set in the world coordinate system when the surgical probe is scribing on the solid bone; The quasi-matrix precisely registers the spatial position of the line point set in the world coordinate system with the 3D model to obtain the registration result.
在粗配准阶段,上位机通过光学导航定位系统中的追踪相机追踪手术探针在世界坐标系的空间位置,从而根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行采点操作时的术中标记点在世界坐标系的空间位置。在获取到骨骼的三维模型中骨骼上的术前规划点在三维模型坐标系的空间位置后,再将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准,得到粗配准矩阵。通过将三维模型坐标下的术中标记点与世界坐标系下的术前规划点进行粗配准,实现将医生的术中操作与三维模型进行对齐,得到世界坐标系与三维模型坐标系的初始转换关系,即粗配准矩阵。在粗配准的基础上,可以进行第二阶段的精配准。在精配准阶段,不需要进行术前规划,在术中利用手术探针等标定设备在实体的骨骼表面进行划线操作。当手术探针在实体的骨骼上进行划线操作时,根据跟踪相机获取的划线操作过程中,手术探针上的示踪器在世界坐标系下的空间位置,确定实体的膝关节的骨骼上的划线点集在世界坐标下的空间位置。In the rough registration stage, the host computer tracks the spatial position of the surgical probe in the world coordinate system through the tracking camera in the optical navigation positioning system, so that the surgical probe can be determined according to the spatial position of the surgical probe in the world coordinate system obtained by the tracking camera. The spatial position of the intraoperative marker point in the world coordinate system when performing the point collection operation on the solid bone. After obtaining the spatial position of the preoperative planning point on the bone in the 3D model coordinate system in the 3D model of the bone, the spatial position of the preoperative planning point in the 3D model coordinate system is compared with the intraoperative marker point in the world coordinate system Coarse registration is carried out at the spatial position of , and the coarse registration matrix is obtained. By roughly registering the intraoperative marker points under the 3D model coordinates with the preoperative planning points under the world coordinate system, the doctor’s intraoperative operation can be aligned with the 3D model, and the initial coordinate system between the world coordinate system and the 3D model coordinate system can be obtained. Transformation relationship, that is, the coarse registration matrix. On the basis of coarse registration, the second stage of fine registration can be performed. In the fine registration stage, no preoperative planning is required, and calibration equipment such as surgical probes are used to draw lines on the solid bone surface during the operation. When the surgical probe is performing a marking operation on the solid bone, according to the spatial position of the tracer on the surgical probe in the world coordinate system during the marking operation process acquired by the tracking camera, determine the bone of the solid knee joint The spatial position of the dashed point set on the world coordinates.
在本实施例的一种可选方式中,将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准包括:通过预设三维空间点云搜索方式,根据术前规划点在三维模型坐标系下的空间位置将术前规划点进行三角化处理,以及根据术中标记点在世界坐标系下的空间位置对术中标记点进行三角化处理,得到术中标记点对应的实操三角形序列和术前规划点对应的规划三角形序列;通过预设三维空间点云搜索方式,根据规划三角形序列对术前规划点在三维模型坐标系下的空间位置进行修正,得到修正后的术前规划点;将实操三角形序列对应的术中标记点与修正后的术前规划点进行配准。In an optional manner of this embodiment, roughly registering the spatial position of the preoperative planning point in the three-dimensional model coordinate system with the spatial position of the intraoperative marker point in the world coordinate system includes: by presetting the three-dimensional space The point cloud search method triangulates the preoperative planning points according to the spatial position of the preoperative planning points in the 3D model coordinate system, and triangulates the intraoperative marking points according to the spatial position of the intraoperative marking points in the world coordinate system Through the processing, the actual operation triangle sequence corresponding to the intraoperative marker points and the planning triangle sequence corresponding to the preoperative planning points are obtained; through the preset 3D space point cloud search method, the preoperative planning points are compared in the 3D model coordinate system according to the planning triangle sequence Correct the spatial position of the corrected preoperative planning point to obtain the corrected preoperative planning point; register the intraoperative marker point corresponding to the actual operation triangle sequence with the corrected preoperative planning point.
术中标记点和术前规划点均为点集。三角化处理是指将每三个点组成一个三角形,三角形的组成原则是周长最大,三角形之间中的点可以存在重合,从而得到术中标记点对应的实操三角形序列和术前规划点对应的规划三角形序列。Both intraoperative marker points and preoperative planning points are point sets. Triangulation processing refers to forming a triangle from every three points. The principle of triangle composition is that the perimeter is the largest, and the points in the triangles can overlap, so as to obtain the practical triangle sequence and preoperative planning points corresponding to the marked points in the operation. The corresponding planning triangle sequence.
通过预设三维空间点云搜索方式,根据术前规划点在三维模型坐标系下的空间位置将术前规划点进行三角化处理,以及根据术中标记点在世界坐标系下的空间位置对术中标 记点进行三角化处理,得到术中标记点对应的实操三角形序列和术前规划点对应的规划三角形序列包括:根据术前规划点在三维模型坐标系下的空间位置将术前规划点的前三个点组成三角形,以及根据术中标记点在世界坐标系下的空间位置将术中标记点的前三个点组成三角形;从第四个点开始,分别在之前的点中选取两个点,与当前点组成三角形,得到术中标记点对应的实操三角形序列和术前规划点对应的规划三角形序列;实操三角形序列与规划三角形序列的三角形组成顺序是相同的。术中标记点与术前规划点的三角化方式是相同的。Through the preset 3D space point cloud search method, the preoperative planning point is triangulated according to the spatial position of the preoperative planning point in the 3D model coordinate system, and the surgical operation is performed according to the spatial position of the intraoperative marker point in the world coordinate system. Triangulation processing is performed on the marked points in the operation to obtain the practical triangle sequence corresponding to the intraoperative marked point and the planning triangle sequence corresponding to the preoperative planning point, including: according to the spatial position of the preoperative planning point in the three-dimensional model coordinate system, the preoperative planning point The first three points of the intraoperative marker point form a triangle, and the first three points of the intraoperative marker point form a triangle according to the spatial position of the intraoperative marker point in the world coordinate system; starting from the fourth point, select two points from the previous points Points form a triangle with the current point to obtain the actual operation triangle sequence corresponding to the marked point in the operation and the planning triangle sequence corresponding to the preoperative planning point; the triangle formation sequence of the actual operation triangle sequence and the planning triangle sequence are the same. The triangulation method of intraoperative marker points and preoperative planning points is the same.
示例性的,对于术前规划点,假设术前规划点中的点云排列顺序为P1、P2、P3...Pn,则前三个点自动组成三角形,从第四个点开始可以从之前的点中挑选出两个点与当前点组成三角形,挑选原则为挑选过后所组成的三角形周长最大。根据这一原则得到若干三角形序列。术中标记点生成三角形序列的方式与术前规划点的方式是相同的。Exemplarily, for the preoperative planning points, assuming that the order of the point clouds in the preoperative planning points is P1, P2, P3...Pn, the first three points automatically form a triangle, and starting from the fourth point can be Select two points from the points to form a triangle with the current point. The selection principle is that the perimeter of the triangle formed after selection is the largest. According to this principle, several triangle sequences are obtained. The way of generating the triangular sequence of marked points during operation is the same as the way of planning points before operation.
在本可选的实现方式中,通过预设三维空间点云搜索方式,根据规划三角形序列对术前规划点在三维模型坐标系下的空间位置进行修正包括:通过预设三维空间点云搜索方式,根据术前规划点在三维模型坐标系下的空间位置在三维模型上确定第二邻域空间点集;在第二邻域空间点集中筛选出第二目标点集;根据规划三角形序列将术前规划点在三维模型坐标下的空间位置修正至第二目标点集的位置。In this optional implementation mode, by presetting the three-dimensional space point cloud search method, correcting the spatial position of the preoperative planning point in the three-dimensional model coordinate system according to the planning triangle sequence includes: through the preset three-dimensional space point cloud search method , determine the second neighborhood space point set on the 3D model according to the spatial position of the preoperative planning point in the 3D model coordinate system; filter out the second target point set from the second neighborhood space point set; The spatial position of the previous planning point under the coordinates of the three-dimensional model is corrected to the position of the second target point set.
可选的,通过预设三维空间点云搜索方式,确定三维模型坐标下系的术前规划点在三维模型上的第二邻域空间点集。第二邻域空间点集中包括大量的点。Optionally, a second neighborhood space point set on the three-dimensional model of the preoperative planning point in the coordinate system of the three-dimensional model is determined through a preset three-dimensional space point cloud search method. The second neighborhood space point set includes a large number of points.
规划三角形序列中包括多个三角形,每个三角形包括三个三角点,对于当前三角形,根据预设筛选策略在第二邻域空间点集中筛选当前三角形的每个三角点对应的目标点,得到第一目标点集。预设筛选策略为筛选出的三个目标点组成的三角形与实操三角形序列中的三角形为全等三角形。由于全等三角形误差极小,可以当前三角形的三个三角点在三维模型坐标下的空间位置分别修正至对应目标点的位置,重复该修正过程,实现通过规划三角形序列中的大量三角形不断对术前规划点在三维模型坐标下的空间位置进行修正,进而得到与术中标记点最相近的修正后的术前规划点。The planning triangle sequence includes multiple triangles, and each triangle includes three triangle points. For the current triangle, according to the preset screening strategy, the target point corresponding to each triangle point of the current triangle is screened in the second neighborhood space point set to obtain the first A set of target points. The default screening strategy is that the triangle formed by the screened three target points is congruent with the triangle in the practical triangle sequence. Since the error of congruent triangles is extremely small, the spatial positions of the three triangular points of the current triangle under the coordinates of the 3D model can be corrected to the positions of the corresponding target points, and the correction process can be repeated to achieve continuous alignment of a large number of triangles in the planned triangle sequence. The spatial position of the pre-planning point under the coordinates of the three-dimensional model is corrected, and then the corrected pre-operative planning point closest to the intraoperative marker point is obtained.
之后,通过配准算法将实操三角形序列对应的术中标记点与修正后的术前规划点进行配准,得到配准结果。例如,配准算法可以是ICP(Iterative Closest Point,迭代最近点算法)。当配准完成后,术前规划点可以变为透明。例如,术前规划的三维模型可以包括三维股骨模型和三维胫骨模型,三维股骨模型可以如图3所示,图中的点为股骨标记点。三维胫骨模型可以如图4所示,图中的点为胫骨标记点。将术中标记点中的股骨标记点与股骨规划点进行配准,配准完成后,该配准点变为透明。相应的,将术中标记点中的 胫骨标记点与胫骨规划点进行配准,配准完成后,该配准点变为透明。Afterwards, the intraoperative marker points corresponding to the actual operation triangle sequence are registered with the corrected preoperative planning points through the registration algorithm to obtain the registration results. For example, the registration algorithm may be ICP (Iterative Closest Point, iterative closest point algorithm). When the registration is completed, the preoperative planning points can become transparent. For example, the preoperatively planned three-dimensional model may include a three-dimensional femoral model and a three-dimensional tibial model, and the three-dimensional femoral model may be as shown in FIG. 3 , and the points in the figure are femoral marker points. The three-dimensional tibial model can be shown in Figure 4, and the points in the figure are tibial marker points. Register the femoral marker points in the intraoperative marker points with the femoral planning points. After the registration is completed, the registration points become transparent. Correspondingly, the tibial marker points in the intraoperative marker points are registered with the tibial planning points. After the registration is completed, the registration points become transparent.
在本实施例中,通过对术中标记点以及术前规划点进行三角化处理,根据规划三角形序列对术前规划点进行修正,得到修正后的术前规划点,由于三角形具有唯一性和足够的稳定性,且预先对术前规划点进行了修正,有效提高了配准的准确性。In this embodiment, by triangulating the intraoperative marker points and preoperative planning points, the preoperative planning points are corrected according to the planning triangle sequence, and the corrected preoperative planning points are obtained. Since the triangle is unique and sufficient The stability of the registration is improved, and the preoperative planning points are corrected in advance, which effectively improves the accuracy of registration.
在完成粗配准后,可以进行第二阶段的精配准。在精配准阶段,可以进行划线操作的划线区域是骨骼表面的关键骨骼区域,即包含关键骨骼点的区域。After the rough registration is completed, the second stage of fine registration can be performed. In the fine registration stage, the scribing area where the scribing operation can be performed is the key bone area on the bone surface, that is, the area containing key bone points.
示例性的,通过光学导航定位系统中的跟踪相机追踪手术探针上的示踪器的位置,根据跟踪相机获取的手术探针在划线过程中,其上的示踪器在世界坐标系下的空间位置,确定实体的骨骼上的划线点集在世界坐标下的空间位置,以得到划线点集。进一步的,在划线操作中,可以通过手术探针以频率S进行采样,在线上进行采点操作,将整条线段细分为若干点集,从而得到划线点集。Exemplarily, the position of the tracer on the surgical probe is tracked by the tracking camera in the optical navigation positioning system, and the tracer on the surgical probe is in the world coordinate system during the marking process obtained according to the tracking camera Determine the spatial position of the line point set on the skeleton of the entity in world coordinates to obtain the line point set. Further, in the line marking operation, the surgical probe can be used to sample at the frequency S, and the point collection operation can be performed online, and the entire line segment can be subdivided into several point sets, so as to obtain the line point set.
在精配准过程中,可以先确定划线点集在三维模型上的邻域空间点集,从而根据邻域空间点集中以及划线点集在世界坐标系下的空间位置对划线点集在三维模型坐标系下的空间位置进行修正,进而将修正后的划线点集与划线点集在世界坐标系下的空间位置进行配准。In the fine registration process, the neighborhood space point set of the dashed point set on the 3D model can be determined first, and then the dashed point set can be aligned according to the neighborhood space point set and the spatial position of the dashed point set in the world coordinate system. The spatial position in the three-dimensional model coordinate system is corrected, and then the corrected line point set is registered with the space position of the line point set in the world coordinate system.
在本实施例的一种可选方式中,根据粗配准矩阵将划线点集在世界坐标系下的空间位置与三维模型进行精配准包括:根据粗配准矩阵将划线点集在世界坐标系下的空间位置反射回三维模型坐标系中,得到划线点集在三维模型坐标系下的位置;根据划线点集在三维模型坐标系下的位置在三维模型上进行邻域空间搜索,得到第一邻域空间点集;根据第一邻域空间点集以及划线点集在世界坐标系下的空间位置对划线点集在三维模型坐标系下的空间位置进行修正,得到修正后的划线点集;将修正后的划线点集与划线点集在世界坐标系下的空间位置进行配准。In an optional manner of this embodiment, finely registering the spatial position of the lined point set in the world coordinate system with the 3D model according to the coarse registration matrix includes: The spatial position in the world coordinate system is reflected back to the 3D model coordinate system, and the position of the dashed point set in the 3D model coordinate system is obtained; according to the position of the dashed point set in the 3D model coordinate system, the neighborhood space is performed on the 3D model Search to obtain the first neighborhood space point set; according to the space position of the first neighborhood space point set and the dashed point set in the world coordinate system, the spatial position of the dashed point set in the 3D model coordinate system is corrected to obtain The corrected set of dashed points; register the corrected set of dashed points with the spatial position of the set of dashed points in the world coordinate system.
粗配准矩阵表示粗配准得到的世界坐标系与三维模型坐标系的转换关系。根据粗配准矩阵可以将划线点集在世界坐标系下的空间位置反射回三维模型坐标系中,从而得到划线点集在三维模型坐标系下的位置。由于三维模型对应三维模型坐标系,由此可根据划线点集在三维模型坐标系下的位置在三维模型上进行邻域空间搜索,得到第一邻域空间点集。第一邻域空间点集为三维模型坐标系下的划线点集对应的邻域空间点集。The coarse registration matrix represents the conversion relationship between the world coordinate system and the 3D model coordinate system obtained by coarse registration. According to the coarse registration matrix, the spatial position of the dashed point set in the world coordinate system can be reflected back to the 3D model coordinate system, so as to obtain the position of the dashed point set in the 3D model coordinate system. Since the 3D model corresponds to the 3D model coordinate system, the neighborhood space search can be performed on the 3D model according to the position of the dashed point set in the 3D model coordinate system to obtain the first neighborhood space point set. The first neighborhood space point set is a neighborhood space point set corresponding to the dashed line point set in the three-dimensional model coordinate system.
在可选方式中,根据第一邻域空间点集以及划线点集在世界坐标系下的空间位置对划线点集在三维模型坐标系下的空间位置进行修正包括:根据划线点集在世界坐标系下的空间位置将划线点集中的点进行三角形配对,得到配对三角形序列;根据第一 邻域空间点集以及配对三角形序列对划线点集在三维模型坐标系下的空间位置进行修正。In an optional manner, correcting the spatial position of the dashed point set in the three-dimensional model coordinate system according to the spatial position of the first neighborhood spatial point set and the dashed point set in the world coordinate system includes: according to the dashed point set At the spatial position in the world coordinate system, triangular pairing is performed on the points in the dashed point set to obtain a paired triangle sequence; according to the first neighborhood space point set and the paired triangle sequence, the spatial position of the dashed point set in the 3D model coordinate system is obtained. Make corrections.
划线点集是由多条线段上的点所组成的,例如,可以包括三条线段中的点。将划线点集中的点进行三角形配对,分别在每条线段中选取一个点,每三个点组成一个三角形,组成原则为三角形周长最大,按照该三角形配对方式,得到配对三角形序列。配对三角形序列包括多个三角形。The dashed point set is composed of points on multiple line segments, for example, may include points in three line segments. The points in the set of dashed points are paired in triangles, and a point is selected in each line segment, and every three points form a triangle. The principle of composition is that the perimeter of the triangle is the largest. According to the triangle pairing method, a sequence of paired triangles is obtained. The sequence of paired triangles includes a plurality of triangles.
采用粗配准中通过第二邻域空间点集对术前规划点在三维模型坐标下的空间位置进行修正的方式,根据第一邻域空间点集以及配对三角形序列对划线点集在三维模型坐标系下的空间位置进行修正。In the coarse registration, the spatial position of the preoperative planning points under the coordinates of the 3D model is corrected through the second neighborhood space point set, and the line point set is adjusted in 3D according to the first neighborhood space point set and the paired triangle sequence. The spatial position in the model coordinate system is corrected.
进一步的,根据第一邻域空间点集以及配对三角形序列对划线点集在三维模型坐标系下的空间位置进行修正包括:在第一邻域空间点集中筛选出第一目标点集;根据配对三角形序列将所述划线点集在三维模型坐标系下的空间位置修正至第一目标点集的位置。Further, correcting the spatial position of the dashed point set in the three-dimensional model coordinate system according to the first neighborhood space point set and the paired triangle sequence includes: filtering out the first target point set from the first neighborhood space point set; The paired triangle sequence corrects the spatial position of the dashed point set in the three-dimensional model coordinate system to the position of the first target point set.
第一邻域空间点集中包括大量的点。配对三角形序列中包括多个三角形,每个三角形包括三个三角点,对于当前三角形,可以根据配对三角形序列在第二邻域空间点集中筛选当前三角形的每个三角点对应的目标点,得到第一目标点集。预设筛选策略为筛选出的三个目标点组成的三角形与配对三角形序列中的三角形为全等三角形。由于全等三角形误差极小,可以当前三角形的三个三角点在三维模型坐标下的空间位置分别修正至第一目标点集中对应目标点的位置,重复该修正过程,实现通过配对三角形序列中的大量三角形不断对划线点集在三维模型坐标下的空间位置进行修正,使得划线点集反射到三维模型坐标系中的空间位置更为准确。The first neighborhood space point set includes a large number of points. The paired triangle sequence includes multiple triangles, and each triangle includes three triangle points. For the current triangle, the target point corresponding to each triangle point of the current triangle can be screened in the second neighborhood space point set according to the paired triangle sequence to obtain the first A set of target points. The default screening strategy is that the triangle formed by the screened three target points is congruent with the triangle in the paired triangle sequence. Since the error of congruent triangles is extremely small, the spatial positions of the three triangle points of the current triangle under the coordinates of the three-dimensional model can be respectively corrected to the positions corresponding to the target points in the first target point set, and the correction process can be repeated to achieve A large number of triangles continuously correct the spatial position of the dashed point set in the 3D model coordinates, making the spatial position of the dashed point set reflected in the 3D model coordinate system more accurate.
之后,通过配准算法将修正后的划线点集与划线点集在世界坐标系下的空间位置进行配准,得到配准结果。例如,配准算法可以是ICP(Iterative Closest Point,迭代最近点算法)。配准结果可以是最终得到的世界坐标系与三维坐标下的转换关系,通过配准结果可以提高术中操作的精准性。Afterwards, the corrected set of dashed points and the spatial position of the set of dashed points in the world coordinate system are registered through a registration algorithm to obtain a registration result. For example, the registration algorithm may be ICP (Iterative Closest Point, iterative closest point algorithm). The registration result can be the transformation relationship between the final world coordinate system and the three-dimensional coordinates, and the accuracy of the intraoperative operation can be improved through the registration result.
在本实施例中,通过划线操作获取实体的骨骼上的划线点集在世界坐标系下的空间位置,从而根据粗配准矩阵将所述划线点集在世界坐标系下的空间位置与所述三维模型进行精配准,与传统的取点配准算法相比,配准效率有利极大的提高,配准精准度也有较大提高。In this embodiment, the spatial position of the scribed point set on the skeleton of the entity in the world coordinate system is obtained through the scribe operation, so that the spatial position of the scribed point set in the world coordinate system is calculated according to the coarse registration matrix Performing fine registration with the 3D model, compared with the traditional point-taking registration algorithm, greatly improves the registration efficiency, and the registration accuracy is also greatly improved.
步骤204,根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结 果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置。 Step 204, according to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time, determine the real-time spatial position of the actuator at the end of the manipulator and the bone in the world coordinate system, and according to the registration result Transform it into the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system.
在术中,光学导航定位系统通过追踪相机实时追踪机械臂末端上的示踪器以及骨骼上的示踪器,根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系中的实时空间位置。例如,机械臂末端的执行器可以是锯片。During the operation, the optical navigation and positioning system tracks the tracker on the end of the manipulator and the tracker on the bone in real time through the tracking camera, and the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time , determine the real-time spatial position of the actuator and bones at the end of the manipulator in the world coordinate system. For example, the actuator at the end of a robotic arm could be a saw blade.
根据配准结果将机械臂末端的执行器在世界坐标系下的实时空间位置,以及骨骼在世界坐标系下的实时空间位置分别转换到三维模型坐标系中,进而得到机械臂末端的执行器在三维模型坐标系下的空间位置,以及骨骼在三维模型坐标系下的空间位置。According to the registration results, the real-time spatial position of the actuator at the end of the manipulator in the world coordinate system and the real-time spatial position of the bone in the world coordinate system are respectively transformed into the three-dimensional model coordinate system, and then the actuator at the end of the manipulator is obtained. The spatial position in the 3D model coordinate system, and the spatial position of the bone in the 3D model coordinate system.
步骤206,根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。 Step 206, according to the spatial position of the current target area and the spatial position of the actuator at the end of the robotic arm, control the robotic arm to limit the movement of the actuator within the current target area.
在一种实施方式中,在机械臂末端的执行器运行过程中,根据执行器的当前空间位置与骨骼的当前目标区域的空间位置确定执行器相对于当前目标区域的偏移量;根据偏移量,对机械臂进行控制,以将执行器的运动限定在当前目标区域内。In one embodiment, during the operation of the actuator at the end of the mechanical arm, the offset of the actuator relative to the current target area is determined according to the current spatial position of the actuator and the spatial position of the current target area of the bone; according to the offset The amount is used to control the robotic arm to limit the movement of the actuator to the current target area.
通过追踪相机(双目红外相机)实时追踪械臂末端设置的光学小球的位置、股骨区域上的光学小球的位置、胫骨区域上的光学小球的位置,可以确定机械臂末端的执行器的当前空间位置、各个目标区域的当前空间位置,从而可以实时确定执行器的空间位置、当前目标区域的空间位置,进而可基于执行器的空间位置,当前目标区域的空间位置确定执行器相对于当前目标区域的偏移量。By tracking the camera (binocular infrared camera) in real time to track the position of the optical ball at the end of the manipulator, the position of the optical ball on the femur area, and the position of the optical ball on the tibial area, the actuator at the end of the manipulator can be determined The current spatial position of the actuator and the current spatial position of each target area, so that the spatial position of the actuator and the spatial position of the current target area can be determined in real time, and then based on the spatial position of the actuator, the spatial position of the current target area can be used to determine the relative position of the actuator. The offset of the current target area.
在手术机器人应用于膝关节置换手术时,目标区域可以为预先规划的膝关节的截骨平面。截骨平面可以包括股骨截骨平面和胫骨截骨平面,对于胫骨截骨平面,其数量可以是1个平面区域,即胫骨截骨平面。对于股骨截骨平面,其数量可以包括5个平面区域,分别包括股骨前端截骨平面、股骨前斜截骨平面、股骨后髁截骨平面、股骨后斜截骨平面、股骨远端截骨平面。When the surgical robot is applied to knee joint replacement surgery, the target area may be a pre-planned osteotomy plane of the knee joint. The osteotomy planes may include femoral osteotomy planes and tibial osteotomy planes, and the number of tibial osteotomy planes may be one plane area, namely the tibial osteotomy planes. For the femoral osteotomy plane, the number can include 5 plane areas, including the frontal femoral osteotomy plane, femoral anterior oblique osteotomy plane, femoral posterior condyle osteotomy plane, femoral posterior oblique osteotomy plane, and distal femoral osteotomy plane .
三维模型中会显示多个目标区域被预先规划的操作顺序,当前目标区域则为响应于操作者从多个目标区域中选定的一个目标区域。The pre-planned operation sequence of multiple target areas is displayed in the three-dimensional model, and the current target area is a target area selected from the multiple target areas in response to the operator.
在本实施方式中,在执行器运行前,当机械臂被操作至骨骼处时,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定执行器的调节路径;在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。In this embodiment, before the actuator runs, when the mechanical arm is operated to the bone, the adjustment path of the actuator is determined according to the current spatial position of the actuator and the current target area in the three-dimensional solid model; it is displayed in the three-dimensional model The path is adjusted so that the operator operates the mechanical arm according to the adjusted path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
可选的,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定当前目标区域的空间位置与执行器的当前空间位置之间的位置差量,根据位置差量确定机 械臂被操作的位移量,从而根据位移量确定执行器的调节路径。在三维模型中显示调节路径,以引导医生托着机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。Optionally, determine the position difference between the spatial position of the current target area and the current spatial position of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model, and determine that the manipulator is operated according to the position difference The displacement amount, so as to determine the adjustment path of the actuator according to the displacement amount. The adjustment path is displayed in the three-dimensional model to guide the doctor to support the robotic arm, so that the robotic arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
在一种实施方式中,根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内的步骤包括:In one embodiment, according to the spatial position of the current target area and the spatial position of the actuator at the end of the robotic arm, the step of controlling the mechanical arm to limit the movement of the actuator within the current target area includes:
在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行执行器;启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。After the plane of the actuator is coplanar with the current target area, when the manipulator is operated, run the actuator; start the Cartesian damping control mode modeled on the virtual spring and damper, and the manipulator is based on the direction of multiple degrees of freedom The preset stiffness value C of each virtual spring and the offset Δx of the actuator relative to the current target area in the direction of multiple degrees of freedom will output the feedback force F opposite to the operated direction, F=Δx*C, so that the actuator The movement of is limited to the current target area.
在本实施方式中,虚拟弹簧的刚度-阻尼模型,也称为笛卡尔阻尼控制模式(Cartesian Impedance Control Mode,CICM)。在阻尼控制模式下,机器人的行为是顺从敏感的,并能对外部影响作出反应,外部影响比如,可以为障碍物或过程力。施加外力可使机器人离开计划的轨道路径。In this embodiment, the stiffness-damping model of the virtual spring is also called Cartesian Impedance Control Mode (CICM). In damped control mode, the behavior of the robot is compliance-sensitive and reacts to external influences such as obstacles or process forces. Applying an external force can cause the robot to deviate from the planned orbital path.
示例性地,在任意的一个目标区域中,在垂直于当前目标区域的方向上,设置比较大的刚度值,刚度值大于预定的阈值,以限制执行器在垂直于当前目标区域的方向上运动,从而有效避免执行器偏离当前目标区域。For example, in any target area, in a direction perpendicular to the current target area, a relatively large stiffness value is set, and the stiffness value is greater than a predetermined threshold, so as to limit the movement of the actuator in a direction perpendicular to the current target area , so as to effectively prevent the actuator from deviating from the current target area.
在执行器与当前目标区域对齐后,运行执行器,此时,控制机器人就进入到虚拟弹簧阻尼模型的状态,此状态下,整个机械臂可以看成一个近似的虚拟弹簧,在任何的方向上施加力后,虚拟弹簧都会遵从胡克定律。在垂直于当前目标区域的方向上,如果该方向的刚度很大,则执行器在该方向上的偏移就会很小,如此就可以将执行器稳定地限定在当前目标区域上,而避免执行器在垂直于当前目标区域的方向上运动,从而尽最大限度的降低执行器超出目标区域,减少给患者带来的误伤害。After the actuator is aligned with the current target area, run the actuator. At this time, the control robot enters the state of the virtual spring damping model. In this state, the entire mechanical arm can be regarded as an approximate virtual spring, in any direction When a force is applied, the virtual spring obeys Hooke's law. In the direction perpendicular to the current target area, if the stiffness in this direction is large, the actuator will have a small deviation in this direction, so that the actuator can be stably limited to the current target area, and avoid The actuator moves in a direction perpendicular to the current target area, so as to minimize the actuator beyond the target area and reduce accidental injuries to patients.
在一种实施方式中,执行器向当前目标区域切入的方向记为进深方向,在当前目标区域内并与切入方向垂直的方向记为横向,与当前目标区域垂直的方向记为垂直方向;偏移量方向包括进深方向的偏移值、横向的偏移值、垂直方向的偏移值;绕进深方向旋转的偏移值、绕横向旋转的偏移值、绕垂直方向旋转的偏移值。In one embodiment, the direction in which the actuator cuts into the current target area is recorded as the depth direction, the direction within the current target area and perpendicular to the cutting direction is recorded as the horizontal direction, and the direction perpendicular to the current target area is recorded as the vertical direction; The displacement direction includes the offset value in the depth direction, the offset value in the horizontal direction, and the offset value in the vertical direction; the offset value for rotating around the depth direction, the offset value for rotating around the horizontal direction, and the offset value for rotating around the vertical direction.
进深方向上的虚拟弹簧的预设刚度值、横向上的虚拟弹簧的预设刚度值的取值范围均为0N/m~500N/m。根据胡克定律,当力一定时,刚度越小,弹簧形变量越大。所以,在进深方向的刚度设置的尽可能小,可以有助于执行器在该方向发生位移。在横向上,设置的刚度也比较小,也有助于执行器在该方向移动,进行切割。The preset stiffness values of the virtual springs in the depth direction and the preset stiffness values of the virtual springs in the transverse direction both range from 0 N/m to 500 N/m. According to Hooke's law, when the force is constant, the smaller the stiffness, the larger the spring deformation. Therefore, setting the stiffness in the depth direction as small as possible can help the displacement of the actuator in this direction. In the transverse direction, the stiffness of the setting is also relatively small, which also helps the actuator to move in this direction for cutting.
垂直方向上的虚拟弹簧的预设刚度值的取值范围为4000N/m~5000N/m。根据胡克定律,当力一定时,刚度越大,弹簧形变量越小。所以,在Z方向的刚度设置的尽可能大,可以有助于避免执行器在Z方向发生位移,因为如果在Z方向发生位移之后,直接造成锯片脱离当前目标区域,容易给患者带来伤害,这是不允许的。The preset stiffness value of the virtual spring in the vertical direction ranges from 4000N/m to 5000N/m. According to Hooke's law, when the force is constant, the greater the stiffness, the smaller the spring deformation. Therefore, setting the stiffness in the Z direction as large as possible can help to avoid the displacement of the actuator in the Z direction, because if the displacement occurs in the Z direction, the saw blade will be directly separated from the current target area, which is likely to cause harm to the patient , which is not allowed.
以垂直方向为轴旋转方向上的虚拟弹簧的预设刚度值的取值范围为0Nm/rad~20Nm/rad,使得锯片可以垂直方向Z为轴在当前目标区域内旋转。The preset stiffness value of the virtual spring in the rotation direction with the vertical direction as the axis ranges from 0Nm/rad to 20Nm/rad, so that the saw blade can rotate in the current target area with the vertical direction Z as the axis.
以进深方向为轴旋转方向上的虚拟弹簧的预设刚度值、以横向为轴旋转方向上的虚拟弹簧的刚度值的取值范围均为200Nm/rad~300Nm/rad,限制锯片以进深方向为轴旋转、以横向为轴旋转的位移,进一步避免了锯片脱离当前目标区域,保证截骨的安全性。The preset stiffness value of the virtual spring in the rotation direction with the depth direction as the axis and the stiffness value of the virtual spring in the rotation direction with the horizontal axis as the axis are both 200Nm/rad to 300Nm/rad, limiting the saw blade to the depth direction The displacement of rotating on the axis and rotating on the horizontal axis further prevents the saw blade from breaking away from the current target area and ensures the safety of osteotomy.
当然,取值范围也可以是其他的范围值。可选的,设定不同自由度方向的弹簧刚度时,可以使用函数setStiffness(…)(type:double)进行设定。Certainly, the value range may also be other range values. Optionally, when setting the spring stiffness in different degrees of freedom directions, you can use the function setStiffness(…)(type:double) to set.
在本实施例中,对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,实现将世界坐标系配准到三维模型坐标系,可以在截骨过程中,实时追踪执行器、骨骼骨骼在世界坐标系下的实时空间位置,从而根据配准结果获得执行器、骨骼在三维模型坐标系中的空间位置,进而确定当前目标区域的空间位置、机械臂末端的执行器的空间位置,以控制机械臂,将执行器的运动限定在当前目标区域内。既可以防止对韧带、血管、神经等组织的损伤,也可以有效避免截骨量过多,提高了手术的精准度和安全性。In this embodiment, the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located is registered with the world coordinate system where the skeleton of the entity is located, so as to realize the registration of the world coordinate system to the 3D model coordinate system, which can be performed during osteotomy. During the process, the real-time spatial positions of the actuators and bones in the world coordinate system are tracked in real time, so that the spatial positions of the actuators and bones in the 3D model coordinate system are obtained according to the registration results, and then the spatial position of the current target area, mechanical The spatial position of the actuator at the end of the arm to control the robotic arm confines the motion of the actuator to the current target region. It can not only prevent damage to ligaments, blood vessels, nerves and other tissues, but also effectively avoid excessive osteotomy, improving the accuracy and safety of surgery.
在一个实施例中,上述方法还包括术前规划的步骤,该步骤包括:In one embodiment, the above method also includes the step of preoperative planning, which step includes:
在获取到骨骼的医学图像后,对医学图像进行分割和三维重建,得到骨骼的三维模型;基于三维模型,确定骨骼关键参数;基于骨骼关键参数确定三维假体模型的类型和型号;将选择的三维假体模型植入三维模型;基于骨骼关键参数和三维假体模型的类型和型号调整三维假体模型的安放位置和安放角度。After the medical image of the bone is obtained, the medical image is segmented and three-dimensionally reconstructed to obtain the three-dimensional model of the bone; based on the three-dimensional model, the key parameters of the bone are determined; the type and model of the three-dimensional prosthesis model are determined based on the key parameters of the bone; the selected The three-dimensional prosthesis model is implanted into the three-dimensional model; the placement position and placement angle of the three-dimensional prosthesis model are adjusted based on the key bone parameters and the type and model of the three-dimensional prosthesis model.
在本实施例中,在获取到目标用户的骨骼CT或者核磁图像数据后,可以通过神经网络模型对扫描图像进行图像分割,可以按需分割成不同粒度的区域,例如股骨区域和胫骨区域,或者还可以按需分割成股骨区域、胫骨区域、腓骨区域和髌骨区域;而后可以对分割后各个区域图像进行三维重建,得到各个骨骼区域的三维模型。In this embodiment, after obtaining the bone CT or MRI image data of the target user, the scanned image can be segmented through the neural network model, and can be segmented into regions of different granularities, such as the femoral region and the tibial region, or It can also be divided into femoral region, tibial region, fibula region and patella region as required; and then 3D reconstruction can be performed on the images of each segmented region to obtain a 3D model of each bone region.
在本实施例中,在获取到目标用户的骨骼CT或者核磁图像数据后,可以通过神经网络模型对扫描图像进行图像分割,可以按需分割成不同粒度的区域,例如在获取的膝关节CT图像后,分割成股骨区域和胫骨区域,或者还可以按需分割成股骨区域、胫 骨区域、腓骨区域和髌骨区域;而后可以对分割后各个区域图像进行三维重建,得到各个骨骼区域的三维模型。In this embodiment, after acquiring the bone CT or nuclear magnetic image data of the target user, the scanned image can be segmented through the neural network model, and can be segmented into regions of different granularities as required, for example, in the acquired CT image of the knee joint Finally, it can be divided into femoral region and tibial region, or can be divided into femoral region, tibial region, fibula region and patella region according to needs; and then 3D reconstruction can be performed on the images of each segmented region to obtain a 3D model of each bone region.
骨骼关键参数可包括骨骼关键解剖点、骨骼关键轴线和骨骼尺寸参数,骨骼关键解剖点可基于深度学习算法,例如神经网络模型,进行识别,并在三维模型上将识别的骨骼关键解剖点进行标记。Key bone parameters can include key anatomical points of bones, key axes of bones and bone size parameters. Key anatomical points of bones can be identified based on deep learning algorithms, such as neural network models, and the identified key anatomical points of bones can be marked on the 3D model .
骨骼尺寸可包括股骨左右径、股骨前后径、胫骨左右径和胫骨前后径,股骨左右径根据股骨内外侧缘连线,股骨前后径根据股骨前皮质切线和股骨后髁切线确定,胫骨左右径根据胫骨内外侧缘连线确定,胫骨前后径根据胫骨前后缘连线确定。Bone size can include left and right femur diameter, femur anteroposterior diameter, tibial left and right diameter and tibial anteroposterior diameter. The line connecting the medial and lateral borders of the tibia is determined, and the anteroposterior diameter of the tibia is determined according to the line connecting the anterior and posterior borders of the tibia.
骨骼关键轴线基于骨骼关键解剖点确定,基于骨骼关键轴线确定骨骼关键角度。而基于骨骼关键轴线、骨骼关键角度有助于确定三维假体模型的类型和型号。膝关节的三维假体模型一般性地包括三维股骨假体模型、三维胫骨假体和连接三维胫骨假体模型和三维股骨假体模型的垫片模型。The key axis of the bone is determined based on the key anatomical points of the bone, and the key angle of the bone is determined based on the key axis of the bone. However, based on the key axis of the bone and the key angle of the bone, it is helpful to determine the type and model of the three-dimensional prosthesis model. The three-dimensional prosthesis model of the knee joint generally includes a three-dimensional femoral prosthesis model, a three-dimensional tibial prosthesis model, and a spacer model connecting the three-dimensional tibial prosthesis model and the three-dimensional femoral prosthesis model.
三维假体模型可为目前市场上已有的全膝关节置换用的假体模型,三维假体模型有多种类型,每种类型的三维假体模型有多种型号。例如,三维股骨假体模型的类型有ATTUNE-PS、ATTUNE-CR、SIGMA-PS150等,ATTUNE-PS的型号有1、2、3、3N、4、4N、5、5N、6、6N。The three-dimensional prosthesis model can be a prosthesis model for total knee joint replacement currently available on the market. There are many types of three-dimensional prosthesis models, and each type of three-dimensional prosthesis model has multiple models. For example, the types of three-dimensional femoral prosthesis models include ATTUNE-PS, ATTUNE-CR, SIGMA-PS150, etc., and the models of ATTUNE-PS include 1, 2, 3, 3N, 4, 4N, 5, 5N, 6, 6N.
在本实施例中,将选择的三维骨骼假体模型植入所述三维骨骼模型,基于骨骼关键参数和三维假体模型的类型和型号调整三维假体模型的安放位置和安放角度。在本实施例中,实现了三维可视化显示骨骼与假体的匹配调节过程、匹配效果。在得到植入三维假体模型后的三维模型后,可以基于股骨外翻角、股骨内翻角、股骨外旋角、股骨内旋角、股骨左右径、股骨前后径确定股骨假体模型是否与三维股骨模型已安装适配。In this embodiment, the selected three-dimensional skeletal prosthesis model is implanted into the three-dimensional bone model, and the placement position and placement angle of the three-dimensional prosthesis model are adjusted based on the key parameters of the skeleton and the type and model of the three-dimensional prosthesis model. In this embodiment, the three-dimensional visual display of the matching adjustment process and matching effect of the bone and the prosthesis is realized. After obtaining the three-dimensional model after implanting the three-dimensional prosthesis model, it can be determined whether the femoral prosthesis model is consistent with The 3D femur model has been fitted.
可以基于胫骨内翻角、股骨外翻角、胫骨左右径、胫骨前后径确定胫骨假体模型是否与三维胫骨模型已安装适配。Based on the tibial varus angle, femoral valgus angle, left and right tibial diameter, and tibial anteroposterior diameter, it can be determined whether the tibial prosthetic model has been installed and adapted to the three-dimensional tibial model.
作为本实施例一种可选的实现方式,三维模型包括三维股骨模型,所述三维假体模型包括三维股骨假体模型,骨骼关键参数包括股骨关键参数,股骨关键参数包括股骨机械轴、股骨通髁线、后髁连线、股骨左右径和股骨前后径;基于骨骼关键参数和三维假体模型的类型和型号调整三维假体模型的安放位置和安放角度的步骤包括:基于股骨左右径和股骨前后径,调整三维股骨假体模型的放置位置;调整三维股骨假体模型的内翻角或外翻角,使三维股骨假体模型的横截面与股骨机械轴垂直;调整三维股骨假体的内旋角或外旋角,使股骨后髁角PCA(股骨通髁线与后髁连线在横断面的 投影线之间的夹角)在预设范围内。As an optional implementation of this embodiment, the three-dimensional model includes a three-dimensional femoral model, the three-dimensional prosthesis model includes a three-dimensional femoral prosthesis model, key bone parameters include femoral key parameters, femoral key parameters include femoral mechanical axis, femoral Condyle line, posterior condyle connection line, left and right femoral diameter and femoral anteroposterior diameter; the steps of adjusting the placement position and placement angle of the three-dimensional prosthesis model based on the key bone parameters and the type and model of the three-dimensional prosthesis model include: based on the left and right femoral diameter and the femoral Adjust the placement position of the 3D femoral prosthesis model; adjust the varus or valgus angle of the 3D femoral prosthesis model so that the cross section of the 3D femoral prosthesis model is perpendicular to the mechanical axis of the femur; adjust the varus angle of the 3D femoral prosthesis Rotation angle or external rotation angle, make the femoral posterior condyle angle PCA (the included angle between the femoral condyle line and the posterior condyle line on the cross-sectional projection line) within the preset range.
在本可选的实现方式中,当股骨假体模型的放置位置满足股骨假体模型能覆盖股骨左右、股骨前后,则安装位置合适。In this optional implementation manner, when the placement position of the femoral prosthesis model satisfies that the femoral prosthesis model can cover the left and right sides of the femur and the front and back of the femur, the installation position is appropriate.
可以基于股骨假体模型的当前位置,根据股骨假体模型在冠状面上下方向上的中轴线与股骨力线的相对角度确定股骨外翻角和股骨内翻角,根据股骨假体模型的横轴和通髁线的相对角度确定外旋角和内旋角;通过股骨机械轴和股骨假体模型在矢状面前后方向上的中轴线的角度确定股骨屈曲角。通过调整上述角度,可以确定三维股骨假体模型的安装角度是否合适,例如,当内/外翻角调整为0°时,PCA一般调整为3°时,则认定为股骨假体模型的安放位置和安放角度调整到合适的位置。Based on the current position of the femoral prosthesis model, the femoral valgus angle and femoral varus angle can be determined according to the relative angle between the central axis of the femoral prosthesis model in the upper and lower direction of the coronal plane and the femoral force line, and according to the transverse axis of the femoral prosthesis model The external rotation angle and internal rotation angle are determined by the relative angle to the condylar line; the femoral flexion angle is determined by the angle between the femoral mechanical axis and the central axis of the femoral prosthesis model in the sagittal front-posterior direction. By adjusting the above angles, it can be determined whether the installation angle of the three-dimensional femoral prosthesis model is appropriate. For example, when the varus/valgus angle is adjusted to 0°, and the PCA is generally adjusted to 3°, then it is determined as the placement position of the femoral prosthesis model Adjust the installation angle to a suitable position.
作为本实施例一种可选的实现方式,三维模型还包括三维胫骨模型,三维股骨假体模型还包括三维胫骨假体模型;骨骼关键参数还包括胫骨关键参数,胫骨关键参数包括胫骨机械轴、胫骨左右径和胫骨前后径;基于骨骼关键参数和三维假体模型的类型和型号调整三维假体模型的安放位置和安放角度的步骤包括:基于胫骨左右径和胫骨前后径,调整三维胫骨假体模型的安放位置;调整三维胫骨假体的内翻角或外翻角,使胫骨机械轴与三维胫骨假体的横截面垂直。As an optional implementation of this embodiment, the three-dimensional model also includes a three-dimensional tibial model, and the three-dimensional femoral prosthesis model also includes a three-dimensional tibial prosthesis model; the bone key parameters also include tibial key parameters, and the tibial key parameters include tibial mechanical axis, Left and right tibial diameter and tibial anteroposterior diameter; the steps of adjusting the placement position and placement angle of the 3D prosthesis model based on the key bone parameters and the type and model of the 3D prosthesis model include: adjusting the 3D tibial prosthesis based on the left and right tibial diameter and the tibial anteroposterior diameter Placement of the model; adjust the varus or valgus angle of the three-dimensional tibial prosthesis so that the mechanical axis of the tibia is perpendicular to the cross-section of the three-dimensional tibial prosthesis.
在本可选的实现方式中,除通过上述方式确定安装位置和角度外,还可以根据胫骨假体的设计原则确定胫骨假体的后倾角,胫骨假体的屈曲角的调整大小可以基于患者生理特性确定,调整为0°或其他,避免出现notch(缺口)、Over。In this optional implementation, in addition to determining the installation position and angle through the above methods, the posterior inclination angle of the tibial prosthesis can also be determined according to the design principles of the tibial prosthesis, and the adjustment of the flexion angle of the tibial prosthesis can be based on the patient's physiological The characteristics are determined, adjust to 0° or other, avoid notch (gap), Over.
作为本实施例一种可选的实现方式,在调整三维假体模型的安放位置和安放角度的步骤之后,方法还包括:基于三维假体模型与三维假体模型的匹配关系进行模拟截骨,得到三维骨骼术后模拟模型;对三维股骨术后模拟模型进行包括伸直位和屈曲位的运动模拟;在伸直位状态确定伸直间隙,在屈曲状态确定屈曲间隙;对比伸直间隙与屈曲间隙,对三维假体模型进行匹配性验证。As an optional implementation of this embodiment, after the step of adjusting the placement position and placement angle of the 3D prosthesis model, the method further includes: performing a simulated osteotomy based on the matching relationship between the 3D prosthesis model and the 3D prosthesis model, Obtain a three-dimensional bone postoperative simulation model; perform motion simulation on the three-dimensional femoral postoperative simulation model including extension and flexion; determine the extension gap in the extension state, and determine the flexion gap in the flexion state; compare the extension gap and flexion Gap, to verify the matching of the three-dimensional prosthesis model.
在本可选的实现方式中,根据股骨假体设计原则确定股骨截骨厚度,不同的股骨假体可能对应不同的截骨厚度;基于假体确定截骨厚度、假体与骨骼匹配后,便可确定截骨平面。In this optional implementation, the femoral osteotomy thickness is determined according to the design principles of the femoral prosthesis, and different femoral prostheses may correspond to different osteotomy thicknesses; The osteotomy plane can be determined.
在调整好三维假体模型的安放位置和安放角度之后,基于三维假体模型与三维模型的匹配关系进行模拟截骨,得到三维骨骼术后模拟模型。After adjusting the placement position and placement angle of the three-dimensional prosthesis model, the simulated osteotomy is performed based on the matching relationship between the three-dimensional prosthesis model and the three-dimensional model, and a three-dimensional skeleton postoperative simulation model is obtained.
在得到三维骨骼术后模拟模型后,进行运动模拟,可以通过伸直位模拟图、屈曲位模拟图,确定伸直间隙、屈曲间隙。基于伸直间隙和屈曲间隙,确定三维假体模型是否与截骨后的三维模型适配。通过对假体的安装效果进行模拟可从不同角度观察假 体大小、位置是否合适,是否出现假体碰撞、异位,进而能够精确地确定假体与骨骼是否适配。用户可通过该最终的模拟图像,确定是否可以对三维假体模型进行调整,如果更换骨骼假体的类型和型号,则可重新调用假体库,基于新的骨骼假体模型生成置换后的三维骨骼术后模拟模型。通过对术后的预期效果进行模拟,可以使最终得到的三维假体模型与患者的膝关节更加匹配。在一种实施方式中,术前规划方法还包括:基于所述三维股骨模型确定股骨髓腔中心点的三维坐标;通过圆形拟合法创建髓内定位模拟杆;由所述髓内定位模拟杆确定股骨开髓点。After obtaining the three-dimensional skeletal postoperative simulation model, the motion simulation can be performed, and the extension gap and flexion gap can be determined through the extension position simulation map and the flexion position simulation map. Based on the extension gap and the flexion gap, it is determined whether the three-dimensional prosthesis model fits the osteotomized three-dimensional model. By simulating the installation effect of the prosthesis, it can be observed from different angles whether the size and position of the prosthesis are appropriate, whether there is collision or dislocation of the prosthesis, and then it is possible to accurately determine whether the prosthesis fits the bone. The user can determine whether the 3D prosthesis model can be adjusted through the final simulation image. If the type and model of the bone prosthesis are changed, the prosthesis library can be called again to generate the replaced 3D prosthesis model based on the new bone prosthesis model. Skeletal postoperative simulation model. By simulating the expected postoperative effect, the final 3D prosthetic model can be more closely matched to the patient's knee joint. In one embodiment, the preoperative planning method further includes: determining the three-dimensional coordinates of the center point of the femoral medullary cavity based on the three-dimensional femoral model; creating an intramedullary positioning analog rod by a circular fitting method; using the intramedullary positioning analog rod Determine the opening point of the femur.
在可选的实现方式中,在膝关节置换术中还可以确定股骨髓内定位模拟杆入针点的位置,其中髁间窝的顶点可作为髓内定位模拟杆的入针点位置,入针点的位置即可作为股骨开髓点。在术中,三维骨骼模型上可视化显示髓内定位模拟杆和股骨开髓点,引导医生开髓。In an optional implementation, the position of the needle entry point of the intramedullary positioning analog rod in the femur can also be determined during knee arthroplasty, wherein the apex of the intercondylar notch can be used as the needle entry point of the intramedullary positioning analog rod, The position of the point can be used as the opening point of the femur. During the operation, the intramedullary locating analog rod and the opening point of the femur are visualized on the three-dimensional bone model to guide the doctor to open the pulp.
在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。The steps shown in the flow diagrams of the figures may be implemented in a computer system, such as a set of computer-executable instructions, and, although a logical order is shown in the flow diagrams, in some cases, may be executed differently from this The steps shown or described are performed in the order shown or described.
在一个实施例中,如图3所示,提供了一种手术机器人追踪和移动的系统,包括:配准模块302、追踪模块304和运动控制模块306;其中:In one embodiment, as shown in FIG. 3 , a system for tracking and moving a surgical robot is provided, including: a registration module 302, a tracking module 304 and a motion control module 306; wherein:
配准模块302,被配置为对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果。The registration module 302 is configured to register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain the registration quasi-result.
追踪模块304,被配置为根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置。The tracking module 304 is configured to determine the real-time spatial position of the actuator at the end of the mechanical arm and the bone in the world coordinate system according to the positions of the tracker on the end of the mechanical arm and the tracker on the bone acquired by the tracking camera in real time, Transform it into the 3D model coordinate system according to the registration result, and obtain the spatial positions of the actuator and bones in the 3D model coordinate system.
运动控制模块306,被配置为根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。The motion control module 306 is configured to control the mechanical arm according to the spatial position of the current target area and the spatial position of the actuator at the end of the robotic arm, so as to limit the movement of the actuator within the current target area.
在一个实施例中,上述系统还包括:In one embodiment, the above system also includes:
平面对齐模块,被配置为在执行器运行前,当机械臂运动至骨骼处时,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定执行器的调节路径;The plane alignment module is configured to determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model when the mechanical arm moves to the bone before the actuator runs;
显示模块,被配置为在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。The display module is configured to display the adjustment path in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is the same as the current target area. noodle.
在一个实施例中,运动控制模块306还被配置为:In one embodiment, the motion control module 306 is further configured to:
在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行所述执行器;After the plane of the actuator is coplanar with the current target area, when the mechanical arm is operated, the actuator is operated;
启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。Start the Cartesian damping control mode modeled on virtual springs and dampers. The manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions The offset Δx outputs the feedback force F opposite to the operated direction, F=Δx*C, so as to limit the movement of the actuator within the current target area.
在一个实施例中,执行器向当前目标区域切入的方向记为进深方向,在当前目标区域内并与切入方向垂直的方向记为横向,与当前目标区域垂直的方向记为垂直方向;偏移量包括进深方向的偏移值、横向的偏移值、垂直方向的偏移值、绕进深方向旋转的偏移值、绕横向旋转的偏移值、绕垂直方向旋转的偏移值;进深方向上的虚拟弹簧的预设刚度值、横向上的虚拟弹簧的预设刚度值的取值范围均为0N/m~500N/m;垂直方向上的虚拟弹簧的预设刚度值的取值范围为4000N/m~5000N/m;以垂直方向为轴旋转方向上的虚拟弹簧的预设刚度值的取值范围为0Nm/rad~20Nm/rad;以进深方向为轴旋转方向上的虚拟弹簧的预设刚度值、以横向为轴旋转方向上的虚拟弹簧的刚度值的取值范围均为200Nm/rad~300Nm/rad。In one embodiment, the direction in which the actuator cuts into the current target area is marked as the depth direction, the direction within the current target area and perpendicular to the cutting direction is marked as the horizontal direction, and the direction perpendicular to the current target area is marked as the vertical direction; The amount includes the offset value in the depth direction, the offset value in the horizontal direction, the offset value in the vertical direction, the offset value rotated around the depth direction, the offset value rotated around the horizontal direction, the offset value rotated around the vertical direction; the depth direction The value ranges of the preset stiffness values of the virtual springs above and the preset stiffness values of the virtual springs in the lateral direction are both 0N/m~500N/m; the value ranges of the preset stiffness values of the virtual springs in the vertical direction are 4000N/m~5000N/m; the value range of the preset stiffness value of the virtual spring in the rotation direction with the vertical axis as the axis is 0Nm/rad~20Nm/rad; the preset stiffness value of the virtual spring in the rotation direction with the depth direction as the axis It is assumed that the value ranges of the stiffness value and the stiffness value of the virtual spring in the rotation direction with the horizontal axis as the axis are both 200Nm/rad-300Nm/rad.
在一个实施例中,配准模块304还被配置为:In one embodiment, registration module 304 is further configured to:
根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行采点操作时的术中标记点在世界坐标系的空间位置;获取骨骼的三维模型中骨骼上的术前规划点在三维模型坐标系的空间位置;将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准,得到粗配准矩阵;根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行划线操作时的划线点集在世界坐标系的空间位置;根据粗配准矩阵将划线点集在世界坐标系下的空间位置与三维模型进行精配准,得到配准结果。According to the spatial position of the surgical probe in the world coordinate system obtained by the tracking camera, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone; obtain the bone in the 3D model of the bone The spatial position of the preoperative planning point on the 3D model coordinate system; the spatial position of the preoperative planning point in the 3D model coordinate system is roughly registered with the spatial position of the intraoperative marker point in the world coordinate system, and the rough registration is obtained. Quasi-matrix; according to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the marking point set in the world coordinate system when the surgical probe is scribing on the solid bone; according to the coarse registration The matrix precisely registers the spatial position of the line point set in the world coordinate system with the 3D model to obtain the registration result.
关于手术机器人追踪和移动的系统的限定可以参见上文中对于手术机器人追踪和移动的方法的限定,在此不再赘述。上述手术机器人追踪和移动的系统中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。For the definition of the system for tracking and moving the surgical robot, please refer to the definition of the method for tracking and moving the surgical robot above, and details will not be repeated here. Each module in the above-mentioned system for tracking and moving the surgical robot can be fully or partially realized by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, and can also be stored in the memory of the computer device in the form of software, so that the processor can invoke and execute the corresponding operations of the above-mentioned modules.
一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器被配置为提供计算和控制能力。该计算 机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库被配置为存储一种手术机器人追踪和移动的方法的数据。该计算机设备的网络接口被配置为与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种手术机器人追踪和移动的方法。In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure may be as shown in FIG. 4 . The computer device includes a processor, memory, network interface and database connected by a system bus. Wherein, the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer equipment includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, computer programs and databases. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database of the computer device is configured to store data of a method of tracking and movement of a surgical robot. The network interface of the computer device is configured to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a method for tracking and moving a surgical robot is realized.
本领域技术人员可以理解,图4中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。Those skilled in the art can understand that the structure shown in Figure 4 is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation to the computer equipment on which the solution of the application is applied. The specific computer equipment can be More or fewer components than shown in the figures may be included, or some components may be combined, or have a different arrangement of components.
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,该存储器存储有计算机程序,该处理器执行计算机程序时实现上述各个实施例中的步骤。In one embodiment, a computer device is provided, including a memory and a processor, the memory stores a computer program, and the processor implements the steps in the foregoing embodiments when executing the computer program.
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现上述各个实施例中的步骤。In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the steps in the foregoing embodiments are implemented.
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be implemented through computer programs to instruct related hardware, and the computer programs can be stored in a non-volatile computer-readable memory In the medium, when the computer program is executed, it may include the processes of the embodiments of the above-mentioned methods. Wherein, any references to memory, storage, database or other media used in the various embodiments provided in the present application may include non-volatile and/or volatile memory. Nonvolatile memory can include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Chain Synchlink DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above embodiments can be combined arbitrarily. To make the description concise, all possible combinations of the technical features in the above embodiments are not described. However, as long as there is no contradiction in the combination of these technical features, they should be It is considered to be within the range described in this specification.
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对公开专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation modes of the present application, and the descriptions thereof are relatively specific and detailed, but should not be construed as limiting the scope of the disclosed patents. It should be noted that those skilled in the art can make several modifications and improvements without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the scope of protection of the patent application should be based on the appended claims.

Claims (10)

  1. 一种手术机器人追踪和移动的方法,包括:A method for tracking and moving a surgical robot, comprising:
    对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果;Register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located with the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain a registration result;
    根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置;According to the position of the tracker on the end of the manipulator and the tracker on the bone acquired by the tracking camera in real time, determine the real-time spatial position of the actuator at the end of the manipulator and the bone in the world coordinate system, and convert them according to the registration result Go to the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
    根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。According to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, the mechanical arm is controlled to limit the movement of the actuator within the current target area.
  2. 根据权利要求1所述的方法,其中,所述方法还包括:The method according to claim 1, wherein the method further comprises:
    在执行器运行前,当机械臂运动至骨骼处时,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定执行器的调节路径;Before the actuator runs, when the mechanical arm moves to the bone, the adjustment path of the actuator is determined according to the current spatial position of the actuator and the current target area in the three-dimensional solid model;
    在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。The adjustment path is displayed in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is coplanar with the current target area.
  3. 根据权利要求2所述的方法,其中,根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内的步骤包括:The method according to claim 2, wherein, according to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, the step of controlling the mechanical arm to limit the movement of the actuator within the current target area comprises:
    在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行执行器;After the plane of the actuator is coplanar with the current target area, when the robot arm is operated, run the actuator;
    启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。Start the Cartesian damping control mode modeled on virtual springs and dampers. The manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions The offset Δx outputs the feedback force F opposite to the operated direction, F=Δx*C, so as to limit the movement of the actuator within the current target area.
  4. 根据权利要求3所述的方法,其中,执行器向当前目标区域切入的方向记为进深方向,在当前目标区域内并与切入方向垂直的方向记为横向,与当前目标区域垂直的方向记为垂直方向;偏移量包括进深方向的偏移值、横向的偏移值、垂直方向的偏移值、绕进深方向旋转的偏移值、绕横向旋转的偏移值、绕垂直方向旋转的偏移值;The method according to claim 3, wherein the direction in which the actuator cuts into the current target area is marked as the depth direction, the direction within the current target area and perpendicular to the cutting direction is marked as the horizontal direction, and the direction perpendicular to the current target area is marked as Vertical direction; the offset includes the offset value in the depth direction, the offset value in the horizontal direction, the offset value in the vertical direction, the offset value rotated around the depth direction, the offset value rotated around the horizontal direction, the offset value rotated around the vertical direction transfer value;
    进深方向上的虚拟弹簧的预设刚度值、横向上的虚拟弹簧的预设刚度值的取值范围均为0N/m~500N/m;The preset stiffness value of the virtual spring in the depth direction and the preset stiffness value of the virtual spring in the lateral direction are all in the range of 0N/m to 500N/m;
    垂直方向上的虚拟弹簧的预设刚度值的取值范围为4000N/m~5000N/m;The preset stiffness value of the virtual spring in the vertical direction ranges from 4000N/m to 5000N/m;
    以垂直方向为轴旋转方向上的虚拟弹簧的预设刚度值的取值范围为0Nm/rad~20 Nm/rad;The preset stiffness value of the virtual spring in the direction of rotation with the vertical direction as the axis ranges from 0Nm/rad to 20 Nm/rad;
    以进深方向为轴旋转方向上的虚拟弹簧的预设刚度值、以横向为轴旋转方向上的虚拟弹簧的刚度值的取值范围均为200Nm/rad~300Nm/rad。The preset stiffness value of the virtual spring in the direction of rotation with the depth direction as the axis and the stiffness value of the virtual spring in the direction of rotation with the horizontal axis as the axis both range from 200Nm/rad to 300Nm/rad.
  5. 根据权利要求1所述的方法,其中,对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系的步骤包括:The method according to claim 1, wherein the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located is registered with the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system The steps include:
    根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行采点操作时的术中标记点在世界坐标系的空间位置;According to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the intraoperative marker point in the world coordinate system when the surgical probe is used to collect points on the solid bone;
    获取骨骼的三维模型中骨骼上的术前规划点在三维模型坐标系的空间位置;Obtain the spatial position of the preoperative planning point on the bone in the three-dimensional model of the bone in the three-dimensional model coordinate system;
    将术前规划点在三维模型坐标系下的空间位置与术中标记点在世界坐标系下的空间位置进行粗配准,得到粗配准矩阵;Roughly register the spatial position of the preoperative planning point in the 3D model coordinate system with the spatial position of the intraoperative marker point in the world coordinate system to obtain a coarse registration matrix;
    根据追踪相机获取的手术探针在世界坐标系的空间位置,确定手术探针在实体的骨骼上进行划线操作时的划线点集在世界坐标系的空间位置;According to the spatial position of the surgical probe acquired by the tracking camera in the world coordinate system, determine the spatial position of the marking point set in the world coordinate system when the surgical probe performs a marking operation on the solid bone;
    根据粗配准矩阵将划线点集在世界坐标系下的空间位置与三维模型进行精配准,得到配准结果。According to the coarse registration matrix, the spatial position of the line point set in the world coordinate system is finely registered with the 3D model, and the registration result is obtained.
  6. 一种手术机器人追踪和移动系统,包括:A surgical robot tracking and movement system comprising:
    配准模块,被配置为对预先获取的骨骼的三维模型所在的三维模型坐标系与实体的骨骼所在的世界坐标系进行配准,以将世界坐标系配准到三维模型坐标系,得到配准结果;The registration module is configured to register the 3D model coordinate system where the pre-acquired 3D model of the skeleton is located and the world coordinate system where the skeleton of the entity is located, so as to register the world coordinate system to the 3D model coordinate system to obtain registration result;
    追踪模块,被配置为根据追踪相机实时获取的机械臂末端上的示踪器、骨骼上的示踪器的位置,确定机械臂末端的执行器、骨骼在世界坐标系下的实时空间位置,根据配准结果将其转换到三维模型坐标系中,获得执行器、骨骼在三维模型坐标系中的空间位置;The tracking module is configured to determine the real-time spatial position of the actuator at the end of the mechanical arm and the bone in the world coordinate system according to the position of the tracker on the end of the mechanical arm and the tracker on the bone acquired by the tracking camera in real time, according to The registration result is transformed into the 3D model coordinate system to obtain the spatial position of the actuator and bones in the 3D model coordinate system;
    运动控制模块,被配置为根据当前目标区域的空间位置、机械臂末端的执行器的空间位置,控制机械臂,以将执行器的运动限定在当前目标区域内。The motion control module is configured to control the mechanical arm according to the spatial position of the current target area and the spatial position of the actuator at the end of the mechanical arm, so as to limit the movement of the actuator to the current target area.
  7. 根据权利要求6所述的系统,其中,所述系统还包括:The system according to claim 6, wherein said system further comprises:
    平面对齐模块,被配置为在执行器运行前,当机械臂运动至骨骼处时,根据执行器、当前目标区域在三维实体模型中的当前空间位置确定执行器的调节路径;The plane alignment module is configured to determine the adjustment path of the actuator according to the current spatial position of the actuator and the current target area in the three-dimensional solid model when the mechanical arm moves to the bone before the actuator runs;
    显示模块,被配置为在三维模型中显示调节路径,以使操作者根据调节路径操作机械臂,使机械臂带动执行器运动至当前目标区域的外缘,使执行器的平面与当前目标区域共面。The display module is configured to display the adjustment path in the three-dimensional model, so that the operator operates the mechanical arm according to the adjustment path, so that the mechanical arm drives the actuator to move to the outer edge of the current target area, so that the plane of the actuator is the same as the current target area. noodle.
  8. 根据权利要求6所述的系统,其中,所述运动控制模块还被配置为:The system of claim 6, wherein the motion control module is further configured to:
    在执行器的平面与当前目标区域共面之后,在机械臂被操作时,运行所述执行器;After the plane of the actuator is coplanar with the current target area, when the mechanical arm is operated, the actuator is operated;
    启动以虚拟弹簧和阻尼器为模型的笛卡尔阻尼控制模式,机械臂基于多个自由度方向上的各个虚拟弹簧的预设刚度值C和多个自由度方向上执行器相对于当前目标区域的偏移量Δx,输出与被操作方向相反的反馈力F,F=Δx*C,从而将执行器的运动限定在当前目标区域内。Start the Cartesian damping control mode modeled on virtual springs and dampers. The manipulator is based on the preset stiffness value C of each virtual spring in multiple degrees of freedom directions and the actuator relative to the current target area in multiple degrees of freedom directions The offset Δx outputs the feedback force F opposite to the operated direction, F=Δx*C, so as to limit the movement of the actuator within the current target area.
  9. 一种计算机设备,包括存储器和处理器,所述存储器存储有可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的方法的步骤。A computer device, comprising a memory and a processor, the memory stores a computer program that can run on the processor, and the processor implements the method according to any one of claims 1 to 5 when executing the computer program A step of.
  10. 一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。A computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are realized.
PCT/CN2022/106427 2021-09-03 2022-07-19 Surgery robot tracking and moving method and system WO2023029784A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111035623.9 2021-09-03
CN202111035623.9A CN113855236B (en) 2021-09-03 2021-09-03 Method and system for tracking and moving surgical robot

Publications (1)

Publication Number Publication Date
WO2023029784A1 true WO2023029784A1 (en) 2023-03-09

Family

ID=78989727

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/106427 WO2023029784A1 (en) 2021-09-03 2022-07-19 Surgery robot tracking and moving method and system

Country Status (2)

Country Link
CN (1) CN113855236B (en)
WO (1) WO2023029784A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117464691A (en) * 2023-12-27 2024-01-30 中国科学院自动化研究所 Rebound prevention control method, rebound prevention control device, electronic equipment and storage medium
CN117770958A (en) * 2023-12-26 2024-03-29 杭州三坛医疗科技有限公司 Tracer orientation positioning method and device, electronic equipment and medium

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113855236B (en) * 2021-09-03 2022-05-31 北京长木谷医疗科技有限公司 Method and system for tracking and moving surgical robot
CN113842214B (en) * 2021-09-03 2022-10-11 北京长木谷医疗科技有限公司 Surgical robot navigation positioning system and method
CN113842211B (en) * 2021-09-03 2022-10-21 北京长木谷医疗科技有限公司 Three-dimensional preoperative planning system for knee joint replacement and prosthesis model matching method
CN113842213B (en) * 2021-09-03 2022-10-11 北京长木谷医疗科技有限公司 Surgical robot navigation positioning method and system
CN113842217B (en) * 2021-09-03 2022-07-01 北京长木谷医疗科技有限公司 Method and system for limiting motion area of robot
CN114533277B (en) * 2022-02-24 2024-07-05 上海电气集团股份有限公司 Control method of mechanical arm, electronic equipment and storage medium
CN115381554B (en) * 2022-08-02 2023-11-21 北京长木谷医疗科技股份有限公司 Intelligent position adjustment system and method for orthopedic surgery robot
CN115444555B (en) * 2022-08-31 2024-09-06 北京长木谷医疗科技股份有限公司 Intelligent positioning tracking method and device for orthopedic operation robot and electronic equipment
CN116077182B (en) * 2022-12-23 2024-05-28 北京纳通医用机器人科技有限公司 Medical surgical robot control method, device, equipment and medium
CN117159144B (en) * 2023-09-06 2024-06-28 北京长木谷医疗科技股份有限公司 Navigation and positioning system and method for knee joint ligament reconstruction operation robot
CN116983086B (en) * 2023-09-26 2024-01-09 北京长木谷医疗科技股份有限公司 Autonomous joint replacement surgery robot navigation positioning system

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104146767A (en) * 2014-04-24 2014-11-19 薛青 Intraoperative navigation method and system for assisting in surgery
CN104463832A (en) * 2013-09-25 2015-03-25 中国科学院深圳先进技术研究院 Registration method and system of neurosurgery navigation
US20180235715A1 (en) * 2017-02-22 2018-08-23 Orthosoft Inc. Bone and tool tracking in robotized computer-assisted surgery
CN109688963A (en) * 2016-07-15 2019-04-26 马科外科公司 System for robot assisted revision program
CN110353806A (en) * 2019-06-18 2019-10-22 北京航空航天大学 Augmented reality navigation methods and systems for the operation of minimally invasive total knee replacement
CN111084661A (en) * 2018-10-23 2020-05-01 株式会社卓越牵引力 Operation support device, method of controlling the same, and recording medium
CN111345896A (en) * 2020-03-13 2020-06-30 北京天智航医疗科技股份有限公司 Osteotomy execution system, positioning, control and simulation execution method and electronic equipment
CN113855236A (en) * 2021-09-03 2021-12-31 北京长木谷医疗科技有限公司 Method and system for tracking and moving surgical robot

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110306873A1 (en) * 2010-05-07 2011-12-15 Krishna Shenai System for performing highly accurate surgery
CN110946654B (en) * 2019-12-23 2022-02-08 中国科学院合肥物质科学研究院 Bone surgery navigation system based on multimode image fusion
CN111388087A (en) * 2020-04-26 2020-07-10 深圳市鑫君特智能医疗器械有限公司 Surgical navigation system, computer and storage medium for performing surgical navigation method
CN112603538A (en) * 2021-02-23 2021-04-06 南京佗道医疗科技有限公司 Orthopedic navigation positioning system and method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104463832A (en) * 2013-09-25 2015-03-25 中国科学院深圳先进技术研究院 Registration method and system of neurosurgery navigation
CN104146767A (en) * 2014-04-24 2014-11-19 薛青 Intraoperative navigation method and system for assisting in surgery
CN109688963A (en) * 2016-07-15 2019-04-26 马科外科公司 System for robot assisted revision program
US20180235715A1 (en) * 2017-02-22 2018-08-23 Orthosoft Inc. Bone and tool tracking in robotized computer-assisted surgery
CN111084661A (en) * 2018-10-23 2020-05-01 株式会社卓越牵引力 Operation support device, method of controlling the same, and recording medium
CN110353806A (en) * 2019-06-18 2019-10-22 北京航空航天大学 Augmented reality navigation methods and systems for the operation of minimally invasive total knee replacement
CN111345896A (en) * 2020-03-13 2020-06-30 北京天智航医疗科技股份有限公司 Osteotomy execution system, positioning, control and simulation execution method and electronic equipment
CN113855236A (en) * 2021-09-03 2021-12-31 北京长木谷医疗科技有限公司 Method and system for tracking and moving surgical robot

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117770958A (en) * 2023-12-26 2024-03-29 杭州三坛医疗科技有限公司 Tracer orientation positioning method and device, electronic equipment and medium
CN117464691A (en) * 2023-12-27 2024-01-30 中国科学院自动化研究所 Rebound prevention control method, rebound prevention control device, electronic equipment and storage medium
CN117464691B (en) * 2023-12-27 2024-03-22 中国科学院自动化研究所 Rebound prevention control method, rebound prevention control device, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN113855236A (en) 2021-12-31
CN113855236B (en) 2022-05-31

Similar Documents

Publication Publication Date Title
WO2023029784A1 (en) Surgery robot tracking and moving method and system
WO2023029362A1 (en) Surgical robot navigation and positioning method and apparatus, and device and storage medium
CN113842214B (en) Surgical robot navigation positioning system and method
US11771498B2 (en) Lower extremities leg length calculation method
US9916421B2 (en) Implant planning using corrected captured joint motion information
CN113850851B (en) Surgical robot bone registration method and system
US9913692B2 (en) Implant planning using captured joint motion information
US20230024464A1 (en) Systems and methods for adaptive planning and control of a surgical tool
US10194990B2 (en) Method for augmenting a surgical field with virtual guidance content
US20200038112A1 (en) Method for augmenting a surgical field with virtual guidance content
CN109925055A (en) Totally digitilized total knee replacement surgical robot system and its simulation operation method
WO2023029922A1 (en) Method and system for limiting movement area of robot
US20220183768A1 (en) Systems and methods for initial assessment warnings
CN116983086B (en) Autonomous joint replacement surgery robot navigation positioning system
US20240216071A1 (en) Surgical registration tools, systems, and methods of use in computer-assisted surgery
Zhang et al. SLAM-TKA: Real-time Intra-operative Measurement of Tibial Resection Plane in Conventional Total Knee Arthroplasty

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22862925

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE