WO2022141153A1 - Ultrasonic positioning puncture system and storage medium - Google Patents

Ultrasonic positioning puncture system and storage medium Download PDF

Info

Publication number
WO2022141153A1
WO2022141153A1 PCT/CN2020/141231 CN2020141231W WO2022141153A1 WO 2022141153 A1 WO2022141153 A1 WO 2022141153A1 CN 2020141231 W CN2020141231 W CN 2020141231W WO 2022141153 A1 WO2022141153 A1 WO 2022141153A1
Authority
WO
WIPO (PCT)
Prior art keywords
coordinate system
joint
ultrasonic detection
ultrasonic
detection device
Prior art date
Application number
PCT/CN2020/141231
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 诺创智能医疗科技(杭州)有限公司
Priority to PCT/CN2020/141231 priority Critical patent/WO2022141153A1/en
Publication of WO2022141153A1 publication Critical patent/WO2022141153A1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34Trocars; Puncturing needles
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves

Definitions

  • the present application relates to the field of robotic arm control, and in particular, to an ultrasonic positioning puncture system and a storage medium.
  • the doctor manually sends small surgical instruments (such as biopsy needles) into the patient's body to detect or treat the lesion.
  • small surgical instruments such as biopsy needles
  • the nearby 2D or 3D scan images determine the appropriate needle entry point and needle entry direction, and then manually complete the puncture operation based on experience.
  • the working principle of the puncture robot with the surgical navigation system is as follows: first, the computer 3D synthesizes the scanned 2D image to form a 3D image near the lesion, and then the doctor judges the target position and the appropriate needle insertion path through the 3D image and inputs it into the navigation system. , and then the navigation system calculates the current state and target point of the robot manipulator arm and plans a trajectory. Finally, the manipulator arm completes the puncture positioning according to the planned trajectory, and then manually completes the needle insertion through the needle insertion mechanism or the doctor, avoiding pure manual operation. errors caused.
  • the existing navigation technology usually locates and registers the puncture needle and medical imaging equipment through optical positioning, and performs puncture with the aid of medical images, but there are large errors in the process of optical positioning and registration, and in During the operation, the optical ball cannot be blocked, which limits the implementation of the operation.
  • an ultrasonic positioning and puncturing system includes a computer device, a first serial manipulator arm, and an ultrasonic detection device mounted on the end of the first serial manipulator arm.
  • Computer equipment is electrically connected to the first serial robotic arm and the ultrasonic detection device respectively, the computer equipment includes a memory and a processor, the memory stores a computer program, and the processor runs the computer program Perform the following steps: locate the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device; obtain the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial robot arm ; According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located; According to the position conversion relationship and the described position conversion relationship The first position information determines the second position information of the target point in the reference coordinate system.
  • the above ultrasonic positioning puncture system has the following advantages: by directly obtaining the coordinates of the lesion in the ultrasonic image, without the need for indirect registration through registration datums such as optical positioning, it solves the existing positioning and registration technology for puncture surgery under medical image navigation There is also the problem of error, which improves the accuracy of the puncture operation.
  • the locating the first position information of the target point in the ultrasound detection plane projected by the ultrasound detection device includes: acquiring an ultrasound image acquired by the ultrasound detection device, in which the ultrasound image is Locate the position information of the target point; according to the position information of the target point in the ultrasonic image and the calibration parameters of the ultrasonic detection device, determine that the target point is in the ultrasonic detection plane projected by the ultrasonic detection device of the first location information.
  • the ultrasonic detection device is rotatably connected to the end of the first serial manipulator, and the ultrasonic detection plane is formed at the end of the ultrasonic detection device and extends from the ultrasonic detection device.
  • the end point of the ultrasonic detection plane extends along the fan plane; the ultrasonic detection plane rotates with the ultrasonic detection device, and the rotation axis of the ultrasonic detection plane passes through the end point; the first attitude information includes the rotation angle of the ultrasonic detection plane.
  • the first position information includes the distance of the target point relative to the end point, and the included angle of the connecting line between the target point and the end point and the rotation axis.
  • the first serial robotic arm includes a plurality of joints
  • the plurality of joints include: rotation joints and/or movement joints
  • the second posture information includes motion amount information of each of the joints .
  • the determining, according to the first attitude information and the second attitude information, the position conversion relationship between the reference coordinate system of the ultrasound positioning and puncturing system and the coordinate system where the ultrasound detection plane is located includes: : when the first serial manipulator is in the initial state, establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the In the initial state, a probe coordinate system is established with the end point of the ultrasonic detection device as the coordinate origin; the reference coordinate system of the first serial manipulator and the DH of the joint coordinate system of each joint of the first serial manipulator are obtained.
  • the parameters determine the first position conversion relationship between the reference coordinate system of the first serial manipulator and the joint coordinate system of the end joint of the first serial manipulator; obtain the probe coordinates The position information of the coordinate origin of the system in the joint coordinate system of the terminal joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the terminal joint, and according to the position information, attitude angle and the first a position conversion relationship, determining a second position conversion relationship between the reference coordinate system and the probe coordinate system; acquiring the rotation angle of the ultrasonic detection plane, and determining the coordinates of the ultrasonic detection plane according to the rotation angle The third position conversion relationship between the system and the probe coordinate system; according to the second position conversion relationship and the third position conversion relationship, determine the reference coordinate system of the ultrasonic positioning and puncturing system and the ultrasonic detection plane.
  • the first serial manipulator includes a rotating joint and a moving joint
  • the Z axis of the joint coordinate system of the rotating joint is set along the rotating axis
  • the Z axis of the joint coordinate system of the moving joint is set along the moving direction
  • the reference coordinate system and the joint coordinate system of each joint are both left-handed or right-handed systems
  • the origin of the joint coordinate system of the rotary joint is the joint coordinate system of the previous mobile joint. The origin coincides.
  • the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the joint coordinate system of the distal joint, and the Z axis of the probe coordinate system The axis is collinear with the Z-axis of the joint coordinate system of the end joint.
  • the end joint of the first serial manipulator is a parallel motion platform
  • the parallel motion platform includes a static platform, a moving platform, and a plurality of pieces disposed between the static platform and the moving platform
  • the telescopic assembly, the ultrasonic detection device is set on the moving platform.
  • the ultrasonic positioning and puncturing system further includes a second serial robotic arm and a puncturing instrument mounted on the end of the second serial robotic arm; the computer equipment is further connected to the second serial robotic arm, respectively. electrically connected, the processor further executes the following steps when running the computer program: taking the second position information of the target point in the reference coordinate system as the target puncture point, and controlling the second serial robotic arm The puncture instrument punctures the target puncture point.
  • a storage medium in which a computer program is stored, the computer program being configured to execute the following steps when running: in an ultrasound detection plane projected by the ultrasound detection device Locating the first position information of the target point; obtaining the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator; according to the first attitude information and the second attitude information, determine The position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located; according to the position conversion relationship and the first position information, it is determined that the target point is in the reference coordinate system the second location information.
  • FIG. 1 is a structural block diagram of an ultrasonic positioning and puncturing system provided by an embodiment of the present application.
  • FIG. 2 is a flowchart of a working method of the ultrasonic positioning puncture system provided by the embodiment of the present application.
  • FIG. 3 is a flowchart of a coordinate conversion method for a multi-joint robotic arm according to an embodiment of the present application.
  • FIG. 4 is a schematic structural diagram of a multi-joint robotic arm provided by a preferred embodiment of the present application.
  • FIG. 5 is a schematic diagram of a mechanical coordinate system and a joint coordinate system of a multi-joint manipulator provided by a preferred embodiment of the present application.
  • FIG. 6 is a schematic diagram of a probe coordinate system and a Stewart calculation coordinate system of an ultrasonic detection device according to a preferred embodiment of the present application.
  • the ultrasonic positioning puncture system and storage medium provided by the present application will be further described below.
  • FIG. 1 is a structural block diagram of an ultrasonic positioning and puncturing system provided by an embodiment of the present application.
  • the ultrasonic positioning and puncturing system includes a computer device 101 , a first serial robotic arm 102 , and a device mounted on the end of the first serial robotic arm 102 .
  • the ultrasonic detection device 103 and the computer equipment 101 are electrically connected to the first serial manipulator 102 and the ultrasonic detection device 103 respectively.
  • the computer equipment 101 includes a memory 1011 and a processor 1012, and the memory 1011 stores a computer program 10111.
  • FIG. 2 is a flowchart of the working method of the ultrasonic positioning and puncturing system provided by the embodiment of the present application. As shown in FIG. 2 , the processor 1012 executes the following steps S201 to S204 when running a computer program.
  • Step S201 Locating the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device.
  • the ultrasonic detection plane projected by the ultrasonic detection device 103 is usually in the shape of a fan, and the fan center of the fan is the starting point of the ultrasonic generation.
  • the ultrasonic detection is sent at the fan center and the echoes reflected by the substance are received to detect the target position in the ultrasonic detection plane.
  • the first serial machine equipped with the ultrasonic detection device 103 The arm 102 needs to move the ultrasonic detection device 103 to a position near the target point.
  • the movement amount of each joint of the first serial robotic arm 102 will change. It is called the second attitude information.
  • the first serial robotic arm includes a plurality of joints, the plurality of joints include: a rotation joint and/or a movement joint, and the second posture information includes motion amount information of each joint.
  • the ultrasonic detection device is rotatably connected to the end of the first serial robotic arm 102, and the ultrasonic detection plane is formed at the end of the ultrasonic detection device and extends from the end of the ultrasonic detection device along the sector; ultrasonic detection The plane rotates following the ultrasonic detection device, and the rotation axis of the ultrasonic detection plane passes through the end point; the first attitude information includes the rotation angle of the ultrasonic detection plane.
  • the ultrasonic detection plane is irradiated toward the target point.
  • the ultrasonic detection device 103 is rotated along the symmetry axis of the fan plane. , the ultrasonic detection plane rotates accordingly, so as to realize the target detection in the truncated or conical space area, and the ultrasonic image of each ultrasonic detection plane will be displayed on the electronic screen in real time.
  • the position of the target point in the real ultrasound detection plane can be determined from the position of the target point in the ultrasound image, which is the first position information.
  • the calibration parameters of the device determine the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device.
  • the calibration parameters of the ultrasonic detection device include, but are not limited to, at least one of the following: the number of array elements, the radius of curvature, the width of the linear array, the field of view (FOV), and the like.
  • the first position information may be represented in the form of Cartesian coordinates or polar coordinates.
  • the first position information includes the distance of the target point relative to the end point, and the angle between the connecting line between the target point and the end point and the rotation axis.
  • Step S202 Acquire first attitude information of the ultrasonic detection plane and second attitude information of the first serial manipulator.
  • Step S203 According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located.
  • Step S204 Determine the second position information of the target point in the reference coordinate system according to the position conversion relationship and the first position information.
  • the conversion relationship between any point on the ultrasonic detection plane and the reference coordinate system is determined. Therefore, according to the first attitude information and the second attitude information, can determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located, and then according to the position conversion relationship and the first position information, it can be determined that the target point is in the reference coordinate system.
  • the second position information of realizes the positioning of the target point in the reference coordinate system.
  • step S201 may also be performed after step S202 or step S203.
  • determining the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located includes the following steps.
  • Step S203-1 When the first serial manipulator is in the initial state, establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the initial state, The probe coordinate system is established with the end point of the ultrasonic detection device as the coordinate origin.
  • Step S203-2 Obtain the reference coordinate system of the first serial manipulator and the DH parameters of the joint coordinate system of each joint of the first serial manipulator, and determine the reference coordinate system of the first serial manipulator to the first serial manipulator according to the DH parameter.
  • Step S203-3 Obtain the position information of the coordinate origin of the probe coordinate system in the joint coordinate system of the terminal joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the terminal joint, and according to the position information, attitude angle and the first A position conversion relationship, which determines a second position conversion relationship between the reference coordinate system and the probe coordinate system.
  • Step S203-4 Obtain the rotation angle of the ultrasonic detection plane, and determine a third position conversion relationship between the coordinate system where the ultrasonic detection plane is located and the probe coordinate system according to the rotation angle.
  • Step S203-5 According to the second position conversion relationship and the third position conversion relationship, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning and puncturing system and the coordinate system where the ultrasonic detection plane is located.
  • step S203-2 aims to determine the first position conversion relationship between the reference coordinate system of the first serial manipulator and the joint coordinate system of the end joint of the first serial manipulator.
  • this embodiment provides a coordinate conversion method for a multi-joint manipulator, which is applied to a manipulator including a rotating joint and a moving joint, and these rotating joints and moving joints are sequentially connected in series .
  • the number of rotating joints and moving joints can be set according to the required degrees of freedom and motion range, that is, the number of rotating joints and moving joints can be one or more.
  • FIG. 3 is a flowchart of a coordinate conversion method for a multi-joint robotic arm according to an embodiment of the present application. As shown in FIG. 3 , the flowchart includes the following steps S301 to S303 .
  • Step S301 When the robot arm is in an initial state, a mechanical coordinate system of the robot arm and a joint coordinate system of each joint of the robot arm are established.
  • the control of the manipulator is usually based on DH parameters for coordinate system transformation, and DH parameters include traditional DH parameters and improved DH parameters.
  • Two connected joints are adjacent joints.
  • the transformation of the joint coordinate system of two adjacent joints is usually represented by a DH parameter or an improved DH parameter.
  • a DH parameter as an example, two adjacent joint coordinate systems can overlap each other by rotating ⁇ around the Z axis and translating d, and then rotating ⁇ around the X axis and translating a.
  • ⁇ , d, ⁇ and a are DH parameters. It can be seen from this that the simpler the DH parameter is, the simpler the transformation of two adjacent joint coordinate systems will be.
  • the Z axis of the joint coordinate system of the rotating joint is set along the rotation axis
  • the Z axis of the joint coordinate system of the moving joint is set along the moving direction
  • the mechanical coordinate system is the same as the joint coordinate system of each joint. It is a left-handed system or a right-handed system, and when the joint before the rotating joint is a moving joint, the origin of the joint coordinate system of the rotating joint coincides with the origin of the joint coordinate system of the moving joint.
  • the above-mentioned mechanical coordinate system can also be called a reference coordinate system or a world coordinate system.
  • the mechanical coordinate system is usually set at the center of the base of the manipulator, and is located directly below the first joint of the multi-joint arm, which can make the mechanical coordinate system
  • the transformation relationship to the first joint is as simple as possible.
  • a rotary joint In a serial manipulator, in most cases, when a rotary joint receives a command with a rotation angle of 0 or 2 ⁇ , it may not be necessary to distinguish between the two, but it will not rotate according to a preset setting or rotate by 2 ⁇ in a set direction. However, in some cases, it is necessary to distinguish between these two corners. In the case where the 0 or 2 ⁇ rotation angle needs to be distinguished, the rotation angle of the Z axis in the DH parameter of the rotary joint is not 0 or 2 ⁇ , so as to avoid confusing the rotation angle of 0 or 2 ⁇ .
  • Step S302 Obtain the mechanical coordinate system of the manipulator and the DH parameters of the joint coordinate system of each joint of the manipulator, and determine the conversion between the mechanical coordinate system of the manipulator and the joint coordinate system of the end joint of the manipulator according to the DH parameters relation.
  • the DH parameters between adjacent joints can be obtained in turn, and according to the DH parameters, the conversion relationship between the mechanical coordinate system of the manipulator and the joint coordinate system of the end joint of the manipulator can be determined. .
  • the homogeneous transformation from the coordinates of the i-1th joint to the coordinates of the ith joint is constructed as a sequence with two rotations and two transformations, which can be expressed as follows using a matrix:
  • n is the total number of rotating joints and moving joints of the robotic arm.
  • the DH parameter of the first joint of the manipulator is the transformation of the coordinate system between the first joint and the mechanical coordinate system, which is recorded as Then the conversion relationship between the mechanical coordinate system and the joint coordinate system of the first joint:
  • the Nth joint is the end joint.
  • Step S303 Perform coordinate conversion between the mechanical coordinate system and the joint coordinate system of the end joint according to the conversion relationship.
  • step S302 It represents the coordinate conversion relationship between the mechanical coordinate system and the joint coordinate system of the end joint. According to the conversion relationship, the coordinate conversion between the mechanical coordinate system and the joint coordinate system of the end joint can be realized.
  • FIG. 4 is a schematic structural diagram of a multi-joint mechanical arm provided by a preferred embodiment of the present application.
  • the multi-joint mechanical arm shown in FIG. The moving joint 6 , the rotating joint 7 , the moving joint 8 , the rotating joint 9 and the moving joint 10 , in addition, the multi-joint robotic arm further includes a base 11 fixedly connected with the moving joint 1 .
  • the coordinate transformation method for a multi-joint robotic arm in this preferred embodiment includes steps 1 to 3, respectively corresponding to steps S301 to S303 in FIG. 3 .
  • Step 1 According to the rules of the world coordinate system, the mechanical coordinate system and the joint coordinate system of each joint are established at the base of the manipulator.
  • FIG. 5 is a schematic diagram of a mechanical coordinate system and a joint coordinate system of a multi-joint manipulator provided by a preferred embodiment of the present application.
  • the origin of the mechanical coordinate system F 0 -X 0 Y 0 Z 0 is fixed at the base of the robot arm, the Z 0 axis points from F 0 to the moving joint 1 , and the Y 0 axis points from the base F 0 Point to the robotic arm, and the X 0 axis points to conform to the right-handed coordinate system.
  • the origin L 1 of the joint coordinate system L 1 -X 1 Y 1 Z 1 of the mobile joint 1 is fixed on the mobile joint 1 , and the orientation of each coordinate axis is the same as that of the corresponding axis of the mechanical coordinate system.
  • the origin of the joint coordinate system R 2 -X 2 Y 2 Z 2 of the rotary joint 2 is fixed on the rotary joint 2 and coincides with L 1.
  • the Z 2 and Z 1 axes point in the same direction, the X 2 axis and the Y 2 axis The directions are opposite to the X 1 axis and Y 1 axis respectively.
  • the joint coordinate system L 3 -X 3 Y 3 Z 3 of the mobile joint 3 is fixed at the origin L 3 of the mobile joint 3, the Z 3 axis points from the L 1 point to the L 3 point, and the X 3 axis and the Y 3 axis are respectively connected with the X 2 Axis and Z 2 axis point the same.
  • the joint coordinate system R 4 -X 4 Y 4 Z 4 of the rotary joint 4 is fixed at the origin R 4 of the rotary joint 4 and coincides with the point L 3 (Fig. 5 is a clear indication of the joint coordinate system of the mobile joint 3 and the rotary joint 4).
  • the joint coordinate system of L 3 and R 4 are marked separately, the same below), the Z 4 axis points to the opposite direction to the Y 3 axis, initially, the X 4 axis and the Y 4 axis point to the X 3 axis and the Z 3 axis respectively. on the contrary.
  • the joint coordinate system R 5 -X 5 Y 5 Z 5 of the rotary joint 5 is fixed at the origin R 5 of the rotary joint 5.
  • the direction of the Z 5 axis and the Z 4 axis are the same. Initially, the X 5 axis and the Y 5 axis point respectively.
  • the X 4 -axis and Y 4 -axis point opposite.
  • the joint coordinate system L 6 -X 6 Y 6 Z 6 of the mobile joint 6 is fixed at the origin L 6 of the mobile joint 6.
  • the Z 6 axis points from the R 5 point to the L 6 point, and the X 6 axis and the Y 6 axis point to the X axis respectively.
  • the 5 -axis and the Z 5 -axis point the same.
  • the joint coordinate system R 7 -X 7 Y 7 Z 7 of the rotary joint 7 is fixed at the origin R 7 of the rotary joint 7 and coincides with the L 6 point. Initially, the Z 7 axis and the Y 7 axis point to the Y 6 axis respectively. The X 7 axis and the Z 6 axis point in the same direction, as opposed to the X 6 axis.
  • the origin of the joint coordinate system L 8 -X 8 Y 8 Z 8 of the mobile joint 8 is fixed on the mobile joint 8.
  • the Z 8 axis points from the L 8 point to the R 7 point, and the directions of the X 8 and Y 8 axes are respectively the same as The X 7 and Y 7 axes point the same.
  • the origin of the joint coordinate system R9 - X9Y9Z9 of the rotary joint 9 is fixed on the rotary joint 9 and coincides with the L8 point.
  • the direction of the Z9 axis is opposite to the direction of the Y8 axis.
  • the direction of the X 9 axis is opposite to the direction of the Z 8 axis, and the direction of the Y 9 axis is the same as the direction of the X8 axis.
  • the origin of the joint coordinate system L 10 -X 10 Y 10 Z 10 of the mobile joint 10 is fixed on the mobile joint 10.
  • the Z 10 axis points from the R 9 point to the L 10 point, and the X 10 axis and the Y 10 axis point respectively. The same point as the X 9 and Z 9 axes.
  • the length of F 0 R 2 is l 1
  • the length of R 2 R 4 is l 2
  • the length of R 4 R 5 is l 3
  • the length of R 5 R 7 is l 4
  • the length of R 7 R 9 is l 5
  • the length of R 9 L 10 is l 6
  • the points L 1 , R 2 , L 3 , R 4 , R 5 , L 6 and R 7 are located on the same horizontal plane.
  • Step 2 Obtain the DH parameters and calculate the conversion relationship from the mechanical coordinate system to the end joints of the multi-joint robotic arm.
  • DH parameters are shown in Table 1.
  • Table 1 in order to avoid confusing the rotation angles 0 and 2 ⁇ of the rotary joint, a deflection angle is added to the Z-axis rotation angle ⁇ of each rotary joint.
  • the Z-axis rotation angle of the rotary joint 2 is ⁇ + ⁇ 2 , where ⁇ 2 is the deflection angle.
  • a parallel platform can also be mounted on the end joint of the multi-joint manipulator shown in Figure 4.
  • the parallel platform includes a static platform, a moving platform, and a plurality of telescopic components arranged between the static platform and the moving platform.
  • the calculation coordinate system of the parallel platform It is completely coincident with the joint coordinate system of the end joint of the multi-joint mechanical arm, and the ultrasonic detection device is arranged on the moving platform.
  • the parallel platform in this embodiment can realize motion with multiple degrees of freedom.
  • the Stewart parallel platform includes a static platform, a moving platform, and a plurality of telescopic and retractable platforms arranged between the static platform and the moving platform.
  • the component can realize six degrees of freedom movement in space, namely displacement along the X axis, displacement along the Y axis, displacement along the Z axis, rotation around the X axis, rotation around the Y axis, and rotation around the Z axis.
  • the Stewart parallel platform is supported by 6 telescopic components. Compared with the passive arm using the cantilever beam structure in series, the rigidity is higher and the structure is stable.
  • the parallel structure Due to the high rigidity, the parallel structure has a higher self-weight or volume than the series structure. Carrying capacity.
  • the error at the end of the passive arm using the cantilever beam structure in series is the accumulation and amplification of the errors of each joint, so the error is large and the precision is low. Precision surgical operation.
  • the inverse solution of the parallel platform is very easy, and it is easy to obtain the motion posture of each telescopic component of the parallel platform according to the coordinate position.
  • the displacement of the master operator can be mapped to the slave user coordinate system by establishing the slave user coordinate system on the static platform, and then the slave user coordinate system can be used to map the displacement of the master operator into the slave user coordinate system.
  • the position information of the target position in the calculation coordinate system can be obtained by the conversion of the system and the calculation coordinate system. According to the position information, the motion posture of each telescopic assembly of the parallel platform can be easily obtained through the inverse solution.
  • the transformation matrix from the mth joint (containing the m degree of freedom) to the nth joint can be expressed as:
  • T trans_mach_st the transformation matrix from the Stewart calculation coordinate system to the mechanical coordinate system
  • Step 3 According to the transformation matrix and the inverse matrix of the transformation matrix, coordinate transformation between the mechanical coordinate system and the Stewart calculation coordinate system is performed.
  • the directions of the X, Y, and Z axes of the probe coordinate system and the directions of the X, Y, and Z axes of the joint coordinate system of the end joint are The directions are the same, and the Z axis of the probe coordinate system and the Z axis of the joint coordinate system of the end joint are collinear.
  • the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the moving platform coordinate system, and the Z axis of the probe coordinate system is the same.
  • the axis and the Z axis of the moving platform coordinate system are collinear.
  • FIG. 6 is a schematic diagram of a probe coordinate system and a Stewart calculation coordinate system of an ultrasonic detection device according to a preferred embodiment of the present application.
  • a rotary motor is arranged between the ultrasonic detection device and the Stewart moving platform, and the rotation angle of the motor is recorded as ⁇ m .
  • the end point of the ultrasonic detection device is located on the Z axis of the Stewart moving platform coordinate system. No matter what the value of ⁇ m is, the position of the end point of the ultrasonic detection device relative to the Stewart moving platform remains unchanged.
  • a probe coordinate system is established at the end of the ultrasonic detection device.
  • the origin of the probe coordinate system coincides with the end point of the probe, and the motion of the rotating motor between the Stewart moving platform and the ultrasonic detection device is ignored.
  • the XYZ axis of the probe coordinate system is always the same as the Stewart moving platform.
  • the XYZ axes of the coordinate system are parallel. That is, the rotation amount ⁇ m of the motor between the Stewart moving platform and the ultrasonic detection device is regarded as the movement of the ultrasonic detection plane relative to the probe coordinate system.
  • the establishment of the probe coordinate system is not unique, and it is not necessary to convert the rotation ⁇ m of the probe relative to the Stewart moving platform into the rotation of the ultrasonic detection plane relative to the probe coordinate system.
  • the motion of the ultrasound probe can be controlled by the master operator.
  • the amount of movement in the horizontal plane of the master manipulator is mapped to the tandem robotic arm.
  • the position C coord_det_st three-dimensional column vector
  • the attitude angle of the coordinate system of the probe relative to the coordinate system calculated by Stewart can be obtained (in is a fixed value)
  • the transformation matrix T trans_det_st from the probe coordinate system to the Stewart calculation coordinate system can be calculated:
  • R coord_det_st is the rotation matrix, using the attitude angle Expressed as:
  • the operator can control the mechanical arm with the ultrasonic detection device installed at the end through the main operator, so that the puncture target is located on the ultrasonic detection plane (sound beam plane).
  • the ultrasonic detection plane is a plane that passes through the end point of the ultrasonic detection device and is perpendicular to the XY plane of the probe coordinate system. At any point in the ultrasonic detection plane, the distance l c relative to the probe coordinate system and the included angle with the Z axis can be measured. ⁇ c .
  • C coord_tar_det [l c sin ⁇ c cos ⁇ m l c sin ⁇ c sin ⁇ m l c cos ⁇ c 1] T .
  • the coordinate C coord_tar_mach of the target point in the mechanical coordinate system can be obtained by coordinate transformation:
  • C coord_tar_mach T trans_mach_st T trans_st_det C coord_tar_det .
  • the ultrasonic positioning and puncturing system further includes a second serial robotic arm and a puncturing instrument mounted on the end of the second serial robotic arm; the computer equipment is also electrically connected to the second serial robotic arm, and the processor is running
  • the computer program also executes the following steps: taking the second position information of the target in the reference coordinate system as the target puncture point, and controlling the puncture instrument of the second serial robotic arm to puncture the target puncture point.
  • the coordinate C coord_tar_stl of the target point in the coordinate system calculated by Stewart of the puncture manipulator can be obtained by coordinate transformation:
  • T trans_stl_mach is the transformation matrix from the calculated coordinate system of the puncture manipulator Stewart to the calculated coordinate system.
  • the joint motion of the Stewart platform of the puncture manipulator can be calculated through the inverse kinematics of the Stewart parallel platform. Arm for precise puncture.
  • the end of the first serial robotic arm is equipped with an ultrasonic detection device
  • the end of the second serial robotic arm is equipped with a puncture instrument
  • the execution of the puncture and ablation operation is completed by the cooperation of the two robotic arms, and the first serial
  • the arm holds the ultrasound probe, scans the surface of the human body to find the lesions, and locates the puncture target, which can achieve precise positioning of the target.
  • the second serial robotic arm holds the ablation needle and performs puncture and ablation on the target point, which can achieve precise puncture of the target point.
  • the embodiments of the present application further provide a computer-readable storage medium for implementation.
  • a computer program is stored on the computer-readable storage medium; the computer program is configured to perform the following steps when executed.
  • S1 The first position information of the target point is located in the ultrasonic detection plane projected by the ultrasonic detection device.
  • S2 Acquire the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator.
  • S4 Determine the second position information of the target point in the reference coordinate system according to the position conversion relationship and the first position information.
  • locating the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device includes: acquiring an ultrasonic image collected by the ultrasonic detection device, and locating the position information of the target point in the ultrasonic image; The position information of the point in the ultrasonic image and the calibration parameters of the ultrasonic detection device determine the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device.
  • the first position information includes the distance of the target point relative to the end point, and the angle between the connecting line between the target point and the end point and the rotation axis.
  • determining the position conversion relationship between the reference coordinate system of the ultrasonic positioning and puncturing system and the coordinate system where the ultrasonic detection plane is located includes: when the first serial robotic arm is in an initial state , establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the initial state, establish the probe coordinate system with the end point of the ultrasonic detection device as the coordinate origin; Obtain the reference coordinate system of the first serial manipulator and the DH parameters of the joint coordinate system of each joint of the first serial manipulator, and determine the reference coordinate system of the first serial manipulator to the end of the first serial manipulator according to the DH parameters
  • the first position conversion relationship between the joint coordinate systems of the joints obtain the position information of the coordinate origin of the probe coordinate system in the joint coordinate system of the end joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the end joint, and According to the position information, the attitude angle
  • the first serial manipulator includes a rotating joint and a moving joint
  • the Z axis of the joint coordinate system of the rotating joint is set along the rotating axis
  • the Z axis of the joint coordinate system of the moving joint is set along the moving direction
  • the reference coordinate system The joint coordinate system of each joint is the same as the left-handed or right-handed system, and when the rotary joint is connected to the previous mobile joint, the origin of the joint coordinate system of the rotary joint coincides with the origin of the joint coordinate system of the previous mobile joint .
  • the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the joint coordinate system of the end joint, and the Z axis of the probe coordinate system and the joint coordinates of the end joint
  • the Z axis of the system is collinear.
  • the computer program is further configured to execute the following steps when running: using the second position information of the target point in the reference coordinate system as the target puncture point, and controlling the puncture instrument of the second serial robotic arm to puncture the target puncture point Perform puncture.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

