WO2015111308A1 - ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム - Google Patents

３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム Download PDF

Info

Publication number
WO2015111308A1
WO2015111308A1 PCT/JP2014/081977 JP2014081977W WO2015111308A1 WO 2015111308 A1 WO2015111308 A1 WO 2015111308A1 JP 2014081977 W JP2014081977 W JP 2014081977W WO 2015111308 A1 WO2015111308 A1 WO 2015111308A1
Authority
WO
WIPO (PCT)
Prior art keywords
medical image
curved surface
display control
dimensional medical
region
Prior art date
Application number
PCT/JP2014/081977
Other languages
English (en)
French (fr)
Inventor
王　彩華
Original Assignee
富士フイルム株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士フイルム株式会社 filed Critical 富士フイルム株式会社
Publication of WO2015111308A1 publication Critical patent/WO2015111308A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/12—Arrangements for detecting or locating foreign bodies
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/136—Segmentation; Edge detection involving thresholding
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR 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/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR 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/30061—Lung

Definitions

  • the present invention relates to a three-dimensional medical image display control apparatus that displays a three-dimensional medical image of a lung, an operating method thereof, and a three-dimensional medical image display control program, and in particular, a curved surface 3 set inside a lung field. It relates to dimension display.
  • Interstitial pneumonia is one of the diffuse lung diseases, and in many cases, lesions are distributed over a large area of the lung. Therefore, in order to diagnose the degree of progression of interstitial pneumonia, it is necessary to grasp the lesion that spreads throughout the lungs.
  • Patent Document 1 proposes a method of displaying a lung nodule and its peripheral region by extracting a lung nodule, specifying a peripheral region continuous with the lung nodule, and performing volume rendering.
  • Patent Document 2 proposes a method of acquiring a three-dimensional medical image of a lung field region and displaying it three-dimensionally.
  • the CT image can clearly grasp the lesion of interstitial pneumonia, and is considered effective for diagnosis.
  • interstitial pneumonia lesions varies greatly from mild to severe, and is distributed over a wide range as described above, so it is difficult to extract the lesion area of interstitial pneumonia.
  • the lesion of interstitial pneumonia is a small change in concentration or a pattern of small cavities surrounded by a thin wall like a honeycomb, it is a common example as described in Patent Document 1, for example.
  • a lesion of interstitial pneumonia is displayed by simple volume rendering, other thin blood vessels, bronchi and the like are also displayed at the same time, so that useful diagnostic information cannot be obtained.
  • the present invention provides a three-dimensional medical image display control device capable of appropriately displaying the presence or absence of a lesion of interstitial pneumonia and the degree of progression thereof, an operating method thereof, and a three-dimensional medical image display.
  • the object is to provide a control program.
  • the three-dimensional medical image display control device of the present invention includes a curved surface setting unit that sets a curved surface along the surface of the lung field inside the lung field in the three-dimensional medical image obtained by imaging a region including the lung field, and a curved surface setting unit And a display control unit that displays the curved surface in a three-dimensional manner using pixel values on the curved surface set by the above.
  • the curved surface can be a surface equidistant from the surface of the lung field.
  • the curved surface setting unit can be set by changing the position of the curved surface from the surface of the lung field.
  • the curved surface setting unit can be set by switching a plurality of curved surfaces having different distances from the lung field surface.
  • the display control unit can display the outline or surface of the lung field.
  • an alveolar region extracting unit that extracts an alveolar region from a three-dimensional medical image, and an alveolar region having a size greater than or equal to a threshold value among a plurality of alveolar regions extracted by the alveolar region extracting unit
  • a display alveolar region identifying unit that identifies the alveolar region may be provided, and the display control unit may display only the alveolar region to be displayed among the plurality of alveolar regions extracted by the alveolar region extracting unit.
  • an alveolar region extracting unit that extracts an alveolar region from a three-dimensional medical image, and a display lung that identifies an alveolar region having a size greater than or equal to a threshold value among a plurality of alveolar regions extracted by the alveolar region extracting unit
  • An alveolar region specifying unit, and the display control unit includes a plurality of alveolar regions extracted by the alveolar region extracting unit, wherein the specified alveolar region is defined as an alveolar region other than the specified alveolar region. Can be displayed in different display modes.
  • the display control unit can display the curved surface by associating the grayscale or color template with the pixel values on the curved surface.
  • the display control unit can display the curved surface by setting the opacity with respect to the pixel value on the curved surface to 100%.
  • the display control unit displays the curved surface by associating a gray scale or a color template with an average value, a minimum value, or a maximum value of a plurality of pixel values within a predetermined thickness range from the curved surface. I can do it.
  • the curved surface can be displayed by setting the opacity with respect to the average value, minimum value, or maximum value to 100% by the display control unit.
  • a lung field region extraction unit that extracts a lung field region from a three-dimensional medical image may be provided, and the curved surface setting unit may set the curved surface with reference to the lung field region.
  • a thoracic region extracting unit that extracts a thoracic region from a three-dimensional medical image can be provided, and the curved surface setting unit can set the curved surface with the thoracic region as a reference.
  • the curved surface setting unit can receive an instruction to change the position of the curved surface by receiving an input from an input device.
  • the operation method of the three-dimensional medical image display control device is an operation method of the three-dimensional medical image display device including a curved surface setting unit and a display control unit, wherein the curved surface setting unit includes an area including a lung field.
  • the display control unit uses the pixel value on the curved surface set by the curved surface setting unit to convert the curved surface into a three-dimensional image. It is characterized by being displayed.
  • a curved surface setting unit that sets a curved surface along the surface of the lung field inside the lung field; Using the pixel value on the curved surface set by the curved surface setting unit, it functions as a display control unit that displays the curved surface in three dimensions.
  • the surface of the lung field is located inside the lung field.
  • a curved surface along the surface of the lung field is displayed using the pixel values on the set curved surface, so that there is an interstitial pneumonia lesion occurring inside the lung field surface. It is also possible to appropriately grasp the three-dimensional shape and the distribution state of interstitial pneumonia lesions throughout the lung.
  • the block diagram which shows schematic structure of the medical image diagnosis assistance system using one Embodiment of the three-dimensional medical image display control apparatus of this invention.
  • the figure which shows an example which displayed in three dimensions the pixel value on an equidistant surface The figure which shows an example which showed equidistant surface on a tomographic image
  • region with the three-dimensional display of an equidistant surface The block diagram which shows the modification of the medical image diagnosis assistance system shown in FIG.
  • FIG. 1 is a block diagram showing a schematic configuration of the medical image diagnosis support system of the present embodiment.
  • the medical image diagnosis support system of this embodiment includes a medical image display device 1, a three-dimensional medical image storage server 2, a display 3, and an input device 4, as shown in FIG.
  • the medical image display apparatus 1 is configured by installing the three-dimensional medical image display control program of this embodiment in a computer.
  • the medical image display device 1 includes a central processing unit (CPU) and a semiconductor memory, and a storage device such as a hard disk or an SSD (Solid State Drive) in which the above-described 3D medical image display control program is installed.
  • the three-dimensional medical image acquisition unit 10, the lung field region extraction unit 11, the curved surface setting unit 12, and the display control unit 13 as illustrated in FIG. 1 are configured by hardware. Then, the respective units operate by the central processing unit executing the three-dimensional medical image display control program installed on the hard disk.
  • the three-dimensional medical image acquisition unit 10 acquires a three-dimensional medical image 5 of the chest imaged in advance before the examination.
  • the three-dimensional medical image 5 for example, volume data reconstructed from slice data output from a CT apparatus or an MRI (Magnetic Resonance Imaging) apparatus, an MS (Multi-Slice) CT apparatus, or a cone beam CT apparatus is output. There is volume data.
  • the 3D medical image 5 is stored in advance in the 3D medical image storage server 2 together with the identification information of the subject, and the 3D medical image acquisition unit 10 receives the identification information of the subject input in the input device 4. Is read out from the three-dimensional medical image storage server 2.
  • the 3D medical image storage server 2 also stores past 3D medical images 5 of the same subject. The 3D medical image 5 taken this time, the 3D medical image 5 taken in the past, Both are read out and can be comparatively interpreted.
  • the lung field region extracting unit 11 extracts a lung field region from the three-dimensional medical image 5 of the chest acquired by the three-dimensional medical image acquiring unit 10. Specifically, the lung field region extraction unit 11 of the present embodiment first extracts the graph structure of the bronchial region included in the three-dimensional medical image 5 of the chest as a bronchial structure, and based on the extracted bronchial structure The field area is extracted.
  • a graph structure extraction method will be described.
  • the bronchi included in the three-dimensional medical image 5 appears as a region showing a low pixel value because the pixels inside the bronchi correspond to the air region, but the bronchial wall is a cylindrical or linear structure showing a relatively high pixel value. It is thought that. Therefore, a bronchus is extracted by performing a structural analysis of the shape based on the distribution of pixel values for each pixel.
  • bronchus branches in multiple stages, and the diameter of the bronchus decreases as the end approaches.
  • a Gaussian pyramid image obtained by multi-resolution conversion of a three-dimensional medical image, that is, a plurality of three-dimensional medical images having different resolutions, is generated. Different size tubular structures are detected by scanning a detection algorithm for each image of the Gaussian pyramid.
  • the Hessian matrix of each pixel of the three-dimensional image of each resolution is calculated, and it is determined whether the pixel is in the tubular structure from the magnitude relationship of the eigenvalues of the Hessian matrix.
  • the Hessian matrix is a matrix whose elements are second-order partial differential coefficients of density values in the directions of the respective axes (the x-axis, y-axis, and z-axis of the three-dimensional image), and is a 3 ⁇ 3 matrix as shown in the following equation. .
  • the eigenvalues of the Hessian matrix in an arbitrary pixel are ⁇ 1, ⁇ 2, and ⁇ 3, when two eigenvalues are large and one eigenvalue is close to 0, for example, when ⁇ 3, ⁇ 2 >> ⁇ 1, ⁇ 1 ⁇ 0 is satisfied
  • the pixel is known to be a tubular structure.
  • the eigenvector corresponding to the minimum eigenvalue ( ⁇ 1 ⁇ 0) of the Hessian matrix coincides with the principal axis direction of the tubular structure.
  • the tubular structure extracted in this way is not always detected as one graph structure in which all the tubular structures are connected. Therefore, after the determination of the entire three-dimensional medical image is completed, the detected tubular structure is within a certain distance, and the direction of the basic line connecting any point on the two extracted tubular structures and each tubular structure Connection relationship of tubular structures extracted by determining whether or not a plurality of tubular structures are connected by evaluating whether the angle formed by the main axis direction of the structure is within a certain angle To rebuild. By this reconstruction, extraction of the bronchial graph structure is completed (for details, refer to JP2010-220A). 742 etc.).
  • the extracted graph structure is classified into a start point, an end point, a branch point, and a side, and a graph structure representing the bronchi can be obtained by connecting the start point, the end point, and the branch point with the sides.
  • the lung field region extraction unit 11 performs anatomical labeling on the graph structure based on the branch of the graph structure.
  • This labeling is a process performed based on the control region of the bronchus.
  • Each control region of the bronchus divided by this labeling is the lung lobe region (upper right lobe region, right middle lobe region, right lower lobe region, upper left lobe). Leaf region and lower left leaf region). Since bronchial labeling is a well-known technique, a detailed description thereof will be omitted.
  • the lung field region extraction unit 11 acquires a lung lobe region (upper right lobe region, right middle lobe region, right lower lobe region, left upper lobe region, and left lower lobe region) based on the bronchial labeling result.
  • Each lobe region can be obtained by applying a Quickhull algorithm to a bronchial graph structure (divided bronchial structure) with the same label, for example, to obtain a convex region.
  • a bronchial graph structure divided bronchial structure
  • the lung field region extraction method is not limited to the above method, and other known methods may be used.
  • the curved surface setting unit 12 sets a curved surface along the surface of the lung field region inside the lung field region extracted as described above.
  • the “curved surface along the surface of the lung field” does not necessarily mean that the entire curved surface is equidistant from the surface of the lung field, and the distance from a part of the lung field of the curved surface is different from the surroundings.
  • the curved surface may be a curved surface substantially along the lung field surface when the entire curved surface is viewed.
  • a plane equidistant from the surface of the lung field region is set.
  • a method of setting a plane equidistant from the surface of the lung field region for example, by performing a distance conversion process on the three-dimensional medical image of the lung field region extracted as described above, A distance conversion image representing the distance of the distance may be generated, and a plane equidistant from the surface of the lung field region may be specified and set based on the distance conversion image. Since the distance conversion process is a known process, detailed description thereof is omitted here. Further, the method for setting the equidistant surface is not limited to the method using the distance conversion image, and other known methods such as a method using a mask may be used.
  • the distance of the equidistant surface from the surface of the lung field is preferably such that the equidistant surface includes a cross section of the alveoli so that the progress of interstitial pneumonia can be observed, for example.
  • the distance is preferably 5 mm to 10 mm.
  • the present invention is not limited to this, and the distance from the lung field surface of the equidistant surface may be changed by the user performing a predetermined input operation using the input device 4, for example.
  • the distance may be sequentially changed in increments of 1 mm from the lung field surface in accordance with the wheel operation of the user's mouse device, or a plurality of preset different distances, for example, 5 mm, 10 mm, and 15 mm.
  • Etc. may be sequentially switched according to the user's operation input. Or you may make it change automatically the distance from the lung field surface of an equidistant surface irrespective of the user's operation input as mentioned above.
  • the range for setting the equidistant surface is not necessarily the entire lung field region, but may be a partial range.
  • interstitial pneumonia often progresses from the back side of the lung field region, and therefore, along the surface of a part of the entire lung field region including at least the surface of the lung field region on the back side, etc.
  • a distance plane may be set.
  • the display control unit 13 specifies the pixel value on the curved surface set by the curved surface setting unit 12, and uses the pixel value to display the curved surface in a three-dimensional manner on the display 3. Specifically, the display control unit 13 of the present embodiment specifies the pixel value on the curved surface from the three-dimensional medical image 5 acquired by the three-dimensional medical image acquisition unit 10, and grayscales the pixel value. And the curved surface is displayed with the opacity of 100% set for the pixel value on the curved surface.
  • a gray scale from white to black is assigned to the CT value of window level ⁇ 600 and window width 1600.
  • This CT value is a value of each pixel of slice data reconstructed from a plurality of radiographic images taken by the CT apparatus. That is, in the present embodiment, a gray scale is directly assigned to each pixel value of slice data reconstructed from a plurality of radiation images.
  • the gray scale is used.
  • the present invention is not limited to this.
  • a color template in which the brightness or saturation of a predetermined color such as green is changed may be used.
  • the opacity is set to 100%.
  • the present invention is not limited to this, and may be set to 10% to 50%.
  • the identification information of the subject is input by the user using the input device 4, and the three-dimensional medical image acquisition unit 10 converts the three-dimensional medical image 5 of the chest corresponding to the input identification information into the three-dimensional medical image. It is read from the storage server 2 and acquired (S10).
  • the three-dimensional medical image 5 of the chest acquired by the three-dimensional medical image acquisition unit 10 is output to the lung field region extraction unit 11, and the lung field region extraction unit 11 uses the three-dimensional medical image 5 input by the method described above.
  • a lung field region is generated (S12).
  • the lung field region extracted by the lung field region extraction unit 11 is input to the curved surface setting unit 12, and the curved surface setting unit 12 sets an equidistant surface along the surface of the lung field region inside the lung field region (S14). ).
  • the display control unit 13 specifies a pixel value on the equidistant surface set by the curved surface setting unit 12 among the pixel values of the three-dimensional medical image of the lung field region, and the gray scale described above for the pixel value. And an opacity of 100% are set, thereby causing the equidistant surface to be displayed three-dimensionally on the display 3 (S16).
  • FIG. 4 shows an example of the three-dimensional display of the equidistant surface displayed as described above
  • FIG. 5 shows the position of the equidistant surface (indicated by a white solid line in FIG. 5). This is shown on a tomographic image.
  • an image showing the equidistant surface on the tomographic image as shown in FIG. 5 may be displayed.
  • the equidistant surface is displayed three-dimensionally as described above, but an image representing the lung field region may be displayed together with the equidistant surface.
  • an image representing the lung field region for example, an opacity of 10% to 50% is assigned to a pixel value on the surface of the lung field region and displayed, and an equidistant surface is displayed via an image representing the lung field region. What is necessary is just to make it display so that a three-dimensional display can be seen through. For example, as shown in FIG. 6, an image of a line representing the outline of the lung field region (shown by a white solid line in FIG. 6) may be displayed.
  • an image representing the alveolar region may be further displayed.
  • an alveolar region extracting unit 14 that further extracts an alveolar region from the three-dimensional medical image 5, and a plurality of alveolar regions extracted by the alveolar region extracting unit 14 Among them, a display alveolar region specifying unit 15 for specifying an alveolar region having a size equal to or larger than a preset threshold as a display target alveolar region is provided, and the display control unit 13 extracts the alveolar region by the alveolar region extracting unit 14. Of the plurality of alveolar regions, only the image representing the alveolar region to be displayed specified by the display alveolar region specifying unit 15 may be displayed.
  • the alveolar region extraction processing in the alveolar region extraction unit 14 a known method can be used.
  • the alveolar wall is inflamed, so-called fibrosis, and shows a higher CT value than the air region.
  • the alveolar region may be extracted by extracting the alveolar wall.
  • a pixel whose CT value is within a predetermined range ( ⁇ 900 to ⁇ 750) is extracted using a threshold processing method. By doing so, the alveolar region may be extracted.
  • the alveolar region may be extracted by removing the blood vessel region or the bronchial region from the lung field region.
  • Various known methods may be used for extracting the lung field region, the blood vessel region, and the bronchial region.
  • the display alveolar region specifying unit 15 specifies an alveolar region having a size greater than or equal to a threshold value from among a plurality of alveolar regions as a display target alveolar region. This is because, in the case of interstitial pneumonia, the alveoli become larger as it gets worse, and the degree of deterioration of interstitial pneumonia can be diagnosed by displaying the enlarged alveolar region.
  • the display alveolar region specifying unit 15 calculates, for example, the volume of each alveolar region, and specifies the alveolar region whose volume is equal to or greater than a preset threshold as a display target alveolar region. Instead of calculating the volume of each alveolar region, for example, the maximum cross-sectional area or maximum length of each alveolar region may be calculated to determine the threshold value.
  • the volume of the alveolar region extracted from the 3D medical image 5 taken in the past may be used as the threshold.
  • the alveolar region extracted from the 3D medical image 5 taken in the past and the alveolar region extracted from the 3D medical image 5 taken this time are associated with each other by alignment or the like. It is desirable to set the size of the past alveolar region to be a threshold value.
  • the display control unit 13 displays only an image representing the alveolar region to be displayed among the plurality of alveolar regions as described above. These may be displayed simultaneously with the three-dimensional display, or may be displayed separately. When displaying an image representing the alveolar region simultaneously with the three-dimensional display of the equidistant surface, for example, an opacity of 10% to 50% is assigned to the pixel value on the equidistant surface and displayed. What is necessary is just to make it display so that the image showing the alveolar region can be seen through through the three-dimensional display.
  • the display control unit 13 may display the alveolar region specifying unit 15.
  • the identified alveolar region may be displayed in a display mode different from the alveolar region other than the identified alveolar region.
  • an image representing the alveolar region identified by the display alveolar region identifying unit 15 is displayed in color, and an image representing the alveolar region other than the identified alveolar region is displayed in monochrome. Also good.
  • the opacity of the image representing the alveolar region specified by the display alveolar region specifying unit 15 is relatively increased, and the opacity of the image representing the alveolar region other than the specified alveolar region is relatively set. May be displayed lower.
  • these displays may be switched in accordance with, for example, an operation input on the input device 4 by the user.
  • three-dimensional display can be performed by assigning 100% opacity to the pixel values on the equidistant surface.
  • the user may switch between display and non-display of the three-dimensional display in which 100% opacity is assigned to the pixel value on the equidistant surface.
  • the grayscale or color template is associated with the pixel value on the equidistant plane and displayed three-dimensionally.
  • the display control unit 13 receives thickness information from an equidistant surface, and a grayscale or color template for an average value, a minimum value, or a maximum value of a plurality of pixel values within the thickness range. May be associated with each other to perform three-dimensional display of equidistant surfaces. For example, when the three-dimensional display is performed using the average value, it is possible to grasp the state of the lesion over a wider range in the thickness direction as compared with the case where only the pixel values on the equidistant surface are used.
  • the fiberization (high CT value) distribution information within the thickness information range can be grasped, and the three-dimensional display is performed using the minimum value.
  • emphysema (low CT value) distribution information within the thickness information range can be grasped.
  • the thickness information is information on the thickness range from the equidistant surface to the inner side of the lung field.
  • the thickness information may be set in advance in the apparatus, or may be set and input by the user using the input device 4. Further, selection of any one of an average value, a minimum value, and a maximum value of a plurality of pixel values may be accepted.
  • the average value, minimum value, or maximum value is also set. It is desirable to assign a opacity of 100% to the three-dimensional display.
  • the present invention is not limited to this, and may be set to 10% to 50% as described above.
  • the lung field region extraction unit 11 extracts the lung field region, and the equidistant surface is set on the basis of the surface of the extracted lung field region.
  • a thoracic region extraction unit 16 is provided instead of the lung field extraction unit 11, and the thoracic region extraction unit 16 performs three-dimensional medical use.
  • a ribcage region may be extracted from the image 5 and an equidistant surface may be set with the extracted ribcage region as a reference.
  • the distance relationship between the surface of the thorax area and the surface of the lung field area is a medically determined relationship, so consider the distance relationship and equidistant inside the lung field area based on the surface of the rib cage area
  • the surface may be set.
  • a method for extracting a rib cage region a known method can be used.
  • a method of extracting a rib and determining a rib cage region from the rib region can be used.
  • “Zhou Hui, et al.,“ Extraction of thorax from 3D torso multi-slice CT images using chest skeletal information ” IEICE paper D-II Vol. J88-D-II No.9 pp. 1999-2002 "can be used.
  • the diaphragm region may be extracted, and the thoracic region may be extracted from the rib region and the diaphragm region.
  • a known method can be used. For example, “Takashi Takahashi, et al.“ Development of diaphragm extraction processing for abdominal extraction from non-contrast 3D abdominal CT image ”, P14 / JAMIT2009JAProceedings Can be used.
  • both the 3D medical image 5 taken this time and the 3D medical image 5 taken in the past are read, and these 3D medical images are read out.
  • the pixel values on the equidistant surface described above are acquired for each of the images 5, and the three-dimensional display of the pixel values on the equidistant surface of the three-dimensional medical image 5 photographed this time and the three-dimensional medical image 5 photographed in the past are obtained. It is also possible to simultaneously perform three-dimensional display of pixel values on the equidistant surface so that they can be comparatively interpreted.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pulmonology (AREA)
  • Physiology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Human Computer Interaction (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

間質性肺炎の病変の有無やその進行度を適切に把握可能な表示を行うことができる３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラムを提供する。 肺野を含む領域を撮影した３次元医用画像において、肺野の内側に肺野の表面に沿った曲面を設定する曲面設定部１２と、曲面設定部１２によって設定された曲面上の画素値を用いてその曲面を立体的に表示させる表示制御部１３とを備える。

Description

３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム
　本発明は、肺の３次元医用画像を表示する３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラムに関するものであり、特に、肺野の内側に設定された曲面の３次元表示に関するものである。
　間質性肺炎は瀰漫性肺疾患のひとつであり、多くの場合、肺の大きな範囲に渡って病変が分布している。したがって、間質性肺炎の進行度を診断するためには、肺全体に広がる病変を把握する必要がある。
　一方、近年、ＣＴ（Computed Tomography）やＭＲＩ（magnetic resonance imaging）
などによって撮影された３次元医用画像を用いた種々の画像診断が行われている。たとえば、特許文献１においては、肺結節を抽出し、その肺結節に連続する周辺領域を特定し、ボリュームレンダリングを行うことによって肺結節およびその周辺領域を表示する方法が提案されている。
　また、特許文献２においては、肺野領域の３次元医用画像を取得し、これを３次元表示する方法が提案されている。
特開２００８－１２２９１号公報 国際公開第２００７/０９４４１２
　ここで、間質性肺炎の画像診断を行う場合、ＣＴ画像は、間質性肺炎の病変を明瞭に捉えることができため、診断に有効と考えられる。
　しかしながら、間質性肺炎の病変は複数の断層画像に存在しているため、通常の２次元断層画像の読影では、間質性肺炎の病変の全体的な広がりを把握することができない。間質性肺炎の病変の有無を診断するためには、３次元な形状を把握することが重要であるとともに、間質性肺炎の進行度を診断するために、間質性肺炎の病変の肺全体における分布を把握することが重要である。
　また、間質性肺炎の病変は、軽症から重症まで見え方が大きく異なり、さらに上述したように広範囲に分布しているため、間質性肺炎の病変領域の抽出は困難である。また、間質性肺炎の病変は、濃度の微小の変化や、蜂窩状のような薄い壁に囲まれた小さな空洞のパターンであるため、たとえば、特許文献１に記載されているような一般的なボリュームレンダリングによって間質性肺炎の病変を表示した場合、他の細い血管や気管支なども同時に表示されるため、有用な診断情報を得ることができない。
　また、特許文献２に記載のように単純に肺野領域を３次元表示しただけでは、肺野領域の表面しか観察することができず、肺野領域の内側において発生している間質性肺炎の病変の状態を観察することができない。
　本発明は、上記事情に鑑み、間質性肺炎の病変の有無やその進行度を適切に把握可能な表示を行うことができる３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラムを提供することを目的とするものである。
　本発明の３次元医用画像表示制御装置は、肺野を含む領域を撮影した３次元医用画像において、肺野の内側に肺野の表面に沿った曲面を設定する曲面設定部と、曲面設定部によって設定された曲面上の画素値を用いてその曲面を立体的に表示させる表示制御部とを備えたことを特徴とする。
　また、上記本発明の３次元医用画像表示制御装置において、上記曲面を、肺野の表面から等距離の面とすることができる。
　また、曲面設定部を、肺野の表面からの上記曲面の位置を変更することにより設定するものとできる。
　また、曲面設定部を、肺野の表面からの距離が異なる複数の上記曲面を切り替えることにより設定するものとできる。
　また、表示制御部を、肺野の輪郭線または表面を表示させるものとできる。
　また、３次元医用画像から肺胞領域を抽出する肺胞領域抽出部と、肺胞領域抽出部によって抽出された複数の肺胞領域のうち閾値以上の大きさの肺胞領域を表示対象の肺胞領域として特定する表示肺胞領域特定部とを設け、表示制御部を、肺胞領域抽出部によって抽出された複数の肺胞領域のうち表示対象の肺胞領域のみを表示させるものとできる。
　また、３次元医用画像から肺胞領域を抽出する肺胞領域抽出部と、肺胞領域抽出部によって抽出された複数の肺胞領域のうち閾値以上の大きさの肺胞領域を特定する表示肺胞領域特定部とを設け、表示制御部を、肺胞領域抽出部によって抽出された複数の肺胞領域のうち上記特定された肺胞領域をその特定された肺胞領域以外の肺胞領域とは異なる表示態様により表示させるものとできる。
　また、表示制御部を、曲面上の画素値に対してグレースケールまたはカラーテンプレートを対応づけて上記曲面を表示させるものとできる。
　また、表示制御部を、曲面上の画素値に対する不透明度を１００％に設定することにより上記曲面を表示させるものとできる。
　また、表示制御部を、上記曲面から予め設定された厚さの範囲内の複数の画素値の平均値、最小値または最大値に対してグレースケールまたはカラーテンプレートを対応づけて上記曲面を表示させるものとできる。
　また、表示制御部を、上記平均値、最小値または最大値に対する不透明度を１００％に設定することにより上記曲面を表示させるものとできる。
　また、３次元医用画像から肺野領域を抽出する肺野領域抽出部を設け、曲面設定部を、肺野領域を基準として上記曲面を設定するものとできる。
　また、３次元医用画像から胸郭領域を抽出する胸郭領域抽出部を設け、曲面設定部を、胸郭領域を基準として上記曲面を設定するものとできる。
　また、曲面設定部を、入力装置からの入力を受け付けることによって上記曲面の位置の変更指示を受け付けるものとできる。
　本発明の３次元医用画像表示制御装置の作動方法は、曲面設定部と、表示制御部とを備えた３次元医用画像表示装置の作動方法であって、曲面設定部が、肺野を含む領域を撮影した３次元医用画像において、肺野の内側に肺野の表面に沿った曲面を設定し、表示制御部が、曲面設定部によって設定された曲面上の画素値を用いてその曲面を立体的に表示させることを特徴とする。
　本発明の３次元医用画像表示制御プログラムは、コンピュータを、肺野を含む領域を撮影した３次元医用画像において、肺野の内側に肺野の表面に沿った曲面を設定する曲面設定部と、曲面設定部によって設定された曲面上の画素値を用いてその曲面を立体的に表示させる表示制御部として機能させることを特徴とする。
　本発明の３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラムによれば、肺野を含む領域を撮影した３次元医用画像において、肺野の内側に肺野の表面に沿った曲面を設定し、その設定した曲面上の画素値を用いてその曲面を立体的に表示させるようにしたので、肺野の表面より内側において発生している間質性肺炎の病変の有無およびその３次元的形状と、肺全体に亘る間質性肺炎の病変の分布状態を適切に把握することができる。
本発明の３次元医用画像表示制御装置の一実施形態を用いた医用画像診断支援システムの概略構成を示すブロック図 曲面設定部において設定された曲面上の画素値に対してグレースケールを対応づけた場合を説明するための図 本発明の３次元医用画像表示制御装置の一実施形態を用いた医用画像診断支援システムの作用を説明するためのフローチャート 等距離面上の画素値を立体表示した一例を示す図 断層画像上に等距離面を示した一例を示す図 等距離面の３次元表示とともに、肺野領域の輪郭を表す線の画像を表示した一例を示す図 図１に示す医用画像診断支援システムの変形例を示すブロック図 図１に示す医用画像診断支援システムの変形例を示すブロック図
　以下、本発明の３次元医用画像表示制御装置および方法並びにプログラムの一実施形態を用いた医用画像診断支援システムについて、図面を参照しながら詳細に説明する。図１は、本実施形態の医用画像診断支援システムの概略構成を示すブロック図である。
　本実施形態の医用画像診断支援システムは、図１に示すように、医用画像表示装置１と、３次元医用画像保管サーバ２と、ディスプレイ３と、入力装置４とを備えている。
　医用画像表示装置１は、コンピュータに本実施形態の３次元医用画像表示制御プログラムをインストールすることによって構成されたものである。医用画像表示装置１は、中央処理装置（ＣＰＵ）および半導体メモリや、上述した３次元医用画像表示制御プログラムがインストールされたハードディスクやＳＳＤ(Solid State Drive)等のストレージデバイスを備えており、これらのハードウェアによって、図１に示すような３次元医用画像取得部１０、肺野領域抽出部１１、曲面設定部１２および表示制御部１３が構成されている。そして、ハードディスクにインストールされた３次元医用画像表示制御プログラムが中央処理装置によって実行されることによって上記各部がそれぞれ動作する。
　３次元医用画像取得部１０は、検査前などに予め撮影された胸部の３次元医用画像５を取得するものである。３次元医用画像５としては、たとえばＣＴ装置やＭＲＩ（MagneticResonance Imaging）装置などから出力されたスライスデータから再構成されたボリュームデータや、ＭＳ（Multi Slice）ＣＴ装置やコーンビームＣＴ装置から出力されたボリュームデータなどがある。
　３次元医用画像５は、３次元医用画像保管サーバ２に被検者の識別情報とともに予め保管されており、３次元医用画像取得部１０は、入力装置４において入力された被検者の識別情報に対応する３次元医用画像５を３次元医用画像保管サーバ２から読み出すものである。また、３次元医用画像保管サーバ２には、同一の被検者の過去の３次元医用画像５も保管されており、今回撮影した３次元医用画像５と過去に撮影した３次元医用画像５との両方が読み出されて比較読影可能なように構成されている。 
　肺野領域抽出部１１は、３次元医用画像取得部１０によって取得された胸部の３次元医用画像５から肺野領域を抽出するものである。具体的には、本実施形態の肺野領域抽出部１１は、まず、胸部の３次元医用画像５に含まれる気管支領域のグラフ構造を気管支構造として抽出し、その抽出した気管支構造に基づいて肺野領域を抽出するものである。以下、グラフ構造の抽出方法の一例を説明する。
　３次元医用画像５に含まれる気管支は、気管支の内部の画素は空気領域に相当するため低い画素値を示す領域として表れるが、気管支壁は比較的高い画素値を示す円柱あるいは線状の構造物であると考えられる。そこで、各画素ごとに画素値の分布に基づく形状の構造解析を行なって気管支を抽出する。
　気管支は多段階に分岐し末端に近づくほど気管支の径は小さくなっていく。異なるサイズの気管支（管状構造物）を検出することができるように、予め、三次元医用画像を多重解像度変換したガウシアンピラミッド画像、つまり異なる解像度の複数の三次元医用画像を生成し、その生成したガウシアンピラミッドの各画像ごとに検出アルゴリズムを走査することにより異なるサイズの管状構造物を検出する。
　まず、各解像度の三次元画像の各画素のヘッセ行列を算出し、ヘッセ行列の固有値の大小関係から管状構造物内の画素であるかを判定する。ヘッセ行列は、各軸（三次元画像のｘ軸、ｙ軸、ｚ軸）方向における濃度値の２階の偏微分係数を要素とする行列であり、下式のように３×３行列となる。
Figure JPOXMLDOC01-appb-M000001
 　任意の画素におけるヘッセ行列の固有値をλ１、λ２、λ３としたとき、固有値のうち２つの固有値が大きく、１つの固有値が０に近い場合、例えば、λ３、λ２≫λ１、λ１≒０を満たすとき、その画素は管状構造物であることが知られている。また、ヘッセ行列の最小の固有値（λ１≒０）に対応する固有ベクトルが管状構造物の主軸方向に一致する。
　気管支はグラフ構造によって表すことができるが、このようにして抽出された管状構造物は、全ての管状構造物が繋がった１つのグラフ構造として検出されるとは限らない。そこで、三次元医用画像全体の判別が終了後、検出された管状構造物が一定の距離内にあり、かつ抽出された二つの管状構造物上の任意の点を結ぶ基本線の向きと各管状構造物の主軸方向とが成す角が一定角度以内であるかについて評価することにより、複数の管状構造物が接続されるものであるか否かを判定して抽出された管状構造物の接続関係を再構築する。この再構築により、気管支のグラフ構造の抽出が完了する（詳細は、特開２０１０-２２０
７４２号公報などを参照）。
　そして、抽出したグラフ構造を、開始点、端点、分岐点及び辺に分類し、開始点、端点及び分岐点を辺で連結することによって、気管支を表すグラフ構造を得ることができる。
　次に、肺野領域抽出部１１は、グラフ構造の分岐に基づいて、このグラフ構造に対して解剖学的なラベリングを行う。このラベリングは、気管支の支配領域に基づいて行われる処理であり、このラベリングによって区分けされた気管支の各支配領域は、それぞれ肺における肺葉領域（右上葉領域、右中葉領域、右下葉領域、左上葉領域および左下葉領域）に対応するものとなる。なお、気管支ラベリングについては、既に公知な技術であるため詳細な説明を省略するが、たとえば「Pechin Lo et al. “A Bottom-up approach for labeling of human airway trees”, Fourth INTERNATIONAL WORKSHOP ON PULMONARY IMAGE　ANALYSIS, P23-P34」や、「Jaesung LEE et al. “Bronchial Segment Matching in Low-dose Lung CT Scan Pairs”, Medical Imaging 2009: Computer-Aided Diagnosis, Proceeding of SPIE Vol. 7260, 72600A1-72600A8」などの技術を用いることができる。
　そして、肺野領域抽出部１１は、気管支ラベリング結果に基づいて、肺葉領域（右上葉領域、右中葉領域、右下葉領域、左上葉領域および左下葉領域）を取得するものである。各肺葉領域は、たとえば同じラベルの付された気管支のグラフ構造（分割気管支構造）に対してQuickhullアルゴリズムを適用して凸領域を取得することによって取得することができる。このように各肺葉領域を抽出することによって肺野全体の領域が抽出される。なお、肺野領域の抽出方法としては、上記のような方法に限らず、その他の公知な方法を用いるようにしてもよい。たとえば「北坂 孝幸,他、“可変ベジエ曲面による形状モデルを用いた３次元胸部Ｘ線ＣＴ像からの肺野領域抽出”電子情報通信学会論文紙 D-II Vol. J83-D-II No.1 pp.165-174」に記載の方法を用いることができる。
　曲面設定部１２は、上述したようにして抽出された肺野領域の内側に、肺野領域の表面に沿った曲面を設定するものである。
　ここで、上記「肺野の表面に沿った曲面」とは、必ずしも曲面全体が肺野の表面から等距離でなくてもよく、曲面の一部の肺野表面からの距離が周辺と異なっていてもよく、曲面全体をみたときに、実質的に肺野表面に沿った曲面であればよい。
　ただし、本実施形態においては、特に、肺野領域の表面から等距離の面を設定する場合について説明する。肺野領域の表面から等距離の面を設定する方法としては、たとえば、上述したようにして抽出された肺野領域の３次元医用画像に対して距離変換処理を施すことによって肺野領域表面からの距離を表す距離変換画像を生成し、その距離変換画像に基づいて肺野領域の表面から等距離の面を特定して設定するようにすればよい。なお、距離変換処理については公知な処理であるので、ここでは詳細な説明は省略する。また、等距離の面の設定方法としては距離変換画像を用いた方法に限らず、たとえばマスクを利用した方法などその他の公知な方法を用いるようにしてもよい。
　等距離面の肺野表面からの距離については、たとえば間質性肺炎の進行度が観察できるように、等距離面が肺胞の断面を含むような距離であることが望ましく、肺野の表面から５ｍｍ～１０ｍｍの距離であることが望ましい。
　また、これに限らず、たとえば、ユーザが入力装置４を用いて所定の入力操作を行うことによって、等距離面の肺野表面からの距離を変更できるようにしてもよい。この場合、たとえばユーザのマウス装置のホイール操作に応じて肺野表面から１ｍｍ刻みなどで細かく距離を順次変更できるようにしてもよいし、予め設定された複数の異なる距離、たとえば５ｍｍと１０ｍｍと１５ｍｍなどをユーザの操作入力に応じて順次切り替えるようにしてもよい。または、上述したようなユーザの操作入力によらずに、等距離面の肺野表面からの距離を自動的に変更するようにしてもよい。
　また、等距離面を設定する範囲としては、必ずしも肺野領域全体でなくてもよく、一部の範囲でもよい。たとえば、間質性肺炎は、肺野領域の背中側から進行することが多いので、肺野領域全体のうち、少なくとも背中側の肺野領域の表面を含む一部の領域の表面に沿った等距離面を設定するようにしてもよい。また、肺野領域全体のうち、横隔膜側の肺野領域の表面を除く一部の領域の表面に沿った等距離面を設定するようにしてもよい。
　表示制御部１３は、曲面設定部１２において設定された曲面上の画素値を特定し、その画素値を用いて上記曲面をディスプレイ３において立体的に表示させるものである。具体的には、本実施形態の表示制御部１３は、３次元医用画像取得部１０によって取得された３次元医用画像５から上記曲面上の画素値を特定し、その画素値に対してグレースケールを対応づけ、かつ曲面上の画素値に対して不透明度１００％を設定して上記曲面を表示させる。
　胸部の３次元医用画像５がＣＴ装置によって撮影されたものである場合、たとえば図２に示すようにウィンドウレベル－６００、ウィンドウ幅１６００のＣＴ値に対して白から黒までのグレースケールが割り当てられる。このＣＴ値は、ＣＴ装置によって撮影された複数の放射線画像から再構成したスライスデータの各画素の値である。すなわち、本実施形態においては、複数の放射線画像から再構成したスライスデータの各画素値に対して直接グレースケールが割り当てられる。
　ただし、本実施形態においては、グレースケールを用いるようにしたが、これに限らず、たとえば緑などの所定の色の明度や彩度などを変化させたカラーテンプレートを用いるようにしてもよい。また、本実施形態においては、不透明度を１００％に設定するようにしたが、これに限らず、１０％～５０％に設定してもよい。
　次に、本実施形態の医用画像診断支援システムの作用について、図３に示すフローチャートを参照しながら説明する。
　まず、ユーザによって入力装置４を用いて被検者の識別情報が入力され、３次元医用画像取得部１０は、その入力された識別情報に対応する胸部の３次元医用画像５を３次元医用画像保管サーバ２から読み出して取得する（Ｓ１０）。
　３次元医用画像取得部１０によって取得された胸部の３次元医用画像５は肺野領域抽出部１１に出力され、肺野領域抽出部１１は、上述した方法によって入力された３次元医用画像５から肺野領域を生成する（Ｓ１２）。
　肺野領域抽出部１１によって抽出された肺野領域は曲面設定部１２に入力され、曲面設定部１２は、肺野領域の内側に肺野領域の表面に沿った等距離面を設定する（Ｓ１４）。
　そして、表示制御部１３は、肺野領域の３次元医用画像の画素値のうち曲面設定部１２において設定された等距離面上の画素値を特定し、その画素値に対して上述したグレースケールを対応づけるとともに、不透明度１００％を設定し、これにより等距離面をディスプレイ３において立体的に表示させる（Ｓ１６）。図４は、上述したようにして表示された等距離面の３次元表示の一例を示すものであり、図５は、等距離面の位置（図５においては白い実線によって示されている）を断層画像上に示したものである。なお、図４に示すような等距離面の３次元表示だけでなく、図５に示すような断層画像上に等距離面を示したものも表示するようにしてもよい。
　また、上記実施形態の医用画像診断支援システムにおいては、上述したように等距離面を３次元表示するようにしたが、さらにその等距離面とともに肺野領域を表す画像を表示するようにしてもよい。肺野領域を表す画像の表示方法としては、たとえば肺野領域の表面の画素値に対して不透明度１０％～５０％を割り当てて表示し、肺野領域を表す画像を介して等距離面の３次元表示が透けて見えるような表示とすればよい。また、これに限らず、たとえば、図６に示すように肺野領域の輪郭を表す線の画像（図６では白い実線で示されている）を表示するようにしてもよい。このように等距離面の３次元表示とともに、肺野領域を表す画像を表示することによって、肺野領域の表面から等距離面までの距離を視覚的に確認することができ、間質性肺炎の診断の精度を向上させることができる。
　また、上記実施形態の医用画像診断支援システムにおいて、さらに肺胞領域を表す画像を表示するようにしてもよい。具体的には、図７に示すように、さらに３次元医用画像５から肺胞領域を抽出する肺胞領域抽出部１４と、その肺胞領域抽出部１４によって抽出された複数の肺胞領域のうち、予め設定された閾値以上の大きさの肺胞領域を表示対象の肺胞領域として特定する表示肺胞領域特定部１５とを設け、表示制御部１３が、肺胞領域抽出部１４によって抽出された複数の肺胞領域のうち、表示肺胞領域特定部１５によって特定された表示対象の肺胞領域を表す画像のみを表示させるようにしてもよい。
　肺胞領域抽出部１４における肺胞領域抽出処理については、既に公知は方法を用いることができる。なお、たとえば間質性肺炎の場合、肺胞壁が炎症をおこして、いわゆる線維化して空気領域よりも高いＣＴ値を示すので、クラスター解析などを行うことによって空気領域よりも高いＣＴ値で構成される肺胞壁を抽出することにより肺胞領域を抽出するようにしてもよい。また、肺胞は上述したように空気領域よりも高く、血管などよりは低いＣＴ値を持つため、閾値処理手法を用いてＣＴ値が所定の範囲（－９００から－７５０）にある画素を抽出することによって肺胞領域を抽出するようにしてもよい。また、肺野領域から血管領域や気管支領域を取り除くことによって肺胞領域を抽出するようにしてもよい。肺野領域、血管領域および気管支領域の抽出方法については、種々の公知な手法を用いるようにすればよい。
　表示肺胞領域特定部１５は、上述したように複数の肺胞領域の中から閾値以上の大きさの肺胞領域を表示対象の肺胞領域として特定するものであるが、これは、たとえば間質性肺炎の場合、悪化するほど肺胞が大きくなるため、大きくなった肺胞領域を表示することによって間質性肺炎の悪化度が診断できるからである。表示肺胞領域特定部１５は、たとえば各肺胞領域の体積を算出し、この体積が予め設定された閾値以上の肺胞領域を表示対象の肺胞領域として特定するものである。なお、各肺胞領域の体積を算出するのではなく、たとえば各肺胞領域の最大断面積や最大長などを算出して閾値判定するようにしてもよい。
　また、上述した閾値判定を行う際に用いられる閾値については、過去に撮影された３次元医用画像５から抽出された肺胞領域の体積などを閾値としてもよい。この場合、過去に撮影された３次元医用画像５から抽出された肺胞領域と、今回撮影された３次元医用画像５から抽出された肺胞領域とで位置合わせなどによって対応付けを行い、対応する過去の肺胞領域の大きさを閾値とすることが望ましい。
　そして、表示制御部１３は、上述したように複数の肺胞領域のうち表示対象の肺胞領域を表す画像のみを表示させるものであるが、この肺胞領域を表す画像については、等距離面の３次元表示と同時に表示するようにしてもよいし、別々に表示してもよい。等距離面の３次元表示と同時に肺胞領域を表す画像を表示させる場合には、たとえば等距離面上の画素値に対して不透明度１０％～５０％を割り当てて表示し、等距離面の３次元表示を介して肺胞領域を表す画像が透けて見えるような表示とすればよい。
　また、上記説明では、肺胞領域のうち表示対象の肺胞領域を表す画像のみを表示させるようにしたが、これに限らず、たとえば、表示制御部１３が、表示肺胞領域特定部１５によって特定された肺胞領域をその特定された肺胞領域以外の肺胞領域とは異なる表示態様で表示するようにしてもよい。具体的には、たとえば表示肺胞領域特定部１５によって特定された肺胞領域を表す画像をカラー表示し、上記特定された肺胞領域以外の肺胞領域を表す画像をモノクロ表示するようにしてもよい。また、たとえば表示肺胞領域特定部１５によって特定された肺胞領域を表す画像の不透明度を相対的に高くし、上記特定された肺胞領域以外の肺胞領域を表す画像の不透明度を相対的に低くして表示するようにしてもよい。
　また、肺胞領域を表す画像の表示と等距離面の３次元表示とを別々に表示させる場合には、たとえばユーザによる入力装置４における操作入力に応じてこれらの表示を切り替えるようにすればよい。この場合、等距離面の画素値に対して不透明度１００％を割り当てて３次元表示を行うことができる。
　また、肺胞領域を表す画像を表示させた状態において、等距離面の画素値に対して不透明度１００％を割り当てた３次元表示の表示と非表示とをユーザが切り替えるようにしてもよい。
　また、上記実施形態の医用画像診断支援システムにおいては、等距離面上における画素値に対してグレースケールまたはカラーテンプレートを対応づけて３次元表示するようにしたが、等距離面上の画素値に限らず、たとえば、表示制御部１３において、等距離面からの厚さ情報を受け付け、その厚さの範囲内の複数の画素値の平均値、最小値または最大値に対してグレースケールまたはカラーテンプレートを対応づけて等距離面の３次元表示を行うようにしてもよい。たとえば平均値を用いて３次元表示を行った場合には、等距離面上の画素値のみを用いた場合と比較すると、厚さ方向についてより広範囲にわたる病変の状況を把握することができる。また、最大値を用いて３次元表示を行った場合には、厚さ情報の範囲内の繊維化（高ＣＴ値）分布情報を把握することができ、最小値を用いて３次元表示を行った場合には、厚さ情報の範囲内の肺気腫（低ＣＴ値）分布情報を把握することができる。
　なお、上記厚さ情報とは、等距離面から肺野の内側方向への厚さの範囲の情報である。厚さ情報は、装置に予め設定しておくようにしてもよいし、ユーザが入力装置４を用いて設定入力するようにしてもよい。また、複数の画素値の平均値、最小値および最大値のいずれかの選択を受け付けるようにしてもよい。
　また、上述したように設定された厚さの範囲内の複数の画素値の平均値、最小値または最大値に対してグレースケールなどを対応づける場合においても、平均値、最小値または最大値に対して不透明度１００％を割り当てて３次元表示することが望ましい。ただし、これに限らず、上述したように１０％～５０％に設定してもよい。
　また、上記実施形態の医用画像診断支援システムにおいては、肺野領域抽出部１１において肺野領域を抽出し、その抽出した肺野領域の表面を基準として等距離面を設定するようにしたが、必ずしも肺野領域の表面を基準としなくてもよく、たとえば、図８に示すように、肺野領域抽出部１１の代わりに胸郭領域抽出部１６を設け、胸郭領域抽出部１６において、３次元医用画像５から胸郭領域を抽出し、その抽出した胸郭領域を基準として等距離面を設定するようにしてもよい。胸郭領域の表面と肺野領域の表面との距離関係は、医学的にある程度決まった関係であるので、その距離関係を考慮して、胸郭領域の表面を基準として肺野領域の内側の等距離面を設定するようにすればよい。なお、胸郭領域の抽出方法については、既に公知は方法を用いることができ、たとえば肋骨を抽出し、その肋骨領域から胸郭領域を決定する方法を用いることができる。たとえば「周 向栄,他、“胸部骨格情報を利用した三次元体幹部マルチスライスＣＴ画像からの胸郭の抽出”電子情報通信学会論文紙 D-II Vol. J88-D-II No.9 pp.1999-2002」に記載の方法を用いることができる。
　また、胸郭領域を抽出する際、横隔膜の領域を抽出し、肋骨領域と横隔膜の領域とから胸郭領域を抽出するようにしてもよい。横隔膜の抽出方法については、既に公知な方法を用いることができ、たとえば「高橋 領，他、“非造影３次元腹部ＣＴ像からの腹腔抽出のための横隔膜抽出処理の開発”、P14/JAMIT2009 Proceedings」に記載の方法を用いることができる。
　また、上記実施形態の医用画像診断支援システムにおいては、同一の被検者について、今回撮影した３次元医用画像５と過去に撮影した３次元医用画像５との両方を読出し、これらの３次元医用画像５のそれぞれについて上述した等距離面上の画素値を取得し、今回撮影した３次元医用画像５の等距離面上の画素値の３次元表示と、過去に撮影した３次元医用画像５の等距離面上の画素値の３次元表示とを同時に行い、これらを比較読影可能なようにしてもよい。
１　　 医用画像表示装置
２　　 ３次元医用画像保管サーバ
３　　 ディスプレイ
４　　 入力装置
１０　 ３次元医用画像取得部
１１　 肺野領域抽出部
１２　 曲面設定部
１３　 表示制御部
１４　 肺胞領域抽出部
１５　 表示肺胞領域特定部
１６　 胸郭領域抽出部

Claims (16)

  1. 　肺野を含む領域を撮影した３次元医用画像において、前記肺野の内側に該肺野の表面に沿った曲面を設定する曲面設定部と、
    　該曲面設定部によって設定された曲面上の画素値を用いて該曲面を立体的に表示させる表示制御部とを備えたことを特徴とする３次元医用画像表示制御装置。
  2. 　前記曲面が、前記肺野の表面から等距離の面である請求項１記載の３次元医用画像表示制御装置。
  3. 　前記曲面設定部が、前記肺野の表面からの前記曲面の位置を変更することにより設定するものである請求項１または２記載の３次元医用画像表示制御装置。
  4. 　前記曲面設定部が、前記肺野の表面からの距離が異なる複数の前記曲面を切り替えることにより設定するものである請求項３記載の３次元医用画像表示制御装置。
  5. 　前記表示制御部が、前記肺野の輪郭線または表面を表示させるものである請求項１から４いずれか１項記載の３次元医用画像表示制御装置。
  6. 　前記３次元医用画像から肺胞領域を抽出する肺胞領域抽出部と、
    　該肺胞領域抽出部によって抽出された複数の肺胞領域のうち閾値以上の大きさの肺胞領域を表示対象の肺胞領域として特定する表示肺胞領域特定部とを備え、
    　前記表示制御部が、前記肺胞領域抽出部によって抽出された複数の肺胞領域のうち前記表示対象の肺胞領域のみを表示させるものである請求項１から５いずれか１項記載の３次元医用画像表示制御装置。
  7. 　前記３次元医用画像から肺胞領域を抽出する肺胞領域抽出部と、
    　該肺胞領域抽出部によって抽出された複数の肺胞領域のうち閾値以上の大きさの肺胞領域を特定する表示肺胞領域特定部とを備え、
    　前記表示制御部が、前記肺胞領域抽出部によって抽出された複数の肺胞領域のうち前記特定された肺胞領域を該特定された肺胞領域以外の肺胞領域とは異なる表示態様で表示させるものである請求項１から５いずれか１項記載の３次元医用画像表示制御装置。
  8. 　前記表示制御部が、前記曲面上の画素値に対してグレースケールまたはカラーテンプレートを対応づけることにより前記曲面を表示させるものである請求項１から７いずれか１項記載の３次元医用画像表示制御装置。
  9. 　前記表示制御部が、前記曲面上の画素値に対する不透明度を１００％に設定することにより前記曲面を表示させるものである請求項８記載の３次元医用画像表示制御装置。
  10. 　前記表示制御部が、前記曲面から予め設定された厚さの範囲内の複数の画素値の平均値、最小値または最大値に対してグレースケールまたはカラーテンプレートを対応づけることにより前記曲面を表示させるものである請求項１から７いずれか１項記載の３次元医用画像表示制御装置。
  11. 　前記表示制御部が、前記平均値、最小値または最大値に対する不透明度を１００％に設定することにより前記曲面を表示させるものである請求項１０記載の３次元医用画像表示制御装置。
  12. 　前記３次元医用画像から肺野領域を抽出する肺野領域抽出部を備え、
    　前記曲面設定部が、前記肺野領域を基準として前記曲面を設定するものである請求項１から１１いずれか１項記載の３次元医用画像表示制御装置。
  13. 　前記３次元医用画像から胸郭領域を抽出する胸郭領域抽出部を備え、
    　前記曲面設定部が、前記胸郭領域を基準として前記曲面を設定するものである請求項１から１１いずれか１項記載の３次元医用画像表示制御装置。
  14. 　前記曲面設定部が、入力装置からの入力を受け付けることによって前記曲面の位置の変更指示を受け付けるものである請求項３または４記載の３次元医用画像表示制御装置。
  15. 　曲面設定部と、表示制御部とを備えた３次元医用画像表示装置の作動方法であって、
    　前記曲面設定部が、肺野を含む領域を撮影した３次元医用画像において、前記肺野の内側に該肺野の表面に沿った曲面を設定し、
    　前記表示制御部が、前記曲面設定部によって設定された曲面上の画素値を用いて該曲面を立体的に表示させることを特徴とする３次元医用画像表示装置の作動方法。
  16. 　コンピュータを、肺野を含む領域を撮影した３次元医用画像において、前記肺野の内側に該肺野の表面に沿った曲面を設定する曲面設定部と、
    　該曲面設定部によって設定された曲面上の画素値を用いて該曲面を立体的に表示させる表示制御部として機能させることを特徴とする３次元医用画像表示制御プログラム。
PCT/JP2014/081977 2014-01-22 2014-12-03 ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム WO2015111308A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014-009466 2014-01-22
JP2014009466A JP2015136480A (ja) 2014-01-22 2014-01-22 ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム

Publications (1)

Publication Number Publication Date
WO2015111308A1 true WO2015111308A1 (ja) 2015-07-30

Family

ID=53681128

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/081977 WO2015111308A1 (ja) 2014-01-22 2014-12-03 ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム

Country Status (2)

Country Link
JP (1) JP2015136480A (ja)
WO (1) WO2015111308A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017150497A1 (ja) * 2016-03-02 2017-09-08 地方独立行政法人神奈川県立病院機構 肺野病変の診断支援装置、該装置の制御方法及びプログラム
JP6855850B2 (ja) 2017-03-10 2021-04-07 富士通株式会社 類似症例画像検索プログラム、類似症例画像検索装置及び類似症例画像検索方法
JP7334900B2 (ja) 2019-05-20 2023-08-29 国立研究開発法人理化学研究所 判別装置、学習装置、方法、プログラム、学習済みモデルおよび記憶媒体

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002099925A (ja) * 2000-09-12 2002-04-05 Ge Medical Systems Global Technology Co Llc 画像表示方法および画像表示装置
JP2007097671A (ja) * 2005-09-30 2007-04-19 Toshiba Corp 画像診断支援装置、及び画像診断支援プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002099925A (ja) * 2000-09-12 2002-04-05 Ge Medical Systems Global Technology Co Llc 画像表示方法および画像表示装置
JP2007097671A (ja) * 2005-09-30 2007-04-19 Toshiba Corp 画像診断支援装置、及び画像診断支援プログラム

Also Published As

Publication number Publication date
JP2015136480A (ja) 2015-07-30

Similar Documents

Publication Publication Date Title
JP5643304B2 (ja) 胸部トモシンセシスイメージングにおけるコンピュータ支援肺結節検出システムおよび方法並びに肺画像セグメント化システムおよび方法
US8676298B2 (en) Medical image alignment apparatus, method, and program
JP5486197B2 (ja) 椎骨中心検出装置および方法並びにプログラム
JP5777070B2 (ja) 領域抽出装置、領域抽出方法および領域抽出プログラム
US9449381B2 (en) Methods, systems, and media for generating and analyzing medical images having elongated structures
US20160292864A1 (en) Medical image data processing apparatus and method for determining the presence of an abnormality
JP5415068B2 (ja) カーブした細長い構造の切断面の視覚化
US9675317B2 (en) Interface identification apparatus and method
JP5785120B2 (ja) 医用画像診断支援装置および方法並びにプログラム
JP5826082B2 (ja) 医用画像診断支援装置および方法並びにプログラム
US10970875B2 (en) Examination support device, examination support method, and examination support program
JP2008509773A (ja) フレキシブル３次元回転血管造影−コンピュータ断層撮影融合方法
US10939800B2 (en) Examination support device, examination support method, and examination support program
JP6580446B2 (ja) 画像処理装置及び画像処理方法
JP6458166B2 (ja) 医用画像処理方法及び装置及びシステム及びプログラム
Akkoul et al. 3D Reconstruction of the proximal femur shape from few pairs of x-ray radiographs
WO2015111308A1 (ja) ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム
US20120188240A1 (en) Medical image display apparatus, method and program
US10896501B2 (en) Rib developed image generation apparatus using a core line, method, and program
JP5624350B2 (ja) 医用画像処理装置
US9307948B2 (en) Image processing apparatus, method, and program
US11003946B2 (en) Examination support device, examination support method, and examination support program
KR102588306B1 (ko) 환자의 두개골의 크기와 형태, 촬영 방향 등에 따라 생기는 다양성을 최소화할 수 있는 인공지능 기반 장치 및 방법
EP4258207A1 (en) Automatic rib fracture detection from unfolded scan images
CN113506277B (zh) 图像处理方法及装置

Legal Events

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

Ref document number: 14880177

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14880177

Country of ref document: EP

Kind code of ref document: A1