[go: up one dir, main page]

CN111134842A - Robot-assisted fracture reduction path planning method - Google Patents

Robot-assisted fracture reduction path planning method Download PDF

Info

Publication number
CN111134842A
CN111134842A CN202010061853.1A CN202010061853A CN111134842A CN 111134842 A CN111134842 A CN 111134842A CN 202010061853 A CN202010061853 A CN 202010061853A CN 111134842 A CN111134842 A CN 111134842A
Authority
CN
China
Prior art keywords
fracture
reduction
robot
path planning
path
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010061853.1A
Other languages
Chinese (zh)
Inventor
雷静桃
程利亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Transpacific Technology Development Ltd
Original Assignee
Beijing Transpacific Technology Development Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Transpacific Technology Development Ltd filed Critical Beijing Transpacific Technology Development Ltd
Priority to CN202010061853.1A priority Critical patent/CN111134842A/en
Publication of CN111134842A publication Critical patent/CN111134842A/en
Pending legal-status Critical Current

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/10Computer-aided planning, simulation or modelling of surgical operations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • 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/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

Landscapes

  • Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Robotics (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Manipulator (AREA)

Abstract

本发明涉及一种机器人辅助骨折复位手术的路径规划方法,采用骨折远端姿态调整及基于A*算法的双平面避碰复位路径规划方法。在mimics和Geomagic studio软件中进行基于健侧镜像匹配方法确定骨折偏差参数,在两个不同坐标平面内分别进行基于A*算法的骨折复位路径规划。本发明所述的骨折复位路径规划方法,能有效避免复位中远端的大幅度牵拉以及骨折远端与周围组织碰撞,实现短耗时、高精度、安全避碰的骨折复位操作。具有方法简单、安全有效等优点,适用于机器人辅助股骨、胫骨等长骨骨折复位手术,能够提高骨折复位操作的精度和安全性,提高复位效率,降低医生工作强度。

Figure 202010061853

The invention relates to a path planning method for robot-assisted fracture reduction surgery, which adopts the attitude adjustment of the distal end of the fracture and a biplane collision avoidance reduction path planning method based on the A* algorithm. In mimics and Geomagic studio software, the fracture deviation parameters were determined based on the mirror matching method of the unaffected side, and the fracture reduction path planning based on the A* algorithm was carried out in two different coordinate planes respectively. The fracture reduction path planning method of the invention can effectively avoid the large pulling of the middle and distal end of the reduction and the collision of the distal end of the fracture with the surrounding tissue, and realize the fracture reduction operation with short time consumption, high precision and safe collision avoidance. It has the advantages of simple, safe and effective method, and is suitable for robot-assisted reduction of long bone fractures such as femur and tibia, which can improve the accuracy and safety of fracture reduction operation, improve reduction efficiency, and reduce the work intensity of doctors.

Figure 202010061853

Description

Robot-assisted fracture reduction path planning method
Technical Field
The invention belongs to the technical field of robots, and relates to a path planning method for a robot-assisted fracture reduction operation.
Background
The operation reduction is a means of fracture surgical treatment, and the traditional fracture reduction operation is performed under the guidance of an X-ray machine by a doctor and is difficult to achieve high precision. The minimally invasive internal fixation for reduction of closed fracture treats femoral shaft fracture, and has the advantages of small operative wound, quick fracture healing, high healing rate, low infection incidence and the like. However, reduction of fractured ends, particularly rotational reduction, has been a surgical challenge because the fractured ends are not directly exposed.
With the development of the robot technology, the robot-assisted fracture reduction research is developed at home and abroad. The supplementary fracture operation that resets of robot, doctor operation robot observe the condition of resetting of patient's fracture position through the X ray machine, and the accurate reasonable route planning that resets before the art guarantees the prerequisite of operation safety, has the significance to improving the operation precision, reducing patient's operation damage.
In the process of robot-assisted fracture reduction, the following points need to be considered:
(1) during reduction, the distal end of the fracture cannot collide with the tissue surrounding the fracture. After the fracture, traction forces caused by the impingement of the musculoskeletal tissues and surrounding soft tissues impede the desired reduction motion. If no rational path planning is performed before surgery, the nervous system, soft tissues and muscle tissues are easily stretched and distorted, causing secondary injuries to the patient, wherein the damage to the nervous system is the most serious and the damage is likely to have irreversible consequences. In order to avoid secondary damage possibly caused by the robot reset operation, improve the safety and the precision of the robot auxiliary reset operation and shorten the operation time, a reasonable and safe collision-preventing reset path is necessary to be planned before the operation.
(2) For closed fractures, there is usually a deviation in displacement as well as a deviation in angular rotation. Under the condition that the internal part around the fracture is difficult to observe, preoperative path planning is needed, and accurate and safe reduction operation on the fracture is realized.
(3) In robotic-assisted surgery, the limited working space and the possibility of collisions between the end of the robotic manipulator arm and the surrounding environment present challenges in robotic-assisted surgery where the end effector of the surgical robot must reach the anatomical target of the patient and not collide with the patient or surrounding instruments.
The path planning refers to searching and finding a collision-free path from a starting point to a target point in an environment with obstacles according to certain evaluation criteria (such as distance, time, cost and the like), and the core of the path planning is a path planning algorithm. The traditional path planning algorithm comprises an artificial potential field method, a fuzzy logic algorithm, a grid method, a free space method and the like. With the development of artificial intelligence technology, some intelligent algorithms for path planning appear, such as ant colony algorithm, genetic algorithm, neural network algorithm, and the like.
The preoperative collision avoidance reduction path planning is to determine a moving path when a robot assists the fracture reduction operation. On the basis, an operation navigation system is utilized to carry out accurate resetting operation, the broken and dislocated bones after the fracture are restored to the normal anatomical positions, and dislocation does not occur in the bone healing process.
For long bone fracture types of femur and tibia, after fracture, axial offset, longitudinal offset, radial (or transverse) offset and rotational offset of a far end relative to a near end occur, and 6 fracture deviation parameters exist. Aiming at the closed reduction of the long bone fracture, a robot is adopted to assist the operation of the fracture reduction operation, the near end of the fracture is fixed on an operating bed, and the robot assists the reduction operation of aligning the far end of the fracture to the near end of the fracture. In the prior art, when the robot-assisted operation is used for resetting, medical accidents are easy to happen due to the large-amplitude traction of the far end and the collision of the fracture far end and surrounding tissues, and the precision and the safety are insufficient.
Disclosure of Invention
The invention relates to a collision-prevention path planning method for a robot-assisted fracture reduction operation, which is characterized in that a fracture near end is fixed, the tail end of a robot is firmly connected with a fracture far end, and the robot drives the fracture far end to move in a collision-prevention manner relative to the fracture near end to perform reduction operation until fracture sections of the far end and the near end are butted to a correct anatomical position, so that robot-assisted fracture reduction is completed. The invention relates to a path planning method for robot-assisted fracture reduction operation, in particular to a fracture distal attitude adjustment and biplane collision prevention reduction path planning method based on an A-algorithm. The method for planning the fracture reduction path can effectively avoid the large-amplitude traction of the far end in the reduction and the collision of the fracture far end and surrounding tissues, and realize the fracture reduction operation with short time consumption, high precision and safe collision avoidance. The preoperative path planning method for robot-assisted fracture reduction has the advantages of simplicity, safety, reliability and the like, and is suitable for not only robot-assisted femoral fracture reduction operations, but also tibial and long-bone fracture reduction operations. The fracture reduction robot can effectively improve the precision and the safety of the robot for assisting fracture reduction operation, improve the reduction efficiency and reduce the damage of ray radiation to doctors and patients.
In order to achieve the above purpose, the technical scheme of the path planning method of the invention comprises the following steps:
a path planning method for a robot-assisted fracture reduction operation comprises the following steps:
1) aiming at the closed reduction of the long bone fracture, when a robot is adopted to assist the operation of the fracture reduction operation, the fracture near end is fixed, and the robot assists to align the fracture far end to the near end for the reduction operation;
2) before an operation, according to CT scanning data of a fracture side and a healthy side bone of a patient before the operation provided by a doctor, a reverse modeling method is adopted to reconstruct and obtain a three-dimensional digital model of the fracture side and the healthy side bone;
3) based on the principle of reduction of healthy side mirror image registration, determining axial, longitudinal, radial or transverse displacement deviation and rotation deviation of the fracture in mimics and Geomagic studio software by an image registration method, and analyzing the displacement and rotation deviation values of the fracture; determining a fracture deviation value as a basis for planning a reduction path;
4) adjusting the fracture far-end space posture:
the posture of the far end of the fracture relative to the near end is adjusted by a robot, so that the far end of the fracture achieves a reasonable anatomical posture relative to the near end of the fracture, and meanwhile, 3 posture angles of the far end of the fracture are tested in real time by an NDI optical tracking array fixed on the far end of the fracture, so that the posture adjustment of the far end of the fracture and the requirement on the reduction precision after the reduction is finished are met;
5) reset path planning preparation:
first, a grid map of the fracture model is given, the number and spacing of the grids obeying two conditions:
(1) the grid spacing is smaller than the reset precision;
(2) accurately representing the shape of the obstacle on the premise of meeting the precision;
then, the lower half parts of the fracture far end and the fracture near end are similar to cylinders, and the parts of the near end and the far end which do not influence the reset path are abandoned;
6) planning collision avoidance paths:
determining an optimal collision avoidance reduction path by using a fracture reduction path planning algorithm based on an A-x algorithm;
7) planning a reset path:
respectively projecting the three-dimensional digital model of the long bone fracture into two rectangular coordinate planes, and respectively planning a reduction path in the two planes by adopting an A-x algorithm;
8) coordinates of data points obtained by planning the reset path based on the A-x algorithm are stored as data files and transmitted to a reset robot controller, and the robot controls the fracture distal end to perform reset operation according to the planned reset path;
9) fixing NDI optical tracking system navigation frames at the far end and the near end of the fracture respectively, and testing the spatial poses of the near end and the far end of the fracture; and the near end and the far end of the fracture are respectively fixed with an optical tracking array of an NDI positioning system, and the relative spatial postures of the near end and the far end of the fracture are tested.
As a preferred technical scheme of the invention, in the step 6), a biplane collision avoidance reduction path planning method based on an a-x algorithm is adopted to determine a fracture reduction path; and projecting the proximal end and the distal end of the fracture into two rectangular coordinate planes, and planning a reduction path in the two coordinate planes by respectively adopting an A-x algorithm.
As a preferred technical solution of the present invention, in the step 6), the a-algorithm determines the search direction through an estimation function, expands from the starting point to the periphery, obtains the cost value of each peripheral node through calculation of a valuation function, selects the minimum cost node as the next expansion node, repeats this process until reaching the target point, and generates the final reset path.
As a preferable technical solution of the present invention, in the step 1), when performing the closed reduction operation with the robot assistance for the long bone fracture of femur or tibia, the proximal end of the fracture is fixed to the operating bed, and the robot assistance is performed for the reduction operation in which the distal end of the fracture is aligned with the proximal end of the fracture.
As a preferable technical scheme of the present invention, in the step 4), the navigation frame is fixed to the fracture distal end, the robot adjusts the posture of the fracture distal end according to the fracture rotation deviation value, and observes the posture angle between the navigation frame and the NDI reference coordinate system in real time until the posture of the fracture distal end is adjusted in place.
As a preferred technical solution of the present invention, in the step 5), in the fracture model projection diagram projected into the rectangular coordinate plane, a series of feature points including a start point and a stop point for planning a fracture reduction path and a plurality of obstacle points on the path are determined, where the obstacle points are used as obstacle avoidance constraint points for planning the reduction path, so as to ensure that collision between the distal end and the proximal end of the femur during reduction operation should be avoided.
Before the path planning is carried out by adopting the A-star algorithm, a grid graph of a fracture model needs to be given. The number and spacing of the grids follows two conditions: the grid spacing is smaller than the reset precision; the shape of the obstacle is accurately represented on the premise of meeting the precision. The lower halves of the fracture distal and proximal ends are approximated to be cylindrical, leaving away portions of the proximal and distal ends that do not affect the reduction path. Determining the starting point and the ending point of the fracture reduction path planning and the characteristic points of a plurality of obstacle points on the path. The soft tissue balance and the collision avoidance of the tissues around the fractured bone during the reduction operation are considered. Because the tissues around the fracture are complex, in order to avoid secondary injury in the reduction operation, the collision between the far end and the near end of the femur is avoided during the reduction operation, and a series of obstacle points are set in the projection drawing of the fracture model and are used as obstacle avoidance constraint points for planning the reduction path.
As a preferred embodiment of the present invention, in the step 6), a fracture reduction path planning algorithm based on the a-algorithm is used. The A-star algorithm based on the grid method can determine the optimal collision avoidance reset path more quickly and effectively. When the path planning is performed on the grid graph by the algorithm A, the center of each grid is assumed to be a node, so that the number of adjacent nodes of each node is limited to 8, and the moving direction angle of each node is also limited to integral multiple.
As a preferred technical scheme of the invention, the A-algorithm determines the search direction through an estimation function, expands from a starting point to the periphery, obtains the cost value of each node around through calculation of the evaluation function, selects the minimum cost node as the next expansion node, repeats the process until reaching the target point, and generates the final reset path.
As a preferred technical solution of the present invention, in the search process, each node on the path is a node with the minimum cost, and therefore the obtained path cost is the minimum. The valuation function of the a-algorithm is: . Wherein, the evaluation function which represents the arrival of the target point from the starting point through any node represents the actual cost from the starting point to the node and represents the nodenEstimated cost to target point.
Manhattan distance: the sum of the distances of the projections generated by the line segment formed by two points on the fixed rectangular coordinate system of the Euclidean space to the axis is the result of summing the distances in multiple dimensions. For example, the manhattan distance between nodes on a plane is as follows. Euclidean distance: the euclidean distance in two-dimensional space is the distance between two points. Such as euclidean distance between nodes on a plane.
As a further preferable technical solution of the present invention, in the step 6), when calculating the valuation function, a calculation formula of manhattan distance and euclidean distance is used; when the nodes are searched along the side length direction of the grid graph, a Manhattan distance formula is adopted for calculation; when the nodes are searched along the diagonal direction of the grid graph, the Euclidean distance formula is adopted for calculation.
Through repeated tests, the fracture reduction collision avoidance path planning method successfully solves the problem of a robot-assisted fracture reduction operation path.
In summary, the path planning method for the robot-assisted fracture reduction surgery operation, namely the method for adjusting the posture of the far end of the fracture and planning the biplane collision prevention reduction path based on the a-algorithm, determines the fracture deviation parameters based on the healthy side mirror image matching method in the mimics and the Geomagic studio software, and plans the fracture reduction path based on the a-algorithm in two different coordinate planes respectively. The method for planning the fracture reduction path can effectively avoid the large-amplitude traction of the far end in the reduction and the collision of the fracture far end and surrounding tissues, and realize the fracture reduction operation with short time consumption, high precision and safe collision avoidance. The method has the advantages of simplicity, safety, effectiveness and the like, is suitable for the robot-assisted femoral fracture reduction operation, can improve the precision and safety of the fracture reduction operation, improves the reduction efficiency, and reduces the damage of ray radiation to doctors and patients.
Compared with the prior art, the invention has the following obvious prominent substantive characteristics and obvious advantages:
1. the invention relates to a path planning method for robot-assisted fracture reduction operation, in particular to a method for adjusting the posture of a fracture far end and planning a biplane collision prevention reduction path based on an A-x algorithm.
2. The invention is based on the principle of reduction of healthy side mirror images, and adopts a fracture deviation determination method based on image matching. And performing healthy side mirror image matching in mimics and Geomagic studio software to determine fracture deviation parameters.
3. According to the biplane reset path planning method, the reset path planning is carried out in two rectangular coordinate planes by adopting an A-algorithm respectively, and the optimal collision avoidance reset path can be determined more quickly and effectively by the A-algorithm based on the grid method.
4. The invention comprehensively considers the collision avoidance of the fracture far end and the surrounding tissues in the fracture reduction operation, and effectively improves the safety and the accuracy of the robot-assisted reduction operation.
5. The method has certain universality, and is suitable for femoral fractures and other long bone fracture types such as tibia and the like.
6. The fracture reduction path planning method has the advantages of being simple, safe, effective and the like, can effectively improve the precision and safety of fracture reduction operation, improve the reduction efficiency and reduce the damage of ray radiation to doctors and patients.
Drawings
FIG. 1 is a diagram of a bone model of a fractured side and a healthy side reconstructed based on a reverse modeling technique according to the present invention.
FIG. 2 is a schematic diagram of fracture deviation determination based on healthy side mirror image registration according to the present invention.
Fig. 3 is a projection schematic view of the fracture model of the invention in the yoz plane.
FIG. 4 is a schematic representation of the collision threshold of the proximal and distal fractures of the present invention.
Fig. 5 is a schematic diagram of a grid of the fracture model of the present invention in the yoz plane.
Fig. 6 is a schematic diagram of the planning result of the reset path of the yoz plane according to the present invention.
FIG. 7 is a schematic diagram of the biplane reduction path planning result of the fracture model of the present invention.
Detailed Description
The invention is further described with reference to specific embodiments and the accompanying drawings.
The embodiment is a patient femoral fracture. Reasonably arranging an operating bed, a reset robot, a fracture model, an NDI optical tracking system and the like. Two optical tracking arrays of the NDI optical tracking system are respectively fixed at the near end and the far end of the fracture and used for testing the spatial position and the posture of the near end and the far end of the fracture. And adjusting a reference coordinate system of the NDI positioning system, and adjusting the initial pose of the resetting robot to prepare for resetting operation.
In this embodiment, a path planning method for a robot-assisted fracture reduction operation includes the following steps:
1) aiming at the closed reduction of the long bone fracture, when a robot is adopted to assist the operation of the fracture reduction operation, the fracture near end is fixed, and the robot assists to align the fracture far end to the near end for the reduction operation; when the robot assists in fracture reduction operation, the near end of the fracture is fixed on an operating bed, the far end of the fracture is fixedly connected with the tail end of the robot, and the robot moves along with the robot, namely the robot moves the far end of the fracture relative to the near end of the fracture in a collision-preventing manner to perform reduction operation until the fracture sections of the far end and the near end are butted to a correct anatomical position;
2) before an operation, according to CT scanning data of a fracture side and a healthy side bone of a patient before the operation provided by a doctor, a reverse modeling method is adopted to reconstruct and obtain a three-dimensional digital model of the fracture side and the healthy side bone; as shown in fig. 1;
3) and (3) determining fracture deviation by software registration:
according to the bone registration principle of the fracture side and the healthy side based on the healthy side mirror image, in the mimics and the Geomagic student software, the axial, longitudinal, radial or transverse displacement deviation and rotation deviation of the fracture are determined by an image registration method, the displacement and rotation deviation values of the fracture are analyzed, and the registration is carried out to determine the fracture deviation, as shown in fig. 2;
4) the robot assists in finishing the posture adjustment of the far end of the fracture relative to the near end of the fracture, so that an ideal anatomical posture can be achieved after the fracture is reset; the posture of the far end of the fracture relative to the near end is adjusted by a robot, so that the far end of the fracture achieves a reasonable anatomical posture relative to the near end of the fracture, and meanwhile, 3 posture angles of the far end of the fracture are tested in real time by an NDI optical tracking array fixed on the far end of the fracture, so that the posture adjustment of the far end of the fracture and the requirement on the reduction precision after the reduction is finished are met;
5) reset path planning preparation:
and (4) planning a biplane collision avoidance reset path based on an A-star algorithm. The algorithm determines the searching direction through an evaluation function, expands from the starting point to the periphery, calculates the cost value of each node around through the evaluation function, selects the minimum cost node as the next expansion node, repeats the process until the end point is reached, and generates the final reset path;
planning a reset path based on an A-x algorithm in a rectangular coordinate plane; the lower half parts of the distal femur and the proximal femur are similar to cylinders, and the parts of the proximal femur and the distal femur, which do not influence the reduction path, are abandoned and projected to a certain plane; the proximal and distal fracture ends are projected in a rectangular coordinate plane as shown in FIG. 3;
6) planning collision avoidance paths:
soft tissue balance and collision avoidance of tissues around broken bones during reduction operation need to be considered; because the tissues around the fracture are complex, in order to avoid secondary damage in the reduction operation, the collision between the far end and the near end of the femur is avoided during the reduction operation, and then obstacle points are required to be set in the actual situation in the fracture model projection drawing; determining characteristic points such as a reset starting point, a reset ending point, a plurality of obstacle points and the like, wherein the obstacle points are used as obstacle avoidance constraint points for reset path planning; the fracture distal end slowly approaches the proximal end and is in a critical state when collision is about to occur, as shown in fig. 4, the left side is the fracture proximal end, and the right side is the fracture distal end; setting the obstacle range to the left part of the dotted line shown in fig. 4, and the obstacle cannot cross the green dotted line to the left during resetting; fracture reduction in the plane, i.e. moving from the starting point to the end point;
before planning the fracture reduction path by adopting an A-x algorithm, providing a grid diagram of a fracture model, as shown in figure 5; the number and spacing of the grids should follow two conditions: the distance is smaller than the reset precision; the shape of the obstacle can be accurately represented as much as possible under the condition of meeting the precision; in the figure, the symbol "x" represents an obstacle point. The fracture proximal characteristic point is the farthest boundary point; the far-end feature point is, the inner boundary point thereof. The initial reset point is, and the end point is; the symbol @ denotes the planned reset path point, as shown in fig. 6;
7) and during planning of the reduction path, the fracture far ends are respectively projected into two rectangular coordinate planes, and the reduction path planning is respectively carried out by adopting an A-x algorithm in the two coordinate planes. Determining a fracture reduction path by adopting a biplane collision avoidance reduction path planning method based on an A-x algorithm, wherein the plane coordinate value of a reduction path point is a coordinate relative to a characteristic point of the proximal end of the fracture, as shown in FIG. 7; respectively projecting the three-dimensional digital model of the long bone fracture into two rectangular coordinate planes, and respectively planning a reduction path in the two planes by adopting an A-x algorithm;
8) storing coordinate values of a series of path points generated by planning a reset path into a data file, transmitting the data file to a controller of the reset robot, starting the robot to move according to the planned reset path, and testing the reset path in real time by an NDI optical tracking system;
9) fixing NDI optical tracking system navigation frames at the far end and the near end of the fracture respectively, and testing the spatial poses of the near end and the far end of the fracture; and the near end and the far end of the fracture are respectively fixed with an optical tracking array of an NDI positioning system, and the relative spatial postures of the near end and the far end of the fracture are tested.
The embodiment relates to a path planning method for a robot-assisted fracture reduction operation. For long bone fracture types of femur, tibia, etc., after fracture, the distal end is axially, longitudinally, radially or laterally, and rotationally offset relative to the proximal end. Aiming at the closed reduction of the long bone fracture, a robot is adopted to assist the operation of the fracture reduction operation, the near end of the fracture is fixed on an operating bed, and the robot assists the reduction operation of aligning the far end of the fracture to the near end of the fracture. The invention relates to a path planning method for robot-assisted fracture reduction operation, which is used for planning a reduction path based on a C-arm perspective two-dimensional image, namely a fracture distal end posture adjustment and biplane collision prevention reduction path planning method based on an A-x algorithm. The square fracture reduction path planning method of the embodiment of the invention can effectively avoid the large-amplitude traction of the far end in the reduction and the collision of the fracture far end and surrounding tissues, and realize the fracture reduction operation with short time consumption, high precision and safe collision avoidance. The method has the advantages of simplicity, safety, effectiveness and the like, is suitable for robot-assisted fracture reduction operations of long bones such as thighbones and shinbones, can improve the precision and safety of fracture reduction operations, improves reduction efficiency, and reduces the damage of ray radiation to doctors and patients.
In conclusion, the invention provides a path planning method for a robot-assisted fracture reduction operation. Aiming at the fracture types of long bones such as thighbone and shinbone, the fracture near end is fixed on an operating bed, and the robot assists the reduction operation of aligning the fracture far end to the fracture near end. After fracture, longitudinal, transverse and rotational displacement of the distal end relative to the proximal end occurred, with 6 fracture deviation parameters. The invention relates to a path planning method for robot-assisted fracture reduction operation, in particular to a method for adjusting the posture of a fracture far end and planning a biplane collision prevention reduction path based on an A-x algorithm. Determining fracture deviation parameters based on a healthy side mirror image matching method in mimics and Geomagic studio software, and respectively planning fracture reduction paths based on an A-star algorithm in two different coordinate planes. The method for planning the fracture reduction path can effectively avoid the large-amplitude traction of the far end in the reduction and the collision of the fracture far end and surrounding tissues, and realize the fracture reduction operation with short time consumption, high precision and safe collision avoidance. The method has the advantages of simplicity, safety, effectiveness and the like, is suitable for robot-assisted fracture reduction operations of long bones such as thighbones and shinbones, and can improve the precision and safety of fracture reduction operations, improve the reduction efficiency and reduce the working strength of doctors.
The embodiments of the present invention have been described with reference to the accompanying drawings, but the present invention is not limited to the embodiments, and various changes and modifications can be made according to the purpose of the invention, and any changes, modifications, substitutions, combinations or simplifications made according to the spirit and principle of the technical solution of the present invention shall be equivalent replacement means, so long as the technical principle and the inventive concept of the path planning method for robot-assisted fracture reduction surgery of the present invention are met, and the present invention shall fall within the protection scope of the present invention.

