WO2004107982A1 - 超音波内視鏡装置 - Google Patents

超音波内視鏡装置 Download PDF

Info

Publication number
WO2004107982A1
WO2004107982A1 PCT/JP2004/008151 JP2004008151W WO2004107982A1 WO 2004107982 A1 WO2004107982 A1 WO 2004107982A1 JP 2004008151 W JP2004008151 W JP 2004008151W WO 2004107982 A1 WO2004107982 A1 WO 2004107982A1
Authority
WO
WIPO (PCT)
Prior art keywords
ultrasonic
image
surface shape
circuit
optical image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2004/008151
Other languages
English (en)
French (fr)
Inventor
Tomonao Kawashima
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Olympus Corp
Original Assignee
Olympus Corp
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 Olympus Corp filed Critical Olympus Corp
Priority to EP04736120A priority Critical patent/EP1632184B1/en
Publication of WO2004107982A1 publication Critical patent/WO2004107982A1/ja
Priority to US11/294,985 priority patent/US20060183992A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/06Devices, other than using radiation, for detecting or locating foreign bodies ; Determining position of diagnostic devices within or on the body of the patient
    • A61B5/061Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
    • A61B5/062Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body using magnetic field
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/12Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4245Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • A61B8/4254Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient using sensors mounted on the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4461Features of the scanning mechanism, e.g. for moving the transducer within the housing of the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/445Details of catheter construction

