WO2006057304A1 - X線ct装置及び画像処理装置 - Google Patents

X線ct装置及び画像処理装置 Download PDF

Info

Publication number
WO2006057304A1
WO2006057304A1 PCT/JP2005/021588 JP2005021588W WO2006057304A1 WO 2006057304 A1 WO2006057304 A1 WO 2006057304A1 JP 2005021588 W JP2005021588 W JP 2005021588W WO 2006057304 A1 WO2006057304 A1 WO 2006057304A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
region
ray
unit
contrast
Prior art date
Application number
PCT/JP2005/021588
Other languages
English (en)
French (fr)
Inventor
Yosuke Okamoto
Satoru Nakanishi
Original Assignee
Kabushiki Kaisha Toshiba
Toshiba Medical Systems Corporation
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 Kabushiki Kaisha Toshiba, Toshiba Medical Systems Corporation filed Critical Kabushiki Kaisha Toshiba
Priority to EP05809487A priority Critical patent/EP1825810A4/en
Priority to JP2006547829A priority patent/JP4874810B2/ja
Publication of WO2006057304A1 publication Critical patent/WO2006057304A1/ja
Priority to US11/753,783 priority patent/US7953263B2/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus 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
    • A61B6/504Apparatus 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 for diagnosis of blood vessels, e.g. by angiography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • G06T5/75Unsharp masking
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/404Angiography