Claims (8)

1.一种机器人辅助骨折复位手术的路径规划方法,其特征在于,包括以下步骤:1. a path planning method for robot-assisted fracture reduction surgery, is characterized in that, comprises the following steps: 1)针对长骨骨折闭合复位,采用机器人辅助骨折复位手术操作时,骨折近端固定不动,机器人辅助将骨折远端对准近端进行复位操作;1) For closed reduction of long bone fractures, when a robot-assisted fracture reduction operation is used, the proximal end of the fracture is fixed, and the robot assists the distal end of the fracture to align with the proximal end for reduction operation; 2)术前,根据医生提供的患者术前的骨折侧和健侧骨的CT扫描数据,采用逆向建模方法,进行重建获得骨折侧和健侧骨的三维数字化模型;2) Before the operation, according to the CT scan data of the fractured side and the healthy side of the patient provided by the doctor, the inverse modeling method is used to reconstruct the three-dimensional digital model of the fractured side and the healthy side bone; 3)基于健侧镜像配准的复位原则,在mimics和Geomagic studio软件中,通过图像配准方法,确定骨折轴向、纵向、径向或横向的移位偏差和旋转偏差,并分析骨折移位和旋转偏差值;3) Based on the reduction principle of mirror registration on the unaffected side, in mimics and Geomagic studio software, through the image registration method, determine the axial, longitudinal, radial or lateral displacement deviation and rotation deviation of the fracture, and analyze the fracture displacement and rotation deviation value; 4)骨折远端空间姿态调整:4) Spatial posture adjustment of the distal end of the fracture: 由机器人调整骨折远端相对于近端姿态,使其相对于骨折近端达到合理的解剖姿态,同时由固定在骨折远端的NDI光学跟踪阵列实时测试骨折远端的3个姿态角,保证骨折远端的姿态调整以及完成复位后满足复位精度要求;The robot adjusts the posture of the distal end of the fracture relative to the proximal end to achieve a reasonable anatomical posture relative to the proximal end of the fracture. At the same time, the NDI optical tracking array fixed at the distal end of the fracture tests the three posture angles of the distal end of the fracture in real time to ensure the fracture. The attitude adjustment of the remote end and the reset accuracy requirements are met after the reset is completed; 5)复位路径规划准备:5) Reset path planning preparation: 首先,给出骨折模型的栅格图,栅格的数量和间距遵循两个条件:First, a grid diagram of the fracture model is given. The number and spacing of grids obey two conditions: (1)栅格间距小于复位精度;(1) The grid spacing is smaller than the reset accuracy; (2)满足精度前提下准确的表示障碍物的形状;(2) Accurately represent the shape of the obstacle under the premise of satisfying the accuracy; 然后,将骨折远端和骨折近端各自的下半部分近似为圆柱形,舍弃近端和远端不影响复位路径的部分;Then, the lower half of the distal end and the proximal end of the fracture are approximately cylindrical, and the parts of the proximal and distal ends that do not affect the reduction path are discarded; 6)避碰路径规划:6) Collision avoidance path planning: 利用基于A*算法的骨折复位路径规划算法,确定最佳避碰复位路径;Use the fracture reduction path planning algorithm based on the A* algorithm to determine the optimal collision avoidance reduction path; 7)复位路径规划:7) Reset path planning: 将长骨骨折三维数字化模型,分别投影到两个直角坐标平面内,在两个平面内分别采用A*算法进行复位路径规划;Project the 3D digital model of long bone fractures into two rectangular coordinate planes respectively, and use the A* algorithm to plan the reduction path in the two planes respectively; 8)将基于A*算法规划的复位路径规划所获得的数据点的坐标保存为数据文件,并传输给复位机器人控制器,机器人控制骨折远端按照所规划的复位路径进行复位手术操作;8) Save the coordinates of the data points obtained by the reduction path planning based on the A* algorithm as a data file and transmit it to the reduction robot controller, and the robot controls the distal end of the fracture to perform reduction operations according to the planned reduction path; 9)在骨折远端和近端分别固定NDI光学跟踪系统导航架,并测试骨折近端和远端的空间位姿;骨折近端和远端分别固定一个NDI定位系统的光学跟踪阵列,测试骨折近端和远端相对空间姿态。9) Fix the NDI optical tracking system navigation frame at the distal and proximal ends of the fracture, and test the spatial pose of the proximal and distal ends of the fracture; fix an optical tracking array of the NDI positioning system at the proximal and distal ends of the fracture, respectively, to test the fracture Relative spatial poses of the proximal and distal ends. 2.根据权利要求书1所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤6)中,采用基于A*算法的双平面避碰复位路径规划方法,确定骨折复位路径;将骨折近端和远端投影到两个直角坐标平面内,在两个坐标平面内分别采用A*算法进行复位路径规划。2. The path planning method for robot-assisted fracture reduction surgery according to claim 1, characterized in that: in step 6), a biplane collision avoidance reduction path planning method based on A* algorithm is used to determine fracture reduction Path; project the proximal and distal ends of the fracture into two rectangular coordinate planes, and use the A* algorithm to plan the reduction path in the two coordinate planes respectively. 3.根据权利要求书1所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤6)中,A*算法通过估计函数来确定搜索方向,从起点开始向周围扩展,通过估价函数计算得到周围每个节点的代价值,选择最小代价节点作为下一个扩展节点,重复这一过程直到达到目标点,生成最终复位路径。3. The path planning method for robot-assisted fracture reduction surgery according to claim 1, characterized in that: in step 6), the A* algorithm determines the search direction by estimating the function, and expands from the starting point to the surrounding area, The cost value of each surrounding node is calculated by the evaluation function, the minimum cost node is selected as the next expansion node, and this process is repeated until the target point is reached, and the final reset path is generated. 4.根据权利要求书3所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤6)中,在搜索过程中,由于路径上的每个节点都是具有最小代价的节点,因此得到的路径代价是最小的;A*算法的估价函数为:;其中,表示从起始点经过任意节点到达目标点的估价函数,表示起始点到节点的实际代价,表示节点n到目标点的估计代价。4. The path planning method for robot-assisted fracture reduction surgery according to claim 3, characterized in that: in step 6), in the search process, since each node on the path has the minimum cost node, so the obtained path cost is the smallest; the evaluation function of the A* algorithm is: ; where, represents the evaluation function from the starting point to the target point through any node, represents the actual cost from the starting point to the node, and represents the node n to the goal. Estimated cost of points. 5.根据权利要求书3所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤6)中,计算估价函数时,使用了曼哈顿距离和欧氏距离的计算公式;当节点沿栅格图的边长方向搜索时,采用曼哈顿距离公式计算;当节点沿栅格图的对角线方向搜索时,采用欧几里得距离公式计算。5. The path planning method for robot-assisted fracture reduction surgery according to claim 3, characterized in that: in step 6), when calculating the evaluation function, the calculation formulas of Manhattan distance and Euclidean distance are used; when When a node is searched along the side length of the grid graph, it is calculated by the Manhattan distance formula; when the node is searched along the diagonal direction of the grid graph, it is calculated by the Euclidean distance formula. 6.根据权利要求书1所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤1)中,针对股骨或胫骨的长骨骨折,采用机器人辅助进行闭合式复位手术操作时,将骨折近端固定于手术床上,机器人辅助将骨折远端对准骨折近端的复位操作。6. The path planning method for robot-assisted fracture reduction surgery according to claim 1, wherein in step 1), for long bone fractures of the femur or tibia, when a robot-assisted closed reduction surgery is used , fix the proximal end of the fracture on the operating table, and the robot assists the reduction operation of aligning the distal end of the fracture with the proximal end of the fracture. 7.根据权利要求书1所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤4)中,导航架固定于骨折远端,机器人根据骨折旋转偏差值进行骨折远端的姿态调整,并实时观察导航架与NDI基准坐标系间的姿态角,直至骨折远端姿态调整到位。7. The path planning method for robot-assisted fracture reduction surgery according to claim 1, characterized in that: in step 4), the navigation frame is fixed on the distal end of the fracture, and the robot performs the distal end of the fracture according to the rotation deviation value of the fracture. Attitude adjustment, and observe the attitude angle between the navigation frame and the NDI reference coordinate system in real time until the distal end of the fracture is adjusted in place. 8.根据权利要求书1所述机器人辅助骨折复位手术的路径规划方法,其特征在于:在所述在步骤5)中,投影到直角坐标平面内的骨折模型投影图中,确定进行骨折复位路径规划的起始点、终止点以及路径上的若干障碍点一系列特征点,其中障碍点作为复位路径规划的避障约束点,以保证复位操作时应避免股骨远端与近端发生碰撞。8. The path planning method for robot-assisted fracture reduction surgery according to claim 1, characterized in that: in step 5), the fracture reduction path is determined by projecting it onto the fracture model projection diagram in the rectangular coordinate plane The planned start point, end point, and a series of characteristic points of several obstacle points on the path, among which the obstacle point is used as the obstacle avoidance constraint point of the reduction path planning to ensure that the collision between the distal end and the proximal end of the femur should be avoided during the reduction operation.
CN202010061853.1A 2020-01-20 2020-01-20 Robot-assisted fracture reduction path planning method Pending CN111134842A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010061853.1A CN111134842A (en) 2020-01-20 2020-01-20 Robot-assisted fracture reduction path planning method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010061853.1A CN111134842A (en) 2020-01-20 2020-01-20 Robot-assisted fracture reduction path planning method