Definitions

  • the present invention relates to an ultrasonic endoscope apparatus for acquiring an optical image and an ultrasonic image of a subject.
  • an ultrasonic transducer that scans with ultrasound and receives the echo is usually observed, and an optical image of the luminal surface of the esophagus, stomach, duodenum, colon, etc. is observed. There is an optical observation window for this.
  • the ultrasound scan plane is not visible in the optical image.
  • the ultrasound scan plane does not always fall within the field of view of the optical image.
  • the lumen is filled with degassed water, which is an ultrasonic medium, or a balloon is used to make the ultrasonic waves reach the affected area. Obstruction of the field of view was disturbed, and after all, observation using ultrasonic images was different from observation using optical images with a certain time interval.
  • diagnosis of the spread of a lesion in the horizontal direction and diagnosis of the depth of the lesion in the vertical direction are performed. Can be performed at the same time, and the diagnostic performance of ultrasonic endoscopy can be improved.
  • the first method is a three-dimensional measurement method based on triangulation, which is based on triangulation and is based on the light-section 'method' using slit light projection. This is the following method.
  • the second method is to use a mechanism that corrects the distortion caused by the wide-angle lens attached to the endoscope.
  • a three-dimensional measurement method based on a light sectioning method using a slit light projection involves guiding a laser beam. It is necessary to provide a special optical system for the endoscope. In the method of correcting distortion, it is necessary to attach a wide-angle lens to the endoscope. A special endoscope is required for the endoscope measurement system, and it must be separated from the endoscope for the three-dimensional ultrasonic image reconstruction system.
  • the three-dimensional measurement method is based on the premise that the object surface is a plain surface, and when projecting onto a surface with uneven gradation such as a lumen surface, the three-dimensional position of the object surface can be accurately grasped. It is hard to say that it can be done.
  • the present invention has been made in view of the above circumstances, and it is possible to display an optical image and an ultrasonic image in an accurate correspondence with each other's position and direction without having to replace an endoscope from a patient. It is an object of the present invention to provide an ultrasonic endoscope apparatus capable of performing the above. Disclosure of the invention
  • An ultrasonic endoscope apparatus is an optical endoscope that acquires an optical image of a subject.
  • a matching circuit that matches the optical image obtained by the optical image obtaining device with the ultrasonic image obtained by the ultrasonic image obtaining device based on the position information obtained by the position information obtaining device. Characterized by having
  • FIG. 1 is a block diagram showing the overall configuration of the ultrasonic endoscope apparatus according to the first embodiment of the present invention.
  • FIG. 2 is an enlarged cross-sectional view of the insertion-side tip of the insertion portion of the endoscope according to the first embodiment of the present invention.
  • FIG. 3 is a block diagram showing the image processing device according to the first embodiment of the present invention.
  • FIG. 4 is a block diagram showing a shape matching circuit according to the first embodiment of the present invention.
  • FIG. 5 is an explanatory diagram showing the operation of the manual scanning according to the first embodiment of the present invention.
  • FIG. 6 is a conceptual diagram showing the three-dimensional image data according to the first embodiment of the present invention.
  • FIG. 7 is an explanatory diagram showing reshuffling of three-dimensional image data according to the first embodiment of the present invention.
  • FIG. 8 is an explanatory diagram showing surface shape data according to the first embodiment of the present invention.
  • FIG. 9 is an explanatory diagram showing imaging of a region of interest by the endoscope according to the first embodiment of the present invention.
  • FIG. 10 is an explanatory diagram showing the synthesis of the optical image and the ultrasonic image according to the first embodiment of the present invention.
  • FIG. 11 is an explanatory diagram showing an image displayed on the monitor according to the first embodiment of the present invention.
  • FIG. 12 is a diagram for explaining an electronic radial scanning ultrasonic endoscope applicable to the first embodiment of the present invention and scanning thereof.
  • FIG. 13 is a diagram for explaining a COMPEX-scanning ultrasonic endoscope applicable to the first embodiment of the present invention and scanning thereof.
  • FIG. 14 is a diagram for explaining a two-dimensional array type ultrasonic endoscope applicable to the first embodiment of the present invention and its scanning.
  • FIG. 15 is a block diagram showing a shape matching circuit according to the second embodiment of the present invention.
  • FIG. 16 shows a shape matching circuit according to the third embodiment of the present invention. It is a block diagram showing a road.
  • FIG. 17 is a block diagram illustrating an image processing device according to the fourth embodiment of the present invention.
  • FIG. 18 is a block diagram illustrating an image processing device according to the fifth embodiment of the present invention.
  • FIG. 19 is an explanatory diagram showing an image displayed on the monitor according to the fifth embodiment of the present invention.
  • FIGS. 1 to 11 relate to a first embodiment of the present invention
  • FIG. 1 is a block diagram showing an overall configuration of an ultrasonic endoscope apparatus.
  • FIG. 2 is an enlarged cross-sectional view of the insertion-side tip of the insertion section of the endoscope.
  • FIG. 3 is a block diagram showing the image processing apparatus.
  • FIG. 4 is a block diagram showing a shape matching circuit.
  • FIG. 5 is an explanatory diagram showing the operation of manual scanning.
  • FIG. 6 is a conceptual diagram showing three-dimensional image data.
  • FIG. 7 is an explanatory diagram showing reshuffling of a three-dimensional image.
  • FIG. 8 is an explanatory view showing the surface shape.
  • FIG. 9 is an explanatory diagram showing imaging of a region of interest by an endoscope.
  • FIG. 10 is an explanatory diagram showing the synthesis of an optical image and an ultrasonic image.
  • FIG. 11 is an explanatory view showing an image displayed on the monitor. You.
  • an ultrasonic endoscope apparatus 1 includes an ultrasonic endoscope 2, an ultrasonic observation apparatus 3, an optical observation apparatus 4, a position detection apparatus 5, an image processing apparatus, It has a device 6, a monitor 7 as display means or a display device, a keyboard 8, and a mouse 9.
  • the ultrasonic endoscope 2 has an insertion section 21 and an operation section 25 connected to each other.
  • the insertion portion 21 is made of a flexible material, and can be inserted into the body cavity of the subject.
  • the operation unit 25 includes a motor 23 that drives the ultrasonic vibrator 22 at the tip of the input unit 21.
  • a transmission coil 24 for exciting a magnetic field in space is provided at the tip of the insertion section 21.
  • the position detection device 5 as a position information acquisition unit or a position information acquisition device includes a coil drive circuit 11, a plurality of reception coils 12, and a position calculation circuit 13.
  • the coil drive circuit 11 outputs a coil excitation signal to the transmission coil 24.
  • the plurality of receiving coils 12 are fixed to specific positions by a predetermined arrangement method, sequentially detect a magnetic field generated by the transmitting coil 24, and output an electrical reception signal.
  • the position calculation circuit 13 is based on the reception signal output from the reception coil 12 Calculate the data indicating the position and direction of the tip of the input section 21 (hereinafter referred to as position and direction).
  • the plurality of receiving coils 12 are integrally fixed to a rectangular parallelepiped casing.
  • this housing and the receiving coil 12 are collectively referred to as a receiving coil 10.
  • the receiving coil 12 is shown fixedly arranged in a straight line in the receiving coil unit 10 for the sake of space, but actually, it is placed on a two-dimensional plane or three-dimensional space. It shall be fixed side by side.
  • the insertion-side tip of the insertion section 21 will be described in detail with reference to FIG.
  • an acoustically translucent tip cap 27 made of a material such as polymethylpentene is provided at the tip 26 of the insertion portion 21.
  • An ultrasonic transducer 22 as an ultrasonic image acquisition means or an ultrasonic image acquisition device is provided inside the tip cap 27, and an ultrasonic transmission medium 28 is filled in the tip cap 27.
  • the ultrasonic vibrator 22 is connected to a flexible shaft 29 made of a flexible material.
  • the flexible shaft 29 is connected to the rotating shaft of the motor 23 in the operation section 25 of the endoscope 2 shown in FIG. 1, and is provided so as to rotate in the direction of the arrow in FIG. .
  • the ultrasonic transducer 22 is connected to the ultrasonic observation device 3 via the endoscope operation unit 25 via a signal line (not shown) in the flexible shaft 29. Outputs an echo signal to
  • the two transmission coils 24 are connected to a coil drive circuit 11 in the position detection device 5 via a signal line 30.
  • One of the transmission coils 24 has an axis with the direction marked "1 o'clock direction” in Fig. 2 and the other has a coiled conductor with the direction marked "normal direction” as the axis.
  • the “normal direction” that is wound is the direction of the insertion axis of the insertion portion 21, and the “12 o'clock direction” is the direction orthogonal to this.
  • the “normal direction” corresponds to the normal direction of the ultrasonic image obtained by the ultrasonic transducer 22 performing radial scanning.
  • the transmitting coil 24 wound in the 12 o'clock direction is provided so that the direction in which the conductor is wound coincides with the 12 o'clock direction of the ultrasonic image among the directions orthogonal to the normal direction. It is assumed that The operation when the ultrasonic transducer 22 performs radial scanning will be described later.
  • the tip portion 26 of the insertion portion 21 has a charge-coupled device solid-state imaging device camera (hereinafter referred to as a CCD camera) 31 for capturing an optical image in color, and a CCD camera 3.
  • a CCD camera charge-coupled device solid-state imaging device camera
  • An imaging light irradiation window 32 for irradiating the body cavity with light necessary for imaging in 1 is provided.
  • Reference numeral 31 is connected to the optical observation device 4, and outputs an imaging signal to the optical observation device 4 via the endoscope operation unit 25 via a signal line (not shown) in the input unit 21.
  • the optical observation device 4 creates an optical image in the body cavity based on the image signal.
  • the imaging light irradiation window 32 is configured so that imaging light from a light source device (not shown) reaches through a light guide path (not shown) such as an optical fiber provided in the inlet 21. Thus, the imaging light required for imaging by the CCD camera 31 can be emitted into the body cavity.
  • a rigid frame 33 for integrally holding the above-described portions of the distal end of the insertion portion 21 is provided as shown in FIG.
  • the image processing apparatus 6 includes an ultrasonic image memory 41, a three-dimensional data construction circuit 42, and a large-capacity first three-dimensional image memory (hereinafter simply referred to as a three-dimensional image memory 4). 3), cross-section extraction circuit 44, cross-section image memory 45, surface extraction circuit 46 as surface shape calculation means or surface shape calculation circuit, optical image memory 47, surface shape calculation means or surface A surface shape estimation circuit 48 as a shape calculation circuit, a large-capacity second 3D image memory (hereinafter simply referred to as a 3D image memory 49), and a shape matching circuit as a matching means. Circuit 50, a coordinate conversion circuit 51, a surface image memory 52, a synthesizing circuit 53, a display circuit 54, a switch 55, and a controller 56 for controlling these. are doing.
  • the switch 55 switches the output destination of the data from the position detection device 5 to one of the ultrasonic image memory 41 and the optical image memory 47.
  • the controller 56 controls these units and circuits according to the input from the keyboard 8 and the mouse 9. '-Next, the configuration of the shape matching circuit 50 shown in FIG. 3 will be described in detail with reference to FIG.
  • the shape matching circuit 50 includes a first surface shape memory 57, a second surface shape memory 58, and a cross-correlation circuit 59.
  • the position detecting device 5 of the present embodiment is configured to acquire information on the position and direction of the ultrasonic image acquiring means with respect to the subject at the time of acquiring the ultrasonic image.
  • the solid line is the signal or data flow related to the optical image
  • the dashed line is the signal or data flow related to the ultrasonic image
  • the two-dot chain line is the insertion part 2.
  • the thick line is the signal or data related to the 3D image Indicates the flow of data
  • the dotted line indicates the flow of matching information
  • the curved arrow indicates the flow of other signals or data.
  • the ultrasonic transducer 22 receives the excitation signal in the form of a pulse voltage from within the ultrasonic observation device 3 and converts it into an ultrasonic beam that is a compression wave of the medium.
  • the ultrasonic beam propagates through the ultrasonic transmission medium 28 and the distal end tip 27 and is radiated to the outside of the ultrasonic endoscope 2, and the reflection echo from the inside of the subject is opposite to the ultrasonic beam. Return to ultrasonic transducer 2 2 along the path o
  • the ultrasonic transducer 22 converts the reflected echo into an electric echo signal and transmits it to the ultrasonic observation device 3 through a path opposite to the excitation signal. Furthermore, while this action is repeated repeatedly, the flexible shaft 29 and the ultrasonic vibrator 22 are rotated by the motor 23 in the operation unit 25, respectively, so that the arrow shown by the arrow in FIG. Rotate in the direction. Therefore, the ultrasonic beam is sequentially and radially irradiated on a plane (hereinafter referred to as a radial scanning plane) orthogonal to the axial direction of the insertion section 21 of the ultrasonic endoscope 2, thereby realizing a so-called mechanical radial scanning. .
  • a radial scanning plane orthogonal to the axial direction of the insertion section 21 of the ultrasonic endoscope 2
  • the ultrasonic observation device 3 performs envelope detection, logarithmic amplification, AZD conversion, scan conversion (radio (A process of converting an image in the polar coordinate system generated by scanning in the horizontal coordinate system into an image in the orthogonal coordinate system), and the like. To build). This ultrasonic image is output to the ultrasonic image memory 41 in the image processing device 6.
  • the CCD camera 31 generates an imaging signal based on imaging information of the inner surface of the body cavity. Specifically, the received light is converted into an electric imaging signal. — Then, the imaging signal is output to the optical observation device 4.
  • the optical observation device 4 constructs image data of an optical image (hereinafter, simply referred to as an optical image) based on the imaging signal. This optical image is output to an optical image memory 47 in the image processing device 6.
  • the coil drive circuit 11 sequentially outputs a coil excitation signal to the transmission coil 24.
  • the transmission coil 24 generates a magnetic field in space.
  • the receiving coil 12 sequentially detects the magnetic field and outputs an electric reception signal to the position calculation circuit 13.
  • the position calculation circuit 13 calculates a position direction error based on the received signal and outputs the calculated position direction data to the image processing device 6.
  • This position / direction data is a data including the position and direction of the transmission coil 24 with respect to the reception coil unit 10. More specifically, the position of the transmission coil 24
  • the direction of the insertion axis of the ultrasonic endoscope 2 (the direction indicated as “normal direction” in FIG. 2) and a specific direction parallel to the ultrasonic image (“ Direction ”).
  • the direction of the insertion axis of the ultrasonic endoscope 2 is the normal direction of the ultrasonic image.
  • the ultrasonic observation apparatus 3 of the present embodiment creates an ultrasonic image such that the 12 o'clock direction in FIG. 2 is the 12 o'clock direction of the ultrasonic image. Therefore, in the end, the position direction-de-time includes the direction indicating the normal direction of the ultrasonic image and the direction of ⁇ 2 o'clock.
  • the operator causes the controller 56 to switch the switch 55 via the keyboard 8 or the mouse 9.
  • the output destination of the position / direction data is set in the ultrasonic image memory 41.
  • this scanning method is referred to as manual scanning.
  • this scanning method is referred to as manual scanning.
  • this scanning method is referred to as manual scanning.
  • a plurality of continuous ultrasonic images 62 are obtained.
  • the region of interest 61 is included in most ultrasonic images.
  • the ultrasonic observation device 3 sequentially outputs the ultrasonic images generated in this way to the ultrasonic image memory 41. Controllers 5 and 6
  • the controller 56 stores the position and direction data as a header of the image data of the ultrasonic image or as a whole data. Due to recent advances in digital technology, the ultrasonic observation device 3 can construct an ultrasonic image with almost no delay with respect to radial scanning, and the recent position detection device 5 has a position direction with little delay with respect to transmission of a magnetic field. Since data can be calculated, the ultrasonic image memory 41 stores, in effect, each ultrasonic image and the position / direction data at the moment when the echo signal is acquired.
  • the three-dimensional data construction circuit 42 reads a plurality of continuous ultrasonic images from the ultrasonic image memory 41 and averages overlapping portions of each, or performs interpolation processing between the ultrasonic images to obtain an address.
  • a three-dimensional image data expressed in three-dimensional rectangular coordinates is created and output to the three-dimensional image memory 43.
  • the three-dimensional image data 63 consists of cells 64 whose addresses are represented by three-dimensional rectangular coordinates, and each cell 64 represents the luminance obtained based on the echo signal. Have value overnight.
  • the cross-section extraction circuit 44 extracts a large number of cells 64 corresponding to an appropriate plurality of cross-sections from the three-dimensional image data 63, and generates cross-section image data (hereinafter referred to as cross-section image data). I do.
  • the cross-sectional image data is output to the cross-sectional image memory 45 and stored. It should be noted that the position and direction of the cross section are set in advance by the operator via the keyboard 8 and the mouse 9, and in the present embodiment, a plurality of cross sections perpendicular to each other are set for convenience of explanation. I do.
  • the surface extraction circuit 46 cuts the three-dimensional image data 63 into a parallel cross-sectional image (hereinafter, referred to as a parallel slice image data 65). Then, cells corresponding to the luminal surface are extracted from each of the parallel slice images 65.
  • a known processing method as described in detail in Japanese Patent Application Laid-Open No. H10-192 by the present applicant is used.
  • the surface extraction circuit 46 creates binarized surface shape data separately from the three-dimensional image data 63 by setting cells corresponding to the surface to 1 and cells corresponding to other than the surface to 0. Then, it outputs to the surface shape memory 57 in the shape matching circuit 50.
  • the operator causes the controller 56 to switch the switch 55 via the keyboard 8 or the mouse 9.
  • the output destination of the position / direction data is set in the optical image memory 47.
  • the operator moves the insertion section 21 of the ultrasonic endoscope 2 while capturing an optical image so that the region of interest 61 is captured at various angles as shown in FIG.
  • the observation device 4 outputs the optical images thus generated to the optical image memory 47 one after another.
  • the controller 56 stores each optical image in the optical image memory 47 in association with the position and orientation at the moment when the optical image is input. For example, the controller 56 stores the position and orientation data as a header of the image data of the optical image or data of the image data.
  • the optical observation device 4 can construct an optical image with almost no delay from the image captured by the CCD camera 31 and the position detection device 5 transmits the magnetic field. Since the position and orientation data can be calculated with almost no delay for the optical image, the optical image memory 47 effectively stores each optical image and the position and orientation data at the moment when it was captured. become.
  • the surface shape estimating circuit 48 reads a plurality of continuous optical images from the optical image memory 47 and estimates the surface shape.
  • 'A method for estimating the surface shape uses a known processing method as described in detail in Japanese Patent Application Laid-Open No. 11-295618 by the present applicant.
  • the position and direction of the tip of the input unit 21 are detected using the position detection device 5 in the same manner as in the present application, and further, the subject is detected using the optical image from the CCD camera 31. It is a method to estimate the surface shape with high accuracy o
  • the surface shape estimating circuit 48 creates a binarized surface shape data by setting the cells corresponding to the surface to 1 and the cells corresponding to other than the surface to 0, and generates the surface shape data in the shape matching circuit 50. Output to shape memory 58.
  • the conceptual diagram of this surface profile is the same as that of Fig. 8 described above.
  • the surface shape estimating circuit 48 creates three-dimensional image data of the lumen surface separately from the surface shape data by mapping the original brightness value of the optical image to the surface shape, 3D image Output to memory 49.
  • the conceptual diagram of the three-dimensional image data is the same as FIG. 6 described above.
  • the three-dimensional image data is composed of cells whose addresses are represented by three-dimensional orthogonal coordinates. Each cell represents the luminance value of the lumen surface obtained from the imaging signal as R (red), G (green), B O with (blue) as data
  • the shape matching circuit 50 compares the surface shape data obtained from the ultrasonic image in the surface shape memory 57 with the surface shape data obtained from the optical image in the surface shape memory 58 to obtain the optical data. Calculate how the surface shape data obtained from the image is rotated, translated, and enlarged / reduced to best match the surface shape data obtained from the ultrasonic image. This is shown in Figure 10.
  • the cross-correlation circuit 59 applies rotation, translation, and enlargement / Z reduction to the surface shape data 72 obtained from the optical image, Calculate the cross-correlation value F with the obtained surface shape data 7 1 and repeat this while changing the degree of rotation, translation, and scaling conversion minutely, thereby obtaining the rotation when the cross-correlation value is maximized.
  • Output to circuit 5 1 The following describes the analytical model. If f (X, y, z) is the surface shape data 71 obtained from the ultrasonic image, and g (X, y, ⁇ ) is the surface shape data obtained from the optical image 7 2, the respective functions are as follows. Take the value of
  • T x ( ⁇ ), ⁇ y ( ⁇ ), and ⁇ ⁇ ( ⁇ ) are rotation matrices around the ⁇ , y, and z axes, respectively.
  • x ', y ,, z 5 are each (x, y, z), ( ⁇ , ⁇ , ⁇ ), ( ⁇ ⁇ , ⁇ y, ⁇ ⁇ ), Ru written as Hino function Will be.
  • the cross-correlation value F is given by the following equation as a function of (, ⁇ , ⁇ ), ( ⁇ X, ⁇ y, ⁇ ), and ⁇ .
  • the value to be obtained satisfies the following and maximizes the cross-correlation value F ⁇ 0, ⁇ , ⁇ , (6 ⁇ ) ⁇ , ( ⁇ y) ⁇ , (dz) o, ⁇ ⁇
  • the circuit 59 repeats the calculations of the equations (3) and (4) while slightly changing the values of ⁇ , ⁇ , ⁇ ), ( ⁇ , ⁇ y, ⁇ ), and Hi in order to obtain these.
  • the cross-correlation circuit 59 outputs 0, 0o, 00, ( ⁇ x) o, dy) ⁇ , ( ⁇ ) ⁇ , ⁇ as the matching information to the coordinate conversion circuit 51.
  • These values are coordinate transformation parameters that give the coincidence between the surface shape data obtained from the ultrasonic image and the surface shape data obtained from the optical image. Next, the synthesis and display of a three-dimensional image will be described.
  • the coordinate transformation circuit 51 performs coordinate transformation on the three-dimensional image data 63 obtained from the optical image in the three-dimensional image memory 49.
  • the coordinate transformation circuit 51 converts the surface shape data 72 obtained from the optical image to the surface shape data 71 obtained from the ultrasonic image. Rotating, translating, enlarging / reducing the surface shape data obtained from the optical image in the same way as rotating and translating the 3D image data obtained from the optical image. Scale up / down.
  • the coordinate conversion circuit 51 converts the values of 0, ⁇ 0, ⁇ 0, ( ⁇ ⁇ ) ⁇ , ( ⁇ y) ⁇ , ( ⁇ )) ⁇ , and h ⁇ from the cross-correlation circuit 59.
  • the coordinate transformation process expressed by the equation (3) is applied to the three-dimensional image data 63 obtained from the optical image in the three-dimensional image memory 49.
  • the coordinate-transformed three-dimensional image data (hereinafter referred to as surface image data) is output to the surface image memory 52 and stored. By performing such processing, the coordinate system and the surface image data of the cross-sectional image data in the cross-sectional image memory 45 are obtained.
  • the coordinate system of the surface image in the memory 52 coincides.
  • the synthesizing circuit 53 synthesizes the cross-sectional image data and the surface image data, performs processing such as hidden surface elimination, and converts the surface 73 obtained from the optical image shown in FIG. 11 and the cross-section 7 obtained from the ultrasonic image. By compositing, a three-dimensional image in which the region of interest 61 is displayed is constructed.
  • This three-dimensional image is converted into a signal such as a video signal and output to the monitor 7 so that it can be output to the screen 7 a of the monitor 7 by the display circuit 54.
  • Monitor 7 Displays a three-dimensional image.
  • the position and direction of the optical image and the ultrasonic image can be accurately correlated and displayed without the need to replace the endoscope from the patient.
  • the examination time can be shortened, and the labor for maintenance of the endoscope before and after the examination such as cleaning and disinfection can be reduced, and the burden on the patient can be reduced.
  • the surface extraction circuit 46 cuts the three-dimensional image data into a parallel slice image 65 and extracts the surface, and then creates the surface shape data 66.
  • it may be configured to extract the surface from each ultrasonic image as it is and to interpolate the surface to create the surface shape data 66.
  • the surface extraction method described in detail in Japanese Patent Application Publication No. Can be Further, the surface extraction method may be any other known method.
  • the surface shape estimating circuit 48 is configured to estimate the surface shape by the method disclosed in Japanese Patent Application Laid-Open No. H11-29556 / 18.
  • the position and direction of the distal end of the insertion section 21 are detected, and optical images of the same subject from the CCD camera 31 at different times are used to detect a time difference.
  • Other methods for estimating the surface shape of the subject by teleo-vision may be used.
  • Fig. 12 is an electronic radial scanning type ultrasonic endoscope applicable to the first embodiment and a diagram for explaining the scanning thereof.
  • FIG. FIG. 13 is a diagram for explaining a convex scanning ultrasonic endoscope applicable to the first embodiment and its scanning.
  • FIG. 14 is a diagram for explaining a two-dimensional array type ultrasonic endoscope applicable to the first embodiment and scanning thereof.
  • a mechanical radial scanning type ultrasonic endoscope 2 that mechanically rotates an ultrasonic transducer 22 and performs radial scanning is used as an ultrasonic endoscope.
  • the present invention is not limited to this.
  • an electronic radial scanning ultrasonic endoscope 81 in which a strip-shaped ultrasonic transducer 82 shown in FIG. 12 is provided in an annular array with respect to the insertion axis, Place the ultrasonic vibrator 92 along the insertion axis shown in Fig. 13.
  • a Compex-scanning ultrasonic endoscope 91 provided in a rectangular shape may be used.
  • scanning in which the insertion portion 93 is twisted about the insertion axis as shown in FIG. 13 is performed instead of manual scanning.
  • the mirror 101 may be applied to the present invention.
  • two transmission coils 24 are provided independently, but this is wound around two axes as shown in FIGS. 12, 13, and 14.
  • the coil 84 may be integrated, and the form of the transmission coil may be variously applied. Furthermore, even if the transmission coil 24 and the reception coil 12 shown in FIG. Since the position and direction data of the tip can be calculated, there is no problem.
  • FIG. 15 is a block diagram showing a shape matching circuit according to the second embodiment of the present invention.
  • the second embodiment differs from the first embodiment in the configuration and operation of the shape matching circuit 110.
  • the shape matching circuit 110 includes a surface shape memory 57, a surface shape memory 58, a center of gravity calculation circuit 111 as a center of gravity calculation means, a center of gravity comparison circuit 112, and a translation circuit 113.
  • An inertia spindle calculation circuit 114 as inertia spindle calculation means, an inertia spindle comparison circuit 115, a rotation circuit 116, and a cross-correlation circuit 117 are provided.
  • the solid line in Fig. 15 shows the flow of signals or data related to the optical image
  • the broken line shows the flow of signals or data related to the ultrasonic image
  • the dotted line shows the flow of matching information.
  • a cross-correlation circuit performs rotation, translation, and scale conversion on a surface shape data obtained from an optical image, and performs a surface shape data obtained from a supersonic image. Calculates the cross-correlation value F, and repeats it while changing the degree of rotation, translation, and scaling in small increments, so that the Euler angles of rotation when the cross-correlation value is maximized, ⁇ , ⁇ ), The displacement ( ⁇ ⁇ , ⁇ y 5 ⁇ ⁇ ) and the enlargement / reduction ratio were calculated.
  • the rotation is based on the positional relationship of the principal axes of inertia of both surface shape data, and only the dog expansion / reduction uses the cross-correlation circuit 117, and ( ⁇ X, ⁇ y, d), calculate.
  • the center-of-gravity calculating circuit 111 reads out the surface shape data stored in both surface shape memories and calculates the position vector of the center of gravity of each.
  • i is a number assigned to a cell constituting the surface shape memory
  • r i is a position vector of each cell
  • I i is data of each cell (the surface is 1; otherwise, 0.).
  • the center-of-gravity comparison circuit 1 1 2 calculates the difference vector between the center-of-gravity position vector G calculated in the both-surface shape data and the position shift between the two-surface shape data, that is, the translational displacement ( ⁇ , ⁇ y, ⁇ ) are calculated. After that, the center of gravity comparison circuit 1 1 2
  • the translation circuit 1 13 applies translation (translation) processing to the surface shape data 7 obtained from the optical image in the surface shape memory 58 to obtain the surface shape data 7 obtained from the ultrasonic image. Align the center of gravity with 1 and output it to the rotation circuit 1 16 and the inertia spindle calculation circuit 1 14.
  • the main shaft inertia calculating circuit 114 reads out the surface shape data from the ultrasonic image stored in the surface shape memory 57 and calculates the unit vector of the main shaft.
  • the inertia spindle calculation circuit 1 14 calculates the inertia of the surface shape data from the optical image centered by the translation circuit 113.
  • the unit vector of the sex axis is also calculated.
  • the principal axis of inertia is a set of three orthogonal axes that are inherent in any rigid body and are usually treated in classical mechanics.
  • the main axis of inertia calculation circuit 114 regards the surface shape data as a set of cells represented by the luminance value Ii and the position vector ri, and further reads the surface shape data by reading the luminance value as mass. Assuming that it is a rigid body, the principal axis of inertia is calculated from the surface shape data just like a rigid body.
  • the principal axis of inertia calculation circuit 114 calculates the right-handed unit vector of three axes perpendicular to the principal axis of inertia for each of the surface shape data from the ultrasonic image and the surface shape data from the optical image.
  • Methods for calculating the principal axes of inertia are known in classical mechanics and linear algebra.
  • the inertia spindle comparison circuit 1 15 calculates the relationship between the unit vectors of the inertia spindle calculated from both surface shape data.
  • the relationship between the two is represented by a 3-by-3 orthogonal matrix, from which the Euler angles of rotation (, ⁇ , ⁇ ) are calculated. This value is the rotational deviation between the two surface shape data.
  • the inertia spindle comparison circuit 115 outputs this value to the coordinate conversion circuit 51 and the rotation circuit 116.
  • the rotation circuit 1 16 performs rotation processing on the surface shape data obtained from the optical image output from the translation circuit 1 13 and aligns the direction with the surface shape data obtained from the ultrasonic image.
  • Output to the cross-correlation circuit 1 17 The cross-correlation circuit 117 reads the surface shape data 71 obtained from the ultrasonic image from the surface shape memory 57. Then, the cross-correlation circuit 117 enlarges or reduces the surface shape data obtained from the optical image output from the rotation circuit 116, and takes a cross-correlation between the two surface shape data. Further, by repeating the enlargement or reduction while changing the magnification, a peak that maximizes the cross-correlation value is obtained and output to the coordinate conversion circuit 51.
  • the rotation boiler angles (, ⁇ , ⁇ ), the translational displacements (dx, ⁇ y, ⁇ ⁇ ), and the enlargement / reduction rates that maximize the cross-correlation value are all independent.
  • the translational displacement (dx, ⁇ y, ⁇ ) is calculated from the center of gravity calculation circuit 11 1 and the center of gravity comparison circuit 1 12 from the rotation Euler.
  • the angles ( ⁇ , ⁇ , ⁇ ) are calculated from the principal axis of inertia calculation circuit 114 and the principal axis of inertia comparison circuit 115, and the cross-correlation circuit 117 only calculates the enlargement / reduction ratio. Since cross-correlation processing is generally heavy, according to the second embodiment, processing can be performed at a higher speed than in the first embodiment.
  • FIG. 16 is a block diagram showing a shape matching circuit according to the third embodiment of the present invention.
  • the third embodiment differs from the second embodiment in the configuration and operation of the shape matching circuit 120. Only different points from the second embodiment will be described.
  • the shape matching circuit 120 has a new adjusting circuit 121.
  • the other configuration of the third embodiment is the same as that of the second embodiment.
  • the translational displacement (dX, ⁇ y, ⁇ ) output from the center-of-gravity comparison circuit 112 and the rotation Euler angle (,, ⁇ , ⁇ ) are directly output to the coordinate conversion circuit 51.
  • these outputs are output to the adjustment circuit 121 as coarse adjustment values.
  • the cross-correlation circuit 1 17 separates the Euler angles of rotation (3 ⁇ 4 ⁇ , ⁇ , ⁇ ) and the translational displacements (5 x, ⁇ y, ⁇ ⁇ ⁇ ) separately from the coarse adjustment values as in the first embodiment. Calculate the enlargement / reduction ratio. However, at this time, the cross-correlation.
  • the adjustment circuit 1 2 1 calculates an accurate value from the coarse adjustment value and the fine adjustment value for the rotation coiler angles ⁇ , ⁇ ) and the translational displacement (5 x, y, ⁇ ⁇ ), and calculates the coordinate conversion circuit 5. Output to 1. Further, the adjustment circuit 12 21 outputs the enlargement / reduction ratio ⁇ from the cross-correlation circuit 1 17 as it is.
  • the rough adjustment value is output to the coordinate conversion circuit 51 as matching information.
  • the surface shape data obtained from the optical image and the surface shape data obtained from the ultrasonic image are compared.
  • the range imaged in the body cavity may be slightly different, and these coarse adjustment values may not accurately represent the Euler angle of rotation and the translational displacement ⁇
  • the shape matching circuit 120 is configured in the second embodiment. In this case, more accurate matching information can be output. Further, by performing the coarse adjustment before the cross-correlation circuit 117 calculates the fine adjustment value, it is possible to perform the cross-correlation processing while changing the independent variable while limiting the change range of the independent variable. Processing can be performed faster than in the first embodiment.
  • FIG. 17 is a block diagram showing an image processing device according to the fourth embodiment of the present invention.
  • a mapping circuit 25 1 is provided in the image processing device 206 in place of the coordinate conversion circuit 51 of the first embodiment.
  • the three-dimensional image data 63 Although the surface of the three-dimensional image shown in FIG. 11 was used as it is, in the fourth embodiment, the surface shape data obtained from the ultrasonic image was replaced with the luminance of the three-dimensional image obtained from the optical image. The value is mapped to represent the surface. Specifically, it is as follows.
  • the mapping circuit 25 1 receives the surface shape data from the surface extraction circuit 46, the matching information from the shape matching circuit 50, and the 3D image data 6-3 from the 3D image memory 49. I do.
  • the surface shape data of these is obtained from an ultrasonic image, and the R (red), G (green) and B (blue) of the lumen surface obtained from the optical image Have in the data.
  • the mapping circuit 251 associates each cell of the surface shape data obtained from the ultrasonic image with the cell of the three-dimensional image data obtained from the optical image based on the matching information. Then, the luminance value of the three-dimensional image data obtained from the optical image is mapped to the surface shape data obtained from the ultrasonic image and output to the surface image memory 52.
  • FIGS. 18 and 19 relate to the fifth embodiment of the present invention.
  • FIG. 18 is a block diagram showing an image processing device
  • FIG. 19 is an explanatory diagram showing an image displayed on a monitor.
  • an associating circuit 353 is provided in place of the synthesizing circuit 53 of the first embodiment.
  • mapping circuit 3 5 3 the coordinates of the mouse cursor on the screen 7 a of the monitor 7 are changed from the controller 5 6 sequentially as the mouse force-sol coordinate value data by the operator's operation of the mouse 9. input.
  • the image processing apparatus 303 includes the cross-sectional image memory 45, the three-dimensional image memory 49, the surface image memory 52, and the coordinate conversion circuit 51 of the first embodiment. Instead, a parallel slice image memory 360 is provided.
  • the parallel slice image memory 360 stores all of the parallel slice image data created by the surface extraction circuit 46 o
  • the three-dimensional image shown in FIG. 11 is synthesized and displayed.
  • the original image of the optical image and the ultrasonic image is simultaneously displayed as it is, and both the images are displayed. Display corresponding points. This is shown in Figure 19.
  • the mapping circuit 35 3 shown in FIG. 18 displays an appropriate optical image 3 71 on the left side of the screen 7 a of the monitor 7 as shown in FIG.
  • the operator selects a desired one using the keyboard 8 or the mouse 9, and the operator selects the desired one using the mouse 9 and the mouse cursor 3 7 2 on the screen 7a of the monitor 7.
  • the controller 56 outputs the mouse cursor coordinate value data to the association circuit 353.
  • the operator designates one point on the optical image 371, by operating the mouse click or the like.
  • the mapping circuit 353 adds this point to this point.
  • the associating circuit 353 converts the parallel slice image data including the corresponding point to the marker 373 on the optical image 371 into parallel slices. Select and read from the image memory 360. Thereafter, the associating circuit 353 attaches a marker 375 to the corresponding point on the parallel slice image data, and displays the parallel slice image 374 on the right side of the monitor 7 screen. Other operations are the same as those of the first embodiment.
  • the surface created from the optical image 37 1 shown in FIG. 11 is combined with a three-dimensional image and displayed.
  • the resolution is the same as that of the optical image 37 1. It may be lower than the original image.
  • the optical image 371 can be observed as it is with the original image having a good resolution by the image processing device 36, and it can be compared with the section having the ultrasonic brightness value. can do.
  • the parallel slice image data is used as the ultrasonic data to be displayed.
  • the associating circuit 353 is connected to the marker 373 on the optical image 371.
  • the ultrasonic image of the original image closest to the corresponding point may be selected from the ultrasonic image memory 41.
  • the position and direction of the optical image and the ultrasonic image can be accurately correlated and displayed without replacing the endoscope from the patient. This shortens the time required for maintenance of the endoscope before and after examinations such as cleaning and disinfection, and reduces the burden on the patient.
  • an ultrasonic endoscope apparatus capable of displaying an optical image and an ultrasonic image in an accurate correspondence with each other in position and direction.

Landscapes

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

Abstract

本発明の超音波内視鏡装置は、被検体の光学像を取得する光学像取得装置と、前記被検体の超音波像を取得する超音波像取得装置と、前記光学像取得時の前記被検体に対する前記光学像取得装置の位置情報を取得する位置情報取得装置と、前記位置情報取得装置により得られた位置情報を基に前記光学像取得装置により得られた光学像と前記超音波像取得装置により得られた超音波像とをマッチングさせるマッチング回路とを備えたことを特徴とする。

Description

明 細 超音波内視鏡装置 技術分野
本発明は、 被検体の光学像と超音波像を取得する超音波内視鏡 装置に関する。 背景技術
従来、 患者の開口部から挿入して体腔内を超音波により走査す る超音波内視鏡が種々実用化されている。
この超音波内視鏡の先端には、 通常、 超音波による走査をして そのエコーを受信する超音波振動子と、 食道、 胃、 十二指腸、 大 腸などの管腔表面の光学像を観察するための光学観察窓が設けら れている。
ところで、 一般に、 この超音波内視鏡による超音波像は、 光学 像との空間的な位置関係を把握しにくいという問題がある。 この ため、 術者には超音波像として管腔のどの部位がモニタに表示さ れているのかわかりにく く、 超音波像により正しく患部を診断す ることには熟練を要していた。
超音波内視鏡で空間的な位置関係を把握しにくい理由は以下の 通りである。 .
第 1に超音波の走査面は光学像では目に見えない。
第 2に超音波の走査面は常に光学像の視野内におさまる訳では ない。
第 3に超音波により走査する際、 超音波を患部まで到達させる ために超音波媒体である脱気水を管腔内に充満させたり、 バル一 ンを用いたりするが、 これらが光学像の視野を妨げるだめ、 結局 、 超音波像による観察は光学像による観察とは時間を置いた別の 観察となっていた。
このことに対応して、 近年では、 別々に撮像された光学像と超 音波像とを合成して表示することで、 光学像と超音波像との空間 的な位置関係を正確に把握できる超音波内視鏡装置が、 特開 2 0
0 2 - 1 7 7 2 9号公報に提案されている。
さらに、 特開 2 0 0 2 — 1 7 7 2 9号公報に提案されている技 術を用いることにより、 水平方向への病変の広がりの診断と、 垂 直方向への病変の深達度診断とを同時に行うことができ、 超音波 内視鏡検査の診断能を向上させることができる。
光学像と超音波像とを合成する際には、 あらかじめ双方の画像 デ一夕の位置関係を対応づける必要がある。 そこで、 この特開 2 0 0 2 - 1 7 7 2 9号公報の技術では、 1つの方法として、 光学 像から被検体病変部表面の 3次元立体把握が可能な、 立体 3次元 内視鏡計測システムと、 3次元超音波像再構成システムとを組合 せ、 両者によって得られた病変部表面の 3次元形態画像情報をパ ターンマッチングさせる方法を提案している。
さらに、 この中で、 光学像から病変部表面の 3次元形態画像情 報を把握する方法として、 次の 2つの方法を提案している。
第 1の方法は、 三角測量を基本原理とした、 スリ ッ ト光投影を 用いた光切断'法 'による立体計測方法である。 これぼ'次のような手 法である。
まず、 レーザーを光源とする投影装置を用いてス リ ッ ト光を発 生させ、 物体面に照射すると、 その物体面にその形状のまま投影 され、 この物体面での光強度分布をスコープで撮影すると、 その 物体形状に応じて変化する輝度パターンが得られる。 この輝度パ 夕一ンを解析することで、 物体までの距離 Z、 そして、 X、 Y座 標を加えた 3次元位置を把握することができる。
第 2の方法は、 内視鏡に装着される広角レンズによる歪曲収差 を補正する仕組みを使う方法である。
このような従来の光学像から病変部表面の 3次元形態画像情報 を把握する技術において、 ス リ ッ ト光投影を用いた光切断法によ る立体計測方法では、 レーザー光を導光するための特別な光学系 を内視鏡に設ける必要がある。 また、 歪曲収差を補正する方法で は、 内視鏡に広角レンズを装着する必要があるため、 立体 3次元 内視鏡計測システムのための内視鏡には特別のものが必要になつ てしまい、 3次元超音波像再構成システムのための内視鏡とは別 体にならざるを得ない。
さらに、 立体計測方法は、 物体面が無地の面であることが前提 であり、 管腔表面のように不均一な階調をもつ面への投影では、 正確に物体面の 3次元位置を把握できるとは言いがたい。
従って、 このような技術による超音波内視鏡装置を用いて光学 像と超音波像とを合成して表示するためには、 一方のシステムで 患者を検査した後、 内視鏡を抜去し、 もう一度別の内視鏡を揷入 して検査し直さなければならない。 その結果、 この超音波内視鏡 装置による検査は検査時間を長引かせ、 洗浄消毒等の検査前後の 内視鏡メンテナンスの手間を倍加させるだけでなく、 患者への負 担を強いるというデメ リ ッ トが生じてしまう。 また、 光学像と超 音波像の位置関係を対応づける際の正確さに欠けていた。
本発明は、 前記事情に鑑みてなされたものであり、 患者から内 視鏡を差し替える必要なく、 光学像と超音波像の互いの位置と方 向とを正確に対応づけて表示するこ とのできる超音波内視鏡装置 を提供することを目的とする。 発明の開示
本発明の超音波内視鏡装置は、 被検体の光学像を取得する光学 像取得装置と、 前記被検体の超音波像を取得する超音波像取得装 置と、 前記光学像取得時の前記被検体に対する前記光学像取得装 置の位置情報を取得する位置情報取得装置と、 前記位置情報取得 装置により得られた位置情報を基に前記光学像取得装置により得 られた光学像と前記超音波像取得装置によ り得られた超音波像と をマッチングさせるマッチング回路とを備えたことを特徴とする
図面の簡単な説明
図 1は、 本発明の第 1の実施の形態に係る超音波内視鏡装置の 全体構成を示すプロック図である。
図 2は、 本発明の第 1の実施の形態に係る内視鏡の揷入部の挿 入側先端の拡大して示す断面図である。
図 3は、 本発明の第 1の実施の形態に係る画像処理装置を示す ブロック図である。
図 4は、 本発明の第 1の実施の形態に係るは形状マッチング回 路を示すプロック図である。
図 5は、 本発明の第 1の実施の形態に係る手引き走査の作用を 示す説明図である。
図 6は、 本発明の第 1の実施の形態に係る 3次元画像デ一夕を 示す概念図である。 図 7は、 本発明の第 1の実施の形態に係る 3次元画像データの 切り直しを示す説明図である。
図 8は、 本発明の第 1 の実施の形態に係る表面形状データを示 す説明図である。
図 9は、 本発明の第 1の実施の形態に係る内視鏡による関心領 域の撮像を示す説明図である。
図 1 0は、 本発明の第 1の実施の形態に係る光学像と超音波像 の合成を示す説明図である。
図 1 1は、 本発明の第 1の実施の形態に係るモニタに表示され る画像を示す説明図である。
図 1 2は、 本発明の第 1の実施の形態第 1の実施の形態に適用 可能な電子ラジアル走査型超音波内視鏡とその走査を説明するた めの図である。
図 1 3は、 本発明の第 1の実施の形態に適用可能なコンペック ス走査型超音波内視鏡とその走査を説明するための図である。
図 1 4は、 本発明の第 1の実施の形態第 1の実施の形態に適用 可能な 2次元アレイ型超音波内視鏡とその走査を説明するための 図である。
図 1 5は、 本発明の第 2の実施の形態に係る形状マッチング回 路を示すブロック図である。
図 1 6は、 本発明の第 3の実施の形態に係る形状マッチング回 路を示すブロック図である。
図 1 7は、 本発明の第 4の実施の形態に係る画像処理装置を示 すブロック図である。
図 1 8は、 本発明の第 5の実施の形態に係る画像処理装置を示 すブロック図である。
図 1 9は、 本発明の第 5の実施の形態に係るモニ夕に表示され る画像を示す説明図である。 発明を実施するための最良の形態
以下、 図面を参照して本発明の実施の形態を説明する。
(第 1の実施の形態)
図 1乃至図 1 1は本発明の第 1の実施の形態に係り、 図 1は超 音波内視鏡装置の全体構成を示すプロック図である。 図 2は内視 鏡の挿入部の挿入側先端の拡大して示す断面図である。 図 3は画 像処理装置を示すブロック図である。 図 4は形状マッチング回路 を示すブロック図である。 図 5は手引き走査の作用を示す説明図 である。 図 6は 3次元画像データを示す概念図である。 図 7は 3 次元画像デ一夕の切り直しを示す説明図である。 図 8は表面形状 デ一夕を示す説明図である。 図 9は内視鏡による関心領域の撮像 を示す説明図である。 図 1 0は光学像と超音波像の合成を示す説 明図である。 図 1 1はモニタに表示される画像を示す説明図であ る。
(構成)
図 1 に示すように、 本実施の形態の超音波内視鏡装置 1は、 超 音波内視鏡 2 と、 超音波観測装置 3 と、 光学観察装置 4 と、 位置 検出装置 5 と、 画像処理装置 6 と、 表示手段あるいは表示装置と してのモニ夕 7 と、 キ一ボード 8 と、 マウス 9 とを有している。 超音波内視鏡 2は、 挿入部 2 1 と、 操作部 2 5 とを連設したも のである。
挿入部 2 1は、 可撓性のある材質で構成され、 被検体の体腔内 へ揷入可能になっている。 操作部 2 5は、 揷入部 2 1の先端部の 超音波振動子 2 2を駆動するモ一夕 2 3を備えている。
挿入部 2 1の先端部には空間に磁場を励起する送信コィル 2 4 を設けている。
位置情報取得手段あるいは位置情報取得装置としての位置検出 装置 5は、 コイル駆動回路 1 1 と、 複数の受信コイル 1 2 と、 位 置算出回路 1 3 とを有している。
コイル駆動回路 1 1は、 送信コイル 2 4にコイル励振信号を出 力する。 複数の受信コィル 1 2は、 所定の配置方法で特定の位置 に固定され、 送信コイル 2 4が作り出す磁場を逐次検知して電気 的な受信信号を出力する。
位置算出回路 1 3は、 受信コイル 1 2が出力する受信信号から 揷入部 2 1の先端部の位置と方向を示すデータ (以下、 位置方向 デ一夕という) を算出する。
なお、 複数の受信コイル 1 2は直方体の筐体に一体的に固定さ れている。 以下、 この筐体と受信コイル 1 2 とを合わせて受.信コ ィルュニヅ ト 1 0 と呼ぶ。
図 1の中では受信コイル 1 2は紙面の都合上、 受信コイルュニ ッ ト 1 0の中で直線上に並べて固定されて表されているが、 実際 には 2次元平面上あるいは 3次元空間上に並べて固定されている ものとする。
図 2を用いて挿入部 2 1の挿入側先端を詳細に説明する。
図 2に示すように、 挿入部 2 1の先端部 2 6 にはポリメチルぺ ンテン等の材質でできた音響的に半透明な先端キヤップ 2 7が設 けられている。 先端キヤップ 2 7の内部には、 超音波像取得手段 あるいは超音波像取得装置としての超音波振動子 2 2が設けられ ており、 先端キヤップ 2 7内には超音波伝達媒体 2 8が充填され ている。 超音波振動子 2 2は可撓性のある材質で作られたフレキ シブルシャフ ト 2 9 に接続されている。 フレキシブルシャフ ト 2 9は図 1に示した内視鏡 2の操作部 2 5内のモ一夕 2 3の回転軸 へ接続されており、 図 2の矢印方向に回転するように設けられて いる。 超音波振動子 2 2はフレキシブルシャフ ト 2 9内の信号線 (図示せず) を介して内視鏡操作部 2 5経由で超音波観測装置 3 へエコー信号を出力する。
挿入部 2 1の先端には空間に磁場を作り出すソレノィ ドコイル である 2個の送信コイル 2 4が設けられている。 2個の送信コィ ル 2 4は、 信号線 3 0を介して位置検出装置 5内のコイル駆動回 路 1 1 と接続している。 この送信コィル 2 4のうちの一つは図 2 の 「 1 2時方向」 と書かれた方向を軸として、 もう一方は 「法線 方向」 と書かれた方向を軸として導線がコィル状に巻かれている この 「法線方向」 は挿入部 2 1 の挿入軸方向、 「 1 2時方向」 はこれに直交する方向である。 「法線方向」 は超音波振動子 2 2 がラジアル走査して得られる超音波像の法線方向と一致する。 ま た、 1 2時方向に卷かれた送信コイル 2 4は、 導線の卷かれてい る方向が、 法線方向に直交する方向のうち超音波像の 1 2時方向 に一致するよう設けられているものとする。 なお、 超音波振動子 2 2がラジアル走査を行う場合の作用は後述する。
これらの他に、 挿入部 2 1の先端部 2 6 には、 光学像をカラ一 で撮像するための電荷結合素子型固体撮像素子カメラ (以下、 C C Dカメラと呼ぶ) 3 1 と、 C C Dカメラ 3 1での撮像に必要な 光を体腔内に照射するための撮像光照射窓 3 2 とが設けられてい る。
光学像取得手段あるいは光学像取得装置としての C C Dカメラ 3 1は、 光学観察装置 4 と接続しており、 揷入部 2 1内の信号線 (図示せず) を介して内視鏡操作部 2 5経由で光学観察装置 4へ 撮像信号を出力する。 光学観察装置 4はこの撮像信号を基に体腔 内の光学像を作成する。 また、 撮像光照射窓 3 2は揷入部 2 1内 に設けられた光ファイバ等の導光路 (図示せず) を経由して光源 装置 (図示せず) からの撮像光が届く よう構成されており、 C C Dカメラ 3 1での撮像に必要な撮像光を体腔内に照射するごとが できる。
そして、 挿入部 2 1の先端部には、 上記した揷入部 2 1の先端 部の各部を一体的に保持するための硬性の硬性フレーム 3 3が図 2のように設けられている。
図 3を参照して図 1 に示した画像処理装置 6 の構成を詳細に説 明する。
図 3に示すように、 画像処理装置 6は、 超音波像メモリ 4 1 と 、 3次元データ構築回路 4 2 と、 大容量の第 1の 3次元画像メモ リ (以下、 単に 3次元画像メモリ 4 3 ) と、 断面抽出回路 4 4 と 、 断面画像メモリ 4 5 と、 表面形状演算手段あるいは表面形状演 算回路としての表面抽出回路 4 6 と、 光学像メモリ 4 7 と、 表面 形状演算手段あるいは表面形状演算回路としての表面形状推定回 路 4 8 と、 大容量の第 2の 3次元画像メモリ (以下、 単に 3次元 画像メモリ 4 9 ) と、 マッチング手段としての形状マッチング回 路 5 0 と、 座標変換回路 5 1 と、 表面画像メモリ 5 2 と、 合成回 路 5 3 と、 表示回路 5 4 と、 スイ ッチ 5 5 と、 これらを制御する コン トローラ 5 6 とを有している。
スィ ッチ 5 5は、 位置検出装置 5からのデ一夕の出力先を超音 波像メモリ 4 1 と光学像メモリ 4 7のうち一方に切り替える。 コ ン トロ一ラ 5 6は、 キ一ボード 8やマウス 9からの入力に応じて これら各部お'よび各回路を制御する。 ' - 次に、 図 4を用いて図 3に示した形状マヅチング回路 5 0の構 成を詳細に説明する。
図 4に示すように、 形状マッチング回路 5 0は、 第 1の表面形 状メモリ 5 7 と、 第 2の表面形状メモリ 5 8 と、 相互相関回路 5 9 とを設けている。
また、 本実施形態の位置検出装置 5は、 前記超音波像取得時の 前記被検体に対する前記超音波像取得手段の位置と方向の情報を 取得するようになっている。
(作用)
以下、 第 1の実施の形態の作用を説明する。
図 1、 図 3及び図 4において、 実線は光学像に関連する信号ま たはデ一夕の流れ、 破線は超音波像に関連する信号またはデ一夕 の流れ、 2点鎖線は挿入部 2 1の先端部の位置と方向に関連する 信号またはデ一夕の流れ、 太線は 3次元画像に関連する信号また はデ一夕の流れ、 点線はマッチング情報の流れ、 曲線の矢印はそ の他の信号またはデ一夕の流れを示している。
超音波像を構築する作用について説明する。
超音波振動子 2 2は、 超音波観測装置 3内からのパルス電圧状 の励起信号を受け取って媒体の疎密波である超音波ビームに変換 する。
超音波ビームは超音波伝達媒体 2 8 と先端キ ップ 2 7 とを伝 わり超音波内視鏡 2外部へと放射され、 被検体内からの反射ェコ 一が超音波ビームとは逆の経路を迪つて超音波振動子 2 2へ戻る o
超音波振動子 2 2は反射エコーを電気的なエコー信号に変換し て励起信号とは逆の経路で超音波観測装置 3へ伝達する。 さらに 、 この作用を反復的に繰り返す一方で、 操作部 2 5内のモー夕 2 3が回転することによ り フレキシブルシャフ ト 2 9 と超音波振動 子 2 2がそれそれ図 2のプロヅク矢印の方向へ回転する。 このた め超音波ビームが超音波内視鏡 2の挿入部 2 1の軸方向に直交す る平面 (以下、 ラジアル走査面という) 内を順次放射状に照射さ れ、 いわゆるメカニカルラジアル走査が実現する。
以下、 メカニカルラジアル走査は単にラジアル走査と呼ぶ。 超音波観測装置 3は、 超音波振動子 2 2からのエコー信号に包 絡線検波、 対数増幅、 A Z D変換、 スキャンコンバート (ラジア ル走査で生成された極座標系の画像デ一夕を直交座標系の画像デ 一夕に変換する処理) 等の公知の処理を施して超音波像の画像デ —夕 (以下、 単に超音波像という) を構築する。 この超音波像は 画像処理装置 6内の超音波像メモリ 4 1 に出力される。
次に、 光学像を構築する作用について説明する。
C C Dカメラ 3 1は、 体腔内表面の撮像情報を基に撮像信号を 生成する。 具体的には受けた光を電気的な撮像信号へ変換する。— そして、 撮像信号を光学観察装置 4へ出力する。 光学観察装置 4 はこの撮像信号を基に光学像の画像データ (以下、 単に光学像と いう) を構築する。 この光学像は画像処理装置 6内の光学像メモ リ 4 7に出力される。
次に、 挿入部 2 1の先端部の位置方向デ一夕を算出する作用に ついて説明する。
コィル駆動回路 1 1は、 送信コィル 2 4にコィル励振信号を逐 次出力する。 送信コイル 2 4は、 空間に磁場を生成する。
受信コイル 1 2は、 磁場を逐次検知して位置算出回路 1 3に電 気的な受信信号を出力する。
位置算出回路 1 3は、 受信信号を基に位置方向デ一夕を算出し 、 画像処理装置 6へ出力する。 この位置方向デ一夕は送信コイル 2 4の受信コィルュニッ ト 1 0に対する位置と方向とを含んだデ —夕とする。 具体的には、 位置方向デ一夕は送信コイル 2 4の位 置だけでなく、 超音波内視鏡 2の挿入軸方向 (図 2の 「法線方向 」 と示された方向) と、 超音波像に平行な特定の方向 (図 2の 「 1 2時方向」 と書かれた方向) とを含んでいるものとする。 ここ で、 超音波内視鏡 2の挿入軸方向は超音波像の法線方向である。
さらに、 本実施の形態の超音波観測装置 3は、 図 2の 1 2時方 向が、 超音波像の 1 2時方向になるよう超音波像を作成する。 従 つて、 結局、 位置方向—デ一夕は、 超音波像の法線方向と Γ 2時方 向を示すデ一夕を含むことになる。
次に、 画像処理装置 6の作用を説明する。
第 1 に、 超音波像に関連する信号あるいはデ一夕の流れについ て説明する。
まず、 術者はキ一ボード 8やマウス 9を介し、 コン トローラ 5 6にスイ ッチ 5 5を切り替えさせる。 ここでは、 位置方向デ一夕 の出力先が超音波.像メモリ 4 1 に設定される。
この後、 術者はラジアル走査をしながら超音波内視鏡 2の揷入 部 2 1をゆつ く り引き抜く。 すると図 5に示すように挿入部 2 1 の先端部がゆっ く り移動する (以下、 この走査方法を手引き走査 と呼ぶことにする) 。 手引き走査に伴い、 連続した複数の超音波 像 6 2が得られる。 図 5のように、 挿入部 2 1の先端部が常に閧 心領域 6 1 の近くにいるよう手引き走査を実行すれば、 ほとんど の超音波像内に関心領域 6 1が含まれる。 超音波観測装置 3 はこのように生成される超音波像を次々と超 音波像メモリ 4 1に出力する。 コン トローラ 5 6は各超音波像と
Ϊ
それが入力された瞬間の位置方向データとを関連付けて超音波像 メモリ 4 1 に記憶させるようにする。 例えば、 コン トローラ 5 6 は位置方向データを超音波像の画像データのヘッダもしくはフッ 夕のデ一夕として記憶させるようにする。 近年のデジタル技術の 進歩によ 、 超音波観測装置 3はラジアル走査に対して'ほとんど 遅延なく超音波像を構築でき、 また近年の位置検出装置 5は磁場 の送信に対してほとんど遅延なく位置方向データを算出できるの で、 超音波像メモリ 4 1 には、 事実上、 各超音波像とそのエコー 信号が取得された瞬間の位置方向デ一夕とが記憶されることにな る。
このようにして、 超音波像メモリ 4 1 には連続した複数の超音 波像が、 おのおのの位置方向データ と関連付けられて記憶される o
3次元データ構築回路 4 2は、 超音波像メモリ 4 1から連続す る複数の超音波像を読み出し、 各々が重複する部分を平均化した り、 超音波像間に補間処理を施してアドレスが 3次元の直交座標 で表現される 3次元画像デ一夕を作成し、 3次元画像メモリ 4 3 に出力する。
3次元画像デ一夕の概念を図 6を参照して説明する。 図 6に示すように、 3次元画像デ一夕 6 3はァドレスが 3次元 の直交座標で表現されるセル 6 4からなり、 各々のセル 6 4はェ コ一信号を基にして得た輝度値をデ一夕に持つ。
断面抽出回路 4 4は、 3次元画像データ 6 3のうち、 適当な複 数の断面に相当する多数のセル 6 4を抽出し、 断面の画像データ (以下、 断面画像デ一夕という) を作成する。
断面画像デ一夕は断面画像メモリ 4 5 に出力され、 記憶される 。 なお、 断面の位置や方向はあらかじめキ一ボード 8やマウス 9 を介して術者が設定するものとし、 本実施の形態では説明の都合 上、 互いに垂直な複数の断面が設定されているものとする。
表面抽出回路 4 6は、 図 7に示すように、 3次元画像デ一夕 6 3 を平行な断面像 (以下、 平行スライス像デ一夕 6 5 という) に 切り直す。 そして、 平行スライス像デ一夕 6 5の各々から管腔表 面に相当するセルを抽出する。 各平行スライス像デ一夕 6 5から 表面を抽出する方法は本出願人による特開平 1 0— 1 9 2号公報 に詳述されているような公知の処理方法を用いる。 この後、 表面 抽出回路 4 6は、 3次元画像デ一夕 6 3 とは別に、 表面に相当す るセルを 1、 表面以外に相当するセルを 0にして 2値化しだ表面 形状データを作成し、 形状マッチング回路 5 0内の表面形状メモ リ 5 7に出力する。
この表面形状デ一夕の概念を図 8を用いて詳細に説明する。 な お、 図 8では説明の都合上、 表面形状デ一夕 6 6の各セル 6 7は メ ッシュが粗く表現されているが、 実際には抽出された表面 6 8 (セルのデータが 1のセルから構成される) が図 8 に表現されて いる程度に滑らかに表現できるようメ ッシュが細かく切られてい る。
次に、 光学像に関連する信号及びデ一夕の流れについて説明す まず、 術者はキ一ボー ド 8やマウス 9を介し、 コン トローラ 5 6 にスイ ッチ 5 5を切り替えさせる。 ここでは、 位置方向デ一夕 の出力先が光学像メモリ 4 7に設定される。 この後、 術者は光学 像を撮像しながら超音波内視鏡 2の挿入部 2 1 を移動させ、 図 9 に示すように、 関心領域 6 1が様々な角度で撮像されるようにす 光学観察装置 4はこのように生成される光学像を次々と光学像 メモリ 4 7に出力する。
コン トローラ 5 6は各光学像とそれが入力された瞬間の位置方 向デ一夕とを関連付けて光学像メモリ 4 7に記憶させるようにす る。 例えば、 コン トローラ 5 6は位置方向データを光学像の画像 デ一夕のヘッダもしくはフッ夕のデータとして記憶させるように する。 光学観察装置 4は C C Dカメラ 3 1 による撮像に対してほ とんど遅延なく光学像を構築でき、 位置検出装置 5は磁場の送信 に対してほとんど遅延なく位置方向デ一夕を算出できるので、 光 学像メモリ 4 7には、 事実上、 各光学像とそれが撮像された瞬間 の位置方向デ一夕とが記憶されることになる。
このようにして、 光学像メモリ 4 7には、 連続した複数の光学 像が、 おのおのの位置方向デ一夕と関連付けられて記憶される。
表面形状推定回路 4 8は、 光学像メモリ 4 7から連続する複数 の光学像を読み出し、 表面形状を推定する。' この表面形状を推定 する方法は本出願人による特開平 1 1 — 2 9 5 6 1 8号公報に詳 述されているような公知の処理方法を用いる。 ここに開示されて いる処理方法は本願と同様に位置検出装置 5を使って揷入部 2 1 の先端部の位置及び方向を検出し、 さらに C C Dカメラ 3 1から の光学像を用いて被検体の表面形状を精度良く推定する方法であ o
この後、 表面形状推定回路 4 8は、 表面に相当するセルを 1、 表面以外に相当するセルを 0にして 2値化した表面形状デ一夕を 作成し、 形状マツチング回路 5 0内の表面形状メモリ 5 8に出力 する。 この表面形状デ一夕の概念図は先に述べた図 8 と同じであ る。
さらに、 表面形状推定回路 4 8は、 この表面形状にもともとの 光学像のカラ一の輝度値をマッピングすることで、 表面形状デ一 夕とは別に管腔表面の 3次元画像データを作成し、 3次元画像メ モリ 4 9 に出力する。 3次元画像データの概念図は先に述べた図 6 と同じである。 3次元画像デ一夕はアドレスが 3次元の直交座 標で表現されるセルからなり、 各々のセルは撮像信号から得た管 腔表面の輝度値を R (赤) 、 G (緑) 、 B (青) をデータに持つ o
次に、 形状マッチング回路 5 0の作用とマッチング情報の流れ について説明する。
形状マッチング回路 5 0は表面形状メモリ 5 7内の、 超音波像 から得た表面形状データ と、 表面形状メモリ 5 8内の、 光学像か ら得た表面形状デ一夕とを比較し、 光学像から得た表面形状デ一 夕をどのように回転、 並進、 拡大/縮小させたら最も良く超音波 像から得た表面形状データと一致するかを演算する。 この様子を 図 1 0に示す。
具体的には、 図 1 0に示すように、 相互相関回路 5 9が光学像 から得た表面形状デ一夕 7 2 に対し、 回転、 並進、 拡大 Z縮小の 変換を施して超音波像から得た表面形状デ一夕 7 1 との相互相関 値 Fを算出し、 回転、 並進、 拡大縮小の変換の程度を微小に変え ながらこれを繰り返すことで相互相関値が最大となる時の回転の オイラー角 ( , θ , φ ) 、 並進の変位 ((ί χ , δ y , δ ζ ) 、 拡大 縮小率 を算出する。 そして相互相関値が最大となる時の これらの値をマッチング情報として座標変換回路 5 1に出力する 以下に、 解析的なモデルを説明する。 f ( X , y , z ) を超音波像から得た表面形状データ 7 1、 g ( X , y , ζ ) を光学像から得た表面形状デ一夕 7 2 とすると各 々の関数は以下の値を取る。
1 ;(x,y,z) =表!)
f (x,y,z) = (1 )
,0 ;(x,y,z) = 人外
1 ; (Χ,Υ,Ζ) =^RE)
g (x,y,z) (2)
0 ;(x,y,z) タ卜
光学像から得た表面形状デ一夕 7 2に対し、 回転、 並進、 拡大 /縮小を施すと、 表面形状デ一夕 7 2上の点 ( x, y , ζ ) は以 下の式で点 (x, , y, , ζ, ) に座標変換される。
X' X
y = αΤχ(Ψ) Τγ(β)Τζ(Φ) y + sy (3) ζ' ζ
ン こ こで T x ( φ ) , Τ y ( θ ) , Τ ζ ( φ ) はそれそれ χ軸、 y軸、 z軸のまわりの回転行列とする。
式 ( 3 ) から、 x ' 、 y, 、 z 5 はそれぞれ (x , y, z ) 、 ( φ , θ , φ ) , ( δ χ , δ y , δ ζ ) , ひの関数として書かれ ることになる。
すると、 相互相関値 Fは、 ( , θ , φ ) 、 ( δ X , δ y , δ ζ ) 、 ひの関数として以下の式で与えられる。
F( »,Sy,8za) = ζ) f(x,y,z)-g(x'y'z') (4)
求める値は以下を満足して、 相互相関値 Fを最大にする ^ 0, θο , φο, ( 6 χ ) ο , ( δ y ) ο , ( d z ) o , ひ οであり、 ネ目 互相関回路 5 9はこれらを求めるために 、ψ、 θ , ø ) 、 ( δ χ , δ y , δ ζ ) 、 ひの値を微小に変えながら式 ( 3 ) と式 ( 4 ) の演算を繰り返す。
F ( φο , Θ 0 , ο , ( δ χ ) ο ( δ y ) ο , ( δ ζ ) ο , a ο;
= m a χ F ( 5 ) 結局、 相互相関回路 5 9は座標変換回路 5 1へマヅチング情報 として 0, 0 o, 0 0 , ( δ x ) o , d y ) ο , { δ ζ ) ο , ο を出力することになる。
これらの値が超音波像から得た表面形状データと光学像から得 た表面形状デ一夕との一致を与える座標変換パラメ一夕である。 次に、 3次元画像の合成と表示について説明する。
座標変換回路 5 1は 3次元画像メモリ 4 9内の光学像から得た 3次元画像データ 6 3に座標変換を施す。
この際、 座標変換回路 5 1は、 形状マッチング回路 5 0からの マッチング情報を基にして、 光学像から得た表面形状デ一夕 7 2 が超音波像から得た表面形状データ 7 1 と最も良く一致するよう 、 光学像から得た表面形状デ一夕 7 2を回転、 並進、 拡大 Ζ縮小 させるのと同じ方法で、 光学像から得た 3次元画像デ一夕 6 3を 回転、 並進、 拡大/縮小させる。
具体的には、 座標変換回路 5 1は相互相関回路 5 9からの 0 , Θ 0 , φ 0 , ( δ χ ) ο , ( δ y ) ο , ( δ ζ ) ο , ひ οの各値を 基に、 式 ( 3 ) で書かれる座標変換処理を 3次元画像メモリ 4 9 内の光学像から得た 3次元画像デ一夕 6 3 に施す。 座標変換され た 3次元画像デ一夕 (以下、 表面画像デ一夕という) は表面画像 メモリ 5 2へ出力され、 記憶される。 このように処理することで 、 断面画像メモリ 4 5内の断面画像デ一夕の座標系と表面画像メ モリ 5 2内の表面画像デ一夕の座標系が一致するのである。
合成回路 5 3は断面画像データと表面画像データを合成し、 陰 面消去等の処理を施して、 図 1 1 に示す光学像から得た表面 7 3 と超音波像から得た断面 7 とを合成して、 関心領域 6 1が表示 された 3次元画像を構築する。
この 3次元画像は表示回路 5 4でモニタ 7の画面 7 aに出力で きるよう、 ビデオ信号等の信号に変換され、 モニタ 7に出ガされ る。
モニタ 7ほこの 3次元画像を表示する。
(効果)
このような第 1の実施の形態によれば、 患者から内視鏡を差し 替える必要なく、 光学像と超音波像の互いの位置と方向とを正確 に対応づけて表示することができるので、 検査時間を短縮し、 洗 浄消毒等の検査前後の内視鏡メンテナンスの手間を減ら し、 患者 への負担を低減できる。
尚、 第 1の実施の形態では、 表面抽出回路 4 6が 3次元画像バ —夕を平行スライス像デ一夕 6 5に切り直して表面を抽出した後 、 表面形状デ一夕 6 6を作成するよう構成したが、 各超音波像か らそのまま表面を抽出してそれらを補間することで表面形状デ一 夕 6 6を作成するよう構成しても良い。 このように構成しても特 開平 1 0— 1 9 2号公報に詳述されている表面抽出の方法を用い ることができる。 また、 表面抽出の方法は他の公知のいかなる方 法でも良い。
また、 第 1の実施の形態では、 表面形状推定回路 4 8が特開平 1 1 - 2 9 5 6 1 8号公報に開示されている方法で表面形状を推 定するよう構成したが、 本方法のほかに、 挿入部 2 1の先端部の 位置及び方向を検出し、 さらに C C Dカメラ 3 1からの同一被検 体-に対する時刻の異なつ""た光学像を用いて、 いわば時差のあるス テレオ視で被検体の表面形状を推定する他の方法であつても良い 図 1 2は第 1の実施の形態に適用可能な電子ラジアル走査型超 音波内視鏡とその走査を説明するための図である。 図 1 3は第 1 の実施の形態に適用可能なコンべックス走査型超音波内視鏡とそ の走査を説明するための図である。 図 1 4は第 1の実施の形態に 適用可能な 2次元アレイ型超音波内視鏡とその走査を説明するた めの図である。
第 1の実施の形態では、 用いられる超音波内視鏡として超音波 振動子 2 2を機械的に回転してラジアル走査するメカニカルラジ アル走査型超音波内視鏡 2を用いたが、 本発明はこれに限定され るものではなく、 図 1 2に示す短冊状の超音波振動子 8 2を揷入 軸に対し環状にァレィ状に設けた電子ラジアル走査型超音波内視 鏡 8 1や、 図 1 3に示す挿入軸に沿って超音波振動子 9 2をァレ ィ状に設けたコンペックス走査型超音波内視鏡 9 1 を用いてもよ い。
コンベックス走査型超音波内視鏡 9 1 を用いた場合には、 手引 き走査に代えて、 図 1 3に示すよう挿入部 9 3を揷入軸を中心に ねじる走査 (ねじり走査) になる。
さらに、 図 1 4に示すように、 近年期待が高まっている超音波 振動子を平面にアレイ状に配した 2次元アレイ—超音波振動子 1 0 2を用いた 2次元アレイ型超音波内視鏡 1 0 1 を本発明に適用し てもよい。
2次元アレイ型超音波内視鏡 1 0 1を用いた場合には、 手引き 走査やねじり走査に代えて、 超音波振動子 1 0 2による走査のみ で、 平面状でなく 3次元的な走査 (ボリュ一ム走査) ができ、 一 度に 3次元の超音波像を取得することになる。 つま り、 このよう に構成することで、 図 2 に示した送信コイル 2 4を用いずとも超 音波による走査のみから 3次元画像デ一夕を構築できるので、 超 音波の走査の際には送信コィル 2 4が不要になる。
また、 図 2に示した第 1の実施の形態では、 送信コィル 2 4を 2個独立に設けたが、 これは図 1 2、 図 1 3、 図 1 4に示すよう 、 2軸に巻かれたコイル 8 4が一体になつたものであっても良く 、 送信コイルの形態は各種適用可能である。 さらには図 1 に示す 送信コィル 2 4 と受信コィル 1 2が逆であっても、 揷入部 2 1の 先端部の位置方向データは算出できるので、 一向に差し支えない 。 なお、 図 1 2、 図 1 3及び図 1 4では C C Dカメラ 3 1 と撮像 光照射窓 3 2 'とは省略されているが、 実際にはこれらの超音波内 視鏡 8 1 , 9 1 , 1 0 1でも設けられている。
(第 2の実施の形態)
以下、 図 1 5 を用いて、 第 2の実施の形態の超音波内視鏡装置 'の構成と作用とを説明する。
図 1 5は本発明の第 2の実施の形態に係る形状マツチング回路 を示すブロック図である。
図 1 5を用いた第 2の実施の形態の説明において、 図 1乃至図 1 1に示した第 1の実施の形態と同様の構成要素には同じ符号を 付して説明を省略している。
(構成)
図 1 5に示すように、 第 2の実施の形態は、 形状マッチング回 路 1 1 0の構成と作用が第 1の実施の形態と異なる。
形状マツチング回路 1 1 0は、 表面形状メモリ 5 7 と、 表面形 状メモリ 5 8 と、 重心算出手段としての重心算出回路 1 1 1 と、 重心比較回路 1 1 2 と、 並進回路 1 1 3 と、 慣性主軸算出手段と しての慣性主軸算出回路 1 1 4 と、 慣性主軸比較回路 1 1 5 と、 回転回路 1 1 6 と、 相互相関回路 1 1 7 とを設けている。
その他の構成は第 1の実施の形態と同じである。 (作用)
図 1 5の実線は光学像に関連する信号またはデ一夕の流れ、 破 線は超音波像に関連する信号またはデータの流れ、 点線はマツチ ング情報の流れを示している。
第 1の実施の形態では、 相互相関回路で、 光学像から得た表面 形状デ一夕に対し、 回転、 並進、 '拡大 縮小の変換を施して超音 波像から得た表面形状デ一夕との相互相関値 Fを算出し、 回転、 並進、 拡大縮小の変換の程度を微小に変えながらこれを繰り返す ことで相互相関値が最大となる時の回転のオイラー角 θ , φ ) 、 並進の変位 ( ί χ , δ y 5 δ ζ ) 、 拡大 Ζ縮小率ひを算出 した。
一方、 第 2の実施の形態では、 光学像から得た表面形状デ一夕 と超音波像から得た表面形状デ一夕の形状はほとんど同じである ことに着目し、 並進は両表面形状データの重心の位置関係から、 回転は両表面形状データの慣性主軸の位置関係から、 そして、 拡 犬/縮小のみ相互相関回路 1 1 7を使い、 ( δ X , δ y , d ζ ) 、 ひを算出する。
まず、 重心算出回路 1 1 1は、 両表面形状メモリに記憶された 表面形状デ一夕を読み出し、 それそれの重心の位置べク トルを算 出する。 重心の位置ベク トル Gの演算は以下の式で与えられる。 G =Zi Ii- ri
ここで iは表面形状メモリを構成するセルにつけられた番号で 、 r iは各セルの位置ベク トル、 I iば各セルのデータ (表面は 1 、 それ以外は 0.) である。 重心比較回路 1 1 2は、 両表面形状デ一夕で算出された重心の 位置ベク トル Gの差ベク トルを算出することで、 両表面形状デ一 夕間の位置ずれ、 すなわち並進の変位 ( χ, δ y , δ ζ ) を算 出する。 その後、 重心比較回路 1 1 2は、 この値を座標変換回路
■ /、 5 1 と並進回路 1 1 3に出力する。 J 並進回路 1 1 3は、 表面形状メモリ 5 8内の光学像から得た表 面形状データ 7 2に対し、 並進 (平行移動) の処理を施し、 超音 波像から得た表面形状データ 7 1 と重心を合わせ、 回転回路 1 1 6 と慣性主軸算出回路 1 1 4に出力する。 慣性主軸算出回路 1 1 4は、 表面形状メモリ 5 7に記憶された 超音波像からの表面形状デ一夕を読み出し、 その慣性主軸の単位 ベク トルを算出する。 また、 慣性主軸算出回路 1 1 4は、 並進回 路 1 1 3で重心を合わせられた光学像からの表面形状デ一夕の慣 性主軸の単位ベク トルをも算出する。 慣性主軸とは、 通常古典力 学で扱われる、 どのような剛体にも固有に存在する 1組の直交 3 軸である。
そこで、 慣性主軸算出回路 1 1 4は、 表面形状デ一夕を輝度値 I i、 位置ベク トル r iで表されるセルの集合と見なし、 さらに輝 度値を質量と読み替えることで表面形状データを剛体と見なして 、 表面形状デ一夕から剛体と同様に慣性主軸を—算出する。
ここで、 慣性主軸算出回路 1 1 4は、 超音波像からの表面形状 データと、 光学像からの表面形状データのそれそれについて慣性 主軸の直交する 3軸の右手系の単位ベク トルを算出する。 慣性主 軸の算出方法は古典力学及び線形代数で公知である。
慣性主軸比較回路 1 1 5は、 両表面形状デ一夕で算出された慣 性主軸の単位べク トル間の関係を算出する。 両者間の関係は 3行 3列の直交行列で表現され、 ここから回転のオイラー角 ( , Θ , φ ) が算出される。 この値が両表面形状データ間の回転ずれで ある。 その後、 慣性主軸比較回路 1 1 5は、 この値を座標変換回 路 5 1 と回転回路 1 1 6 に出力する。
回転回路 1 1 6は、 並進回路 1 1 3から出力された光学像から 得た表面形状デ一夕に対し、 回転の処理を施し、 超音波像から得 た表面形状デ一夕と方向を合わせ、 相互相関回路 1 1 7へ出力す o 相互相関回路 1 1 7は、 表面形状メモリ 5 7から超音波像から 得た表面形状データ 7 1を読み出す。 そして、 相互相関回路 1 1 7は、 回転回路 1 1 6から出力された光学像から得た表面形状デ 一夕を拡大または縮小し、 両表面形状データの相互相関を取る。 さらに拡大または縮小を倍率ひを変えながら繰り返すことで、 相 互相関値が最大になるひを求め、 座標変換回路 5 1へ出力する。
その他め作用は第 1の実施の形態と同じである。
(効果)
第 1の実施の形態では、 マッチング情報として相互相関値が最 大になる回転のォイラ一角 ( , θ , φ ) 、 並進の変位 ( d x , δ y , δ ζ ) 、 拡大/縮小率 を全て独立変数と見なして算出し ていたが、 第 2の実施の形態では、 並進の変位 ( d x , δ y , δ ζ ) は重心算出回路 1 1 1 と重心比較回路 1 1 2 とから、 回転の オイラー角 (^ , θ , φ ) は慣性主軸算出回路 1 1 4 と慣性主軸 比較回路 1 1 5 とから算出され、 相互相関回路 1 1 7が算出する のは拡大ノ縮小率ひのみである。 相互相関は一般に処理が重いた め、 第 2の実施の形態によれば、 第 1の実施の形態に比べて処理 を高速に実施できる。
第 2の実施の形態のその他の効果は第 1の実施の形態と同じで ある。
(第 3の実施の形態) 以下、 図 1 6を用いて、 第 3の実施の形態の超音波内視鏡装置 の構成と作用とを説明する。
図 1 6は本発明の第 3の実施の形態に係る形状マツチング回路 を示すブロック図である。
図 1 6を用いた第 3の実施の形態の説明において、 図 1 5 に示 した第 2の実施の形態と同様の構成要素には同じ符号を付して説 明を省略している。
(構成)
図 1 6に示すように、 第 3の実施の形態は、 第 2の実施の形態 と形状マツチング回路 1 2 0の構成と作用とが異なる。 第 2の実 施の形態とは異なる個所のみ説明する。
形状マッチング回路 1 2 0は、 調整回路 1 2 1を新たに設けて いる。
その他の第 3の実施の形態の構成は第 2の実施の形態と同じで ある。
(作用)
第 2の実施の形態は、 重心比較回路 1 1 2の出力である並進の 変位 ( d X , δ y , δ ζ ) と、 慣性主軸比較回路 1 1 5の出力で ある回転のオイラー角 ( , θ , φ ) を直接座標変換回路 5 1へ 出力するよう構成していた。 しかし、 第 3の実施の形態ではこれ らの出力を粗調整値として調整回路 1 2 1へ出力させる。 相互相関回路 1 1 7は、 この粗調整値とは別に第 1の実施の形 態と同じく回転のオイラー角 (¾^ , θ , φ ) 、 並進の変位 ( 5 x , δ y , δ ζ ) 、 拡大/縮小率ひを算出する。 ただし、 このとき 相互相関.回路 1 1 7で算出される回転のオイラー角 ( , θ , φ ) と並進の変位 ( 5 χ, δ y , δ ζ ) に関しては、 重心比較回路 1 1 2、 慣性主軸比較回路 1 1 5での調整に加えた再度の調整値 となり微調整値である 9
調整回路 1 2 1は、 回転のォイラ一角 θ , φ ) 、 並進の 変位 ( 5 x , y , δ ζ ) については粗調整値と微調整値とから 、 正確な値を算出し座標変換回路 5 1へ出力する。 また、 調整回 路 1 2 1は、 相互相関回路 1 1 7からの拡大 縮小率 αをそのま ま出力する。
(効果)
第 2の実施の形態では、 マツチング情報として粗調整値を座標 変換回路 5 1へ出力していたが、 光学像から得られた表面形状デ —夕と超音波像から得られた表面形状データとでは、 体腔内で撮 像された範囲が微妙に異なる場合があり、 これらの粗調整値が回 転のオイラー角と並進の変位を正確に表していない可能性がある ο
第 3の実施の形態では、 相互相関回路 1 1 7で微調整値を算出 するよう構成したので、 形状マツチング回路 1 2 0は第 2の実施 の形態に比べて正確なマツチング情報を出力できる。 さらに相互 相関回路 1 1 7が微調整値を算出する前に粗調整を実施すること で、 独立変数を変えながら行う相互相関処理を、 独立変数の変化 範囲を限定して行うことができ、 第 1の実施の形態に比べて処理 を高速に実施できる。
その他の効果は第 1の実施の形態と同じである。
(第 4の実施の形態) '
以下、 図 1 7を用いて、 第 4の実施の形態の超音波内視鏡装置 の構成と作用とを説明する。
図 1 7は本発明の第 4の実施の形態に係る画像処理装置を示す ブロック図である。
図 1 7を用いた第 4の実施の形態の説明において、 図 1乃至図 1 1 に示した第 1の実施の形態と同様の構成要素には同じ符号を 付して説明を省略している。
(構成)
図 1 7に示すように、 第 4の実施の形態は、 画像処理装置 2 0 6 において、 第 1の実施の形態の座標変換回路 5 1の代わりにマ ヅピング回路 2 5 1 を設けている。
その他の構成は第 1の実施の形態と同じである。
(作用)
第 1の実施の形態では光学像から得た 3次元画像デ一夕 6 3を そのまま用いて図 1 1に示す 3次元画像の表面を表現していたが 、 第 4の実施の形態では超音波像から得た表面形状データに光学 像から得た 3次元画像デ一夕の輝度値をマッビングして表面を表 現する。 具体的には以下の通りである。
マッピング回路 2 5 1には、 表面抽出回路 4 6からの表面形状 デ一夕、 形状マツチング回路 5 0からのマヅチング情報、 3次元 画像メモリ 4 9からの 3次元画像デ一夕 6— 3が入力する。 このう ちの表面形状データは超音波像から得られており、 3次元画像デ 一夕は光学像から得た管腔表面の R (赤) 、 G (緑) 、 B (青) を輝度値のデータに持つ。
マッピング回路 2 5 1は、 超音波像から得た表面形状デ一夕の 各セルに対し、 マッチング情報を基にして、 光学像から得た 3次 元画像データのセルを対応づける。 そして、 光学像から得た 3次 元画像データの輝度値を超音波像から得た表面形状デ一夕にマッ ビングして、 表面画像メモリ 5 2へ出力する。
その他の作用は第 1の実施の形態と同じである。
(効果)
このような第 4の実施の形態では、 第 1の実施の形態と同様の 効果が得られる。
(第 5の実施の形態)
以下、 図 1 8及び図 1 9 を用いて、 第 5の実施の形態の超音波 内視鏡装置の構成と作用とを説明する。
図 1 8及ぴ図 1 9は本発明の第 5の実施の形態に係り、 図 1 8 は画像処理装置を示すブロック図、 図 1 9はモニタに表示される 画像を示す説明図である。
図 1 8及び図 1 9を用いた第 5の実施の形態の説明において、 図 1乃至図 1 1 に示した第 1の実施の形態と同様の構成要素には 同じ符号を付して説明を省略している。
(構成)
図 1 8に示すように、 第 5の実施の形態は、 画像処理装置 3 0 6において、 第 1の実施の形態の合成回路 5 3の代わり に対応付 け回路 3 5 3を設けている。
対応付け回路 3 5 3には、 術者のマウス 9操作により、 モニタ 7の画面 7 a上のマウスカーソルの座標値が、 マウス力一ソル座 標値デ一夕として、 逐次コン トローラ 5 6から入力する。
また、 第 5の実施の形態は、 画像処理装置 3 0 6 において、 第 1の実施の形態の断面画像メモリ 4 5、 3次元画像メモリ 4 9、 表面画像メモリ 5 2、 座標変換回路 5 1 の代わりに平行スライス 像メモリ 3 6 0を設けている。 平行スライス像メモリ 3 6 0は表 面抽出回路 4 6が作成する平行スライス像デ一夕を全て保存する o
その他の構成は第 1の実施の形態と同じである。 (作用)
第 1の実施の形態では図 1 1に示す 3次元画像を合成して表示 していたが、 第 5の実施の形態では光学像と超音波像の原画像を ' そのまま同時表示し、 双方の対応点を表示する。 この様子を図 1 9に示す。
具体的には以下の通りである。
図 1 8に示す対応付け回路 3 5 3は、 図 Γ 9に示ずように 適 当な光学像 3 7 1をモニタ 7の画面 7 aの左側に表示させる。 こ の光学像 3 7 1 に対して、 術者がキ一ボード 8やマウス 9で所望 のものを選んで、 術者がマウス 9でモニタ 7の画面 7 a上のマウ スカ一ソル 3 7 2を動かすと、 コン トロ一ラ 5 6はマウスカーソ ル座標値データを対応付け回路 3 5 3へ出力する。
次に、 術者は、 マウスク リ ック等の操作により、 光学像 3 7 1 上で 1点を指定する。 対応付け回路 3 5 3はこの点にマ一力 3 7 3を付す。
次に対応付け回路 3 5 3は形状マッチング回路 5 0からのマツ チング情報を基に、 光学像 3 7 1上のマーカ 3 7 3への対応点を 含んだ平行スライス像デ一夕を平行スライス像メモリ 3 6 0の中 から選択して読み出す。 その後、 対応付け回路 3 5 3は平行スラ ィス像デ一夕上の対応点にマーカ 3 7 5を付して、 平行スライス 像 3 7 4をモニタ 7の画面右側に表示する。 その他の作用は第 1の実施の形態と同じである。
(効果)
第 1の実施の形態によれば、 図 1 1 に示す光学像 3 7 1から作 成した表面を 3次元画像に合成して表示していたが、 この場合、 解像度は光学像 3 7 1の原画像よ り も落ちる可能性がある。 第 5 の実施の形態は、 画像処理装置 3 0 6により、 光学像 3 7 1を解 像度の良い原画像のままで観察'でき、 なおかゥ超音波の'輝度値を 持つ断面と対比することができる。
(変形例)
第 5の実施の形態では、 表示する超音波のデ一夕として平行ス ライス像デ一夕を用いるよう構成したが、 対応付け回路 3 5 3が 光学像 3 7 1上のマーカ 3 7 3への対応点に最も近い原画の超音 波像を超音波像メモリ 4 1から選び出すよう構成しても良い。 こ のように構成することで、 モニタ 7の画面上に表示される画像は どちらも原画像になり、 さらに同原画像の位置関係を対応づけな がら全く劣化なく観察することができる。
以上述べた様に本発明によれば、 患者から内視鏡を差し替える ことなく、 光学像と超音波像の互いの位置と方向とを正確に対応 づけて表示することができるので、 検査時間を短縮し、 洗浄消毒 等の検査前後の内視鏡メンテナンスの手間を減ら し、 患者への負 担を低減できる。 以上、 本発明の実施形態について説明したが、 上記実施形態に 限定されるものではなく、 本発明の精神を逸脱しない範囲で幾多 の変化がなしえることは勿論である。 産業上の利用可能性
以上説明したように本発明によれば、 光学像と超音波像の互い の位置と方向とを正確に対応づけて表示するごとのできる超音波 内視鏡装置を提供することができる。

