WO2007102296A1 - 医療用画像処理装置及び医療用画像処理方法 - Google Patents
医療用画像処理装置及び医療用画像処理方法 Download PDFInfo
- Publication number
- WO2007102296A1 WO2007102296A1 PCT/JP2007/052347 JP2007052347W WO2007102296A1 WO 2007102296 A1 WO2007102296 A1 WO 2007102296A1 JP 2007052347 W JP2007052347 W JP 2007052347W WO 2007102296 A1 WO2007102296 A1 WO 2007102296A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- edge
- information
- shape
- dimensional model
- image processing
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/04—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor combined with photographic or television appliances
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10068—Endoscopic image
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30028—Colon; Small intestine
- G06T2207/30032—Colon polyp
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Definitions
- the present invention relates to a medical image processing apparatus and a medical image processing method, and more particularly to a medical image processing apparatus and a medical image processing method capable of detecting a raised shape in a three-dimensional model of a living tissue in a body cavity. Is.
- An endoscope system configured to include an endoscope, a medical image processing apparatus, and the like is widely used in the medical field and the like.
- the endoscope system is formed by, for example, an insertion portion that is inserted into a body cavity as a living body, an objective optical system that is disposed at a distal end portion of the insertion portion, and an image formed by the objective optical system.
- An endoscope including an imaging unit that captures an image of a body cavity and outputs it as an imaging signal, and displays the image of the body cavity on a monitor as a display unit based on the imaging signal
- a medical image processing apparatus for performing processing for the purpose. Then, based on the image in the body cavity displayed on the monitor or the like as the display unit, the user observes the organ or the like as the subject in the body cavity, for example.
- an endoscope system having the above-described configuration can also take an image of a digestive tract mucous membrane such as a large intestine as a subject in a body cavity. Therefore, the user can comprehensively observe various findings such as the color of the mucous membrane, the shape of the lesion, and the fine structure of the mucosal surface.
- a two-dimensional image including an image of a biological tissue such as a fold of the large intestine or a polyp as an image of a biological tissue having a raised shape imaged by an endoscope is a field of view of the endoscope.
- the invisible region in the above-described two-dimensional image is generally called occlusion and is a region where it is difficult to estimate an accurate three-dimensional model. Therefore, in the part where the occlusion of the two-dimensional image occurs, for example, there is no estimation result at the corresponding position of the three-dimensional model based on the two-dimensional image, or the three-dimensional image based on the two-dimensional image. An unreliable estimation result may be calculated at the corresponding position in the model.
- US Patent Application No. 20030223627 A method for detecting a lesion such as a polyp is obtained by observation using, for example, CT (Computed Tomog rapy). This method is assumed to be applied to a 3D model that does not generate occlusion. Therefore, the method for detecting a lesion site such as a polyp proposed in US Patent Application No. 20030223627 is, for example, a two-dimensional image having an occlusion portion obtained by observation using an endoscope. When applied to a 3D model estimated on the basis of this, the detection accuracy of a lesion site such as a polyp is lowered. As a result, when an object is observed while using an endoscope system that adopts the method power proposed in US Patent Application No. 20030223627, the burden on the user to find a diseased site such as a polyp There is a problem that becomes large.
- CT Computer Tomog rapy
- the present invention has been made in view of the above-described points, and in the case where observation is performed using a three-dimensional model estimated based on a two-dimensional image having occlusion, a polyp by a user, etc.
- the burden on the user by preventing oversight of the lesion site It is an object of the present invention to provide a medical image processing apparatus and a medical image processing method capable of reducing the above-mentioned problem.
- the first medical image processing apparatus extracts edges of the two-dimensional image based on the two-dimensional image of the biological tissue image in the body cavity input from the medical imaging apparatus.
- a three-dimensional model estimation unit that estimates a three-dimensional model of the living tissue based on the two-dimensional image, and a position of each botacel where the edge estimated as part of the three-dimensional model exists.
- the shape feature amount of at least a part of the botacell extraction unit that extracts a predetermined botacell group that is a target of calculation of the shape feature amount and the predetermined battercell group that is extracted by the botacell extraction unit
- a shape feature amount calculating unit for calculating a shape, and extracting, from the at least some of the at least a part of the botacells, a group of botasels in which a three-dimensional model is estimated as a predetermined shape.
- the voxel group characterized by having a a raised shape detecting section for detecting as a Botaseru group constituting the raised shape in the three-dimensional model of the living tissue.
- the second medical image processing apparatus is the first medical image processing apparatus, wherein the predetermined botacell group includes a predetermined one of the botacells included in each of the botacells where the edge exists. Is a group of botasels existing on the near side of the visual field direction of the medical imaging apparatus.
- the predetermined one of the button cells is the one in the field-of-view direction among the button cells in which the edge exists. It is characterized in that it is one bottom cell that exists at the innermost side.
- a fourth medical image processing apparatus is characterized in that, in the first to third medical image processing apparatuses, the predetermined shape is a convex shape.
- a fifth medical image processing apparatus is based on a two-dimensional image of an image of a living tissue in a body cavity that is inputted, and extracts an edge of the two-dimensional image by fine line display. For each edge extracted by the edge extraction unit, position information about the position of the end of one edge, and information about the shape of the one edge calculated based on the position information
- An edge information acquisition unit that acquires shape information as a report, a three-dimensional model estimation unit that estimates a three-dimensional model of the biological tissue based on the two-dimensional image, and the position information and the shape information.
- a depth information that is information related to the depth of the three-dimensional model of the biological tissue and a height information that is information related to the height of the three-dimensional model of the biological tissue; and the depth And a raised shape detecting unit for detecting whether or not the biological tissue has a raised shape based on the information and the height information.
- the shape information includes information on a position of a midpoint in the thin line of the one edge, and the one Information on the position of the midpoint of a straight line connecting two end points at the end of the edge.
- a first medical image processing method extracts an edge of a two-dimensional image based on a two-dimensional image of a biological tissue image in a body cavity input from a medical imaging device.
- a three-dimensional model estimation step for estimating a three-dimensional model of the living tissue based on the two-dimensional image, and for each of the botasels where the edges estimated as part of the three-dimensional model exist. Based on the position, the step of extracting a predetermined botacell group for calculation of the shape feature quantity, and the shape feature quantity of at least a part of the boutelle cells extracted from the boutelle cell extraction step.
- a three-dimensional model is estimated as a predetermined shape among the at least some of the botasels. Extracts the Taseru group, the voxel group, characterized by having a a raised shape detecting step of detecting as Botaseru group forming the raised shape in the three-dimensional model of the living tissue.
- a second medical image processing method is the first medical image processing method according to the first medical image processing method, in which the predetermined botacell group is included in each of the botacells where the edge exists. Is a group of botasels existing on the near side of the visual field direction of the medical imaging apparatus.
- the third medical image processing method in the present invention is the second medical image processing method.
- a fourth medical image processing method is characterized in that, in the first to third medical image processing methods, the predetermined shape is a convex shape.
- a fifth medical image processing method is based on an input two-dimensional image of a body tissue image in a body cavity, and extracts an edge of the two-dimensional image by thin lines. For each edge extracted in the step and the edge extraction step, the position information of the position that becomes the end of the one edge, and the shape that is information on the shape of the one edge calculated based on the position information
- An edge information acquisition step for acquiring information
- a three-dimensional model estimation step for estimating a three-dimensional model of the living tissue based on the two-dimensional image
- the position information and the shape information based on the position information and the shape information.
- a sixth medical image processing method is the fifth medical image processing method, wherein the shape information includes information on a position of a midpoint in the thin line of the one edge, and the one Information on the position of the midpoint of a straight line connecting two end points at the end of the edge.
- FIG. 1 is a diagram showing an example of the overall configuration of an endoscope system in which a medical image processing apparatus according to each embodiment of the present invention is used.
- FIG. 2 is a schematic diagram showing a state where the endoscope of FIG. 1 is inserted into a tubular organ.
- FIG. 3 is a schematic diagram showing images of tubular organs and living tissue imaged by an endoscope in the state of FIG. 2.
- FIG. 4 is a flowchart showing a procedure of processing performed by the medical image processing apparatus of FIG. 1 in the first embodiment when an edge portion in a two-dimensional image of a biological tissue image is extracted.
- FIG. 5 is a flowchart showing a procedure of processing performed by the medical image processing apparatus of FIG. 1 in the first embodiment when detecting a raised shape in a three-dimensional model of biological tissue.
- FIG. 6 is a diagram showing an example of a three-dimensional model of biological tissue estimated by the medical image processing apparatus of FIG.
- FIG. 7 is a diagram showing an example of a region to be processed in the flowchart of FIG. 5 for detecting a raised shape in a three-dimensional model of biological tissue.
- FIG. 8 is a diagram showing an example different from FIG. 7 of a region to be processed in the flowchart of FIG. 5 for detecting a raised shape in a three-dimensional model of biological tissue.
- FIG. 9 is a flowchart showing an outline of a processing procedure performed by the medical image processing apparatus of FIG. 1 in the second embodiment.
- FIG. 10 is a flowchart showing a specific procedure of attention point extraction processing in FIG.
- FIG. 11 is a flowchart showing a specific procedure of the raised shape detection process in FIG.
- FIG. 12 is a diagram showing an example of a case where a two-dimensional image including an edge portion of a biological tissue is thinned by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 13 is a diagram showing an example of the positional relationship of each point of interest in the two-dimensional image shown in FIG.
- FIG. 14 is a diagram showing an example of the positional relationship of each point of interest in a three-dimensional model of biological tissue obtained by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 15 shows an example different from FIG. 14 of the positional relationship of each point of interest in the three-dimensional model of biological tissue obtained by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 1 is a diagram illustrating an example of the overall configuration of an endoscope system in which the medical image processing apparatus according to the present embodiment is used.
- FIG. 2 is a schematic diagram showing a state when the endoscope of FIG. 1 is inserted into a tubular organ.
- FIG. 3 shows a tubular organ and a living body image captured by an endoscope in the state of FIG. It is a schematic diagram which shows the image of a weaving.
- FIG. 4 is a flowchart showing a procedure of processing performed by the medical image processing apparatus of FIG. 1 when extracting an edge portion in a two-dimensional image of a biological tissue image.
- FIG. 1 is a diagram illustrating an example of the overall configuration of an endoscope system in which the medical image processing apparatus according to the present embodiment is used.
- FIG. 2 is a schematic diagram showing a state when the endoscope of FIG. 1 is inserted into a tubular organ.
- FIG. 3 shows a tubular organ and a living body image captured by an endo
- FIG. 5 is a flowchart showing a procedure of processing performed by the medical image processing apparatus of FIG. 1 when detecting a raised shape in a three-dimensional model of biological tissue.
- FIG. 6 is a diagram illustrating an example of a three-dimensional model of biological tissue estimated by the medical image processing apparatus of FIG.
- FIG. 7 is a diagram showing an example of a region to be processed in the flowchart of FIG. 5 for detecting a raised shape in a three-dimensional model of a living body and weaving.
- FIG. 8 is a diagram showing an example different from FIG. 7 of the region to be processed in the flow chart of FIG. 5 for detecting the raised shape in the three-dimensional model of living body and weaving.
- the endoscope system 1 includes a medical observation apparatus 2 that captures an image of a subject and outputs a two-dimensional image of the image of the subject, a personal computer, and the like.
- a medical image processing device 3 that performs image processing on the video signal of the two-dimensional image output from the observation device 2 and outputs the video signal after the image processing as an image signal, and medical image processing And a monitor 4 for displaying an image based on the image signal output from the device 3.
- the medical observation device 2 is inserted into a body cavity, images the subject in the body cavity and outputs it as an imaging signal, and the subject imaged by the endoscope 6
- the signal processing is performed on the imaging signal output from the endoscope 6 to generate a two-dimensional image.
- a camera control unit (hereinafter abbreviated as CCU) 8 that outputs as a video signal, and a monitor that displays the image of the subject imaged by the endoscope 6 based on the video signal of the two-dimensional image that also outputs CCU 8 power 9
- the principal part is comprised.
- the endoscope 6 as a medical imaging device includes an insertion portion 11 that is inserted into a body cavity and an operation portion 12 that is provided on the proximal end side of the insertion portion 11.
- a light guide 13 for transmitting illumination light supplied from the light source device 7 is inserted into a portion from the proximal end side in the insertion portion 11 to the distal end portion 14 on the distal end side in the insertion portion 11. Yes.
- the light guide 13 has a distal end side disposed at the distal end portion 14 of the endoscope 6 and a rear end side connected to the light source device 7. Since the light guide 13 has such a configuration, The illumination light supplied from the device 7 is transmitted by the light guide 13 and then the illumination window force (not shown) provided on the distal end surface of the distal end portion 14 of the insertion portion 11 is also emitted. Then, illumination light is emitted from an illumination window (not shown) to illuminate a living tissue or the like as a subject.
- an objective optical system 15 attached to an observation window (not shown) adjacent to an illumination window (not shown) and an imaging position of the objective optical system 15 are arranged.
- An imaging unit 17 having an imaging device 16 constituted by a CCD (charge coupled device) or the like is provided. With such a configuration, the subject image formed by the objective optical system 15 is captured by the imaging element 16 and then output as an imaging signal.
- the image sensor 16 is connected to the CCU 8 via a signal line.
- the image sensor 16 is driven based on the drive signal output from the CCU 8, and outputs an image signal to the CCU 8.
- the imaging signal input to the CCU 8 is subjected to signal processing in a signal processing circuit (not shown) provided in the CCU 8, thereby being converted and output as a video signal of a two-dimensional image.
- the video signal of the two-dimensional image output from the CCU 8 is output to the monitor 9 and the medical image processing device 3.
- an image of the subject based on the video signal output from the CCU 8 is displayed on the monitor 9 as a two-dimensional image.
- the medical image processing apparatus 3 performs an AZD conversion on the video signal of the two-dimensional image output from the medical observation apparatus 2, and outputs the image input section 21 from the image input section 21.
- CPU 22 as a central processing unit that performs image processing on a video signal to be processed, a processing program storage unit 23 in which a processing program related to the image processing is written, a video signal output from the image input unit 21, and the like
- an information storage unit 25 for storing image data as an image processing result of the CPU 22 and the like.
- the medical image processing apparatus 3 also includes a storage device interface 26, a hard disk 27 as a storage device that stores image data and the like as image processing results of the CPU 22 via the storage device interface 26, Based on the image data as the image processing result of the CPU 22, a display processing unit 28 performs a display process for displaying the image data on the monitor 4 and outputs the image data after the display process as an image signal.
- C It has an input operation unit 29 constituted by a keyboard or the like that allows a user to input parameters for image processing performed by the PU 22 and operation instructions for the medical image processing apparatus 3.
- the monitor 4 displays an image based on the image signal output from the display processing unit 28.
- the image input unit 21, the CPU 22, the processing program storage unit 23, the image storage unit 24, the information storage unit 25, the storage device interface 26, the display processing unit 28, and the input operation unit 29 of the medical image processing apparatus 3 are used. Are connected to each other via a data bus 30.
- the user inserts the insertion portion 11 of the endoscope 6 into a tubular organ 31 such as the large intestine, for example.
- a tubular organ 31 such as the large intestine
- an image of the living tissue 31A which is a lesion site existing on the inner wall of the tubular organ 31
- the imaging unit 17 provided at the distal end portion 14. Images are taken as shown in Fig. 3.
- the images of the tubular organ 31 and the living tissue 31A imaged by the imaging unit 17 as an image as shown in FIG. 3 are output to the CCU 8 as an imaging signal.
- the CCU 8 performs signal processing on the imaging signal output from the imaging device 16 of the imaging unit 17 in a signal processing circuit (not shown), thereby converting the imaging signal as a video signal of a two-dimensional image. Output.
- the monitor 9 displays the images of the tubular organ 31 and the living tissue 31A as a two-dimensional image as shown in FIG. 3, for example, based on the video signal output from the CCU 8. Further, the CCU 8 outputs a video signal of a two-dimensional image obtained by performing signal processing on the imaging signal output from the imaging device 16 of the imaging unit 17 to the medical image processing apparatus 3.
- the video signal of the two-dimensional image output to the medical image processing device 3 is AZD converted by the image input unit 21 and then input to the CPU 22.
- the CPU 22 performs, for example, the following processing based on the video signal of the two-dimensional image output from the image input unit 21 and the processing program written in the processing program storage unit 23. Extracting the edge portion of the living tissue 31A in the two-dimensional image
- the CPU 22 having a function as an edge extraction unit is output from the image input unit 21.
- the CPU 22 having a function as an edge extraction unit is output from the image input unit 21.
- the CPU 22 having a function as an edge extraction unit is output from the image input unit 21.
- the video signal of the two-dimensional image for example, by applying a bandpass filter to the red component of the two-dimensional image, all edge portions included in the two-dimensional image are extracted (FIG. 4). Step Sl).
- the CPU 22 having a function as an edge extracting unit fine-line all the extracted edge portions (step S2 in FIG. 4), and at the same time, out of all the edge portions that have been thin-lined,
- the length L of the edge part E is calculated (step S3 in FIG. 4). Further, the CPU 22 determines whether or not the length L of one edge portion E is longer than the threshold value thLl and shorter than the threshold value thL2.
- step S4 in FIG. 4 when the CPU 22 detects that the length L force of one edge portion E is a length equal to or shorter than a predetermined threshold thLl, or is equal to or larger than the threshold thL2 (step S4 in FIG. 4), The process shown in step S11 of FIG. 4 described later is performed.
- the CPU 22 obtains the normal NCc subtracted from the midpoint Cc of one edge portion E, and obtains N normal NCn subtracted from each control point Cn (step in FIG. 4). S6). Thereafter, the CPU 22 detects the number Na of the N normals NCn that intersect the normal NCc (step S7 in FIG. 4).
- the CPU 22 determines whether or not the number Na of N normals NCn that intersects the normal NCc is greater than the threshold value tha. If the CPU 22 detects that the number Na intersects the normal NCc is greater than the threshold value tha (step S8 in FIG. 4), the pixel group ip included in one edge portion E is converted to the edge of the living tissue 31A. It is determined that the pixel group is included in the part, and the value of the variable edge (i) in each pixel of the pixel group ip is set to ON (step S9 in FIG. 4). Furthermore, when the CPU 22 detects that the number Na of crossing the normal NCc is equal to or less than the threshold value tha (step S8 in FIG.
- the pixel group ip included in one edge portion E is converted to the living tissue 31. It is determined that the pixel group is not included in the edge portion of A, and the value of the variable edge (i) in each pixel of the pixel group ip is set to OFF (step S10 in FIG. 4).
- the CPU 22 determines whether or not the processing has been completed for all the extracted edge portions. When the CPU 22 detects that the processing has not been completed for all the extracted edge portions (step S11 in FIG. 4), the CPU 22 performs the steps in FIG. Steps S3 to S10 in Fig. 4 are performed. In addition, when the CPU 22 detects that the processing has been completed for all the extracted edge portions (step S11 in FIG. 4), a series of steps for extracting the edge portions of the living tissue 31A in the two-dimensional image are performed. The process ends.
- the CPU 22 obtains each edge of the two-dimensional image as a processing result obtained by performing the above-described series of processing for extracting the edge portion of the living tissue 31A in the two-dimensional image.
- the value of the variable edge (i) in the pixel group ip included in the section is temporarily stored in the information storage section 25.
- the CPU 22 having a function as a three-dimensional model estimation unit converts, for example, luminance information of the video signal of the two-dimensional image output from the image input unit 21 by a process such as geometric conversion. Based on this, image processing for acquiring image data necessary for estimating the three-dimensional model of the tubular organ 31 and the living tissue 31A is performed. In other words, the CPU 22 generates, for example, a voxel corresponding to each pixel in the two-dimensional image by processing such as geometric transformation, and acquires the voxel as image data for estimating the three-dimensional model. To do. In other words, the pixel group ip is converted as the botacell group ib by the above-described processing.
- the CPU 22 uses a plane including the botacell group ib whose variable edge (i) is ON as image data for estimating the three-dimensional model of the living tissue 31A shown in FIG. Obtain 3 lb of data for a tissue boundary.
- the living tissue 31A is estimated as a three-dimensional model having a shape as shown in FIG. 6, for example, when the z-axis direction is the visual field direction during observation with the endoscope 6.
- the CPU 22 detects the raised shape in the three-dimensional model of the living tissue 31A by performing the processing described below based on the data of the living tissue boundary 3 lb.
- the CPU 22 selects an endoscope from among the botacell group ib whose variable edge (i) is ON.
- the z-coordinate is the most One large botacell is extracted, and the z coordinate of the one botacell is set as Maxz (step S 21 in FIG. 5).
- the CPU 22 determines that the endoscope 6 has a position more than the position of the one of all the botasels obtained as image data for estimating the three-dimensional model of the tubular organ 31 and the living tissue 31A.
- a botacell group rb whose z coordinate is smaller than Maxz is detected as a botacell existing in front of the viewing direction (step S22 in FIG. 5).
- botacell group rb is a botacell group included in the raised shape detection target area AR shown in FIG.
- Botacell group has R number of Botacells.
- the CPU 22 having the function as the botacell extraction unit is configured to determine the position of each botacell included in the botacell group ib whose variable edge (i) is ON. Based on the above, the botacell group rb having R voxels contained in the raised shape detection target area AR is extracted as a predetermined botacell group that is a target for calculating the shape feature quantity.
- the CPU 22 having a function as a shape feature quantity calculation unit sets the variable a to 1.
- Shapelndex value and Curvedness value can be calculated by using a method similar to the method described in US Patent Application No. 20030223627, for example. For this reason, in the present embodiment, the description of the calculation method of the Shapelndex value and the Curvedness value in one botasel Ba is omitted.
- the CPU 22 compares the Shapelndex value SBa with a preset Shapelndex value threshold value Sth, and compares the Curvedness value CBa with a preset Curvedness value threshold value Cth.
- the CPU 22 having a function as a raised shape detecting unit performs the above-described processing, and the 3D model is projected as a process for detecting whether or not the living tissue 31A has a raised shape. Extract the estimated Botacel group Perform the process.
- the threshold value Sth is set to 0.9 and the threshold value Cth is set to 0 in order to detect the living tissue 31A having a convex shape as the raised shape. Assume that it is set as 2.
- One bocellel Ba is determined as a bocellel constituting a part of the raised shape, and the value of the variable ryuuki (Ba) in the one bocellel Ba is set to ON (step S27 in FIG. 5).
- the CPU 22 If the CPU 22 detects that the Shapelndex value SBa is less than or equal to the threshold value Sth or that the Curve dness value CBa is less than or equal to the threshold value Cth (step S26 in FIG. 5), the CPU 22 raises one botel cell Ba. It is determined that the cell is not a part of the shape, and the value of the variable ryuuki (Ba) in the one voxel Ba is turned OFF (step S28 in FIG. 5).
- the CPU 22 applies a character string or coloring for indicating that the living tissue 31A has a raised shape to a three-dimensional model having a button cell whose ryuuki (Ba) value is ON.
- Control for superimposing the three-dimensional model on the display processing unit 28 is performed.
- the monitor 4 displays a three-dimensional model of the tubular organ 31 and the living tissue 31A so that the user can easily find the living tissue 31A having the raised shape.
- the CPU 22 detects each pixel on the two-dimensional image that exists at a position corresponding to the position of each of the botasels based on the position of each of the botasels whose ryuuki (Ba) value is ON.
- the display processing unit 28 is controlled to superimpose a character string or coloring to indicate that the tissue 31A has a raised shape on the two-dimensional image having the pixels. It may be performed.
- the raised shape detection target area AR is not limited to an area that includes a botacell group whose z coordinate is smaller than Maxz, for example, as shown in FIG. Thus, it may be determined as a region including only the botacell group existing in the vicinity of 3 lb of the biological tissue boundary.
- the medical image processing apparatus 3 of the present embodiment can prevent the user from overlooking a lesion site such as a polyp by performing the series of processes described above, thereby reducing the burden on the user. Can be reduced.
- the medical image processing apparatus 3 of the present embodiment performs the above-described processing for detecting the raised shape of the living tissue 31A in the botell cell included in the raised shape detection target area AR in the three-dimensional model. Only for groups. Therefore, the medical image processing apparatus 3 of the present embodiment can detect the raised shape of the living tissue 31A at high speed.
- FIG. 9 is a flowchart showing an outline of a processing procedure performed by the medical image processing apparatus of FIG. 1 in the second embodiment.
- FIG. 10 is a flowchart showing a specific procedure of the attention point extraction process in FIG.
- FIG. 11 is a flowchart showing a specific procedure of the raised shape detection process in FIG.
- FIG. 12 is a diagram illustrating an example of a case where a two-dimensional image including an edge portion of a living tissue is thinned by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 13 is a diagram showing an example of the positional relationship between the points of interest in the two-dimensional image shown in FIG. FIG.
- FIG. 14 is a diagram showing an example of the positional relationship of each point of interest in the three-dimensional model of biological tissue obtained by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 15 is an example different from FIG. 14 of the positional relationship of each point of interest in the three-dimensional model of biological tissue obtained by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- FIG. 14 is a diagram showing an example of the positional relationship of each point of interest in the three-dimensional model of biological tissue obtained by the medical image processing apparatus of FIG. 1 performing the processing of the second embodiment.
- the CPU 22 having a function as an edge extraction unit is based on a video signal of a two-dimensional image obtained by imaging the tubular organ 31 and the living tissue 31A output from the image input unit 21, for example, the two-dimensional image
- a bandpass filter to the red component of the image
- the CPU 22 having a function as an edge extraction unit binarizes and thins all the extracted edge portions (step S32 in FIG. 9), and generates a thinned two-dimensional image. Attention point extraction processing is performed (step S33 in FIG. 9).
- the CPU 22 having a function as an edge information acquisition unit acquires the end points K1 and K2 as the position information of the position of the end of one edge in the thin line of the one edge portion E1 (FIG. 10). Step S331).
- the CPU 22 having a function as an edge information acquisition unit calculates a line segment length LE1 of the thin line of the one edge portion E1 as shape information regarding the shape of the one edge portion E1 (FIG. 10). Along with step S332), acquire the midpoint KE that becomes LEI X 1Z2 (step S333 in Fig. 10) o
- the CPU 22 having a function as an edge information acquisition unit extracts a straight line F1 connecting the end point K1 and the end point K2 as shape information regarding the shape of the one edge part E1, and calculates the length LF1 of the straight line. Along with the calculation (step S334 in FIG. 10), the midpoint KF that becomes LF1 X 1Z2 is obtained (step S335 in FIG. 10).
- the CPU 22 extracts the end point Kl, the end point ⁇ 2, the midpoint KE, and the midpoint KF as attention points in one edge portion E1 (step S336 in FIG. 10).
- the end point Kl, the end point ⁇ 2, the midpoint KE, and the midpoint KF in one edge portion E1 are extracted as attention points each having a positional relationship as shown in FIG. 13, for example.
- the CPU 22 determines whether or not the processing has been completed for all the extracted edge portions.
- step S337 in FIG. 10 If the CPU 22 detects that the processing has not been completed for all the extracted edge portions (step S337 in FIG. 10), the CPU 22 performs the steps in FIG. The processing from S331 to step S336 in FIG. 10 is performed. CPU22 also extracts When it is detected that the processing for all the edge portions has been completed (step S337 in FIG. 10), the attention point extraction processing in the two-dimensional image is terminated.
- the CPU 22 obtains the end point K1 and the end point that are the attention points of each edge portion of the two-dimensional image.
- Data of K2, midpoint 2 and midpoint KF are temporarily stored in the information storage unit 25
- the CPU 22 having a function as a three-dimensional model estimation unit after performing the attention point extraction process in the two-dimensional image, is output from the image input unit 21 by a process such as geometric transformation, for example. Based on the luminance information of the video signal of the two-dimensional image, image processing is performed to obtain image data necessary for estimating the three-dimensional model of the tubular organ 31 and the living tissue 31 ((step S34 in FIG. 9). ). In other words, the CPU 22 generates, for example, a votacell corresponding to each pixel in the two-dimensional image by processing such as geometric conversion, and acquires the voxel as image data for estimating the three-dimensional model.
- the CPU 22 uses one of the living tissue 31A based on the data of the end point Kl, the end point ⁇ 2, the midpoint ⁇ , and the midpoint KF, which are data of each attention point stored in the information storage unit 25.
- Estimate the coordinates of each point of interest in the 3D model of the edge E1 (Step S35 in Fig. 9).
- the endpoint Kl, the endpoint ⁇ 2, the midpoint ⁇ and the midpoint in the 3D model of one edge E1 KF is estimated by the above-described processing, for example, as having a positional relationship as shown in FIG.
- the CPU 22 performs a ridge shape detection process as described below, whereby a ridge shape in the three-dimensional model of the living tissue 31A. Is detected (step S36 in FIG. 9).
- the CPU 22 having a function as a three-dimensional model information acquisition unit, as depth information in the three-dimensional model of the living tissue 31 ⁇ , the three-dimensional model of one edge portion E1.
- the z coordinate KFz of the midpoint KF, the z coordinate Klz of the end point K1, and the z coordinate K2z of the end point K2 are extracted (step S361 in FIG. 11).
- the CPU 22 having a function as a raised shape detecting unit determines whether the z coordinate KFz of the midpoint KF The value is compared with the z-coordinate Klz value of the end point Kl and the z-coordinate K2z value of the end point K2. Then, the CPU 22 determines that the value of the ⁇ coordinate KFz of the midpoint KF is smaller than the value of the ⁇ coordinate Klz of the end point K1, and the value of the z coordinate KFz of the midpoint KF is smaller than the value of the z coordinate K2z of the end point K2. If small is detected (step S362 in FIG. 11), the processing shown in step S363 in FIG. 11 described later is further performed.
- the CPU 22 determines that the z-coordinate KFz value of the midpoint KF is greater than or equal to the z-coordinate Klz value of the end point K1, or the z-coordinate KFz value of the midpoint KF is greater than or equal to the z-coordinate K2z value of the end point K2. If it is detected (step S362 in FIG. 11), it is determined that one edge portion ⁇ 1 of the living tissue 31 ⁇ ⁇ is not an edge portion constituting the raised shape, and step S367 in FIG. Perform the process shown.
- the CPU 22 having a function as a three-dimensional model information acquisition unit determines that the z-coordinate KFz value of the midpoint KF is smaller than the z-coordinate Klz value of the end point Kl and the z-coordinate KFz value of the midpoint KF. Is detected to be smaller than the value of the z-coordinate K2z of the end point K2 (step S362 in Fig. 11), the straight line F2 connecting the end point K1 and the end point ⁇ 2 in the 3D model is extracted, and A vertical line F3 that is a straight line connecting the dot ⁇ and the straight line F2 is acquired (step S363 in FIG. 11).
- the CPU 22 acquires the line segment length LF3 of the vertical line F3 as the height information in the three-dimensional model of the living body yarn and weave 31A (step S364 in FIG. 11), and the line segment length LF3 and the line segment Compare with the threshold value thLF3 of minute length LF3.
- step S365 in FIG. 11 when the CPU 22 having the function as the raised shape detecting unit detects that the value of the line segment length LF3 is larger than the threshold value thLF3 (step S365 in FIG. 11), one of the living tissue 31A is detected.
- the edge part E1 is determined to be an edge part constituting the raised shape, and the value of the variable ryuuki (E1) in all the botasel groups of the one edge part E1 is set to ON (step in FIG. 11). S366).
- the edge E1 of the living tissue 31A is not an edge constituting the raised shape.
- the value of the variable ryuuki (El) in all voxel groups possessed by the one edge portion E1 is set to OFF (step S367 in FIG. 11).
- the CPU 22 determines whether or not the processing has been completed for all the extracted edge portions.
- step S368 in FIG. 11 When the CPU 22 detects that the processing has not been completed for all the extracted edge portions (step S368 in FIG. 11), the CPU 22 in FIG. The processing from step S361 to step S367 in FIG. 11 is performed. Further, when the CPU 22 detects that the processing for all the extracted edge portions is completed (step S368 in FIG. 11), the CPU 22 ends the raised shape detection processing in the three-dimensional model of the living tissue 31A.
- the CPU 22 completes the series of image processing shown in FIG. 9 when the above-described raised shape detection processing shown in FIG. 11 is completed.
- the CPU 22 is not limited to the comparison and determination based on the value of the z coordinate of each point of interest.
- the comparison and determination based on the red component value of each attention point as described in (1) may be performed.
- the CPU 22 extracts the red component value RKF at the midpoint KF, the red component value RK1 at the end point K1, and the red component value RK2 at the end point K2. Then, the CPU 22 detects that the red component value RKF at the midpoint KF is greater than the red component value RK1 at the end point K1, and that the red component value RKF at the midpoint KF is greater than the red component value RK2 at the end point K2. In this case, the processing in step S363 in FIG. 11 is further performed.
- the CPU 22 confirms that the red component value RKF at the midpoint KF is less than or equal to the red component value RK1 at the end point K1, or the red component value RKF at the midpoint KF is less than or equal to the red component value RK2 at the end point K2. If it is detected, it is determined that one edge portion E1 of the living tissue 31A is not an edge portion constituting the raised shape, and the processing of step S367 in FIG. 11 described above is performed.
- the CPU 22 applies a character string or coloring for indicating that the living tissue 31A has a raised shape to a three-dimensional model having a button cell whose ryuuki (El) value is ON.
- Control for superimposing the three-dimensional model on the display processing unit 28 is performed.
- the monitor 4 displays a three-dimensional model of the tubular organ 31 and the living tissue 31A so that the user can easily find the living tissue 31A having the raised shape.
- the CPU 22 detects each pixel on the two-dimensional image existing at a position corresponding to the position of each of the botasels based on the position of each of the botasels whose ryuuki (E1) value is ON.
- the display processing unit 28 is controlled to superimpose a character string or coloring to indicate that the tissue 31A has a raised shape on the two-dimensional image having the pixels. It may be performed.
- the medical image processing apparatus 3 can prevent the user from overlooking a lesion site such as a polyp by performing the series of processes described above, thereby reducing the burden on the user. Can be reduced.
- the medical image processing apparatus 3 quantitatively detects whether or not the above-described living tissue 31A has a raised shape while calculating a value specific to the raised shape. Therefore, the medical image processing apparatus 3 of the present embodiment can accurately detect the raised shape of the living tissue 31A.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Endoscopes (AREA)
Abstract
本発明の医療用画像処理装置は、入力される生体組織の像の二次元画像のエッジを抽出するエッジ抽出部と、前記二次元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデル推定部と、前記三次元モデルの一部として推定された前記エッジが存在する各ボクセルの位置に基づき、形状特徴量の算出対象となる所定のボクセル群を抽出するボクセル抽出部と、前記所定のボクセル群のうち、少なくとも一部のボクセルについて前記形状特徴量を算出する形状特徴量算出部と、前記形状特徴量に基づき、三次元モデルが所定の形状として推定されたボクセル群を抽出する三次元形状抽出部と、前記ボクセル群を前記生体組織の三次元モデルにおける隆起形状を構成するボクセル群とする隆起形状検出部とを有する。
Description
明 細 書
医療用画像処理装置及び医療用画像処理方法
技術分野
[0001] 本発明は、医療用画像処理装置及び医療用画像処理方法に関し、特に、体腔内 の生体組織の三次元モデルにおける隆起形状を検出可能な医療用画像処理装置 及び医療用画像処理方法に関するものである。
背景技術
[0002] 内視鏡及び医療用画像処理装置等を具備して構成される内視鏡システムは、医療 分野等において広く用いられている。具体的には、内視鏡システムは、例えば、生体 としての体腔内に挿入される挿入部と、該揷入部の先端部に配置された対物光学系 と、該対物光学系により結像された体腔内の像を撮像して撮像信号として出力する 撮像部とを有して構成される内視鏡と、該撮像信号に基づき、表示部としてのモニタ 等に該体腔内の像を画像表示させるための処理を行う医療用画像処理装置とを具 備して構成されている。そして、ユーザは、表示部としてのモニタ等に画像表示され た体腔内の像に基づき、例えば、体腔内における被写体としての臓器等の観察を行
[0003] また、前述した構成を具備する内視鏡システムは、体腔内における被写体として、 例えば、大腸等の消化管粘膜の像を撮像することもまた可能である。そのため、ユー ザは、例えば、粘膜の色調、病変の形状及び粘膜表面の微細な構造等の様々な所 見を総合的に観察することができる。
[0004] さらに、近年においては、内視鏡により撮像された被写体の像の撮像信号に応じた 二次元画像のデータに基づき、該被写体の三次元モデルを生成することが可能であ る、例えば、日本国特開平 11— 337845号公報に記載されているような内視鏡装置 が提案されている。
[0005] 一方、三次元モデルにおいてポリープ等の病変部位を検出するための方法として 、例えば、 Shapelndex及び Curvednessといった、曲率に基づく形状特徴量を用い つつ該三次元モデルの形状の評価を行うことにより、該三次元モデルにおける病変
部位を検出可能である、 US Patent Application No. 20030223627に記載 されて ヽる方法が提案されて 、る。
[0006] 特に、内視鏡により撮像された隆起形状を有する生体組織の像としての、例えば、 大腸のひだまたはポリープ等の生体組織の像を含む二次元画像は、該内視鏡の視 野方向に対して不可視領域を含む場合が多い。そして、前述した二次元画像におけ る不可視領域は、一般的にォクルージョンと呼ばれており、正確な三次元モデルの 推定が困難な領域である。そのため、二次元画像のォクルージョンが発生している部 分においては、例えば、該二次元画像に基づく三次元モデルの対応する位置に推 定結果が存在しない、または、該二次元画像に基づく三次元モデルの対応する位置 に信頼性の低い推定結果が算出されるといったことが起こり得る。
[0007] その結果、従来の画像処理としての、例えば、日本国特開平 11— 337845号公報 の内視鏡装置において行われている画像処理により推定された被写体の三次元モ デルが用いられつつ該被写体の観察が行われた場合、ユーザがポリープ等の病変 部位を見つける際の負担が大きくなつてしまうと 、う課題が生じて 、る。
[0008] また、 US Patent Application No. 20030223627【こお!ヽて提案されて!ヽる ポリープ等の病変部位を検出するための方法は、例えば、 CT (Computed Tomog rapy)を用いた観察により得られる、ォクルージョンとなる部分が発生しない三次元モ デルに対する適用を想定した方法である。そのため、 US Patent Application No. 20030223627において提案されているポリープ等の病変部位を検出するた めの方法は、例えば、内視鏡を用いた観察により得られる、ォクルージョンとなる部分 を有する二次元画像に基づいて推定された三次元モデルに対して適用された場合 に、ポリープ等の病変部位の検出精度が低下してしまう。その結果、 US Patent A pplication No. 20030223627において提案されている方法力採用された内視 鏡システムが用いられつつ被写体の観察が行われた場合、ユーザがポリープ等の病 変部位を見つける際の負担が大きくなつてしまうという課題が生じている。
[0009] 本発明は、前述した点に鑑みてなされたものであり、ォクルージョンを有する二次元 画像に基づ 、て推定された三次元モデルを用いて観察が行われる場合に、ユーザ によるポリープ等の病変部位の見落としを防止可能であることにより、ユーザの負担
を軽減させることのできる医療用画像処理装置及び医療用画像処理方法を提供す ることを目的としている。
発明の開示
課題を解決するための手段
[0010] 本発明における第 1の医療用画像処理装置は、医療用撮像装置から入力される体 腔内の生体組織の像の二次元画像に基づき、該二次元画像のエッジを抽出するェ ッジ抽出部と、前記二次元画像に基づき、前記生体組織の三次元モデルを推定する 三次元モデル推定部と、前記三次元モデルの一部として推定された前記エッジが存 在する各ボタセルの位置に基づき、形状特徴量の算出対象となる所定のボタセル群 を抽出するボタセル抽出部と、前記ボタセル抽出部により抽出された前記所定のボタ セル群のうち、少なくとも一部のボタセルについて前記形状特徴量を算出する形状 特徴量算出部と、前記形状特徴量に基づき、前記少なくとも一部のボタセルのうち、 三次元モデルが所定の形状として推定されたボタセル群を抽出するとともに、該ボク セル群を、前記生体組織の三次元モデルにおける隆起形状を構成するボタセル群と して検出する隆起形状検出部と、を有することを特徴とする。
[0011] 本発明における第 2の医療用画像処理装置は、前記第 1の医療用画像処理装置 において、前記所定のボタセル群は、前記エッジが存在する各ボタセルに含まれる 所定の一のボタセルよりも、前記医療用撮像装置の視野方向手前側に存在するボタ セル群であることを特徴とする。
[0012] 本発明における第 3の医療用画像処理装置は、前記第 2の医療用画像処理装置 において、前記所定の一のボタセルは、前記エッジが存在する各ボタセルのうち、前 記視野方向の最も奥側に存在する一のボタセルであることを特徴とする。
[0013] 本発明における第 4の医療用画像処理装置は、前記第 1乃至第 3の医療用画像処 理装置において、前記所定の形状は、凸型形状であることを特徴とする。
[0014] 本発明における第 5の医療用画像処理装置は、入力される体腔内の生体組織の像 の二次元画像に基づき、該二次元画像のエッジを細線ィ匕して抽出するエッジ抽出部 と、前記エッジ抽出部が抽出したエッジ各々について、一のエッジの端部となる位置 の位置情報と、該位置情報に基づいて算出される、該一のエッジの形状に関する情
報である形状情報とを取得するエッジ情報取得部と、前記二次元画像に基づき、前 記生体組織の三次元モデルを推定する三次元モデル推定部と、前記位置情報と前 記形状情報とに基づき、前記生体組織の三次元モデルの奥行に関する情報である 奥行情報と、前記生体組織の三次元モデルの高さに関する情報である高さ情報とを 取得する三次元モデル情報取得部と、前記奥行情報と前記高さ情報とに基づき、前 記生体組織が隆起形状を有するか否かを検出する隆起形状検出部と、を有すること を特徴とする。
[0015] 本発明における第 6の医療用画像処理装置は、前記第 5の医療用画像処理装置 において、前記形状情報は、前記一のエッジの細線における中点の位置に関する情 報と、前記一のエッジの端部における 2つの端点を結ぶ直線の中点の位置に関する 情報とを有することを特徴とする。
[0016] 本発明における第 1の医療用画像処理方法は、医療用撮像装置から入力される体 腔内の生体組織の像の二次元画像に基づき、該二次元画像のエッジを抽出するェ ッジ抽出ステップと、前記二次元画像に基づき、前記生体組織の三次元モデルを推 定する三次元モデル推定ステップと、前記三次元モデルの一部として推定された前 記エッジが存在する各ボタセルの位置に基づき、形状特徴量の算出対象となる所定 のボタセル群を抽出するボタセル抽出ステップと、前記ボタセル抽出ステップにより抽 出された前記ボタセル群のうち、少なくとも一部のボタセルについて前記形状特徴量 を算出する形状特徴量算出ステップと、前記形状特徴量に基づき、前記少なくとも一 部のボタセルのうち、三次元モデルが所定の形状として推定されたボタセル群を抽出 するとともに、該ボクセル群を、前記生体組織の三次元モデルにおける隆起形状を 構成するボタセル群として検出する隆起形状検出ステップと、を有することを特徴と する。
[0017] 本発明における第 2の医療用画像処理方法は、前記第 1の医療用画像処理方法 において、前記所定のボタセル群は、前記エッジが存在する各ボタセルに含まれる 所定の一のボタセルよりも、前記医療用撮像装置の視野方向手前側に存在するボタ セル群であることを特徴とする。
[0018] 本発明における第 3の医療用画像処理方法は、前記第 2の医療用画像処理方法
において、前記所定の一のボタセルは、前記エッジが存在する各ボタセルのうち、前 記視野方向の最も奥側に存在する一のボタセルであることを特徴とする請求項 8に 記載の医療用画像処理方法。
[0019] 本発明における第 4の医療用画像処理方法は、前記第 1乃至第 3の医療用画像処 理方法において、前記所定の形状は、凸型形状であることを特徴とする。
[0020] 本発明における第 5の医療用画像処理方法は、入力される体腔内の生体組織の像 の二次元画像に基づき、該二次元画像のエッジを細線ィ匕して抽出するエッジ抽出ス テツプと、前記エッジ抽出ステップが抽出したエッジ各々について、一のエッジの端 部となる位置の位置情報と、該位置情報に基づいて算出される、該一のエッジの形 状に関する情報である形状情報とを取得するエッジ情報取得ステップと、前記二次 元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデル推定ステ ップと、前記位置情報と前記形状情報とに基づき、前記生体組織の三次元モデルの 奥行に関する情報である奥行情報と、前記生体組織の三次元モデルの高さに関す る情報である高さ情報とを取得する三次元モデル情報取得ステップと、前記奥行情 報と前記高さ情報とに基づき、前記生体糸且織が隆起形状を有する力否かを検出する 隆起形状検出ステップと、を有することを特徴とする。
[0021] 本発明における第 6の医療用画像処理方法は、前記第 5の医療用画像処理方法 において、前記形状情報は、前記一のエッジの細線における中点の位置に関する情 報と、前記一のエッジの端部における 2つの端点を結ぶ直線の中点の位置に関する 情報とを有することを特徴とする。
図面の簡単な説明
[0022] [図 1]本発明の各実施形態に係る医療用画像処理装置が用いられる内視鏡システム の全体構成の一例を示す図。
[図 2]図 1の内視鏡が管状器官内に挿入された場合の状態を示す模式図。
[図 3]図 2の状態において、内視鏡により撮像される管状器官及び生体組織の像を示 す模式図。
[図 4]生体組織の像の二次元画像におけるエッジ部を抽出する場合に、図 1の医療 用画像処理装置が第 1の実施形態において行う処理の手順を示すフローチャート。
[図 5]生体組織の三次元モデルにおける隆起形状を検出する場合に、図 1の医療用 画像処理装置が第 1の実施形態において行う処理の手順を示すフローチャート。
[図 6]図 1の医療用画像処理装置により推定された、生体組織の三次元モデルの一 例を示す図。
[図 7]図 5のフローチャートにおいて、生体組織の三次元モデルにおける隆起形状を 検出するための処理の対象となる領域の一例を示す図。
[図 8]図 5のフローチャートにおいて、生体組織の三次元モデルにおける隆起形状を 検出するための処理の対象となる領域の、図 7とは異なる例を示す図。
[図 9]図 1の医療用画像処理装置が第 2の実施形態において行う処理の手順の概要 を示すフローチャート。
[図 10]図 9における注目点抽出処理の具体的な手順を示すフローチャート。
[図 11]図 9における隆起形状検出処理の具体的な手順を示すフローチャート。
[図 12]図 1の医療用画像処理装置が第 2の実施形態の処理を行うことにより、生体組 織のエッジ部を含む二次元画像が細線ィ匕された場合の一例を示す図。
[図 13]図 12の細線ィ匕された二次元画像における、各注目点の位置関係の一例を示 す図。
[図 14]図 1の医療用画像処理装置が第 2の実施形態の処理を行うことにより得る、生 体組織の三次元モデルにおける各注目点の位置関係の一例を示す図。
[図 15]図 1の医療用画像処理装置が第 2の実施形態の処理を行うことにより得る、生 体組織の三次元モデルにおける各注目点の位置関係の、図 14とは異なる一例を示 す図。
発明を実施するための最良の形態
[0023] 以下、図面を参照して本発明の実施の形態を説明する。
[0024] (第 1の実施形態)
図 1から図 8は、本発明の第 1の実施形態に係るものである。図 1は、本実施形態に 係る医療用画像処理装置が用いられる内視鏡システムの全体構成の一例を示す図 である。図 2は、図 1の内視鏡が管状器官内に挿入された場合の状態を示す模式図 である。図 3は、図 2の状態において、内視鏡により撮像される管状器官及び生体組
織の像を示す模式図である。図 4は、生体組織の像の二次元画像におけるエッジ部 を抽出する場合に、図 1の医療用画像処理装置が行う処理の手順を示すフローチヤ ートである。図 5は、生体組織の三次元モデルにおける隆起形状を検出する場合に、 図 1の医療用画像処理装置が行う処理の手順を示すフローチャートである。図 6は、 図 1の医療用画像処理装置により推定された、生体組織の三次元モデルの一例を示 す図である。図 7は、図 5のフローチャートにおいて、生体糸且織の三次元モデルにお ける隆起形状を検出するための処理の対象となる領域の一例を示す図である。図 8 は、図 5のフローチャートにおいて、生体糸且織の三次元モデルにおける隆起形状を 検出するための処理の対象となる領域の、図 7とは異なる例を示す図である。
[0025] 内視鏡システム 1は、図 1に示すように、被写体を撮像するとともに、該被写体の像 の二次元画像を出力する医療用観察装置 2と、パーソナルコンピュータ等により構成 され、医療用観察装置 2から出力される二次元画像の映像信号に対して画像処理を 行うとともに、該画像処理を行った後の映像信号を画像信号として出力する医療用 画像処理装置 3と、医療用画像処理装置 3から出力される画像信号に基づく画像を 表示するモニタ 4とを有して要部が構成されて 、る。
[0026] また、医療用観察装置 2は、体腔内に挿入されるとともに、該体腔内の被写体を撮 像して撮像信号として出力する内視鏡 6と、内視鏡 6により撮像される被写体を照明 するための照明光を供給する光源装置 7と、内視鏡 6に対する各種制御を行うととも に、内視鏡 6から出力される撮像信号に対して信号処理を行い、二次元画像の映像 信号として出力するカメラコントロールユニット(以降、 CCUと略記する) 8と、 CCU8 力も出力される二次元画像の映像信号に基づき、内視鏡 6により撮像された被写体 の像を画像表示するモニタ 9とを有して要部が構成されている。
[0027] 医療用撮像装置としての内視鏡 6は、体腔内に挿入される挿入部 11と、挿入部 11 の基端側に設けられた操作部 12とを有して構成されている。また、挿入部 11内の基 端側から、挿入部 11内の先端側の先端部 14にかけての部分には、光源装置 7から 供給される照明光を伝送するためのライトガイド 13が挿通されている。
[0028] ライトガイド 13は、先端側が内視鏡 6の先端部 14に配置されるとともに、後端側が 光源装置 7に接続される。ライトガイド 13がこのような構成を有することにより、光源装
置 7から供給される照明光は、ライトガイド 13により伝送された後、挿入部 11の先端 部 14の先端面に設けられた、図示しない照明窓力も出射される。そして、図示しない 照明窓から照明光が出射されることにより、被写体としての生体組織等が照明される
[0029] 内視鏡 6の先端部 14には、図示しない照明窓に隣接する図示しない観察窓に取り 付けられた対物光学系 15と、対物光学系 15の結像位置に配置され、例えば、 CCD (電荷結合素子)等により構成される撮像素子 16とを有する撮像部 17が設けられて いる。このような構成により、対物光学系 15により結像された被写体の像は、撮像素 子 16により撮像された後、撮像信号として出力される。
[0030] 撮像素子 16は、信号線を介して CCU8に接続されている。そして、撮像素子 16は 、 CCU8から出力される駆動信号に基づいて駆動するとともに、 CCU8に対して撮像 信号を出力する。
[0031] また、 CCU8に入力された撮像信号は、 CCU8の内部に設けられた図示しない信 号処理回路において信号処理されることにより、二次元画像の映像信号として変換さ れて出力される。 CCU8から出力された二次元画像の映像信号は、モニタ 9及び医 療用画像処理装置 3に対して出力される。これにより、モニタ 9には、 CCU8から出力 される映像信号に基づく被写体の像が二次元の画像として表示される。
[0032] 医療用画像処理装置 3は、医療用観察装置 2から出力される二次元画像の映像信 号に対し、 AZD変換を行って出力する画像入力部 21と、画像入力部 21から出力さ れる映像信号に対して画像処理を行う、中央演算処理装置としての CPU22と、該画 像処理に関する処理プログラムが書き込まれた処理プログラム記憶部 23と、画像入 力部 21から出力される映像信号等を記憶する画像記憶部 24と、 CPU22の画像処 理結果としての画像データ等を記憶する情報記憶部 25とを有する。
[0033] また、医療用画像処理装置 3は、記憶装置インターフェース 26と、記憶装置インタ 一フェース 26を介して CPU22の画像処理結果としての画像データ等を記憶する、 記憶装置としてのハードディスク 27と、 CPU22の画像処理結果としての画像データ に基づき、該画像データをモニタ 4に画像表示するための表示処理を行うとともに、 該表示処理を行った後の画像データを画像信号として出力する表示処理部 28と、 C
PU22が行う画像処理におけるパラメータ及び医療用画像処理装置 3に対する操作 指示をユーザが入力可能な、キーボード等により構成される入力操作部 29とを有す る。そして、モニタ 4は、表示処理部 28から出力される画像信号に基づく画像を表示 する。
[0034] なお、医療用画像処理装置 3の画像入力部 21、 CPU22、処理プログラム記憶部 2 3、画像記憶部 24、情報記憶部 25、記憶装置インターフェース 26、表示処理部 28 及び入力操作部 29は、データバス 30を介して相互に接続されている。
[0035] 次に、内視鏡システム 1の作用について説明を行う。
[0036] まず、ユーザは、図 2に示すように、例えば、大腸等である管状器官 31内に内視鏡 6の挿入部 11を挿入する。そして、ユーザにより挿入部 11が管状器官 31に挿入され ると、例えば、管状器官 31の内壁に存在する病変部位である生体組織 31Aの像が、 先端部 14に設けられた撮像部 17により、図 3に示すような像として撮像される。そし て、図 3に示すような像として撮像部 17により撮像された、管状器官 31及び生体組 織 31Aの像は、撮像信号として CCU8に対して出力される。
[0037] CCU8は、図示しない信号処理回路において、撮像部 17の撮像素子 16から出力 される撮像信号に対して信号処理を行うことにより、該撮像信号を二次元画像の映像 信号として変換して出力する。そして、モニタ 9は、 CCU8から出力される映像信号に 基づき、管状器官 31及び生体組織 31 Aの像を、例えば、図 3に示すような二次元の 画像として表示する。また、 CCU8は、撮像部 17の撮像素子 16から出力される撮像 信号に対して信号処理を行うことにより得られた二次元画像の映像信号を、医療用 画像処理装置 3に対して出力する。
[0038] 医療用画像処理装置 3に対して出力された二次元画像の映像信号は、画像入力 部 21において AZD変換された後、 CPU22に入力される。
[0039] CPU22は、画像入力部 21から出力された二次元画像の映像信号と、処理プログ ラム記憶部 23に書き込まれた処理プログラムとに基づき、例えば、以降に述べるよう な処理を行うことにより、該二次元画像における生体組織 31Aのエッジ部を抽出する
[0040] まず、エッジ抽出部としての機能を有する CPU22は、画像入力部 21から出力され
た二次元画像の映像信号に基づき、例えば、該二次元画像の赤色成分に対してバ ンドパスフィルタを適用することにより、該二次元画像に含まれる全てのエッジ部を抽 出する(図 4のステップ Sl)。
[0041] その後、エッジ抽出部としての機能を有する CPU22は、抽出した全てのエッジ部を 細線ィ匕する(図 4のステップ S2)とともに、細線ィ匕した該全てのエッジ部のうち、一の エッジ部 Eの長さ Lを算出する(図 4のステップ S3)。さらに、 CPU22は、一のエッジ 部 Eの長さ Lが、閾値 thLlより長ぐかつ、閾値 thL2より短いか否かの判断を行う。
[0042] そして、 CPU22は、一のエッジ部 Eの長さ L力 所定の閾値 thLl以下の長さである こと、または、閾値 thL2以上であることを検出した場合(図 4のステップ S4)、後述す る図 4のステップ S 11に示す処理を行う。また、 CPU22は、一のエッジ部 Eの長さ が 、閾値 thLlより長ぐかつ、閾値 thL2より短いことを検出した場合(図 4のステップ S4 )、該一のエッジ部 Eを制御点 Cn(n= l, 2, · ··, N)により N等分する(図 4のステップ S5)。
[0043] さらに、 CPU22は、一のエッジ部 Eの中点 Ccから引いた法線 NCcを取得するととも に、各制御点 Cnから引いた N本の法線 NCnを取得する(図 4のステップ S6)。その 後、 CPU22は、 N本の法線 NCnのうち、法線 NCcと交わるものの本数 Naを検出す る(図 4のステップ S 7)。
[0044] また、 CPU22は、 N本の法線 NCnのうち、法線 NCcと交わるものの本数 Naが、閾 値 thaより多いか否かの判断を行う。そして、 CPU22は、法線 NCcと交わるものの本 数 Naが閾値 thaより多いことを検出した場合(図 4のステップ S8)、一のエッジ部 Eに 含まれるピクセル群 ipを生体組織 31 Aのエッジ部に含まれるピクセル群であると判断 し、該ピクセル群 ipが有する各ピクセルにおける変数 edge (i)の値を ONとする(図 4 のステップ S9)。さらに、そして、 CPU22は、法線 NCcと交わるものの本数 Naが閾値 tha以下であることを検出した場合(図 4のステップ S8)、一のエッジ部 Eに含まれるピ クセル群 ipを生体組織 31 Aのエッジ部に含まれるピクセル群ではな 、と判断し、該ピ クセル群 ipが有する各ピクセルにおける変数 edge (i)の値を OFFとする(図 4のステツ プ S10)。
[0045] CPU22は、抽出した全てのエッジ部に対し、処理が完了したか否かを判断する。
そして、 CPU22は、抽出した全てのエッジ部に対しての処理が完了していないことを 検出した場合(図 4のステップ S11)、他の一のエッジ部に対し、前述した、図 4のステ ップ S3から図 4のステップ S10までの処理を行う。また、 CPU22は、抽出した全ての エッジ部に対しての処理が完了したことを検出した場合(図 4のステップ S11)、二次 元画像における生体組織 31Aのエッジ部を抽出するための一連の処理を終了する。
[0046] その後、 CPU22は、前述した、二次元画像における生体組織 31 Aのエッジ部を抽 出するための一連の処理を行うことにより得た処理結果としての、二次元画像の各ェ ッジ部に含まれるピクセル群 ipにおける変数 edge (i)の値を、情報記憶部 25に一時 的に記憶させる。
[0047] そして、三次元モデル推定部としての機能を有する CPU22は、例えば、幾何学的 な変換等の処理により、画像入力部 21から出力された二次元画像の映像信号の輝 度情報等に基づき、管状器官 31及び生体組織 31Aの三次元モデルを推定する際 に必要となる画像データを取得するための画像処理を行う。換言すると、 CPU22は 、例えば、幾何学的な変換等の処理により、二次元画像における各ピクセルに対応 するボタセルを生成するとともに、該ボクセルを、三次元モデルを推定するための画 像データとして取得する。すなわち、ピクセル群 ipは、前述した処理によりボタセル群 ibとして変換される。
[0048] CPU22は、前述した処理により、図 3に示した生体組織 31 Aの三次元モデルを推 定するための画像データとして、変数 edge (i)が ONであるボタセル群 ibを含む平面 である生体組織境界部 3 lbのデータを得る。これにより、生体組織 31Aは、例えば、 z軸方向を内視鏡 6による観察時の視野方向とした場合、図 6に示すような形状を有 する三次元モデルとして推定される。
[0049] その後、 CPU22は、生体組織境界部 3 lbのデータに基づき、さらに、以降に述べ るような処理を行うことにより、生体組織 31 Aの三次元モデルにおける隆起形状を検 出する。
[0050] まず、 CPU22は、処理プログラム記憶部 23に書き込まれた処理プログラムと、生体 組織境界部 31bのデータとに基づき、変数 edge (i)が ONであるボタセル群 ibのうち 、内視鏡 6の視野方向の最も奥側に存在する所定の一のボタセルとして、 z座標が最
大である一のボタセルを抽出するとともに、該一のボタセルの z座標を Maxzとして設 定する(図 5のステップ S 21 )。
[0051] そして、 CPU22は、管状器官 31及び生体組織 31 Aの三次元モデルを推定するた めの画像データとして得た全てのボタセルのうち、前記一のボタセルの位置よりも内 視鏡 6の視野方向手前側に存在するボタセルとして、 z座標が Maxzよりも小さいボタ セル群 rbを検出する(図 5のステップ S 22)。
[0052] なお、前述したボタセル群 rbは、図 7に示す、隆起形状検出対象領域 AR内に含ま れるボタセル群、すなわち、形状特徴量の算出対象となるボタセル群であるとともに、
R個のボタセルを有するボタセル群であるとする。
[0053] また、前述した、図 5のステップ S21及び S22の処理を換言すると、ボタセル抽出部 としての機能を有する CPU22は、変数 edge (i)が ONであるボタセル群 ibが有する 各ボタセルの位置に基づき、隆起形状検出対象領域 AR内に含まれる、 R個のボクセ ルを有するボタセル群 rbを、形状特徴量の算出対象となる所定のボタセル群として 抽出する。
[0054] さらに、形状特徴量算出部としての機能を有する CPU22は、変数 aを 1に設定した
(図 5のステップ S23)後、ボタセル群 rbが有する R個のボタセルのうち、一のボタセル である Ba (a= l, 2, · ··, R— 1, R)を抽出するとともに(図 5のステップ S24)、該一の ボタセル Baにおける形状特徴量として、 Shapelndex値 SBa及び Curvedness値 C Baを算出する(図 5のステップ S25)。
[0055] なお、前述した Shapelndex値及び Curvedness値は、例えば、 US Patent Ap plication No. 20030223627に記載されている方法と同様の方法を用いることに より算出可能である。そのため、本実施形態においては、一のボタセル Baにおける S hapelndex値及び Curvedness値の算出方法に関しては、説明を省略する。
[0056] また、 CPU22は、 Shapelndex値 SBaと、予め設定された Shapelndex値の閾値 Sthとの比較を行うとともに、 Curvedness値 CBaと、予め設定された Curvedness値 の閾値 Cthとの比較を行う。換言すると、隆起形状検出部としての機能を有する CP U22は、前述した処理を行うことにより、生体組織 31 Aが隆起形状であるか否かを検 出するための処理として、三次元モデルが凸型形状と推定されたボタセル群を抽出
する処理を行う。なお、図 5に示す一連の処理においては、隆起形状として凸型形状 を有する生体組織 31Aを検出するために、例えば、閾値 Sthは 0. 9として設定され ているとし、また、閾値 Cthは 0. 2として設定されているとする。
[0057] そして、隆起形状検出部としての機能を有する CPU22は、 Shapelndex値 SBa力 S 閾値 Sthより大きく、かつ、 Curvedness値 CBaが閾値 Cthより大きいことを検出した 場合(図 5のステップ S26)、一のボタセル Baを隆起形状の一部を構成するボタセル であると判断し、該一のボタセル Baにおける変数 ryuuki(Ba)の値を ONとする(図 5 のステップ S27)。
[0058] また、 CPU22は、 Shapelndex値 SBaが閾値 Sth以下であること、または、 Curve dness値 CBaが閾値 Cth以下であることを検出した場合(図 5のステップ S26)、一の ボタセル Baを隆起形状の一部を構成するボタセルではな 、と判断し、該一のボクセ ル Baにおける変数 ryuuki (Ba)の値を OFFとする(図 5のステップ S28)。
[0059] その後、 CPU22は、 R個のボタセル全てにおいて、前述した処理が行われたかどう 力 すなわち、変数 a =Rである力否かの判定を行う。
[0060] そして、 CPU22は、 a=Rではないことを検知した場合(図 5のステップ S29)、変数 iに 1を加える処理を行った(図 5のステップ S30)後、前述した、図 5のステップ S24か らステップ S29に示す処理を再度行う。
[0061] また、 CPU22は、 a=Rであることを検知した場合(図 5のステップ S29)、生体組織 31Aの三次元モデルにおける隆起形状を検出するための一連の処理を完了する。
[0062] そして、 CPU22は、例えば、 ryuuki (Ba)の値が ONであるボタセルを有する三次 元モデルにぉ 、て、生体組織 31Aが隆起形状であることを示すための文字列または 着色等を該三次元モデルに重畳させるような制御を表示処理部 28に対して行う。こ れにより、モニタ 4には、隆起形状を有する生体組織 31 Aをユーザが容易に発見可 能であるような、管状器官 31及び生体組織 31 Aの三次元モデルが画像表示される。
[0063] また、 CPU22は、 ryuuki (Ba)の値が ONである各ボタセルの位置に基づき、該各 ボタセルの位置に対応する位置に存在する二次元画像上の各ピクセルを検出すると ともに、生体組織 31Aが隆起形状であることを示すための文字列または着色等を、 該各ピクセルを有する該二次元画像に重畳させるような制御を表示処理部 28に対し
て行うものであっても良い。
[0064] なお、前述した一連の処理にぉ 、て、隆起形状検出対象領域 ARは、 z座標が Ma xzよりも小さいボタセル群を含む領域として定められるものに限らず、例えば、図 8に 示すように、生体組織境界部 3 lb近傍に存在するボタセル群のみを含む領域として 定められるものであっても良い。
[0065] 本実施形態の医療用画像処理装置 3は、以上に述べた一連の処理を行うことによ り、ユーザによるポリープ等の病変部位の見落としを防止可能であることにより、ユー ザの負担を軽減させることができる。
[0066] また、本実施形態の医療用画像処理装置 3は、前述した、生体組織 31Aの隆起形 状を検出するための処理を、三次元モデルにおける隆起形状検出対象領域 AR内 に含まれるボタセル群に対してのみ行う。そのため、本実施形態の医療用画像処理 装置 3は、生体組織 31Aの隆起形状の検出を高速に行うことができる。
[0067] (第 2の実施形態)
図 9から図 15は、本発明の第 2の実施形態に係るものである。なお、第 1の実施形 態と同様の構成を持つ部分については、詳細説明は省略する。また、第 1の実施形 態と同様の構成要素については、同一の符号を用いて説明は省略する。さらに、本 実施形態に用いる内視鏡システム 1の構成は、第 1の実施形態と同様である。
[0068] 図 9は、図 1の医療用画像処理装置が第 2の実施形態において行う処理の手順の 概要を示すフローチャートである。図 10は、図 9における注目点抽出処理の具体的 な手順を示すフローチャートである。図 11は、図 9における隆起形状検出処理の具 体的な手順を示すフローチャートである。図 12は、図 1の医療用画像処理装置が第 2の実施形態の処理を行うことにより、生体組織のエッジ部を含む二次元画像が細線 化された場合の一例を示す図である。図 13は、図 12の細線ィ匕された二次元画像に おける、各注目点の位置関係の一例を示す図である。図 14は、図 1の医療用画像処 理装置が第 2の実施形態の処理を行うことにより得る、生体組織の三次元モデルに おける各注目点の位置関係の一例を示す図である。図 15は、図 1の医療用画像処 理装置が第 2の実施形態の処理を行うことにより得る、生体組織の三次元モデルに おける各注目点の位置関係の、図 14とは異なる一例を示す図である。
[0069] まず、エッジ抽出部としての機能を有する CPU22は、画像入力部 21から出力され る、管状器官 31及び生体組織 31Aが撮像された二次元画像の映像信号に基づき、 例えば、該二次元画像の赤色成分に対してバンドパスフィルタを適用することにより、 該二次元画像に含まれる全てのエッジ部を抽出する(図 9のステップ S31)。
[0070] その後、エッジ抽出部としての機能を有する CPU22は、抽出した全てのエッジ部を 2値化及び細線化する(図 9のステップ S32)ととも〖こ、細線化された二次元画像にお ける注目点抽出処理を行う(図 9のステップ S33)。
[0071] 具体的には、エッジ部が細線ィ匕された二次元画像として、例えば、図 12に示すよう な生体組織 31 Aの一のエッジ部 E1が含まれている画像を取得した場合、まず、エツ ジ情報取得部としての機能を有する CPU22は、該一のエッジ部 E1の細線における 、一のエッジの端部となる位置の位置情報として、端点 K1及び K2を取得する(図 10 のステップ S331)。
[0072] 次に、エッジ情報取得部としての機能を有する CPU22は、一のエッジ部 E1の形状 に関する形状情報として、該一のエッジ部 E1の細線の線分長 LE1を算出する(図 10 のステップ S332)とともに、 LEI X 1Z2となる中点 KEを取得する(図 10のステップ S 333) o
[0073] また、エッジ情報取得部としての機能を有する CPU22は、一のエッジ部 E1の形状 に関する形状情報として、端点 K1と端点 K2とを結ぶ直線 F1を抽出し、該直線の長 さ LF1を算出する(図 10のステップ S334)とともに、 LF1 X 1Z2となる中点 KFを取 得する(図 10のステップ S335)。
[0074] そして、 CPU22は、端点 Kl、端点 Κ2、中点 KE及び中点 KFを、一のエッジ部 E1 における注目点として抽出する(図 10のステップ S336)。なお、一のエッジ部 E1に おける端点 Kl、端点 Κ2、中点 KE及び中点 KFは、例えば、各々が図 13に示すよう な位置関係を有する注目点として抽出される。
[0075] CPU22は、抽出した全てのエッジ部に対し、処理が完了したか否かを判断する。
そして、 CPU22は、抽出した全てのエッジ部に対しての処理が完了していないことを 検出した場合(図 10のステップ S337)、他の一のエッジ部に対し、前述した、図 10の ステップ S331から図 10のステップ S336までの処理を行う。また、 CPU22は、抽出
した全てのエッジ部に対しての処理が完了したことを検出した場合 (図 10のステップ S337)、二次元画像における注目点抽出処理を終了する。
[0076] そして、 CPU22は、前述した、二次元画像における注目点抽出処理を行うことによ り得た処理結果として、該二次元画像の各エッジ部が有する注目点である、端点 K1 、端点 K2、中点 ΚΕ及び中点 KFのデータを、情報記憶部 25に一時的に記憶させる
[0077] 三次元モデル推定部としての機能を有する CPU22は、二次元画像における注目 点抽出処理を行った後、例えば、幾何学的な変換等の処理により、画像入力部 21か ら出力された二次元画像の映像信号の輝度情報等に基づき、管状器官 31及び生体 組織 31 Αの三次元モデルを推定する際に必要となる画像データを取得するための 画像処理を行う(図 9のステップ S34)。換言すると、 CPU22は、例えば、幾何学的な 変換等の処理により、二次元画像における各ピクセルに対応するボタセルを生成す るとともに、該ボクセルを、三次元モデルを推定するための画像データとして取得す る。
[0078] さらに、 CPU22は、情報記憶部 25に記憶されている各注目点のデータである、端 点 Kl、端点 Κ2、中点 ΚΕ及び中点 KFのデータに基づき、生体組織 31Aの一のエツ ジ部 E1の三次元モデルにおける該各注目点の座標を推定する(図 9のステップ S35 ) oなお、一のエッジ部 E1の三次元モデルにおける端点 Kl、端点 Κ2、中点 ΚΕ及び 中点 KFは、前述した処理により、例えば、各々が図 14に示すような位置関係を有す るとして推定される。
[0079] そして、 CPU22は、前述した処理により推定された各注目点の座標に基づき、さら に、以降に述べるような隆起形状検出処理を行うことにより、生体組織 31Aの三次元 モデルにおける隆起形状を検出する(図 9のステップ S 36)。
[0080] 具体的には、まず、三次元モデル情報取得部としての機能を有する CPU22は、生 体組織 31 Αの三次元モデルにおける奥行情報として、一のエッジ部 E1の三次元モ デル〖こおける、中点 KFの z座標 KFzと、端点 K1の z座標 Klzと、端点 K2の z座標 K 2zとを抽出する(図 11のステップ S361)。
[0081] 次に、隆起形状検出部としての機能を有する CPU22は、中点 KFの z座標 KFzの
値と、端点 Klの z座標 Klzの値及び端点 K2の z座標 K2zの値との比較を行う。そし て、 CPU22は、中点 KFの ζ座標 KFzの値が端点 K1の ζ座標 Klzの値よりも小さく、 かつ、中点 KFの z座標 KFzの値が端点 K2の z座標 K2zの値よりも小さ!/、ことを検出 した場合(図 11のステップ S362)、後述する、図 11のステップ S363に示す処理をさ らに行う。また、 CPU22は、中点 KFの z座標 KFzの値が端点 K1の z座標 Klzの値 以上であるか、または、中点 KFの z座標 KFzの値が端点 K2の z座標 K2zの値以上 であることを検出した場合(図 11のステップ S362)、生体組織 31 Αの一のエッジ部 Ε 1を、隆起形状を構成するエッジ部ではないと判断し、後述する、図 11のステップ S3 67に示す処理を行う。
[0082] 三次元モデル情報取得部としての機能を有する CPU22は、中点 KFの z座標 KFz の値が端点 Klの z座標 Klzの値よりも小さぐかつ、中点 KFの z座標 KFzの値が端 点 K2の z座標 K2zの値よりも小さいことを検出した場合(図 11のステップ S362)、三 次元モデルにおける端点 K1と端点 Κ2とを結ぶ直線 F2を抽出するとともに、三次元 モデルにおける中点 ΚΕと該直線 F2とを結ぶ直線である垂直線 F3を取得する(図 1 1のステップ S363)。その後、 CPU22は、生体糸且織 31 Aの三次元モデルにおける高 さ情報として、垂直線 F3の線分長 LF3を取得する(図 11のステップ S364)とともに、 該線分長 LF3と、該線分長 LF3の閾値 thLF3との比較を行う。
[0083] そして、隆起形状検出部としての機能を有する CPU22は、線分長 LF3の値が閾値 thLF3よりも大き ヽことを検出した場合(図 11のステップ S 365)、生体組織 31 Aの一 のエッジ部 E1を、隆起形状を構成するエッジ部であると判断するとともに、該一のェ ッジ部 E1が有するボタセル群全てにおける変数 ryuuki (E1)の値を ONとする(図 11 のステップ S366)。また、 CPU22は、線分長 LF3の値が閾値 thLF3以下であること を検出した場合(図 11のステップ S365)、生体組織 31Aの一のエッジ部 E1を、隆起 形状を構成するエッジ部ではないと判断するとともに、該一のエッジ部 E1が有するボ クセル群全てにおける変数 ryuuki (El)の値を OFFとする(図 11のステップ S 367)。
[0084] CPU22は、抽出した全てのエッジ部に対し、処理が完了したか否かを判断する。
そして、 CPU22は、抽出した全てのエッジ部に対しての処理が完了していないことを 検出した場合(図 11のステップ S368)、他の一のエッジ部に対し、前述した、図 11の
ステップ S361から図 11のステップ S367までの処理を行う。また、 CPU22は、抽出 した全てのエッジ部に対しての処理が完了したことを検出した場合(図 11のステップ S368)、生体組織 31Aの三次元モデルにおける隆起形状検出処理を終了する。
[0085] そして、 CPU22は、前述した、図 11に示す隆起形状検出処理が終了した時点に おいて、図 9に示す一連の画像処理を完了する。
[0086] なお、前述した隆起形状検出処理の、図 11のステップ S362に示す処理において 、 CPU22は、各注目点の z座標の値に基づく比較及び判断を行うものに限らず、例 えば、以降に記すような、各注目点の赤色成分値に基づく比較及び判断を行うもの であっても良い。
[0087] 具体的には、まず、 CPU22は、中点 KFにおける赤色成分値 RKFと、端点 K1に おける赤色成分値 RK1と、端点 K2における赤色成分値 RK2とを抽出する。そして、 CPU22は、中点 KFにおける赤色成分値 RKFが端点 K1における赤色成分値 RK1 よりも大きぐかつ、中点 KFにおける赤色成分値 RKFが端点 K2における赤色成分 値 RK2よりも大きいことを検出した場合、前述した、図 11のステップ S363の処理をさ らに行う。また、 CPU22は、中点 KFにおける赤色成分値 RKFが端点 K1における赤 色成分値 RK1以下であるか、または、中点 KFにおける赤色成分値 RKFが端点 K2 における赤色成分値 RK2以下であることを検出した場合、生体組織 31Aの一のエツ ジ部 E1を、隆起形状を構成するエッジ部ではないと判断し、前述した、図 11のステツ プ S367の処理を行う。
[0088] そして、 CPU22は、例えば、 ryuuki (El)の値が ONであるボタセルを有する三次 元モデルにぉ 、て、生体組織 31Aが隆起形状であることを示すための文字列または 着色等を該三次元モデルに重畳させるような制御を表示処理部 28に対して行う。こ れにより、モニタ 4には、隆起形状を有する生体組織 31 Aをユーザが容易に発見可 能であるような、管状器官 31及び生体組織 31 Aの三次元モデルが画像表示される。
[0089] また、 CPU22は、 ryuuki (E1)の値が ONである各ボタセルの位置に基づき、該各 ボタセルの位置に対応する位置に存在する二次元画像上の各ピクセルを検出すると ともに、生体組織 31Aが隆起形状であることを示すための文字列または着色等を、 該各ピクセルを有する該二次元画像に重畳させるような制御を表示処理部 28に対し
て行うものであっても良い。
[0090] 本実施形態の医療用画像処理装置 3は、以上に述べた一連の処理を行うことによ り、ユーザによるポリープ等の病変部位の見落としを防止可能であることにより、ユー ザの負担を軽減させることができる。
[0091] また、本実施形態の医療用画像処理装置 3は、前述した、生体組織 31Aが隆起形状 であるか否かを、隆起形状特有の値を算出しつつ定量的に検出する。そのため、本 実施形態の医療用画像処理装置 3は、生体組織 31Aの隆起形状の検出を精度良く 行うことができる。
[0092] なお、本発明は、上述した各実施形態に限定されるものではなぐ発明の趣旨を逸 脱しな 、範囲内にお 、て種々の変更や応用が可能であることは勿論である。
[0093] 本出願は、 2006年 3月 8日に日本国に出願された特願 2006— 63234号を優先 権主張の基礎として出願するものであり、上記の開示内容は、本願明細書、請求の 範囲、図面に引用されたものとする。
Claims
[1] 医療用撮像装置力も入力される体腔内の生体組織の像の二次元画像に基づき、 該二次元画像のエッジを抽出するエッジ抽出部と、
前記二次元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデ ル推定部と、
前記三次元モデルの一部として推定された前記エッジが存在する各ボタセルの位 置に基づき、形状特徴量の算出対象となる所定のボタセル群を抽出するボタセル抽 出部と、
前記ボタセル抽出部により抽出された前記所定のボタセル群のうち、少なくとも一部 のボタセルについて前記形状特徴量を算出する形状特徴量算出部と、
前記形状特徴量に基づき、前記少なくとも一部のボタセルのうち、三次元モデルが 所定の形状として推定されたボタセル群を抽出するとともに、該ボクセル群を、前記 生体組織の三次元モデルにおける隆起形状を構成するボタセル群として検出する隆 起形状検出部と、
を有することを特徴とする医療用画像処理装置。
[2] 前記所定のボタセル群は、前記エッジが存在する各ボタセルに含まれる所定の一 のボタセルよりも、前記医療用撮像装置の視野方向手前側に存在するボタセル群で あることを特徴とする請求項 1に記載の医療用画像処理装置。
[3] 前記所定の一のボタセルは、前記エッジが存在する各ボタセルのうち、前記視野方 向の最も奥側に存在する一のボタセルであることを特徴とする請求項 2に記載の医療 用画像処理装置。
[4] 前記所定の形状は、凸型形状であることを特徴とする請求項 1乃至請求項 3のいず れかーに記載の医療用画像処理装置。
[5] 入力される体腔内の生体組織の像の二次元画像に基づき、該二次元画像のエツ ジを細線ィ匕して抽出するエッジ抽出部と、
前記エッジ抽出部が抽出したエッジ各々について、一のエッジの端部となる位置の 位置情報と、該位置情報に基づいて算出される、該一のエッジの形状に関する情報 である形状情報とを取得するエッジ情報取得部と、
前記二次元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデ ル推定部と、
前記位置情報と前記形状情報とに基づき、前記生体組織の三次元モデルの奥行 に関する情報である奥行情報と、前記生体組織の三次元モデルの高さに関する情 報である高さ情報とを取得する三次元モデル情報取得部と、
前記奥行情報と前記高さ情報とに基づき、前記生体組織が隆起形状を有するか否 かを検出する隆起形状検出部と、
を有することを特徴とする医療用画像処理装置。
[6] 前記形状情報は、前記一のエッジの細線における中点の位置に関する情報と、前 記一のエッジの端部における 2つの端点を結ぶ直線の中点の位置に関する情報とを 有することを特徴とする請求項 5に記載の医療用画像処理装置。
[7] 医療用撮像装置力 入力される体腔内の生体組織の像の二次元画像に基づき、 該二次元画像のエッジを抽出するエッジ抽出ステップと、
前記二次元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデ ル推定ステップと、
前記三次元モデルの一部として推定された前記エッジが存在する各ボタセルの位 置に基づき、形状特徴量の算出対象となる所定のボタセル群を抽出するボタセル抽 出ステップと、
前記ボタセル抽出ステップにより抽出された前記ボタセル群のうち、少なくとも一部 のボタセルについて前記形状特徴量を算出する形状特徴量算出ステップと、 前記形状特徴量に基づき、前記少なくとも一部のボタセルのうち、三次元モデルが 所定の形状として推定されたボタセル群を抽出するとともに、該ボクセル群を、前記 生体組織の三次元モデルにおける隆起形状を構成するボタセル群として検出する隆 起形状検出ステップと、
を有することを特徴とする医療用画像処理方法。
[8] 前記所定のボタセル群は、前記エッジが存在する各ボタセルに含まれる所定の一 のボタセルよりも、前記医療用撮像装置の視野方向手前側に存在するボタセル群で あることを特徴とする請求項 7に記載の医療用画像処理方法。
[9] 前記所定の一のボタセルは、前記エッジが存在する各ボタセルのうち、前記視野方 向の最も奥側に存在する一のボタセルであることを特徴とする請求項 8に記載の医療 用画像処理方法。
[10] 前記所定の形状は、凸型形状であることを特徴とする請求項 7乃至請求項 9の 、ず れかーに記載の医療用画像処理方法。
[11] 入力される体腔内の生体組織の像の二次元画像に基づき、該二次元画像のエツ ジを細線ィ匕して抽出するエッジ抽出ステップと、
前記エッジ抽出ステップが抽出したエッジ各々について、一のエッジの端部となる 位置の位置情報と、該位置情報に基づいて算出される、該一のエッジの形状に関す る情報である形状情報とを取得するエッジ情報取得ステップと、
前記二次元画像に基づき、前記生体組織の三次元モデルを推定する三次元モデ ル推定ステップと、
前記位置情報と前記形状情報とに基づき、前記生体組織の三次元モデルの奥行 に関する情報である奥行情報と、前記生体組織の三次元モデルの高さに関する情 報である高さ情報とを取得する三次元モデル情報取得ステップと、
前記奥行情報と前記高さ情報とに基づき、前記生体組織が隆起形状を有するか否 かを検出する隆起形状検出ステップと、
を有することを特徴とする医療用画像処理方法。
[12] 前記形状情報は、前記一のエッジの細線における中点の位置に関する情報と、前 記一のエッジの端部における 2つの端点を結ぶ直線の中点の位置に関する情報とを 有することを特徴とする請求項 11に記載の医療用画像処理方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP07714003.6A EP1992274B1 (en) | 2006-03-08 | 2007-02-09 | Medical image processing device and medical image processing method |
| CN2007800075843A CN101394780B (zh) | 2006-03-08 | 2007-02-09 | 医疗用图像处理装置以及医疗用图像处理方法 |
| US12/204,330 US8086005B2 (en) | 2006-03-08 | 2008-09-04 | Medical image processing apparatus and medical image processing method |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006063234A JP4981336B2 (ja) | 2006-03-08 | 2006-03-08 | 医療用画像処理装置及び医療用画像処理方法 |
| JP2006-063234 | 2006-03-08 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US12/204,330 Continuation US8086005B2 (en) | 2006-03-08 | 2008-09-04 | Medical image processing apparatus and medical image processing method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2007102296A1 true WO2007102296A1 (ja) | 2007-09-13 |
Family
ID=38474738
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2007/052347 Ceased WO2007102296A1 (ja) | 2006-03-08 | 2007-02-09 | 医療用画像処理装置及び医療用画像処理方法 |
Country Status (4)
| Country | Link |
|---|---|
| EP (1) | EP1992274B1 (ja) |
| JP (1) | JP4981336B2 (ja) |
| CN (1) | CN101394780B (ja) |
| WO (1) | WO2007102296A1 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101996329B (zh) * | 2010-11-17 | 2012-10-31 | 沈阳东软医疗系统有限公司 | 一种对血管形变区域的检测装置和方法 |
| JP6298035B2 (ja) | 2015-12-17 | 2018-03-20 | ファナック株式会社 | モデル生成装置、位置姿勢算出装置、およびハンドリングロボット装置 |
| CN106909770B (zh) * | 2015-12-21 | 2020-11-03 | 佳能株式会社 | 医疗图像处理装置、其控制方法和存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0742536A1 (en) | 1995-05-11 | 1996-11-13 | Agfa-Gevaert N.V. | Method of recognising one or more irradiation fields |
| JPH11337845A (ja) * | 1998-05-25 | 1999-12-10 | Mitsubishi Electric Corp | 内視鏡装置 |
| JP2005506140A (ja) * | 2001-10-16 | 2005-03-03 | ザ・ユニバーシティー・オブ・シカゴ | コンピュータ支援の3次元病変検出方法 |
| JP2005177477A (ja) * | 2003-12-15 | 2005-07-07 | Siemens Ag | カテーテル装置 |
| WO2005107601A2 (en) | 2004-05-06 | 2005-11-17 | Focus Surgery, Inc. | Method and apparatus for the selective treatment of tissue |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2324602C (en) * | 1998-03-20 | 2012-07-24 | Barbara Ann Karmanos Cancer Institute | Multidimensional detection and characterization of pathologic tissues |
| US7043064B2 (en) * | 2001-05-04 | 2006-05-09 | The Board Of Trustees Of The Leland Stanford Junior University | Method for characterizing shapes in medical images |
| US20020164061A1 (en) * | 2001-05-04 | 2002-11-07 | Paik David S. | Method for detecting shapes in medical images |
| US6726849B2 (en) * | 2002-02-01 | 2004-04-27 | Construction Research & Technology Gmbh | Method of dusting coal mine surfaces |
| US20040101184A1 (en) * | 2002-11-26 | 2004-05-27 | Radhika Sivaramakrishna | Automatic contouring of tissues in CT images |
-
2006
- 2006-03-08 JP JP2006063234A patent/JP4981336B2/ja not_active Expired - Fee Related
-
2007
- 2007-02-09 EP EP07714003.6A patent/EP1992274B1/en not_active Not-in-force
- 2007-02-09 CN CN2007800075843A patent/CN101394780B/zh not_active Expired - Fee Related
- 2007-02-09 WO PCT/JP2007/052347 patent/WO2007102296A1/ja not_active Ceased
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0742536A1 (en) | 1995-05-11 | 1996-11-13 | Agfa-Gevaert N.V. | Method of recognising one or more irradiation fields |
| JPH11337845A (ja) * | 1998-05-25 | 1999-12-10 | Mitsubishi Electric Corp | 内視鏡装置 |
| JP2005506140A (ja) * | 2001-10-16 | 2005-03-03 | ザ・ユニバーシティー・オブ・シカゴ | コンピュータ支援の3次元病変検出方法 |
| JP2005177477A (ja) * | 2003-12-15 | 2005-07-07 | Siemens Ag | カテーテル装置 |
| WO2005107601A2 (en) | 2004-05-06 | 2005-11-17 | Focus Surgery, Inc. | Method and apparatus for the selective treatment of tissue |
Non-Patent Citations (2)
| Title |
|---|
| PENG Y.-J. ET AL.: "A new tissue segmentation algorithm in 3d data based on boundary model and local character structure", PROCEEDINGS. SPIE-THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING, vol. 6044, 2005, pages 60441J.1 - 60441J.6, XP003017616 * |
| See also references of EP1992274A4 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101394780B (zh) | 2010-09-29 |
| CN101394780A (zh) | 2009-03-25 |
| JP2007236630A (ja) | 2007-09-20 |
| EP1992274A4 (en) | 2013-05-01 |
| EP1992274B1 (en) | 2016-04-13 |
| EP1992274A1 (en) | 2008-11-19 |
| JP4981336B2 (ja) | 2012-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7830378B2 (en) | Medical image processing apparatus and medical image processing method | |
| JP5918548B2 (ja) | 内視鏡画像診断支援装置およびその作動方法並びに内視鏡画像診断支援プログラム | |
| US8086005B2 (en) | Medical image processing apparatus and medical image processing method | |
| JP5676058B1 (ja) | 内視鏡システム及び内視鏡システムの作動方法 | |
| JP5276225B2 (ja) | 医用画像処理装置及び医用画像処理装置の作動方法 | |
| JP5771757B2 (ja) | 内視鏡システム及び内視鏡システムの作動方法 | |
| JP4902735B2 (ja) | 医療用画像処理装置及び医療用画像処理方法 | |
| US10939800B2 (en) | Examination support device, examination support method, and examination support program | |
| JP5078486B2 (ja) | 医療用画像処理装置及び医療用画像処理装置の作動方法 | |
| US12433478B2 (en) | Processing device, endoscope system, and method for processing captured image | |
| US20220358750A1 (en) | Learning device, depth information acquisition device, endoscope system, learning method, and program | |
| WO2019130868A1 (ja) | 画像処理装置、プロセッサ装置、内視鏡システム、画像処理方法、及びプログラム | |
| EP1992273A1 (en) | Medical image processing device and medical image processing method | |
| EP1992274B1 (en) | Medical image processing device and medical image processing method | |
| JP5148096B2 (ja) | 医療用画像処理装置及び医療用画像処理装置の作動方法 | |
| JP2008093213A (ja) | 医療用画像処理装置及び医療用画像処理方法 | |
| JP2008023266A (ja) | 医療用画像処理装置及び医療用画像処理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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: 2007714003 Country of ref document: EP |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 200780007584.3 Country of ref document: CN |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |