WO2013183276A1 - カテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラム - Google Patents

カテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラム Download PDF

Info

Publication number
WO2013183276A1
WO2013183276A1 PCT/JP2013/003485 JP2013003485W WO2013183276A1 WO 2013183276 A1 WO2013183276 A1 WO 2013183276A1 JP 2013003485 W JP2013003485 W JP 2013003485W WO 2013183276 A1 WO2013183276 A1 WO 2013183276A1
Authority
WO
WIPO (PCT)
Prior art keywords
catheter
rotation angle
unit
tip
hole
Prior art date
Application number
PCT/JP2013/003485
Other languages
English (en)
French (fr)
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 JP2014506652A priority Critical patent/JP5526297B1/ja
Publication of WO2013183276A1 publication Critical patent/WO2013183276A1/ja
Priority to US14/197,802 priority patent/US9918682B2/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/12Arrangements for detecting or locating foreign bodies
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • GPHYSICS
    • G09EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09BEDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/28Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
    • G09B23/285Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine for injections, endoscopy, bronchoscopy, sigmoidscopy, insertion of contraceptive devices or enemas
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M25/00Catheters; Hollow probes
    • A61M25/01Introducing, guiding, advancing, emplacing or holding catheters
    • A61M25/0105Steering means as part of the catheter or advancing means; Markers for positioning