Claims

請求の範囲
1 . 被検体の光学像を取得する光学像取得装置と、
前記被検体の超音波像を取得する超音波像取得装置と、 前記光学像取得時の前記被検体に対する前記光学像取得装置の 位置情報を取得する位置情報取得装置と、
前記位置情報取得装置によ り得られた位置情報を基に前記光学 像取得装置により得られた光学像ど前記超音波像取得装置により 得られた超音波像とをマッチングさせるマッチング回路と、
を備えたことを特徴とする超音波内視鏡装置。
2 . 前記位置情報を基に前記光学像から第 1の表面形状を算出す る第 1の表面形状演算回路と、
前記超音波像から第 2の表面形状を算出する第 2の表面形状演 算回路と、
を更に有し、
前記マッチング回路は、 前記第 1および第 2の表面形状を用い てマッチングすることを特徴とする請求項 1 に記載の超音波内視
3 . 前記位置情報取得装置が、 前記超音波像取得時の前記被検体 に対する前記超音波像取得装置の位置情報を取得し、
前記第 2の表面形状演算回路が、 前記位置情報を基に第 2の表 面形状を算出すること、 を特徴とする請求項 2に記載の超音波内視鏡装置。
4 . 前記光学像の輝度値を表面画像デ一夕に用い、 前記超音波像 の輝度値を断面画像デ一夕に用いて、 前記表面画像データと前記 断面画像デ一夕とを合成して 3次元画像を構築する合成回路を設 けたことを特徴とする請求項 1 に記載の超音波内視鏡装置。
5 . 前記光学像から得られた画像と前記超音波像から得られた画 像を同時表示す-る表示装-置と、'
一方の画像上の任意点に対する、 他方への対応点を算出する対 応付け制御回路と、
を設けたことを特徴とする請求項 1に記載の超音波内視鏡装置
6 . 前記マッチング回路は、 相互相関処理を用いて前記マツチン グを行う ことを特徴とする請求項 2 に記載の超音波内視鏡装置。
7 . 前記マッチング回路には、 前記第 1の表面形状と前記第 2の 表面形状の重心を算出する重心算出回路が設けられていることを 特徴とする請求項 2に記載の超音波内視鏡装置。
8 . 前記マッチング回路には、 前記第 1の表面形状と前記第 2の 表面形状の慣性主軸を算出する慣性生軸算出回路が設けられてい ることを特徴とする請求項 2 に記載の超音波内視鏡装置。
9 . 前記マッチング回路には、 前記第 1の表面形状と前記第 2の 表面形状の重心を算出する重心算出回路、 も しくは前記第 1の表 面形状と前記第 2の表面形状の慣性主軸を算出する慣性主軸算出 回路を設け、
前記相互相関処理の前に、 前記重心も しくは前記慣性主軸を算 出したことを特徴とする請求項 6に記載の超音波内視鏡装置。
1 0 . 前記超音波像取得装置は、 それ自身でボリューム走査を行 うことを特徴とする請求項 1 に記載の超音波内視鏡装置。
1 1 . 前記超音波像取得装置は、 超音波振動子を 2次元アレイ状 に配した 2次元アレイ超音波振動子であることを特徴とする請求 項 1 0に記載の超音波内視鏡装置。
1 2 . 被検体の光学像を取得する光学像取得手段と、
前記被検体の超音波像を取得する超音波像取得手段と、 前記光学像取得時の前記被検体に対する前記光学像取得手段の 位置情報を取得する位置情報取得手段と、
前記位置情報取得手段により得られた位置情報を基に前記光学 像取得手段により得られた光学像と前記超音波像取得手段によ り 得られた超音波像とをマッチングさせるマッチング手段と、
を備えたことを特徴とする超音波内視鏡装置。
1 3 . 前記位置情報を基に前記光学像から第 1の表面形状を算出 する第 1の表面形状演算手段と、
前記超音波像から第 2の表面形状を算出する第 2の表面形状演 算手段と、 を更に有し、
前記マツチング手段は、 前記第 1および第 2の表面形状を用い 'てマッチングすることを特徴とする請求項 1 2に記載の超音波内 視鏡装置。
1 4 . 前記位置情報取得手段が、 前記超音波像取得時の前記被検 体に対する前記超音波像取得手段の位置情報を取得し、
前記第 2の表面形状'演算手段が、 前記位置情報を基に第 2'の表 面形状を算出することを特徴とする請求項 1 3 に記載の超音波内 1 5 . 前記光学像の輝度値を表面画像データに用い、 前記超音波 像の輝度値を断面画像データに用いて、 前記表面画像デ一夕と前 記断面画像デ一夕とを合成して 3次元画像を構築する合成手段を 設けたことを特徴とする請求項 1 2 に記載の超音波内視鏡装置。 1 6 . 前記光学像から得られた画像と前記超音波像から得られた 画像を同時表示する表示手段と、
一方の画像上の任意点に対する、 他方への対応点を算出する対 応付け制御手段と、
を設けたことを特徴とする請求項 1 2に記載の超音波内視鏡装 置。
1 7 . 前記マッチング手段は、 相互相関処理を用いて前記マヅチ ングを行うことを特徴とする請求項 1 3 に記載の超音波内視鏡装 置 o
1 8 . 前記マッチング手段には、 前記第 1の表面形状と前記第 2 の表面形状の重心を算出する重心算出手段が設けられていること を特徴とする請求項 1 3に記載の超音波内視鏡装置。
1 9 . 前記マツチング手段には、 前記第 1の表面形状と前記第 2 の表面形状の慣性主軸を算出する慣性主軸算出手段が設けられて いることを特徴とずる請求項 1 3に記載の超音波内視鏡装置。 2 0 . 前記マッチング手段には、 前記第 1の表面形状と前記第 2 の表面形状の重心を算出する重心算出手段、 も しくは前記第 1の 表面形状と前記第 2の表面形状の慣性主軸を算出する慣性主軸算 出手段を設け、
前記相互相関処理の前に、 前記重心も しくは前記慣性主軸を算 出したことを特徴とする請求項 1 7に記載の超音波内視鏡装置。
PCT/JP2004/008151 2003-06-06 2004-06-04 超音波内視鏡装置 Ceased WO2004107982A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP04736120A EP1632184B1 (en) 2003-06-06 2004-06-04 Ultrasonic endoscope
US11/294,985 US20060183992A1 (en) 2003-06-06 2005-12-06 Ultrasonic endoscope device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2003162845A JP4414682B2 (ja) 2003-06-06 2003-06-06 超音波内視鏡装置
JP2003-162845 2003-06-06

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/294,985 Continuation US20060183992A1 (en) 2003-06-06 2005-12-06 Ultrasonic endoscope device