Publications (1)

Publication Number Publication Date
CN111134842A true CN111134842A (en) 2020-05-12

Family

ID=70526228

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010061853.1A Pending CN111134842A (en) 2020-01-20 2020-01-20 Robot-assisted fracture reduction path planning method

Country Status (1)

Country Link
CN (1) CN111134842A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111590584A (en) * 2020-05-27 2020-08-28 京东方科技集团股份有限公司 Determination method and device for safety limit area, reset method and medical robot
CN112043384A (en) * 2020-07-29 2020-12-08 上海大学 An external force prediction method for fracture reduction robot
CN112370154A (en) * 2020-11-13 2021-02-19 山东中医药大学附属医院 Fracture distal end reduction trajectory planning method and reduction system based on mirror image method
CN112741688A (en) * 2020-12-11 2021-05-04 北京信息科技大学 Path planning method for minimally invasive surgery software manipulator
CN112932662A (en) * 2021-02-01 2021-06-11 威海威高骨科手术机器人有限公司 Reset position detection method with patient safety protection function
CN114451991A (en) * 2021-01-30 2022-05-10 威海威高骨科手术机器人有限公司 Fracture reduction path planning method based on virtual reduction collision detection
CN114974509A (en) * 2022-05-26 2022-08-30 哈尔滨工业大学 Fracture reduction path planning method, fracture reduction method and electronic equipment
WO2023000560A1 (en) * 2021-07-21 2023-01-26 天津大学 Reduction trajectory automatic planning method for parallel fracture surgical robot
CN117237242A (en) * 2023-11-16 2023-12-15 深圳爱递医药科技有限公司 Oral and maxillofacial surgery postoperative care system based on structured light data
CN118787444A (en) * 2024-06-14 2024-10-18 中国科学院深圳先进技术研究院 A robot-assisted fracture reduction path planning method and system
CN119791842A (en) * 2025-03-12 2025-04-11 锦恒科技(大连)有限公司 A method and system for reducing long bone fractures based on image processing

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021238876A1 (en) * 2020-05-27 2021-12-02 京东方科技集团股份有限公司 Safety limit region determination method and device, resetting method, and medical robot
US12172324B2 (en) 2020-05-27 2024-12-24 Boe Technology Group Co., Ltd. Method for determining safety-limit zone, and device, reset method and medical robot using the same
CN111590584A (en) * 2020-05-27 2020-08-28 京东方科技集团股份有限公司 Determination method and device for safety limit area, reset method and medical robot
CN111590584B (en) * 2020-05-27 2021-12-10 京东方科技集团股份有限公司 Determination method and device for safety limit area, reset method and medical robot
CN112043384B (en) * 2020-07-29 2023-07-18 上海大学 A fracture reduction robot external force prediction system
CN112043384A (en) * 2020-07-29 2020-12-08 上海大学 An external force prediction method for fracture reduction robot
CN112370154A (en) * 2020-11-13 2021-02-19 山东中医药大学附属医院 Fracture distal end reduction trajectory planning method and reduction system based on mirror image method
CN112741688A (en) * 2020-12-11 2021-05-04 北京信息科技大学 Path planning method for minimally invasive surgery software manipulator
CN112741688B (en) * 2020-12-11 2023-09-29 北京信息科技大学 Path planning method for minimally invasive surgery software manipulator
CN114451991B (en) * 2021-01-30 2024-02-27 威海威高骨科手术机器人有限公司 Fracture reduction path planning method based on virtual reduction collision detection
CN114451991A (en) * 2021-01-30 2022-05-10 威海威高骨科手术机器人有限公司 Fracture reduction path planning method based on virtual reduction collision detection
CN112932662A (en) * 2021-02-01 2021-06-11 威海威高骨科手术机器人有限公司 Reset position detection method with patient safety protection function
WO2023000560A1 (en) * 2021-07-21 2023-01-26 天津大学 Reduction trajectory automatic planning method for parallel fracture surgical robot
CN114974509A (en) * 2022-05-26 2022-08-30 哈尔滨工业大学 Fracture reduction path planning method, fracture reduction method and electronic equipment
CN114974509B (en) * 2022-05-26 2024-05-10 哈尔滨工业大学 Fracture reduction path planning method, fracture reduction method and electronic equipment
CN117237242A (en) * 2023-11-16 2023-12-15 深圳爱递医药科技有限公司 Oral and maxillofacial surgery postoperative care system based on structured light data
CN117237242B (en) * 2023-11-16 2024-02-27 深圳爱递医药科技有限公司 Oral maxillofacial surgery postoperative care system based on structured light data
CN118787444A (en) * 2024-06-14 2024-10-18 中国科学院深圳先进技术研究院 A robot-assisted fracture reduction path planning method and system
CN118787444B (en) * 2024-06-14 2025-08-15 中国科学院深圳先进技术研究院 Robot-assisted fracture reduction path planning method and system
CN119791842A (en) * 2025-03-12 2025-04-11 锦恒科技(大连)有限公司 A method and system for reducing long bone fractures based on image processing

