JP2015112481A - 二重エネルギーctスキャン画像における物体識別方法 - Google Patents

二重エネルギーctスキャン画像における物体識別方法 Download PDF

Info

Publication number
JP2015112481A
JP2015112481A JP2014246381A JP2014246381A JP2015112481A JP 2015112481 A JP2015112481 A JP 2015112481A JP 2014246381 A JP2014246381 A JP 2014246381A JP 2014246381 A JP2014246381 A JP 2014246381A JP 2015112481 A JP2015112481 A JP 2015112481A
Authority
JP
Japan
Prior art keywords
pixel
pixel group
pixels
value
contrast agent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2014246381A
Other languages
English (en)
Other versions
JP6438752B2 (ja
Inventor
チャンイン・リ
Jianying Li
シュオ・リ
Shuo Li
ジーフイ・スン
Zhihui Sun
ジー・イェ
Zhi Ye
シャオヤン・ユ
Xiaoyan Yu
ケ・スン
Ke Sun
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of JP2015112481A publication Critical patent/JP2015112481A/ja
Application granted granted Critical
Publication of JP6438752B2 publication Critical patent/JP6438752B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/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/48Diagnostic techniques
    • A61B6/482Diagnostic techniques involving multiple energy imaging
    • 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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices 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
    • 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/10Segmentation; Edge detection
    • G06T7/174Segmentation; Edge detection involving the use of two or more images
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray 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
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

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

Abstract

【課題】二重エネルギーCTスキャン画像を処理するための方法を提供する。【解決手段】その方法は、二重エネルギーCTスキャン画像内のピクセルをフィルタリングし、グループ化されることになるピクセルを得ることと、二重エネルギーCTスキャン画像内のグループ化されることになるピクセルの位置に基づいて、グループ化されることになるピクセルを複数のピクセルグループにグループ化することと、各ピクセルグループ内のピクセルに関する物質分解を実行することと、物質分解の結果に基づいて各ピクセルグループに対応する物体を特定することとを含む。その方法を用いることによって、スキャン時間が短縮され、対象物(例えば、診断されることになるユーザ)に照射されるX線の線量が削減される。【選択図】図1

Description

本発明の実施形態は包括的にはコンピュータ断層撮影(CT)の分野に関し、より具体的には、二重エネルギーCTスキャン画像における物体識別方法に関する。
コンピュータ断層撮影(CT)技術を通して得られた対象物体のスキャン画像は、異なるグレースケールからなる複数のピクセルを含む。ピクセルのグレースケールは、スキャンされることになる対象物体の密度に比例する。スキャンされることになる対象物内の密度が異なるエリアでは、X線の吸収が異なる。それゆえ、CTスキャン画像では、低いグレースケールを有するピクセルを用いて、X線吸収レベルが低い低密度のエリアを表し、高いグレースケールを有するピクセルを用いて、X線吸収レベルが高い高密度のエリアを表す。通常、CT値(単位:Hu)を用いて、CTスキャン画像内のピクセルを表す。例えば、X線吸収レベルがより高い物体又は組織に対応するピクセルのCT値ほど、X線吸収レベルがより低い物体又は組織に対応するピクセルのCT値より高く設定される。例えば、骨組織に対応するピクセルのCT値は+1000Huと設定され、水に対応するピクセルのCT値は0Huと設定され、空気に対応するピクセルのCT値は−1000Huと設定されるなどである。
現在、二重エネルギーCT技術が提案されている。この技術では、2つの異なるエネルギーのX線を用いて、スキャン対象物をスキャンし、それにより、対象物の二重エネルギーCTスキャン画像を得る。そのような二重エネルギーCTスキャン画像は、単一エネルギーのX線を利用して対象物をスキャンする従来の単一エネルギーCT技術を用いることによって得られたスキャン画像より多くの情報を含む。
しかしながら、単一エネルギーCTスキャン画像及び二重エネルギーCTスキャン画像のいずれにおいても、対象物内の異なる物体が同じ、又は類似のX線吸収レベルを有する場合がある。例えば、骨組織と、病変に起因して発生した血管又は他の組織のような石灰化した部分とは同じ、又は類似のX線吸収レベルを有し、それゆえ、同じ、又は類似のCT値を有する場合がある。それゆえ、CTスキャン画像内の同じ、又は類似のグレースケール値のピクセルによって表される物体を区別することは難しい。
さらに、対象物内の異なる物体をハイライトするために、CTスキャンを実行する前に、対象物に造影剤が投与(例えば、注入)される。物体に造影剤が投与されたCTスキャンは、CT造影スキャンとも呼ばれ、得られたCTスキャン画像は、CT造影スキャン画像とも呼ばれる。
しかしながら、CT造影スキャン画像では、造影剤に対応するピクセルが相対的に高いCT値を有するので、造影剤に対応するピクセルと、相対的に高い密度及び相対的に高いX線吸収レベルを有する骨組織及び他の組織に対応するピクセルとを区別するのが難しい。従来技術において、造影剤を投与する前のCTスキャン(非造影スキャン)によって得られたCTスキャン画像(非造影画像)と、CT造影スキャン画像(造影スキャン画像)とを比較することによって、造影剤に対応するピクセルを識別するための方法が存在する。しかしながら、そのような方法では、対象物を2回CTスキャンする必要があるので、対象物によるX線吸収線量が増加する。
それゆえ、CTスキャン画像内の異なる物体を識別するための方法が必要とされている。
本発明の例示的な実施形態の目的は、従来技術における上記の問題及び/又は他の問題を克服することである。それゆえ、本発明の例示的な実施形態は、異なる物体を容易に区別することができる、二重エネルギーCTスキャン画像を処理するための方法を提供することである。
本発明の例示的な実施形態によれば、二重エネルギーCT造影スキャン画像内の造影剤を識別するための方法が提供され、その方法は、二重エネルギーCT造影スキャン画像内のピクセルをフィルタリングし、グループ化されることになるピクセルを得ることと、二重エネルギーCT造影スキャン画像内のグループ化されることになるピクセルの位置に基づいて、グループ化されることになるピクセルを複数のピクセルグループにグループ化することと、各ピクセルグループ内のピクセルに関して物質分解を実行することと、物質分解の結果に基づいて、複数のピクセルグループ内の造影剤に対応するピクセルグループを特定することとを含む。
本発明の例示的な実施形態によれば、二重エネルギーCTスキャン画像を処理するための方法が提供され、その方法は、二重エネルギーCTスキャン画像内のピクセルをフィルタリングし、グループ化されることになるピクセルを得ることと、二重エネルギーCTスキャン画像内のグループ化されることになるピクセルの位置に基づいて、グループ化されることになるピクセルを複数のピクセルグループにグループ化することと、各ピクセルグループ内のピクセルに関して物質分解を実行することと、物質分解の結果に基づいて各ピクセルグループに対応する物体を特定することとを含む。
本発明は、添付の図面を参照しながら、本発明の例示的な実施形態の以下の説明を読むことから更に深く理解される。
例示的な実施形態による、二重エネルギーCT造影スキャン画像における造影剤識別方法の流れ図である。 例示的な実施形態による、二重エネルギーCT造影スキャン画像の一例を示す図である。 例示的な実施形態による、フィルタリングされたスキャン画像の一例を示す図である。 例示的な実施形態による、異なるピクセルグループに対応する分解散布図の一例を示す図である。 例示的な実施形態による、異なるピクセルグループに対応する分解散布図の一例を示す図である。 例示的な実施形態による、CT値を下げたピクセルを含む、処理されたスキャン画像の一例を示す図である。 別の例示的な実施形態による、二重エネルギーCTスキャン画像における物体識別方法の流れ図である。
以下の詳細な説明において、本発明の実施形態が説明される。これらの実施形態の詳細な説明では、簡潔にするために、実際の実施態様の全ての特徴が詳細に説明されるとは限らないことに留意されたい。1つの実施態様を実際に実施する際に、任意の工学又は設計プロジェクトと同様に、開発者の具体的な目標を達成するために、かつシステム関連又はビジネス関連の制約を満たすために、種々の具体的な決定が行われることになるので、実施される実施形態が状況に応じて異なることは理解されたい。さらに、そのような開発のプロセスにおいて行われる努力は、本発明の開示に関連する分野の当業者にとって複雑で、非常に長い場合があるが、本開示において開示される技術的内容に基づいて行われる製造又は生産の設計及び変形は、単に従来の技術的手段に関連するものであり、それゆえ、本開示が不十分あると見なされるべきではないことも理解されよう。
別に規定されない限り、特許請求の範囲及び説明において用いられるような技術用語又は科学用語は、当業者によって一般的に理解される意味を有するように解釈されるべきである。説明及び特許請求の範囲において、「第1の(first)」、「第2の(second)」などの用語は、任意の順序、量又は重要性を表すこと意図するものではなく、単に、異なる構成要素を区別するために用いられる。「1つの(a、an)」などの用語は、数量制限を与えるものではなく、少なくとも1つものが存在することを指示する。「備える(comprise)」、「含む(include)」などの用語は、その用語の後に列挙されない他の要素又は物体を除外しない。「接続する(connect)」、「接続(connection)」などの用語は、物理的若しくは機械的な接続、又は直接的若しくは間接的な接続には限定されない。
図1は、例示的な実施形態による、二重エネルギーCT造影スキャン画像内の造影剤識別方法の流れ図である。
図1に示されるように、ステップS110において、対象物(例えば、診断されることになるユーザ)の二重エネルギーCT造影スキャン画像をフィルタリングすることができる。図2は、例示的な実施形態による、二重エネルギーCT造影スキャン画像の一例を示す。画像のフィルタリングを実行する際に、図2に示されるような二重エネルギーCT造影スキャン画像内のピクセルのCT値が基準CT値より大きいか否かに関する判断を行うことができ、基準CT値より大きいCT値を有するピクセルをフィルタリングされたピクセルと判断することができる。フィルタリングされたピクセルはその後、グループ化されるので(後に詳細に説明される)、フィルタリングされたピクセルは、グループ化されることになるピクセルと呼ぶこともできる。ここで、基準CT値はあらかじめ選択された値とすることができる。この例示的な実施形態では、造影剤に対応するピクセルを骨組織のような他の物体に対応するピクセルと区別するために、基準CT値は100Huとして選択することができる。しかしながら、例示的な実施形態はそれには限定されず、二重エネルギーCT造影スキャン画像は異なる基準CT値を用いてフィルタリングすることができることは理解されたい。例えば、基準CT値より小さいCT値を有するピクセルをグループ化されることになるピクセルと見なすことができるか、又は第1の基準CT値より大きく、かつ第2の基準CT値より小さいCT値を有するピクセルを、グループ化されることになるピクセルと見なすことができる(第1の基準CT値は第2の基準CT値より小さい)。
図3は、例示的な実施形態による、フィルタリングされたスキャン画像の一例を示す。
図2に示されるような二重エネルギーCT造影スキャン画像が、100Huである基準CT値を用いてフィルタリングされた後に、図3に示されるように、グループ化されることになるピクセルを単に含むフィルタリングされたスキャン画像を得ることができる。図3において、部分Aは、骨組織に対応するピクセルとすることができ、部分Bは、造影剤(血液中に含まれる造影剤など)に対応するピクセルとすることができる。
再び図1を参照すると、二重エネルギーCT造影スキャン画像がフィルタリングされた後に、結果として生成されたグループ化されることになるピクセルがグループ化される(S130)。グループ化されることになるピクセルは、二重エネルギーCT造影スキャン画像内のグループ化されることになるピクセルの位置に基づいて、複数のピクセルグループに自動的に、又は手作業でグループ化することができる。例えば、自動グループ化の場合、グループ化されることになるピクセルを最初に選択することができ(例えば、任意に選択される)、その後、前記選択されたピクセルに隣接するグループ化されることになるピクセルのうちの1つのピクセルを、選択されたピクセルと同じグループの中にグループ化することができる。そのプロセスは、グループ化されることになる全てのピクセルが対応するピクセルグループにグループ化されるまで繰り返すことができる。このようにして、グループ化されることになるピクセル内で隣接して位置するピクセルを、同じピクセルグループの中に自動的にグループ化することができる。
図1を参照すると、グループ化されることになるピクセルがグループ化された後に、各ピクセルグループ内のピクセルを物質分解にかけることができる(S150)。二重エネルギーCTスキャン画像は単一エネルギーCTスキャン画像より多くの情報を提供することができるので、物質分解は、ピクセルごとに少なくとも2つの分解値を得ることができるように実行することができる。ここで、各ピクセルの2つの分解値はそれぞれ同等の密度値の2つの異なる基底物質とすることができるか、又は二重エネルギーCTスキャンを実行するために用いられる2つの異なるエネルギーに対応するCT値とすることができる。以下では、各ピクセルの2つの分解値がそれぞれ同等の密度値の水、及び造影剤内に含まれる同等の密度値のヨウ素である例示的な実施形態が説明される。
図1に示されるように、各ピクセルグループ内のピクセルのCT値の物質分解後に、物質分解の結果に基づいて、対象物(例えば、骨組織、造影剤など)以外の物体に対応する、複数のピクセルグループ内のピクセルグループを特定することができる(S170)。
例えば、ピクセルグループ内の各ピクセルの2つの分解値間の関係を指示する分解値散布図を各ピクセルグループ内のピクセルの分解値に基づいて確立することができる。図4及び図5は、例示的な実施形態による異なるピクセルグループに対応する分解散布図の例であり、図4は図2の部分A(グループA)に対応するピクセルグループの分解散布図であり、図5は図2の部分B(グループB)に対応するピクセルグループの分解散布図である。
グループA及びグループBは異なる物体(すなわち、それぞれ骨組織及び造影剤)を表すので、図4及び図5に示されるように、分解散布図において、グループA内のピクセルの分布はグループB内のピクセルの分布とは異なる可能性がある。それゆえ、各ピクセルグループに対応する物体を、分解散布図内のピクセルグループのピクセルの位置に基づいて手作業で、又は自動的に特定することができる(S170)。
例えば、各ピクセルグループに対応する物体を自動的に特定する場合、分解値散布図内のピクセルグループのピクセルの位置と、異なる物体を表す特性曲線との間の平均距離に関する計算を行うことができ、前者と特定の特性曲線との間の平均距離が基準値以下である場合には、前記ピクセルグループは前記特性曲線によって表される物体に対応すると判断することができる。そのような特性曲線及び基準値は、あらかじめ求めることができる。しかしながら、例示的な実施形態は、それには限定されないことは理解されたい。他の例示的な実施形態では、分解値散布図内のピクセルグループのピクセルの位置と、特定の物体(例えば、造影剤)を表す特性曲線(例えば、造影剤特性曲線)との間の平均距離、及び分解値散布図内のピクセルグループのピクセルの位置と、二重エネルギーCT造影スキャン画像内の前記特定の物体(例えば、造影剤)以外の物体の特性曲線との間の平均距離に関する計算を行うことができる。分解値散布図内のピクセルグループのピクセルの位置と、特定の物体(例えば、造影剤)の特性曲線との間の平均距離が、分解値散布図内のピクセルグループのピクセルの位置と、他の物体を表す特性曲線との間の平均距離より短い場合には、そのピクセルグループは、特定の物体(例えば、造影剤)に対応すると判断することができる。
図4及び図5は、骨組織を表す特性曲線C1及び造影剤を表す特性曲線C2の例を示す。図4に示されるように、グループA内のピクセルと骨組織を表す特性曲線C1との間の平均距離は、基準値より小さい場合があり、グループA内のピクセルと造影剤を表す特性曲線C2との間の平均距離は基準値より大きい場合がある。それゆえ、グループAは骨組織に対応すると判断することができる。図5に示されるように、グループB内のピクセルと骨組織を表す特性曲線C1との間の平均距離は、基準値より大きい場合があり、グループB内のピクセルと造影剤を表す特性曲線C2との間の平均距離は基準値より小さい場合がある。それゆえ、グループBは造影剤に対応すると判断することができる。
代替的には、グループA内のピクセルと、骨組織を表す特性曲線C1との間の平均距離は、グループA内のピクセルと造影剤を表す特性曲線C2との間の平均距離に比べて短い可能性があるので、グループA内のピクセルを骨組織に対応すると判断することができる。同様に、グループB内のピクセルと造影剤を表す特性曲線C2との間の平均距離は、グループB内のピクセルと骨組織を表す特性曲線C1との間の平均距離に比べて短い可能性があるので、グループB内のピクセルは造影剤に対応すると判断することができる。
それゆえ、例示的な実施形態によれば、二重エネルギーCTスキャン画像内の異なる物体を自動的に、又は手作業で区別し、識別することができる。例えば、二重エネルギーCT造影スキャン画像内の造影剤を正確に識別することができる。
さらに、オプションの例示的な実施形態では、ピクセルグループに対応する物体が特定された後に、二重エネルギーCTスキャン画像内のピクセルのCT値を変更することができる。例えば、複数のピクセルグループから特定の物体(例えば、造影剤)に対応するピクセルグループを選択することができ、前記選択されたピクセルグループ(例えば、造影剤に対応する)内のピクセルのCT値を下げることができる。1つの例示的な実施形態では、そのようにしてピクセルのCT値を下げることができ、すなわち、選択されたピクセルグループ内のピクセルの、2つの分解値(例えば、1つが造影剤に対応する)のうちの1つとしての同等密度のヨウ素が、例えば、2mg/ccより大きいときに、その値を2mg/ccに置き換えることができる。
図6は、例示的な実施形態による、CT値を下げたピクセルを含む、処理されたスキャン画像の一例を示す。図6に示されるように、造影剤に対応すると判断されたエリアB内のピクセルのグレースケールは、CT値を下げることに起因して下げられる。このようにして、図6に示される処理された画像では、造影剤に対応するエリアBは、もはやハイライトされず、骨組織に対応するエリアAのみがハイライトされる。言い換えると、上記の処理を通して、造影されないスキャン画像(非造影スキャン画像)を得ることができる。言い換えると、二重エネルギーCT造影スキャン画像(造影スキャン画像)に基づいて非造影スキャン画像を得ることができる。それゆえ、対象物上で二重エネルギーCT造影スキャンを一度だけ実行することによって、造影スキャン画像及び非造影スキャン画像を得ることができる。それゆえ、スキャン時間が短縮され、対象物(例えば、診断されることになるユーザ)に照射されるX線の線量が削減される。
造影スキャン画像に基づいて非造影スキャン画像を得ることが説明されてきたが、本発明の例示的な実施形態がそれには限定されないことは理解されよう。以下に説明されるように、別の例示的な実施形態では、二重エネルギーCTスキャン画像内の異なる物体を識別し、識別結果に基づいて二重エネルギーCTスキャン画像を処理することができる。
図7は、別の例示的な実施形態による、二重エネルギーCTスキャン画像内の物体識別方法の流れ図を示す。この実施形態では、上記の例示的な実施形態と同一、又は類似の要素又は機構は、繰り返しを避けるために説明されない。
最初に、二重エネルギーCTスキャン画像内のピクセルをフィルタリングして、グループ化されることになるピクセルを得る(S710)。例えば、ピクセルのCT値を基準CT値と比較することができ、基準CT値より大きいCT値を有するピクセルをグループ化されることになるピクセルと判断することができる。
その後、グループ化されることになるピクセルは、二重エネルギーCTスキャン画像内のグループ化されることになるピクセルの位置に基づいて、複数のピクセルグループに自動的に、又は手作業でグループ化することができる(S730)。例えば、自動グループ化の場合、グループ化されることになるピクセル内で隣接して位置するピクセルを同じグループにグループ化することができる。
次に、各ピクセルグループ内のピクセルを物質分解にかけることができる(S750)。ピクセルの物質分解の結果は、同等の密度値の2つの異なる基底物質、又は二重エネルギーCTスキャンを実行するために用いられる2つの異なるエネルギーに対応するCT値とすることができる。
その後、各ピクセルグループ内のピクセルの分解結果に基づいて、各ピクセルグループに対応する物体を自動的に、又は手作業で判断することができる(S770)。例えば、ピクセルグループ内の各ピクセルの物質分解結果間の関係を指示する分解値散布図を確立することができ、その後、分解値散布図内のピクセルグループのピクセルの位置に基づいて、ピクセルグループに対応する物体を自動的に、又は手作業で判断することができる。ピクセルグループに対応する物体を自動的に判断する場合、分解値散布図内のピクセルグループのピクセルの位置と、異なる物体を表す特性曲線との間の平均距離に関する計算を行うことができ、最小平均距離を有する特性曲線によって表される物体を、前記ピクセルグループに対応する物体と判断することができる。代替的には、分解値散布図内のピクセルグループのピクセルの位置と、異なる物体を表す特性曲線との間の平均距離に関する計算を行うことができ、ピクセルグループのピクセルと特定の物体を表す特性曲線との間の平均距離が基準値以下である場合には、前記ピクセルグループを物体に対応すると判断することができる。ここで、異なる物体を表す特性曲線は、骨組織、造影剤又は他の物体を表す特性曲線を含むことができる。
さらに、オプションの例示的な実施形態では、二重エネルギーCTスキャン画像内のピクセルのCT値を、各ピクセルグループに対応する物体の判断結果に従って変更することができる。例えば、複数のピクセルグループから特定の物体に対応するピクセルグループを選択することができ、前記選択されたピクセルグループ内のピクセルのCT値を変更することができる。1つの例示的な実施形態では、選択されたピクセルグループ内のピクセルのCT値を下げることができる。
1つの例示的な実施形態によれば、二重エネルギーCTスキャン画像内の各物体に対応するピクセルを識別し、それにより、造影剤に対応するピクセルのCT値を変更することができる。それゆえ、二重エネルギーCTスキャン画像内の各物体に対応するピクセルのグレースケールを所望により変更することができる。それゆえ、対象物上で二重エネルギーCTスキャンを一度だけ実行することにより、異なる物体がそれぞれハイライトされた画像を得ることができ、それにより、スキャン時間が短縮され、対象物(例えば、診断されることになるユーザ)に照射されるX線の線量が削減される。
本発明の例示的な実施形態が上記で説明されてきたが、種々の変更を加えることができることは理解されたい。例えば、説明された技法が異なる順序において実施される場合でも、及び/又は説明されたシステム、アーキテクチャ、装置又は回路内の構成要素が異なるやり方で組み合わせられ、及び/又は異なる構成要素若しくは同等物によって置き換えられるか、又は補われる場合であっても、適切な結果が達成されることになることは理解されよう。したがって、他の実施態様も特許請求の範囲の保護範囲内に入る。
この書面による説明は、複数の例を用いて、最良の態様を含む、本発明を開示し、それにより、任意のデバイス又はシステムを作製し、使用すること、及び任意の組み込まれた方法を実行することを含む、当業者が本発明を実践できるようにする。本発明の特許可能な範囲は、特許請求の範囲によって規定され、当業者が思い浮かぶ他の例も含むことができる。そのような他の例は、特許請求の範囲の文字通りの言葉と異ならない構造的要素を有する場合には、又は特許請求の範囲の文字通りの言葉とわずかに異なるが、同等の構造的要素を含む場合には、特許請求の範囲に入ることを意図している。
S110 ステップ
S130 ステップ
S150 ステップ
S170 ステップ
S710 ステップ
S730 ステップ
S750 ステップ
S770 ステップ
A 部分、グループ
B 部分、グループ
C1 特性曲線
C2 特性曲線

Claims (18)

  1. 二重エネルギーCT造影スキャン画像内の造影剤を識別するための方法であって、前記方法は、
    前記二重エネルギーCT造影スキャン画像内のピクセルをフィルタリングし、グループ化されることになるピクセルを得ることと、
    前記二重エネルギーCT造影スキャン画像内の前記グループ化されることになるピクセルの位置に基づいて、前記グループ化されることになるピクセルを複数のピクセルグループにグループ化することと、
    各ピクセルグループ内の前記ピクセルに関する物質分解を実行することと、
    物質分解の結果に基づいて、前記複数のピクセルグループ内の前記造影剤に対応するピクセルグループを特定することとを含む、二重エネルギーCT造影スキャン画像内の造影剤を識別するための方法。
  2. フィルタリングする前記ステップは、
    前記ピクセルのCT値が基準CT値より大きいか否かを判断することと、
    そのCT値が前記基準CT値より大きい前記ピクセルを、グループ化されることになるピクセルと判断することとを含む、請求項1記載の方法。
  3. グループ化されることになるピクセルをグループ化する前記ステップは、
    前記グループ化されることになるピクセル内で隣接して位置するピクセルを同じピクセルグループにグループ化することを含む、請求項1記載の方法。
  4. 前記造影剤に対応する前記ピクセルグループを判断する前記ステップは、
    各ピクセルグループ内の前記ピクセルの分解結果に基づいて、前記ピクセルグループ内の各ピクセルの物質分解結果間の関係を指示する分解値散布図を確立することと、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置に基づいて、ピクセルグループが前記造影剤に対応するか否かを判断することとを含む、請求項1記載の方法。
  5. ピクセルグループが前記造影剤に対応するか否かを判断する前記ステップは、
    前記分解値散布図内のピクセルグループのピクセルの位置と前記造影剤を表す造影剤特性曲線との間の平均距離を計算することと、
    計算された距離が基準値以下である場合には、前記ピクセルグループを前記造影剤に対応すると判断することとを含む、請求項4記載の方法。
  6. ピクセルグループが前記造影剤に対応するか否かを判断する前記ステップは、
    前記分解値散布図内のピクセルグループの前記ピクセルの位置と前記造影剤を表す造影剤特性曲線との間の平均距離と、前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と前記二重エネルギーCT造影スキャン画像内の前記造影剤以外の物体を表す特性曲線との間の平均距離とを計算することと、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と前記造影剤特性曲線との間の平均距離が、前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と他の物体を表す前記特性曲線との間の平均距離より小さい場合には、前記ピクセルグループを前記造影剤に対応すると判断することとを含む、請求項4記載の方法。
  7. 他の物体を表す前記特性曲線は、骨組織を表す特性曲線を含む、請求項6記載の方法。
  8. 前記ピクセルの2つの特性値は、それぞれ同等の密度値の2つの異なる基底物質であるか、又は前記二重エネルギーCT造影スキャン画像を得るために、前記二重エネルギーCTスキャンを実行するために用いられる2つの異なるエネルギーに対応するCT値である、請求項3記載の方法。
  9. 前記基準CT値は100Huである、請求項1記載の方法。
  10. 前記造影剤に対応する前記ピクセルグループを判断した結果に基づいて、前記造影剤に対応する前記ピクセルグループのCT値を変更することを更に含む、請求項1記載の方法。
  11. 二重エネルギーCTスキャン画像内の物体を識別するための方法であって、前記方法は、
    前記二重エネルギーCTスキャン画像をフィルタリングし、グループ化されることになるピクセルを得ることと、
    前記二重エネルギーCTスキャン画像内の前記グループ化されることになるピクセルの位置に基づいて、前記グループ化されることになるピクセルを複数のピクセルグループにグループ化することと、
    各ピクセルグループ内の前記ピクセルに関する物質分解を実行することと、
    物質分解の結果に基づいて、各ピクセルグループに対応する物体を特定することとを含む、二重エネルギーCTスキャン画像内の物体を識別するための方法。
  12. フィルタリングする前記ステップは、
    前記ピクセルのCT値が基準CT値より大きいか否かを判断することと、
    そのCT値が前記基準CT値より大きい前記ピクセルを前記グループ化されることになるピクセルと判断することとを含む、請求項11記載の方法。
  13. 前記グループ化されることになるピクセルをグループ化する前記ステップは、
    前記グループ化されることになるピクセル内で隣接して位置するピクセルを同じピクセルグループにグループ化することとを含む、請求項11記載の方法。
  14. 各ピクセルグループに対応する前記物体を判断する前記ステップは、
    各ピクセルグループ内の前記ピクセルの分解結果に基づいて、前記ピクセルグループ内の各ピクセルの物質分解結果間の関係を指示する分解値散布図を確立することと、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの前記位置に基づいて前記ピクセルグループに対応する前記物体を特定することとを含む、請求項11記載の方法。
  15. 前記ピクセルグループに対応する前記物体を判断する前記ステップは、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と、異なる物体を表す特性曲線との間の平均距離を計算することと、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と、特定の物体を表す特性曲線との間の平均距離が、基準値以下である場合には、前記ピクセルグループを前記特定の物体に対応すると判断することとを含む、請求項14記載の方法。
  16. ピクセルグループに対応する前記物体を判断する前記ステップは、
    前記分解値散布図内の前記ピクセルグループの前記ピクセルの位置と、異なる物体を表す特性曲線との間の平均距離を計算することと、
    最小平均距離に対応する前記特性曲線によって表される前記物体を前記ピクセルグループに対応する前記物体と判断することとを含む、請求項14記載の方法。
  17. 前記ピクセルの2つの特性値は、それぞれ同等の密度値の2つの異なる基底物質であるか、又は前記二重エネルギーCT造影スキャン画像を得るために、前記二重エネルギーCTスキャンを実行するために用いられる2つの異なるエネルギーに対応するCT値である、請求項13記載の方法。
  18. 各ピクセルグループに対応する前記物体を判断した結果に基づいて、前記二重エネルギーCTスキャン画像内の前記ピクセルのCT値を変更することを更に含む、請求項11記載の方法。
JP2014246381A 2013-12-09 2014-12-05 二重エネルギーctスキャン画像における物体識別方法 Active JP6438752B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201310660019.4A CN104700389B (zh) 2013-12-09 2013-12-09 双能ct扫描图像中的对象识别方法
CN201310660019.4 2013-12-09

Publications (2)

Publication Number Publication Date
JP2015112481A true JP2015112481A (ja) 2015-06-22
JP6438752B2 JP6438752B2 (ja) 2018-12-19

Family

ID=53271694

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014246381A Active JP6438752B2 (ja) 2013-12-09 2014-12-05 二重エネルギーctスキャン画像における物体識別方法

Country Status (3)

Country Link
US (1) US10452948B2 (ja)
JP (1) JP6438752B2 (ja)
CN (1) CN104700389B (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2016158234A1 (ja) * 2015-03-30 2017-12-28 株式会社日立製作所 画像生成装置、画像生成方法及びx線ct装置

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10820865B2 (en) 2016-07-13 2020-11-03 Koninklijke Philiops N.V. Spectral computed tomography fingerprinting
CN106473761B (zh) * 2016-10-14 2019-04-02 山东大学 一种医学双能ct电子密度图像的重建与数值校准方法
US10769822B2 (en) * 2017-04-13 2020-09-08 Canon Medical Systems Corporation X-ray CT apparatus
CN107527057B (zh) * 2017-09-07 2020-03-31 国能日新科技股份有限公司 一种风速功率异常数据剔除方法及装置
CN109959669B (zh) 2017-12-26 2021-02-23 清华大学 基材料分解方法和装置及计算机可读存储介质
CN110547821B (zh) * 2018-05-31 2024-04-26 通用电气公司 追踪目标介入物的ct系统、方法及相关计算机可读介质
US11295487B2 (en) * 2018-09-07 2022-04-05 Canon Medical Systems Corporation X-ray CT apparatus, medical image processing apparatus, and X-ray CT system
WO2021194358A1 (en) * 2020-03-23 2021-09-30 Mars Bioimaging Limited Segmentation in multi-energy ct data
CN114916950B (zh) * 2022-07-21 2022-11-01 中国科学院深圳先进技术研究院 基于多层平板探测器的高空间分辨能谱ct图像重建方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040101086A1 (en) * 2002-11-27 2004-05-27 Sabol John Michael Method and apparatus for quantifying tissue fat content
JP2004174253A (ja) * 2002-11-27 2004-06-24 Ge Medical Systems Global Technology Co Llc 軟組織空間の視覚化の方法及び装置
US20070030944A1 (en) * 2005-08-08 2007-02-08 Michael Grasruck Method for an X-ray machine
JP2007111525A (ja) * 2005-10-17 2007-05-10 Siemens Ag X線コンピュータ断層撮影におけるコンピュータ断層撮影画像作成方法
JP2009142518A (ja) * 2007-12-17 2009-07-02 Ge Medical Systems Global Technology Co Llc X線ct装置
JP2010500901A (ja) * 2006-08-15 2010-01-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ エネルギー感受性コンピュータ断層撮影における動き補償
US20100328313A1 (en) * 2009-06-26 2010-12-30 Kabushiki Kaisha Toshiba Advanced clustering method for material separation in dual energy ct

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1622649A2 (en) * 2003-05-06 2006-02-08 Philips Intellectual Property & Standards GmbH Apparatus and method for examining an object by means of elastically scattered x-ray radiation and contrast agent
US20060285737A1 (en) * 2005-06-17 2006-12-21 Hamill James J Image-based artifact reduction in PET/CT imaging
CN1940992A (zh) * 2005-06-17 2007-04-04 Cti分子成像公司 Pet/ct成像中基于图像的伪影降低
US8837677B2 (en) * 2007-04-11 2014-09-16 The Invention Science Fund I Llc Method and system for compton scattered X-ray depth visualization, imaging, or information provider

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040101086A1 (en) * 2002-11-27 2004-05-27 Sabol John Michael Method and apparatus for quantifying tissue fat content
JP2004174261A (ja) * 2002-11-27 2004-06-24 Ge Medical Systems Global Technology Co Llc 組織脂肪含有量を定量化する方法及び装置
JP2004174253A (ja) * 2002-11-27 2004-06-24 Ge Medical Systems Global Technology Co Llc 軟組織空間の視覚化の方法及び装置
US20070030944A1 (en) * 2005-08-08 2007-02-08 Michael Grasruck Method for an X-ray machine
JP2007044520A (ja) * 2005-08-08 2007-02-22 Siemens Ag X線装置用の方法
JP2007111525A (ja) * 2005-10-17 2007-05-10 Siemens Ag X線コンピュータ断層撮影におけるコンピュータ断層撮影画像作成方法
JP2010500901A (ja) * 2006-08-15 2010-01-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ エネルギー感受性コンピュータ断層撮影における動き補償
JP2009142518A (ja) * 2007-12-17 2009-07-02 Ge Medical Systems Global Technology Co Llc X線ct装置
US20100328313A1 (en) * 2009-06-26 2010-12-30 Kabushiki Kaisha Toshiba Advanced clustering method for material separation in dual energy ct

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2016158234A1 (ja) * 2015-03-30 2017-12-28 株式会社日立製作所 画像生成装置、画像生成方法及びx線ct装置

Also Published As

Publication number Publication date
JP6438752B2 (ja) 2018-12-19
US20150161787A1 (en) 2015-06-11
US10452948B2 (en) 2019-10-22
CN104700389A (zh) 2015-06-10
CN104700389B (zh) 2019-08-13

Similar Documents

Publication Publication Date Title
JP6438752B2 (ja) 二重エネルギーctスキャン画像における物体識別方法
Heuvelmans et al. Optimisation of volume-doubling time cutoff for fast-growing lung nodules in CT lung cancer screening reduces false-positive referrals
Ko et al. Dual-energy computed tomography: concepts, performance, and thoracic applications
Jeong et al. Usefulness of a metal artifact reduction algorithm for orthopedic implants in abdominal CT: phantom and clinical study results
JP2019516460A (ja) 空間とスペクトル情報に基づく複数エネルギーのct画像におけるノイズ制御のためのシステムと方法
JP6548208B2 (ja) 二重エネルギーct造影剤強調されたスキャン画像内の石灰化部分を識別する装置の作動方法
Ebner et al. Lung nodule detection by microdose CT versus chest radiography (standard and dual-energy subtracted)
Wang et al. Optimal virtual monoenergetic image in “TwinBeam” dual‐energy CT for organs‐at‐risk delineation based on contrast‐noise‐ratio in head‐and‐neck radiotherapy
Zbijewski et al. Dual‐energy cone‐beam CT with a flat‐panel detector: Effect of reconstruction algorithm on material classification
CN102737375B (zh) 确定坏死组织的方法、图像处理装置和系统
Hamill et al. Respiratory‐gated CT as a tool for the simulation of breathing artifacts in PET and PET/CT
Altunbas et al. A unified scatter rejection and correction method for cone beam computed tomography
JP2018500061A (ja) 画像データ用のテクスチャ解析マップ
JP2016539704A (ja) 画像データからの骨の区分化
Li et al. Accuracy of dual‐energy computed tomography for the quantification of iodine in a soft tissue‐mimicking phantom
Kim et al. The optimal energy level of virtual monochromatic images from spectral CT for reducing beam-hardening artifacts due to contrast media in the thorax
Patel et al. Effects on image quality of a 2D antiscatter grid in x‐ray digital breast tomosynthesis: initial experience using the dual modality (x‐ray and molecular) breast tomosynthesis scanner
JP2012228404A (ja) 放射線画像処理方法および放射線画像処理装置
Srinivas-Rao et al. Dual-energy computed tomography to photon counting computed tomography: emerging technological innovations
Silva et al. Pre-surgical assessment of mandibular bone invasion from oral cancer: comparison between different imaging techniques and relevance of radiologist expertise
CN109389577B (zh) X射线图像处理方法和系统及计算机存储介质
Tamura et al. Feasibility of thin-slice abdominal CT in overweight patients using a vendor neutral image-based denoising algorithm: Assessment of image noise, contrast, and quality
JP2016137050A (ja) 医用画像処理装置、医用画像処理方法および医用画像診断装置
JP7346546B2 (ja) Ct仮想単色撮像のための自動的適応エネルギ設定を行う方法
Chaber et al. Hounsfield units and fractal dimension (test HUFRA) for determining PET positive/negative lymph nodes in pediatric Hodgkin’s lymphoma patients

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171205

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180328

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180403

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180703

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20181023

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20181119

R150 Certificate of patent or registration of utility model

Ref document number: 6438752

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250