Publications (1)

Publication Number Publication Date
WO2004107982A1 true WO2004107982A1 (ja) 2004-12-16

Family

ID=33508679

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2004/008151 Ceased WO2004107982A1 (ja) 2003-06-06 2004-06-04 超音波内視鏡装置

Country Status (4)

Country Link
US (1) US20060183992A1 (ja)
EP (1) EP1632184B1 (ja)
JP (1) JP4414682B2 (ja)
WO (1) WO2004107982A1 (ja)

Families Citing this family (72)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9510732B2 (en) 2005-10-25 2016-12-06 Intuitive Surgical Operations, Inc. Methods and apparatus for efficient purging
US10064540B2 (en) 2005-02-02 2018-09-04 Intuitive Surgical Operations, Inc. Visualization apparatus for transseptal access
US8078266B2 (en) 2005-10-25 2011-12-13 Voyage Medical, Inc. Flow reduction hood systems
US7860555B2 (en) 2005-02-02 2010-12-28 Voyage Medical, Inc. Tissue visualization and manipulation system
US11478152B2 (en) 2005-02-02 2022-10-25 Intuitive Surgical Operations, Inc. Electrophysiology mapping and visualization system
US8137333B2 (en) 2005-10-25 2012-03-20 Voyage Medical, Inc. Delivery of biological compounds to ischemic and/or infarcted tissue
US8050746B2 (en) 2005-02-02 2011-11-01 Voyage Medical, Inc. Tissue visualization device and method variations
US20080015569A1 (en) 2005-02-02 2008-01-17 Voyage Medical, Inc. Methods and apparatus for treatment of atrial fibrillation
CA2789262C (en) 2005-04-28 2016-10-04 Proteus Digital Health, Inc. Pharma-informatics system
US9198608B2 (en) 2005-04-28 2015-12-01 Proteus Digital Health, Inc. Communication system incorporated in a container
WO2007028035A2 (en) 2005-09-01 2007-03-08 Proteus Biomedical, Inc. Implantable zero-wire communications system
US8221310B2 (en) 2005-10-25 2012-07-17 Voyage Medical, Inc. Tissue visualization device and method variations
US8956287B2 (en) 2006-05-02 2015-02-17 Proteus Digital Health, Inc. Patient customized therapeutic regimens
US8446410B2 (en) * 2006-05-11 2013-05-21 Anatomage Inc. Apparatus for generating volumetric image and matching color textured external surface
US9055906B2 (en) 2006-06-14 2015-06-16 Intuitive Surgical Operations, Inc. In-vivo visualization systems
US10004388B2 (en) * 2006-09-01 2018-06-26 Intuitive Surgical Operations, Inc. Coronary sinus cannulation
US20080097476A1 (en) 2006-09-01 2008-04-24 Voyage Medical, Inc. Precision control systems for tissue visualization and manipulation assemblies
WO2008028149A2 (en) 2006-09-01 2008-03-06 Voyage Medical, Inc. Electrophysiology mapping and visualization system
US10335131B2 (en) 2006-10-23 2019-07-02 Intuitive Surgical Operations, Inc. Methods for preventing tissue migration
SG175681A1 (en) 2006-10-25 2011-11-28 Proteus Biomedical Inc Controlled activation ingestible identifier
US8718193B2 (en) 2006-11-20 2014-05-06 Proteus Digital Health, Inc. Active signal processing personal health signal receivers
US20080183036A1 (en) 2006-12-18 2008-07-31 Voyage Medical, Inc. Systems and methods for unobstructed visualization and ablation
US8131350B2 (en) 2006-12-21 2012-03-06 Voyage Medical, Inc. Stabilization of visualization catheters
US9226648B2 (en) 2006-12-21 2016-01-05 Intuitive Surgical Operations, Inc. Off-axis visualization systems
CN101686827B (zh) * 2007-01-19 2014-08-13 桑尼布鲁克健康科学中心 具有组合的超声和光学成像装置的成像探头
EP3785599B1 (en) 2007-02-01 2022-08-03 Otsuka Pharmaceutical Co., Ltd. Ingestible event marker systems
WO2008101107A1 (en) 2007-02-14 2008-08-21 Proteus Biomedical, Inc. In-body power source having high surface area electrode
WO2008112577A1 (en) 2007-03-09 2008-09-18 Proteus Biomedical, Inc. In-body device having a multi-directional transmitter
EP2148608A4 (en) 2007-04-27 2010-04-28 Voyage Medical Inc CATHETER FOR VISUALIZATION AND HANDLING OF MANUFACTURED FABRICS OF COMPLEX SHAPES
US8657805B2 (en) 2007-05-08 2014-02-25 Intuitive Surgical Operations, Inc. Complex shape steerable tissue visualization and manipulation catheter
EP3025636B1 (en) 2007-05-11 2017-11-01 Intuitive Surgical Operations, Inc. Visual electrode ablation systems
US8540632B2 (en) 2007-05-24 2013-09-24 Proteus Digital Health, Inc. Low profile antenna for in body device
US8235985B2 (en) 2007-08-31 2012-08-07 Voyage Medical, Inc. Visualization and ablation system variations
US8961412B2 (en) 2007-09-25 2015-02-24 Proteus Digital Health, Inc. In-body device with virtual dipole signal amplification
DK2215726T3 (en) 2007-11-27 2018-04-09 Proteus Digital Health Inc Transbody communication modules with communication channels
JP2009180607A (ja) * 2008-01-30 2009-08-13 Olympus Corp 自動分析装置
US8858609B2 (en) 2008-02-07 2014-10-14 Intuitive Surgical Operations, Inc. Stent delivery under direct visualization
EP3235491B1 (en) 2008-03-05 2020-11-04 Proteus Digital Health, Inc. Multi-mode communication ingestible event markers and systems
US9101735B2 (en) 2008-07-07 2015-08-11 Intuitive Surgical Operations, Inc. Catheter control systems
EP2313002B1 (en) 2008-07-08 2018-08-29 Proteus Digital Health, Inc. Ingestible event marker data framework
US8894643B2 (en) 2008-10-10 2014-11-25 Intuitive Surgical Operations, Inc. Integral electrode placement and connection systems
US8333012B2 (en) 2008-10-10 2012-12-18 Voyage Medical, Inc. Method of forming electrode placement and connection systems
US9468364B2 (en) 2008-11-14 2016-10-18 Intuitive Surgical Operations, Inc. Intravascular catheter with hood and image processing systems
CA2746650A1 (en) 2008-12-11 2010-06-17 Proteus Biomedical, Inc. Evaluation of gastrointestinal function using portable electroviscerography systems and methods of using the same
JP2012514799A (ja) 2009-01-06 2012-06-28 プロテウス バイオメディカル インコーポレイテッド 摂取に関連するバイオフィードバックおよび個別薬物療法の方法およびシステム
TWI517050B (zh) 2009-11-04 2016-01-11 普羅托斯數位健康公司 供應鏈管理之系統
US20110144576A1 (en) * 2009-12-14 2011-06-16 Voyage Medical, Inc. Catheter orientation control system mechanisms
US8694071B2 (en) 2010-02-12 2014-04-08 Intuitive Surgical Operations, Inc. Image stabilization techniques and methods
US9814522B2 (en) 2010-04-06 2017-11-14 Intuitive Surgical Operations, Inc. Apparatus and methods for ablation efficacy
TWI557672B (zh) 2010-05-19 2016-11-11 波提亞斯數位康健公司 用於從製造商跟蹤藥物直到患者之電腦系統及電腦實施之方法、用於確認將藥物給予患者的設備及方法、患者介面裝置
JP5414637B2 (ja) * 2010-08-24 2014-02-12 富士フイルム株式会社 撮像装置および内視鏡システム
WO2015112603A1 (en) 2014-01-21 2015-07-30 Proteus Digital Health, Inc. Masticable ingestible product and communication system therefor
US9756874B2 (en) 2011-07-11 2017-09-12 Proteus Digital Health, Inc. Masticable ingestible product and communication system therefor
KR101898964B1 (ko) 2011-07-21 2018-09-14 프로테우스 디지털 헬스, 인코포레이티드 모바일 통신 장치, 시스템, 및 방법
KR20130015146A (ko) * 2011-08-02 2013-02-13 삼성전자주식회사 의료 영상 처리 방법 및 장치, 영상 유도를 이용한 로봇 수술 시스템
US9235683B2 (en) 2011-11-09 2016-01-12 Proteus Digital Health, Inc. Apparatus, system, and method for managing adherence to a regimen
EP2856098B1 (en) 2012-05-25 2019-10-16 Vascular Imaging Corporation Optical fiber pressure sensor
JP6511439B2 (ja) 2013-06-04 2019-05-15 プロテウス デジタル ヘルス, インコーポレイテッド データ収集および転帰の査定のためのシステム、装置、および方法
JP6236909B2 (ja) * 2013-06-24 2017-11-29 株式会社リコー 浸透過程計測装置及び浸透過程計測方法
CN105324082B (zh) * 2013-10-01 2017-12-12 奥林巴斯株式会社 超声波内窥镜系统及其通信方法
US10327645B2 (en) 2013-10-04 2019-06-25 Vascular Imaging Corporation Imaging techniques using an imaging guidewire
US10084880B2 (en) 2013-11-04 2018-09-25 Proteus Digital Health, Inc. Social media networking based on physiologic information
US10537255B2 (en) 2013-11-21 2020-01-21 Phyzhon Health Inc. Optical fiber pressure sensor
JP2017518147A (ja) 2014-03-28 2017-07-06 インテュイティブ サージカル オペレーションズ, インコーポレイテッド 手術シーンの定量的三次元撮像
KR102397670B1 (ko) 2014-03-28 2022-05-16 인튜어티브 서지컬 오퍼레이션즈 인코포레이티드 정량적 3차원 영상화에 기초한 햅틱 피드백을 갖는 수술 시스템
US10334227B2 (en) 2014-03-28 2019-06-25 Intuitive Surgical Operations, Inc. Quantitative three-dimensional imaging of surgical scenes from multiport perspectives
EP3122232B1 (en) * 2014-03-28 2020-10-21 Intuitive Surgical Operations Inc. Alignment of q3d models with 3d images
CN110251047B (zh) 2014-03-28 2022-01-18 直观外科手术操作公司 手术植入物的定量三维成像和打印
KR102387096B1 (ko) * 2014-03-28 2022-04-15 인튜어티브 서지컬 오퍼레이션즈 인코포레이티드 시계 내의 기기들의 정량적 3차원 시각화
US11564660B2 (en) * 2016-03-04 2023-01-31 Canon Medical Systems Corporation Ultrasonic diagnostic apparatus and method for generating ultrasonic image
CN109843149B (zh) 2016-07-22 2020-07-07 普罗秋斯数字健康公司 可摄入事件标记的电磁感测和检测
JP6832165B2 (ja) * 2017-01-16 2021-02-24 オリンパス株式会社 観察システム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05324835A (ja) * 1992-05-18 1993-12-10 Kobe Steel Ltd 二次元図形のパターンマッチング方法及びその装置
JP2000116655A (ja) 1998-10-14 2000-04-25 Olympus Optical Co Ltd 診断装置
JP2002017729A (ja) * 2000-07-11 2002-01-22 Toshiba Corp 超音波内視鏡診断装置
JP2003038492A (ja) * 2001-07-30 2003-02-12 Pentax Corp 超音波内視鏡装置

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5054491A (en) * 1988-10-17 1991-10-08 Olympus Optical Co., Ltd. Ultrasonic endoscope apparatus
JPH08299260A (ja) * 1995-04-28 1996-11-19 Fuji Photo Optical Co Ltd 超音波内視鏡
JPH10192A (ja) * 1996-04-15 1998-01-06 Olympus Optical Co Ltd 超音波画像診断装置
JPH10262973A (ja) * 1997-03-25 1998-10-06 Olympus Optical Co Ltd 超音波画像診断装置
US6248074B1 (en) * 1997-09-30 2001-06-19 Olympus Optical Co., Ltd. Ultrasonic diagnosis system in which periphery of magnetic sensor included in distal part of ultrasonic endoscope is made of non-conductive material
JP4248615B2 (ja) * 1997-10-23 2009-04-02 オリンパス株式会社 超音波画像診断装置
JP4054104B2 (ja) * 1998-04-10 2008-02-27 オリンパス株式会社 内視鏡画像処理装置
US6625316B1 (en) * 1998-06-01 2003-09-23 Canon Kabushiki Kaisha Image processing apparatus and method, and image processing system
US6306091B1 (en) * 1999-08-06 2001-10-23 Acuson Corporation Diagnostic medical ultrasound systems and methods utilizing estimation of 3-dimensional rigid body transformation
US6443894B1 (en) * 1999-09-29 2002-09-03 Acuson Corporation Medical diagnostic ultrasound system and method for mapping surface data for three dimensional imaging
US7027650B2 (en) * 1999-12-10 2006-04-11 Christian Williame Dynamic computing imagery, especially for visceral osteopathy and for articular kinetics
DE10015826A1 (de) * 2000-03-30 2001-10-11 Siemens Ag System und Verfahren zur Erzeugung eines Bildes

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05324835A (ja) * 1992-05-18 1993-12-10 Kobe Steel Ltd 二次元図形のパターンマッチング方法及びその装置
JP2000116655A (ja) 1998-10-14 2000-04-25 Olympus Optical Co Ltd 診断装置
JP2002017729A (ja) * 2000-07-11 2002-01-22 Toshiba Corp 超音波内視鏡診断装置
JP2003038492A (ja) * 2001-07-30 2003-02-12 Pentax Corp 超音波内視鏡装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP1632184A4 *