An ultrasonic positioning puncture system and a storage medium. The ultrasonic positioning puncture system comprises a computer device (101), a first series mechanical arm (102) and an ultrasonic detection apparatus (103). The computer device (101) comprises a processor (1012), and when running a computer program, the processor (1012) executes the following steps: acquiring first attitude information of an ultrasonic detection plane, second attitude information of the first series mechanical arm (102), and first position information of a target point in the ultrasonic detection plane; determining, according to the first attitude information and the second attitude information, the position conversion relationship between a reference coordinate system of the ultrasonic positioning puncture system and a coordinate system in which the ultrasonic detection plane is located (S203); and determining, according to the position conversion relationship and the first position information, second position information of the target point in the reference coordinate system (S204), thereby improving the accuracy of a puncture surgery.

Description

超声定位穿刺系统和存储介质Ultrasound positioning puncture system and storage medium 技术领域technical field
本申请涉及机械臂控制领域,特别是涉及超声定位穿刺系统和存储介质。The present application relates to the field of robotic arm control, and in particular, to an ultrasonic positioning puncture system and a storage medium.
背景技术Background technique
传统的穿刺手术是医生在超声、CT、MRI等图像扫描设备的引导下,由医生手动将小型手术器械(如活检穿刺针等)送入患者体内,对病变部位进行检测或治疗,医生通过病灶附近的二维或三维扫描图像判断合适的入针点和入针方向,然后凭借经验手动完成穿刺操作。In the traditional puncture operation, under the guidance of image scanning equipment such as ultrasound, CT, and MRI, the doctor manually sends small surgical instruments (such as biopsy needles) into the patient's body to detect or treat the lesion. The nearby 2D or 3D scan images determine the appropriate needle entry point and needle entry direction, and then manually complete the puncture operation based on experience.
计算机辅助导航技术和机器人技术被引入到穿刺手术中。具有手术导航系统的穿刺机器人工作原理如下:首先由计算机将扫描的二维图像进行三维合成形成病灶附近的三维图像,然后由医生通过三维图像判断靶点位置和合适的入针路径并输入导航系统,再由导航系统计算出机器人操作臂的当前状态和目标点并规划出一条轨迹,最后操作臂按照规划的轨迹完成穿刺定位,再通过进针机构或医生手动完成进针,避免了纯手动操作带来的误差。然而,现有的导航技术通常通过光学定位对穿刺针和医学影像设备进行定位和配准,在医学影像的辅助下进行穿刺,但光学定位和配准的过程中存在较大的误差,并且在操作的过程中不能使光学小球被遮挡,限制了手术的实施。Computer-aided navigation techniques and robotics were introduced into puncture procedures. The working principle of the puncture robot with the surgical navigation system is as follows: first, the computer 3D synthesizes the scanned 2D image to form a 3D image near the lesion, and then the doctor judges the target position and the appropriate needle insertion path through the 3D image and inputs it into the navigation system. , and then the navigation system calculates the current state and target point of the robot manipulator arm and plans a trajectory. Finally, the manipulator arm completes the puncture positioning according to the planned trajectory, and then manually completes the needle insertion through the needle insertion mechanism or the doctor, avoiding pure manual operation. errors caused. However, the existing navigation technology usually locates and registers the puncture needle and medical imaging equipment through optical positioning, and performs puncture with the aid of medical images, but there are large errors in the process of optical positioning and registration, and in During the operation, the optical ball cannot be blocked, which limits the implementation of the operation.
发明内容SUMMARY OF THE INVENTION
根据本申请的各种实施例,提供一种超声定位穿刺系统,所述超声定位穿刺系统包括计算机设备、第一串联机械臂和搭载在所述第一串联机械臂末端的超声探测装置,所述计算机设备分别与所述第一串联机械臂和所述超声探测装置电性连接,所述计算机设备包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器在运行所述计算机程序时执行下列步骤:在所述超声探测装置投射的超声探测平面内定位靶点的第一位置信息;获取所述超声探测平面的第一姿态信息以及所述第一串联机械臂的第二姿态信息;根据所述第一姿态信息和所述第二姿态信息,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系;根据所述位置转换关系和所述第一位置信息,确定所述靶点在所述参考坐标系内的第二位置信息。According to various embodiments of the present application, an ultrasonic positioning and puncturing system is provided. The ultrasonic positioning and puncturing system includes a computer device, a first serial manipulator arm, and an ultrasonic detection device mounted on the end of the first serial manipulator arm. Computer equipment is electrically connected to the first serial robotic arm and the ultrasonic detection device respectively, the computer equipment includes a memory and a processor, the memory stores a computer program, and the processor runs the computer program Perform the following steps: locate the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device; obtain the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial robot arm ; According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located; According to the position conversion relationship and the described position conversion relationship The first position information determines the second position information of the target point in the reference coordinate system.
上述超声定位穿刺系统具有以下优点:通过直接在超声影像中获取病灶坐标,而不需 要通过光学定位等配准基准间接配准,解决了现有的医学影像导航下穿刺手术的定位和配准技术也存在误差的问题,提高了穿刺手术的精度。The above ultrasonic positioning puncture system has the following advantages: by directly obtaining the coordinates of the lesion in the ultrasonic image, without the need for indirect registration through registration datums such as optical positioning, it solves the existing positioning and registration technology for puncture surgery under medical image navigation There is also the problem of error, which improves the accuracy of the puncture operation.
在其中的一些实施例中,所述在所述超声探测装置投射的超声探测平面内定位靶点的第一位置信息包括:获取由所述超声探测装置采集的超声图像,在所述超声图像中定位所述靶点的位置信息;根据所述靶点在所述超声图像中的位置信息以及所述超声探测装置的标定参数,确定所述靶点在所述超声探测装置投射的超声探测平面内的所述第一位置信息。In some of these embodiments, the locating the first position information of the target point in the ultrasound detection plane projected by the ultrasound detection device includes: acquiring an ultrasound image acquired by the ultrasound detection device, in which the ultrasound image is Locate the position information of the target point; according to the position information of the target point in the ultrasonic image and the calibration parameters of the ultrasonic detection device, determine that the target point is in the ultrasonic detection plane projected by the ultrasonic detection device of the first location information.
在其中的一些实施例中,所述超声探测装置与所述第一串联机械臂的末端可转动地连接,所述超声探测平面形成于所述超声探测装置的末端端点并自所述超声探测装置的末端端点沿扇面延伸;所述超声探测平面跟随所述超声探测装置旋转,所述超声探测平面的旋转轴经过所述末端端点;所述第一姿态信息包括所述超声探测平面的旋转角度。In some of these embodiments, the ultrasonic detection device is rotatably connected to the end of the first serial manipulator, and the ultrasonic detection plane is formed at the end of the ultrasonic detection device and extends from the ultrasonic detection device. The end point of the ultrasonic detection plane extends along the fan plane; the ultrasonic detection plane rotates with the ultrasonic detection device, and the rotation axis of the ultrasonic detection plane passes through the end point; the first attitude information includes the rotation angle of the ultrasonic detection plane.
在其中的一些实施例中,所述第一位置信息包括所述靶点相对于所述末端端点的距离,以及所述靶点和所述末端端点的连接线与所述旋转轴的夹角。In some of the embodiments, the first position information includes the distance of the target point relative to the end point, and the included angle of the connecting line between the target point and the end point and the rotation axis.
在其中的一些实施例中,所述第一串联机械臂包括多个关节,多个所述关节包括:旋转关节和/或移动关节,所述第二姿态信息包括每个所述关节的运动量信息。In some of these embodiments, the first serial robotic arm includes a plurality of joints, the plurality of joints include: rotation joints and/or movement joints, and the second posture information includes motion amount information of each of the joints .
在其中的一些实施例中,所述根据所述第一姿态信息和所述第二姿态信息,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系包括:在所述第一串联机械臂处于初始状态时,建立所述第一串联机械臂的参考坐标系和所述第一串联机械臂的各关节的关节坐标系,以及在所述超声探测装置处于初始状态时,以所述超声探测装置的末端端点为坐标原点建立探头坐标系;获取所述第一串联机械臂的参考坐标系和所述第一串联机械臂的各关节的关节坐标系的DH参数,并根据所述DH参数,确定所述第一串联机械臂的参考坐标系到所述第一串联机械臂的末端关节的关节坐标系之间的第一位置转换关系;获取所述探头坐标系的坐标原点在所述末端关节的关节坐标系内的位置信息,以及所述探头坐标系相对于所述末端关节的关节坐标系的姿态角,并根据该位置信息、姿态角以及所述第一位置转换关系,确定所述参考坐标系与所述探头坐标系之间的第二位置转换关系;获取所述超声探测平面的旋转角度,并根据所述旋转角度确定所述超声探测平面所在坐标系与所述探头坐标系之间的第三位置转换关系;根据所述第二位置转换关系和所述第三位置转换关系,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系。In some of the embodiments, the determining, according to the first attitude information and the second attitude information, the position conversion relationship between the reference coordinate system of the ultrasound positioning and puncturing system and the coordinate system where the ultrasound detection plane is located includes: : when the first serial manipulator is in the initial state, establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the In the initial state, a probe coordinate system is established with the end point of the ultrasonic detection device as the coordinate origin; the reference coordinate system of the first serial manipulator and the DH of the joint coordinate system of each joint of the first serial manipulator are obtained. parameters, and according to the DH parameters, determine the first position conversion relationship between the reference coordinate system of the first serial manipulator and the joint coordinate system of the end joint of the first serial manipulator; obtain the probe coordinates The position information of the coordinate origin of the system in the joint coordinate system of the terminal joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the terminal joint, and according to the position information, attitude angle and the first a position conversion relationship, determining a second position conversion relationship between the reference coordinate system and the probe coordinate system; acquiring the rotation angle of the ultrasonic detection plane, and determining the coordinates of the ultrasonic detection plane according to the rotation angle The third position conversion relationship between the system and the probe coordinate system; according to the second position conversion relationship and the third position conversion relationship, determine the reference coordinate system of the ultrasonic positioning and puncturing system and the ultrasonic detection plane. The position conversion relationship of the coordinate system where it is located.
在其中的一些实施例中,所述第一串联机械臂包括旋转关节和移动关节,旋转关节的关节坐标系的Z轴沿旋转轴设置,移动关节的关节坐标系的Z轴沿移动方向设置,参考坐标系与各关节的关节坐标系同为左手系或右手系,且当旋转关节与在前的移动关节连接时,旋 转关节的关节坐标系的原点与该在前的移动关节的关节坐标系的原点重合。In some of the embodiments, the first serial manipulator includes a rotating joint and a moving joint, the Z axis of the joint coordinate system of the rotating joint is set along the rotating axis, and the Z axis of the joint coordinate system of the moving joint is set along the moving direction, The reference coordinate system and the joint coordinate system of each joint are both left-handed or right-handed systems, and when the rotary joint is connected to the previous mobile joint, the origin of the joint coordinate system of the rotary joint is the joint coordinate system of the previous mobile joint. The origin coincides.
在其中的一些实施例中,所述探头坐标系的X、Y、Z轴的方向与所述末端关节的关节坐标系的X、Y、Z轴的方向相同,且所述探头坐标系的Z轴和所述末端关节的关节坐标系的Z轴共线。In some of these embodiments, the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the joint coordinate system of the distal joint, and the Z axis of the probe coordinate system The axis is collinear with the Z-axis of the joint coordinate system of the end joint.
在其中的一些实施例中,所述第一串联机械臂的末端关节为并联运动平台,所述并联运动平台包括静平台、动平台和多个设置在所述静平台与所述动平台之间的伸缩组件,所述超声探测装置设于所述动平台。In some of the embodiments, the end joint of the first serial manipulator is a parallel motion platform, and the parallel motion platform includes a static platform, a moving platform, and a plurality of pieces disposed between the static platform and the moving platform The telescopic assembly, the ultrasonic detection device is set on the moving platform.
在其中的一些实施例中,所述超声定位穿刺系统还包括第二串联机械臂和搭载在所述第二串联机械臂末端的穿刺器械;所述计算机设备还分别与所述第二串联机械臂电性连接,所述处理器在运行所述计算机程序时还执行下列步骤:将所述靶点在所述参考坐标系内的第二位置信息作为目标穿刺点,控制所述第二串联机械臂的穿刺器械对所述目标穿刺点进行穿刺。In some of the embodiments, the ultrasonic positioning and puncturing system further includes a second serial robotic arm and a puncturing instrument mounted on the end of the second serial robotic arm; the computer equipment is further connected to the second serial robotic arm, respectively. electrically connected, the processor further executes the following steps when running the computer program: taking the second position information of the target point in the reference coordinate system as the target puncture point, and controlling the second serial robotic arm The puncture instrument punctures the target puncture point.
根据本申请的各种实施例,提供一种存储介质,所述存储介质中存储有计算机程序,所述计算机程序被设置为运行时执行下列步骤:在所述超声探测装置投射的超声探测平面内定位靶点的第一位置信息;获取所述超声探测平面的第一姿态信息以及所述第一串联机械臂的第二姿态信息;根据所述第一姿态信息和所述第二姿态信息,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系;根据所述位置转换关系和所述第一位置信息,确定所述靶点在所述参考坐标系内的第二位置信息。According to various embodiments of the present application, there is provided a storage medium in which a computer program is stored, the computer program being configured to execute the following steps when running: in an ultrasound detection plane projected by the ultrasound detection device Locating the first position information of the target point; obtaining the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator; according to the first attitude information and the second attitude information, determine The position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located; according to the position conversion relationship and the first position information, it is determined that the target point is in the reference coordinate system the second location information.
附图说明Description of drawings
为了更好地描述和说明这里公开的那些发明的实施例和/或示例,可以参考一幅或多幅附图。用于描述附图的附加细节或示例不应当被认为是对所公开的发明、目前描述的实施例和/或示例以及目前理解的这些发明的最佳模式中的任何一者的范围的限制。In order to better describe and illustrate embodiments and/or examples of those inventions disclosed herein, reference may be made to one or more of the accompanying drawings. The additional details or examples used to describe the drawings should not be construed as limiting the scope of any of the disclosed inventions, the presently described embodiments and/or examples, and the best mode presently understood of these inventions.
图1是本申请实施例提供的超声定位穿刺系统的结构框图。FIG. 1 is a structural block diagram of an ultrasonic positioning and puncturing system provided by an embodiment of the present application.
图2是本申请实施例提供的超声定位穿刺系统的工作方法的流程图。FIG. 2 is a flowchart of a working method of the ultrasonic positioning puncture system provided by the embodiment of the present application.
图3是本申请实施例的多关节机械臂的坐标转换方法的流程图。FIG. 3 is a flowchart of a coordinate conversion method for a multi-joint robotic arm according to an embodiment of the present application.
图4是本申请优选实施例提供的多关节机械臂的结构示意图。FIG. 4 is a schematic structural diagram of a multi-joint robotic arm provided by a preferred embodiment of the present application.
图5是本申请优选实施例提供的多关节机械臂的机械坐标系和关节坐标系的示意图。FIG. 5 is a schematic diagram of a mechanical coordinate system and a joint coordinate system of a multi-joint manipulator provided by a preferred embodiment of the present application.
图6是本申请优选实施例的超声探测装置的探头坐标系与Stewart计算坐标系的示意图。FIG. 6 is a schematic diagram of a probe coordinate system and a Stewart calculation coordinate system of an ultrasonic detection device according to a preferred embodiment of the present application.
具体实施方式Detailed ways
以下将对本申请提供的超声定位穿刺系统和存储介质作进一步说明。The ultrasonic positioning puncture system and storage medium provided by the present application will be further described below.
图1是本申请实施例提供的超声定位穿刺系统的结构框图,如图1所示,该超声定位穿刺系统包括计算机设备101、第一串联机械臂102和搭载在第一串联机械臂102末端的超声探测装置103,计算机设备101分别与第一串联机械臂102和超声探测装置103电性连接,计算机设备101包括存储器1011和处理器1012,存储器1011中存储有计算机程序10111。FIG. 1 is a structural block diagram of an ultrasonic positioning and puncturing system provided by an embodiment of the present application. As shown in FIG. 1 , the ultrasonic positioning and puncturing system includes a computer device 101 , a first serial robotic arm 102 , and a device mounted on the end of the first serial robotic arm 102 . The ultrasonic detection device 103 and the computer equipment 101 are electrically connected to the first serial manipulator 102 and the ultrasonic detection device 103 respectively. The computer equipment 101 includes a memory 1011 and a processor 1012, and the memory 1011 stores a computer program 10111.
图2是本申请实施例提供的超声定位穿刺系统的工作方法的流程图,如图2所示,处理器1012在运行计算机程序时执行下列步骤S201至步骤S204。FIG. 2 is a flowchart of the working method of the ultrasonic positioning and puncturing system provided by the embodiment of the present application. As shown in FIG. 2 , the processor 1012 executes the following steps S201 to S204 when running a computer program.
步骤S201:在超声探测装置投射的超声探测平面内定位靶点的第一位置信息。Step S201: Locating the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device.
超声探测装置103投射的超声探测平面通常呈扇面状,扇面的扇心即为超声发生的起点位置,在扇心位置发出超声并接收物质反射的回波来探测该超声探测平面内的目标位置。The ultrasonic detection plane projected by the ultrasonic detection device 103 is usually in the shape of a fan, and the fan center of the fan is the starting point of the ultrasonic generation. The ultrasonic detection is sent at the fan center and the echoes reflected by the substance are received to detect the target position in the ultrasonic detection plane.
由于超声探测平面只能够探测到一个扇面的信息,而待穿刺的靶点位置位于人体内的一个空间位置内,因此若要探测到准确的靶点位置,搭载超声探测装置103的第一串联机械臂102需要将超声探测装置103移动到靶点附近位置,在移动到靶点附近位置时,第一串联机械臂102的各关节的运动量会发生变化,在本申请实施例中这些关节运动量的集合称为第二姿态信息。在其中的一些实施例中,第一串联机械臂包括多个关节,多个关节包括:旋转关节和/或移动关节,第二姿态信息包括每个关节的运动量信息。Since the ultrasonic detection plane can only detect the information of one sector, and the position of the target to be punctured is located in a spatial position in the human body, to detect the accurate target position, the first serial machine equipped with the ultrasonic detection device 103 The arm 102 needs to move the ultrasonic detection device 103 to a position near the target point. When moving to a position near the target point, the movement amount of each joint of the first serial robotic arm 102 will change. It is called the second attitude information. In some of the embodiments, the first serial robotic arm includes a plurality of joints, the plurality of joints include: a rotation joint and/or a movement joint, and the second posture information includes motion amount information of each joint.
在其中的一些实施例中,超声探测装置与第一串联机械臂102的末端可转动地连接,超声探测平面形成于超声探测装置的末端端点并自超声探测装置的末端端点沿扇面延伸;超声探测平面跟随超声探测装置旋转,超声探测平面的旋转轴经过末端端点;第一姿态信息包括超声探测平面的旋转角度。In some of these embodiments, the ultrasonic detection device is rotatably connected to the end of the first serial robotic arm 102, and the ultrasonic detection plane is formed at the end of the ultrasonic detection device and extends from the end of the ultrasonic detection device along the sector; ultrasonic detection The plane rotates following the ultrasonic detection device, and the rotation axis of the ultrasonic detection plane passes through the end point; the first attitude information includes the rotation angle of the ultrasonic detection plane.
当超声探测装置103被移动到靶点附近位置之后,超声探测平面朝向靶点方向照射,为了在超声探测平面中显示岀靶点,在本申请实施例中沿扇面的对称轴旋转超声探测装置103,超声探测平面随之旋转,从而实现圆台形或圆锥形空间区域内的靶点探测,每个超声探测平面的超声图像都会实时显示在电子屏幕上。一旦在电子屏幕上观察到与靶点特征相符的图像,则可以该靶点在超声图像内的位置确定靶点在真实的超声探测平面内的位置,即为第一位置信息。After the ultrasonic detection device 103 is moved to a position near the target point, the ultrasonic detection plane is irradiated toward the target point. In order to display the target point in the ultrasonic detection plane, in the embodiment of the present application, the ultrasonic detection device 103 is rotated along the symmetry axis of the fan plane. , the ultrasonic detection plane rotates accordingly, so as to realize the target detection in the truncated or conical space area, and the ultrasonic image of each ultrasonic detection plane will be displayed on the electronic screen in real time. Once an image consistent with the characteristics of the target point is observed on the electronic screen, the position of the target point in the real ultrasound detection plane can be determined from the position of the target point in the ultrasound image, which is the first position information.
例如,确定靶点在超声探测平面内的位置,通过实时获取由超声探测装置采集的超声图像,在超声图像中定位靶点的位置信息;然后根据靶点在超声图像中的位置信息以及超声探测装置的标定参数,确定靶点在超声探测装置投射的超声探测平面内的第一位置信息。超声探测装置的标定参数包括但不限于以下至少之一:阵元数、曲率半径、线阵宽度、查扫视 野(FOV)等。For example, determine the position of the target point in the ultrasonic detection plane, and locate the position information of the target point in the ultrasonic image by acquiring the ultrasonic image collected by the ultrasonic detection device in real time; then according to the position information of the target point in the ultrasonic image and the ultrasonic detection The calibration parameters of the device determine the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device. The calibration parameters of the ultrasonic detection device include, but are not limited to, at least one of the following: the number of array elements, the radius of curvature, the width of the linear array, the field of view (FOV), and the like.
第一位置信息可以采用直角坐标或极坐标的方式表示。例如,在采用极坐标方式表示的情况下,第一位置信息包括靶点相对于末端端点的距离,以及靶点和末端端点的连接线与旋转轴的夹角。The first position information may be represented in the form of Cartesian coordinates or polar coordinates. For example, in the case of using polar coordinates to represent, the first position information includes the distance of the target point relative to the end point, and the angle between the connecting line between the target point and the end point and the rotation axis.
步骤S202:获取超声探测平面的第一姿态信息以及第一串联机械臂的第二姿态信息。Step S202: Acquire first attitude information of the ultrasonic detection plane and second attitude information of the first serial manipulator.
步骤S203:根据第一姿态信息和第二姿态信息,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系。Step S203: According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located.
步骤S204:根据位置转换关系和第一位置信息,确定靶点在参考坐标系内的第二位置信息。Step S204: Determine the second position information of the target point in the reference coordinate system according to the position conversion relationship and the first position information.
在超声探测平面的第一姿态信息和第一串联机械臂的第二姿态信息确定的情况下,则超声探测平面上的任意一点与参考坐标系的转换关系就是确定的,因此,根据第一姿态信息和第二姿态信息,可以确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系,再根据该位置转换关系和第一位置信息就能够确定靶点在参考坐标系内的第二位置信息,实现靶点在参考坐标系下的定位。When the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator are determined, the conversion relationship between any point on the ultrasonic detection plane and the reference coordinate system is determined. Therefore, according to the first attitude information and the second attitude information, can determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located, and then according to the position conversion relationship and the first position information, it can be determined that the target point is in the reference coordinate system. The second position information of , realizes the positioning of the target point in the reference coordinate system.
需要说明的是,在上述流程中或者附图的流程图中示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。例如,步骤S201也可以在步骤S202或步骤S203之后执行。It should be noted that the steps shown in the above flow or the flow chart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical sequence is shown in the flow chart, in the In some cases, steps shown or described may be performed in an order different from that herein. For example, step S201 may also be performed after step S202 or step S203.
在其中的一些实施例中,根据第一姿态信息和第二姿态信息,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系包括如下步骤。In some of the embodiments, according to the first attitude information and the second attitude information, determining the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located includes the following steps.
步骤S203-1:在第一串联机械臂处于初始状态时,建立第一串联机械臂的参考坐标系和第一串联机械臂的各关节的关节坐标系,以及在超声探测装置处于初始状态时,以超声探测装置的末端端点为坐标原点建立探头坐标系。Step S203-1: When the first serial manipulator is in the initial state, establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the initial state, The probe coordinate system is established with the end point of the ultrasonic detection device as the coordinate origin.
步骤S203-2:获取第一串联机械臂的参考坐标系和第一串联机械臂的各关节的关节坐标系的DH参数,并根据DH参数,确定第一串联机械臂的参考坐标系到第一串联机械臂的末端关节的关节坐标系之间的第一位置转换关系。Step S203-2: Obtain the reference coordinate system of the first serial manipulator and the DH parameters of the joint coordinate system of each joint of the first serial manipulator, and determine the reference coordinate system of the first serial manipulator to the first serial manipulator according to the DH parameter. The first position transformation relationship between the joint coordinate systems of the end joints of the serial manipulator.
步骤S203-3:获取探头坐标系的坐标原点在末端关节的关节坐标系内的位置信息,以及探头坐标系相对于末端关节的关节坐标系的姿态角,并根据该位置信息、姿态角以及第一位置转换关系,确定参考坐标系与探头坐标系之间的第二位置转换关系。Step S203-3: Obtain the position information of the coordinate origin of the probe coordinate system in the joint coordinate system of the terminal joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the terminal joint, and according to the position information, attitude angle and the first A position conversion relationship, which determines a second position conversion relationship between the reference coordinate system and the probe coordinate system.
步骤S203-4:获取超声探测平面的旋转角度,并根据旋转角度确定超声探测平面所在坐标系与探头坐标系之间的第三位置转换关系。Step S203-4: Obtain the rotation angle of the ultrasonic detection plane, and determine a third position conversion relationship between the coordinate system where the ultrasonic detection plane is located and the probe coordinate system according to the rotation angle.
步骤S203-5:根据第二位置转换关系和第三位置转换关系,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系。Step S203-5: According to the second position conversion relationship and the third position conversion relationship, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning and puncturing system and the coordinate system where the ultrasonic detection plane is located.
上述的步骤S203-2旨在确定第一串联机械臂的参考坐标系到第一串联机械臂的末端关节的关节坐标系之间的第一位置转换关系。The above-mentioned step S203-2 aims to determine the first position conversion relationship between the reference coordinate system of the first serial manipulator and the joint coordinate system of the end joint of the first serial manipulator.
针对相关技术中机械臂的控制复杂的问题,在本实施例提供了一种多关节机械臂的坐标转换方法,应用于包括旋转关节和移动关节的机械臂,这些旋转关节和移动关节依次串联连接。为了实现机械臂在空间中多个自由度的运动,旋转关节和移动关节的数量可以根据需要的自由度和运动范围进行设置,即旋转关节和移动关节的数量都可以是一个或者多个。Aiming at the problem of complex control of the manipulator in the related art, this embodiment provides a coordinate conversion method for a multi-joint manipulator, which is applied to a manipulator including a rotating joint and a moving joint, and these rotating joints and moving joints are sequentially connected in series . In order to realize the movement of the manipulator with multiple degrees of freedom in space, the number of rotating joints and moving joints can be set according to the required degrees of freedom and motion range, that is, the number of rotating joints and moving joints can be one or more.
图3是本申请实施例的多关节机械臂的坐标转换方法的流程图,如图3所示,该流程包括如下步骤S301至步骤S303。FIG. 3 is a flowchart of a coordinate conversion method for a multi-joint robotic arm according to an embodiment of the present application. As shown in FIG. 3 , the flowchart includes the following steps S301 to S303 .
步骤S301:在机械臂处于初始状态时,建立机械臂的机械坐标系和机械臂的各关节的关节坐标系。Step S301: When the robot arm is in an initial state, a mechanical coordinate system of the robot arm and a joint coordinate system of each joint of the robot arm are established.
机械臂的控制通常基于DH参数进行坐标系的转换,DH参数包括传统的DH参数和改进的DH参数。相连接的两个关节为相邻关节。相邻两个关节的关节坐标系的变换通常采用DH参数或改进的DH参数来表示。以DH参数为例,相邻两个关节坐标系通过绕Z轴旋转θ、平移d后,再绕X轴旋转α、平移a,就能够相互重合,上述的θ、d、α和a即为DH参数。由此可知,如果DH参数越简洁,则相邻两个关节坐标系的转换越简单。The control of the manipulator is usually based on DH parameters for coordinate system transformation, and DH parameters include traditional DH parameters and improved DH parameters. Two connected joints are adjacent joints. The transformation of the joint coordinate system of two adjacent joints is usually represented by a DH parameter or an improved DH parameter. Taking the DH parameter as an example, two adjacent joint coordinate systems can overlap each other by rotating θ around the Z axis and translating d, and then rotating α around the X axis and translating a. The above θ, d, α and a are DH parameters. It can be seen from this that the simpler the DH parameter is, the simpler the transformation of two adjacent joint coordinate systems will be.
为了简化DH参数,在本实施例中,旋转关节的关节坐标系的Z轴沿旋转轴设置,移动关节的关节坐标系的Z轴沿移动方向设置,机械坐标系与各关节的关节坐标系同为左手系或右手系,且当旋转关节之前的关节为移动关节时,旋转关节的关节坐标系的原点与该移动关节的关节坐标系的原点重合。In order to simplify the DH parameters, in this embodiment, the Z axis of the joint coordinate system of the rotating joint is set along the rotation axis, the Z axis of the joint coordinate system of the moving joint is set along the moving direction, and the mechanical coordinate system is the same as the joint coordinate system of each joint. It is a left-handed system or a right-handed system, and when the joint before the rotating joint is a moving joint, the origin of the joint coordinate system of the rotating joint coincides with the origin of the joint coordinate system of the moving joint.
上述的机械坐标系又可以称为参考坐标系或者世界坐标系,该机械坐标系通常设置于机械臂的底座中心,并且位于多关节臂的第一个关节的正下方,这样可以使得机械坐标系与第一个关节的转换关系尽可能地简单。The above-mentioned mechanical coordinate system can also be called a reference coordinate system or a world coordinate system. The mechanical coordinate system is usually set at the center of the base of the manipulator, and is located directly below the first joint of the multi-joint arm, which can make the mechanical coordinate system The transformation relationship to the first joint is as simple as possible.
在串联机械臂中,在大多数情况下旋转关节在接收到转角为0或2π的指令时,可能不需要区分二者,而是根据预先的设定保持不旋转或者以设定方向旋转2π。但在一些情况下,需要对这两个转角进行区分。在需要对0或2π转角进行区分的情形下,旋转关节的DH参数中Z轴的旋转角度不为0或2π,以避免混淆0或2π的旋转角度。In a serial manipulator, in most cases, when a rotary joint receives a command with a rotation angle of 0 or 2π, it may not be necessary to distinguish between the two, but it will not rotate according to a preset setting or rotate by 2π in a set direction. However, in some cases, it is necessary to distinguish between these two corners. In the case where the 0 or 2π rotation angle needs to be distinguished, the rotation angle of the Z axis in the DH parameter of the rotary joint is not 0 or 2π, so as to avoid confusing the rotation angle of 0 or 2π.
步骤S302:获取机械臂的机械坐标系和机械臂的各关节的关节坐标系的DH参数,并根据DH参数,确定机械臂的机械坐标系到机械臂的末端关节的关节坐标系之间的转换关系。Step S302: Obtain the mechanical coordinate system of the manipulator and the DH parameters of the joint coordinate system of each joint of the manipulator, and determine the conversion between the mechanical coordinate system of the manipulator and the joint coordinate system of the end joint of the manipulator according to the DH parameters relation.
在建立各关节的关节坐标系之后,就可以依次获得相邻关节之间的DH参数,并根据DH参数,确定机械臂的机械坐标系到机械臂的末端关节的关节坐标系之间的转换关系。After the joint coordinate system of each joint is established, the DH parameters between adjacent joints can be obtained in turn, and according to the DH parameters, the conversion relationship between the mechanical coordinate system of the manipulator and the joint coordinate system of the end joint of the manipulator can be determined. .
其中,根据DH法则,第i-1关节的坐标到第i关节的坐标的齐次变换被构造为具有两个旋转和两个变换的序列,采用矩阵可表示如下:Among them, according to the DH rule, the homogeneous transformation from the coordinates of the i-1th joint to the coordinates of the ith joint is constructed as a sequence with two rotations and two transformations, which can be expressed as follows using a matrix:
Figure PCTCN2020141231-appb-000001
Figure PCTCN2020141231-appb-000001
其中,i=2,3,4,…,n;n为机械臂的旋转关节和移动关节的总数量。机械臂的第1个关节的DH参数是表示的该第一个关节与机械坐标系之间的坐标系的转换,记为
Figure PCTCN2020141231-appb-000002
则机械坐标系与第1个关节的关节坐标系的转换关系:
Among them, i=2,3,4,...,n; n is the total number of rotating joints and moving joints of the robotic arm. The DH parameter of the first joint of the manipulator is the transformation of the coordinate system between the first joint and the mechanical coordinate system, which is recorded as
Figure PCTCN2020141231-appb-000002
Then the conversion relationship between the mechanical coordinate system and the joint coordinate system of the first joint:
Figure PCTCN2020141231-appb-000003
Figure PCTCN2020141231-appb-000003
可见,该变换关系与第i-1关节的坐标到第i关节的坐标的齐次变换形式完全相同。It can be seen that the transformation relationship is exactly the same as the homogeneous transformation form from the coordinates of the i-1th joint to the coordinates of the ith joint.
在得到
Figure PCTCN2020141231-appb-000004
Figure PCTCN2020141231-appb-000005
之后,从第0个坐标系(机械坐标系)到第N关节的关节坐标系的转换矩阵可表示为:
in getting
Figure PCTCN2020141231-appb-000004
and
Figure PCTCN2020141231-appb-000005
After that, the transformation matrix from the 0th coordinate system (mechanical coordinate system) to the joint coordinate system of the Nth joint can be expressed as:
Figure PCTCN2020141231-appb-000006
Figure PCTCN2020141231-appb-000006
其中,第N个关节为末端关节。Among them, the Nth joint is the end joint.
步骤S303:根据转换关系进行机械坐标系与末端关节的关节坐标系之间的坐标转换。Step S303: Perform coordinate conversion between the mechanical coordinate system and the joint coordinate system of the end joint according to the conversion relationship.
上述步骤S302中确定的
Figure PCTCN2020141231-appb-000007
表示的就是机械坐标系与末端关节的关节坐标系之间的坐标转换关系,根据该转换关系就能够实现机械坐标系与末端关节的关节坐标系之间的坐标转换。
Determined in the above step S302
Figure PCTCN2020141231-appb-000007
It represents the coordinate conversion relationship between the mechanical coordinate system and the joint coordinate system of the end joint. According to the conversion relationship, the coordinate conversion between the mechanical coordinate system and the joint coordinate system of the end joint can be realized.
下面通过优选实施例对本申请实施例进行描述和说明。The embodiments of the present application will be described and illustrated below through preferred embodiments.
图4是本申请优选实施例提供的多关节机械臂的结构示意图,如图4所示的多关节机械臂依次包括移动关节1、旋转关节2、移动关节3、旋转关节4、旋转关节5、移动关节6、旋转关节7、移动关节8、旋转关节9和移动关节10,此外,该多关节机械臂还包括与移动关节1固连的基座11。FIG. 4 is a schematic structural diagram of a multi-joint mechanical arm provided by a preferred embodiment of the present application. The multi-joint mechanical arm shown in FIG. The moving joint 6 , the rotating joint 7 , the moving joint 8 , the rotating joint 9 and the moving joint 10 , in addition, the multi-joint robotic arm further includes a base 11 fixedly connected with the moving joint 1 .
该优选实施例的多关节机械臂的坐标转换方法包括步骤1至步骤3,分别对应于图3中的步骤S301至步骤S303。The coordinate transformation method for a multi-joint robotic arm in this preferred embodiment includes steps 1 to 3, respectively corresponding to steps S301 to S303 in FIG. 3 .
步骤1:依照世界坐标系的规则在机械臂的基座处建立机械坐标系,以及各关节的关节坐标系。Step 1: According to the rules of the world coordinate system, the mechanical coordinate system and the joint coordinate system of each joint are established at the base of the manipulator.
图5是本申请优选实施例提供的多关节机械臂的机械坐标系和关节坐标系的示意图。FIG. 5 is a schematic diagram of a mechanical coordinate system and a joint coordinate system of a multi-joint manipulator provided by a preferred embodiment of the present application.
参考图5,机械坐标系F 0-X 0Y 0Z 0坐标系原点F 0固连于机械臂基座处,Z 0轴由F 0指向移动关节1,Y 0轴由基座F 0点指向机械臂,X 0轴指向符合右手坐标系。 Referring to Figure 5, the origin of the mechanical coordinate system F 0 -X 0 Y 0 Z 0 is fixed at the base of the robot arm, the Z 0 axis points from F 0 to the moving joint 1 , and the Y 0 axis points from the base F 0 Point to the robotic arm, and the X 0 axis points to conform to the right-handed coordinate system.
移动关节1的关节坐标系L 1-X 1Y 1Z 1的原点L 1固连在移动关节 1上,各坐标轴的指向和机械坐标系的对应轴指向相同。 The origin L 1 of the joint coordinate system L 1 -X 1 Y 1 Z 1 of the mobile joint 1 is fixed on the mobile joint 1 , and the orientation of each coordinate axis is the same as that of the corresponding axis of the mechanical coordinate system.
旋转关节2的关节坐标系R 2-X 2Y 2Z 2的原点R 2固连在旋转关节2上且与L 1重合,Z 2指向和Z 1轴指向相同,X 2轴和Y 2轴的指向分别与X 1轴和Y 1轴指向相反。 The origin of the joint coordinate system R 2 -X 2 Y 2 Z 2 of the rotary joint 2 is fixed on the rotary joint 2 and coincides with L 1. The Z 2 and Z 1 axes point in the same direction, the X 2 axis and the Y 2 axis The directions are opposite to the X 1 axis and Y 1 axis respectively.
移动关节3的关节坐标系L 3-X 3Y 3Z 3原点L 3固连在移动关节3上,Z 3轴由L 1点指向L 3点,X 3轴和Y 3轴分别与X 2轴和Z 2轴指向相同。 The joint coordinate system L 3 -X 3 Y 3 Z 3 of the mobile joint 3 is fixed at the origin L 3 of the mobile joint 3, the Z 3 axis points from the L 1 point to the L 3 point, and the X 3 axis and the Y 3 axis are respectively connected with the X 2 Axis and Z 2 axis point the same.
旋转关节4的关节坐标系R 4-X 4Y 4Z 4原点R 4固连在旋转关节4上且与L 3点重合(图5中为清晰标明移动关节3的关节坐标系和旋转关节4的关节坐标系将L 3和R 4分开标出,下同),Z 4轴指向和Y 3轴指向相反,初始时,X 4轴和Y 4轴指向分别与X 3轴和Z 3轴指向相反。 The joint coordinate system R 4 -X 4 Y 4 Z 4 of the rotary joint 4 is fixed at the origin R 4 of the rotary joint 4 and coincides with the point L 3 (Fig. 5 is a clear indication of the joint coordinate system of the mobile joint 3 and the rotary joint 4). The joint coordinate system of L 3 and R 4 are marked separately, the same below), the Z 4 axis points to the opposite direction to the Y 3 axis, initially, the X 4 axis and the Y 4 axis point to the X 3 axis and the Z 3 axis respectively. on the contrary.
旋转关节5的关节坐标系R 5-X 5Y 5Z 5原点R 5固连在旋转关节5上,Z 5轴方向和Z 4轴方向相同,初始时,X 5轴和Y 5轴指向分别与X 4轴和Y 4轴指向相反。 The joint coordinate system R 5 -X 5 Y 5 Z 5 of the rotary joint 5 is fixed at the origin R 5 of the rotary joint 5. The direction of the Z 5 axis and the Z 4 axis are the same. Initially, the X 5 axis and the Y 5 axis point respectively. The X 4 -axis and Y 4 -axis point opposite.
移动关节6的关节坐标系L 6-X 6Y 6Z 6原点L 6固连在移动关节6上,Z 6轴由R 5点指向L 6点,X 6轴和Y 6轴指向分别与X 5轴和Z 5轴指向相同。 The joint coordinate system L 6 -X 6 Y 6 Z 6 of the mobile joint 6 is fixed at the origin L 6 of the mobile joint 6. The Z 6 axis points from the R 5 point to the L 6 point, and the X 6 axis and the Y 6 axis point to the X axis respectively. The 5 -axis and the Z 5 -axis point the same.
旋转关节7的关节坐标系R 7-X 7Y 7Z 7原点R 7固连在旋转关节7上且和L 6点重合,初始时,Z 7轴和Y 7轴的指向分别与Y 6轴和X 6轴的指向相反,X 7轴和Z 6轴的指向相同。 The joint coordinate system R 7 -X 7 Y 7 Z 7 of the rotary joint 7 is fixed at the origin R 7 of the rotary joint 7 and coincides with the L 6 point. Initially, the Z 7 axis and the Y 7 axis point to the Y 6 axis respectively. The X 7 axis and the Z 6 axis point in the same direction, as opposed to the X 6 axis.
移动关节8的关节坐标系L 8-X 8Y 8Z 8的原点L 8固连在移动关节8上,Z 8轴由L 8点指向R 7点,X 8和Y 8轴的指向分别与X 7和Y 7轴的指向相同。 The origin of the joint coordinate system L 8 -X 8 Y 8 Z 8 of the mobile joint 8 is fixed on the mobile joint 8. The Z 8 axis points from the L 8 point to the R 7 point, and the directions of the X 8 and Y 8 axes are respectively the same as The X 7 and Y 7 axes point the same.
旋转关节9的关节坐标系R 9-X 9Y 9Z 9的原点R 9点固连在旋转关节9上且和L 8点重合,Z 9轴的方向与Y 8轴的方向相反,初始时,X 9轴的方向与Z 8轴的方向相反,Y 9轴的方向与X8轴的方向相同。 The origin of the joint coordinate system R9 - X9Y9Z9 of the rotary joint 9 is fixed on the rotary joint 9 and coincides with the L8 point. The direction of the Z9 axis is opposite to the direction of the Y8 axis. Initially , the direction of the X 9 axis is opposite to the direction of the Z 8 axis, and the direction of the Y 9 axis is the same as the direction of the X8 axis.
移动关节10的关节坐标系L 10-X 10Y 10Z 10的原点L 10固连在移动关节10上,Z 10轴由R 9点指向L 10点,X 10轴和Y 10轴的指向分别与X 9轴和Z 9轴的指向相同。 The origin of the joint coordinate system L 10 -X 10 Y 10 Z 10 of the mobile joint 10 is fixed on the mobile joint 10. The Z 10 axis points from the R 9 point to the L 10 point, and the X 10 axis and the Y 10 axis point respectively. The same point as the X 9 and Z 9 axes.
F 0R 2的长度为l 1,R 2R 4的长度为l 2,R 4R 5的长度为l 3,R 5R 7的长度为l 4,R 7R 9的长度为l 5,R 9L 10的长度为l 6,点L 1、R 2、L 3、R 4、R 5、L 6和R 7位于同一水平面上。 The length of F 0 R 2 is l 1 , the length of R 2 R 4 is l 2 , the length of R 4 R 5 is l 3 , the length of R 5 R 7 is l 4 , the length of R 7 R 9 is l 5 , The length of R 9 L 10 is l 6 , and the points L 1 , R 2 , L 3 , R 4 , R 5 , L 6 and R 7 are located on the same horizontal plane.
步骤2:获取DH参数,并计算机械坐标系到多关节机械臂的末端关节的转换关系。Step 2: Obtain the DH parameters and calculate the conversion relationship from the mechanical coordinate system to the end joints of the multi-joint robotic arm.
DH参数如表1所示。在表1中,为了避免混淆旋转关节的旋转角度0和2π,在每个旋转关节的Z轴旋转角θ上分别增加了一个偏转角度,例如旋转关节2的Z轴旋转角度为π+θ 2,其中θ 2即为偏转角度。 DH parameters are shown in Table 1. In Table 1, in order to avoid confusing the rotation angles 0 and 2π of the rotary joint, a deflection angle is added to the Z-axis rotation angle θ of each rotary joint. For example, the Z-axis rotation angle of the rotary joint 2 is π+θ 2 , where θ 2 is the deflection angle.
表1图5所示的多关节机械臂的DH参数表Table 1. DH parameter table of the multi-joint manipulator shown in Fig. 5
Figure PCTCN2020141231-appb-000008
Figure PCTCN2020141231-appb-000008
在图4所示的多关节机械臂的末端关节上还可以搭载并联平台,并联平台包括静平台、动平台和多个设置在静平台与动平台之间的伸缩组件,并联平台的计算坐标系与多关节机械臂的末端关节的关节坐标系完全重合,超声探测装置设于动平台。A parallel platform can also be mounted on the end joint of the multi-joint manipulator shown in Figure 4. The parallel platform includes a static platform, a moving platform, and a plurality of telescopic components arranged between the static platform and the moving platform. The calculation coordinate system of the parallel platform It is completely coincident with the joint coordinate system of the end joint of the multi-joint mechanical arm, and the ultrasonic detection device is arranged on the moving platform.
本实施例的并联平台可以实现多个自由度的运动,以具有六自由度的Stewart并联平台为例,Stewart并联平台包括静平台、动平台和多个设置在静平台与动平台之间的伸缩组件,可以实现空间上六个自由度的运动,分别为沿X轴位移、沿Y轴位移、沿Z轴位移、绕X轴转动、绕Y轴转动和绕Z轴转动。Stewart并联平台由6根伸缩组件支撑,与采用串联的悬臂梁结构的被动臂相比刚度大,结构稳定,并且由于刚度大,并联结构较串联结构在相同的自重或体积下,有较高的承载能力。采用串联的悬臂梁结构的被动臂末端的误差是各个关节误差的积累和放大,因而误差大、精度低,而并联平台则没有那样的误差积累和放大关系,微动精度高,更适宜执行高精度的手术操作。此外,在位置求解上,并联平台的逆解非常容易,容易根据坐标位置求得并联平台各伸缩组件的运动姿态。The parallel platform in this embodiment can realize motion with multiple degrees of freedom. Taking the Stewart parallel platform with six degrees of freedom as an example, the Stewart parallel platform includes a static platform, a moving platform, and a plurality of telescopic and retractable platforms arranged between the static platform and the moving platform. The component can realize six degrees of freedom movement in space, namely displacement along the X axis, displacement along the Y axis, displacement along the Z axis, rotation around the X axis, rotation around the Y axis, and rotation around the Z axis. The Stewart parallel platform is supported by 6 telescopic components. Compared with the passive arm using the cantilever beam structure in series, the rigidity is higher and the structure is stable. Due to the high rigidity, the parallel structure has a higher self-weight or volume than the series structure. Carrying capacity. The error at the end of the passive arm using the cantilever beam structure in series is the accumulation and amplification of the errors of each joint, so the error is large and the precision is low. Precision surgical operation. In addition, in terms of position solution, the inverse solution of the parallel platform is very easy, and it is easy to obtain the motion posture of each telescopic component of the parallel platform according to the coordinate position.
正是基于并联平台逆解非常容易的特点,在上述实施例中,可以通过在静平台上建立从用户坐标系,将主操作手的位移量映射到从用户坐标系内,再通过从用户坐标系与计算坐标系的转换从而得到目标位置在计算坐标系内的位置信息,根据该位置信息通过逆解就能够很容易地得到并联平台各伸缩组件的运动姿态,相对于相关技术中通过正运动学求解并联平台末端的位置信息,上述方式大大降低了运算复杂度,提高了控制效率,节约了运算资源。Based on the fact that the inverse solution of the parallel platform is very easy, in the above embodiment, the displacement of the master operator can be mapped to the slave user coordinate system by establishing the slave user coordinate system on the static platform, and then the slave user coordinate system can be used to map the displacement of the master operator into the slave user coordinate system. The position information of the target position in the calculation coordinate system can be obtained by the conversion of the system and the calculation coordinate system. According to the position information, the motion posture of each telescopic assembly of the parallel platform can be easily obtained through the inverse solution. The above method greatly reduces the computational complexity, improves the control efficiency, and saves computational resources.
解算Stewart计算坐标系到机械坐标系的转换矩阵,即移动关节坐标系L 10-X 10Y 10Z 10 到机械坐标系F 0-X 0Y 0Z 0的转换矩阵与逆矩阵。 Solve the transformation matrix from the Stewart-calculated coordinate system to the mechanical coordinate system, that is, the transformation matrix and inverse matrix from the moving joint coordinate system L 10 -X 10 Y 10 Z 10 to the mechanical coordinate system F 0 -X 0 Y 0 Z 0 .
根据机器人正运动学可知两相邻关节之间的转换矩阵为:According to the forward kinematics of the robot, the transformation matrix between two adjacent joints is:
Figure PCTCN2020141231-appb-000009
Figure PCTCN2020141231-appb-000009
从第m关节(包含m关节的自由度)到第n关节的转换矩阵可表示为:The transformation matrix from the mth joint (containing the m degree of freedom) to the nth joint can be expressed as:
Figure PCTCN2020141231-appb-000010
Figure PCTCN2020141231-appb-000010
可解算得单条机械臂从基座到Stewart平台静平台的转换矩阵
Figure PCTCN2020141231-appb-000011
即机械坐标系到Stewart计算坐标系的转换矩阵,命名为T trans_mach_st,而Stewart计算坐标系到机械坐标系的转换矩阵则为逆矩阵
Figure PCTCN2020141231-appb-000012
The transformation matrix of a single robotic arm from the base to the static platform of the Stewart platform can be solved
Figure PCTCN2020141231-appb-000011
That is, the transformation matrix from the mechanical coordinate system to the Stewart calculation coordinate system is named T trans_mach_st , and the transformation matrix from the Stewart calculation coordinate system to the mechanical coordinate system is the inverse matrix
Figure PCTCN2020141231-appb-000012
步骤3:根据转换矩阵和转换矩阵的逆矩阵进行机械坐标系与Stewart计算坐标系的坐标转换。Step 3: According to the transformation matrix and the inverse matrix of the transformation matrix, coordinate transformation between the mechanical coordinate system and the Stewart calculation coordinate system is performed.
为了简化探头坐标系与参考坐标系之间的转换关系,在其中的一些实施例中,探头坐标系的X、Y、Z轴的方向与末端关节的关节坐标系的X、Y、Z轴的方向相同,且探头坐标系的Z轴和末端关节的关节坐标系的Z轴共线。在第一串联机械臂的末端关节为并联平台的情况下,探头坐标系的X、Y、Z轴的方向与动平台坐标系的X、Y、Z轴的方向相同,且探头坐标系的Z轴和动平台坐标系的Z轴共线。In order to simplify the conversion relationship between the probe coordinate system and the reference coordinate system, in some of the embodiments, the directions of the X, Y, and Z axes of the probe coordinate system and the directions of the X, Y, and Z axes of the joint coordinate system of the end joint are The directions are the same, and the Z axis of the probe coordinate system and the Z axis of the joint coordinate system of the end joint are collinear. When the end joint of the first serial manipulator is a parallel platform, the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the moving platform coordinate system, and the Z axis of the probe coordinate system is the same. The axis and the Z axis of the moving platform coordinate system are collinear.
下面通过优选实施例对本申请实施例的探头坐标系以及探头坐标系与参考坐标系的转换关系进行描述和说明。The following describes and illustrates the probe coordinate system and the conversion relationship between the probe coordinate system and the reference coordinate system in the embodiments of the present application by using preferred embodiments.
图6是本申请优选实施例的超声探测装置的探头坐标系与Stewart计算坐标系的示意图。在超声探测装置与Stewart动平台之间设有一个旋转电机,电机的旋转角度记为θ m。超声探测装置的末端点位于Stewart动平台坐标系的Z轴上,无论θ m取何值,超声探测装置的末端点相对Stewart动平台的位置都是不变。 FIG. 6 is a schematic diagram of a probe coordinate system and a Stewart calculation coordinate system of an ultrasonic detection device according to a preferred embodiment of the present application. A rotary motor is arranged between the ultrasonic detection device and the Stewart moving platform, and the rotation angle of the motor is recorded as θ m . The end point of the ultrasonic detection device is located on the Z axis of the Stewart moving platform coordinate system. No matter what the value of θ m is, the position of the end point of the ultrasonic detection device relative to the Stewart moving platform remains unchanged.
在超声探测装置末端建立探头坐标系,探头坐标系的原点与探头末端点重合,且忽略Stewart动平台与超声探测装置之间的旋转电机的运动,默认探头坐标系的XYZ轴始终与Stewart动平台坐标系的XYZ轴平行。也就是将Stewart动平台与超声探测装置之间的电机的 旋转量θ m视为超声探测平面相对探头坐标系的运动。 A probe coordinate system is established at the end of the ultrasonic detection device. The origin of the probe coordinate system coincides with the end point of the probe, and the motion of the rotating motor between the Stewart moving platform and the ultrasonic detection device is ignored. By default, the XYZ axis of the probe coordinate system is always the same as the Stewart moving platform. The XYZ axes of the coordinate system are parallel. That is, the rotation amount θ m of the motor between the Stewart moving platform and the ultrasonic detection device is regarded as the movement of the ultrasonic detection plane relative to the probe coordinate system.
需要说明的是,探头坐标系的建立方式并不是唯一的,并不一定要将探头相对于Stewart动平台的转动θ m转化为超声探测平面相对于探头坐标系的转动。 It should be noted that the establishment of the probe coordinate system is not unique, and it is not necessary to convert the rotation θ m of the probe relative to the Stewart moving platform into the rotation of the ultrasonic detection plane relative to the probe coordinate system.
在通过主从映射方法控制Stewart并联平台的情况下,超声探头的运动可以由主操作手控制,例如通过混联主从映射算法将主操作手在竖直方向上的运动量映射到并联平台,将主操作手的水平平面内的运动量映射到串联机械臂。通过混联主从映射算法可得超声探测装置末端点在Stewart计算坐标系的位置C coord_det_st(三维列向量),和探头坐标系相对于Stewart计算坐标系的姿态角
Figure PCTCN2020141231-appb-000013
(其中
Figure PCTCN2020141231-appb-000014
为固定值),即可解算得探头坐标系到Stewart计算坐标系的转换矩阵T trans_det_st
In the case of controlling the Stewart parallel platform through the master-slave mapping method, the motion of the ultrasound probe can be controlled by the master operator. The amount of movement in the horizontal plane of the master manipulator is mapped to the tandem robotic arm. Through the hybrid master-slave mapping algorithm, the position C coord_det_st (three-dimensional column vector) of the end point of the ultrasonic detection device in the coordinate system calculated by Stewart, and the attitude angle of the coordinate system of the probe relative to the coordinate system calculated by Stewart can be obtained
Figure PCTCN2020141231-appb-000013
(in
Figure PCTCN2020141231-appb-000014
is a fixed value), the transformation matrix T trans_det_st from the probe coordinate system to the Stewart calculation coordinate system can be calculated:
Figure PCTCN2020141231-appb-000015
Figure PCTCN2020141231-appb-000015
其中,R coord_det_st为旋转矩阵,用姿态角
Figure PCTCN2020141231-appb-000016
表示为:
Among them, R coord_det_st is the rotation matrix, using the attitude angle
Figure PCTCN2020141231-appb-000016
Expressed as:
Figure PCTCN2020141231-appb-000017
Figure PCTCN2020141231-appb-000017
操作者可以通过主操作手操控末端安装有超声探测装置的机械臂,使穿刺目标靶点位于超声探测平面(声束平面)上。超声探测平面为经过超声探测装置末端点且垂直于探头坐标系XY平面的一个平面,在超声探测平面内的任意一点可测得其相对于探头坐标系的距离l c和与Z轴的夹角θ cThe operator can control the mechanical arm with the ultrasonic detection device installed at the end through the main operator, so that the puncture target is located on the ultrasonic detection plane (sound beam plane). The ultrasonic detection plane is a plane that passes through the end point of the ultrasonic detection device and is perpendicular to the XY plane of the probe coordinate system. At any point in the ultrasonic detection plane, the distance l c relative to the probe coordinate system and the included angle with the Z axis can be measured. θ c .
通过l c,θ c,θ m这三个已知量,可知穿刺目标靶点在探头坐标系下的坐标Sst(Xst,Yst,Zst),记为C coord_tar_detThrough the three known quantities of l c , θ c , and θ m , the coordinates Sst (Xst, Yst, Zst) of the puncture target in the probe coordinate system can be known, which is denoted as C coord_tar_det :
C coord_tar_det=[l csinθ ccosθ m l csinθ csinθ m l ccosθ c 1] TC coord_tar_det =[l c sinθ c cosθ m l c sinθ c sinθ m l c cosθ c 1] T .
通过坐标转换可得目标靶点在机械坐标系下的坐标C coord_tar_machThe coordinate C coord_tar_mach of the target point in the mechanical coordinate system can be obtained by coordinate transformation:
C coord_tar_mach=T trans_mach_stT trans_st_detC coord_tar_detC coord_tar_mach =T trans_mach_st T trans_st_det C coord_tar_det .
在其中的一些实施例中,超声定位穿刺系统还包括第二串联机械臂和搭载在第二串联机械臂末端的穿刺器械;计算机设备还分别与第二串联机械臂电性连接,处理器在运行计算机程序时还执行下列步骤:将靶点在参考坐标系内的第二位置信息作为目标穿刺点,控制第二串联机械臂的穿刺器械对目标穿刺点进行穿刺。In some of the embodiments, the ultrasonic positioning and puncturing system further includes a second serial robotic arm and a puncturing instrument mounted on the end of the second serial robotic arm; the computer equipment is also electrically connected to the second serial robotic arm, and the processor is running The computer program also executes the following steps: taking the second position information of the target in the reference coordinate system as the target puncture point, and controlling the puncture instrument of the second serial robotic arm to puncture the target puncture point.
通过坐标转换可得目标靶点在穿刺机械臂Stewart计算坐标系下的坐标C coord_tar_stlThe coordinate C coord_tar_stl of the target point in the coordinate system calculated by Stewart of the puncture manipulator can be obtained by coordinate transformation:
C coord_tar_stl=T trans_stl_machC coord_tar_machC coord_tar_stl =T trans_stl_mach C coord_tar_mach .
其中,T trans_stl_mach为穿刺机械臂Stewart计算坐标系到计算坐标系的转换矩阵。 Among them, T trans_stl_mach is the transformation matrix from the calculated coordinate system of the puncture manipulator Stewart to the calculated coordinate system.
已知目标靶点在穿刺机械臂Stewart计算坐标系下的坐标以及预设的穿刺路径和穿刺速度,即可通过Stewart并联平台的逆运动学解算穿刺机械臂Stewart平台的关节运动量,用穿刺机械臂进行精准的穿刺。Knowing the coordinates of the target point in the coordinate system calculated by the puncture manipulator Stewart, as well as the preset puncture path and puncture speed, the joint motion of the Stewart platform of the puncture manipulator can be calculated through the inverse kinematics of the Stewart parallel platform. Arm for precise puncture.
通过上述具有两个机械臂的系统,第一串联机械臂的末端搭载超声探测装置,第二串联机械臂的末端搭载穿刺器械,穿刺消融手术的执行由两个机械臂协作完成,第一串联机械臂持超声探头,在人体表面扫描寻找病灶,并定位穿刺目标靶点,可以实现靶点的精准定位。第二串联机械臂持消融针,对目标靶点进行穿刺消融,可以实现对靶点的精准穿刺。整个定位和穿刺的过程几乎对人为操作的依赖性较小,更加安全,并且上述靶点定位和穿刺方法兼容现有的具有多机械臂且机械臂末端设置有并联平台的机器人,不需要独立开发复杂的系统。Through the above system with two robotic arms, the end of the first serial robotic arm is equipped with an ultrasonic detection device, and the end of the second serial robotic arm is equipped with a puncture instrument, the execution of the puncture and ablation operation is completed by the cooperation of the two robotic arms, and the first serial The arm holds the ultrasound probe, scans the surface of the human body to find the lesions, and locates the puncture target, which can achieve precise positioning of the target. The second serial robotic arm holds the ablation needle and performs puncture and ablation on the target point, which can achieve precise puncture of the target point. The whole process of positioning and puncturing is almost less dependent on human operation and is safer, and the above-mentioned target positioning and puncturing methods are compatible with existing robots with multiple manipulators and a parallel platform at the end of the manipulator, and no independent development is required. complex system.
另外,结合上述实施例中的方法,本申请实施例还提供一种计算机可读存储介质来实现。该计算机可读存储介质上存储有计算机程序;计算机程序被设置为运行时执行下列步骤。In addition, in combination with the methods in the foregoing embodiments, the embodiments of the present application further provide a computer-readable storage medium for implementation. A computer program is stored on the computer-readable storage medium; the computer program is configured to perform the following steps when executed.
S1:在超声探测装置投射的超声探测平面内定位靶点的第一位置信息。S1: The first position information of the target point is located in the ultrasonic detection plane projected by the ultrasonic detection device.
S2:获取超声探测平面的第一姿态信息以及第一串联机械臂的第二姿态信息。S2: Acquire the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator.
S3:根据第一姿态信息和第二姿态信息,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系。S3: According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located.
S4:根据位置转换关系和第一位置信息,确定靶点在参考坐标系内的第二位置信息。S4: Determine the second position information of the target point in the reference coordinate system according to the position conversion relationship and the first position information.
在其中一些实施例中,在超声探测装置投射的超声探测平面内定位靶点的第一位置信息包括:获取由超声探测装置采集的超声图像,在超声图像中定位靶点的位置信息;根据靶点在超声图像中的位置信息以及超声探测装置的标定参数,确定靶点在超声探测装置投射的 超声探测平面内的第一位置信息。In some of these embodiments, locating the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device includes: acquiring an ultrasonic image collected by the ultrasonic detection device, and locating the position information of the target point in the ultrasonic image; The position information of the point in the ultrasonic image and the calibration parameters of the ultrasonic detection device determine the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device.
在其中一些实施例中,第一位置信息包括靶点相对于末端端点的距离,以及靶点和末端端点的连接线与旋转轴的夹角。In some of the embodiments, the first position information includes the distance of the target point relative to the end point, and the angle between the connecting line between the target point and the end point and the rotation axis.
在其中一些实施例中,根据第一姿态信息和第二姿态信息,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系包括:在第一串联机械臂处于初始状态时,建立第一串联机械臂的参考坐标系和第一串联机械臂的各关节的关节坐标系,以及在超声探测装置处于初始状态时,以超声探测装置的末端端点为坐标原点建立探头坐标系;获取第一串联机械臂的参考坐标系和第一串联机械臂的各关节的关节坐标系的DH参数,并根据DH参数,确定第一串联机械臂的参考坐标系到第一串联机械臂的末端关节的关节坐标系之间的第一位置转换关系;获取探头坐标系的坐标原点在末端关节的关节坐标系内的位置信息,以及探头坐标系相对于末端关节的关节坐标系的姿态角,并根据该位置信息、姿态角以及第一位置转换关系,确定参考坐标系与探头坐标系之间的第二位置转换关系;获取超声探测平面的旋转角度,并根据旋转角度确定超声探测平面所在坐标系与探头坐标系之间的第三位置转换关系;根据第二位置转换关系和第三位置转换关系,确定超声定位穿刺系统的参考坐标系与超声探测平面所在坐标系的位置转换关系。In some of the embodiments, according to the first attitude information and the second attitude information, determining the position conversion relationship between the reference coordinate system of the ultrasonic positioning and puncturing system and the coordinate system where the ultrasonic detection plane is located includes: when the first serial robotic arm is in an initial state , establish the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator, and when the ultrasonic detection device is in the initial state, establish the probe coordinate system with the end point of the ultrasonic detection device as the coordinate origin; Obtain the reference coordinate system of the first serial manipulator and the DH parameters of the joint coordinate system of each joint of the first serial manipulator, and determine the reference coordinate system of the first serial manipulator to the end of the first serial manipulator according to the DH parameters The first position conversion relationship between the joint coordinate systems of the joints; obtain the position information of the coordinate origin of the probe coordinate system in the joint coordinate system of the end joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the end joint, and According to the position information, the attitude angle and the first position conversion relationship, determine the second position conversion relationship between the reference coordinate system and the probe coordinate system; obtain the rotation angle of the ultrasonic detection plane, and determine the coordinate system where the ultrasonic detection plane is located according to the rotation angle The third position conversion relationship with the probe coordinate system; according to the second position conversion relationship and the third position conversion relationship, the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located is determined.
在其中一些实施例中,第一串联机械臂包括旋转关节和移动关节,旋转关节的关节坐标系的Z轴沿旋转轴设置,移动关节的关节坐标系的Z轴沿移动方向设置,参考坐标系与各关节的关节坐标系同为左手系或右手系,且当旋转关节与在前的移动关节连接时,旋转关节的关节坐标系的原点与该在前的移动关节的关节坐标系的原点重合。In some of these embodiments, the first serial manipulator includes a rotating joint and a moving joint, the Z axis of the joint coordinate system of the rotating joint is set along the rotating axis, the Z axis of the joint coordinate system of the moving joint is set along the moving direction, and the reference coordinate system The joint coordinate system of each joint is the same as the left-handed or right-handed system, and when the rotary joint is connected to the previous mobile joint, the origin of the joint coordinate system of the rotary joint coincides with the origin of the joint coordinate system of the previous mobile joint .
在其中一些实施例中,探头坐标系的X、Y、Z轴的方向与末端关节的关节坐标系的X、Y、Z轴的方向相同,且探头坐标系的Z轴和末端关节的关节坐标系的Z轴共线。In some of these embodiments, the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the joint coordinate system of the end joint, and the Z axis of the probe coordinate system and the joint coordinates of the end joint The Z axis of the system is collinear.
在其中一些实施例中,计算机程序还被设置为运行时执行下列步骤:将靶点在参考坐标系内的第二位置信息作为目标穿刺点,控制第二串联机械臂的穿刺器械对目标穿刺点进行穿刺。In some of these embodiments, the computer program is further configured to execute the following steps when running: using the second position information of the target point in the reference coordinate system as the target puncture point, and controlling the puncture instrument of the second serial robotic arm to puncture the target puncture point Perform puncture.
需要说明的是,本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例,本实施例在此不再赘述。It should be noted that, for specific examples in this embodiment, reference may be made to the examples described in the foregoing embodiments and optional implementation manners, and details are not described herein again in this embodiment.
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above-described embodiments can be combined arbitrarily. For the sake of brevity, all possible combinations of the technical features in the above-described embodiments are not described. However, as long as there is no contradiction between the combinations of these technical features, All should be regarded as the scope described in this specification.
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对申请专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在 不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only represent several embodiments of the present application, and the descriptions thereof are relatively specific and detailed, but should not be construed as a limitation on the scope of the patent application. It should be pointed out that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can also be made, which all belong to the protection scope of the present application. Therefore, the scope of protection of the patent of the present application shall be subject to the appended claims.