Similar Documents

Publication Publication Date Title
CN111134842A (en) Robot-assisted fracture reduction path planning method
US12059215B2 (en) Navigation system and method for joint replacement surgery
CN112641510B (en) Joint replacement surgery robot navigation and positioning system and method
JP4815054B2 (en) System and method for performing an image-directed robotic orthopedic procedure without utilizing a trust criteria system
CN109984843B (en) Fracture closed reduction navigation system and method
CN111887989A (en) A robot-assisted fracture space collision avoidance reduction path planning method
Wu et al. Motion planning of continuum tubular robots based on centerlines extracted from statistical atlas
CN109124835A (en) The localization method and system of late-segmental collapse point
Burghart et al. Robot assisted craniofacial surgery: first clinical evaluation
CN116784975B (en) Flexible puncture needle path planning method based on improved Bi-RRT algorithm
CN218784458U (en) End effector and surgical system
CN112370158B (en) Robot fracture reduction path planning method and system and robot
CN112370156B (en) Memory alloy reset external fixing system based on virtual reset and use method
CN119112369A (en) Surgical treatment methods, devices, surgical robot systems and related equipment
CN116725680B (en) End effector and surgical system
CN219579025U (en) Full-featured orthopedic surgery control system
CN218784457U (en) End effector and surgical system
CN115153836B (en) A 3D osteotomy preoperative intelligent planning method and system
US20220233256A1 (en) Assistive surgical robot for distal hole localization in intramedullary nail
Zhu et al. Development of a surgeon and patient-friendly orthopedic surgical robot
CN115444523A (en) An automatic puncture control method and system
CN115300041A (en) Acetabular osteotomy orthopedic surgery positioning tool, system and computer equipment
Liu et al. A novel computer-assisted system for long bone fracture reduction with a hexapod external fixator
CN119587155B (en) Remotely controlled pelvis fracture reduction system
JP2022535772A (en) lockable surgical system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200512