Definitions

  • the present invention relates to a catheter tip rotation angle measurement device, a catheter rotation angle measurement method, and a catheter rotation angle measurement program for measuring the rotation angle of a catheter tip.
  • a catheter has been inserted into a patient's blood vessel, etc. to diagnose or treat a lesion or stenosis (for example, PCI (Percutaneous Coronary Intervention) or TAE (Transcatheter Arterial Embolization)). Etc.) are being carried out.
  • PCI Percutaneous Coronary Intervention
  • TAE Transcatheter Arterial Embolization
  • the operator operates the catheter while monitoring the fluoroscopic image of the tip of the catheter during the above diagnosis or treatment. At that time, in order to adjust the direction of the distal end of the catheter in the branching direction of the target blood vessel, for example, in the branching of the blood vessel, the operator needs to correctly grasp the direction of the distal end of the catheter.
  • Patent Document 1 As a conventional technique for detecting the orientation of a catheter, a method for detecting the rotation angle of a catheter at the hand of an operator is disclosed (see Patent Document 1). In addition, a method is disclosed in which a rotation detection sensor is attached to the distal end portion of an electronic endoscope, and a display image on a monitor caused by twisting of the electronic endoscope is automatically corrected (see Patent Document 2).
  • JP 2009-162920 A JP-A-7-246183 Special table 2009-522016
  • Patent Document 1 cannot correctly grasp the rotation angle of the catheter tip when the rotation angle at the operator's hand does not match the rotation angle of the catheter tip.
  • Patent Document 2 it is necessary to attach a rotation detection sensor to the distal end of the catheter.
  • an object of the present invention is to provide a catheter tip rotation angle measurement device, a catheter rotation angle measurement method, and a catheter rotation angle measurement program for solving the above-described problems and measuring the rotation angle of the catheter tip more accurately.
  • a catheter tip rotation angle measuring device is a device that measures the rotation angle of the tip of a catheter inserted into a living body tube, and is a second through a first hole.
  • a movement restricting portion that restricts movement of the catheter in the insertion direction by passing the catheter from the first hole to the second hole;
  • a first rotation angle measurement unit that measures a first rotation angle that is the rotation angle of the catheter in the hole, and a second rotation angle measurement unit that measures a second rotation angle that is the rotation angle of the catheter in the second hole.
  • a catheter insertion length measuring unit that measures a catheter insertion length that is a length from the first hole or the second hole to the distal end portion of the catheter, and a first measured by the first rotation angle measuring unit.
  • One rotation angle and the second Using the angle difference from the second rotation angle measured by the turning angle measurement unit, the catheter insertion length measured by the catheter insertion length measurement unit, and the distance between the first hole and the second hole,
  • a catheter tip rotation angle measuring device comprising: a catheter tip rotation angle calculation unit that calculates a rotation angle of the catheter tip.
  • the rotation angle of the distal end of the catheter can be measured more accurately.
  • FIG. 1 is a block diagram showing a functional configuration of a catheter tip rotation angle measuring device in the first embodiment
  • FIG. 2 is a diagram illustrating a configuration of the movement restriction unit in the first embodiment
  • FIG. 3 is a diagram illustrating an example of the first rotation angle measurement unit and the second rotation angle measurement unit in the first embodiment.
  • FIG. 4 is a diagram illustrating regions measured by the first rotation angle measurement unit and the second rotation angle measurement unit in the first embodiment.
  • FIG. 5 is a diagram showing the correspondence between line sensors and pixels.
  • FIG. 6 is an explanatory diagram showing the relationship between the cross section of the catheter, the imaging unit, and the pixels, FIG.
  • FIG. 7 is a cross-sectional view showing the relationship between the rotary encoder and the catheter.
  • FIG. 8 is an explanatory diagram showing the configuration of the catheter insertion length measurement unit in the first embodiment
  • FIG. 9 is an overall view of a catheter having a curved portion
  • FIG. 10 is a diagram illustrating a relationship between the first rotation angle, the second rotation angle, the catheter insertion length, the length of the movement restricting portion, and the rotation angle of the catheter tip in the first embodiment
  • FIG. 11A is a diagram illustrating a state in which an example of the rotation angle of the distal end of the catheter presented to the user in the first embodiment is displayed on a display;
  • FIG. 11B is a diagram showing a state in which an example of the rotation angle of the distal end of the catheter presented to the user in the first embodiment is displayed on the display together with the display “correcting”;
  • FIG. 12 is a flowchart showing an example of the processing operation of the surgeon's catheterization in the first embodiment.
  • FIG. 13 is a flowchart illustrating an example of a processing operation of the catheter tip rotation angle measuring device according to the first embodiment.
  • FIG. 14 is a block diagram showing a functional configuration of the catheter tip rotation angle measuring device in the second embodiment,
  • FIG. 15 is a schematic diagram illustrating a functional configuration of the X-ray imaging apparatus according to the second embodiment.
  • FIG. 16 is a diagram illustrating an example of an X-ray image captured by the X-ray imaging apparatus according to the second embodiment.
  • FIG. 17 is a diagram illustrating an example of the distal end of the catheter extracted from the catheter distal end extraction unit in the second embodiment.
  • FIG. 18 is an explanatory diagram illustrating an example of a processing operation of the operation detection unit according to the second embodiment.
  • FIG. 19 is a flowchart showing an example of the processing operation of the catheter tip rotation angle measuring device according to the second embodiment.
  • FIG. 20 is a block diagram showing a functional configuration of the catheter tip rotation angle measuring device in the third embodiment
  • FIG. 21 is an explanatory diagram showing an example of a method for determining the rotation axis of the catheter in the third embodiment, FIG.
  • FIG. 22 is a diagram illustrating an example of a three-dimensional shape of the catheter restored from the catheter rotation axis determination unit in the third embodiment.
  • FIG. 23 is a block diagram showing a functional configuration of the catheter tip rotation angle correction unit in the third embodiment
  • FIG. 24 is a diagram showing a table showing the relationship between the calibration image and the rotation angle of the catheter tip in the third embodiment
  • FIG. 25 is an explanatory diagram showing a comparison between an X-ray image and a calibration image in the third embodiment.
  • FIG. 26 is an explanatory diagram showing a method of calculating an epipolar line when determining the rotation axis of the catheter in the third embodiment
  • FIG. 27 is a flowchart showing an example of the processing operation of the catheter tip rotation angle measuring device according to the third embodiment.
  • FIG. 23 is a block diagram showing a functional configuration of the catheter tip rotation angle correction unit in the third embodiment
  • FIG. 24 is a diagram showing a table showing the relationship between the calibration image and the rotation angle of the catheter tip in the third
  • FIG. 28 is a flowchart showing an example of the processing operation of the catheter tip rotation angle correction unit in the third embodiment.
  • FIG. 29 is an explanatory view showing an example of a guide wire inserted prior to the catheter in the modification of the first embodiment
  • FIG. 30 is a block diagram when automatic determination of stop of the tip rotation angle calculation is performed in the catheter tip rotation angle measuring device in the first modification example
  • FIG. 31A is a diagram illustrating an example of an X-ray imaging image of a distal state of a catheter (or guide wire);
  • FIG. 31B is a diagram showing an example of an X-ray imaging image of the distal end state of the catheter (or guide wire);
  • FIG. 31C is a diagram illustrating an example of an X-ray imaging image of the distal end state of the catheter (or guide wire);
  • FIG. 32 is a flowchart of the operation of the catheter tip rotation angle measuring device in the first modification
  • FIG. 33 is a block diagram of the catheter tip rotation angle measuring device in the second modification
  • FIG. 34 is a block diagram of the catheter tip rotation angle measuring device in the second modification
  • FIG. 35 is a flowchart of the operation of the catheter tip rotation angle measuring device in the second modification
  • FIG. 36 is a block diagram of a catheter tip rotation angle measuring device in a third modification.
  • the operator performs an operation of adjusting the tip of the catheter or the tip of the guide wire to the branch direction of the target blood vessel or the like.
  • the operator needs to accurately grasp the orientation of the distal end of the catheter or the like.
  • the X-ray fluoroscopic image becomes a two-dimensional image, so it is difficult to grasp the direction of the distal end of the catheter or the like, particularly the depth direction.
  • drawbacks That is, it is difficult for the operator to know whether the tip of the catheter or the like is in the front direction or the rear direction of the monitor, which may hinder the smooth operation of the catheter or the like.
  • a method for X-ray imaging of a catheter or the like from two directions using an X-ray fluoroscopy device called a biplane type having two sets of an X-ray generator and an X-ray detector.
  • X-ray imaging by the biplane type is usually performed with an angle difference of 90 degrees, while observing two X-ray fluoroscopic images with an angle difference of 90 degrees alternately, It is a great burden on the operator to instantly and intuitively grasp the shape of the three-dimensional catheter tip while taking into account the imaging angles of the X-ray generator and the X-ray detector.
  • Patent Document 1 and Patent Document 2 have been proposed as conventional methods for detecting the direction of the distal end of the catheter.
  • Patent Document 1 discloses a method of measuring the rotation angle of a catheter at the operator's hand using a rotary encoder or the like in an apparatus for simulating the three-dimensional movement of the catheter.
  • the catheter may be twisted.
  • the rotation angle of the catheter at the operator's hand and the rotation angle of the distal end of the catheter do not match. Therefore, in the method of Patent Document 1, when the catheter is twisted, the direction of the distal end of the catheter cannot be accurately detected.
  • Patent Document 2 discloses a method for correcting rotation of the axis of the monitor display image due to twisting of the electronic endoscope by mounting a rotation detection sensor at the tip of the electronic endoscope in the electronic endoscope system. ing.
  • the rotation detection sensor is attached to the distal end of the catheter to be inserted into a blood vessel having a diameter of several millimeters, there is a risk of damaging the blood vessel when the rotation detection sensor contacts the blood vessel during the catheter operation.
  • the present invention solves the above-described conventional problems, and an object thereof is to provide an apparatus for measuring a rotation angle of a catheter tip in which twist is generated without attaching a rotation detection sensor to the tip of the catheter. .
  • a catheter tip rotation angle measuring device is a device for measuring the rotation angle of the tip of a catheter inserted into a living body tube, and has a cylindrical shape penetrating from a first hole to a second hole.
  • a movement restricting portion that restricts movement of the catheter in the insertion direction by passing the catheter from the first hole to the second hole, and a rotation angle of the catheter in the first hole.
  • a first rotation angle measurement unit that measures a first rotation angle
  • a second rotation angle measurement unit that measures a second rotation angle that is the rotation angle of the catheter in the second hole, and the first hole
  • a catheter insertion length measuring unit for measuring a catheter insertion length which is a length from the second hole to the distal end portion of the catheter; a first rotation angle measured by the first rotation angle measuring unit; and a second rotation. Measured from the angle measurement unit The rotation angle of the tip of the catheter is calculated using the angle difference from the two rotation angles, the catheter insertion length measured by the catheter insertion length measurement unit, and the distance between the first hole and the second hole.
  • a catheter tip rotation angle measuring unit is used to measure a first rotation angle
  • a second rotation angle measurement unit measures a second rotation angle that is the rotation angle of the catheter in the second hole
  • a catheter tip rotation angle measuring unit for measuring a catheter insertion length which is a length from the second hole to the distal end portion of the catheter.
  • the rotation angle of the catheter tip can be measured without mounting a rotation detection sensor or the like on the tip of the catheter. Therefore, since there is no need to attach a rotation detection sensor or the like to the distal end of the catheter, it is possible to reduce the possibility of damaging the blood vessel in the catheter operation that requires a delicate and delicate operation.
  • the catheter tip rotation angle measuring device may further include a presentation unit that presents the rotation angle of the catheter tip calculated by the catheter tip rotation angle calculation unit to the user.
  • the measured rotation angle of the distal end of the catheter can be presented to the user. Therefore, the surgeon can accurately grasp the direction of the distal end of the catheter from the measured rotation angle of the distal end of the catheter. Thereby, the operator can easily adjust the distal end of the catheter to the branching direction of the target blood vessel or the like at the branching portion of the blood vessel, and can perform a smooth catheter operation.
  • the operation time is shortened by smooth catheter operation, the irradiation time of X-rays irradiated to the patient is shortened. Therefore, the amount of exposure caused by X-ray irradiation to the patient can be reduced.
  • the catheter tip rotation angle calculation unit sets the first rotation angle measured by the first rotation angle measurement unit as ⁇ 1 and the second rotation angle measured by the second rotation angle measurement unit as ⁇ 2.
  • the catheter tip rotation angle calculation unit calculates the difference between the first rotation angle and the rotation angle of the catheter tip ( ⁇ 1 ) and the first rotation angle.
  • the relationship between the rotation angle and the angle difference ( ⁇ 1 ⁇ 2 ) between the second rotation angle is related to the catheter insertion length (L) and the distance (L 12 ) between the first hole and the second hole.
  • the proportionality can be used to measure the rotation angle of the catheter tip.
  • the catheter tip rotation angle measuring device further includes a first image obtaining unit that obtains an X-ray image captured by an X-ray imaging device that captures an area including the distal end portion of the catheter inserted into a living body tube. Detecting a distal end portion of the catheter from the X-ray image acquired from the first image acquiring portion, and an operation of the distal end portion of the catheter extracted from the distal end portion extracting portion of the catheter.
  • the catheter tip rotation angle calculator calculates the rotation angle of the catheter tip when the motion detector detects a motion, and the motion detector does not detect a motion. In this case, the rotation angle of the distal end of the catheter may not be calculated.
  • the catheter tip rotation angle measuring device is further provided with a plurality of imaging parameter information acquisition units that acquire imaging parameter information that is information related to imaging conditions and postures of the respective X-ray imaging devices at the time of imaging.
  • a second image acquisition unit that acquires a plurality of X-ray images including the distal end portion of the catheter imaged by the X-ray imaging apparatus, and a plurality of X-ray images acquired from the second image acquisition unit in each X-ray image Extracting the straight line portion of the catheter, using the extracted straight line portion of the catheter and the imaging parameter information when each X-ray image acquired from the imaging parameter information acquisition unit is captured, Catheter rotation which is a rotation axis with respect to the insertion direction of the catheter from the restored three-dimensional shape of the catheter.
  • a catheter rotation axis determination unit for determining a plurality of types of catheters, a catheter shape storage unit for storing the shape of the catheter tip for each rotation angle of the catheter tip, and a catheter tip rotation angle correction unit
  • the catheter tip rotation angle correction unit receives the first X-ray image from the plurality of X-ray images from the second image acquisition unit, and the catheter tip rotation angle calculation unit and the catheter tip rotation angle and the first X-ray image from the catheter tip rotation angle calculation unit.
  • the catheter A plurality of the shape of the distal end of the catheter, and a calibration image generating unit that generates a plurality of calibration images in which the acquired shape of the distal end of the catheter and the rotation angle are associated with each other, and the acquired from the acquiring unit
  • a comparison unit that compares the first X-ray image and a plurality of calibration images generated by the calibration image generation unit, and determines a calibration image having the highest degree of similarity with the first X-ray image;
  • a rotation angle determination unit that determines the rotation angle corresponding to the calibration image determined by the comparison unit as a rotation angle of the distal end of the catheter; and the presentation unit is determined by the rotation angle determination unit
  • FIG. 1 shows a functional configuration of a catheter tip rotation angle measuring apparatus 100 according to the first embodiment.
  • the catheter tip rotation angle measuring device 100 includes a first rotation angle measurement unit 102, a second rotation angle measurement unit 103, a catheter insertion length measurement unit 104, a catheter tip rotation angle calculation unit 105, Have Moreover, it has the movement control part 101 which is not shown in figure. Below, each component which the catheter tip rotation angle measuring apparatus 100 has is demonstrated.
  • the movement restricting portion 101 has a cylindrical shape penetrating from the first hole 1010 to the second hole 1011 and moves the catheter 10 by passing the catheter 10 from the first hole 1010 to the second hole 1011. Is regulated in the insertion direction.
  • FIG. 2 shows an example of an outline of the movement restriction unit 101.
  • the movement restricting portion 101 has a cylindrical shape penetrating from the first hole 1010 to the second hole 1011.
  • the first hole 1010 is the rear end of the movement restricting portion 101 when viewed from the insertion direction
  • the second hole 1011 is the front end of the movement restricting portion 101 when viewed from the insertion direction. That is, the first hole 1010 is a hole on the operator side, and the second hole 1011 is a hole on the patient side.
  • the movement of the catheter 10 is restricted in the insertion direction.
  • the insertion direction includes rotation around the axial direction in addition to the axial direction of the catheter 10. That is, that the movement of the catheter 10 is restricted in the insertion direction means that the catheter 10 cannot move in the radial direction orthogonal to the axial direction of the catheter 10.
  • the first hole 1010 and the second hole 1011 have a size that allows the catheter 10 to pass therethrough and are larger than the diameter of the catheter 10 by a predetermined size.
  • the first hole 1010 and the second hole 1011 are approximately the same size as the diameter of the catheter 10, and there is almost no gap in the radial direction perpendicular to the axial direction of the catheter 10 when the catheter 10 is inserted. .
  • the size of the first hole 1010 and the second hole 1011 is, for example, about 1 to 3 mm larger than the diameter of the catheter 10.
  • the sizes of the first hole 1010 and the second hole 1011 may differ depending on the diameter or type of the catheter 10.
  • the first hole 1010 and the second hole 1011 are about 3 mm larger than the diameter of the catheter 10 if the catheter 10 is thicker than a normal catheter, and if the catheter 10 is thinner than a normal catheter, It may be about 1 mm larger than the diameter of the catheter 10.
  • the cylindrical inner surface has little frictional resistance. Therefore, the catheter operation can be performed smoothly, and the twist of the catheter 10 can be prevented from occurring inside the movement restricting portion 101.
  • a distance L 12 from the first hole 1010 to the second hole 1011 is equal to the length of the movement restricting portion 101.
  • the length of the movement restricting portion 101 is about 10 cm.
  • the movement restricting unit 101 may be fixed to the base 1013 by the legs 1012. By fixing the movement restricting portion 101 to the base 1013, a first rotation angle and a second rotation angle described later can be measured with high accuracy.
  • the first rotation angle measurement unit 102 measures the first rotation angle ⁇ 1 that is the rotation angle of the catheter 10 in the first hole 1010.
  • the measurement location of the first rotation angle measurement unit 102 is within a predetermined range including the position of the first hole 1010.
  • the measurement location of the first rotation angle measurement unit 102 is preferably as close to the first hole 1010 as possible, for example, within a range of 3 mm from the first hole.
  • FIG. 3 shows an example of an outline of the first rotation angle measurement unit 102.
  • the first rotation angle measurement unit 102 includes an imaging unit 1020 and a calculation unit 1021.
  • FIG. 4 shows the movement restricting portion 101 as viewed from above.
  • a line 11 is drawn in advance on the surface of the catheter 10 in a direction parallel to the central axis of the catheter 10.
  • the imaging unit 1020 captures the range 1022.
  • the calculation unit 1021 extracts the position of the line 11 from the image captured by the imaging unit 1020, and acquires the rotation angle in the range 1022 from the position displacement of the line 11 with respect to the center of the catheter 10.
  • a specific method in which the first rotation angle measurement unit 102 measures the first rotation angle will be described.
  • the photographing unit 1020 is a camera, for example, and the calculation unit 1021 is a CPU of a PC, for example.
  • the imaging unit 1020 starts imaging upon receiving a detection start input from the operator.
  • the imaging unit 1020 continues to image the catheter 10 at regular intervals.
  • the calculation unit 1021 extracts the position of the line 11 in each image for a plurality of captured images. For example, when the catheter 10 is white and the line 11 is black, the black position is extracted.
  • FIG. 5 shows a photographed image when the photographing unit 1020 is a line sensor including 15 pixels.
  • the imaging unit 1020 is arranged so that the entire width of the catheter 10 (in the direction perpendicular to the axial direction) is imaged by these 15 pixels.
  • the calculation unit 1021 performs a process of acquiring “5” that is the number of the pixel 5.
  • the calculation unit 1021 determines whether the brightness of each pixel is equal to or less than a threshold value.
  • the brightness of the pixel takes a value from 0 to 255, where 0 is the darkest and 255 is the brightest, and the threshold is 50.
  • the number of the pixel whose brightness is equal to or less than the threshold value is set as the position of the line 11.
  • FIG. 6 shows the relationship between the position of the line 11 and the rotation angle of the catheter.
  • a circle on the right side of FIG. 6 is a cross-sectional view of the catheter 10.
  • the left side of FIG. 6 shows the photographing unit 1020.
  • the photographing unit 1020 has an element for photographing pixels 0 to 14, and an enlarged view of this is shown in the middle of FIG.
  • the boundary of the range captured by each pixel is indicated by a dotted line.
  • the pixel 7 is a pixel in which the center of the catheter range (here, pixel 0 to pixel 14) in the image photographed by the photographing unit 1020 is photographed.
  • the pixel 7 is referred to as a center pixel 1022.
  • a surface portion of the catheter 10 photographed by the center pixel 1022 is set as a center position 12.
  • Pixel 5 is the position of line 11 drawn on the surface of catheter 10.
  • the pixel 5 is referred to as a line position pixel 1023.
  • the angle ⁇ 1 is the rotation angle of the line 11 from the center position 12, that is, the first rotation angle that is the rotation angle of the catheter 10 in the first hole 1010.
  • the angle ⁇ 1 is expressed by the relational expression (1) from the definition of the trigonometric function.
  • ⁇ 1 arcsin (h / r) (2)
  • the first rotation angle measurement unit 102 does not need to be configured by the photographing unit 1020 and the calculation unit 1021, and may be configured by a rotary encoder or the like.
  • FIG. 7 the block diagram of the 1st rotation angle measurement part 102 at the time of measuring the rotation angle of the catheter 10 using a rotary encoder is shown.
  • FIG. 7 shows a state seen from the insertion direction of the catheter 10. Since the groove portion 13 provided in advance in the catheter 10 and the protrusion portion 1014 provided in advance in the movement restricting portion 101 are engaged with each other, the rotation of the catheter 10 and the rotation of the rotary encoder are synchronized. Ten rotation angles can be measured.
  • a rotary encoder an absolute value encoder or the like (Soichi Ogawa et al., Basic Robotics, 1998, p98, Tokyo Denki University Press) can be used.
  • the second rotation angle measurement unit 103 measures a second rotation angle ⁇ 2 that is the rotation angle of the catheter 10 in the second hole 1011.
  • the second rotation angle measurement unit 103 can be realized by the same configuration as the first rotation angle measurement unit 102.
  • FIG. 3 shows a schematic example of the second rotation angle measurement unit 103.
  • the second rotation angle measurement unit 103 includes an imaging unit 1030 and a calculation unit 1031.
  • a specific method in which the second rotation angle measurement unit 103 measures the second rotation angle ⁇ 2 using the imaging unit 1030 and the calculation unit 1031 is that the first rotation angle measurement unit 102 uses the first rotation angle ⁇ . This is the same as the method of measuring 1 .
  • the second rotation angle measurement unit 103 does not need to be configured by the imaging unit 1030 and the calculation unit 1031, and may be configured by a rotary encoder or the like.
  • a specific method in which the second rotation angle measurement unit 103 measures the second rotation angle ⁇ 2 using a rotary encoder is a method in which the first rotation angle measurement unit 102 measures the first rotation angle ⁇ 1. It is the same.
  • the catheter insertion length measurement unit 104 measures the catheter insertion length L, which is the length from the first hole 1010 or the second hole 1011 to the distal end portion of the catheter 10.
  • FIG. 8 shows an example of the outline of the catheter insertion length measurement unit 104.
  • FIG. 8 shows a case where the length from the first hole 1010 to the distal end portion 15 of the catheter 10 is the catheter insertion length L.
  • the catheter insertion length measurement unit 104 is configured to change, for example, a roller 1040 that interlocks with advancement and retraction of the catheter 10 by pressing the catheter 10, an encoder 1041 that measures the amount of rotation of the roller 1040, and a rotation angle value of the encoder 1041.
  • the calculation unit 1042 calculates the insertion length from the quantity.
  • the calculation unit 1042 sets the initial value “0” before starting the operation of the catheter 10.
  • the catheter 10 is passed through the roller 1040, and from that time, the calculation unit 1042 starts calculating the insertion length of the catheter 10 and calculates until the operation of the catheter 10 ends. to continue.
  • the catheter insertion length measurement unit 104 can continue to measure the insertion length L of the catheter 10 in real time during the operation of the catheter 10.
  • the distal end of the catheter 10 has a curved portion 15.
  • the bending portion 15 (hereinafter, also referred to as “catheter distal end portion 15”) is different from other catheters 10 in material and the like, and is easily bent.
  • the catheter insertion length measurement unit 104 measures the distance from the first hole 1010 to the bending portion 15 as the catheter insertion length L. That is, the length not including the distal end portion 15 of the catheter corresponds to the catheter insertion length L. By not including the distal end portion 15 of the catheter in the catheter insertion length L, the rotation angle ⁇ of the distal end of the catheter 10 can be calculated more accurately.
  • the catheter insertion length measuring unit 104 does not necessarily need to have a length that does not include the distal end portion 15 of the catheter as the catheter insertion length L, and measures the length that includes the distal end portion 15 of the catheter as the catheter insertion length L. Also good.
  • the catheter insertion length measuring unit 104 may be configured to take an image of a scale relating to the insertion position described on the surface of the catheter 10 by an imaging means such as a camera and calculate the insertion length based on the value of the scale. .
  • the catheter insertion length measurement unit 104 does not need to be provided independently of the movement restricting unit 101 and may be provided in the movement restricting unit 101.
  • the catheter insertion length measurement unit 104 may set the length from the second hole 1011 to the distal end portion of the catheter 10 as the catheter insertion length L.
  • the catheter tip rotation angle calculation unit 105 calculates an angle difference ( ⁇ between the first rotation angle ⁇ 1 measured by the first rotation angle measurement unit 102 and the second rotation angle ⁇ 2 measured by the second rotation angle measurement unit 103. 1 - [theta] 2), the catheter insertion length L measured from the catheter insertion length measurement unit 104, and the first hole 1010 with a distance L 12 between the second hole 1011, the rotation angle ⁇ of the tip of the catheter calculate.
  • the catheter tip rotation angle calculation unit 105 receives the first rotation angle ⁇ 1 from the first rotation angle measurement unit 102, the second rotation angle ⁇ 2 from the second rotation angle measurement unit 103, and the catheter insertion from the catheter insertion length measurement unit 104. Each of the lengths L is acquired. The catheter tip rotation angle calculation unit 105 calculates an angle difference ( ⁇ 1 ⁇ 2 ) between the acquired first rotation angle ⁇ 1 and second rotation angle ⁇ 2 . In addition, the operator sets in advance the distance L 12 between the first hole 1010 and the second hole 1011 in the catheter tip rotation angle calculation unit 105 before starting the operation of the catheter 10.
  • FIG. 10 shows the relationship between parameters when the first hole 1010 is on the operator side.
  • the first hole 1010 is on the operator side and the catheter insertion length L is the length from the first hole 1010 to the distal end portion 15 of the catheter 10 will be described.
  • the relationship of the expression (3) is established between the parameters.
  • Equation (3) represents the overall twist angle ( ⁇ 1 ⁇ ) of the catheter 10 and the first The relationship between the twist angle ( ⁇ 1 ⁇ 2 ) between the first hole 1010 and the second hole 1011 is proportional to the relationship between the catheter insertion length L and the distance L 12 between the first hole 1010 and the second hole 1011. Represents that. Therefore, the rotation angle ⁇ of the catheter tip can be calculated from the equation (4).
  • the catheter insertion length L is the second hole.
  • the rotation angle ⁇ of the catheter tip when the length from 1011 to the tip 15 of the catheter 10 can be calculated from the equation (7).
  • the presentation unit 106 presents the rotation angle ⁇ of the distal end of the catheter 10 calculated by the catheter tip rotation angle calculation unit 105 to the user.
  • the presentation unit 106 is not an essential component of the catheter tip rotation angle measuring device 100, but is a user by presenting the rotation angle ⁇ of the tip of the catheter 10 calculated by the catheter tip rotation angle calculation unit 105 to the user. Helps the operator operate the catheter.
  • the presentation unit 106 displays, for example, a numerical value on the display of the rotation angle ⁇ of the distal end of the catheter 10 calculated by the catheter tip rotation angle calculation unit 105. Moreover, you may display in three dimensions using CG etc. of the catheter 10 which reflected the rotation angle (theta) of the front-end
  • FIG. 11A shows an example in which the rotation angle ⁇ of the distal end of the catheter 10 is displayed at the lower right of the display using an arrow.
  • Step S100> The operator sets initial values of the first rotation angle measurement unit 102, the second rotation angle measurement unit 103, and the catheter insertion length measurement unit 104 to 0 as advance preparations for performing the treatment. In addition, the operator inputs the distance L 12 between the first hole 1010 and the second hole 1011 to the catheter tip rotation angle calculation unit 105. Then, the operator selects the catheter 10 suitable for the diagnosis site or the treatment site, passes the selected catheter 10 from the first hole 1010 to the second hole 1011 of the movement restricting portion 101, and further measures the catheter insertion length. Pass through part 104.
  • “0” in the first rotation angle measurement unit 102 and the second rotation angle measurement unit 103 may be, for example, 0 in the direction orthogonal to the insertion direction of the catheter 10 and from the floor to the ceiling.
  • the initial value of the 1st rotation angle measurement part 102, the 2nd rotation angle measurement part 103, and the catheter insertion length measurement part 104 was set to 0, it may not be 0.
  • Each value may be measured from the numerical value at the start of measurement and the numerical value during measurement.
  • ⁇ Step S200> The surgeon inserts the catheter 10 into the patient.
  • Step S300> The surgeon operates the catheter 10 while monitoring the X-ray image of the catheter 10 projected on the display, and inserts the catheter 10 to the diagnosis site or treatment site.
  • Step S400> When the operator wants to recognize the direction of the distal end of the catheter 10, the operator performs an input for starting measurement of the rotation angle ⁇ of the distal end of the catheter 10.
  • Examples of the device that accepts an input include a keyboard and a push button.
  • the catheter tip rotation angle measuring apparatus 100 executes steps S410 to S450 described later, and presents the rotation angle of the tip of the catheter 10 to the operator.
  • the measurement of the rotation angle of the distal end of the catheter 10 is not limited to being started after waiting for input, and may be measured at a predetermined time interval.
  • ⁇ Step S500> The surgeon continues to operate the catheter 10 with reference to the rotation angle ⁇ of the tip of the catheter presented by the presentation unit 106.
  • Step S600> The surgeon executes step S300 when the diagnosis or treatment is not completed, and executes step S700 when the diagnosis or treatment is completed.
  • Step S700> When completing the diagnosis or treatment, the operator removes the catheter 10 from the patient and finishes the treatment.
  • the method of measuring the rotation angle ⁇ of the distal end of the catheter 10 is simply to measure the rotation angle ⁇ of the distal end of the catheter 10. That is, the method does not include the operator performing the catheter operation in step S500 based on the rotation angle ⁇ of the distal end of the catheter measured in step S400. Therefore, the method for measuring the rotation angle ⁇ of the distal end of the catheter 10 is not a method for performing a surgical process on the human body and does not correspond to a method for operating a human.
  • FIG. 13 shows a flow of processing operations of the catheter tip rotation angle measuring apparatus 100 in the first embodiment.
  • the first rotation angle measurement unit 102 measures a first rotation angle that is the rotation angle of the catheter 10 in the first hole 1010. Specific methods of the first rotation angle theta 1 of the measurement is as described above. Further, the first rotation angle measurement unit 102 outputs the measured first rotation angle ⁇ 1 to the catheter tip rotation angle calculation unit 105.
  • the second rotation angle measurement unit 103 measures a second rotation angle that is the rotation angle of the catheter 10 in the second hole 1011. Specific method of the second rotation angle theta 2 of the measurement are as described above. Further, the second rotation angle measurement unit 103 outputs the measured second rotation angle ⁇ 2 to the catheter tip rotation angle calculation unit 105.
  • the catheter insertion length measurement unit 104 measures the catheter insertion length, which is the length from the first hole 1010 to the distal end portion 15 of the catheter 10.
  • the specific method for measuring the catheter insertion length L is as described above. Further, the catheter insertion length measurement unit 104 outputs the measured catheter insertion length L to the catheter tip rotation angle calculation unit 105.
  • the catheter tip rotation angle calculation unit 105 calculates an angle difference ( ⁇ between the first rotation angle ⁇ 1 measured by the first rotation angle measurement unit 102 and the second rotation angle ⁇ 2 measured by the second rotation angle measurement unit 103. 1 - [theta] 2), the catheter insertion length L measured from the catheter insertion length measurement unit 104, and the first hole 1010 with a distance L 12 between the second hole 1011, the rotation angle ⁇ of the tip of the catheter calculate.
  • the catheter tip rotation angle calculation unit 105 calculates the first rotation angle ⁇ 1 measured in step S410, the second rotation angle ⁇ 2 measured in step S420, and the catheter insertion length L measured in step S430. Get each.
  • the catheter tip rotation angle calculation unit 105 an angular difference from the first rotation angle theta 1 and the second rotation angle theta 2 which the acquired first rotation angle theta 1 and the second rotation angle theta 2 (theta 1 - ⁇ 2 ) is calculated.
  • the catheter tip rotation angle calculator 105 calculates the angle difference ( ⁇ 1 ⁇ 2 ) between the calculated first rotation angle ⁇ 1 and the second rotation angle ⁇ 2 , the acquired catheter insertion length L, and step S100. in using the first hole 1010, which is set a distance L 12 between the second hole 1011, (4) from equation to calculate the rotation angle ⁇ of the tip of the catheter 10.
  • the catheter tip rotation angle calculation unit 105 outputs the calculated rotation angle ⁇ of the tip of the catheter 10 to the presentation unit 106.
  • Step S450> The presentation unit 106 acquires the rotation angle ⁇ of the distal end of the catheter 10 calculated in step S440, and presents the acquired rotation angle ⁇ of the distal end of the catheter 10 to the user.
  • the specific method of presentation by the presentation unit 106 is as described above.
  • ⁇ Effect> As described above, according to the catheter tip rotation angle measuring device 100 in the first embodiment, when the catheter 10 is twisted uniformly, the catheter 10 is not attached to the tip of the catheter 10 without attaching a rotation detection sensor or the like. The rotation angle of the tip of 10 can be measured. Thereby, since there is no need to attach a rotation detection sensor or the like to the distal end of the catheter 10, the possibility of damaging a blood vessel in the operation of the catheter 10 requiring a delicate and delicate operation can be reduced.
  • the operator can accurately grasp the orientation of the distal end of the catheter 10 from the measured rotation angle ⁇ of the distal end of the catheter 10. I can do it.
  • the operator can easily adjust the distal end of the catheter 10 in the branching direction of the target blood vessel or the like at the branching portion of the blood vessel, and can perform a smooth catheter operation.
  • the operation time is shortened by smooth catheter operation, the irradiation time of X-rays irradiated to the patient is shortened. Therefore, the amount of exposure caused by X-ray irradiation to the patient can be reduced.
  • FIG. 29 shows an example of an image when the guide wire 16 is advanced into the catheter 10 through the guide wire 16 through the lumen of the catheter 10 inserted into the blood vessel.
  • the guide wire 16 is preceded before insertion of the catheter 10 to make it difficult to damage the blood vessel wall and to smoothly carry the catheter 10 to the lesioned part.
  • the surface of the lumen of the catheter 10 has a small frictional resistance, or the surgeon injects a small amount of water into the lumen of the catheter 10 during treatment to facilitate passage through the guide wire. Therefore, when the guide wire 16 is twisted, a load is applied to the rotation of the distal end of the guide wire 16, and the twist of the guide wire 16 caused by the twist is uniformly generated in the lumen of the catheter 10. It is done.
  • the parent catheter makes it easy to guide the latter child catheter, the thickness is about 2 mm, and the child catheter is for actually treating the lesioned part, and the thickness is 1 mm.
  • the parent catheter makes it easy to guide the latter child catheter, the thickness is about 2 mm, and the child catheter is for actually treating the lesioned part, and the thickness is 1 mm.
  • the child catheter By passing the child catheter through the lumen of the parent catheter inserted into the blood vessel, the child catheter can be smoothly carried to the lesion, and the child catheter is used to treat the lesion, such as expanding the narrowed blood vessel with a balloon. be able to. Even in this case, when twisting occurs in the child catheter, a load is applied to the rotation of the distal end of the child catheter, and it is considered that the twisting of the child catheter is caused uniformly in the lumen of the parent catheter. It is done.
  • the catheter tip rotation angle measuring apparatus 100 can be used particularly effectively in a guide wire inserted prior to the catheter or a child catheter inserted beyond the parent catheter. .
  • the catheter tip rotation angle measuring device 200 detects the movement of the tip 15 of the catheter 10 and, when motion is detected, calculates the rotation angle ⁇ of the tip of the catheter 10 and no motion is detected. In this case, the point of rotation angle ⁇ of the catheter 10 is not calculated, which is different from the catheter tip rotation angle measuring apparatus 100 of the first embodiment. Below, it demonstrates focusing on a different point from 1st Embodiment.
  • FIG. 14 shows a functional configuration of the catheter tip rotation angle measuring apparatus 200 in the second embodiment.
  • the same components as those in FIG. 1 are denoted by the same reference numerals, and description thereof will be omitted as appropriate.
  • the catheter tip rotation angle measuring device 200 includes a first rotation angle measurement unit 102, a second rotation angle measurement unit 103, a catheter insertion length measurement unit 104, a catheter tip rotation angle calculation unit 205, , A presentation unit 106, a first image acquisition unit 207, a catheter tip portion extraction unit 208, and a motion detection unit 209.
  • a first rotation angle measurement unit 102 a second rotation angle measurement unit 103
  • a catheter insertion length measurement unit 104 the catheter tip rotation angle calculation unit 205
  • a presentation unit 106 a first image acquisition unit 207
  • motion detection unit 209 a motion detection unit
  • the first image acquisition unit 207 acquires the X-ray image 30 captured by the X-ray imaging apparatus 20 that captures an area including the distal end portion of the catheter 10 inserted into the living body tube.
  • FIG. 15 shows the configuration of the first image acquisition unit 207 and the X-ray imaging apparatus 20.
  • the X-ray imaging apparatus 20 is an apparatus that acquires an X-ray image captured by irradiating X-rays to an imaging target region of a patient or a blood vessel image captured when a contrast agent is injected.
  • the X-ray imaging apparatus 20 is an apparatus that acquires an X-ray image captured by irradiating X-rays to an imaging target region of a patient or a blood vessel image captured when a contrast agent is injected.
  • it is called an X-ray angiography apparatus or angiography.
  • the X-ray imaging apparatus 20 includes an X-ray generation unit 21, an X-ray detection unit 22, a mechanism unit 23, a mechanism control unit 24, an arm 26, and a database 25.
  • the X-ray generation unit 21 includes an X-ray tube that generates X-rays using a high voltage and an X-ray diaphragm that controls an irradiation field by shielding a part of the X-rays.
  • the patient 28 is irradiated with X-rays.
  • the X-ray detection unit 22 receives X-rays transmitted through the patient 28, records image information, and outputs the recorded image information. For example, an X-ray sensing layer is disposed, and the X-rays are converted into digital data.
  • the X-ray detection unit 22 outputs image information indicating the irradiated X-ray image 30 to the database 25.
  • the mechanism unit 23 moves the arm 26 and the bed 27 based on the instruction of the mechanism control unit 24 that has received the operator's operation instruction.
  • the first image acquisition unit 207 reads an X-ray image or an X-ray image stored in the database 25 and outputs it to the catheter tip portion extraction unit 208.
  • FIG. 16 shows an example of the X-ray image 30 of the catheter 10 acquired by the first image acquisition unit 207.
  • the distal end portion 15 of the catheter is a portion including the bending portion 15 shown in FIG. 9, for example, a portion from the distal end of the catheter 10 to 1 to 2 cm in addition to the length La of the bending portion 15.
  • the distal end portion 15 of the catheter may be the same portion as the bending portion 15.
  • the catheter tip extraction unit 208 extracts the catheter tip 15 from the X-ray image 30 acquired from the first image acquisition unit 207.
  • the catheter tip portion extraction unit 208 compares the acquired X-ray image 30 with the X-ray image of the previous frame, and extracts only the portion where the lightness is changed for each pixel. Thereby, only when the catheter 10 moves, the distal end portion 15 of the catheter can be extracted. For example, when the brightness of a pixel takes a value of “0 to 255”, “0” is the darkest and “255” is the brightest, the pixel of the current frame is subtracted from the pixel of the previous frame of the X-ray image 30.
  • FIG. 17 shows an image of the distal end 15 of the catheter 10 extracted by the catheter distal end extracting unit 208.
  • the catheter tip portion extraction unit 208 may compare the frame with a certain time before (for example, one second before).
  • the catheter distal end extraction unit 208 is not limited to the above method, and determines the linearity of the catheter 10, recognizes a portion having a low linearity value as the bending portion 15, and includes a predetermined bending portion 15 including the recognized bending portion 15. The region may be extracted as the distal end portion 15 of the catheter 10.
  • the motion detection unit 209 detects the motion of the distal end portion 15 of the catheter 10 extracted by the catheter distal end portion extraction unit 208.
  • the motion detection unit 209 can indirectly detect the insertion or rotation of the proximal portion by detecting the motion of the distal end portion 15 of the catheter 10.
  • the hand portion is a portion including the movement restricting portion 101 and the catheter insertion length measuring portion 104.
  • the motion detection unit 209 detects the motion of the distal end portion 15 when the number of pixels whose brightness has changed in the extraction of the distal end portion 15 of the catheter 10 is greater than or equal to a predetermined number of pixels. For example, the motion detection unit 209 detects the motion of the tip 15 when the brightness is “0” and the number of pixels is “10” or more. In addition, when extracting the distal end portion 15 of the catheter 10, when the distal end portion 15 of the catheter 10 does not move, the motion detection unit 209 extracts the brightness of all the pixels as “255” and detects the motion of the distal end portion 15. Not.
  • FIG. 18 shows whether the distal end portion 15 of the catheter 10 is detected by the motion detection unit 209 by comparing the previous frame and the current frame.
  • the motion detection unit 209 is not limited to the above method, and compares the tip portion 15 drawn from the X-ray image 30 with a marginal edge with a predetermined time interval (for example, an interval of 1 ms) to determine a predetermined value (for example, , 0.03 mm) or more may be detected, and the operation may not be detected if there is no change greater than a predetermined value.
  • a predetermined time interval for example, an interval of 1 ms
  • the catheter tip rotation angle calculation unit 205 calculates an angular difference ( ⁇ between the first rotation angle ⁇ 1 measured by the first rotation angle measurement unit 102 and the second rotation angle ⁇ 2 measured by the second rotation angle measurement unit 103. 1 - [theta] 2), the catheter insertion length L measured from the catheter insertion length measurement unit 104, and the first hole 1010 with a distance L 12 between the second hole 1011, the rotation angle ⁇ of the tip of the catheter calculate.
  • the catheter tip rotation angle calculation unit 205 calculates the rotation angle ⁇ of the tip of the catheter 10 when the motion detection unit 209 detects motion, and when the motion detection unit 209 does not detect motion, The rotation angle ⁇ of the tip is not calculated.
  • the specific calculation method of the distal end rotation angle ⁇ of the catheter 10 calculated by the catheter distal end rotation angle calculation unit 205 is the same as the calculation method of the catheter distal end rotation angle calculation unit 105 in the first embodiment.
  • FIG. 19 shows a flow of processing operations of the catheter tip rotation angle measuring apparatus 200 in the second embodiment.
  • the same steps as those in FIG. 13 are denoted by the same reference numerals, and description thereof will be omitted as appropriate.
  • the first image acquisition unit 207 acquires the X-ray image 30 captured by the X-ray imaging apparatus 20 from the database 25.
  • the specific method for acquiring the image is as described above.
  • the first image acquisition unit 207 outputs the acquired X-ray image 30 to the catheter tip portion extraction unit 208.
  • the catheter tip portion extraction unit 208 extracts the tip portion 15 of the catheter 10 from the X-ray image 30 acquired from the first image acquisition unit 207. First, the catheter tip extraction unit 208 acquires the X-ray image 30 acquired from step S460. Next, the catheter tip portion extraction unit 208 extracts the tip portion 15 of the catheter 10 from the acquired X-ray image 30. The specific method for extracting the tip 15 is as described above. Next, the catheter distal end extraction unit 208 outputs the extracted distal end 15 to the motion detection unit 209.
  • Step S480> The motion detector 209 detects the motion of the catheter tip 15 extracted by the catheter tip extractor 208. First, the motion detection unit 209 acquires the tip 15 extracted from step S470. Next, the motion detection unit 209 detects the motion of the acquired tip portion 15. A specific method for detecting the operation of the distal end portion 15 is as described above. Next, the motion detection unit 209 executes step S440 when the motion of the tip portion 15 is detected, and executes step S450 when the motion of the tip portion 15 is not detected. That is, when the motion detection unit 209 does not detect the motion of the distal end portion 15, the presentation unit 105 continues to present the rotation angle ⁇ of the distal end of the catheter 10 currently being presented to the user.
  • the catheter tip rotation angle measuring apparatus 200 in the second embodiment when the operation of the tip portion 15 of the catheter 10 is detected, the rotation angle ⁇ of the tip of the catheter 10 is measured. . Therefore, the surgeon can grasp the rotation angle ⁇ of the distal end of the catheter 10 in real time in accordance with his / her own catheter operation. Thus, the operator can easily adjust the distal end of the catheter 10 in the branching direction of the target blood vessel or the like at the branching portion of the blood vessel, and can perform a smooth catheter operation. In addition, since the operation time is shortened by smooth catheter operation, the irradiation time of X-rays irradiated to the patient is shortened. Therefore, the amount of exposure caused by X-ray irradiation to the patient can be reduced.
  • the catheter tip rotation angle measuring apparatus 300 is capable of accurately measuring the rotation angle of the tip of the catheter 10 even when a part of the catheter 10 is twisted. Different from the catheter tip rotation angle measuring apparatus 100 of FIG. below, it demonstrates focusing on a different point from 1st Embodiment.
  • FIG. 20 shows a functional configuration of the catheter tip rotation angle measuring apparatus 300 in the third embodiment.
  • the same components as those in FIG. 1 are denoted by the same reference numerals, and description thereof will be omitted as appropriate.
  • the catheter tip rotation angle measuring device 300 includes a first rotation angle measurement unit 102, a second rotation angle measurement unit 103, a catheter insertion length measurement unit 104, a catheter tip rotation angle calculation unit 305, A presentation unit 306, an imaging parameter information acquisition unit 310, a second image acquisition unit 311, a catheter rotation axis determination unit 312, a catheter shape storage unit 313, and a catheter tip rotation angle correction unit 314. Moreover, it has the movement control part 101 which is not shown in figure. Below, each component which the catheter front-end
  • the imaging parameter information acquisition unit 310 acquires imaging parameter information that is information relating to imaging conditions and postures of the respective X-ray imaging devices 20 at the time of imaging.
  • the imaging parameter information is an apparatus parameter (internal parameter) of the X-ray imaging apparatus 20 and a position / orientation parameter (external parameter) at the time of imaging of the X-ray imaging apparatus 20.
  • the internal parameter is a parameter related to imaging of the X-ray imaging apparatus 20, and for example, shutter speed, aperture value, exposure correction value, ISO sensitivity value, white balance setting value, recording image size, and recording It is a parameter that affects image quality and other image creation.
  • the external parameter is a parameter related to the position or orientation of the X-ray imaging apparatus 20 at the time of imaging. For example, when the origin and the directions of the X, Y, and Z axes are determined in space, the X-ray imaging apparatus 20 The position is represented by (x, y, z).
  • the imaging parameter information acquisition unit 310 outputs the acquired imaging parameter information to the outside as a perspective projection row example.
  • the perspective projection matrix will be described. Assuming that the point X in the three-dimensional space is projected onto the point m in the two-dimensional perspective image, this projection can be expressed as in equation (8) using homogeneous coordinates.
  • the image point m in the equation (8) is a 3 ⁇ 1 vector
  • the three-dimensional point X is a 4 ⁇ 1 vector
  • the matrix P of the 3 ⁇ 4 matrix is called a perspective projection matrix or a projection matrix.
  • is a real number indicating indefiniteness of a constant multiple.
  • the projection matrix P indicates a geometric relationship such as a focal length at the time of fluoroscopic image shooting, a three-dimensional rotation and a parallel movement from the three-dimensional space coordinate system to the fluoroscopic image coordinate system.
  • the 2nd image acquisition part 311 images the area
  • the second image acquisition unit 311 has substantially the same configuration as the first image acquisition unit 207.
  • the second image acquisition unit 311 acquires a plurality of X-ray images 30 captured by a plurality of X-ray imaging devices 20 arranged at different positions.
  • the plurality of X-ray imaging devices 20 arranged at different positions include, for example, two sets of an X-ray generation unit 21 and an X-ray detection unit 22, and the imaging angle difference between them is 90 degrees.
  • This is a biplane type X-ray imaging apparatus capable of simultaneously capturing X-ray images.
  • the different positions mean that the X-ray detection units 22 are installed at the same location and the cases where the X-ray detection units 22 are installed at opposite locations.
  • the angle difference should be 5 degrees or more.
  • the catheter rotation axis determination unit 312 extracts the straight line portion of the catheter 10 in each X-ray image from the plurality of X-ray images acquired from the second image acquisition unit 311, and extracts the straight line portion of the catheter 10 and the parameters Using the parameter information obtained when each X-ray image acquired from the information acquisition unit 310 is captured, the three-dimensional shape of the catheter 10 is restored. From the restored three-dimensional shape of the catheter 10, the insertion direction of the catheter 10 is restored. Determine the axis of rotation. The catheter rotation axis determination unit 312 outputs the determined rotation axis and the X-ray image acquired from the second image acquisition unit 311 to the outside.
  • the catheter rotation axis determination unit 312 extracts a straight portion of the catheter 10 in each X-ray image 30 from the plurality of X-ray images acquired from the second image acquisition unit 311.
  • a specific extraction method will be described using an example in which two X-ray images 30 are acquired from the second image acquisition unit 311.
  • the catheter rotation axis determination unit 312 acquires two X-ray images 30 having an angle difference from the second image acquisition unit 311.
  • the catheter rotation axis determination unit 312 extracts a straight portion of the catheter 10 for each of the two acquired X-ray images 30.
  • the straight portion of the catheter 10 is a portion not including the bending portion 15 of the catheter 10 as shown in FIG.
  • the catheter rotation axis determination unit 312 determines the linearity of the catheter 10, recognizes a portion having a predetermined linearity or more as the straight portion 14, and sets a predetermined region including the recognized straight portion 14 as a straight portion of the catheter 10. Extract.
  • FIG. 21 shows a schematic diagram for determining the rotation axis of the catheter 10.
  • the catheter rotation axis determination unit 312 acquires a straight line image 140 and a straight line image 141 that are images obtained by extracting straight portions of the catheter 10 from the two X-ray images 30. Further, the catheter rotation axis determination unit 312 includes the X-ray generation unit 210 and the X-ray generation unit 211, the X-ray detection unit 220, and the X-ray detection unit 221 when the linear part image 140 and the linear part image 141 are acquired. Is acquired from the imaging parameter information acquisition unit 310.
  • the catheter rotation axis determination unit 312 performs the three-dimensional shape restoration of the straight portion of the catheter 10 using the two straight portion images 140 and the straight portion image 141 and two sets of imaging parameter information.
  • a specific method for restoring the three-dimensional shape is as follows. First, in order to extract corresponding points between two images, a Fundamental matrix indicating a geometrical relationship between the two images is calculated from the imaging parameter information using equation (9).
  • FIG. 26 shows an outline of the epipolar line. As shown in FIG. 26, (a) in the figure is a point on a straight line in the straight line portion image 140, and a dotted line in (b) in the figure is an epipolar line with respect to a straight point.
  • the epipolar line indicates the existence range in the straight line image 140. Therefore, the corresponding point between the straight line portion image 140 and the straight line portion image 141 can be calculated by calculating the intersection of the epipolar line and the extracted straight line. Since the projection formula like Formula (8) is established at the corresponding points between the imaging parameter information (projection matrix of each X-ray tube) and the straight line portion image 140 and the straight line portion image 141, each projection formula is expressed as a three-dimensional point. Summarizing X, it can be transformed as shown in equation (11).
  • a three-dimensional point can be restored from corresponding points between the imaging parameter information (each projection matrix) and the straight line portion image 140 and the straight line portion image 141 by using the expression (12) obtained by modifying the expression (11). it can.
  • M is an abbreviation of the 4 ⁇ 3 matrix on the left side of Equation (11)
  • b is an abbreviation of the 4 ⁇ 1 matrix on the right side of Equation (11).
  • FIG. 22 shows the restored three-dimensional straight line of the catheter 10.
  • the catheter rotation axis determination unit 312 determines the rotation axis with respect to the insertion direction of the catheter 10 from the three-dimensional shape of the catheter 10 restored by the above method.
  • the catheter shape storage unit 313 stores the shape of the distal end of the catheter 10 for each of the rotation angles of the distal end of the catheter 10 for a plurality of types of catheters 10.
  • catheters 10 or guidewires There are many types of catheters 10 or guidewires inserted prior to them, and there are various shapes depending on the application. For example, there is a catheter called Judkins for the right coronary artery or the left coronary artery, and the tip is bent according to each anatomical position, and the shape is easy to insert.
  • the catheter shape storage unit 313 stores various three-dimensional shape models corresponding to them in advance. Then, a surgeon or a radiologist who handles radiation equipment sets the type in the catheter shape storage unit 313 in advance before inserting the catheter 10 or the guide wire.
  • the catheter tip rotation angle calculation unit 305 includes an angle difference ( ⁇ between the first rotation angle ⁇ 1 measured by the first rotation angle measurement unit 102 and the second rotation angle ⁇ 2 measured by the second rotation angle measurement unit 103. 1 - [theta] 2), the catheter insertion length L measured from the catheter insertion length measurement unit 104, and the first hole 1010 with a distance L 12 between the second hole 1011, the rotation angle ⁇ of the tip of the catheter calculate.
  • the catheter tip rotation angle calculation unit 305 outputs the second rotation angle ⁇ 2 measured by the second rotation angle measurement unit 103 to the outside in addition to the calculated catheter tip rotation angle ⁇ .
  • FIG. 23 shows a functional configuration of the catheter tip rotation angle correction unit 314.
  • the catheter tip rotation angle correction unit 314 includes an acquisition unit 3141, a calibration image generation unit 3142, a comparison unit 3143, and a rotation angle determination unit 3144. Below, each component which the catheter front-end
  • the acquisition unit 3141 obtains the first X-ray image from the plurality of X-ray images from the second image acquisition unit 311, and the rotation angle ⁇ and the second rotation angle ⁇ 2 of the distal end of the catheter 10 from the catheter tip rotation angle calculation unit 305.
  • the catheter rotation axis determination unit 312 acquires the catheter rotation axis, and the imaging parameter information acquisition unit 310 acquires the imaging parameter information when the first X-ray image is captured.
  • the calibration image generation unit 3142 receives the shape of the tip of the catheter 10 for each rotation angle from the rotation angle ⁇ of the tip of the catheter 10 acquired from the acquisition unit 3141 from the catheter shape storage unit 313 to the second rotation angle ⁇ 2. And a plurality of calibration images in which the acquired tip shape of the catheter 10 is associated with the rotation angle ⁇ a are generated.
  • FIG. 24 shows an example of a calibration image. A specific method by which the calibration image generation unit 3142 generates a calibration image will be described.
  • the calibration image generation unit 3142 matches the rotation axis of the three-dimensional shape model acquired from the catheter shape storage unit 313 with the rotation axis acquired from the catheter rotation axis determination unit 312.
  • the calibration image generation unit 3142 uses either one of the two-direction imaging parameter information acquired from the imaging parameter information acquisition unit 310 to generate a three-dimensional shape model in which the rotation axes are aligned using the imaging parameter information.
  • an image in which the rotation angle ⁇ of the tip of the catheter 10 corresponds to the rotation angle ⁇ calculated by the catheter tip rotation angle calculation unit 305 to the second rotation angle ⁇ 2 is constant.
  • the reason why the generation range is changed from ⁇ to ⁇ 2 is that the catheter is calculated on the assumption that there is a uniform twist when some or all of the twist is eliminated and a partial twist exists. This is because an actual rotation angle exists between the rotation angle ⁇ of the tip of 10 and the second rotation angle ⁇ 2 .
  • the comparison unit 3143 compares the first X-ray image acquired from the acquisition unit 3141 with the plurality of calibration images generated by the calibration image generation unit 3142, and the calibration having the highest degree of similarity with the first X-ray image. Determine the image for the action.
  • FIG. 25 shows a schematic diagram for comparing the X-ray image and the calibration image.
  • the comparison unit 3143 extracts feature points such as edges or corners from the first X-ray image, and calculates a feature amount from the surrounding area. Similarly, the comparison unit 3143 calculates feature amounts for a plurality of calibration images. The comparison unit 3143 compares the feature amount of the first X-ray image with each feature amount of the calibration image, and selects the calibration image having the highest feature amount coincidence with the first X-ray image. It is determined that the image is a large calibration image.
  • the method for calculating the similarity between images performed by the comparison unit 3143 may be any method for calculating the similarity between monochrome images.
  • the rotation angle determination unit 3144 determines the rotation angle ⁇ a corresponding to the calibration image determined by the comparison unit 3143 as the rotation angle ⁇ of the distal end of the catheter 10.
  • ⁇ Presentation unit 306> The presentation unit 306 presents the rotation angle ⁇ of the distal end of the catheter 10 determined by the rotation angle determination unit 3144 to the user. At this time, as shown in FIG. 11B, in order to indicate that calibration (correction) is being performed to a user such as a doctor or an engineer, “being corrected” may be displayed.
  • FIG. 27 shows a flow of processing operations of the catheter tip rotation angle measuring apparatus 300 in the third embodiment.
  • steps similar to those in FIG. 13 are denoted by the same reference numerals, and description thereof will be omitted as appropriate.
  • the imaging parameter information acquisition unit 310 acquires imaging parameter information that is information relating to imaging conditions and postures of the respective X-ray imaging devices 20 at the time of imaging. Further, the imaging parameter information acquisition unit 310 outputs the acquired imaging parameter information to the catheter tip rotation angle correction unit 314.
  • the 2nd image acquisition part 311 images the area
  • the second image acquisition unit 311 outputs the acquired plurality of X-ray images 30 to the catheter rotation axis determination unit 312.
  • the catheter rotation axis determination unit 312 extracts the straight line portion of the catheter 10 in each X-ray image from the plurality of X-ray images acquired from the second image acquisition unit 311, and extracts the straight line portion of the catheter 10 and the parameters Using the parameter information obtained when each X-ray image acquired from the information acquisition unit 310 is captured, the three-dimensional shape of the catheter 10 is restored. From the restored three-dimensional shape of the catheter 10, the insertion direction of the catheter 10 is restored. Determine the axis of rotation. The catheter rotation axis determination unit 312 outputs the determined rotation axis and the X-ray image acquired from the second image acquisition unit 311 to the outside. The catheter rotation axis determination unit 312 outputs the determined rotation axis and X-ray image to the catheter tip rotation angle correction unit 314.
  • the catheter tip rotation angle correction unit 314 corrects the rotation angle of the tip of the catheter 10. Further, the catheter tip rotation angle correction unit 314 outputs the corrected result to the presentation unit 306. A specific description of the operation of the catheter tip rotation angle correction unit 314 will be described later.
  • the presentation unit 306 presents the rotation angle ⁇ of the distal end of the catheter 10 determined by the rotation angle determination unit 3144 to the user.
  • the above operation is the processing operation of the catheter tip rotation angle measuring apparatus 300 in the third embodiment.
  • FIG. 28 shows a flow of processing operations of the catheter tip rotation angle correction unit 314.
  • the acquisition unit 3141 obtains the first X-ray image from the plurality of X-ray images from the second image acquisition unit 311, and the rotation angle ⁇ and the second rotation angle ⁇ 2 of the distal end of the catheter 10 from the catheter tip rotation angle calculation unit 305.
  • the catheter rotation axis determination unit 312 acquires the catheter rotation axis
  • the imaging parameter information acquisition unit 310 acquires the imaging parameter information when the first X-ray image is captured.
  • the acquisition unit 3141 outputs the acquired information and the like to the calibration image generation unit 3142 and the comparison unit 3143.
  • the timing at which the acquisition unit 3141 acquires each piece of information is every predetermined time interval.
  • the acquisition unit 3141 may acquire the catheter insertion length L from the catheter insertion length measurement unit 104 and may acquire the catheter insertion length L when there is a change.
  • the calibration image generation unit 3142 receives the shape of the tip of the catheter 10 for each rotation angle from the rotation angle ⁇ of the tip of the catheter 10 acquired from the acquisition unit 3141 from the catheter shape storage unit 313 to the second rotation angle ⁇ 2. And a plurality of calibration images in which the acquired tip shape of the catheter 10 is associated with the rotation angle ⁇ a are generated. Further, the calibration image generation unit 3142 outputs the generated plurality of calibration images to the comparison unit 3143.
  • the comparison unit 3143 compares the first X-ray image acquired from the acquisition unit 3141 with the plurality of calibration images generated by the calibration image generation unit 3142, and the calibration having the highest degree of similarity with the first X-ray image. Determine the image for the action.
  • the comparison unit 3143 also outputs the determined calibration image to the rotation angle determination unit 3144.
  • the rotation angle determination unit 3144 determines the rotation angle ⁇ a corresponding to the calibration image determined by the comparison unit 3143 as the rotation angle ⁇ of the distal end of the catheter 10. In addition, the rotation angle determination unit 3144 outputs the determined rotation angle ⁇ of the distal end of the catheter 10 to the presentation unit 106.
  • the tip rotation angle measuring apparatus 300 in the third embodiment even when a part of the catheter 10 is twisted, the tip rotation angle of the catheter 10 can be accurately measured. it can. Therefore, the surgeon can grasp the rotation angle ⁇ of the distal end of the catheter 10 even when a part of the catheter 10 is twisted. Thereby, the operator can easily adjust the distal end of the catheter to the branching direction of the target blood vessel or the like at the branching portion of the blood vessel, and can perform a smooth catheter operation. In addition, since the operation time is shortened by smooth catheter operation, the irradiation time of X-rays irradiated to the patient is shortened. Therefore, the amount of exposure caused by X-ray irradiation to the patient can be reduced.
  • a catheter tip rotation angle measuring device 400 includes a first rotation angle measurement unit 102, a second rotation angle measurement unit 103, a catheter insertion length measurement unit 104, a catheter tip rotation angle calculation unit 205, The presentation unit 106, the first image acquisition unit 207, the catheter tip extraction unit 208, and the determination unit 230.
  • the second embodiment is the same as the first embodiment. Therefore, detailed description is omitted, and only different configurations will be described below.
  • the first image acquisition unit 207 acquires an X-ray image captured by the X-ray imaging apparatus 20 that captures an area including the distal end portion of the catheter 10 inserted into the living body tube.
  • the catheter tip extraction unit 208 extracts the tip of the catheter 10 from the X-ray image acquired from the first image acquisition unit 207.
  • the determination unit 230 performs linear image determination as to whether or not the distal end portion of the catheter 10 extracted by the catheter distal end portion extraction unit 208 is a linear image on the X-ray image. Further, when the determination unit 230 determines that the distal end portion of the catheter 10 is a linear image, the concentration change determination is further performed on whether or not there is a concentration change at a predetermined threshold value or more at the distal end portion. As a result of the determination by the determination unit 230, when it is determined that the image is a linear image and there is no change in the concentration at the distal end, the determination unit 230 determines whether to stop calculating the catheter tip rotation angle.
  • FIG. 31A to 31C show examples of X-ray imaging images of the distal end state of the catheter 10 (or guide wire), respectively.
  • FIG. 31A shows a case where the determination unit 230 determines that the image is not a linear image at the time of determining a linear image (that is, for example, when the tip shape is a saddle shape).
  • FIG. 31B the straight line image and the determination unit 230 determine when the straight line image is determined.
  • the determination unit 230 determines that there is a density change when the tip density change is determined (the saddle-shaped tip part is superimposed and darkened). (When projected).
  • FIG. 31A shows a case where the determination unit 230 determines that the image is not a linear image at the time of determining a linear image (that is, for example, when the tip shape is a saddle shape).
  • the straight line image and the determination unit 230 determine when the straight line image is determined.
  • the determination unit 230 determines that there is a density change when the tip density change is determined (the saddle
  • 31C shows a case where the straight line image is determined by the determination unit 230 at the time of determining the straight line image, and further when the determination unit 230 determines that there is no density change at the time of determining the tip end density change (when the tip shape is a straight line). .
  • FIG. 32 shows the operation of the catheter tip rotation angle measuring apparatus 400. Note that steps S410, S420, S430, and S450 in FIG. 32 are the same as those in the previous embodiment, and thus description thereof will be omitted, and only different parts will be described below.
  • the first image acquisition unit 207 acquires the X-ray image 30 captured by the X-ray imaging apparatus 20 that captures an area including the distal end portion of the catheter 10 inserted into the living body tube ( Step S460 in FIG. 32).
  • the catheter tip extraction unit 208 extracts the catheter tip 15 from the X-ray image 30 acquired from the first image acquisition unit 207 (step S470 in FIG. 32).
  • the determination unit 230 determines whether or not the distal end portion of the catheter 10 is a linear image (step S471 in FIG. 32). If the determination unit 230 determines that the distal end portion of the catheter 10 is not a straight image (see FIG. 31A) (NO in step S471 in FIG. 32), the distal end shape of the catheter 10 is a saddle shape, and therefore the distal end rotation is performed.
  • the calculation of the angle is determined by the determination unit 230 (or without determination by the determination unit 230), and the catheter tip rotation angle calculation unit 205 calculates the catheter tip rotation angle (step S440 in FIG. 32).
  • the concentration change of whether there is a concentration change equal to or greater than a predetermined threshold at the distal end portion is performed by the determination unit 230 (step S472 in FIG. 32). Such a process is performed for the following reason.
  • the distal end portion of the catheter 10 when the distal end portion of the catheter 10 is a straight image, the distal end portion of the catheter 10 has a saddle shape, and the distal end portion of the catheter 10 is captured in an overlapping manner depending on the rotation angle.
  • the shape of the part is a straight line.
  • the tip part receives X-rays, a part of the tip part is attenuated twice because of the X-ray attenuation (see FIG. 31B). X1). That is, there is a part with a large attenuation rate, and if there is a part of a grayscale image that is equal to or greater than the threshold value, the tip is saddle-shaped.
  • the determination unit 230 determines that the distal end portion of the catheter 10 is a linear image, the determination unit 230 further determines whether or not there is a concentration change equal to or greater than a predetermined threshold at the distal end portion. .
  • the determination unit 230 determines that the image is a linear image and there is a change in the concentration of the distal end (YES in step S472 of FIG. 32), the distal end shape of the catheter 10 is It is a type. Therefore, the determination of the tip rotation angle is determined by the determination unit 230, and the catheter tip rotation angle calculation unit 205 calculates the catheter tip rotation angle (step S440 in FIG. 32).
  • the determination unit 230 determines that there is no change in density at the distal end as a result of determination by the determination unit 230 (NO in step S472 of FIG. 32), the distal end of the catheter 10 is determined.
  • the shape is a straight line. Therefore, the determination unit 230 determines that the calculation of the tip rotation angle is unnecessary, and the calculation of the catheter tip rotation angle is stopped by the catheter tip rotation angle calculation unit 205 (in FIG. 32, the flow ends without performing step S440). To do).
  • the distal end portion shape of the catheter 10 is a straight line and calculation of the distal end rotation angle is not necessary, useless calculation processing can be automatically omitted.
  • the catheter tip rotation angle measuring device 100 for calculating the tip rotation angle of the catheter 10 is further provided with a device for calculating the tip rotation angle of the guide wire 16. That is, as shown in FIG. 33, the catheter tip rotation angle measuring device 500 includes a first rotation angle measurement unit 102, a second rotation angle measurement unit 103, a third rotation angle measurement unit 112, and a fourth rotation angle measurement. Unit 113, guide wire insertion length measurement unit 114, guide wire tip rotation angle calculation unit 115, and target switching unit 116. Further, it further includes a movement restricting portion 111 for a guide wire. In FIG. 33, 10 is a catheter and 110 is a guide wire. Generally, the guide wire 110 is longer than the catheter 10.
  • FIG. 34 shows a block diagram of the catheter tip rotation angle measuring device 500.
  • the catheter tip rotation angle measuring apparatus 500 calculates the catheter tip rotation angle and the guide wire tip rotation angle separately.
  • the third rotation angle measurement unit 112 measures the third rotation angle ⁇ 3 of the guide wire 110 that is the rotation angle of the guide wire 110 in the third hole 1110. Similar to the first rotation angle measurement unit 102, the third rotation angle measurement unit 112 includes an imaging unit 1120 (corresponding to the imaging unit 1020) and a calculation unit 1121 (corresponding to the calculation unit 1021).
  • the fourth rotation angle measurement unit 113 measures the fourth rotation angle ⁇ 4 of the guide wire 110 that is the rotation angle of the guide wire 110 in the fourth hole 1111, similarly to the second rotation angle measurement unit 103. Similar to the second rotation angle measurement unit 103, the fourth rotation angle measurement unit 113 includes an imaging unit 1130 (corresponding to the imaging unit 1030) and a calculation unit 1131 (corresponding to the calculation unit 1031).
  • the guide wire insertion length measurement unit 114 calculates a guide wire insertion length L g that is the length from the third hole 1110 or the fourth hole 1111 to the distal end portion of the guide wire 110. taking measurement.
  • the guide wire tip rotation angle calculation unit 105 is measured by the third rotation angle ⁇ 3 measured by the third rotation angle measurement unit 112 and the fourth rotation angle measurement unit 113.
  • the guide restricting portion 111 for the guide wire penetrates from the first hole 1110 (corresponding to the first hole 1010) to the second hole 1111 (corresponding to the second hole 1011). It has a cylindrical shape and restricts the movement of the guide wire 110 in the insertion direction by passing the guide wire 110 from the first hole 1110 to the second hole 1111.
  • the target switching unit 116 performs target switching processing at a timing when either the catheter insertion length (L or L c ) or the guide wire insertion length (L g ) has changed.
  • L c > L g ⁇ L 13 (where L 13 is the distance between the movement restricting portions (101, 111) of the catheter 10 and the guide wire 110 (the first hole 1010 and the third hole 1010)).
  • the distance is set in advance if the catheter tip rotation angle is selected, and if L c ⁇ L g ⁇ L 13 , the guide wire tip is selected.
  • the rotation angle is selected, and the transmission information (presentation target and its tip rotation angle) to the presentation unit 106 is switched.
  • FIG. 35 shows a flowchart of the operation of the catheter tip rotation angle measuring apparatus 500.
  • the first rotation angle measurement part 102 measures the first rotation angle theta 1 is a rotation angle of the catheter 10 in the first bore 1010. Then, in step S420, the second rotation angle measurement part 103 measures the second rotation angle theta 2 is a rotational angle of the catheter 10 in the second hole 1011.
  • the catheter insertion length measurement unit 104 measures the catheter insertion length, which is the length from the first hole 1010 to the distal end portion 15 of the catheter 10.
  • step S440 the catheter tip rotation angle calculation unit 105, a second rotation angle first rotation angle theta 1 and measured from the first rotation angle measurement part 102 measured from the second rotation angle measurement part 103 theta 2 angle difference (theta 1 - [theta] 2) and, using the distance L 12 of the catheter insertion length was measured from a catheter insertion length measurement unit 104 L g, and a first hole 1010 and the second hole 1011, a catheter The rotation angle ⁇ c of the tip of is calculated. Steps S411 to S441 are performed simultaneously or separately with these steps S410 to S440.
  • step S411 the third rotation angle measurement unit 112 measures the third rotation angle ⁇ 3 that is the rotation angle of the guide wire 110 in the third hole 1110.
  • step S421 the fourth rotation angle measurement part 113 measures the fourth rotation angle theta 4 is a rotational angle of the guide wire 110 in the fourth hole 1111.
  • step S431 the guide wire insertion length measurement unit 114 measures a guide wire insertion length L g that is the length from the third hole 1110 to the distal end portion of the guide wire 110.
  • step S441 the guide wire tip rotation angle calculation unit 115 calculates the third rotation angle ⁇ 3 measured by the third rotation angle measurement unit 112 and the fourth rotation angle ⁇ measured by the fourth rotation angle measurement unit 113. 4 ( ⁇ 3 ⁇ 4 ), the guide wire insertion length L g measured by the guide wire insertion length measurement unit 114, and the distance L 34 between the third hole 1110 and the fourth hole 1111.
  • the rotation angle ⁇ g of the tip of the guide wire 110 is calculated.
  • step S442 the target switching section 116 acquires the catheter insertion length L c measured in step S430 from the catheter insertion length measuring unit 104, the guide wire measured guidewire insertion length L g at step S431 Obtained from the insertion length measurement unit 114.
  • the target switching unit 116 acquires the rotation angle ⁇ c of the distal end of the catheter 10 calculated in step S440 from the catheter distal end rotation angle calculation unit 105, and also calculates the rotation angle of the distal end of the guide wire 110 calculated in step S441.
  • ⁇ g is acquired from the guide wire tip rotation angle calculation unit 115. Then, the target switching unit 116 determines whether or not L c > L g ⁇ L 13 is satisfied.
  • step S441 If L c > L g ⁇ L 13 is satisfied (YES in step S441), the catheter tip rotation angle is selected, and the process proceeds to step S450. If L c > L g ⁇ L 13 is not satisfied (NO in step S441), the guide wire tip rotation angle is selected, and the process proceeds to step S450.
  • step S450 the presentation unit 106 acquires the tip rotation angle selected in step S443 or S444 from the target switching unit 116, and presents the acquired tip rotation angle to the user.
  • the distal end rotation angle of the catheter 10 and the distal end rotation angle of the guide wire 110 are separately accurate. It can be measured well. Therefore, for example, since only the guide wire 110 is used up to the vicinity of the diagnosis site or treatment site, the tip rotation angle of the guide wire 110 is accurately measured by the catheter tip rotation angle measuring device, and the diagnosis site or Since the insertion of the guide wire 110 is stopped near the treatment site and the catheter 10 is inserted along the guide wire 110 instead, the distal end rotation angle of the catheter 10 is accurately measured by the catheter distal end rotation angle measuring device. It becomes possible to do.
  • the catheter tip rotation angle measuring apparatus 600 includes a correction switching information storage unit 320 and a correction switching unit 321 as a database, for example, the catheter tip rotation angle measuring apparatus 500 of the second modification and FIG.
  • the correction switching information storage unit 320 is configured to be provided in an apparatus that is combined with the catheter tip rotation angle measuring apparatus 300 in the third embodiment, and stores information (correction switching information) necessary for correction switching.
  • the correction switching information storage unit 320 includes type information as to whether it is a catheter or a guide wire, a predetermined threshold for determining a difference in insertion length, an insertion length to a lesion site, a diagnosis site, or a treatment site (for example, a brain or heart), Information necessary for correction switching is stored.
  • the correction switching unit 321 switches the presence / absence of correction based on the correction switching information stored in the correction switching information storage unit 320 and the information from the target switching unit 116, and the result is converted into a catheter tip rotation angle correction unit 314 and a presentation unit 306. And output. Therefore, for example, the correction switching unit 321 acquires the value of the distal end rotation angle, its type (catheter or guide wire) and the corresponding insertion length from the target switching unit 116, and acquires the correction switching information from the correction switching information storage unit 320. Then, correction switching is performed based on the information.
  • the correction switching unit 321 does not perform the tip rotation angle correction and notifies the presentation unit 306 of the value of the tip rotation angle.
  • the type is output.
  • the catheter tip rotation angle correction unit 314 performs the catheter tip rotation angle correction process.
  • the correction switching unit 321 determines that the catheter 10 is in use and the distal end rotation angle of the catheter 10 is calculated, the distal end rotation considering the calibration is assumed to be corrected. It outputs to the catheter tip rotation angle correction unit 314 and the presentation unit 306 so as to calculate the angle.
  • the catheter tip rotation angle correction unit 314 and the presentation unit are calculated so as to calculate the tip rotation angle without calibration (correction). Output to 306.
  • the insertion length to the brain or heart is stored in advance in the correction switching information storage unit 320, and the correction switching unit 321 determines that the insertion length has been reached, thereby determining blood vessels such as the brain or the heart.
  • the correction switching unit 321 estimates that the structure has arrived at a complicated place, it outputs to the catheter tip rotation angle correction unit 314 and the presentation unit 306 so as to calculate the tip rotation angle in consideration of calibration (correction). To do.
  • the correction switching unit 321 estimates that the lesion length has been reached by preliminarily storing the insertion length to the lesion site in the correction switching information storage unit 320 and determining by the correction switching unit 321 that the insertion length has been reached. Then, it outputs to the catheter tip rotation angle correction unit 314 and the presentation unit 306 so as to calculate the tip rotation angle in consideration of calibration (correction).
  • the tip rotation angle measurement operation can be performed efficiently and quickly with high accuracy by performing correction processing only when necessary, instead of always performing correction.
  • correction operation is performed by inputting correction switching information from the user input unit 322 to the correction switching unit 321. Can be switched at an arbitrary timing.
  • the catheter tip rotation angle measuring device, the catheter tip rotation angle measuring method, and the catheter tip rotation angle measurement program according to the present invention have a function of measuring the rotation angle of the catheter tip and are useful during catheter contrast examination and treatment. It can also be applied to uses such as catheter operation education for residents.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Algebra (AREA)
  • Human Computer Interaction (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Technology (AREA)
  • Educational Administration (AREA)
  • Business, Economics & Management (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Medicinal Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Pulmonology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Endoscopes (AREA)

Abstract

 カテーテル先端回転角測定装置は、カテーテル10の先端側の一端である第1の孔1010から他端である第2の孔1011まで貫通した筒状の形状を有し、第1の孔から第2の孔までカテーテルを通すことによりカテーテルの移動を挿入方向に規制する移動規制部101と、第1の孔におけるカテーテルの回転角である第1回転角を測定する第1回転角測定部102と、第2の孔におけるカテーテルの回転角である第2回転角を測定する第2回転角測定部103と、第1の孔からカテーテルの先端部までの長さであるカテーテル挿入長を測定するカテーテル挿入長測定部104と、第1回転角と第2回転角との角度差、カテーテル挿入長、及び第1の孔と第2の孔との距離を用いて、カテーテルの先端の回転角を算出するカテーテル先端回転角算出部105と、を備える。

Description

カテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラム
 本発明は、カテーテルの先端の回転角を測定するカテーテル先端回転角測定装置、カテーテル回転角測定方法及びカテーテル回転角測定プログラムに関する。
 近年、患者の血管等にカテーテルを挿入し、病変又は狭窄等の診断又はその治療(例えば、PCI(Percutaneous Coronary Intervention:経皮的冠動脈形成術)あるいはTAE(Transcatheter Arterial Embolisation:経カテーテル動脈塞栓術)など)が行なわれている。
 術者は、上記の診断又は治療中、カテーテルの先端のX線透視像をモニタリングしながらカテーテル操作を行なう。その際、血管の分岐等において、目的の血管の分岐方向にカテーテルの先端の向きを合わせるため、術者は、カテーテルの先端の向きを正しく把握する必要がある。
 カテーテルの向きを検出する従来技術として、術者手元のカテーテルの回転角度を検出する方法が開示されている(特許文献1参照)。また、電子内視鏡の先端部に回転検出センサを装着し、電子内視鏡の捩れによるモニタの表示画像を自動的に補正する方法が開示されている(特許文献2参照)。
特開2009-162920号公報 特開平7-246183号公報 特表2009-522016号公報
 しかしながら、特許文献1による方法では、術者手元の回転角度とカテーテルの先端の回転角度が一致しない場合、カテーテル先端の回転角を正しく把握することができない。また、特許文献2による方法では、カテーテルの先端に回転検出センサを装着する必要がある。
 そこで、本発明の目的は、上記課題を解決し、カテーテルの先端の回転角をより正確に測定するカテーテル先端回転角測定装置、カテーテル回転角測定方法及びカテーテル回転角測定プログラムを提供することである。
 上記目的を達成するために、本発明の一態様に係るカテーテル先端回転角測定装置は、生体管内に挿入されたカテーテルの先端の回転角を測定する装置であって、第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部と、前記第1の孔における前記カテーテルの回転角である第1回転角を測定する第1回転角測定部と、前記第2の孔における前記カテーテルの回転角である第2回転角を測定する第2回転角測定部と、前記第1の孔又は前記第2の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長を測定するカテーテル挿入長測定部と、前記第1回転角測定部より測定された第1回転角と前記第2回転角測定部より測定された第2回転角との角度差、前記カテーテル挿入長測定部より測定されたカテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出部と、を備えるカテーテル先端回転角測定装置。
 なお、これらの全般的又は具体的な態様は、システム、方法、集積回路、コンピュータプログラム又は記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
 本発明の一態様によれば、より正確にカテーテルの先端の回転角を測定できる。
 本発明のこれらと他の目的と特徴は、添付された図面についての好ましい実施形態に関連した次の記述から明らかになる。この図面においては、
図1は、第1実施形態におけるカテーテル先端回転角測定装置の機能構成を示すブロック図であり、 図2は、第1実施形態における移動規制部の構成を示す図であり、 図3は、第1実施形態における第1回転角測定部及び第2回転角測定部の一例を示す図であり、 図4は、第1実施形態における第1回転角測定部及び第2回転角測定部が測定する領域を示す図であり、 図5は、ラインセンサと画素の対応を示す図であり、 図6は、カテーテルの断面と撮影部と画素との関係を示す説明図であり、 図7は、ロータリーエンコーダとカテーテルとの関係を示す断面図であり、 図8は、第1実施形態におけるカテーテル挿入長測定部の構成を示す説明図であり、 図9は、湾曲部を有するカテーテルの全体図であり、 図10は、第1実施形態における第1回転角と第2回転角とカテーテル挿入長と移動規制部の長さとカテーテル先端の回転角との関係を示す図であり、 図11Aは、第1実施形態におけるユーザに提示されるカテーテルの先端の回転角の一例をディスプレイ上に表示した状態を示す図であり、 図11Bは、第1実施形態におけるユーザに提示されるカテーテルの先端の回転角の一例を、「補正中」との表示と共に、ディスプレイ上に表示した状態を示す図であり、 図12は、第1実施形態における術者のカテーテル施術の処理動作の一例を示すフローチャートであり、 図13は、第1実施形態におけるカテーテル先端回転角測定装置の処理動作の一例を示すフローチャートであり、 図14は、第2実施形態におけるカテーテル先端回転角測定装置の機能構成を示すブロック図であり、 図15は、第2実施形態におけるX線撮像装置の機能構成を示す概略図であり、 図16は、第2実施形態におけるX線撮像装置により撮像されたX線画像の一例を示す図であり、 図17は、第2実施形態におけるカテーテル先端抽出部より抽出されたカテーテルの先端の一例を示す図であり、 図18は、第2実施形態における動作検出部の処理動作の一例を示す説明図であり、 図19は、第2実施形態におけるカテーテル先端回転角測定装置の処理動作の一例を示すフローチャートであり、 図20は、第3実施形態におけるカテーテル先端回転角測定装置の機能構成を示すブロック図であり、 図21は、第3実施形態におけるカテーテルの回転軸を決定する方法の一例を示す説明図であり、 図22は、第3実施形態におけるカテーテル回転軸決定部より復元されたカテーテルの3次元形状の一例を示す図であり、 図23は、第3実施形態におけるカテーテル先端回転角補正部の機能構成を示すブロック図であり、 図24は、第3実施形態におけるキャリブレーション用画像とカテーテルの先端の回転角との関係を示すテーブルを示す図であり、 図25は、第3実施形態におけるX線画像とキャリブレーション用画像とを比較を示す説明図であり、 図26は、第3実施形態におけるカテーテルの回転軸を求める際のエピポーラ線の算出の仕方を示す説明図であり、 図27は、第3実施形態におけるカテーテル先端回転角測定装置の処理動作の一例を示すフローチャートであり、 図28は、第3実施形態におけるカテーテル先端回転角補正部の処理動作の一例を示すフローチャートであり、 図29は、第1実施形態の変形例におけるカテーテルに先行して挿入されるガイドワイヤの一例を示す説明図であり、 図30は、第1変形例におけるカテーテル先端回転角測定装置において、先端回転角算出停止の自動判定を行うときのブロック図であり、 図31Aは、カテーテル(又はガイドワイヤ)の先端状態のX線撮像画像例を示す図であり、 図31Bは、カテーテル(又はガイドワイヤ)の先端状態のX線撮像画像例を示す図であり、 図31Cは、カテーテル(又はガイドワイヤ)の先端状態のX線撮像画像例を示す図であり、 図32は、第1変形例におけるカテーテル先端回転角測定装置の動作のフローチャートであり、 図33は、第2変形例におけるカテーテル先端回転角測定装置のブロック図であり、 図34は、第2変形例におけるカテーテル先端回転角測定装置のブロック図であり、 図35は、第2変形例におけるカテーテル先端回転角測定装置の動作のフローチャートであり、 図36は、第3変形例におけるカテーテル先端回転角測定装置のブロック図である。
 (本発明の基礎となった知見)
 カテーテル又はカテーテルに先行して挿入されるガイドワイヤには、診断部位又は治療部位まで挿入がスムーズに行なえるよう、多くの種類があり、用途に応じて様々な形態又は機能がある。例えば、右冠動脈用又は左冠動脈用のジャドキンスと呼ばれるカテーテルがあり、それぞれの解剖学的位置に合わせて先端が曲げられることにより、診断部位又は治療部位に挿入しやすい形状となっている。
 また、術者は、血管等の分岐部において、目的の血管等にカテーテル又はガイドワイヤを挿入するため、カテーテルの先端又はガイドワイヤの先端を目的の血管等の分岐方向に合わせる操作を行なう。その際、カテーテル等の先端を目的の血管等の分岐方向に合わせるためには、術者は、カテーテル等の先端の向きを正確に把握しておく必要がある。
 しかし、1方向からX線撮影したカテーテルのX線透視像をモニタに映し出す従来の方法では、X線透視像が2次元画像となるため、カテーテル等の先端の向き、特に奥行き方向を把握し辛い欠点がある。つまり、カテーテル等の先端がモニタ前方向又はモニタ後方向いずれを向いているのか、術者は把握し辛く、スムーズなカテーテル等の操作に支障をきたすことがある。
 一方、X線発生器とX線検出器とを2組有するバイプレーンタイプと呼ばれるX線透視装置を用いて、2方向からカテーテル等をX線撮影する方法がある。しかし、バイプレーンタイプによるX線撮影は、通常90度の角度差にて撮影しているため、撮影された90度の角度差のある2つのX線透視像を交互に観察しながら、かつ各々のX線発生器とX線検出器の撮影角度を考慮しつつ、立体的なカテーテル先端の形状を瞬時にかつ直感的に把握することは術者にとって大きな負担である。
 そこで、カテーテルの先端の向きを検出する従来の方法として、特許文献1及び特許文献2の方法が提案されている。
 特許文献1は、カテーテルの3次元の動きをシミュレーションする装置において、ロータリーエンコーダ等を用いて術者の手元におけるカテーテルの回転角を測定する方法について開示している。しかし、カテーテルが血管の内壁又は分岐部に接触した状態で、術者がカテーテルを挿入方向又は挿入方向を軸とした回転方向に動かした場合、カテーテルに捩れが生じることがある。捩れが生じた場合、術者の手元のカテーテルの回転角とカテーテルの先端の回転角とは一致しなくなる。そのため、特許文献1の方法では、カテーテルに捩れが生じた場合、カテーテルの先端の向きを正確に検出できない。
 また、特許文献2は、電子内視鏡システムにおいて、電子内視鏡の先端部に回転検出センサを装着し、電子内視鏡の捩れによるモニタ表示画像の軸の回転を補正する方法について開示している。しかし、直径が数mmの血管へ挿入するカテーテルの先端に回転検出センサを装着するため、カテーテル操作中に回転検出センサが血管に接触することにより、血管を傷つける危険性がある。
 本発明は、前記従来の課題を解決するもので、カテーテルの先端に回転検出センサを装着することなく、捩れが発生しているカテーテル先端の回転角を測定する装置を提供することを目的とする。
 本発明の記述を続ける前に、添付図面において同じ部品については同じ参照符号を付している。
 以下、図面を参照して本発明における実施形態を詳細に説明する前に、本発明の種々の態様について説明する。
 本発明の一態様に係るカテーテル先端回転角測定装置は、生体管内に挿入されたカテーテルの先端の回転角を測定する装置であって、第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部と、前記第1の孔における前記カテーテルの回転角である第1回転角を測定する第1回転角測定部と、前記第2の孔における前記カテーテルの回転角である第2回転角を測定する第2回転角測定部と、前記第1の孔又は前記第2の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長を測定するカテーテル挿入長測定部と、前記第1回転角測定部より測定された第1回転角と前記第2回転角測定部より測定された第2回転角との角度差、前記カテーテル挿入長測定部より測定されたカテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出部と、を備えるカテーテル先端回転角測定装置。
 この構成によれば、カテーテルの先端に回転検出センサ等を装着することなく、カテーテルの先端の回転角を測定することができる。これにより、カテーテルの先端に回転検出センサ等を装着する必要がないため、繊細かつ微妙な操作を求められるカテーテル操作において、血管を傷つける可能性を低減できる。
 例えば、さらに、前記カテーテル先端回転角測定装置は、前記カテーテル先端回転角算出部より算出された前記カテーテルの先端の回転角をユーザに提示する提示部と、を備えてもよい。
 この構成によれば、測定されたカテーテルの先端の回転角をユーザに提示することができる。したがって、術者は、測定されたカテーテルの先端の回転角からカテーテルの先端の向きを正確に把握することが出来る。これにより、術者は、血管の分岐部において、カテーテルの先端を目的の血管等の分岐方向に合わせることが容易となり、スムーズなカテーテル操作を行なうことが出来る。また、スムーズなカテーテル操作により術時間が短縮されるため、患者に照射されるX線の照射時間が短縮される。そのため、患者へのX線照射による被爆量を低減できる。
 例えば、前記カテーテル先端回転角算出部は、前記第1回転角測定部より測定される第1回転角をθとし、前記第2回転角測定部より測定される第2回転角をθとし、前記第1の孔から測定されるカテーテル挿入長をLとし、前記第1の孔と前記第2の孔との距離をL12とした場合、カテーテルの先端の回転角であるθをθ=θ-(θ-θ)・L/L12より算出するとしてもよい。
 この構成によれば、カテーテルに一様な捩れが生じている場合、カテーテル先端回転角算出部は、第1回転角とカテーテルの先端の回転角との差(θ-θ)と、第1回転角と第2回転角との角度差(θ-θ)との関係が、カテーテル挿入長(L)、第1の孔と第2の孔との距離(L12)との関係に比例していることを用いて、カテーテルの先端の回転角を測定できる。
 例えば、前記カテーテル先端回転角測定装置は、さらに、生体管内に挿入された前記カテーテルの先端部を含む領域を撮像するX線撮像装置により撮像されたX線画像を取得する第1画像取得部と、前記第1画像取得部より取得されたX線画像から前記カテーテルの先端部を抽出するカテーテル先端部抽出部と、前記カテーテル先端部抽出部より抽出された前記カテーテルの先端部の動作を検出する動作検出部と、を備え、前記カテーテル先端回転角算出部は、前記動作検出部が動作を検出した場合、前記カテーテルの先端の回転角を算出し、前記動作検出部が動作を検出しなかった場合、前記カテーテルの先端の回転角を算出しないとしてもよい。
 この構成によれば、カテーテルの先端の動きが検出されると、カテーテルの先端の回転角が測定されることになる。したがって、術者は、自身のカテーテル操作に合わせて、リアルタイムにカテーテルの先端の回転角を把握することができる。
 例えば、前記カテーテル先端回転角測定装置は、さらに、撮像時の各X線撮像装置の撮像条件及び姿勢に関する情報である撮像パラメータ情報を取得する撮像パラメータ情報取得部と、異なる位置に配される複数のX線撮像装置により撮像された前記カテーテルの先端部を含むX線画像を複数取得する第2画像取得部と、前記第2画像取得部より取得した複数のX線画像より各X線画像における前記カテーテルの直線の部分を抽出し、抽出された前記カテーテルの直線の部分と、前記撮像パラメータ情報取得部より取得した前記各X線画像が撮像された際の撮像パラメータ情報とを用いて、前記カテーテルの3次元形状を復元し、復元された前記カテーテルの3次元形状より前記カテーテルの挿入方向に対する回転軸であるカテーテル回転軸を決定するカテーテル回転軸決定部と、複数種類のカテーテルについて、前記カテーテルの先端の回転角毎に前記カテーテルの先端の形状を記憶するカテーテル形状記憶部と、カテーテル先端回転角補正部と、を備え、前記カテーテル先端回転角補正部は、前記第2画像取得部より前記複数のX線画像の中から第1X線画像を、前記カテーテル先端回転角算出部より前記カテーテルの先端の回転角及び前記第2回転角を、前記カテーテル回転軸決定部より前記カテーテル回転軸を、撮像パラメータ情報取得部より前記第1X線画像を撮影した際の撮像パラメータ情報をそれぞれ取得する取得部と、前記カテーテル形状記憶部より、前記取得部より取得した前記カテーテルの先端の回転角から前記第2回転角までの回転角毎に、前記カテーテルの先端の形状を複数取得し、取得した前記カテーテルの先端の形状と前記回転角とが対応付けられたキャリブレーション用画像を複数生成するキャリブレーション用画像生成部と、前記取得部より取得した前記第1X線画像と、前記キャリブレーション用画像生成部より生成された複数のキャリブレーション用画像とを比較し、前記第1X線画像と最も類似度の大きいキャリブレーション用画像を決定する比較部と、前記比較部より決定されたキャリブレーション用画像に対応する前記回転角を、前記カテーテルの先端の回転角と決定する回転角度決定部とを有し、前記提示部は、前記回転角度決定部より決定された前記カテーテルの先端の回転角をユーザに提示するとしてもよい。
 この構成によれば、カテーテルの一部に捩れが生じている場合においても、カテーテルの先端の回転角を正確に測定することができる。
 以下本発明の実施の形態について、図面を参照しながら説明する。なお、以下で説明する実施の形態は、いずれも本発明の好ましい一具体例を示すものである。以下の実施の形態で示される数値、構成要素、構成要素の接続形態、ステップ、ステップの順序などは、一例であり、本発明を限定する主旨ではない。本発明は、請求の範囲だけによって限定される。よって、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素については、本発明の課題を達成するのに必ずしも必要ではないが、より好ましい形態を構成するものとして説明される。
 (第1実施形態)
 図1に、第1実施形態におけるカテーテル先端回転角測定装置100の機能構成を示す。
 図1に示すように、カテーテル先端回転角測定装置100は、第1回転角測定部102と、第2回転角測定部103と、カテーテル挿入長測定部104と、カテーテル先端回転角算出部105とを有する。また、図示せぬ移動規制部101を有する。以下に、カテーテル先端回転角測定装置100が有する各構成要素について説明する。
 <構成>
 <移動規制部101>
 移動規制部101は、第1の孔1010から第2の孔1011まで貫通した筒状の形状を有し、第1の孔1010から第2の孔1011までカテーテル10を通すことによりカテーテル10の移動を挿入方向に規制するものである。
 図2に、移動規制部101の概略の一例を示す。移動規制部101は、第1の孔1010から第2の孔1011まで貫通した筒状の形状である。また、第1の孔1010は、挿入方向からみて移動規制部101の後端であり、第2の孔1011は、挿入方向からみて移動規制部101の前端である。つまり、第1の孔1010は術者側の孔であり、第2の孔1011は患者側の孔である。
 カテーテル10を第1の孔1010から第2の孔1011まで通すことにより、カテーテル10の移動は挿入方向に規制される。ここで、挿入方向とは、カテーテル10の軸方向に加え、軸方向回りの回転も含む。つまり、カテーテル10の移動が挿入方向に規制されるとは、カテーテル10は、カテーテル10の軸方向と直交する半径方向へ移動できないことを意味する。
 第1の孔1010及び第2の孔1011は、カテーテル10が貫通可能な大きさであり、カテーテル10の直径より所定の大きさだけ大きい。第1の孔1010及び第2の孔1011は、カテーテル10の直径とほぼ同じ大きさであり、カテーテル10を挿入した状態では、カテーテル10の軸方向と直交する半径方向には、ほとんど隙間が無い。第1の孔1010及び第2の孔1011の大きさは、例えば、カテーテル10の直径より1~3mm程度大きい。なお、第1の孔1010及び第2の孔1011の大きさは、カテーテル10の直径又は種類により許容される大きさが異なっても良い。つまり、第1の孔1010及び第2の孔1011は、カテーテル10が通常のカテーテルより太いものであれば、カテーテル10の直径より3mm程度大きく、カテーテル10が通常のカテーテルより細いものであれば、カテーテル10の直径より1mm程度大きくあればよい。
 また、筒状の内部表面は、摩擦抵抗が少ないことが望ましい。これにより、カテーテル操作をスムーズに行なうことができ、かつ移動規制部101の内部にてカテーテル10の捩れが生じることを防ぐことができる。
 第1の孔1010から第2の孔1011までの距離L12は、移動規制部101の長さと等しい。また、後述するが、移動規制部101の長さは長いほど、カテーテル10の先端の回転角を高精度に測定することができる。しかし、移動規制部101の長さが長いほど、カテーテル操作は難しくなるデメリットもある。ここでは、例えば、移動規制部101の長さは10cm程度とする。
 なお、移動規制部101は、脚1012によって、台1013に固定されていてもよい。移動規制部101が台1013に固定されることにより、後述の第1回転角及び第2回転角を高精度に測定することができる。
 <第1回転角測定部102>
 第1回転角測定部102は、第1の孔1010におけるカテーテル10の回転角である第1回転角θを測定する。
 第1回転角測定部102の測定箇所は、第1の孔1010の位置を含む所定の範囲内である。第1回転角測定部102の測定箇所は、第1の孔1010に近いほど望ましく、例えば、第1の孔より3mm以内の範囲である。
 図3に、第1回転角測定部102の概略の一例を示す。第1回転角測定部102は、撮影部1020及び計算部1021より構成される。図4に、移動規制部101を上から下にみたものを示す。カテーテル10の表面には予めカテーテル10の中心軸に平行な方向に線11が引かれている。撮影部1020は、範囲1022を撮影する。計算部1021は、撮影部1020が撮影した画像から線11の位置を抽出し、カテーテル10の中心に対する線11の位置変位より、範囲1022における回転角を取得する。以下に、第1回転角測定部102が第1回転角を測定する具体的な方法について説明する。
 撮影部1020は例えばカメラであり、計算部1021は例えばPCのCPUである。撮影部1020は、術者からの検出スタートの入力を受けて撮影を開始する。撮影部1020は一定間隔でカテーテル10を撮影し続ける。
 計算部1021は、撮影された複数の画像について、各画像における線11の位置をそれぞれ抽出する。抽出は、例えば、カテーテル10が白で、線11が黒の場合、黒の位置を抽出する。図5に、撮影部1020が15画素からなるラインセンサの場合、撮影される画像を示す。ここでは、簡略化のため、カテーテル10の(軸方向に垂直な向きである)幅全体がこの15画素に撮影されるように撮影部1020が配置されているものとして説明する。図5に示すように、各画素には上から下に順番に番号(n=0,1,2,…,14)が付与されている。図5の場合では、計算部1021は、画素5の番号である「5」を取得する処理を行う。具体的には、計算部1021は各画素の明度が閾値以下か否かを判定する。例えば、画素の明度は0~255の値を取り、0が最も暗く、255が最も明るいとして、閾値を50とする。明度が閾値以下であった画素の番号を線11の位置とする。
 次に、計算部1021は、抽出された線11の位置から範囲1022における回転角を測定する。図6に、線11の位置とカテーテルの回転角度との関係を示す。図6の右側に示す円は、カテーテル10の断面図である。図6の左側は撮影部1020を示している。撮影部1020には画素0~14までを撮影する素子があり、これを拡大したものを図6の真ん中に示している。また、各画素に撮影される範囲の境界を点線で示している。画素7は、撮影部1020が撮影した画像におけるカテーテルの範囲(ここでは画素0から画素14)の中心が撮影される画素である。以降、画素7を中心画素1022と呼ぶ。中心画素1022に撮影されるカテーテル10の表面部分を中心位置12とする。半径rは、撮影部1020が撮影した画像におけるカテーテル10の半径の長さである。半径rは、図6では、中心画素1022から、カテーテル10の端が撮影される画素14までの長さと等しく、14-7=7より7画素となる。画素5は、カテーテル10の表面上に引かれている線11の位置である。以下では、画素5を線位置画素1023と呼ぶ。距離hは、線位置画素1023から中心画素1022までの画素上での距離であり、ここでは7-5=2である。角θは、中心位置12からの線11の回転角、つまり第1の孔1010におけるカテーテル10の回転角である第1回転角である。
 角θは、三角関数の定義より、(1)式の関係式で表される。
  sinθ=h/r     ・・・・・(1)式
 よって、角θは、(1)式より(2)式となる。
  θ=arcsin(h/r)     ・・・・・(2)式
 ここで、関数arcsinはsin関数の逆関数である。よって、図6の場合では、θ=arcsin(2/7)となり、これにより、第1の孔1010におけるカテーテル10の回転角である第1回転角θを測定することができる。
 なお、第1回転角測定部102は、撮影部1020と計算部1021とで構成されている必要はなく、ロータリーエンコーダ等より構成されていても良い。図7に、ロータリーエンコーダを用いてカテーテル10の回転角を測定する際の第1回転角測定部102の構成図を示す。図7は、カテーテル10の挿入方向から見た様子を示している。カテーテル10に予め施されている溝部13と、移動規制部101の内部に予め施されている突起部1014とが噛み合わされることにより、カテーテル10の回転とロータリーエンコーダの回転とが同期し、カテーテル10の回転角が測定できる。ここで、ロータリーエンコーダとして、絶対値形エンコーダ等(小川鑛一ら、基礎ロボット工学、1998、p98、東京電機大学出版局)を用いることができる。
 <第2回転角測定部103>
 第2回転角測定部103は、第2の孔1011におけるカテーテル10の回転角である第2回転角θを測定する。
 第2回転角測定部103は、第1回転角測定部102と同様の構成により実現できる。図3に、第2回転角測定部103の概略の一例を示す。第2回転角測定部103は、撮影部1030及び計算部1031より構成される。また、撮影部1030と計算部1031とを用いて、第2回転角測定部103が第2回転角θを測定する具体的な方法は、第1回転角測定部102が第1回転角θを測定する方法と同様である。
 なお、第2回転角測定部103は、撮影部1030と計算部1031とで構成されている必要はなく、ロータリーエンコーダ等より構成されていても良い。また、ロータリーエンコーダを用いて、第2回転角測定部103が第2回転角θを測定する具体的な方法は、第1回転角測定部102が第1回転角θを測定する方法と同様である。
 <カテーテル挿入長測定部104>
 カテーテル挿入長測定部104は、第1の孔1010又は第2の孔1011からカテーテル10の先端部までの長さであるカテーテル挿入長Lを測定する。
 図8に、カテーテル挿入長測定部104の概略の一例を示す。図8は、第1の孔1010からカテーテル10の先端部15までの長さをカテーテル挿入長Lとしている場合である。カテーテル挿入長測定部104は、例えば、カテーテル10を押圧することにより、カテーテル10の進退に連動するローラ1040と、ローラ1040の回転量を測定するエンコーダ1041と、エンコーダ1041の回転角の値の変化量から挿入長を算出する計算部1042とからなる。
 計算部1042は、カテーテル10の操作を開始する前に初期値「0」に設定しておく。カテーテル10の操作を開始する際、カテーテル10をローラ1040に通し、その時点より、計算部1042は、カテーテル10の挿入長の算出を開始し、カテーテル10の操作を終了するまでの間、算出し続ける。これにより、カテーテル挿入長測定部104は、カテーテル10の操作中、リアルタイムにカテーテル10の挿入長Lを測定し続けることができる。
 また、図9に示すように、カテーテル10の先端は湾曲部15を有している。湾曲部15(以後、「カテーテルの先端部15」ともいう。)は、その他のカテーテル10の部分と材質等が異なり、曲がり易くなっている。カテーテル挿入長測定部104は、第1の孔1010から当該湾曲部15までの距離をカテーテル挿入長Lとして測定する。つまり、カテーテルの先端部15を含まない長さが、カテーテル挿入長Lに相当する。カテーテル挿入長Lに、カテーテルの先端部15を含めないことにより、カテーテル10の先端の回転角θを、より正確に算出することができる。
 なお、カテーテル挿入長測定部104は、必ずしもカテーテルの先端部15を含まない長さをカテーテル挿入長Lとする必要は無く、カテーテルの先端部15を含む長さをカテーテル挿入長Lとして測定してもよい。
 なお、カテーテル挿入長測定部104は、カテーテル10の表面に記載された挿入位置に関する目盛りをカメラ等の撮影手段により撮影し、当該目盛りの値に基づいて挿入長を算出する構成であっても良い。
 なお、カテーテル挿入長測定部104は、移動規制部101と独立して設けられている必要な無く、移動規制部101に設けられていてもよい。
 なお、カテーテル挿入長測定部104は、第2の孔1011からカテーテル10の先端部までの長さをカテーテル挿入長Lとしても良い。
 <カテーテル先端回転角算出部105>
 カテーテル先端回転角算出部105は、第1回転角測定部102より測定された第1回転角θと第2回転角測定部103より測定された第2回転角θとの角度差(θ-θ)、カテーテル挿入長測定部104より測定されたカテーテル挿入長L、及び第1の孔1010と第2の孔1011との距離L12を用いて、カテーテルの先端の回転角θを算出する。
 カテーテル先端回転角算出部105は、第1回転角測定部102より第1回転角θを、第2回転角測定部103より第2回転角θを、カテーテル挿入長測定部104よりカテーテル挿入長Lをそれぞれ取得する。カテーテル先端回転角算出部105は、取得した第1回転角θと第2回転角θとの角度差(θ-θ)を算出する。また、術者は、カテーテル先端回転角算出部105に、第1の孔1010と第2の孔1011との距離L12を、カテーテル10の操作を開始する前に予め設定しておく。
 図10に、第1の孔1010が術者側である場合の各パラメータの関係を示す。以下、第1の孔1010が術者側であり、カテーテル挿入長Lが第1の孔1010からカテーテル10の先端部15までの長さである場合について説明する。
 捩れがカテーテル10に一様に生じている場合、各パラメータの間には(3)式の関係が成立している。
  (θ-θ)=(θ-θ)・L/L12     ・・・・・(3)式
 (3)式は、カテーテル10の全体の捩れ角度(θ-θ)及び第1の孔1010と第2の孔1011との捩れ角度(θ-θ)の関係が、カテーテル挿入長L及び第1の孔1010と第2の孔1011との距離L12の関係に比例することを表している。よって、カテーテル先端の回転角θは、(4)式より算出することができる。
  θ=θ-(θ-θ)・L/L12     ・・・・・(4)式
 なお、第1の孔1010が術者側として説明を行なってきたが、第1の孔1010が患者側であり、第2の孔1011が術者側である場合も考えられる。この場合、(3)式は(5)式のようになる。
  (θ-θ)=(θ-θ)・L/L12     ・・・・・(5)式
 また、この場合、カテーテル先端の回転角θは(6)式より算出することができる。
  θ=θ-(θ-θ)・L/L12     ・・・・・(6)式
 なお、第2の孔1011が患者側である場合、カテーテル挿入長Lが、第2の孔1011からカテーテル10の先端部15までの長さである際のカテーテル先端の回転角θは(7)式より算出することができる。
  θ=θ-(θ-θ)・(L+L12)/L12     ・・・・・(7)式
 <提示部106>
 提示部106は、カテーテル先端回転角算出部105より算出されたカテーテル10の先端の回転角θをユーザに提示する。
 提示部106は、カテーテル先端回転角測定装置100の必須構成要素ではないが、カテーテル先端回転角算出部105より算出されたカテーテル10の先端の回転角θをユーザに提示することにより、ユーザである術者のカテーテル操作に役立つ。
 提示部106は、例えば、カテーテル先端回転角算出部105より算出されたカテーテル10の先端の回転角θを、ディスプレイ上に数値等を表示する。また、カテーテル10の先端の回転角θを反映させたカテーテル10のCG等を用いて3次元的に表示してもよい。また、提示部106は、術者がカテーテル操作中にモニタリングするディスプレイ上に表示されるX線画像に、カテーテル10の先端の回転角θを重畳表示してもよい。図11Aに、カテーテル10の先端の回転角θを、矢印を用いてディスプレイの右下に表示している例を示す。
 <動作>
 次に、以上のように構成されたカテーテル先端回転角測定装置100の各構成要素の動作に関して、具体的に説明する。
 カテーテル先端回転角測定装置100の各構成要素の動作の説明の前に、カテーテル先端回転角測定装置100が利用される場面であるカテーテル施術の全体の流れについて説明する。図12に、カテーテルによる施術の全体のフローを示す。
 <ステップS100>
 術者は、施術を行なう事前準備として、第1回転角測定部102、第2回転角測定部103、及びカテーテル挿入長測定部104の初期値を0に設定する。また、術者は、カテーテル先端回転角算出部105に、第1の孔1010と第2の孔1011との距離L12を入力する。そして、術者は、診断部位又は治療部位に適したカテーテル10を選定し、選定したカテーテル10を移動規制部101の第1の孔1010から第2の孔1011まで通し、さらに、カテーテル挿入長測定部104に通す。
 なお、ここで、第1回転角測定部102及び第2回転角測定部103における「0」とは、例えば、カテーテル10の挿入方向に直交し、かつ床から天井へ向く方向を0としてもよい。
 なお、第1回転角測定部102、第2回転角測定部103、及びカテーテル挿入長測定部104の初期値を0に設定すると説明したが、0でなくてもよい。測定開始時点の数値と測定中の数値より、それぞれの値が測定しても良い。
 <ステップS200>
 術者は、患者にカテーテル10を挿入する。
 <ステップS300>
 術者は、ディスプレイ上に写し出されたカテーテル10のX線画像をモニタリングしながらカテーテル10を操作し、カテーテル10を診断部位又は治療部位まで挿入する。
 <ステップS400>
 術者は、カテーテル10の先端の向きを認識したい場合、カテーテル10の先端の回転角θの測定を開始する入力を行なう。入力を受付けるデバイスは、例えば、キーボード又はプッシュボタン等が挙げられる。カテーテル先端回転角測定装置100は、入力を受付けると、後述のステップS410からステップS450を実行し、カテーテル10の先端の回転角を術者に提示する。
 なお、カテーテル10の先端の回転角の測定は、入力を待って開始されるものに限らず、所定の時間間隔で測定しても良い。
 <ステップS500>
 術者は、提示部106より提示されたカテーテルの先端の回転角θ等を参考に、カテーテル10の操作を続行する。
 <ステップS600>
 術者は、診断又は治療を完了していない場合、ステップS300を実行し、診断又は治療を完了した場合、ステップS700を実行する。
 <ステップS700>
 術者は、診断又は治療を完了した場合、患者よりカテーテル10を抜去し、施術を終える。
 以上の動作が、カテーテル施術の全体の説明である。次に、第1実施形態におけるカテーテル先端回転角測定装置100の処理動作について説明する。
 カテーテル10の先端の回転角θの測定方法は、単にカテーテル10の先端の回転角θを測定するものである。つまり、当該方法には、術者が、ステップS400で測定されたカテーテルの先端の回転角θに基づいて、ステップS500のカテーテル操作を行なうことは含まれていない。したがって、カテーテル10の先端の回転角θの測定方法は、人体に対して外科的処理を施す方法ではなく、人間を手術する方法には該当しない。
 図13に、第1実施形態におけるカテーテル先端回転角測定装置100の処理動作のフローを示す。
 <ステップS410>
 第1回転角測定部102は、第1の孔1010におけるカテーテル10の回転角である第1回転角を測定する。
 第1回転角θの測定の具体的な方法は、前述の通りである。また、第1回転角測定部102は、測定した第1回転角θをカテーテル先端回転角算出部105へ出力する。
 <ステップS420>
 第2回転角測定部103は、第2の孔1011におけるカテーテル10の回転角である第2回転角を測定する。
 第2回転角θの測定の具体的な方法は、前述の通りである。また、第2回転角測定部103は、測定した第2回転角θをカテーテル先端回転角算出部105へ出力する。
 <ステップS430>
 カテーテル挿入長測定部104は、第1の孔1010からカテーテル10の先端部15までの長さであるカテーテル挿入長を測定する。
 カテーテル挿入長Lの測定の具体的な方法は、前述の通りである。また、カテーテル挿入長測定部104は、測定したカテーテル挿入長Lをカテーテル先端回転角算出部105へ出力する。
 <ステップS440>
 カテーテル先端回転角算出部105は、第1回転角測定部102より測定された第1回転角θと第2回転角測定部103より測定された第2回転角θとの角度差(θ-θ)、カテーテル挿入長測定部104より測定されたカテーテル挿入長L、及び第1の孔1010と第2の孔1011との距離L12を用いて、カテーテルの先端の回転角θを算出する。
 始めに、カテーテル先端回転角算出部105は、ステップS410で測定された第1回転角θを、ステップS420で測定された第2回転角θを、ステップS430で測定されたカテーテル挿入長Lをそれぞれ取得する。
 次に、カテーテル先端回転角算出部105は、取得した第1回転角θ及び第2回転角θより、第1回転角θと第2回転角θとの角度差(θ-θ)を計算する。
 次に、カテーテル先端回転角算出部105は、計算した第1回転角θと第2回転角θとの角度差(θ-θ)と、取得したカテーテル挿入長Lと、ステップS100で設定された第1の孔1010と第2の孔1011との距離L12とを用いて、(4)式より、カテーテル10の先端の回転角θを算出する。
 カテーテル先端回転角算出部105は、算出したカテーテル10の先端の回転角θを、提示部106へ出力する。
 <ステップS450>
 提示部106は、ステップS440で算出されたカテーテル10の先端の回転角θを取得し、取得したカテーテル10の先端の回転角θをユーザに提示する。
 提示部106による提示の具体的な方法は、前述の通りである。
<効果>
 以上のように、第1実施形態におけるカテーテル先端回転角測定装置100によれば、カテーテル10に捩れが一様に生じている場合、カテーテル10の先端に回転検出センサ等を装着することなく、カテーテル10の先端の回転角を測定することができる。これにより、カテーテル10の先端に回転検出センサ等を装着する必要がないため、繊細かつ微妙な操作を求められるカテーテル10の操作において、血管を傷つける可能性を低減できる。また、測定されたカテーテル10の先端の回転角θをユーザに提示することにより、術者は、測定されたカテーテル10の先端の回転角θからカテーテル10の先端の向きを正確に把握することが出来る。これにより、術者は、血管の分岐部において、カテーテル10の先端を目的の血管等の分岐方向に合わせることが容易となり、スムーズなカテーテル操作を行なうことが出来る。また、スムーズなカテーテル操作により術時間が短縮されるため、患者に照射されるX線の照射時間が短縮される。そのため、患者へのX線照射による被爆量を低減できる。
 なお、本第1実施形態は、カテーテルに先行して挿入されるガイドワイヤに対しても適用される。図29に血管に挿入したカテーテル10の内腔にガイドワイヤ16を通して、ガイドワイヤ16をカテーテル10に先行して進入させている際の一画像例を示す。血管カテーテル治療においては、一般的に、カテーテル10を挿入する前にガイドワイヤ16を先行させることにより、血管壁を傷つけにくくし、またカテーテル10を病変部までスムーズに運べるようにしている。また、カテーテル10の内腔の表面は摩擦抵抗を小さくしていたり、あるいは治療中に術者がカテーテル10の内腔に少量の水を注入することによって、ガイドワイヤを通りやすくしている。そのため、ガイドワイヤ16に捩れが発生する場合は、ガイドワイヤ16の先端の回転に対して負荷がかかり、それによるガイドワイヤ16の捩れは、カテーテル10の内腔において一様に生じていると考えられる。
 また、前述の関係は、親カテーテルと子カテーテルとの間にも成立し得る。ここで、親カテーテルとは、後者の子カテーテルを誘導しやすいようにするもので、太さは2mm前後であり、子カテーテルは、実際に病変部を治療するためのもので、太さは1mm以下の場合が多い。血管に挿入した親カテーテルの内腔に子カテーテルを通すことで、子カテーテルを病変部までスムーズに運ぶことができて、子カテーテルにより、狭窄した血管をバルーンで広げるなどの病変部の治療を行うことができる。この場合においても、子カテーテルに捩れが発生する場合は、子カテーテルの先端の回転に対して負荷がかかり、それによる子カテーテルの捩れは、親カテーテルの内腔において一様に生じていると考えられる。
 従って、本第1実施形態におけるカテーテル先端回転角測定装置100は、カテーテルに先行して挿入されるガイドワイヤ、又は、親カテーテルの先へ挿入される子カテーテルにおいて、特に効果的に用いることができる。
 (第2実施形態)
 第2実施形態におけるカテーテル先端回転角測定装置200は、カテーテル10の先端部15の動きを検出し、動きが検出された場合、カテーテル10の先端の回転角θを算出し、動きが検出されなかった場合、カテーテル10の先端の回転角θを算出しない点が、第1実施形態のカテーテル先端回転角測定装置100と異なる。以下に、第1実施形態と異なる点を中心に説明する。
 図14に、第2実施形態におけるカテーテル先端回転角測定装置200の機能構成を示す。なお、図14において、図1と同様の構成要素については同一の符号を付し、適宜説明を省略する。
 図14に示すように、カテーテル先端回転角測定装置200は、第1回転角測定部102と、第2回転角測定部103と、カテーテル挿入長測定部104と、カテーテル先端回転角算出部205と、提示部106と、第1画像取得部207と、カテーテル先端部抽出部208と、動作検出部209とを有する。また、図14には図示しないが、移動規制部101を有する。以下に、カテーテル先端回転角測定装置200が有する各構成要素について説明する。
 <構成>
 <第1画像取得部207>
 第1画像取得部207は、生体管内に挿入されたカテーテル10の先端部を含む領域を撮像するX線撮像装置20により撮像されたX線画像30を取得する。
 図15に、第1画像取得部207とX線撮像装置20との構成を示す。ここで、X線撮像装置20とは、患者の撮影対象部位に対してX線を照射して撮影されたX線画像又は造影剤を注入した際に撮影された血管画像を取得する装置であり、例えば、X線血管造影装置又はアンギオグラフィと称される。
 X線撮像装置20は、X線発生部21と、X線検出部22と、機構部23と、機構制御部24と、アーム26と、データベース25とを有する。X線発生部21は、高電圧を用いてX線を発生するX線管とX線の一部を遮蔽することによって照射野を制御するX線絞り器とを有しており、寝台27上の患者28にX線を照射する。X線検出部22は、患者28を透過したX線を受けて画像情報を記録し、記録した画像情報を出力するものであり、例えば、X線観応層を配置し、X線をデジタルデータに変換して出力するFPD(Flat Panel Detector)として構成されている。また、X線検出部22は、X線が照射されると、照射されたX線画像30を示す画像情報をデータベース25へ出力する。機構部23は、術者の操作指示を受け付けた機構制御部24の指示に基づいてアーム26及び寝台27を移動させる。第1画像取得部207は、データベース25に格納されたX線画像又はX線映像を読み出し、カテーテル先端部抽出部208へ出力する。
 図16に、第1画像取得部207が取得するカテーテル10のX線画像30の一例を示す。カテーテルの先端部15とは、図9に示した湾曲部15を含む部分であり、例えば、カテーテル10の先端から湾曲部15の長さLaに加えて1~2cmまでの部分をいう。なお、カテーテルの先端部15は、湾曲部15と同一の部分であってもよい。
 <カテーテル先端部抽出部208>
 カテーテル先端部抽出部208は、第1画像取得部207より取得されたX線画像30からカテーテルの先端部15を抽出する。
 以下に、カテーテル先端部抽出部208がカテーテルの先端部15を抽出する具体的な方法について説明する。カテーテル先端部抽出部208は、取得したX線画像30と前フレームのX線画像とを比較し、画素ごとに明度が暗くなる方へ変化した部分のみを抽出する。これにより、カテーテル10に動きがあった場合のみ、カテーテルの先端部15を抽出できる。例えば、画素の明度は「0~255」の値を取り、「0」が最も暗く、「255」が最も明るいとした場合、X線画像30の前フレームの画素から現フレームの画素を引いた値が「50」以上の画素をカテーテル10の部分として「0」とし、それ以外の画素をカテーテル10以外の部分として「255」とする。これにより、2値化され、カテーテル10の先端部15を抽出できる。図17に、カテーテル先端部抽出部208より抽出されたカテーテル10の先端部15の画像を示す。
 なお、カテーテル先端部抽出部208は、一定の時間前(例えば1秒前)のフレームと比較しても良い。
 なお、カテーテル先端部抽出部208は、上記の方法に限らず、カテーテル10の直線度を判定し、直線度の値が低い部分を湾曲部15として認識し、認識した湾曲部15を含む所定の領域をカテーテル10の先端部15として抽出しても良い。
 <動作検出部209>
 動作検出部209は、カテーテル先端部抽出部208より抽出されたカテーテル10の先端部15の動作を検出する。動作検出部209は、カテーテル10の先端部15の動作を検出することにより、手元部の挿入又は回転の動作を間接的に検出することができる。ここで、手元部とは、移動規制部101及びカテーテル挿入長測定部104を含めた部分である。
 以下に、動作検出部209が先端部15の動作を検出する具体的な方法について説明する。動作検出部209は、カテーテル10の先端部15の抽出において、明度に変化のあった画素数が所定の画素数以上であった場合、先端部15の動作を検出する。例えば、動作検出部209は、明度が「0」であり画素数が「10」以上の場合、先端部15の動作を検出する。また、動作検出部209は、カテーテル10の先端部15を抽出において、カテーテル10の先端部15に動きが無い場合、全ての画素の明度が「255」として抽出され、先端部15の動作が検出されない。図18に、前フレームと現フレームとを比較することにより、動作検出部209が検出するカテーテル10の先端部15の動作の有無を示す。
 なお、動作検出部209は、上記の方法に限らず、X線画像30から辺縁描出等された先端部15を、所定の時間間隔(例えば、1ms間隔)に比較し、所定の値(例えば、0.03mm)以上の変化があった場合に動作を検出し、所定の値以上の変化がなかった場合に動作を検出しないとしてもよい。
 <カテーテル先端回転角算出部205>
 カテーテル先端回転角算出部205は、第1回転角測定部102より測定された第1回転角θと第2回転角測定部103より測定された第2回転角θとの角度差(θ-θ)、カテーテル挿入長測定部104より測定されたカテーテル挿入長L、及び第1の孔1010と第2の孔1011との距離L12を用いて、カテーテルの先端の回転角θを算出する。
 また、カテーテル先端回転角算出部205は、動作検出部209が動作を検出した場合、カテーテル10の先端の回転角θを算出し、動作検出部209が動作を検出しなかった場合、カテーテル10の先端の回転角θを算出しない。
 なお、カテーテル先端回転角算出部205が算出するカテーテル10の先端の回転角θの具体的な算出方法は、第1実施形態におけるカテーテル先端回転角算出部105の算出方法と同様である。
 <動作>
 次に、以上のように構成されたカテーテル先端回転角測定装置200の各構成要素の動作に関して具体的に説明する。図19に、第2実施形態におけるカテーテル先端回転角測定装置200の処理動作のフローを示す。なお、図19において、図13と同様のステップについては同一の符号を付し、説明を適宜省略する。
 <ステップS460>
 第1画像取得部207は、X線撮像装置20より撮像されたX線画像30をデータベース25より取得する。
 画像取得の具体的な方法は、前述の通りである。また、第1画像取得部207は、取得したX線画像30をカテーテル先端部抽出部208へ出力する。
 <ステップS470>
 カテーテル先端部抽出部208は、第1画像取得部207より取得されたX線画像30からカテーテル10の先端部15を抽出する。
 始めに、カテーテル先端部抽出部208は、ステップS460より取得されたX線画像30を取得する。次に、カテーテル先端部抽出部208は、取得したX線画像30からカテーテル10の先端部15を抽出する。先端部15の具体的な抽出方法は、前述の通りである。次に、カテーテル先端部抽出部208は、抽出した先端部15を動作検出部209へ出力する。
 <ステップS480>
 動作検出部209は、カテーテル先端部抽出部208より抽出されたカテーテルの先端部15の動作を検出する。
 始めに、動作検出部209は、ステップS470より抽出された先端部15を取得する。次に、動作検出部209は、取得した先端部15の動作を検出する。先端部15の動作の具体的な検出方法は、前述の通りである。次に、動作検出部209は、先端部15の動作を検出した場合、ステップS440を実行し、先端部15の動作を検出しなかった場合、ステップS450を実行する。つまり、動作検出部209が先端部15の動作を検出しなかった場合、提示部105は、現在提示しているカテーテル10の先端の回転角θをユーザに提示し続ける。したがって、先端部15に動作が無い場合、第1の回転角θ、第2の回転角θ及びカテーテル挿入長Lに変化がある場合でも、その変化は僅かであるとして、各パラメータの変化量は「0」とみなす。
 <効果>
 以上のように、第2実施形態におけるカテーテル先端回転角測定装置200によれば、カテーテル10の先端部15の動作が検出されると、カテーテル10の先端の回転角θが測定されることになる。したがって、術者は、自身のカテーテル操作に合わせて、リアルタイムにカテーテル10の先端の回転角θを把握することができる。これにより、術者は、血管の分岐部において、カテーテル10の先端を目的の血管等の分岐方向に合わせることが容易となり、スムーズなカテーテル操作を行なうことが出来る。また、スムーズなカテーテル操作により術時間が短縮されるため、患者に照射されるX線の照射時間が短縮される。そのため、患者へのX線照射による被爆量を低減できる。
 (第3実施形態)
 第3実施形態におけるカテーテル先端回転角測定装置300は、カテーテル10の一部に捩れが生じている場合においてもカテーテル10の先端の回転角を正確に測定することができる点が、第1実施形態のカテーテル先端回転角測定装置100と異なる。以下に、第1実施形態と異なる点を中心に説明する。
 図20に、第3実施形態におけるカテーテル先端回転角測定装置300の機能構成を示す。なお、図20において、図1と同様の構成要素については同一の符号を付し、適宜説明を省略する。
 図20に示すように、カテーテル先端回転角測定装置300は、第1回転角測定部102と、第2回転角測定部103と、カテーテル挿入長測定部104と、カテーテル先端回転角算出部305と、提示部306と、撮像パラメータ情報取得部310と、第2画像取得部311と、カテーテル回転軸決定部312と、カテーテル形状記憶部313と、カテーテル先端回転角補正部314とを有する。また、図示せぬ移動規制部101を有する。以下に、カテーテル先端回転角測定装置300が有する各構成要素について説明する。
 <構成>
 <撮像パラメータ情報取得部310>
 撮像パラメータ情報取得部310は、撮像時の各X線撮像装置20の撮像条件及び姿勢に関する情報である撮像パラメータ情報を取得する。
 撮像パラメータ情報とは、X線撮像装置20の装置パラメータ(内部パラメータ)と、X線撮像装置20の撮影時の位置姿勢パラメータ(外部パラメータ)である。より、具体的には、内部パラメータとは、X線撮像装置20の撮影に関するパラメータであり、例えば、シャッター速度、絞り値、露出補正値、ISO感度値、ホワイトバランス設定値、記録画像サイズ及び記録画像画質、その他画像作りに影響するパラメータである。また、外部パラメータとは、X線撮像装置20の撮影時の位置又は向きに関するパラメータであり、例えば、空間中において原点とX、Y、Z軸の方向を定めた場合、X線画像装置20の位置は(x、y、z)で表される。
 撮像パラメータ情報取得部310は、取得した撮像パラメータ情報を透視投影行例として外部に出力する。
 ここで、透視投影行列について説明する。3次元空間中の点Xが2次元透視像中の点mに投影されているとすると、この投影は斉次座標を用いて(8)式のように表せる。
Figure JPOXMLDOC01-appb-M000001
 ここで、(8)式の画像点mは3×1ベクトル、3次元点Xは4×1ベクトルであり、3×4行列の行列Pを透視投影行列又は投影行列と呼ぶ。ここで、λは定数倍の不定性を示す実数である。投影行列Pは、透視像撮影時における焦点距離、3次元空間座標系から透視像座標系への3次元の回転及び並行移動などの幾何学的な関係を示している。
 <第2画像取得部311>
 第2画像取得部311は、生体管内に挿入されたカテーテル10の先端部15を含む領域を撮像し、異なる位置に配される複数のX線撮像装置20により撮像されたX線画像30を複数取得する。
 第2画像取得部311は、第1画像取得部207とほぼ同様の構成である。第2画像取得部311は、異なる位置に配される複数のX線撮像装置20より撮像された複数のX線画像30を取得する。ここで、異なる位置に配される複数のX線撮像装置20とは、例えば、X線発生部21とX線検出部22とを2組揃え、それらの撮影角度差が90度あり、2つのX線画像を同時に撮影できるバイプレーンタイプのX線撮像装置である。また、異なる位置とは、各X線検出部22が同じ箇所に設置される場合及び各X線検出部22が対向する箇所に設置される場合を除く意味であり、各X線検出部22の角度差が5度以上あれば良い。
 <カテーテル回転軸決定部312>
 カテーテル回転軸決定部312は、第2画像取得部311より取得した複数のX線画像より各X線画像におけるカテーテル10の直線の部分を抽出し、抽出されたカテーテル10の直線の部分と、パラメータ情報取得部310より取得した各X線画像が撮像された際のパラメータ情報とを用いて、カテーテル10の3次元形状を復元し、復元されたカテーテル10の3次元形状よりカテーテル10の挿入方向に対する回転軸を決定する。また、カテーテル回転軸決定部312は、決定した回転軸及び第2画像取得部311より取得したX線画像を外部に出力する。
 カテーテル回転軸決定部312は、第2画像取得部311より取得した複数のX線画像から、各X線画像30におけるカテーテル10の直線の部分を抽出する。以下に、具体的な抽出方法について、第2画像取得部311より2つのX線画像30を取得した場合を例に説明する。
 カテーテル回転軸決定部312は、第2画像取得部311より角度差のある2つのX線画像30を取得する。カテーテル回転軸決定部312は、取得した2つのX線画像30のそれぞれについて、カテーテル10の直線の部分を抽出する。ここで、カテーテル10の直線の部分とは、図9に示すようなカテーテル10の湾曲部15を含まない部分、つまり直線部14である。カテーテル回転軸決定部312は、カテーテル10の直線度を判定し、所定の直線度以上の部分を直線部14として認識し、認識した直線部14を含む所定の領域をカテーテル10の直線の部分として抽出する。
 図21に、カテーテル10の回転軸を決定する概略図を示す。
 始めに、カテーテル回転軸決定部312は、2つのX線画像30より、カテーテル10の直線の部分を抽出した画像である直線部画像140及び直線部画像141を取得する。また、カテーテル回転軸決定部312は、直線部画像140及び直線部画像141を取得した際のそれぞれのX線発生部210及びX線発生部211とX線検出部220及びX線検出部221との撮像パラメータ情報を撮像パラメータ情報取得部310より取得する。
 次に、カテーテル回転軸決定部312は、2つの直線部画像140及び直線部画像141と、2組の撮像パラメータ情報とを用いて、カテーテル10の直線の部分の3次元形状復元を行なう。3次元形状復元の具体的な方法は、以下の通りである。
 まず、2枚の画像間で対応点を抽出するために、(9)式を用いて撮像パラメータ情報から2つの画像間の幾何学的な関係を示すFundamental行列を算出する。
Figure JPOXMLDOC01-appb-M000002
 次に、(10)式を用いて、X線発生器210及びX線検出器220によって得られた直線部画像140中の直線上の点に対するエピポーラ線を算出する。図26に、エピポーラ線の概略を示す。図26に示すとおり、図中の(a)が直線部画像140中の直線上の点であり、図中の(b)の点線が直線状の点に対するエピポーラ線である。
Figure JPOXMLDOC01-appb-M000003
 ここで、エピポーラ線は直線部画像140中における存在範囲を示す。従って、エピポーラ線と抽出した直線の交点を算出することで直線部画像140と直線部画像141との間の対応点を算出することができる。撮像パラメータ情報(各X線管の投影行列)と直線部画像140及び直線部画像141との間の対応点には(8)式のような投影式が成り立つため、各投影式を3次元点Xについてまとめると、(11)式のように変形することができる。
Figure JPOXMLDOC01-appb-M000004
 そのため、(11)式を変形した(12)式を用いて、撮像パラメータ情報(各投影行列)と直線部画像140及び直線部画像141との間の対応点から3次元点を復元することができる。ここで、Mは、(11)式の左辺の4×3行列の省略形であり、bは、(11)式の右辺の4×1行列の省略形である。
Figure JPOXMLDOC01-appb-M000005
 この処理を直線部画像140から抽出した直線上の複数の点に行うことにより、3次元直線を算出することができる。図22に、復元されたカテーテル10の3次元直線を示す。
 カテーテル回転軸決定部312は、上記の方法で復元されたカテーテル10の3次元形状より、カテーテル10の挿入方向に対する回転軸を決定する。
 <カテーテル形状記憶部313>
 カテーテル形状記憶部313は、複数種類のカテーテル10について、カテーテル10の先端の回転角毎にカテーテル10の先端の形状を記憶する。
 カテーテル10、又は、それに先行して挿入されるガイドワイヤには多くの種類があり、用途に応じて様々な形状がある。例えば、右冠動脈用又は左冠動脈用のジャドキンスと呼ばれるカテーテルがあり、それぞれの解剖学的位置に合わせて先端が曲げられて挿入しやすい形状になっている。
 カテーテル形状記憶部313は、それらに対応する様々な3次元形状のモデルを予め記憶している。そして、術者又は放射線機器を扱う放射線技師は、カテーテル10又はガイドワイヤを挿入する前に、そのタイプをカテーテル形状記憶部313に予め設定しておく。
 <カテーテル先端回転角算出部305>
 カテーテル先端回転角算出部305は、第1回転角測定部102より測定された第1回転角θと第2回転角測定部103より測定された第2回転角θとの角度差(θ-θ)、カテーテル挿入長測定部104より測定されたカテーテル挿入長L、及び第1の孔1010と第2の孔1011との距離L12を用いて、カテーテルの先端の回転角θを算出する。また、カテーテル先端回転角算出部305は、算出したカテーテルの先端の回転角θに加えて、第2回転角測定部103より測定された第2回転角θを外部へ出力する。
 <カテーテル先端回転角補正部314>
 図23に、カテーテル先端回転角補正部314の機能構成を示す。図23に示すように、カテーテル先端回転角補正部314は、取得部3141と、キャリブレーション用画像生成部3142と、比較部3143と、回転角度決定部3144とを有する。以下に、カテーテル先端回転角補正部314が有する各構成要素について説明する。
 <取得部3141>
 取得部3141は、第2画像取得部311より複数のX線画像の中から第1X線画像を、カテーテル先端回転角算出部305よりカテーテル10の先端の回転角θ及び第2回転角θを、カテーテル回転軸決定部312よりカテーテル回転軸を、撮像パラメータ情報取得部310より第1X線画像を撮影した際の撮像パラメータ情報をそれぞれ取得する。
 <キャリブレーション用画像生成部3142>
 キャリブレーション用画像生成部3142は、カテーテル形状記憶部313より、取得部3141より取得したカテーテル10の先端の回転角θから第2回転角θまでの回転角毎に、カテーテル10の先端の形状を複数取得し、取得したカテーテル10の先端の形状と回転角θaとが対応付けられたキャリブレーション用画像を複数生成する。
 図24に、キャリブレーション用画像の一例を示す。キャリブレーション用画像生成部3142がキャリブレーション用画像を生成する具体的な方法について説明する。
 始めに、キャリブレーション用画像生成部3142は、カテーテル形状記憶部313より取得した3次元形状モデルの回転軸を、カテーテル回転軸決定部312より取得した回転軸に合わせる。
 次に、キャリブレーション用画像生成部3142は、撮像パラメータ情報取得部310より取得した2方向の撮像パラメータ情報のどちらか一方を用いて、回転軸を合わせた3次元形状モデルを当該撮像パラメータ情報で撮像した場合に想定されるキャリブレーション用画像のうち、カテーテル10の先端の回転角θがカテーテル先端回転角算出部305で算出した回転角θから第2回転角θまでに相当する画像を一定の角度間隔(例えば1度毎)に生成する。ここで、生成範囲をθからθにしている理由は、捩れの一部あるいは全部が解消されて部分的に捩れが存在している場合、一様に捩れがあると想定して算出したカテーテル10の先端の回転角θから第2回転角θまでの間に実際の回転角度が存在するからである。
 <比較部3143>
 比較部3143は、取得部3141より取得した第1X線画像と、キャリブレーション用画像生成部3142より生成された複数のキャリブレーション用画像とを比較し、第1X線画像と最も類似度の大きいキャリブレーション用画像を決定する。
 図25に、X線画像とキャリブレーション用画像とを比較する概略図を示す。
 比較部3143が行なう比較の具体的な方法について説明する。比較部3143は、第1X線画像からエッジ又はコーナー等の特徴点を抽出し、その周囲の領域から特徴量を計算する。また、比較部3143は、同様に、複数のキャリブレーション用画像についても特徴量の計算を行なう。比較部3143は、第1X線画像の特徴量と、キャリブレーション用画像の各特長量とを比較し、最も特徴量の一致度が高いキャリブレーション用画像を、第1X線画像と最も類似度が大きいキャリブレーション用画像であると決定する。
 なお、比較部3143が行なう画像同士の類似度の計算方法は、モノクロ画像同士の類似度を計算する方法であれば良い。
 <回転角度決定部3144>
 回転角度決定部3144は、比較部3143より決定されたキャリブレーション用画像に対応する回転角θaを、カテーテル10の先端の回転角θと決定する。
 <提示部306>
 提示部306は、回転角度決定部3144より決定されたカテーテル10の先端の回転角θをユーザに提示する。このとき、図11Bに示すように、医師又は技師などのユーザにキャリブレーション(補正)を行っていることを示すため、「補正中」であることを表示してもよい。
 <動作>
 次に、以上のように構成されたカテーテル先端回転角測定装置300の各構成要素の動作に関して具体的に説明する。図27に、第3実施形態におけるカテーテル先端回転角測定装置300の処理動作のフローを示す。なお、図27において、図13と同様のステップについては同一の符号を付し、説明を適宜省略する。
 <ステップS490>
 撮像パラメータ情報取得部310は、撮像時の各X線撮像装置20の撮像条件及び姿勢に関する情報である撮像パラメータ情報を取得する。また、撮像パラメータ情報取得部310は、取得した撮像パラメータ情報をカテーテル先端回転角補正部314に出力する。
 <ステップS500>
 第2画像取得部311は、生体管内に挿入されたカテーテル10の先端部15を含む領域を撮像し、異なる位置に配される複数のX線撮像装置20により撮像されたX線画像30を複数取得する。また、第2画像取得部311は、取得した複数のX線画像30をカテーテル回転軸決定部312に出力する。
 <ステップS510>
 カテーテル回転軸決定部312は、第2画像取得部311より取得した複数のX線画像より各X線画像におけるカテーテル10の直線の部分を抽出し、抽出されたカテーテル10の直線の部分と、パラメータ情報取得部310より取得した各X線画像が撮像された際のパラメータ情報とを用いて、カテーテル10の3次元形状を復元し、復元されたカテーテル10の3次元形状よりカテーテル10の挿入方向に対する回転軸を決定する。また、カテーテル回転軸決定部312は、決定した回転軸及び第2画像取得部311より取得したX線画像を外部に出力する。また、カテーテル回転軸決定部312は、決定した回転軸及びX線画像をカテーテル先端回転角補正部314に出力する。
 <ステップS520>
 カテーテル先端回転角補正部314は、カテーテル10の先端の回転角を補正する。また、カテーテル先端回転角補正部314は、補正した結果を提示部306へ出力する。
 カテーテル先端回転角補正部314の動作に関する具体的な説明は後述する。
 <ステップS530>
 提示部306は、回転角度決定部3144より決定されたカテーテル10の先端の回転角θをユーザに提示する。
 以上の動作が、第3実施形態のおけるカテーテル先端回転角測定装置300の処理動作である。
 次に、カテーテル先端回転角測定装置300の特徴的な構成要素であるカテーテル先端回転角補正部314の動作に関して詳しく説明する。
 図28に、カテーテル先端回転角補正部314の処理動作のフローを示す。
 <ステップS5201>
 取得部3141は、第2画像取得部311より複数のX線画像の中から第1X線画像を、カテーテル先端回転角算出部305よりカテーテル10の先端の回転角θ及び第2回転角θを、カテーテル回転軸決定部312よりカテーテル回転軸を、撮像パラメータ情報取得部310より第1X線画像を撮影した際の撮像パラメータ情報をそれぞれ取得する。また、取得部3141は、取得したそれぞれの情報等をキャリブレーション用画像生成部3142及び比較部3143へ出力する。
 取得部3141が各情報等を取得するタイミングは、所定の時間間隔毎である。なお、取得部3141は、カテーテル挿入長測定部104よりカテーテル挿入長Lを取得し、カテーテル挿入長Lに変化があった場合に取得しても良い。
 <ステップS5202>
 キャリブレーション用画像生成部3142は、カテーテル形状記憶部313より、取得部3141より取得したカテーテル10の先端の回転角θから第2回転角θまでの回転角毎に、カテーテル10の先端の形状を複数取得し、取得したカテーテル10の先端の形状と回転角θaとが対応付けられたキャリブレーション用画像を複数生成する。また、キャリブレーション用画像生成部3142は、生成した複数のキャリブレーション用画像を比較部3143へ出力する。
 <ステップS5203>
 比較部3143は、取得部3141より取得した第1X線画像と、キャリブレーション用画像生成部3142より生成された複数のキャリブレーション用画像とを比較し、第1X線画像と最も類似度の大きいキャリブレーション用画像を決定する。また、比較部3143は、決定したキャリブレーション用画像を回転角度決定部3144へ出力する。
 <ステップS5204>
 回転角度決定部3144は、比較部3143より決定されたキャリブレーション用画像に対応する回転角θaを、カテーテル10の先端の回転角θと決定する。また、回転角度決定部3144は、決定したカテーテル10の先端の回転角θを提示部106へ出力する。
 <効果>
 以上のように、第3実施形態におけるカテーテル先端回転角測定装置300によれば、カテーテル10の一部に捩れが生じている場合においても、カテーテル10の先端の回転角を正確に測定することができる。したがって、術者は、カテーテル10の一部に捩れが生じている場合においても、カテーテル10の先端の回転角θを把握することができる。これにより、術者は、血管の分岐部において、カテーテルの先端を目的の血管等の分岐方向に合わせることが容易となり、スムーズなカテーテル操作を行なうことが出来る。また、スムーズなカテーテル操作により術時間が短縮されるため、患者に照射されるX線の照射時間が短縮される。そのため、患者へのX線照射による被爆量を低減できる。
 (第1変形例)
 例えば、カテーテル10の先端形状が鉤型の場合と直線の場合との少なくとも2種類の異なる形状のカテーテル10を使用する場合には、鉤型先端形状のカテーテル10の場合には前記した先端回転角の算出が必要であるが、直線の先端形状のカテーテル10の場合には前記した先端回転角の算出は不要となり、算出を停止させる場合がある。このように先端回転角の先端形状に応じて、算出の停止を自動的に判定する例について、第1変形例として、以下に説明する。なお、カテーテル10に代えて又はカテーテル10と共にガイドワイヤ16を使用する場合でも、先端形状が直線のガイドワイヤ16と、先端形状が鉤型のガイドワイヤ16との二種類を使用するときも、同様である。
 図30に示すように、カテーテル先端回転角測定装置400は、第1回転角測定部102と、第2回転角測定部103と、カテーテル挿入長測定部104と、カテーテル先端回転角算出部205と、提示部106と、第1画像取得部207と、カテーテル先端部抽出部208と、判定部230とを有する。また、図30には図示しないが、移動規制部101を有する。第1画像取得部207と、カテーテル先端部抽出部208と、判定部230と以外は、第1実施形態などと同様であるため、詳細な説明は省略し、異なる構成のみ以下に説明する。
 第1画像取得部207は、生体管内に挿入されたカテーテル10の先端部を含む領域を撮像するX線撮像装置20により撮像されたX線画像を取得する。
 カテーテル先端部抽出部208は、第1画像取得部207より取得されたX線画像からカテーテル10の先端部を抽出する。
 判定部230は、カテーテル先端部抽出部208より抽出されたカテーテル10の先端部がX線画像上において直線像であるか否かの直線像判定を行う。さらに、判定部230で、カテーテル10の先端部が直線像であると判定された場合は、更に先端部において所定閾値以上の濃度変化があるか否かの濃度変化判定を行う。判定部230による判定の結果、直線像でかつ先端部の濃度変化が無いと判定された場合は、判定部230でカテーテル先端回転角算出停止の判定を行う。
 図31A~図31Cは、それぞれ、カテーテル10(又はガイドワイヤ)の先端状態のX線撮像画像例を示している。図31Aは、直線像判定時に直線像ではないと判定部230で判定された場合(すなわち、例えば、先端形状が鉤型の場合)を示している。図31Bは、直線像判定時に直線像と判定部230で判定されたが、先端部濃度変化判定時に濃度変化があると判定部230で判定された場合(鉤型の先端部が重畳して濃く投影されている場合)を示している。図31Cは、直線像判定時に直線像と判定部230で判定され、更に先端部濃度変化判定時に濃度変化が無いと判定部230で判定された場合(先端形状が直線の場合)を示している。これらの画像については、以下のフローの説明中で詳しく説明する。
 このようなカテーテル先端回転角測定装置400の動作を図32に示す。なお、図32のステップS410、S420、S430、S450は、先の実施形態と同様であるため、説明を省略し、異なる部分のみ以下に説明する。
 第1画像取得部207では、第2実施形態と同様に、生体管内に挿入されたカテーテル10の先端部を含む領域を撮像するX線撮像装置20により撮像されたX線画像30を取得する(図32のステップS460)。
 次いで、カテーテル先端部抽出部208では、第1画像取得部207より取得されたX線画像30からカテーテルの先端部15を抽出する(図32のステップS470)。このとき、X線画像30において、まず、カテーテル10の先端部が直線像であるか否かを判定部230で判定する(図32のステップS471)。カテーテル10の先端部が直線像ではない(図31A参照)と判定部230で判定すれば(図32のステップS471でNOの場合)、カテーテル10の先端部形状は鉤型であるため、先端回転角の算出は必要と判定部230で判定し(又は、判定部230で何ら判定することなく)、カテーテル先端回転角算出部205でカテーテル先端回転角の算出を行う(図32のステップS440)。一方、X線画像30において、カテーテル10の先端部が直線像であるとき(図32のステップS471でYESの場合)は、先端部において所定閾値以上の濃度変化があるか否かの濃度変化の判定を判定部230で行う(図32のステップS472)。このような処理を行うのは、以下の理由による。X線画像30において、カテーテル10の先端部が直線像であるときは、カテーテル10の先端部形状は鉤型であるにもかかわらず、回転角度により重なって撮像された場合と、カテーテル10の先端部形状が直線の場合との2つの場合が考えられる。前者の場合には、直線像に見えていても、先端部がX線を受けたとき、先端部の一部では、X線が2回減衰するため、濃くなる部分が発生する(図31BのX1参照)。すなわち、一部に減衰率が大きい部分があり、閾値以上の濃淡画像の部分があれば、先端が鉤型であることになる。逆に、そのような濃淡部分が無ければ、先端は直線形状であることになる(図31C参照)。このため、カテーテル10の先端部が直線像であると判定部230で判定したときは、更に、先端部において所定閾値以上の濃度変化があるか否かの濃度変化の判定を判定部230で行う。
 判定部230による判定の結果、直線像でありかつ先端部の濃度変化があると判定部230で判定された場合(図32のステップS472でYESの場合)は、カテーテル10の先端部形状は鉤型である。このため、先端回転角の算出は必要と判定部230で判定し、カテーテル先端回転角算出部205でカテーテル先端回転角の算出を行う(図32のステップS440)。
 一方、判定部230による判定の結果、直線像であるが、先端部の濃度変化は無いと判定部230で判定された場合(図32のステップS472でNOの場合)は、カテーテル10の先端部形状は直線である。このため、先端回転角の算出は不要と判定部230で判定し、カテーテル先端回転角算出部205でカテーテル先端回転角の算出を停止する(図32において、ステップS440を行うことなく、フローを終了する)。
 このような変形例によれば、カテーテル10の先端部形状が直線であり、先端回転角の算出は不要であるにもかかわらず、無駄な算出処理を自動的に省略することができる。
 (第2変形例)
 ガイドワイヤ16とカテーテル10との両方の先端回転角を算出する例を、第2変形例として以下に説明する。カテーテル10の先端回転角を算出するためのカテーテル先端回転角測定装置100において、さらに、ガイドワイヤ16の先端回転角を算出するための装置を追加して備えて構成している。すなわち、図33に示すように、カテーテル先端回転角測定装置500は、第1回転角測定部102と、第2回転角測定部103と、第3回転角測定部112と、第4回転角測定部113と、ガイドワイヤ挿入長測定部114と、ガイドワイヤ先端回転角算出部115と、対象切替部116とを有する。また、ガイドワイヤ用の移動規制部111をさらに有する。図33において、10はカテーテル、110はガイドワイヤである。一般的に、カテーテル10よりもガイドワイヤ110の方が長い。
 図34にカテーテル先端回転角測定装置500のブロック図を示す。カテーテル先端回転角測定装置500では、カテーテル先端回転角とガイドワイヤ先端回転角をそれぞれ別々に算出する。
 第3回転角測定部112は、第1回転角測定部102と同様に、第3の孔1110におけるガイドワイヤ110の回転角であるガイドワイヤ110の第3回転角θを測定する。第3回転角測定部112は、第1回転角測定部102と同様に、撮影部1120(撮影部1020に対応)及び計算部1121(計算部1021に対応)より構成される。
 第4回転角測定部113は、第2回転角測定部103と同様に、第4の孔1111におけるガイドワイヤ110の回転角であるガイドワイヤ110の第4回転角θを測定する。第4回転角測定部113は、第2回転角測定部103と同様に、撮影部1130(撮影部1030に対応)及び計算部1131(計算部1031に対応)より構成される。
 ガイドワイヤ挿入長測定部114は、カテーテル挿入長測定部104と同様に、第3の孔1110又は第4の孔1111からガイドワイヤ110の先端部までの長さであるガイドワイヤ挿入長Lを測定する。
 ガイドワイヤ先端回転角算出部105は、カテーテル先端回転角算出部105と同様に、第3回転角測定部112より測定された第3回転角θと第4回転角測定部113より測定された第4回転角θとの角度差(θ-θ)、ガイドワイヤ挿入長測定部114より測定されたガイドワイヤ挿入長L、及び第3の孔1110と第4の孔1111との距離L34を用いて、ガイドワイヤの先端の回転角θを算出する。
 ガイドワイヤ用の移動規制部111は、移動規制部101と同様に、第1の孔1110(第1の孔1010に対応)から第2の孔1111(第2の孔1011に対応)まで貫通した筒状の形状を有し、第1の孔1110から第2の孔1111までガイドワイヤ110を通すことによりガイドワイヤ110の移動を挿入方向に規制するものである。
 対象切替部116は、カテーテル挿入長(L又はL)とガイドワイヤ挿入長(L)とのどちらかに変化があったタイミングで、対象切替処理を行う。具体的には、L>L-L13(ただし、L13はカテーテル10とガイドワイヤ110とのそれぞれの移動規制部(101、111)の間の距離(第1の孔1010と第3の孔1110との間の距離)である。この距離は、事前に設定される。)の場合にはカテーテル先端回転角を選択し、L≦L-L13の場合にはガイドワイヤ先端回転角を選択して、提示部106への送信情報(提示対象及びその先端回転角)を切り替える。
 図35に、カテーテル先端回転角測定装置500の動作のフローチャートを示す。
 まず、ステップS410で、第1回転角測定部102は、第1の孔1010におけるカテーテル10の回転角である第1回転角θを測定する。
 次いで、ステップS420で、第2回転角測定部103は、第2の孔1011におけるカテーテル10の回転角である第2回転角θを測定する。
 次いで、ステップS430で、カテーテル挿入長測定部104は、第1の孔1010からカテーテル10の先端部15までの長さであるカテーテル挿入長を測定する。
 次いで、ステップS440で、カテーテル先端回転角算出部105は、第1回転角測定部102より測定された第1回転角θと第2回転角測定部103より測定された第2回転角θとの角度差(θ-θ)、カテーテル挿入長測定部104より測定されたカテーテル挿入長L、及び第1の孔1010と第2の孔1011との距離L12を用いて、カテーテルの先端の回転角θを算出する。
 これらのステップS410~S440と同時的に又は別々に、ステップS411~S441を行う。
 すなわち、ステップS411で、第3回転角測定部112は、第3の孔1110におけるガイドワイヤ110の回転角である第3回転角θを測定する。
 次いで、ステップS421で、第4回転角測定部113は、第4の孔1111におけるガイドワイヤ110の回転角である第4回転角θを測定する。
 次いで、ステップS431で、ガイドワイヤ挿入長測定部114は、第3の孔1110からガイドワイヤ110の先端部までの長さであるガイドワイヤ挿入長Lを測定する。
 次いで、ステップS441で、ガイドワイヤ先端回転角算出部115は、第3回転角測定部112より測定された第3回転角θと第4回転角測定部113より測定された第4回転角θとの角度差(θ-θ)、ガイドワイヤ挿入長測定部114より測定されたガイドワイヤ挿入長L、及び第3の孔1110と第4の孔1111との距離L34を用いて、ガイドワイヤ110の先端の回転角θgを算出する。
 次いで、ステップS442で、対象切替部116は、ステップS430で測定されたカテーテル挿入長Lをカテーテル挿入長測定部104から取得するとともに、ステップS431で測定されたガイドワイヤ挿入長Lをガイドワイヤ挿入長測定部114から取得する。また、対象切替部116は、ステップS440で算出されたカテーテル10の先端の回転角θをカテーテル先端回転角算出部105から取得するとともに、ステップS441で算出されたガイドワイヤ110の先端の回転角θgをガイドワイヤ先端回転角算出部115から取得する。そして、対象切替部116は、L>L-L13が成立するか否かを判定する。
 L>L-L13が成立する場合(ステップS441でYESの場合)には、カテーテル先端回転角を選択して、ステップS450に進む。
 L>L-L13が成立しない場合(ステップS441でNOの場合)には、ガイドワイヤ先端回転角を選択して、ステップS450に進む。
 次いで、ステップS450で、提示部106は、ステップS443又はS444で選択された先端回転角を対象切替部116から取得し、取得した先端回転角をユーザに提示する。
 この第2変形例によれば、カテーテル10のみではなく、カテーテル10とガイドワイヤ110との両方を使用する場合にも、カテーテル10の先端回転角とガイドワイヤ110の先端回転角とを別々に精度良く測定することができる。よって、例えば、診断部位又は治療部位の近くまでは、ガイドワイヤ110のみを使用してするため、ガイドワイヤ110の先端回転角をカテーテル先端回転角測定装置により精度良く測定しておき、診断部位又は治療部位の近くになると、ガイドワイヤ110の挿入を停止して、代わりにガイドワイヤ110に沿わせてカテーテル10を挿入するため、カテーテル10の先端回転角をカテーテル先端回転角測定装置により精度良く測定することが可能となる。
 (第3変形例)
 第1実施形態にかかる先端回転角算出と、第3実施形態におけるキャリブレーションを考慮した先端回転角算出とを切り替えることができる例について、第3変形例として以下に説明する。
 この第3変形例にかかるカテーテル先端回転角測定装置600は、データベースとして補正切替情報記憶部320と、補正切替部321とを、例えば、第2変形例のカテーテル先端回転角測定装置500と図20の第3実施形態におけるカテーテル先端回転角測定装置300とを組み合わせた装置に備える構成としている
 補正切替情報記憶部320は、補正切替のために必要な情報(補正切替情報)を記憶している。例えば、補正切替情報記憶部320は、カテーテルか、ガイドワイヤかの種別情報、挿入長の差の判定用の所定閾値、病変箇所又は診断部位又は治療部位(例えば脳又は心臓)までの挿入長、などの補正切替のために必要な情報を記憶している。
 補正切替部321は、補正切替情報記憶部320に記憶された補正切替情報と、対象切替部116からの情報により、補正の有無を切り替えて、結果をカテーテル先端回転角補正部314と提示部306とに出力する。よって、補正切替部321は、例えば、対象切替部116から先端回転角の値とその種別(カテーテル又はガイドワイヤ)及び対応する挿入長を取得し、補正切替情報記憶部320から補正切替情報を取得し、それらの情報に基づき補正切替を行う。例えば、補正切替部321は、先端回転角の種別がガイドワイヤであり、且つ各挿入長の差が所定閾値より小さい場合、先端回転角補正は行わずに提示部306に先端回転角の値とその種別を出力する。それ以外の場合は、カテーテル先端回転角補正処理をカテーテル先端回転角補正部314で実施する。
 このような構成によれば、一例として、カテーテル10を使用中でありかつカテーテル10の先端回転角を算出すると補正切替部321により判定する場合は、補正を行うとして、キャリブレーションを考慮した先端回転角の算出を行うように、カテーテル先端回転角補正部314と提示部306とに出力する。
 別の例として、ガイドワイヤ110の先端回転角を算出する場合でかつ、挿入長の差が所定閾値以上となったと補正切替部321により判定する場合には、キャリブレーション(補正)を考慮した先端回転角の算出を行い、所定閾値より小さいと補正切替部321により判定する場合には、キャリブレーション(補正)無しの先端回転角の算出を行うように、カテーテル先端回転角補正部314と提示部306とに出力する。
 又は、別の例として、脳又は心臓までの挿入長を補正切替情報記憶部320に予め記憶しておき、その挿入長に到達したと補正切替部321により判定することにより脳又は心臓などの血管構造が複雑な場所に到達したと補正切替部321で推定したときには、キャリブレーション(補正)を考慮した先端回転角の算出を行うように、カテーテル先端回転角補正部314と提示部306とに出力する。
 病変箇所までの挿入長を補正切替情報記憶部320に予め記憶しておき、その挿入長に到達したと補正切替部321により判定することにより病変箇所に到達したと補正切替部321で推定したときには、キャリブレーション(補正)を考慮した先端回転角の算出を行うように、カテーテル先端回転角補正部314と提示部306とに出力する。
 この第3変形例によれば、補正を常時行うのではなく、必要なときのみ補正処理を行うことにより、効率良くかつ迅速に精度良く、先端回転角測定動作を行うことができる。
 なお、カテーテル先端回転角測定装置を使用する技師又は医者などのユーザからの入力により、補正を切替える場合には、ユーザ入力部322から補正切替部321に補正切替情報を入力することにより、補正動作を任意のタイミングで切替えることが可能となる。
 なお、本発明を第1~第3実施形態及び変形例に基づいて説明してきたが、本発明は、前記の第1~第3実施形態及び変形例に限定されないのはもちろんである。
 なお、上記様々な実施形態又は変形例のうちの任意の実施形態又は変形例を適宜組み合わせることにより、それぞれの有する効果を奏するようにすることができる。
 本発明にかかるカテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラムは、カテーテルの先端の回転角を測定する機能を有し、カテーテル造影検査及び治療時に有用である。また、研修医のためのカテーテル操作教育等の用途にも応用できる。
 本発明は、添付図面を参照しながら好ましい実施形態に関連して充分に記載されているが、この技術の熟練した人々にとっては種々の変形又は修正は明白である。そのような変形又は修正は、添付した請求の範囲による本発明の範囲から外れない限りにおいて、その中に含まれると理解されるべきである。

Claims (16)

  1.  生体管内に挿入されたカテーテルの先端の回転角を測定する装置であって、
     第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部と、
     前記第1の孔における前記カテーテルの回転角である第1回転角を測定する第1回転角測定部と、
     前記第2の孔における前記カテーテルの回転角である第2回転角を測定する第2回転角測定部と、
     前記第1の孔又は前記第2の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長を測定するカテーテル挿入長測定部と、
     前記第1回転角測定部より測定された第1回転角と前記第2回転角測定部より測定された第2回転角との角度差、前記カテーテル挿入長測定部より測定されたカテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出部と、
     を備えるカテーテル先端回転角測定装置。
  2.  前記カテーテルと、前記カテーテルに先行して挿入されるガイドワイヤとの両方を使用する場合には、前記生体管内において、前記カテーテルの先端の前記回転角を測定することに加えて、前記カテーテルに代えて前記ガイドワイヤの先端の回転角を算出し、前記カテーテルとして親カテーテルと、親カテーテルの先へ挿入される子カテーテルとを使用する場合には、前記子カテーテルの先端の回転角を測定することに加えて、前記親カテーテルの先端の回転角を算出する、請求項1に記載のカテーテル先端回転角測定装置。
  3.  前記生体管内に挿入された前記カテーテルの前記先端部を含む領域を撮像するX線撮像装置と、
     前記X線撮像装置で撮像された画像から前記カテーテルの前記先端部を含む領域を抽出するカテーテル先端部抽出部と、
     前記カテーテル先端部抽出部で抽出された情報を基に、前記カテーテル先端回転角算出部により前記カテーテルの先端の回転角を算出するか否かを判定する判定部とを備える、
    請求項1又は2に記載のカテーテル先端回転角測定装置。
  4.  前記判定部は、前記X線撮像装置で撮像された画像から抽出した前記カテーテルの前記先端部の画像を基に、前記カテーテルの前記先端部が直線像であるか否かの直線像判定を行うとともに、前記判定部で、前記カテーテルの前記先端部が直線像であると判定された場合は、更に前記先端部において所定閾値以上の濃度変化があるか否かの濃度変化の判定を行い、前記所定閾値以上の濃度変化があれば、前記カテーテルの前記先端部が直線ではなく、先端回転角測定が必要であると判定する一方、前記所定閾値以上の濃度変化がなければ、前記カテーテルの前記先端部が直線であり、先端回転角測定が不要であると判定する、
    請求項3に記載のカテーテル先端回転角測定装置。
  5.  さらに、前記カテーテル先端回転角測定装置は、
     前記カテーテル先端回転角算出部より算出された前記カテーテルの先端の回転角をユーザに提示する提示部と、
     を備える請求項1又は2に記載のカテーテル先端回転角測定装置。
  6.  前記カテーテル先端回転角算出部は、前記第1回転角測定部より測定される第1回転角をθとし、前記第2回転角測定部より測定される第2回転角をθとし、前記第1の孔から測定されるカテーテル挿入長をLとし、前記第1の孔と前記第2の孔との距離をL12とした場合、カテーテルの先端の回転角であるθをθ=θ-(θ-θ)・L/L12より算出する、
     請求項1又は2に記載のカテーテル先端回転角測定装置。
  7.  前記カテーテル先端回転角測定装置は、さらに、
     生体管内に挿入された前記カテーテルの先端部を含む領域を撮像するX線撮像装置により撮像されたX線画像を取得する第1画像取得部と、
     前記第1画像取得部より取得されたX線画像から前記カテーテルの先端部を抽出するカテーテル先端部抽出部と、
     前記カテーテル先端部抽出部より抽出された前記カテーテルの先端部の動作を検出する動作検出部と、
     を備え、
     前記カテーテル先端回転角算出部は、前記動作検出部が動作を検出した場合、前記カテーテルの先端の回転角を算出し、前記動作検出部が動作を検出しなかった場合、前記カテーテルの先端の回転角を算出しない、
     請求項1又は2に記載のカテーテル先端回転角測定装置。
  8.  前記カテーテル先端回転角測定装置は、さらに、
     撮像時の各X線撮像装置の撮像条件及び姿勢に関する情報である撮像パラメータ情報を取得する撮像パラメータ情報取得部と、
     異なる位置に配される複数のX線撮像装置により撮像された前記カテーテルの先端部を含むX線画像を複数取得する第2画像取得部と、
     前記第2画像取得部より取得した複数のX線画像より各X線画像における前記カテーテルの直線の部分を抽出し、抽出された前記カテーテルの直線の部分と、前記撮像パラメータ情報取得部より取得した前記各X線画像が撮像された際の撮像パラメータ情報とを用いて、前記カテーテルの3次元形状を復元し、復元された前記カテーテルの3次元形状より前記カテーテルの挿入方向に対する回転軸であるカテーテル回転軸を決定するカテーテル回転軸決定部と、
     複数種類のカテーテルについて、前記カテーテルの先端の回転角毎に前記カテーテルの先端の形状を記憶するカテーテル形状記憶部と、
     カテーテル先端回転角補正部と、
     を備え、
     前記カテーテル先端回転角補正部は、
     前記第2画像取得部より前記複数のX線画像の中から第1X線画像を、前記カテーテル先端回転角算出部より前記カテーテルの先端の回転角及び前記第2回転角を、前記カテーテル回転軸決定部より前記カテーテル回転軸を、撮像パラメータ情報取得部より前記第1X線画像を撮影した際の撮像パラメータ情報をそれぞれ取得する取得部と、
     前記カテーテル形状記憶部より、前記取得部より取得した前記カテーテルの先端の回転角から前記第2回転角までの回転角毎に、前記カテーテルの先端の形状を複数取得し、取得した前記カテーテルの先端の形状と前記回転角とが対応付けられたキャリブレーション用画像を複数生成するキャリブレーション用画像生成部と、
     前記取得部より取得した前記第1X線画像と、前記キャリブレーション用画像生成部より生成された複数のキャリブレーション用画像とを比較し、前記第1X線画像と最も類似度の大きいキャリブレーション用画像を決定する比較部と、
     前記比較部より決定されたキャリブレーション用画像に対応する前記回転角を、前記カテーテルの先端の回転角と決定する回転角度決定部とを有し、
     前記提示部は、前記回転角度決定部より決定された前記カテーテルの先端の回転角をユーザに提示する、
     請求項5に記載のカテーテル先端回転角測定装置。
  9.  前記生体管内に前記カテーテルの先端部が到達すべき治療部位に応じて、前記カテーテル先端回転角補正部で補正を行うか否かを切替える補正切替部をさらに備える、
     請求項8に記載のカテーテル先端回転角測定装置。
  10.  前記生体管内に挿入された前記カテーテルの挿入長さに応じて、前記カテーテル先端回転角補正部で補正を行うか否かを切替える補正切替部をさらに備える、
     請求項8に記載のカテーテル先端回転角測定装置。
  11.  前記生体管内に挿入された前記カテーテルの先端部の画像に応じて、前記カテーテル先端回転角補正部で補正を行うか否かを切替える補正切替部をさらに備える、
     請求項8に記載のカテーテル先端回転角測定装置。
  12.  第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部を用いて、生体管内に挿入されたカテーテルの先端の回転角を測定する装置であって、
     第1回転角測定部より測定された前記第1の孔における前記カテーテルの回転角である第1回転角、第2回転角測定部より測定された前記第2の孔における前記カテーテルの回転角である第2回転角、及びカテーテル挿入長測定部より測定された前記第1の孔又は前記第2の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長をそれぞれ取得し、前記第1回転角と前記第2回転角との角度差、前記カテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出部を備えるカテーテル先端回転角測定装置。
  13.  第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部を用いて、生体管内に挿入された前記カテーテルの先端の回転角を測定する方法であって、
     前記第1の孔における前記カテーテルの回転角である第1回転角を測定する第1回転角測定ステップと、
     前記第2の孔における前記カテーテルの回転角である第2回転角を測定する第2回転角測定ステップと、
     前記第1の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長を測定するカテーテル挿入長測定ステップと、
     前記第1回転角測定ステップより測定された第1回転角と前記第2回転角測定ステップより測定された第2回転角との角度差、前記カテーテル挿入長測定ステップより測定されたカテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出ステップと、
     を含むカテーテル先端回転角測定方法。
  14.  第1の孔から第2の孔まで貫通した筒状の形状を有し、前記第1の孔から前記第2の孔まで前記カテーテルを通すことにより前記カテーテルの移動を挿入方向に規制する移動規制部を用いて、生体管内に挿入されたカテーテルの先端の回転角を測定する方法であって、
     第1回転角測定部より測定された前記第1の孔における前記カテーテルの回転角である第1回転角、第2回転角測定部より測定された前記第2の孔における前記カテーテルの回転角である第2回転角、及びカテーテル挿入長測定部より測定された前記第1の孔又は前記第2の孔から前記カテーテルの先端部までの長さであるカテーテル挿入長をそれぞれ取得し、前記第1回転角と前記第2回転角との角度差、前記カテーテル挿入長、及び前記第1の孔と前記第2の孔との距離を用いて、前記カテーテルの先端の回転角を算出するカテーテル先端回転角算出ステップを含むカテーテル先端回転角測定方法。
  15.  請求項13に記載のカテーテル先端回転角測定方法をコンピュータに実行させるためのカテーテル先端回転角測定プログラム。
  16.  請求項14に記載のカテーテル先端回転角測定方法をコンピュータに実行させるためのカテーテル先端回転角測定プログラム。
PCT/JP2013/003485 2012-06-07 2013-06-03 カテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラム WO2013183276A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2014506652A JP5526297B1 (ja) 2012-06-07 2013-06-03 カテーテル先端回転角測定装置、カテーテル先端回転角測定装置の駆動方法及びカテーテル先端回転角測定プログラム
US14/197,802 US9918682B2 (en) 2012-06-07 2014-03-05 Catheter tip-end rotation angle measurement apparatus, catheter tip-end rotation angle measurement method, and catheter tip-end rotation angle measurement program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012129552 2012-06-07
JP2012-129552 2012-06-07

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/197,802 Continuation US9918682B2 (en) 2012-06-07 2014-03-05 Catheter tip-end rotation angle measurement apparatus, catheter tip-end rotation angle measurement method, and catheter tip-end rotation angle measurement program

Publications (1)

Publication Number Publication Date
WO2013183276A1 true WO2013183276A1 (ja) 2013-12-12

Family

ID=49711682

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/003485 WO2013183276A1 (ja) 2012-06-07 2013-06-03 カテーテル先端回転角測定装置、カテーテル先端回転角測定方法及びカテーテル先端回転角測定プログラム

Country Status (3)

Country Link
US (1) US9918682B2 (ja)
JP (1) JP5526297B1 (ja)
WO (1) WO2013183276A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180092388A (ko) * 2017-02-09 2018-08-20 (주)엘메카 카테터 이동의 정밀 제어 가능한 의료용 석션기

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10004467B2 (en) * 2014-04-25 2018-06-26 Medtronic, Inc. Guidance system for localization and cannulation of the coronary sinus
FR3024651A1 (fr) * 2014-08-11 2016-02-12 Eurl Cornier Dispositif de visualisation de l'uretre feminin en imagerie
WO2016040614A1 (en) * 2014-09-10 2016-03-17 The University Of North Carolina At Chapel Hill Radiation-free simulator system and method for simulating medical procedures
JP6595232B2 (ja) 2015-07-02 2019-10-23 ソニー・オリンパスメディカルソリューションズ株式会社 内視鏡用撮像装置、内視鏡装置、及び内視鏡用ケーブル
JP6981793B2 (ja) * 2017-06-30 2021-12-17 キヤノンメディカルシステムズ株式会社 医用画像処理装置、医用画像診断装置及びプログラム
KR102125228B1 (ko) * 2018-03-28 2020-06-22 강원대학교산학협력단 내시경적 유두괄약근 절개도의 카테터 말단부 굴곡의 동작 특성을 측정하기 위한 시험장치
CN108630049A (zh) * 2018-04-28 2018-10-09 苏州敏行医学信息技术有限公司 置管检测模组及具有该模组的置管术虚拟仿真设备
CN109064861B (zh) * 2018-09-20 2024-02-20 武汉湾流科技股份有限公司 一种用于腹腔镜模拟手术的腹腔镜器械钳子
CN115721408B (zh) * 2022-11-21 2024-05-14 绍兴梅奥心磁医疗科技有限公司 一种可视化消融装置
CN116570818B (zh) * 2023-05-22 2024-04-16 极限人工智能有限公司 导管控制方向与影像动作方向的一致性校准方法及***

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09225023A (ja) * 1996-02-21 1997-09-02 Buaayu:Kk バルーンカテーテル
JP2007306993A (ja) * 2006-05-16 2007-11-29 Olympus Medical Systems Corp 回転自走式内視鏡システム
JP2009507543A (ja) * 2005-09-08 2009-02-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ イメージングシステム用磁気トラッキングシステム
JP2012035080A (ja) * 2010-08-05 2012-02-23 Biosense Webster (Israel) Ltd カテーテルのもつれの指示

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07246183A (ja) 1994-03-11 1995-09-26 Olympus Optical Co Ltd 電子内視鏡システム
JP3872210B2 (ja) 1998-06-17 2007-01-24 株式会社東海理化電機製作所 カテーテル操作シミュレータ及びそれを用いたシミュレーション方法
US8945095B2 (en) 2005-03-30 2015-02-03 Intuitive Surgical Operations, Inc. Force and torque sensing for surgical instruments
US20080146942A1 (en) * 2006-12-13 2008-06-19 Ep Medsystems, Inc. Catheter Position Tracking Methods Using Fluoroscopy and Rotational Sensors
JP2009162920A (ja) 2007-12-28 2009-07-23 Macroware Co Ltd カテーテル動作検出装置およびコンピュータプログラム
US9549707B2 (en) * 2011-11-01 2017-01-24 Siemens Medical Solutions Usa, Inc. System for orientation and display of ultrasound imaging data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09225023A (ja) * 1996-02-21 1997-09-02 Buaayu:Kk バルーンカテーテル
JP2009507543A (ja) * 2005-09-08 2009-02-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ イメージングシステム用磁気トラッキングシステム
JP2007306993A (ja) * 2006-05-16 2007-11-29 Olympus Medical Systems Corp 回転自走式内視鏡システム
JP2012035080A (ja) * 2010-08-05 2012-02-23 Biosense Webster (Israel) Ltd カテーテルのもつれの指示

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180092388A (ko) * 2017-02-09 2018-08-20 (주)엘메카 카테터 이동의 정밀 제어 가능한 의료용 석션기
KR101958286B1 (ko) 2017-02-09 2019-03-15 (주)엘메카 카테터 이동의 정밀 제어 가능한 의료용 석션기

Also Published As

Publication number Publication date
US20140187921A1 (en) 2014-07-03
US9918682B2 (en) 2018-03-20
JPWO2013183276A1 (ja) 2016-01-28
JP5526297B1 (ja) 2014-06-18

Similar Documents

Publication Publication Date Title
JP5526297B1 (ja) カテーテル先端回転角測定装置、カテーテル先端回転角測定装置の駆動方法及びカテーテル先端回転角測定プログラム
US6370417B1 (en) Method for positioning a catheter in a vessel, and device for implementing the method
US8718338B2 (en) System and method to compensate for respiratory motion in acquired radiography images
CA2559308C (en) Method and system for determining a three dimensional representation of an organ
US20090310847A1 (en) Medical image processing apparatus and x-ray diagnosis apparatus
US20050027193A1 (en) Method for automatically merging a 2D fluoroscopic C-arm image with a preoperative 3D image with one-time use of navigation markers
US20080199059A1 (en) Information Enhanced Image Guided Interventions
WO2014024422A1 (ja) カテーテル先端回転角度検出装置、方法、及び、プログラム
JP2012519528A (ja) 関心領域を医療画像上に表示するための医療画像観察システム
US11127153B2 (en) Radiation imaging device, image processing method, and image processing program
JP2007282974A (ja) 循環器用画像診断装置
US11488313B2 (en) Generating a motion-compensated image or video
US20140093030A1 (en) X-ray ct system, image display device, and image display method
US8363901B2 (en) Radiographic tomography apparatus
JP6662428B2 (ja) 動態解析システム
US10765392B2 (en) X-ray imaging apparatus and x-ray image display method
US20100189217A1 (en) X-ray diagnosis apparatus
US11481957B2 (en) Medical image processing apparatus and storage medium
US8348507B2 (en) Object localization in X-ray images
US8175684B2 (en) Method for processing images and associated medical imaging system
KR20140120145A (ko) 방사선 불투과성 반구형 입체 표지를 기반으로 하는 ct 혈관 조영 영상과 x-선 혈관 조영 영상 간의 정합 방법
US8457720B2 (en) Apparatus for assisting with the positioning of an implant
US8467850B2 (en) System and method to determine the position of a medical instrument
JP2020036819A (ja) X線撮影装置およびx線画像処理方法
JP7434008B2 (ja) 医用画像処理装置およびプログラム

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2014506652

Country of ref document: JP

Kind code of ref document: A

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

Ref document number: 13800635

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

Country of ref document: EP

Kind code of ref document: A1