Also Published As

Publication number Publication date
EP1632184B1 (en) 2012-04-04
US20060183992A1 (en) 2006-08-17
JP4414682B2 (ja) 2010-02-10
JP2004358096A (ja) 2004-12-24
EP1632184A4 (en) 2007-11-14
EP1632184A1 (en) 2006-03-08

Similar Documents

Publication Publication Date Title
WO2004107982A1 (ja) 超音波内視鏡装置
CN100518656C (zh) 超声波摄像装置及其摄像数据处理方法以及处理程序
US6416476B1 (en) Three-dimensional ultrasonic diagnosis apparatus
JP5208495B2 (ja) 医療用システム
JP4551051B2 (ja) 超音波診断装置
JP4868959B2 (ja) 体腔内プローブ装置
WO2007114375A1 (ja) 超音波診断装置及び超音波診断装置制御方法
WO2004028375A1 (ja) 超音波診断装置
CN102300505A (zh) 超声波诊断装置以及图像数据显示用控制程序
JP2009089736A (ja) 超音波診断装置
KR101117026B1 (ko) 서로 다른 영상간에 영상 정합을 수행하는 영상 정합 시스템 및 방법
JP4869197B2 (ja) 医用ガイド装置
JP2009061086A (ja) 超音波診断装置、並びに、画像処理方法及びプログラム
JP2002017729A (ja) 超音波内視鏡診断装置
JPH10262963A (ja) 超音波画像診断装置
JP2006055493A (ja) 超音波診断装置および医用画像解析装置
JP2010088699A (ja) 医療画像処理システム
JP2004121488A (ja) 超音波診断装置
JP4869189B2 (ja) 医用ガイドシステム
JP2000116655A (ja) 診断装置
JP4651379B2 (ja) 超音波画像処理装置及び超音波画像処理方法、並びに、超音波画像処理プログラム
JP2004209277A (ja) 超音波画像診断装置
JP2013013452A (ja) 超音波診断装置及び制御プログラム
JPH11113912A (ja) 超音波画像診断装置
JP4847126B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2004736120

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 11294985

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 2004736120

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 11294985

Country of ref document: US