Claims (11)

  1. 一种超声定位穿刺系统,其特征在于,所述超声定位穿刺系统包括计算机设备、第一串联机械臂和搭载在所述第一串联机械臂末端的超声探测装置,所述计算机设备分别与所述第一串联机械臂和所述超声探测装置电性连接,所述计算机设备包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器在运行所述计算机程序时执行下列步骤:An ultrasonic positioning and puncturing system, characterized in that, the ultrasonic positioning and puncturing system comprises computer equipment, a first serial manipulator arm and an ultrasonic detection device mounted on the end of the first serial manipulator arm, and the computer equipment is respectively connected with the The first serial manipulator is electrically connected to the ultrasonic detection device, and the computer device includes a memory and a processor, where a computer program is stored in the memory, and the processor executes the following steps when running the computer program:
    在所述超声探测装置投射的超声探测平面内定位靶点的第一位置信息;Locating the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device;
    获取所述超声探测平面的第一姿态信息以及所述第一串联机械臂的第二姿态信息;acquiring the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator;
    根据所述第一姿态信息和所述第二姿态信息,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系;According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located;
    根据所述位置转换关系和所述第一位置信息,确定所述靶点在所述参考坐标系内的第二位置信息。According to the position conversion relationship and the first position information, second position information of the target point in the reference coordinate system is determined.
  2. 根据权利要求1所述的超声定位穿刺系统,其特征在于,所述在所述超声探测装置投射的超声探测平面内定位靶点的第一位置信息包括:The ultrasonic positioning and puncturing system according to claim 1, wherein the first position information for locating the target point in the ultrasonic detection plane projected by the ultrasonic detection device comprises:
    获取由所述超声探测装置采集的超声图像,在所述超声图像中定位所述靶点的位置信息;acquiring an ultrasonic image collected by the ultrasonic detection device, and locating the position information of the target point in the ultrasonic image;
    根据所述靶点在所述超声图像中的位置信息以及所述超声探测装置的标定参数,确定所述靶点在所述超声探测装置投射的超声探测平面内的所述第一位置信息。According to the position information of the target point in the ultrasonic image and the calibration parameters of the ultrasonic detection device, the first position information of the target point in the ultrasonic detection plane projected by the ultrasonic detection device is determined.
  3. 根据权利要求1所述的超声定位穿刺系统,其特征在于,所述超声探测装置与所述第一串联机械臂的末端可转动地连接,所述超声探测平面形成于所述超声探测装置的末端端点并自所述超声探测装置的末端端点沿扇面延伸;所述超声探测平面跟随所述超声探测装置旋转,所述超声探测平面的旋转轴经过所述末端端点;所述第一姿态信息包括所述超声探测平面的旋转角度。The ultrasonic positioning and puncturing system according to claim 1, wherein the ultrasonic detection device is rotatably connected to the end of the first serial mechanical arm, and the ultrasonic detection plane is formed at the end of the ultrasonic detection device The end point extends along the sector from the end end point of the ultrasonic detection device; the ultrasonic detection plane rotates following the ultrasonic detection device, and the rotation axis of the ultrasonic detection plane passes through the end end point; the first attitude information includes all the The rotation angle of the ultrasonic detection plane.
  4. 根据权利要求3所述的超声定位穿刺系统,其特征在于,所述第一位置信息包括所述靶点相对于所述末端端点的距离,以及所述靶点和所述末端端点的连接线与所述旋转轴的夹角。The ultrasonic positioning puncture system according to claim 3, wherein the first position information includes the distance of the target point relative to the end point, and the connecting line between the target point and the end point and the The included angle of the rotation axis.
  5. 根据权利要求1所述的超声定位穿刺系统,其特征在于,所述第一串联机械臂包括多个关节,多个所述关节包括:旋转关节和/或移动关节,所述第二姿态信息包括每个所述关节的运动量信息。The ultrasonic positioning and puncturing system according to claim 1, wherein the first serial robotic arm includes a plurality of joints, the plurality of joints include: a rotation joint and/or a moving joint, and the second posture information includes Movement amount information for each of the joints.
  6. 根据权利要求3所述的超声定位穿刺系统,其特征在于,所述根据所述第一姿态信息和所述第二姿态信息,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系包括:The ultrasonic positioning and puncturing system according to claim 3, wherein the reference coordinate system of the ultrasonic positioning and puncturing system and the ultrasonic detection plane are determined according to the first attitude information and the second attitude information. The position conversion relationship of the coordinate system where it is located includes:
    在所述第一串联机械臂处于初始状态时,建立所述第一串联机械臂的参考坐标系和所述 第一串联机械臂的各关节的关节坐标系,以及在所述超声探测装置处于初始状态时,以所述超声探测装置的末端端点为坐标原点建立探头坐标系;When the first serial manipulator is in the initial state, the reference coordinate system of the first serial manipulator and the joint coordinate system of each joint of the first serial manipulator are established, and when the ultrasonic detection device is in the initial state In the state, the probe coordinate system is established with the end point of the ultrasonic detection device as the coordinate origin;
    获取所述第一串联机械臂的参考坐标系和所述第一串联机械臂的各关节的关节坐标系的DH参数,并根据所述DH参数,确定所述第一串联机械臂的参考坐标系到所述第一串联机械臂的末端关节的关节坐标系之间的第一位置转换关系;Obtain the reference coordinate system of the first serial manipulator and the DH parameters of the joint coordinate system of each joint of the first serial manipulator, and determine the reference coordinate system of the first serial manipulator according to the DH parameters a first position conversion relationship to the joint coordinate system of the end joint of the first serial manipulator;
    获取所述探头坐标系的坐标原点在所述末端关节的关节坐标系内的位置信息,以及所述探头坐标系相对于所述末端关节的关节坐标系的姿态角,并根据该位置信息、姿态角以及所述第一位置转换关系,确定所述参考坐标系与所述探头坐标系之间的第二位置转换关系;Obtain the position information of the coordinate origin of the probe coordinate system in the joint coordinate system of the terminal joint, and the attitude angle of the probe coordinate system relative to the joint coordinate system of the terminal joint, and according to the position information, attitude angle and the first position conversion relationship, to determine the second position conversion relationship between the reference coordinate system and the probe coordinate system;
    获取所述超声探测平面的旋转角度,并根据所述旋转角度确定所述超声探测平面所在坐标系与所述探头坐标系之间的第三位置转换关系;acquiring the rotation angle of the ultrasonic detection plane, and determining a third position conversion relationship between the coordinate system where the ultrasonic detection plane is located and the probe coordinate system according to the rotation angle;
    根据所述第二位置转换关系和所述第三位置转换关系,确定所述超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系。According to the second position conversion relationship and the third position conversion relationship, the position conversion relationship between the reference coordinate system of the ultrasonic positioning and puncturing system and the coordinate system where the ultrasonic detection plane is located is determined.
  7. 根据权利要求6所述的超声定位穿刺系统,其特征在于,所述第一串联机械臂包括旋转关节和移动关节,所述旋转关节的关节坐标系的Z轴沿旋转轴设置,所述移动关节的关节坐标系的Z轴沿移动方向设置,所述参考坐标系与各关节的关节坐标系同为左手系或右手系,且当旋转关节与在前的移动关节连接时,旋转关节的关节坐标系的原点与该在前的移动关节的关节坐标系的原点重合。The ultrasonic positioning and puncturing system according to claim 6, wherein the first serial mechanical arm comprises a rotation joint and a moving joint, the Z axis of the joint coordinate system of the rotating joint is set along the rotation axis, and the moving joint The Z-axis of the joint coordinate system is set along the moving direction, the reference coordinate system and the joint coordinate system of each joint are both left-handed or right-handed systems, and when the rotating joint is connected to the previous moving joint, the joint coordinates of the rotating joint The origin of the system coincides with the origin of the joint coordinate system of the preceding moving joint.
  8. 根据权利要求6所述的超声定位穿刺系统,其特征在于,所述探头坐标系的X、Y、Z轴的方向与所述末端关节的关节坐标系的X、Y、Z轴的方向相同,且所述探头坐标系的Z轴和所述末端关节的关节坐标系的Z轴共线。The ultrasonic positioning puncture system according to claim 6, wherein the directions of the X, Y, and Z axes of the probe coordinate system are the same as the directions of the X, Y, and Z axes of the joint coordinate system of the distal joint, And the Z axis of the probe coordinate system and the Z axis of the joint coordinate system of the end joint are collinear.
  9. 根据权利要求1所述的超声定位穿刺系统,其特征在于,所述第一串联机械臂的末端关节为并联运动平台,所述并联运动平台包括静平台、动平台和多个设置在所述静平台与所述动平台之间的伸缩组件,所述超声探测装置设于所述动平台。The ultrasonic positioning and puncturing system according to claim 1, wherein the end joint of the first serial manipulator is a parallel motion platform, and the parallel motion platform comprises a static platform, a moving platform and a plurality of A telescopic assembly between a platform and the moving platform, and the ultrasonic detection device is arranged on the moving platform.
  10. 根据权利要求1所述的超声定位穿刺系统,其特征在于,所述超声定位穿刺系统还包括第二串联机械臂和搭载在所述第二串联机械臂末端的穿刺器械;所述计算机设备还分别与所述第二串联机械臂电性连接,所述处理器在运行所述计算机程序时还执行下列步骤:The ultrasonic positioning and puncturing system according to claim 1, wherein the ultrasonic positioning and puncturing system further comprises a second serial robotic arm and a puncturing instrument carried at the end of the second serial mechanical arm; the computer equipment further comprises Electrically connected with the second serial manipulator, the processor also executes the following steps when running the computer program:
    将所述靶点在所述参考坐标系内的第二位置信息作为目标穿刺点,控制所述第二串联机械臂的穿刺器械对所述目标穿刺点进行穿刺。The second position information of the target point in the reference coordinate system is used as the target puncture point, and the puncture instrument of the second serial robotic arm is controlled to puncture the target puncture point.
  11. 一种存储介质,所述存储介质中存储有计算机程序,其特征在于,所述计算机程序被设置为运行时执行下列步骤:A storage medium, wherein a computer program is stored in the storage medium, wherein the computer program is configured to execute the following steps when running:
    在超声探测装置投射的超声探测平面内定位靶点的第一位置信息;Locating the first position information of the target in the ultrasonic detection plane projected by the ultrasonic detection device;
    获取所述超声探测平面的第一姿态信息以及第一串联机械臂的第二姿态信息;acquiring the first attitude information of the ultrasonic detection plane and the second attitude information of the first serial manipulator;
    根据所述第一姿态信息和所述第二姿态信息,确定超声定位穿刺系统的参考坐标系与所述超声探测平面所在坐标系的位置转换关系;According to the first attitude information and the second attitude information, determine the position conversion relationship between the reference coordinate system of the ultrasonic positioning puncture system and the coordinate system where the ultrasonic detection plane is located;
    根据所述位置转换关系和所述第一位置信息,确定所述靶点在所述参考坐标系内的第二位置信息。According to the position conversion relationship and the first position information, second position information of the target point in the reference coordinate system is determined.