Definitions

  • the present invention relates to an X-ray CT apparatus and an image processing apparatus that generate an image inside a subject based on projection data obtained by irradiating the subject with X-rays.
  • multi-slice X-ray CT equipment In response to the strong demand from the medical field that the advancement of X-ray CT equipment has taken higher-definition (high-resolution) and wide-area imaging, multi-slice X-ray CT equipment has been developed in recent years. This has become quite popular.
  • This multi-slice X-ray CT system detects an X-ray source that emits a fan beam X-ray that spreads in the slice direction (longitudinal direction of the bed) and multiple rows (4 rows, 8 rows, 16 rows, etc.)
  • This is a scanner that includes a two-dimensional detector having a structure in which element rows are arranged in the slice direction, and operates this by multi-scan or helical scan. As a result, volume data over a wide range of the subject can be obtained with high accuracy and in a short time compared to a single slice X-ray CT apparatus.
  • An object of the present invention is to reduce blurring generated around an object having a high X-ray attenuation coefficient when an X-ray absorption distribution image inside a subject obtained by an X-ray CT apparatus is displayed.
  • a first aspect of the present invention is an X-ray CT apparatus that collects projection data related to a subject and reconstructs an image inside the subject based on the collected projection data. And said Based on the unit configured to extract a high-contrast region having a relatively high X-ray attenuation coefficient from the image, the position of the extracted high-contrast region, and the point image intensity distribution function unique to the device, A unit configured to generate a blurred image relating to a contrast region and a relatively low X-ray attenuation coefficient! The subtracting the blurred image from the image cover in order to generate a low contrast image relating to a low contrast region.
  • an X-ray CT apparatus comprising a mute configured as described above.
  • a second aspect of the present invention relates to an X-ray CT apparatus that collects projection data related to a subject and reconstructs an image inside the subject based on the collected projection data.
  • a unit that subtracts a blurred image related to a region with a relatively high X-ray attenuation coefficient included in the image from the image In order to generate a low-contrast image related to a region with a relatively low X-ray attenuation coefficient, a unit that subtracts a blurred image related to a region with a relatively high X-ray attenuation coefficient included in the image from the image.
  • a unit that classifies the low-contrast image into a plurality of regions according to CT values, and a unit that replaces pixel values of the low-contrast image with values that are unique to the plurality of classified regions. Provides a CT line.
  • a third aspect of the present invention relates to an X-ray CT apparatus that collects projection data related to a subject and reconstructs an image inside the subject based on the collected projection data.
  • An X-ray CT apparatus comprising: a unit that selects a specific region from the plurality of extracted region candidates based on a distance between the centroids.
  • FIG. 1 is a block diagram of an X-ray CT apparatus that is useful in the present invention.
  • FIG. 2 is a block diagram of an image processing unit useful for the first embodiment of the present invention.
  • FIG. 3 is an explanatory diagram of a transfer function.
  • FIG. 4 is a block diagram of an image processing unit useful for the second embodiment of the present invention.
  • FIG. 5 is a block diagram of a clustering processing unit that works according to the second embodiment of the present invention.
  • FIG. 6 is a diagram closing three regions divided by the clustering process.
  • FIG. 7 is a flowchart showing a region determination processing procedure by the region determination unit of FIG.
  • FIG. 8 is a bottom view of S13 in FIG.
  • FIG. 9 is a diagram showing the regions Rl extracted in S14 of FIG. 7 and the center-of-gravity position Bl calculated in S15.
  • FIG. 10 is a diagram showing the region R2 extracted in S14 of FIG. 7 and the barycentric position B2 calculated in S15.
  • FIG. 11 is a diagram showing the region R3 extracted in S14 of FIG. 7 and the barycentric position B3 calculated in S15.
  • FIG. 1 is a configuration diagram of an X-ray CT apparatus according to the first embodiment.
  • the X-ray CT system 1 includes a gantry 2 configured to collect projection data related to the subject, a bed 3 for moving the subject P, and input and image display for operating the X-ray CT device.
  • An operation console is provided to perform the operation.
  • the gantry 1 includes an X-ray tube 5 and an X-ray detector 6.
  • the X-ray tube 5 and the X-ray detector 6 are mounted on a ring-shaped rotating frame 8 that is rotated by a gantry driving device 7.
  • the bed 3 includes a top plate 8 on which a subject is placed and a top plate driving device 9 for moving the top plate 8.
  • the rotating frame 8 has an opening in its central part. The subject P placed on the top plate 8 is inserted into the opening.
  • the rotation center axis of the rotating frame 8 is defined as the Z axis (slice direction axis), and the plane perpendicular to the Z axis is defined as two orthogonal axes of XY.
  • a tube voltage is applied between the cathode and anode of the X-ray tube 5 from the high voltage generator 10.
  • Filament current is supplied from the high voltage generator 10 to the filament of the X-ray tube 5.
  • X-rays are generated by applying a tube voltage and supplying a filament current.
  • a one-dimensional array type detector or a two-dimensional array type detector (also referred to as a multi-slice type detector) is adopted.
  • the X-ray detection element has a square light receiving surface of 0.5 mm X O. 5 mm, for example.
  • 916 X-ray detection elements are arranged in the channel direction.
  • a two-dimensional array detector has 40 rows arranged in the slice direction.
  • a single row force is a one-dimensional array type detector.
  • the data acquisition device 11 is generally called a DAS (data acquisition system). Yes.
  • the data collection device 11 converts the signal output from the detector 6 for each channel into a voltage signal, amplifies it, and further converts it into a digital signal. This data (raw data) is supplied to the operation console 4 outside the gantry.
  • the preprocessing unit 12 of the operation console 4 performs correction processing such as sensitivity correction on the data (raw data) output from the data collection device 11 and outputs projection data.
  • This projection data is sent to the reconstruction processing unit 13.
  • the reconstruction processing unit 13 reconstructs image data based on projection data collected by, for example, helical scan, volume scan using cone beam X-rays, or a combination thereof, and stores the image data in the image storage unit 14. .
  • the image processing unit 15 generates a display image based on the image data stored in the image storage unit 14. Setting of conditions for displaying an image, setting of a region of interest, and the like are performed based on an input to the input device 16 by an operator. Details of the image processing unit will be described later.
  • the display device 17 displays the image generated by the image processing unit 15.
  • the scan control unit 18 of the operation console 4 is configured so that a scan such as a helical scan is performed based on an input from the operator, the high pressure generator 10, the gantry drive unit 7, the data collection unit 11, the top plate drive unit. 9 is controlled.
  • the operation console 4 may be configured by dedicated hardware, or the same function may be realized by software using a computer!
  • FIG. 2 is a diagram illustrating the configuration of the image processing unit 15 in FIG. 2A is a schematic diagram of a blood vessel cross-sectional structure as an example of an image processing target. It contains vascular wall Pl, stent P2, cosmetic agent (blood flow part) P3, and lipid P4.
  • vascular wall Pl vascular wall
  • stent P2 vascular wall
  • cosmetic agent blood flow part
  • lipid P4 lipid
  • the image processing unit 15 includes a region of interest setting / image extracting unit 19, a high contrast extracting unit 20, a convolution unit 21, a CT value scaling unit 22, a subtractor 23, and an adder 24.
  • the region-of-interest setting 'image extracting unit 19 sets a region of interest based on the input to the input device 16, and extracts and outputs image data of the region of interest from the image storage unit 14.
  • (B) in FIG. 2 represents an image of the region of interest extracted by the region of interest setting / image extracting unit 19.
  • the high contrast extraction unit 20 extracts a high contrast region having a very high X-ray attenuation coefficient, that is, a pixel group having a very high CT value, from the image of the region of interest extracted by the extraction unit 19.
  • High-contrast materials include steps embedded in the subject for treatment. There are metallic instruments such as mint and calcium mineral. Note that contrast agents that are not mainly stents are mainly used as high-contrast substances in low-contrast images, excluding stents, which will be described later. A contrast agent having a CT value between the stent and the vessel wall is selectively used.
  • a differential image of an input image is obtained by obtaining a difference between the value of the pixel of interest and the average value of a plurality of pixels around the pixel of interest, and this differential image is obtained.
  • a binarized image is obtained by applying threshold processing to the edge using the threshold corresponding to the edge of the stent.
  • This binarized image is an image with different pixel values in the area where the high-contrast material will exist and the other areas.
  • this binarized image is an image including information on a position where a high contrast material will exist (hereinafter referred to as a high contrast position image).
  • FIG. 2C shows this high contrast position image.
  • a relatively high-contrast position image can be easily obtained with relatively little arithmetic processing.
  • the convolution unit 21 performs a convolution operation on a high-contrast position image with a point spread intensity function (PSF) unique to the X-ray computed tomography apparatus.
  • PSF point spread intensity function
  • the PSF scans a wire phantom that has a diameter of its heel (0.05 mm) that is smaller than the detector pitch (resolution limit) of the detector 6, for example 0.5 mm, and obtains its projection data force by reconstruction. It is obtained as 2D image (blurred image) data about the wire phantom.
  • the CT value scaling unit 22 standardizes the image (blurred image) output from the convolution unit 21 according to a non-singular value as a CT value, that is, a standard CT value as a contrast agent here. (Scaling). This normalization process is a pixel value level adjustment process for a subtraction process described later. (D) in FIG. 2 represents an image output from the CT value scaling unit 22 (hereinafter referred to as a high contrast image).
  • the subtraction unit 23 performs a subtraction process, and obtains a difference image between the image of the region of interest and the high-contrast image. Ask. As a result, a low-contrast image in which the high-contrast substance and the blur component around it are reduced can be obtained from the image of the region of interest.
  • This low-contrast image is an image showing a substance with a low attenuation coefficient. (E) in Fig. 2 shows this low contrast image.
  • the adding unit 24 adds the position information of the high-contrast substance by adding the high-contrast position image to the low-contrast image.
  • (F) in Fig. 2 shows this output image.
  • This output image is obtained by adding the position information of the high-contrast material to the image obtained by reducing the blur (artifact) due to the high-contrast material from the original image of the region of interest. Is displayed. In this output image, the lipid around the high-contrast substance, which was difficult to identify due to blurring in the original image, can be observed well.
  • the blur component generated around the high-contrast substance can be reduced, so that the substance around the high-contrast substance can be observed well.
  • a stent or calcification is present in a meridian vessel with a diameter of 3 mm to 5 mm, such as the coronary artery of the heart, the surrounding vessel wall, lipid accumulation state, contrast agent state, etc. are observed well. Is possible.
  • the first embodiment may be implemented with various modifications.
  • the high contrast position image is added in the above embodiment, only the low contrast image may be displayed without adding.
  • marks of different colors indicating the position of the high contrast material may be superimposed and displayed.
  • the force described in the X-ray CT apparatus is implemented in a medical image processing apparatus that displays an image based on projection data, CT image, etc. output from the X-ray CT apparatus. You may give it. Further, in the above-described embodiment, the force described for the processing of the two-dimensional image may be performed by performing the above-described processing on a plurality of two-dimensional images.
  • the pixel value of the image when it is a CT value, it may be a value other than the CT value as long as it represents the force X-ray attenuation coefficient described above.
  • the high-contrast position image is obtained by the differential / binary processing.
  • the transfer function representing the blur due to the high-contrast material is inversely convolved with the image, and the high-contrast position image is obtained. Based on this information, the information on the method to obtain and the position of the stent A method for obtaining position information may be used.
  • classification is automatically performed by clustering processing so that similar substances in an image can be easily identified.
  • plaques and lipids are known to have CT values of 100 to 50, blood vessel walls of 50 to 129, and contrast agents in the bloodstream of 130 to 350. ing.
  • CT value fluctuates.
  • the region extraction (region division) method cannot cope with such variations in CT values and cannot be classified well.
  • the region A mainly corresponding to the plaque, the region B mainly corresponding to the blood vessel wall, and the region C mainly corresponding to the contrast agent are extracted from the low contrast image with high accuracy.
  • FIG. 4 is a configuration diagram of the X-ray CT apparatus 1 according to the second embodiment.
  • the clustering processing unit 25 obtains a plurality of predetermined regions based on CT f and performs CT value conversion processing so that each region is displayed with the same luminance or color.
  • the case where the CT value is divided into three regions will be described.
  • the number of these regions may be other numbers. It is also possible to change the number of areas according to the operator's input.
  • FIG. 5 is a configuration diagram of the cluster link processing unit 25.
  • the clustering processing unit 25 includes a histogram generation unit 26, a region determination unit 27, and a CT value conversion unit 28.
  • the histogram generator 26 obtains a histogram of the input low contrast image.
  • the histogram shows the frequency distribution of pixel appearance for each CT value.
  • CT values are distributed in the range between -100 and 537.
  • the region determining unit 27 extracts, for example, three regions A, B, and C from the low contrast image as shown in FIG.
  • Region A mainly corresponds to plaque.
  • Region B mainly corresponds to the vessel wall.
  • Region C mainly corresponds to the contrast agent.
  • FIG. 7 shows three regions A, B, and C extracted from the low-contrast image by the region determination unit 27.
  • the processing procedure is shown.
  • the low contrast image force blood vessel wall region B is extracted (S11).
  • a plurality of tracking lines are set radially from the approximate center of the blood vessel designated by the operator.
  • CT values are tracked along each tracking line. Two positions where the CT value fluctuates relatively large, that is, two positions showing differential values exceeding a predetermined value, are identified as the inner and outer points of the blood vessel wall.
  • the plaque region A is also extracted with a low contrast image force (S12).
  • a pixel group inside the blood vessel wall and having a CT value lower than that of the blood vessel wall is extracted as the plaque region A.
  • a contrast medium region C is extracted (S13-S19).
  • the range of distribution of CT values on the histogram (in Fig. 6, the range from 100 to 537) is divided into N segments SEG to SEG by the k-means method. (S13). Highest CT value
  • SEG is the SEG segment.
  • N is set to 2 or 3 times the number n of regions finally classified.
  • m is a processing variable. m is initialized to 1.
  • the region candidate R corresponding to the segment SEG is extracted from the low-contrast image (S14). Region R is also extracted with low contrast image power with the lowest value of segment SEG as a threshold.
  • the barycentric position B of the extracted region candidate R is calculated (S15).
  • the variable m is incremented by one (S17).
  • "CTmin seg m ⁇ CTmaxZ2" is determined as a stop condition for preventing diffusion.
  • CTm in seg m is the minimum value of segment m
  • CTmax is the maximum CT value of low contrast images.
  • region candidate R corresponding to segment SEG and segment SEG is
  • the value of the region candidate R is also extracted as a value. Extracted region candidate R
  • Region candidate R is identified as region C of the contrast agent (S 19).
  • the displacement of the center of gravity is the threshold Th
  • the largest region candidate R that fits below is selected as a specific region C.
  • the contrast agent region c can be extracted with high accuracy without being greatly affected by the CT value variation factor.
  • the region B of the blood vessel wall can be extracted with high accuracy without being significantly affected by the CT value fluctuation factor.
  • the plaque region A can be extracted with high accuracy without being affected by the factor of CT value fluctuation.
  • the region determination unit 27 may classify the histogram into three regions using the k-mean method.
  • the threshold Th for determining each region is calculated by the k-mean method so that the variance of CT values in each region is equal.
  • FIG. 6 shows an example of area division according to the present embodiment, and CT values from ⁇ 100 to 537 are classified into three areas. Pixels with CT values in the range of region A mainly correspond to the lipid portion. Pixels with CT values in the range of region B mainly correspond to the blood vessel wall. Pixels with CT values in the area C range mainly correspond to the plaque part.
  • the CT value conversion unit 28 replaces the CT values of the pixels included in each of the three regions extracted or classified by the region determination unit 27 with preset eigenvalues for each region. Thereby The same area portion is displayed with the same brightness or color on the display image, and each area is displayed in a form that can be grasped.
  • an image divided into regions according to the CT value is displayed, so that an image that allows a good grasp of different substances such as tissue, lipid, and contrast medium inside the subject are displayed. Can be provided.
  • the threshold value for determining the region is automatically obtained based on the histogram, it is possible to provide a well-classified image even when there is a variation in the CT value.
  • the blur component (artifact component) generated around the high contrast material is reduced and then the clustering process is performed on the low contrast image, the material around the high contrast material Can be observed well.
  • the stents or calcifications in a meridian vessel having a diameter of 3 to 5 mm such as the coronary artery of the heart, the surrounding blood vessel wall, lipid accumulation state, contrast agent state, etc. are observed well. It is possible.
  • the present invention is not limited to the above-described embodiments, and can be embodied by modifying the constituent elements without departing from the spirit of the invention at the stage of implementation.
  • Various inventions can be formed by appropriately combining a plurality of constituent elements disclosed in the embodiments. For example, some components such as all the components shown in the embodiment may be deleted. In addition, the components of different embodiments may be combined as appropriate.
  • the threshold value may be obtained using another clustering method such as a force group averaging method or a Ward method obtained using the k mean method.
  • the clustering process is performed on the image from which the high-contrast material and the surrounding blur are removed. However, such a process is not performed, and the clustering process is directly performed on the original image of the region of interest. It is also possible to generate a display image by processing!
  • the high contrast position image is added, but only the low contrast image may be displayed without adding. Also, different colors indicating the position of high contrast material The mark may be superimposed and displayed.
  • the force is described in the X-ray CT apparatus.
  • the medical image processing apparatus performs image display based on the projection data CT image output from the X-ray CT apparatus. May be.
  • the processing of the three-dimensional image may be performed by performing the above-described processing on the two-dimensional images of the force multiple described for the processing of the two-dimensional image.
  • the case where the pixel value of the image is a CT value has been described, but any value other than the CT value may be used as long as it represents an X-ray attenuation coefficient.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Vascular Medicine (AREA)
  • Optics & Photonics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Physiology (AREA)
  • Quality & Reliability (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

 X線CT装置は、画像から、X線減弱係数の比較的高い高コントラスト領域を抽出するユニット(20)と、抽出した高コントラスト領域の位置と装置固有の点像強度分布関数とに基づいて、高コントラスト領域に関するボケ画像を発生するユニット(21,22)と、X線減弱係数の比較的低い低コントラスト領域に関する低コントラスト画像を発生するために、画像からボケ画像を引き算するユニット(23)とを具備する。

Description

明 細 書
X線 CT装置及び画像処理装置
技術分野
[0001] 本発明は、被検体に X線を曝射して得られた投影データに基づ!/ヽて被検体内部の 画像を生成する X線 CT装置及び画像処理装置に関する。
背景技術
[0002] X線 CT装置の進歩は目覚しぐより高精細(高解像度)且つ広範囲に撮影した 、と いう医療現場からの強い要望に応えて、近年、マルチスライス X線 CT装置が開発さ れ、これがかなり普及してきている。このマルチスライス X線 CT装置は、スライス方向( 寝台の長手方向)に広がり幅を有するファンビーム X線を曝射する X線源と、複数列( 4列、 8列、 16列など)の検出素子列をスライス方向に並べた構造の 2次元検出器と を備え、これをマルチスキャン又はへリカルスキャンで動作させるスキャナである。こ れにより、シングルスライス X線 CT装置に比べて、被検体の広範囲にわたるボリユー ムデータを高精度に且つ短時間のうちに得ることができる。
[0003] このような X線 CT装置により生成された CT画像中の血管は、近傍に高コントラスト 物質 (高い X線減弱係数を有する物質)がある場合、そのボケの影響を受けて、血管 内腔の様子が不鮮明になる。従来、高コントラスト領域の影響を受けずに注目領域の 画質を改善する方法として、高コントラスト領域の CT値を周囲の CT値に合わせた適 当な値に設定して力 画像処理し、その後、高コントラスト領域の値を戻すという方法 が知られている。し力し、この方法では、高コントラスト領域周辺のボケを良好に除去 することはできな力つた (特開平 10— 40372号公報参照)。
発明の開示
[0004] 本発明の目的は、 X線 CT装置により得られた被検体内部の X線吸収分布画像を 表示する場合に、高い X線減弱係数の物体の周辺に生じるボケを低減することにあ る。
[0005] 本発明の第 1局面は、被検体に関する投影データを収集し、前記収集した投影デ ータに基づ ヽて前記被検体の内部の画像を再構成する X線 CT装置にぉ ヽて、前記 画像から、 X線減弱係数の比較的高い高コントラスト領域を抽出するように構成され たユニットと、前記抽出した高コントラスト領域の位置と装置固有の点像強度分布関 数とに基づいて、前記高コントラスト領域に関するボケ画像を発生するように構成され たユニットと、 X線減弱係数の比較的低!、低コントラスト領域に関する低コントラスト画 像を発生するために、前記画像カゝら前記ボケ画像を引き算するように構成されたュ- ットとを具備する X線 CT装置を提供する。
[0006] 本発明の第 2局面は、被検体に関する投影データを収集し、前記収集した投影デ ータに基づ ヽて前記被検体の内部の画像を再構成する X線 CT装置にぉ ヽて、 X線 減弱係数の比較的低 ヽ領域に関する低コントラスト画像を発生するために、前記画 像から、前記画像に含まれる X線減弱係数の比較的高い領域に関するボケ画像を引 き算するユニットと、前記低コントラスト画像を CT値に応じて複数の領域に分類する ユニットと、前記低コントラスト画像の画素値を、前記分類した複数の領域各々に固 有の値に置き換えるユニットとを具備する X線 CT装置を提供する。
[0007] 本発明の第 3局面は、被検体に関する投影データを収集し、前記収集した投影デ ータに基づ ヽて前記被検体の内部の画像を再構成する X線 CT装置にぉ ヽて、前記 画像から、複数の閾値に基づいて、順番に内包するサイズの異なる複数の領域候補 を抽出するユニットと、前記抽出された複数の領域候補各々の重心を計算するュ- ットと、前記重心間の距離に基づいて、前記抽出された複数の領域候補から特定の 領域を選択するユニットとを具備する X線 CT装置を提供する。
図面の簡単な説明
[0008] [図 1]図 1は本発明に力かる X線 CT装置の構成図である。
[図 2]図 2は本発明の第 1実施例に力かる画像処理部の構成図である。
[図 3]図 3は伝達関数の説明図。
[図 4]図 4は本発明の第 2実施例に力かる画像処理部の構成図である。
[図 5]図 5は本発明の第 2実施例に力かるクラスタリング処理部の構成図である。
[図 6]図 6はクラスタリング処理により分割される 3つの領域を閉めす図である。
[図 7]図 7は図 5の領域決定部による領域決定処理手順を示すフローチャートである。
[図 8]図 8は図 7の S13のネ甫足図である。 [図 9]図 9は図 7の S14で抽出される領域 Rlと S 15で計算される重心位置 Blを示す 図である。
[図 10]図 10は図 7の S14で抽出される領域 R2と S 15で計算される重心位置 B2を示 す図である。
[図 11]図 11は図 7の S14で抽出される領域 R3と S 15で計算される重心位置 B3を示 す図である。
発明を実施するための最良の形態
[0009] (第 1実施例)
以下、本発明の第 1実施例について説明する。
図 1は、第 1実施例に係る X線 CT装置の構成図である。 X線 CT装置 1は、被検体 に関する投影データを収集するために構成された架台 2、被検体 Pを載せて移動さ せるための寝台 3、 X線 CT装置の操作するための入力及び画像表示を行うための操 作コンソールを備えて 、る。
[0010] 架台 1は、 X線管 5と X線検出器 6を有する。 X線管 5と X線検出器 6は、架台駆動装 置 7により回転駆動されるリング状の回転フレーム 8に搭載される。寝台 3は、被検体 を載せる天板 8、天板 8を移動させるための天板駆動装置 9を備えている。回転フレ ーム 8はその中央部分に開口を有する。その開口に天板 8に載置された被検体 Pが 挿入される。なお、回転フレーム 8の回転中心軸を Z軸 (スライス方向軸)、 Z軸に直交 する平面を XYの直交 2軸で定義する。
[0011] X線管 5の陰極陽極間には高電圧発生器 10から管電圧が印加される。 X線管 5の フィラメントには高電圧発生器 10からフィラメント電流が供給される。管電圧の印加及 びフィラメント電流の供給により X線が発生される。 X線検出器 6には、 1次元アレイ型 検出器又は 2次元アレイ型検出器 (マルチスライス型検出器ともいう)が採用される。 X線検出素子は例えば 0. 5mm X O. 5mmの正方の受光面を有する。例えば 916個 の X線検出素子がチャンネル方向に配列される。この列がスライス方向に例えば 40 列並設されたものが 2次元アレイ型検出器である。単一の列力 なるものが 1次元ァ レイ型検出器である。
[0012] データ収集装置 11は、一般的に DAS (data acquisition system)と呼ばれて いる。データ収集装置 11は、検出器 6からチャンネルごとに出力される信号を電圧信 号に変換し、増幅し、さらにディジタル信号に変換する。このデータ (生データ)は架 台外部の操作コンソール 4に供給される。
[0013] 操作コンソール 4の前処理部 12は、データ収集装置 11から出力されるデータ(生 データ)に対して感度補正等の補正処理を施して投影データを出力する。この投影 データは再構成処理部 13に送られる。再構成処理部 13は、例えばへリカルスキャン 、コーンビーム X線を使ったボリュームスキャン又はそれらの併用により収集された投 影データに基づいて、画像データを再構成し、画像記憶部 14に記憶する。画像処 理部 15は、画像記憶部 14に記憶された画像データに基づ 、て表示画像を生成する 。画像を表示するための条件の設定、関心領域の設定等は、操作者による入力装置 16への入力に基づいて行われる。画像処理部の詳細は後述する。表示装置 17は、 画像処理部 15で生成された画像を表示する。又、操作コンソール 4のスキャン制御 部 18は、操作者の入力に基づいて、ヘリカルスキャン等のスキャンが行われるように 、高圧発生装置 10、架台駆動装置 7、データ収集装置 11、天板駆動装置 9を制御 する。尚、操作コンソール 4は、専用ハードウェアで構成しても良いし、コンピュータを 用いてソフトウェアで同様の機能を実現しても良!、。
[0014] 図 2は、図 1の画像処理部 15の構成を説明する図である。尚、図 2の (A)は、画像 処理対象の一例としての血管断面構造の模式図である。血管壁 Pl、ステント P2、造 影剤(血流部分) P3、脂質 P4を含んでいる。以下、このような実体の画像を処理した 場合を例として説明する。
[0015] 画像処理部 15は、関心領域設定 ·画像抽出部 19、高コントラスト抽出部 20、コンボ リューシヨン部 21、 CT値スケーリング部 22、減算器 23、加算器 24を備えている。
[0016] 関心領域設定'画像抽出部 19は、入力装置 16への入力に基づいて関心領域を設 定し、その関心領域の画像データを画像記憶部 14から抽出して出力する。図 2の(B )は、関心領域設定 ·画像抽出部 19により抽出された関心領域の画像を表している。
[0017] 高コントラスト抽出部 20は、抽出部 19により抽出された関心領域の画像から、 X線 減弱係数の非常に高い高コントラスト領域、つまり非常に高い CT値を有する画素群 を抽出する。高コントラスト物質としては、治療のために被検体内に埋め込まれたステ ント等の金属性の器具、石灰化したカルシウム等がある。なお、後述するステント等を 除外した低コントラスト画像における高コントラスト物質としては、ステント等ではなぐ 主に造影剤を対象としている。ステントと血管壁との間の CT値を有する造影剤が選 択的に使用される。
[0018] 具体的な抽出処理としては、注目画素の値と、その注目画素周辺の複数画素の値 の平均値との差分を求めることにより入力された画像の微分画像を求め、この微分画 像に対してステントの辺縁に対応する閾値を用いて閾値処理を施すことにより、 2値 化画像を得る。この 2値化画像は、高コントラスト物質が存在するであろう部分と、そ れ以外の部分で画素値が異なる画像となる。つまり、この 2値化画像は、高コントラス ト物質が存在するであろう位置の情報を含んだ画像 (以下、高コントラスト位置画像と いう。)となる。図 2の(C)は、この高コントラスト位置画像を表している。この微分 · 2値 化処理では、比較的少な 、演算処理で容易に高コントラスト位置画像を得ることがで きる。
[0019] コンボリューシヨン部 21は、高コントラスト位置画像に対して、当該 X線コンピュータ 断層撮影装置固有の点像強度分布関数 (PSF:point spread fonction)を畳み込み演 算をする。 PSFは、装置固有のボケ特性を定義する関数として伝達関数又はボケ関 数とも呼ばれる。 PSFは、検出器 6の検出器ピッチ (解像限界)である例えば 0. 5mm より微小な例えばその ΙΖΙΟの径(0. 05mm)を有するワイヤファントムをスキャンし、 その投影データ力も再構成により得られたワイヤファントムに関する 2次元画像 (ボケ 画像)のデータとして得られる。この PSFを高コントラスト位置画像に畳み込むことに より、高コントラスト対象をそれ単独でスキャンして取得したデータ力も再構成した画 像、つまり高コントラスト対象とその周辺に現れるボケとを含んだ画像が得られる。
[0020] CT値スケーリング部 22は、コンボリューシヨン部 21から出力された画像 (ボケ画像) を、 CT値として特異でない値、つまりここでは造影剤としての標準的な CT値に従つ て規格化 (スケール)する。この規格化処理は、後述する減算処理のための画素値の レベル合わせ処理である。図 2の(D)は、 CT値スケーリング部 22から出力される画 像 (以下、高コントラスト画像という)を表している。
[0021] 減算部 23は、減算処理を行い、関心領域の画像と高コントラスト画像の差分画像を 求める。これにより、関心領域の画像から高コントラスト物質及びその周辺のボケ成分 が低減された低コントラスト画像が得られる。この低コントラスト画像は、低い減弱係数 の物質を表した画像である。図 2の (E)は、この低コントラスト画像を表している。
[0022] 加算部 24は、この低コントラスト画像に高コントラスト位置画像を加算することにより 、高コントラスト物質の位置情報を加える。図 2の (F)は、この出力画像を表している。 この出力画像は、元の関心領域の画像から高コントラスト物質によるボケ (アーチファ タト)を低減した画像に、高コントラスト物質の位置情報を加えたものであり、後段の処 理を経て表示装置 17に表示される。この出力画像では、元の画像ではボケの影響 により識別困難だった高コントラスト物質周辺の脂質が良好に観察可能である。
[0023] このような第 1実施例によれば、高コントラスト物質の周辺に生じるボケ成分を低減 することができるため、高コントラスト物質周辺の物質を良好に観察できる。特に、心 臓の冠状動脈のような 3mmから 5mmの径を有する細経血管にステントや石灰化が ある場合において、その周辺の血管壁、脂質の蓄積状態、造影剤の状態等を良好 に観察することが可能である。
[0024] この第 1実施例は種々変形して実施しても良 、。例えば、上記実施例では、高コン トラスト位置画像を加算したが、加算せずに低コントラスト画像のみの表示を行っても 良い。又、高コントラスト物質の位置を示す色の異なるマークを重畳して表示するよう にしても良い。
[0025] さら〖こ、上記実施例では、 X線 CT装置について説明した力 X線 CT装置から出力 された投影データ、 CT画像等に基づ ヽて画像の表示を行う医用画像処理装置で実 施しても良い。さらに、上記実施例では 2次元画像の処理について説明した力 複数 の 2次元画像に上述の処理を行うことにより 3次元画像の処理を行うようにしても良い
[0026] さらに、上記実施例では画像の画素値が CT値である場合にっ 、て説明した力 X 線減弱係数を表すものであれば CT値以外の値であっても良い。又、上記実施例で は微分 · 2値ィ匕処理により高コントラスト位置画像を求めたが、高コントラスト物質によ るボケを表した伝達関数を画像に逆コンボリューシヨンして高コントラスト位置画像を 求める方法や、ステントの位置 '形状の情報を予め記憶しておきこの情報に基づいて 位置情報を求める方法等を用いても良い。
[0027] (第 2実施例)
以下、本発明の第 2実施例について説明する。尚、第 1実施例と同じ部分について は説明を省略する。
[0028] 第 2実施例は、画像中の類似物質が識別し易くなるようにクラスタリング処理により、 自動的に分類を行うものである。理想的な条件下では、プラークや脂質は— 100から 50の CT値、血管壁は 50力ら 129の CT値、血流に含まれる造影剤は 130— 350の CT値を示すことが知られている。し力しながら、実際の検査においては、被検体の体 格、ビームハード-ング、再構成関数、物質のサイズ、関心領域外の物質の状態、造 影剤の濃度等の影響を受けて、 CT値が変動している。このため領域抽出 (領域分割 )を行う方法では、このような CT値の変動に対応できず、良好に分類することができ ない。本実施例では、低コントラスト画像から、主にプラークに対応する領域 A、主に 血管壁に対応する領域 B、主に造影剤に対応する領域 Cを高精度に抽出する。
[0029] 図 4は、第 2実施例に係る X線 CT装置 1の構成図である。図 4では、減算部 23と加 算部 24の間にクラスタリング処理部 25がある。このクラスタリング処理部 25は、 CT f に基づいて所定の複数の領域を求め、各領域が同一の輝度又は色で表示されるよう に CT値の変換処理を行う。尚、本実施例では、 CT値を 3つの領域に分ける場合に ついて説明するが、この領域の数は他の数でも良い。又、操作者の入力に応じて領 域数を変更できるようにしても良 、。
[0030] 図 5は、クラスタリンク処理部 25の構成図である。クラスタリング処理部 25は、ヒスト グラム生成部 26、領域決定部 27、 CT値変換部 28を備えている。ヒストグラム生成部 26は、入力された低コントラスト画像のヒストグラムを求める。ヒストグラムは、 CT値ご との画素出現の頻度分布を表している。ここでは CT値は— 100から 537間での範囲 に分布している。
[0031] 領域決定部 27は、低コントラスト画像から、図 6に示すように、例えば 3つの領域 A, B, Cを抽出する。領域 Aは主にプラークに対応する。領域 Bは主に血管壁に対応す る。領域 Cは主に造影剤に対応する。
[0032] 図 7は領域決定部 27による低コントラスト画像から 3つの領域 A, B, Cを抽出する 処理の手順を示している。まず、低コントラスト画像力 血管壁領域 Bを抽出する(S1 1)。操作者により指定された血管の略中心から放射状に複数の追跡線が設定される 。各追跡線に沿って CT値が追跡される。 CT値の変動が比較的大きな位置、つまり 所定値を超える微分値を示す 2つの位置が血管壁の内点と外点として特定される。 次に、低コントラスト画像力もプラーク領域 Aが抽出される(S12)。血管壁の内側であ つて血管壁よりも低い CT値を有する画素群が、プラーク領域 Aとして抽出される。
[0033] 次に、造影剤の領域 Cが抽出される(S13— S19)。まず、図 8に例示するように、ヒ ストグラム上の CT値が分布する範囲(図 6では— 100から 537までの範囲)が、 k—m eans法により、 N個のセグメント SEG〜SEG に分割される(S13)。 CT値の最も高
1 N
いセグメントを、 SEGとする。 Nは最終的に分類される領域数 n (ここでは n= 3)を超 える整数に設定される。好ましくは、 Nは最終的に分類される領域数 nの 2倍または 3 倍に設定される。 mは処理変数である。 mは 1に初期化される。
[0034] まず、図 9に示すように、セグメント SEGに対応する領域候補 Rが低コントラスト画 像カゝら抽出される(S14)。セグメント SEGの最低値を閾値として領域 Rが低コントラ スト画像力も抽出される。抽出された領域候補 Rの重心位置 Bが計算される(S15) 。 S17において変数 mが 1つインクリメントされる(S17)。なお、 S18において、拡散 防止のための停止条件として、 "CTmin seg m < CTmaxZ2"が判断される。 CTm in seg mはセグメント mの最小値であり、 CTmaxは低コントラスト画像の最大 CT値で ある。つまり、セグメント mの最小値が低コントラスト画像の最大 CT値の 1Z2未満のと き、領域拡大処理が停止される。 S18で Noのとき、 S19が実行される。 S18で Yesの とき、 S14に戻る。
[0035] 図 10に示すように、セグメント SEGとセグメント SEGとに対応する領域候補 Rが
1 2 2 低コントラスト画像カゝら抽出される(S14)。実際には、セグメント SEGの最低値を閾
2
値として領域候補 Rが低コントラスト画像力も抽出される。抽出された領域候補 Rの
2 2 重心位置 Bが計算される(S 15)。この領域候補 Rの重心位置 Bと、一つ前の領域
2 2 2 候補 Rの重心位置 Bとの間の距離 Dが計算され、所定の閾値 Thと比較される(S1
1 1 2
6)。
[0036] S14〜S18の処理は、重心位置の変位量 D が閾値 Thを超過するまでループされ る。つまり、変数 mのインクリメントに従って領域候補 Rは少しずつ拡大される。領域候 補 Rの拡大に伴って重心 Bは変位する。抽出対象が同じときは重心 Bの変位量は比 較的少ない。抽出対象が他の物質を含むように拡大したとき、重心 Bの変位量は比 較的大きくなる。つまり、抽出対象が造影剤であるときは、領域候補 Rが拡大しても、 重心位置はあまり変化しない。しかし、抽出対象に血管壁等が加わったとき、領域候 補 Rは飛躍的に大きくなり、重心位置も大きく変化する。重心位置が大きく変化する セグメントが探索される。
[0037] m= 3において、図 11に示すように、領域候補 Rの重心位置 Bと、一つ前の領域
3 3
候補 Rの重心位置 Bとの間の距離 Dが閾値 Thを超過したとき(S16)、一つ前の領
2 2 3
域候補 Rが造影剤の領域 Cとして特定される(S 19)。つまり重心の変位量が閾値 Th
2
未満に収まる最大の領域候補 Rが特定の領域 Cとして選択される。
2
[0038] このように領域拡大による重心位置の変位に着目することで、 CT値変動要因の影 響をあまり受けずに、造影剤の領域 cを高精度に抽出することができる。また、血管 壁抽出についても同様に、 CT値変化に着目することで、 CT値変動要因の影響をあ まり受けずに、血管壁の領域 Bを高精度に抽出することができる。さらに、プラーク抽 出についても同様に、血管壁内側という条件を設定することにより、 CT値変動要因 の影響をあまり受けずに、プラークの領域 Aを高精度に抽出することができる。
[0039] なお、領域決定部 27は、このヒストグラムを、 k— mean法を用いて 3つの領域に分 類するようにしてもよい。 k— mean法により、各領域内の CT値の分散値が等しくなる ように、各領域を決めるための閾値 Thを求める。各領域内の CT値の分散値が等しく なるよう各領域を決めることにより、類似性の高い物質が同じ領域内に含まれやすく なる。
[0040] 図 6は、本実施例により領域分けされた例を表しており、—100から 537の CT値が 3つの領域に分類されている。領域 Aの範囲の CT値を持つ画素は主に脂質の部分 に対応する。領域 Bの範囲の CT値を持つ画素は主に血管壁の部分に対応する。領 域 Cの範囲の CT値を持つ画素は主にプラークの部分に対応する。
[0041] CT値変換部 28は、領域決定部 27で抽出又は分類した 3つの領域それぞれに含 まれる画素の CT値を、予め設定された領域ごとの固有値に置き換える。それにより 表示画像上で同じ領域部分が同一の輝度又は色で表示され、各領域が把握できる 形態で表示される。
[0042] このような本実施例によれば、 CT値に応じて領域分けされた画像が表示されるた め、被検体内部の組織、脂質、造影剤等の異なる物質を良好に把握できる画像を提 供することができる。又、領域を決めるための閾値をヒストグラムに基づいて自動的に 求めて ヽるため、 CT値の変動がある場合でも良好に分類した画像を提供することが できる。
[0043] 又、この実施例では、高コントラスト物質の周辺に生じるボケ成分 (アーチファクト成 分)を低減させてから、低コントラスト画像を対象としてクラスタリング処理を行っている ため、高コントラスト物質周辺の物質を良好に観察できる。特に、心臓の冠状動脈の ような 3mmから 5mmの径を有する細経血管にステントや石灰化がある場合において 、その周辺の血管壁、脂質の蓄積状態、造影剤の状態等を良好に観察することが可 能である。
[0044] 又、高コントラスト物質が含まれな 、画像に対してクラスタリング処理を行 、、後から 高コントラスト物質の情報を加えているため、低い CT値を有する部分の領域分けが 良好にできる。
[0045] 尚、本発明は上記実施例に限定されるものではなぐ実施段階ではその要旨を逸 脱しない範囲で構成要素を変形して具体ィ匕できる。また、上記実施例に開示されて いる複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば 、実施例に示される全構成要素カゝら幾つかの構成要素を削除してもよい。さらに、異 なる実施例にわたる構成要素を適宜組み合わせてもよ ヽ。
[0046] 例えば、上記実施例では、 k mean法を用いて閾値を求めた力 群平均法、 War d法等の他のクラスタリング手法を用いて閾値を求めるようにしても良い。又、上記実 施例では、高コントラスト物質及びその周辺のボケを除去した画像に対してクラスタリ ング処理を行ったが、このような処理をせず、関心領域の元画像に対して直接クラス タリング処理を行って表示画像を生成するようにしても良!ヽ。
[0047] さらに上記実施例では、高コントラスト位置画像を加算したが、加算せずに低コント ラスト画像のみの表示を行っても良い。又、高コントラスト物質の位置を示す色の異な るマークを重畳して表示するようにしても良い。
[0048] さら〖こ、上記実施例では、 X線 CT装置について説明した力 X線 CT装置から出力 された投影データ CT画像等に基づ ヽて画像の表示を行う医用画像処理装置で実 施しても良い。さらに、上記実施例では、 2次元画像の処理について説明した力 複 数の 2次元画像に上述の処理を行うことにより 3次元画像の処理を行うようにしても良 い。さらに、上記実施例では、画像の画素値が CT値である場合について説明したが 、 X線減弱係数を表すものであれば CT値以外の値であっても良 、。
産業上の利用可能性
[0049] 本発明によれば、 X線 CT装置により得られた被検体内部の X線吸収分布画像を表 示する場合に、高い X線減弱係数の物体の周辺に生じるボケを低減することができる

Claims

請求の範囲
[1] 被検体に関する投影データを収集し、前記収集した投影データに基づ ヽて前記被 検体の内部の画像を再構成する X線 CT装置において、
前記画像から、 X線減弱係数の比較的高い高コントラスト領域を抽出するように構 成されたユニットと、
前記抽出した高コントラスト領域の位置と装置固有の点像強度分布関数とに基づ いて、前記高コントラスト領域に関するボケ画像を発生するように構成されたユニット と、
X線減弱係数の比較的低い低コントラスト領域に関する低コントラスト画像を発生す るために、前記画像力も前記ボケ画像を引き算するように構成されたユニットとを具備 する X線 CT装置。
[2] 前記高コントラスト領域を抽出するユニットは、前記画像から微分画像を発生し、前記 微分画像に対して閾値処理を行うことにより前記高コントラスト領域を抽出する請求 項 1記載の X線 CT装置。
[3] 前記ボケ画像を発生するユニットは、前記抽出した高コントラスト領域の位置画像に 前記点像強度分布関数を畳み込むユニットと、前記畳み込まれた高コントラスト領域 の位置画像を前記高コントラスト領域の標準的 CT値に従って基準化するユニットとを 有する請求項 1記載の X線 CT装置。
[4] 前記高コントラスト領域は、治療のために被検体内に埋め込まれた金属性の器具又 は石灰化したカルシウムである請求項 1記載の X線 CT装置。
[5] 前記低コントラスト画像に、前記抽出した高コントラスト領域の位置画像を加算するュ ニットをさらに備える請求項 1記載の X線 CT装置。
[6] 前記低コントラスト画像を CT値に応じて複数の領域に分類するユニットと、
前記低コントラスト画像の画素値を、前記分類した複数の領域各々に固有の値に 置き換えるユニットとをさらに備える請求項 1記載の X線 CT装置。
[7] 前記置き換えられた低コントラスト画像に、前記抽出した高コントラスト領域の位置画 像を加算するユニットをさらに備える請求項 6記載の X線 CT装置。
[8] 前記分類した複数の領域には、脂質が支配的な領域と、血管壁が支配的な領域と、 造影剤が支配的な領域とが含まれる請求項 6記載の X線 CT装置。
[9] 前記分類ユニットは、前記低コントラスト画像から、複数の閾値に基づいて、サイズの 異なる複数の領域を抽出するユニットと、
前記抽出された複数の領域各々の重心を計算するユニットと、
前記重心間の距離に基づいて、前記抽出された複数の領域から抽出対象物に対 応する特定の領域を選択するユニットとを有する請求項 6記載の X線 CT装置。
[10] 前記複数の閾値は、前記画像に関する CT値ごとの画素出現の頻度分布に基づい て決定される請求項 9記載の X線 CT装置。
[11] 被検体に関する投影データを収集し、前記収集した投影データに基づいて前記被 検体の内部の画像を再構成する X線 CT装置において、
X線減弱係数の比較的低い領域に関する低コントラスト画像を発生するために、前 記画像から、前記画像に含まれる X線減弱係数の比較的高 ヽ領域に関するボケ画 像を引き算するユニットと、
前記低コントラスト画像を CT値に応じて複数の領域に分類するユニットと、 前記低コントラスト画像の画素値を、前記分類した複数の領域各々に固有の値に 置き換えるユニットとを具備する X線 CT装置。
[12] 前記分類ユニットは、前記低コントラスト画像を、造影剤領域、血管壁領域及びブラ ーク領域に分類する請求項 11記載の X線 CT装置。
[13] 前記分類ユニットは、前記造影剤領域を同定するために、
前記低コントラスト画像から、複数の閾値に基づいて、順番に内包するサイズの異 なる複数の領域候補を抽出するユニットと、
前記抽出された複数の領域候補各々の重心を計算するユニットと、
前記重心間の距離に基づ!、て、前記抽出された複数の領域候補から前記造影剤 領域を選択するユニットとを有する請求項 12記載の X線 CT装置。
[14] 前記複数の閾値は、前記画像に関する CT値ごとの画素出現の頻度分布に基づい て決定される請求項 13記載の X線 CT装置。
[15] 前記選択ユニットは、前記重心の変位量が所定量未満に収まる最大の領域を前記 特定の領域として選択する請求項 13記載の X線 CT装置。
[16] 前記分類ユニットは、前記血管壁領域を同定するために、
前記低コントラスト画像の血管略中心カゝら放射状に CT値を追跡し、前記 CT値の変 化が所定値を超過する 2つの位置を血管の内壁と外壁として特定する請求項 13記 載の X線 CT装置。
[17] 前記分類ユニットは、前記プラーク領域を同定するために、
前記血管壁の内側であって、前記血管壁よりも低い CT値を有する画素群を抽出 する請求項 16記載の X線 CT装置。
[18] 被検体に関する投影データを収集し、前記収集した投影データに基づ ヽて前記被 検体の内部の画像を再構成する X線 CT装置において、
前記画像から、複数の閾値に基づいて、順番に内包するサイズの異なる複数の領 域候補を抽出するユニットと、
前記抽出された複数の領域候補各々の重心を計算するユニットと、
前記重心間の距離に基づいて、前記抽出された複数の領域候補から特定の領域 を選択するユニットとを具備する X線 CT装置。
[19] 前記複数の閾値は、前記画像に関する CT値ごとの画素出現の頻度分布に基づい て決定される請求項 18記載の X線 CT装置。
[20] 前記選択ユニットは、前記重心の変位量が所定量未満に収まる最大の領域を前記 特定の領域として選択する請求項 18記載の X線 CT装置。
[21] X線画像から、 X線減弱係数の比較的高い高コントラスト領域を抽出するように構成さ れたユニットと、
前記抽出した高コントラスト領域の位置と装置固有の点像強度分布関数とに基づ いて、前記高コントラスト領域に関するボケ画像を発生するように構成されたユニット と、
X線減弱係数の比較的低い低コントラスト領域に関する低コントラスト画像を発生す るために、前記画像力も前記ボケ画像を引き算するように構成されたユニットとを具備 する画像処理装置。
[22] X線画像から、前記 X線画像に含まれる X線減弱係数の比較的高い領域に関するボ ケ画像を引き算し、 X線減弱係数の比較的低い領域に関する低コントラスト画像を発 生するユニットと、
前記低コントラスト画像を画素値に応じて複数の領域に分類するユニットと、 前記低コントラスト画像の画素値を、前記分類した複数の領域各々に固有の値に 置き換えるユニットとを具備する画像処理装置。
X線画像から、複数の閾値に基づいて、順番に内包するサイズの異なる複数の領域 候補を抽出するユニットと、
前記抽出された複数の領域候補各々の重心を計算するユニットと、
前記重心間の距離に基づいて、前記抽出された複数の領域候補から特定の領域 を選択するユニットとを具備する画像処理装置。
PCT/JP2005/021588 2004-11-26 2005-11-24 X線ct装置及び画像処理装置 WO2006057304A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP05809487A EP1825810A4 (en) 2004-11-26 2005-11-24 X-RAY CT DEVICE AND IMAGE PROCESSING DEVICE
JP2006547829A JP4874810B2 (ja) 2004-11-26 2005-11-24 X線ct装置
US11/753,783 US7953263B2 (en) 2004-11-26 2007-05-25 X-ray CT apparatus and image processing apparatus

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2004342272 2004-11-26
JP2004-342272 2004-11-26
JP2004342273 2004-11-26
JP2004-342273 2004-11-26

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/753,783 Continuation US7953263B2 (en) 2004-11-26 2007-05-25 X-ray CT apparatus and image processing apparatus

Publications (1)

Publication Number Publication Date
WO2006057304A1 true WO2006057304A1 (ja) 2006-06-01

Family

ID=36498040

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/021588 WO2006057304A1 (ja) 2004-11-26 2005-11-24 X線ct装置及び画像処理装置

Country Status (4)

Country Link
US (1) US7953263B2 (ja)
EP (2) EP2284794B1 (ja)
JP (2) JP4874810B2 (ja)
WO (1) WO2006057304A1 (ja)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL1034681C2 (nl) * 2006-11-17 2009-03-04 Gen Electric Werkwijze en apparatuur voor het aangeven van contouren van ten minste één bloedvat.
JP2009118985A (ja) * 2007-11-14 2009-06-04 Ge Medical Systems Global Technology Co Llc X線ct装置
JP2009545800A (ja) * 2006-08-02 2009-12-24 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 画像中のボクセルのクラスター・マップを再構成する方法
JP2010011980A (ja) * 2008-07-02 2010-01-21 Toshiba Corp 画像処理装置及び画像診断装置
JP2010036041A (ja) * 2008-08-05 2010-02-18 Medison Co Ltd 弾性情報のカラーマップを形成する超音波システム及びカラーマップ形成方法
JP2010193975A (ja) * 2009-02-23 2010-09-09 Toshiba Corp 医用画像処理装置及び医用画像処理方法
JP2011115481A (ja) * 2009-12-07 2011-06-16 Hitachi Medical Corp 画像解析装置及びその方法、画像解析プログラム
JP2012200371A (ja) * 2011-03-24 2012-10-22 Toshiba Corp プラーク領域抽出方法及びその装置
JP2015503389A (ja) * 2011-12-27 2015-02-02 コーニンクレッカ フィリップス エヌ ヴェ 3dスキャンからのem場発生器によるアーティファクトの除去
JP2015112487A (ja) * 2013-12-06 2015-06-22 株式会社東芝 医用画像における構造物をセグメンテーションする医用画像処理装置、医用画像をセグメンテーションするための方法及び医用画像をセグメンテーションするコンピュータプログラムを記憶する記憶媒体
JP2016198504A (ja) * 2015-04-08 2016-12-01 株式会社東芝 画像生成装置、x線コンピュータ断層撮影装置及び画像生成方法
JP2017018177A (ja) * 2015-07-07 2017-01-26 キヤノン株式会社 画像処理装置及び画像処理方法
US9706973B2 (en) 2012-08-30 2017-07-18 Toshiba Medical Systems Corporation Medical image processing apparatus and X-ray computed tomography apparatus
JPWO2016104082A1 (ja) * 2014-12-26 2017-10-05 株式会社日立製作所 画像処理装置及び画像処理方法
JPWO2019003506A1 (ja) * 2017-06-30 2020-03-26 株式会社島津製作所 断層像生成方法および放射線撮影装置

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1859412B1 (en) * 2005-02-14 2019-04-10 Mayo Foundation For Medical Education And Research Electronic stool subtraction in virtual colonoscopy
US20080226082A1 (en) * 2007-03-12 2008-09-18 Storage Appliance Corporation Systems and methods for secure data backup
CN101677782B (zh) * 2007-05-09 2013-05-15 株式会社日立医药 磁共振成像装置及图像分类方法
US8090171B2 (en) * 2007-10-19 2012-01-03 Siemens Medical Solutions Usa, Inc. Image data subtraction system suitable for use in angiography
US9760677B2 (en) 2009-04-29 2017-09-12 Onemednet Corporation Methods, systems, and devices for managing medical images and records
US9171344B2 (en) 2007-10-30 2015-10-27 Onemednet Corporation Methods, systems, and devices for managing medical images and records
US8065166B2 (en) 2007-10-30 2011-11-22 Onemednet Corporation Methods, systems, and devices for managing medical images and records
RU2479038C2 (ru) * 2007-12-20 2013-04-10 Конинклейке Филипс Электроникс Н.В. Трехмерная реконструкция тела и контура тела
JP5305821B2 (ja) * 2008-10-10 2013-10-02 株式会社東芝 医用画像処理装置及び医用画像診断装置
US9478050B2 (en) * 2011-01-06 2016-10-25 Koninklijke Philips N.V. Imaging system for imaging an object
JP2012200285A (ja) * 2011-03-23 2012-10-22 Toshiba Corp 画像処理装置、x線ct装置、及び画像処理方法
WO2014181383A1 (ja) * 2013-05-10 2014-11-13 株式会社島津製作所 画像処理装置
US9438771B2 (en) 2013-10-08 2016-09-06 Canon Kabushiki Kaisha Image processing apparatus, image pickup apparatus, image pickup system, image processing method, and non-transitory computer-readable storage medium
CN104580879B (zh) 2013-10-09 2018-01-12 佳能株式会社 图像处理设备、图像拾取设备以及图像处理方法
CN105022719B (zh) * 2014-04-23 2019-06-28 Ge医疗系统环球技术有限公司 医学造影系统及方法
US9462987B2 (en) * 2014-12-04 2016-10-11 Siemens Aktiengesellschaft Determining plaque deposits in blood vessels
WO2017047893A1 (en) * 2015-09-15 2017-03-23 Samsung Electronics Co., Ltd. Tomography apparatus and controlling method for the same
KR20170088681A (ko) 2016-01-25 2017-08-02 삼성전자주식회사 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법
KR101873109B1 (ko) 2016-10-17 2018-06-29 경희대학교 산학협력단 관심 영역의 동적 설정이 가능한 2차원 엑스선 디텍터와 이를 포함하는 콘빔 ct 장치 및 그 동작 방법
US20210049809A1 (en) 2019-08-12 2021-02-18 Canon Medical Systems Corporation Image processing method and apparatus
CN111214252B (zh) * 2020-01-07 2023-05-02 东软医疗系统股份有限公司 Ct数据的采集方法及装置、图像重建方法及系统
JP7306552B1 (ja) 2022-09-26 2023-07-11 東洋インキScホールディングス株式会社 積層体の分離回収方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001137228A (ja) * 1999-09-07 2001-05-22 General Electric Co <Ge> 2回パス式コーン・ビーム画像再構成の方法及び装置
JP2001190550A (ja) * 2000-01-12 2001-07-17 Univ Nihon X線ct撮影方法及びその装置
JP2004065975A (ja) * 2002-07-31 2004-03-04 Ge Medical Systems Global Technology Co Llc プラーク特徴評価の方法、システム及びコンピュータ・プロダクト
JP2004174218A (ja) * 2002-10-01 2004-06-24 Japan Science & Technology Agency 画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1040372A (ja) 1996-07-18 1998-02-13 Shimadzu Corp 画像処理装置
US6920238B1 (en) * 1996-12-03 2005-07-19 Synchrotronics, Co. Precision imaging system
US7236618B1 (en) * 2000-07-07 2007-06-26 Chee-Kong Chui Virtual surgery system with force feedback
DE10249643A1 (de) * 2002-10-24 2004-05-13 Siemens Ag Verfahren zur Unterstützung der Diagnose und/oder Therapie einer krankhaften Veränderung eines Blutgefäßes und hierzu hergerichtete Datenverarbeitungseinrichtung
US20040130546A1 (en) * 2003-01-06 2004-07-08 Porikli Fatih M. Region growing with adaptive thresholds and distance function parameters
AU2003219634A1 (en) * 2003-02-27 2004-09-17 Agency For Science, Technology And Research Method and apparatus for extracting cerebral ventricular system from images
CN1261910C (zh) 2003-12-08 2006-06-28 西安理工大学 彩色多窗ct图像的自动生成方法
US7783096B2 (en) * 2005-10-17 2010-08-24 Siemens Corporation Device systems and methods for imaging
US7840047B2 (en) * 2006-04-06 2010-11-23 Siemens Aktiengesellschaft Method for reproducibly generating views of tomographic image data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001137228A (ja) * 1999-09-07 2001-05-22 General Electric Co <Ge> 2回パス式コーン・ビーム画像再構成の方法及び装置
JP2001190550A (ja) * 2000-01-12 2001-07-17 Univ Nihon X線ct撮影方法及びその装置
JP2004065975A (ja) * 2002-07-31 2004-03-04 Ge Medical Systems Global Technology Co Llc プラーク特徴評価の方法、システム及びコンピュータ・プロダクト
JP2004174218A (ja) * 2002-10-01 2004-06-24 Japan Science & Technology Agency 画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
See also references of EP1825810A4
VAN STRATEN MARCEL ET AL.: "Removal of arterial wall calcifications in CT angiography by local subtraction", MEDICAL PHYSICS, AIP, MELVILLE, NY, US LNKD- DOI: 10.1118/ 1.1567271, vol. 30, no. 5, 1 May 2003 (2003-05-01), pages 761 - 770

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009545800A (ja) * 2006-08-02 2009-12-24 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 画像中のボクセルのクラスター・マップを再構成する方法
NL1034681C2 (nl) * 2006-11-17 2009-03-04 Gen Electric Werkwijze en apparatuur voor het aangeven van contouren van ten minste één bloedvat.
JP2009118985A (ja) * 2007-11-14 2009-06-04 Ge Medical Systems Global Technology Co Llc X線ct装置
JP2010011980A (ja) * 2008-07-02 2010-01-21 Toshiba Corp 画像処理装置及び画像診断装置
JP2010036041A (ja) * 2008-08-05 2010-02-18 Medison Co Ltd 弾性情報のカラーマップを形成する超音波システム及びカラーマップ形成方法
US9629605B2 (en) 2008-08-05 2017-04-25 Samsung Medison Co., Ltd. Formation of a color map for an elastic image
US9330481B2 (en) 2009-02-23 2016-05-03 Kabushiki Kaisha Toshiba Medical image processing apparatus and medical image processing method
JP2010193975A (ja) * 2009-02-23 2010-09-09 Toshiba Corp 医用画像処理装置及び医用画像処理方法
JP2011115481A (ja) * 2009-12-07 2011-06-16 Hitachi Medical Corp 画像解析装置及びその方法、画像解析プログラム
JP2012200371A (ja) * 2011-03-24 2012-10-22 Toshiba Corp プラーク領域抽出方法及びその装置
JP2015503389A (ja) * 2011-12-27 2015-02-02 コーニンクレッカ フィリップス エヌ ヴェ 3dスキャンからのem場発生器によるアーティファクトの除去
US9706973B2 (en) 2012-08-30 2017-07-18 Toshiba Medical Systems Corporation Medical image processing apparatus and X-ray computed tomography apparatus
JP2015112487A (ja) * 2013-12-06 2015-06-22 株式会社東芝 医用画像における構造物をセグメンテーションする医用画像処理装置、医用画像をセグメンテーションするための方法及び医用画像をセグメンテーションするコンピュータプログラムを記憶する記憶媒体
JPWO2016104082A1 (ja) * 2014-12-26 2017-10-05 株式会社日立製作所 画像処理装置及び画像処理方法
JP2016198504A (ja) * 2015-04-08 2016-12-01 株式会社東芝 画像生成装置、x線コンピュータ断層撮影装置及び画像生成方法
JP7187131B2 (ja) 2015-04-08 2022-12-12 キヤノンメディカルシステムズ株式会社 画像生成装置、x線コンピュータ断層撮影装置及び画像生成方法
JP2017018177A (ja) * 2015-07-07 2017-01-26 キヤノン株式会社 画像処理装置及び画像処理方法
JPWO2019003506A1 (ja) * 2017-06-30 2020-03-26 株式会社島津製作所 断層像生成方法および放射線撮影装置

Also Published As

Publication number Publication date
EP1825810A4 (en) 2010-07-14
JPWO2006057304A1 (ja) 2008-06-05
US20070230653A1 (en) 2007-10-04
JP4874810B2 (ja) 2012-02-15
JP2012005891A (ja) 2012-01-12
JP5159937B2 (ja) 2013-03-13
EP2284794B1 (en) 2017-01-18
EP2284794A3 (en) 2014-08-13
US7953263B2 (en) 2011-05-31
EP2284794A2 (en) 2011-02-16
EP1825810A1 (en) 2007-08-29

Similar Documents

Publication Publication Date Title
JP4874810B2 (ja) X線ct装置
JP4402435B2 (ja) 軟組織空間の視覚化の方法及び装置
JP4619781B2 (ja) プラークの成分を検出するためのシステム
US6035012A (en) Artifact correction for highly attenuating objects
US7702141B2 (en) Method for quantifying an object in a larger structure using a reconstructed image
NL1024447C2 (nl) Werkwijzen en inrichtingen voor afknottingscompensatie.
US7340027B2 (en) Metal artifact correction in computed tomography
JP5442530B2 (ja) 画像処理装置、画像表示装置およびプログラム並びにx線ct装置
EP2891455B1 (en) Radiation tomographic image-generating apparatus and radiation tomographic image-generating method
US20020114503A1 (en) Method and apparatus for processing a computed tomography image of a lung obtained using contrast agent
JP2004174264A (ja) 計算機式断層写真法(ct)スカウト画像を形成する方法及び装置
JP2007520300A (ja) ボクセル組織クラスを分散するctにおける高減衰オブジェクトにより生じる画像全体のアーティファクトの縮小
EP1048008A1 (en) Methods and apparatus for calcification scoring
JP2005511128A (ja) 灌流測定のための低レベル信号の補正
CN100589760C (zh) X射线ct装置和图像处理装置
US20070183644A1 (en) Image processing method and computer readable medium for image processing
JP5329103B2 (ja) 画像処理装置及びx線ct装置
US20030083565A1 (en) Method and apparatus of determining and displaying a helical artifact index
CN107680078B (zh) 一种图像处理方法和装置
US7359476B2 (en) Radiation tomographic imaging apparatus and radiation tomographic imaging method, and image producing apparatus
WO2008065394A1 (en) Method and apparatus for reducing distortion in a computed tomography image
JP7258744B2 (ja) スペクトルコンピュータ断層撮影フィンガープリンティング
US20050018889A1 (en) Systems and methods for filtering images
CN112545541B (zh) 医用图像处理装置以及医用图像处理方法
CN114052754A (zh) 用于对图像进行伪影检测的系统和方法

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

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

AL Designated countries for regional patents

Kind code of ref document: A1

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

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

Ref document number: 2006547829

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2005809487

Country of ref document: EP

Ref document number: 11753783

Country of ref document: US

Ref document number: 200580040425.4

Country of ref document: CN

NENP Non-entry into the national phase

Ref country code: DE

WWP Wipo information: published in national office

Ref document number: 2005809487

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 11753783

Country of ref document: US