PCT/CN2020/141231 2020-12-30 2020-12-30 Ultrasonic positioning puncture system and storage medium WO2022141153A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/141231 WO2022141153A1 (en) 2020-12-30 2020-12-30 Ultrasonic positioning puncture system and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/141231 WO2022141153A1 (en) 2020-12-30 2020-12-30 Ultrasonic positioning puncture system and storage medium

Publications (1)

Publication Number Publication Date
WO2022141153A1 true WO2022141153A1 (en) 2022-07-07

Family

ID=82258797

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/141231 WO2022141153A1 (en) 2020-12-30 2020-12-30 Ultrasonic positioning puncture system and storage medium

Country Status (1)

Country Link
WO (1) WO2022141153A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115005991A (en) * 2022-08-03 2022-09-06 北京壹点灵动科技有限公司 Precision detection method of surgical navigation device and surgical navigation precision detection device
CN115252131A (en) * 2022-08-12 2022-11-01 上海精劢医疗科技有限公司 Pose information correction method for surgical tool
CN115556109A (en) * 2022-10-24 2023-01-03 深圳市通用测试系统有限公司 Method and device for positioning mechanical arm in test system
CN116077155A (en) * 2023-04-06 2023-05-09 深圳惟德精准医疗科技有限公司 Puncture method based on optical tracking equipment and mechanical arm and related device
CN117481682A (en) * 2023-09-12 2024-02-02 中山大学 Remote control method, device and system for color ultrasonic probe, electronic equipment and medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030187345A1 (en) * 2002-03-27 2003-10-02 Aloka Co., Ltd. Ultrasonic medical system
US20170055940A1 (en) * 2014-04-28 2017-03-02 Mazor Robotics Ltd. Ultrasound guided hand held robot
CN109701169A (en) * 2018-12-27 2019-05-03 成植温 A kind of system of tumor treatment that mechanical arm punctures
CN110200683A (en) * 2019-07-03 2019-09-06 北京积水潭医院 Ultrasound-guided Biopsy attitude monitoring method and device
CN111227943A (en) * 2020-01-23 2020-06-05 诺创智能医疗科技(杭州)有限公司 Control method of surgical mechanical arm, computer equipment and surgical mechanical arm
CN111437011A (en) * 2020-03-30 2020-07-24 中国科学院深圳先进技术研究院 Puncture surgery robot system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030187345A1 (en) * 2002-03-27 2003-10-02 Aloka Co., Ltd. Ultrasonic medical system
US20170055940A1 (en) * 2014-04-28 2017-03-02 Mazor Robotics Ltd. Ultrasound guided hand held robot
CN109701169A (en) * 2018-12-27 2019-05-03 成植温 A kind of system of tumor treatment that mechanical arm punctures
CN110200683A (en) * 2019-07-03 2019-09-06 北京积水潭医院 Ultrasound-guided Biopsy attitude monitoring method and device
CN111227943A (en) * 2020-01-23 2020-06-05 诺创智能医疗科技(杭州)有限公司 Control method of surgical mechanical arm, computer equipment and surgical mechanical arm
CN111437011A (en) * 2020-03-30 2020-07-24 中国科学院深圳先进技术研究院 Puncture surgery robot system

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115005991A (en) * 2022-08-03 2022-09-06 北京壹点灵动科技有限公司 Precision detection method of surgical navigation device and surgical navigation precision detection device
CN115252131A (en) * 2022-08-12 2022-11-01 上海精劢医疗科技有限公司 Pose information correction method for surgical tool
CN115556109A (en) * 2022-10-24 2023-01-03 深圳市通用测试系统有限公司 Method and device for positioning mechanical arm in test system
CN115556109B (en) * 2022-10-24 2024-06-11 深圳市通用测试系统有限公司 Positioning method and device for mechanical arm in test system
CN116077155A (en) * 2023-04-06 2023-05-09 深圳惟德精准医疗科技有限公司 Puncture method based on optical tracking equipment and mechanical arm and related device
CN116077155B (en) * 2023-04-06 2023-06-27 深圳惟德精准医疗科技有限公司 Surgical navigation method based on optical tracking equipment and mechanical arm and related device
CN117481682A (en) * 2023-09-12 2024-02-02 中山大学 Remote control method, device and system for color ultrasonic probe, electronic equipment and medium
CN117481682B (en) * 2023-09-12 2024-04-09 中山大学 Remote control method, device and system for color ultrasonic probe, electronic equipment and medium

Similar Documents

Publication Publication Date Title
WO2022141153A1 (en) Ultrasonic positioning puncture system and storage medium
CN112754616B (en) Ultrasonic positioning puncture system and storage medium
CN114041875B (en) Integrated operation positioning navigation system
CN113133832B (en) Calibration method and system for double-arm robot puncture system
WO2022218389A1 (en) Ultrasonic and x-ray combined execution operation method, device and system, and computer-readable storage medium
CN113288429A (en) Space registration and real-time navigation method of breast minimally invasive interventional operation robot
JP7383489B2 (en) Integration of robotic device guidance and acoustic probes
CN113940755A (en) Surgical operation planning and navigation method integrating operation and image
CN113180828B (en) Surgical robot constraint motion control method based on rotation theory
CN113400325A (en) Navigation and positioning method for dental implant robot
WO2022141138A1 (en) Hybrid master-slave mapping method, robotic arm system, and computer device
CN111603205A (en) Three-dimensional image reconstruction and positioning analysis system used in CT (computed tomography) cabin of puncture surgical robot
CN113974835A (en) Operation robot motion control method based on telecentric motionless point constraint
CN115500940A (en) Positioning display method of surgical needle and related device
CN113456220B (en) Alignment method, surgical robot, and computer storage medium
CN116269812A (en) Master-slave operation puncture system and planning method
WO2022198615A1 (en) Calibration method and system for dual-arm robot puncture system
Onogi et al. Robotic ultrasound guidance by B-scan plane positioning control
CN112107366B (en) Mixed reality ultrasonic navigation system
Azimi et al. Teleoperative control of intraocular robotic snake: Vision-based angular calibration
WO2021012142A1 (en) Surgical robot system and control method therefor
CN117084790B (en) Puncture azimuth control method and device, computer equipment and storage medium
Tseng et al. Image‐guided robotic navigation system for neurosurgery
CN116919596B (en) Instrument navigation method, system, device, equipment and storage medium
CN117084791B (en) Puncture azimuth resolving method and puncture operation executing system

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: 20967498

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20967498

Country of ref document: EP

Kind code of ref document: A1