JP2024513735A - サブトラクション撮像 - Google Patents
サブトラクション撮像 Download PDFInfo
- Publication number
- JP2024513735A JP2024513735A JP2023557761A JP2023557761A JP2024513735A JP 2024513735 A JP2024513735 A JP 2024513735A JP 2023557761 A JP2023557761 A JP 2023557761A JP 2023557761 A JP2023557761 A JP 2023557761A JP 2024513735 A JP2024513735 A JP 2024513735A
- Authority
- JP
- Japan
- Prior art keywords
- image
- data
- live
- ray
- image data
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 claims description 77
- 238000002601 radiography Methods 0.000 claims description 27
- 238000012937 correction Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 15
- 238000002347 injection Methods 0.000 claims description 14
- 239000007924 injection Substances 0.000 claims description 14
- 230000002792 vascular Effects 0.000 claims description 14
- 210000004204 blood vessel Anatomy 0.000 claims description 13
- 239000002872 contrast media Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 claims description 10
- 210000001519 tissue Anatomy 0.000 claims description 8
- 210000000988 bone and bone Anatomy 0.000 claims description 7
- 210000004872 soft tissue Anatomy 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 210000005166 vasculature Anatomy 0.000 claims description 3
- 239000003795 chemical substances by application Substances 0.000 claims 1
- 230000000873 masking effect Effects 0.000 abstract description 2
- 230000006978 adaptation Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 11
- 230000003044 adaptive effect Effects 0.000 description 9
- 238000002583 angiography Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 6
- 238000012512 characterization method Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000002083 X-ray spectrum Methods 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000002594 fluoroscopy Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000003936 working memory Effects 0.000 description 2
- 206010002329 Aneurysm Diseases 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000010102 embolization Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000926 neurological effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000013515 script Methods 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
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/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
- A61B6/504—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 for diagnosis of blood vessels, e.g. by angiography
-
- 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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- 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/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- 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/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/10116—X-ray image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
画像をマスキングする精度を更に改善するために、画像データ入力、データプロセッサ及び出力インタフェースを有する、デジタルサブトラクション撮像のための装置が提供される。複数の2D投影画像からの再構成に基づく対象の3D表現を有する対象の関心領域の3D画像データと、関心領域の2DライブX線画像とが受信され、位置合わせされ、2DライブX線画像に対応する3D画像データのマッチする姿勢が決定される。デジタル再構成X線撮影が、2Dマスク画像を生成するために、決定されたマッチする姿勢に基づいて3D画像データから計算される。2DライブX線画像に関連する現在のデータは、適合2Dマスク画像を得るために使用され、2Dライブ取得パラメータ及び/又は2DライブX線画像のデータを有する。生成された適合2Dマスクは、2DライブX線画像から減算された画像であり、関心領域の変化を強調するデジタル画像が提供される。
Description
本発明は、サブトラクション撮像に関し、特に、デジタルサブトラクション撮像のための装置、デジタルサブトラクション撮像のための撮像システム、及びデジタルサブトラクション撮像のための方法に関する。
X線撮像は、対象内の関心領域を視覚化するために使用される。しかしながら、X線に基づく撮像技術、例えば、蛍光透視法、コンピュータ断層撮影などは、血液及び隣接組織のX線減衰特性が非常に類似しているので、血管に関して制限される。血管系を視覚化するために、放射線不透過性造影剤が、血管に注入され、血管内のX線減衰を増加させることができる。血管の可視性を低減することができる周囲組織からの散乱なしに血管をより明確に視覚化するために、デジタルサブトラクション血管造影(DSA)画像は、同じ構造のコントラスト強調画像から造影前マスク画像を減算することによって計算される。しかしながら、患者の動きのために、マスク画像(典型的には注入されていない)及びライブ画像(典型的には造影剤で満たされている)のアライメントは、完全ではない。正味の効果は、通常動きアーチファクトと称されるものが存在する準最適減算である。US10682112B2は、3D術前ボリュームを用いた一連の3D X線透視画像における独立した動きの抑制に関し、X線画像から除去されるべき構造の(術前)3次元ボリュームの3D-2D位置合わせのデジタル再構成放射線画像が生成されることを記載している。これは、正しい観察面への構造体の投影、並びにその位置合わせを含む。デジタル再構成放射線画像は、X線画像から減算されて、抑制された構造の望ましくない/干渉する動きのない構造抑制透視画像を生成する。しかしながら、アーチファクトが依然として生じることができることが示されている。
したがって、画像をマスキングする精度を更に向上させる必要があり得る。
本発明の目的は、独立請求項の主題によって解決され、更なる実施形態は、従属請求項に組み込まれる。本発明の以下に説明する態様は、デジタルサブトラクション撮像のための装置、デジタルサブトラクション撮像のための撮像システム、及びデジタルサブトラクション撮像のための方法にも適用されることに留意されたい。
本発明によれば、デジタルサブトラクション撮像のための装置が提供される。装置は、画像データ入力部、データプロセッサ、及び出力インタフェースを有する。画像データ入力部は、対象の関心領域の3D画像データを受信するように構成され、3D画像データは、関心領域の複数の2D投影画像からの再構成に基づく対象の3D表現を有する。画像データ入力部は、また、関心領域の2DライブX線画像を受信するように構成される。更に、データプロセッサは、3D画像データと2DライブX線画像とを位置合わせするように構成され、2DライブX線画像に対応する3D画像データのマッチする姿勢が、決定される。データプロセッサは、また、決定されたマッチする姿勢に基づいて、3D画像データからデジタル再構成放射線撮影(digitally reconstructed radiography)を計算して、2Dマスク画像を生成するように構成される。データプロセッサは、また、2Dマスク画像の生成のために、適合された2Dマスク画像を得るために2DライブX線画像に関連する現在のデータを使用するように構成される。2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び2DライブX線画像のデータのグループのうちの少なくとも1つを有する。データプロセッサは、生成された適合2Dマスク画像を2DライブX線画像から減算するように更に構成される。更に、出力インタフェースは、関心領域の変化を強調するデジタル画像を提供するように構成される。
効果として、適合は、デジタル再構成放射線画像(DRR)画像及びライブ画像が異なるデータ及び形成プロセスから生じることを補償する。適合は、サブトラクション画像におけるアーチファクトを低減し、したがって、改善された精度でサブトラクション画像の形で現在の状況の視覚化を提供する。
一例によれば、画像データ入力部は、関心領域の血管造影注入2DライブX線画像を2DライブX線画像として受信するように構成される。データプロセッサは、生成された適合2Dマスク画像を、血管造影剤注入2DライブX線画像から減算するように構成される。出力インタフェースは、関心領域内の血管構造を表すデジタルサブトラクション血管造影を提供するように構成される。
一例によれば、適合された2Dマスク画像を得るために、画像データ入力部は、複数の2D投影画像を有する関心領域の投影データを受信するように構成される。画像データ入力部は、また、2DライブX線画像の取得から2Dデータ特性を受信するように構成される。更に、データプロセッサは、2Dデータ特性に基づいて、関心領域の投影データの再構成を適合させて、対象の適合3D表現を有する適合3D画像データを生成するように構成される。データプロセッサは、また、デジタル再構成放射線撮影の計算のために、適合された3D画像データを使用するように構成される。
オプションによれば、データプロセッサは、適合された3D画像データを位置合わせのためにも使用するように構成される。
一例によれば、適合された2Dマスク画像を得るために、画像データ入力部は、2DライブX線画像の取得から2Dデータ特性を受信するように構成される。データプロセッサは、2Dデータ特性に基づいて、3D画像データからデジタル再構成放射線撮影の計算を適合させて、適合された2Dマスク画像を生成するように構成される。
一例によれば、2Dデータ特性は、i)スペクトル、分解能、及びX線撮像エネルギのグループのうちの少なくとも1つを有する2Dライブ取得パラメータに基づく。加えて又は代わりに、2Dデータ特性は、ii)解像度及びコントラストのグループのうちの少なくとも1つを有する2DライブX線画像の画像データに基づく。
一例によれば、3D画像データについて、適合された再構成は、2DライブX線画像に対して使用されると予想される所定の平均光子エネルギに対応する再構成された減衰値を有する。
オプションによれば、データプロセッサは、各組織タイプが期待される光子エネルギに対して個別に調整された減衰値を有するように、骨及び軟組織タイプの領域にセグメント化された3D画像をセグメント化するように構成される。
別のオプションによれば、データプロセッサは、デジタル再構成放射線撮影のために異なるスペクトル関連因子を適用するように構成される。
更なるオプションによれば、3D画像データは、血管なしである。
一例によれば、データプロセッサは、2DライブX線画像内の対象の姿勢を推定するように構成される。3D画像データ及び2DライブX線画像の位置合わせのために、データプロセッサは、3D画像データの対象の3D表現を推定された姿勢にアラインするように構成される。加えて又は代わりに、関心領域の投影データの再構成のために、データプロセッサは、再構成のために推定された姿勢を考慮するように構成される。
一例によれば、データプロセッサは、入力として2DライブX線画像を使用する2D散乱推定又は3D散乱推定を提供するように構成される。データプロセッサは、計算されたデジタル再構成放射線撮影に散乱信号を加えるように構成される。
一例によれば、3D画像データの再構成のための補正として、データプロセッサは、i)仮想ヒール効果補正を提供するように構成される。加えて又は代わりに、データプロセッサは、ii)仮想逆検出器利得補正を提供するように構成される。
一例によれば、ディスプレイが、関心領域内の血管構造を表すデジタルを提示するために提供される。
本発明によれば、デジタルサブトラクション撮像のための撮像システムも、提供される。システムは、X線源及びX線検出器を有するX線撮像装置と、前述の例のうちの1つによるデジタルサブトラクション撮像のための装置とを有する。X線撮像装置は、2DライブX線画像を生成するように構成される。
本発明によれば、デジタルサブトラクション撮像のための方法も、提供される。方法は、以下のステップ:対象の関心領域の3D画像データを受信するステップであって、3D画像データは関心領域の複数の2D投影画像からの再構成に基づく対象の3D表現を有する、ステップと、関心領域の2DライブX線画像を受信するステップと、3D画像データと2DライブX線画像とを位置合わせするステップであって、2DライブX線画像に対応する3D画像データのマッチする姿勢が決定される、ステップと、決定されたマッチする姿勢に基づいて、3D画像データからデジタル再構成放射線撮影を計算して、2Dマスク画像を生成するステップとを有する。2Dマスク画像を生成するために、2DライブX線画像に関連する現在のデータが、適合された2Dマスク画像を得るのに使用され、2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び2DライブX線画像のデータのグループのうちの少なくとも1つを有する。更に、生成された適合2Dマスク画像を2DライブX線画像から減算するステップと、関心領域の変化を強調するデジタル画像を提供するステップとが、提供される。
応用の分野は、神経学、特に神経介入である。別の適切な分野は、脊椎関連の介入である。
一態様によれば、3D生成マスク画像は、ライブ画像に適合される。そのために、一例では、ライブ画像特性が、抽出され、生成された3Dボリュームをライブ画像に適合させるように再構成プロセスに導入される。代替的に又は追加的に、マッチングするデジタル再構成放射線画像は、ライブ画像の内容及び特性に適合される。この(潜在的に二重の)機構は、少なくともスペクトル、散乱、及び解像度適合を標的とする。
一態様によれば、適合の2つの主要なソースが、提供される。第1のオプションでは、1つの適合が、3D再構成時に、すなわち3D再構成中に提供され、第2のオプションでは、追加的に又は代替的に、1つの適合が、デジタル再構成放射線撮影生成時に、すなわち、デジタル再構成放射線撮影の計算中に提供される。一態様によれば、これらの2つのソースは、併せて使用される場合には結合されるが、独立して、例えば他方が存在しない状態で、使用されることもできる。
適合は、特に、例えば血管内処置、例えば狭窄検出、動脈瘤塞栓形成、神経学卒中治療などを含むがこれらに限定されないコントラスト強調撮像を使用する様々なX線ベースの処置のための、介入型X線撮像システム、例えば固定式及び可動式の撮像システムにおける使用に適している。
別の態様によれば、オプションとして、現在のライブ状況に対するマスク画像生成の適合は、特にロードマッピング又は装置視覚化における、他のサブトラクティブ撮像技法に対して提供される。例えば、血管系における装置の操縦中に、蛍光透視画像は、背景コンテンツを除去し、装置のみを可視のままにするために減算後に視覚化される。別の例では、以前の介入時間に取得された血管のフットプリントが、ガイダンス目的の装置操作中に重ね合わされる。撮像技術の2つのステップ(以前の瞬間対ライブ瞬間)の間に生じる任意の動きは、減算プロセスを劣化させることに留意されたい。したがって、適合は、これらの場合にも適している。
更なる態様によれば、画像データ入力部と、データプロセッサと、出力インタフェースとを有する、デジタルサブトラクション撮像のための装置が、提供される。複数の2D投影画像からの再構成に基づく対象の3D表現を有する対象の関心領域の3D画像データが、受信される。更に、関心領域の2DライブX線画像が、受信される。3D画像データ及び2DライブX線画像が、位置合わせされ、2DライブX線画像に対応する3D画像データのマッチする姿勢が、決定される。デジタル再構成放射線撮影が、2Dマスク画像を生成するために、決定されたマッチする姿勢に基づいて3D画像データから計算される。2Dマスク画像について、2DライブX線画像に関連する現在のデータは、適合された2Dマスク画像を得るために使用され、2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び/又は2DライブX線画像のデータを有する。生成された適合2Dマスクは、2DライブX線画像から減算された画像であり、関心領域の変化を強調するデジタル画像が提供される。
本発明のこれら及び他の態様は、以下に記載される実施形態から明らかになり、それを参照して説明される。
本発明の例示的な実施形態は、以下の図面を参照して以下に説明される。
ここで、特定の実施形態が、添付の図面を参照して詳細に説明される。以下の説明では、同様の図面参照番号が、異なる図面においても、同様の要素に対して使用される。詳細な構成及び要素のような、本明細書で規定される事項は、例示的な実施形態の包括的な理解を助けるために提供される。また、周知の機能又は構成は、不必要な詳細で実施形態を不明瞭にするので、詳細には説明されない。更に、「のうちの少なくとも1つ」などの表現は、要素のリストに先行する場合、要素のリスト全体を修飾し、リストの個々の要素を修飾しない。
用語「被験者」は、個体とも称され得る。「被験者」は、更に、患者とも称され得るが、この用語は、任意の病気又は疾患が被験者に実際に存在するかどうかを示さないことに留意されたい。
図1は、デジタルサブトラクション撮像のための装置10の一例を概略的に示す。装置10は、画像データ入力部12、データプロセッサ14、及び出力インタフェース16を有する。
オプションとして、ディスプレイ18も、被験者の関心領域内の血管構造を表すデジタル画像を提示するために提供される。点線の矢印20は、ディスプレイ18のためのデータを提供するデータ接続を示す。
画像データ入力部12は、対象の関心領域の3D画像データ22を受信するように構成され、3D画像データは、関心領域の複数の2D投影画像からの再構成に基づく対象の3D表現を有する。画像データ入力部12は、また、関心領域の2DライブX線画像24を受信するように構成される。
データプロセッサ14は、3D画像データと2DライブX線画像とを位置合わせするように構成され、2DライブX線画像に対応する3D画像データのマッチする姿勢が、決定される。データプロセッサ14は、また、決定されたマッチする姿勢に基づいて、3D画像データからデジタル再構成放射線撮影を計算して、2Dマスク画像を生成するように構成される。データプロセッサ14は、2Dマスク画像の生成のために、適合された2Dマスク画像を得るように2DライブX線画像に関連する現在データを使用するように更に構成され、2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び2DライブX線画像のデータのグループのうちの少なくとも1つを有する。データプロセッサ14は、更に、生成された適合2Dマスク画像を2DライブX線画像から減算するように構成される。
出力16インタフェースは、矢印20によって示されるデジタル画像を提供するように構成され、デジタル画像は関心領域の変化を強調する。
周囲のフレーム26は、画像データ入力部12、データプロセッサ14、及び出力インタフェース16がハウジングのような共通の構造内に配置されることができるというオプションを示す。しかしながら、それらは、別個のコンポーネントとして提供されることもできる。
図2は、デジタルサブトラクション撮像のための撮像システム50の一例を示す。システム50は、X線源54及びX線検出器56を有するX線撮像装置52を有する。更に、前述及び以下の例のうちの1つによるデジタルサブトラクション撮像のための装置10の一例が、提供される。X線撮像装置52は、2DライブX線画像を生成するように構成される。
X線撮像装置52は、天井構造から移動可能に吊り下げられる支持構造58を有する。別のオプションでは、移動式X線撮像装置、又は床上に支持された固定式であるが移動可能なシステムが、提供される。図示の例では、X線源54及びX線検出器56を担持する可動Cアーム60が、設けられる。関心領域を有する対象62が、示され、対象は、調整可能な被験者支持テーブル64上に配置される。ベッドサイドコントローラ66も、示される。上記には、照明機器68及びディスプレイ機器70が、設けられる。
デジタルサブトラクション撮像のための装置10の例は、前景に示されるワークステーション又はコンソール72の文脈において提供される。コンソール72は、ディスプレイ74と、キーボード76と、マウス78とを有する。更に、制御コンソール80及びグラフィックタブレット82が、示される。
デジタルサブトラクション撮像のための撮像システムは、デジタルサブトラクション血管造影のための撮像システムとも称される。
オプションでは、X線撮像装置は、また、3D画像データの再構成のために複数の2D投影画像を生成するように構成される。
図3は、デジタルサブトラクション撮像のための方法100の一例の基本ステップを示す。方法100は、以下のステップを有する。
第1のステップ102において、対象の関心領域の3D画像データが、受信される。3D画像データは、関心領域の複数の2D投影画像からの再構成に基づく対象の3D表現を有する。
第2のステップ104では、関心領域の2DライブX線画像が、受信される。
第3のステップ106において、3D画像データ及び2DライブX線画像が、位置合わせされる。2DライブX線画像に対応する3D画像データのマッチする姿勢が、決定される。
第4のステップ108では、デジタル再構成放射線撮影が、2Dマスク画像を生成するように、決定されたマッチする姿勢に基づいて3D画像データから計算される。2Dマスク画像を生成するために、2DライブX線画像に関連する現在のデータが、適合された2Dマスク画像を得るのに使用され、2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び2DライブX線画像のデータのグループのうちの少なくとも1つを有する。
第5のステップ110において、生成された適合2Dマスク画像は、2DライブX線画像から減算される。
第6のステップ112では、関心領域の変化を強調するデジタル画像が、提供される。
第1のステップ102において、対象の関心領域の3D画像データが、受信される。3D画像データは、関心領域の複数の2D投影画像からの再構成に基づく対象の3D表現を有する。
第2のステップ104では、関心領域の2DライブX線画像が、受信される。
第3のステップ106において、3D画像データ及び2DライブX線画像が、位置合わせされる。2DライブX線画像に対応する3D画像データのマッチする姿勢が、決定される。
第4のステップ108では、デジタル再構成放射線撮影が、2Dマスク画像を生成するように、決定されたマッチする姿勢に基づいて3D画像データから計算される。2Dマスク画像を生成するために、2DライブX線画像に関連する現在のデータが、適合された2Dマスク画像を得るのに使用され、2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び2DライブX線画像のデータのグループのうちの少なくとも1つを有する。
第5のステップ110において、生成された適合2Dマスク画像は、2DライブX線画像から減算される。
第6のステップ112では、関心領域の変化を強調するデジタル画像が、提供される。
オプションでは、2DライブX線画像が、関心領域の血管造影剤注入2DライブX線画像として提供される。減算のために、生成された適合2Dマスク画像は、血管造影剤注入2DライブX線画像から減算される。デジタル画像を提供するために、関心領域内の血管構造を表すデジタルサブトラクション血管造影が、提供される。
デジタルサブトラクション撮像のための方法は、デジタルサブトラクション血管造影のための方法とも称されることができる。
この方法の一例では、複数の2D投影画像が、画像のスキャンを形成する画像のセットとして提供される。一例では、複数のスキャン、すなわち、各々複数の2D投影画像を有する画像の2以上のセットが、提供される。
計算されたデジタル再構成放射線撮影は、提供された姿勢の下で3Dボリュームを撮像するときにX線システムが生成する2D画像の近似である。
3D画像データは、3D再構成データとも称される。3D画像データは、2D投影画像、すなわちX線画像のような画像の生成のための投影データを提供するが、この場合、仮想検出器平面上への3D画像データの仮想投影を提供する。したがって、3D画像データは、投影データとも称される。
2DライブX線画像は、2Dライブデータとも称される。
3D画像データは、3D画像データセットとも称される。
3D画像データは、3D再構成データとして、すなわち関心領域を表す3Dデータセットとして提供される。
一例では、3D再構成データが、関心領域を表す3Dデータセットである。例えば、神経介入では、これは、造影剤注入なしで実行され、いわゆる投影データを生成する回転スキャンから再構成された、頭部の3D表現であることができる。この場合、3Dデータは、骨及び軟組織の頭部の解剖学的構造を表すが、血管部分を含まない(X線では不可視である)。造影剤注入なしでは、これらは、X線を用いて周囲の軟組織と区別されることができない。
別のシナリオでは、造影された血管を有する3Dスキャンが提供され、造影剤なしで2Dライブ(蛍光透視)データが提供される。これは、蛍光透視データ上に仮想血管ロードマップを作成するために使用され、更なる造影剤注入なしでナビゲーションを可能にする。
2DライブX線画像は、血管造影取得によって提供される、すなわち取得される。
一例では、2D投影スキャンが、回転Cアームスキャンの複数のスキャンとして提供される。別の例では、2D投影スキャンが、CTスキャンの複数のスキャンとして提供される。
一変形形態では、Cアームスキャンが、造影剤なしで提供される。
一例では、2DライブX線画像に関連する現在のデータが、2Dデータ特性を有する。一例では、ライブ画像又はその解像度を生成するために使用されるX線スペクトルなどの2Dデータの所定の必須の特性が、収集される。2Dデータ特性は、システムからの取得パラメータから、又はライブ画像コンテンツ自体から直接的に、導出されることができる。
一例では、2DライブX線画像が、関心領域の血管造影剤注入2DライブX線画像として提供される。
別の例では、2DライブX線画像が、造影なしで、装置が挿入された蛍光透視画像として提供される。例えば、減算は、次いで、血管内ナビゲーション中に装置を強調するのを助けることができる。
図1に戻って参照すると、詳細には図示されていない装置の例では、画像データ入力部12が、2DライブX線画像として関心領域の血管造影剤注入2DライブX線画像を受信するように構成される。データプロセッサ14は、生成された適合2Dマスク画像を、血管造影剤注入2DライブX線画像から減算するように構成される。出力インタフェース16は、関心領域内の血管構造を表すデジタルサブトラクション血管造影を提供するように構成される。
デジタルサブトラクション撮像のための装置10は、デジタルサブトラクション血管造影のための装置とも称されることができる。
デジタルサブトラクション撮像のための装置10は、デジタルサブトラクション血管造影のための装置とも称されることができる。
図1においてオプションとして示される、装置の一例では、適合された2Dマスク画像を得るために、画像データ入力部12が、複数の2D投影画像を有する関心領域の投影データを受信するように構成される。画像データ入力部12は、また、2DライブX線画像の取得から2Dデータ特性28を受信するように構成される。データプロセッサ14は、2Dデータ特性28に基づいて、関心領域の投影データの再構成を適合させて、対象の適合3D表現を有する適合3D画像データを生成するように構成される。データプロセッサ14は、デジタル再構成放射線撮影の計算のために、適合された3D画像データを使用するように更に構成される。
オプションでは、データプロセッサ14が、適合された3D画像データを位置合わせにも使用するように構成される。
一例では、適合3D再構成が、提供される。これは、2D-ライブ互換性のあるデジタル再構成放射線画像を送達することがより可能であるボリュームを生成するための特定の適合を可能にする。これは、適合されたデジタル再構成放射線撮影アイデアを通して、標的化された解剖学的構造の所与の体積との互換性をもたらす。
図1における追加又は代替のオプションとして示される、装置の一例では、適合された2Dマスク画像を得るために、画像データ入力部12が、2DライブX線画像の取得から2Dデータ特性30を受信するように構成される。データプロセッサ14は、2Dデータ特性30に基づいて、3D画像データからのデジタル再構成放射線撮影の計算を適合させて、適合された2Dマスク画像を生成するように構成される。
一例では、適合されたDRR生成のために、デジタル再構成放射線撮影は、ライブ画像コンテンツ及び特性に対して適合可能である。例えば、3D又は2D散乱推定器が、散乱信号をデジタル再構成放射線撮影に追加するのに使用される。2Dアプローチは、3D再構成画像の視野内にない散乱構造、例えば患者テーブルとも称される対象支持体を含み得るので、DRR(デジタル再構成放射線撮影)又は3D再構成ではなく、血管造影画像を入力として使用する。
デジタル再構成放射線撮影及び血管造影をより一貫させることができる更なる改善は、仮想ヒール効果及び仮想逆検出器ゲイン補正である。両方の効果のために、画像再構成パイプラインからの較正データが、使用される。しかしながら、一例では、デジタル再構成放射線撮影プロセスが、後処理としてスペクトル適合又は解像度回復をも含む。一例では、適合されたデジタル再構成放射線撮影計算が、ライブ画像とその特性との両方を入力として受信する。
オプションでは、適合された2Dマスク画像を得るために、関心領域の投影データの再構成及びデジタル再構成放射線撮影の計算が、2Dデータ特性に基づいて適合される。
詳細には示されていない装置の一例では、2Dデータ特性が、i)スペクトル、分解能、及びX線撮像エネルギのグループのうちの少なくとも1つを有する2Dライブ取得パラメータに基づく。代替的に、又は追加的に、2Dデータ特性は、ii)解像度及びコントラストのグループのうちの少なくとも1つを有する2DライブX線画像の画像データに基づく。
2Dデータ特性を使用する適合は、3Dデータのための(複数の2D画像の)データ取得と2Dライブデータのためのデータ取得との間の較正型アライメントを提供する。
適合された3D再構成の間、スペクトルなどの2Dライブデータの特性は、一例では3D/2D互換性を改善するために使用される。一例として、3D再構成の互換性が、3D減衰値の再構成中に特定のビーム硬化補正方法を使用することによって改善される。例えば、再構成された減衰値が、2Dライブデータ取得によって使用されることが予想される平均光子エネルギに対応する場合に、有益であり得る。この平均光子エネルギは、例えば、3Dデータ及び2Dデータを生成するために使用されるX線スペクトルが異なり得るので、3D再構成のために使用される回転取得中の平均光子エネルギとは異なり得る。
この目的のために、3D画像は、更に、例えば閾値又はモデルベースのセグメント化を使用することによって、骨及び軟組織領域にセグメント化され得、その結果、各組織タイプは、期待される光子エネルギについて個別に調整された減衰値を有する。効果として、投影されたデータとライブデータとの間のX線スペクトルの不一致が、対処される。
更なる例では、散乱補正(散乱効果が再構成時に除去される)及び空間分解能の不一致など、3D再構成の他の適合が、提供される。
一例では、二重エネルギ又はスペクトル3D取得が、提供され、特定のエネルギ及び/又は材料タイプに対する1つ又は複数の再構成が、作成され、最終的な適合3D画像に結合される。
一例では、投影データが、理想的な3D再構成解像度が2Dライブデータにマッチすることを可能にしない場合、ニューラルネットワークベースの超解像方法が、3D画像解像度を更に改善するために使用される。
詳細には示されていない装置の例では、3D画像データのために、適合された再構成が、2DライブX線画像のために使用されると予想される所定の平均光子エネルギに対応する再構成された減衰値を有する。
更に、オプションとして、データプロセッサ14は、各組織タイプが期待される光子エネルギに対して個別に調整された減衰値を有するように、骨及び軟組織タイプの領域にセグメント化された3D画像をセグメント化するように構成される。
更に、別のオプションとして、代替的に又は追加的に、データプロセッサ14は、デジタル再構成放射線撮影のために異なるスペクトル関連因子を適用するように構成される。
更に、別のオプションとして、代替的に又は追加的に、3D画像データは、血管なしである。
その結果、2Dで仮想投影するために、3D画像データは、デジタル再構成放射線撮影を生成するようにデジタル再構成に対して最適化された、適合された3D画像データとして提供される。適合された3D画像データは、標準的な3D再構成と比較して、3Dビューイングの目的にはあまり適していないが、マスク画像を作成する観点から改善され得る。
平均光子エネルギは、典型的には、使用される回転取得中の平均光子エネルギとは異なる(3D及び2Dデータを生成するために使用されるX線スペクトルは通常異なるため)。
オプションでは、セグメンテーションのために、閾値ベース又はモデルベースのセグメンテーションが、提供される。
一例では、所定の平均光子エネルギが、ビーム硬化補正の代わりに提供される。
一例では、関心領域の投影データの再構成のために、再構成中のビーム硬化補正を含む適合された再構成が、提供される。3D画像データは、少なくとも2つの異なる吸収速度材料タイプに分離される。デジタル再構成放射線撮影のために、異なるスペクトル関連因子が、適用される。3D画像データは、血管なしである。「血管なし(vessel-less)」という用語は、血管構造の血管が画像に基づいて識別されることができない2D画像又は3D画像データを指す。したがって、3D画像は、識別可能な血管なしで提供される。
例えば、第1の吸収速度材料は、骨を表す部分に割り当てられ、第2の吸収速度材料は、組織を表す部分に割り当てられる。
一例では、3D画像データが、造影剤注入X線画像に基づく。コントラストを示す画像部分は、ボリューム生成のために除去される。
一例として、造影剤で満たされた血管、又は取り外し可能な医療装置、又は1つもしくはいくつかのECGケーブルのような外部オブジェクトは、それらが2Dライブデータに存在せず、特定の応用に対して関心がない場合、3Dデータから除去され得る。
ビーム硬化は、2D-ライブ互換性のあるデジタル再構成放射線画像を送達することがより可能であるボリュームを生成するために提供される特定の適合の一例である。
詳細には図示されない装置の一例では、データプロセッサ14が、2DライブX線画像内の対象の姿勢を推定するように構成される。3D画像データと2DライブX線画像との位置合わせのために、データプロセッサ14は、3D画像データの対象の3D表現を、推定された姿勢に対してアラインするように構成される。代替的に又は追加的に、関心領域の投影データの再構成のために、データプロセッサ14は、再構成に対して推定された姿勢を考慮するように構成される。
「姿勢」という用語は、撮像される被験者に対する撮像装置の空間的配置に関する。したがって、姿勢という用語は、空間的識別に関する。「姿勢」という用語は、被検体の姿勢、すなわち、閉じた又は開いた顎のような、可動身体部分又はセグメントの互いに対する配置、並びに互いに対する及び被検体に対するX線源及びX線検出器の幾何学的配置に関する。
一例では、2Dライブ画像に対応する3Dデータのマッチする姿勢が、3D再構成データ及びライブ画像から決定される。一例では、これは、3D-2D位置合わせの例として提供され、この動作が、一方又は他方のデータソース(通常、2Dライブデータ)内の注入された構造(血管)の存在を無視し得るという可能な制約を伴う。
姿勢推定のために、多数のコンピュータビジョン又は機械学習方法が、このタスクを解決するために提供される。一例では、再構成動作が、適合され、患者に埋め込まれた外部放射線不透過性マーカが、この動作をより速く又はよりロバストにするために使用される。例えば、マーカは、3D及び2Dライブデータの両方に存在し、3Dからの投影を、2Dにおけるそれらのライフ位置にマッチさせる幾何学的変換が、目標姿勢を決定する。
詳細には示されていない装置の例では、マッチする姿勢を見つけるために、データプロセッサ14が、i)システムベースの幾何学的位置合わせに基づく姿勢推定、ii)セグメント化ベースの姿勢推定、iii)造影注入血管構造の存在が2DライブX線画像及び/又は3D画像データにおいて無視されるという制約、並びにiv)3D画像データの異なる投影方向に対する複数のデジタル再構成放射線画像の生成、及び2Dライブ画像との可能なマッチングの識別に対する複数のデジタル再構成放射線画像の使用のグループのうちの少なくとも1つを提供するように構成される。
造影剤を注入された血管構造は、血管を表す。
一例では、マッチする姿勢に対して、放射線不透過性マーカが、被験者に取り付けられて提供され、これらのマーカは、3D画像データ及び2DライブX線画像の両方に存在する。3Dからの投影を2Dにおける位置にマッチさせる幾何学的変換は、目標姿勢を決定する。
一例では、外部放射線不透過性マーカが、被験者に一時的に埋め込まれる。
一例では、画像ベースの位置合わせが、提供される。
一例では、3D画像データ及び2Dライブ画像の位置合わせが、目標生体構造に依存して提供される。剛性解剖学的構造の場合、剛体変換が、提供される。非剛性構造の場合、弾性変形が、提供される。
一例では、3D画像データが、少なくとも2つの異なるサブボリュームとして提供される。サブボリュームは、互いに対する空間的配向において調整可能である。
例えば、被験者の顎運動を補償するために、3Dデータ内で関節運動させることができる、それぞれのサブボリュームが、識別される。
詳細には示されていない装置の例では、マスク生成のために、データプロセッサ14は、3D再構成プラス投影プロセスの伝達関数を反転するように構成される。言い換えれば、「3D再構成+投影」プロセスの伝達関数の逆関数が、提供される。これは、ライブ画像と可能な限り良好にマッチする再投影が提供されるという効果を有する。
例えば、放射測定と全空間スペクトルの両方が、適合される。
加えて、これは、空間的に変化する基準で行われるべきである。この目標を達成するために、ブラインド又は情報に基づくデコンボリューション/復元技術が、提供される。
詳細には示されない、装置の一例では、データプロセッサ14は、入力として2DライブX線画像を使用する2D散乱推定又は3D散乱推定を提供するように構成される。データプロセッサ14は、散乱信号を、計算されたデジタル再構成放射線撮影に加えるように構成される。
これは、3D再構成画像の視野内にない散乱構造、例えば患者テーブルを含み得るという効果を提供する。
これは、DRR及び血管造影をより一貫したものにすることができる更なる改善を提供する。
一例では、2D散乱推定が、入力として血管造影画像を使用する。
詳細には示されていない装置の例では、3D画像データの再構成のための補正として、データプロセッサ14は、i)仮想ヒール効果補正を提供するように構成される。加えて又は代わりに、データプロセッサ14は、ii)仮想逆検出器ゲイン補正を提供するように構成される。
ヒール効果補正及び検出器ゲイン補正は、仮想的に均一なX線強度を生成するために使用される。ヒール効果により、X線の放射強度は、放射角度ごとにわずかに異なる。検出器のX線検出効率は、内部検出器構成のような様々な要因により、画素ごとに変化する。2Dライブデータでは、これらの効果は、補正されてもされなくてもよい。更に、補正が適用される場合、補正は、実際の効果の単純化であり得、したがって、完全に正確ではないかもしれない。DRR生成のために、これらの補正ステップは、補正が2Dライブデータに対して使用されない場合、逆に適用されることができる。したがって、DRRでは、仮想不均一性が、2Dライブデータと同様に生成されることができる。代わりに、ヒール効果及び検出器ゲインは、システムの物理的現実をシミュレートするDRR生成中に考慮されることができ、次いで、2Dライブデータの簡略化された補正は、DRR上で作成される限られた精度を持つ同じ補正を達成するために適用される。
一例では、両方の効果について、画像再構成パイプラインからの較正データが、使用される。
一例では、デジタル再構成放射線撮影が、復元後に、ライブ画像と幾何学的にマッチしない場合、2D位置合わせプロセスが、提供され得る。一例として、面外運動の少なくとも一部が、適合ステップによって補正される。したがって、後続の2D位置合わせプロセスは、厳しい患者の動きが存在する場合であっても、良好な結果をもたらすことができる。次いで、減算が適用され、クリーンかつ改善されたデジタルサブトラクション血管造影法を生成することができる。
図3に戻って参照すると、詳細には示されていない方法の例において、適合された2Dマスク画像を得るために、以下のステップ、すなわち、
複数の2D投影画像を有する関心領域の投影データを受信するステップと、
2DライブX線画像の取得から2Dデータ特性を受信するステップと、
2Dデータ特性に基づいて、関心領域の投影データの再構成を適合させて、対象の適合3D表現を有する適合3D画像データを生成するステップと、
が、提供される。適合された3D画像データは、デジタル再構成放射線撮影の計算のために使用される。
複数の2D投影画像を有する関心領域の投影データを受信するステップと、
2DライブX線画像の取得から2Dデータ特性を受信するステップと、
2Dデータ特性に基づいて、関心領域の投影データの再構成を適合させて、対象の適合3D表現を有する適合3D画像データを生成するステップと、
が、提供される。適合された3D画像データは、デジタル再構成放射線撮影の計算のために使用される。
オプションでは、適合された3D画像データが、位置合わせにも使用される。
本方法の一例では、適合された2Dマスク画像を得るために、以下のステップ、すなわち、
2DライブX線画像の取得から2Dデータ特性を受信するステップと、
2Dデータ特性に基づいて、3D画像データからのデジタル再構成放射線撮影の計算を適合させて、適合された2Dマスク画像を生成するステップと、
が、提供される。
2DライブX線画像の取得から2Dデータ特性を受信するステップと、
2Dデータ特性に基づいて、3D画像データからのデジタル再構成放射線撮影の計算を適合させて、適合された2Dマスク画像を生成するステップと、
が、提供される。
一例では、2Dデータ特性は、i)スペクトル、解像度、及びX線撮像エネルギのグループのうちの少なくとも1つを有する2Dライブ取得パラメータに基づく。加えて又は代わりに、2Dデータ特性は、ii)解像度及びコントラストのグループのうちの少なくとも1つを有する2DライブX線画像の画像データに基づく。
本方法の一例では、3D画像データについて、適合された再構成は、2DライブX線画像に使用されると予想される所定の平均光子エネルギに対応する再構成された減衰値を有する。
本方法のオプションでは、3D画像は、各組織タイプが期待される光子エネルギに対して個別に調整された減衰値を有するように、骨及び軟組織タイプの領域にセグメント化される。別のオプションでは、デジタル再構成放射線撮影のために、異なるスペクトル関連因子が、適用される。更なるオプションでは、3D画像データは、血管なしである。
本方法の一例では、2DライブX線画像内の対象の姿勢を推定するステップが、提供される。更に、i)3D画像データと2DライブX線画像との位置合わせのために、3D画像データの対象の3D表現が、推定された姿勢にアラインされる。加えて又は代わりに、ii)関心領域の投影データの再構成のために、推定された姿勢が、再構成のために考慮される。
方法の一例では、マッチする姿勢を見つけるために、
i)システムベースの幾何学的位置合わせに基づく姿勢推定、
ii)セグメント化ベースの姿勢推定、
iii)造影剤を注入された血管構造の存在が、2DライブX線画像及び/又は3D画像データにおいて無視されるという制約、及び
iv)3D画像データの異なる投影方向のための複数のデジタル再構成放射線画像を生成し、2Dライブ画像との可能なマッチングを識別するために複数のデジタル再構成放射線画像を使用すること、
のグループの少なくとも1つが、提供される。
本方法の一例では、マスク生成のために、3D再構成プラス投影プロセスの伝達関数が、反転される。
i)システムベースの幾何学的位置合わせに基づく姿勢推定、
ii)セグメント化ベースの姿勢推定、
iii)造影剤を注入された血管構造の存在が、2DライブX線画像及び/又は3D画像データにおいて無視されるという制約、及び
iv)3D画像データの異なる投影方向のための複数のデジタル再構成放射線画像を生成し、2Dライブ画像との可能なマッチングを識別するために複数のデジタル再構成放射線画像を使用すること、
のグループの少なくとも1つが、提供される。
本方法の一例では、マスク生成のために、3D再構成プラス投影プロセスの伝達関数が、反転される。
本方法の一例では、2D散乱推定が、入力として2DライブX線画像、例えば血管造影画像を使用して提供される。別の例では、3D散乱推定が、提供される。散乱信号は、計算されたデジタル再構成放射線撮影に加えられる。
本方法の一例では、3D画像データの再構成のための補正として、i)仮想ヒール効果補正、及び/又はii)仮想逆検出器ゲイン補正が、提供される。
図4は、デジタルサブトラクション撮像のためのワークフロー200の一例を示す。言わば開始として、2Dライブデータ202及び3D再構成データ204が、提供される。2Dライブデータ202は、姿勢推定プロシージャ208のために使用される(206)。3D再構成データ204は、また、姿勢推定プロシージャ208のために使用される(210)。
オプションとして、破線で示されるように、投影データ212が、提供され、適合3D再構成プロシージャ216に対して使用され(214)、適合3D再構成データ204’に至る(218)。
姿勢推定プロシージャ208は、3D再構成データ204又は適合3D再構成データ204’からの入力をも受信する(224)DRR生成222に接続される(220)。DRR生成222は、伝達関数反転プロシージャ228に接続される(226)。
オプションとして、破線で示されるように、2Dライブデータ202も、伝達関数反転プロシージャ228に供給される(230)。
伝達関数反転プロシージャ228は、可能な2D位置合わせを有する減算プロシージャ234に接続され(232)、このプロシージャのために、2Dライブデータ202も、提供される(236)。可能な2D位置合わせを有する減算プロシージャ234の結果として、例えばサブトラクション画像、すなわち、適合されたサブトラクション画像として表示される、結果240が提供される(238)。
図5は、デジタルサブトラクション撮像のためのワークフロー300の別の例を示す。言わば開始として、2Dライブデータ302、2Dライブ取得パラメータ304、及び投影データ306が、提供される。投影データ306は、適合3D再構成プロシージャ310において使用される(308)。2Dライブ取得パラメータ304は、スペクトル、解像度などのデータを提供する2Dデータ特性決定プロシージャ314に提供され(312)、2Dデータ特性決定プロシージャ314は、また、2Dライブデータ302を受信する(316)。2Dデータ特性決定プロシージャ314は、適合3D再構成プロシージャ310に接続される(318)。適合3D再構成プロシージャ310は、姿勢推定プロシージャ326に供給される(324)3D再構成データ322に導く(320)。姿勢推定プロシージャ326は、また、2Dライブデータ302を受信する(328)。
オプションとして、破線によって示されるように、姿勢推定プロシージャは、適合3D再構成プロシージャ310に戻るように接続され(330)、姿勢を決定する際の最適化のためのループ状プロシージャを提供する。姿勢推定326から3D再構成310に向かう接続矢印330は、ライブデータに対する3D再構成の適合が、姿勢が知られると精緻化されることができることを示す。その場合、ライブデータ特性だけでなく、実際のライブデータコンテンツも、適合プロセスを駆動することができる。一例では、反復再構成方法が、追加の入力として位置合わせされたライブデータとともに使用される。再構成パラメータは、位置合わせされたライブデータに対するDRRの類似性を最大化するように反復的に最適化される。例えば、ライブデータの選択された領域における二乗平均平方根差が、この類似性を測定するように使用される。
姿勢推定プロシージャ326は、更に、2Dライブデータ302を受信する(336)適合DRRプロシージャ334に接続され(332)、また、2Dデータ特性決定プロシージャ314にも接続される(338)。適合されたDRRプロシージャ334は、減算+2D位置合わせプロシージャ344に接続され(340)、これは、また、2Dライブデータ302を受信する(342)。減算+2D位置合わせプロシージャ344の結果として、例えばサブトラクション画像、すなわち適合されたサブトラクション画像として表示される、結果348が、提供される(346)。
なお、図4及び図5において、直角の角の四角形は、データを表し、丸い角の四角形は、プロシージャ、すなわち、処理ステップ又はデータ収集/提供ステップを表す。矢印は、出力データを生成するための入力データに対する処理ステップの適用を表す。破線の矢印は、出力データを生成するための入力データに対する処理ステップのオプションの適用を表す。
図5において、フレーム310、314及び334は、基本的に異なる適合オプションを提供していることに更に留意されたい。フレーム306及び322は、3Dコンテキストに関連し、フレーム304及び302は、ライブ2Dコンテキストに関連する。
別の例では、プロセッサによって実行されると、プロセッサが前述の例のうちの1つの方法を、例えば適切なシステム上で、実行することを可能にする、コンピュータ可読コード又は命令を有するコンピュータプログラムが、提供される。
一例では、上記の例のうちの1つによる装置を制御するためのコンピュータプログラム又はプログラム要素が、提供され、このプログラム又はプログラム要素は、処理ユニットによって実行されるとき、上記の方法例のうちの1つの方法ステップを実行するように構成される。
したがって、コンピュータプログラム要素は、コンピュータユニットに記憶されてもよく、又は本発明の実施形態の一部であってもよい1つより多いコンピュータユニットにわたって分散されてもよい。この計算ユニットは、上述の方法のステップを実行する又はその実行を誘導するように構成されてもよい。更に、それは、上述の装置の構成要素を動作させるように構成されてもよい。計算ユニットは、自動的に動作するように、及び/又はユーザのオーダを実行するように、構成されることができる。コンピュータプログラムは、データプロセッサのワーキングメモリにロードされてもよい。従って、データプロセッサは、本発明の方法を実行するように装備されてもよい。
請求されるべき別の例では、コンピュータ可読媒体が前述の例のコンピュータプログラムを記憶した。
本発明の態様は、コンピュータによって実行され得るコンピュータ可読記憶装置上に記憶されたコンピュータプログラム命令の集合であってもよい、コンピュータプログラム製品において実施されてもよい。本発明の命令は、スクリプト、解釈可能プログラム、ダイナミックリンクライブラリ(DLL)又はJavaクラスを含むが、これらに限定されない、任意の解釈可能又は実行可能コードメカニズムであってもよい。命令は、完全な実行可能プログラム、部分的な実行可能プログラム、既存のプログラムに対する修正(例えば、更新)又は既存のプログラムに対する拡張(例えば、プラグイン)として提供されることができる。更に、本発明の処理の一部は、複数のコンピュータ又はプロセッサに分散されてもよい。
上述のように、処理ユニット、例えばコントローラは、制御方法を実施する。このコントローラは、必要とされる様々な機能を実行するために、ソフトウェア及び/又はハードウェアを用いて、様々な方法で実施されることができる。プロセッサは、必要とされる機能を実行するために、ソフトウェア(例えば、マイクロコード)を使用してプログラムされ得る1つ以上のマイクロプロセッサを採用するコントローラの一例である。しかしながら、コントローラは、プロセッサ使用の有無にかかわらず実施されてもよく、幾つかの機能を実行するための専用ハードウェアと、他の機能を実行するためのプロセッサ(例えば、1つ以上のプログラムされるマイクロプロセッサ及び関連する回路)との組み合わせとして実施されてもよい。
本開示の様々な実施形態に用いられ得るコントローラ構成要素の例は、従来のマイクロプロセッサ、特定用途向け集積回路(ASIC)及びフィールドプログラマブルゲートアレイ(FPGA)を含むが、これらに限定されない。
本発明のこの例示的な実施形態は、最初から本発明を使用するコンピュータプログラムと、アップデートの手段によって既存のプログラムを本発明を使用するプログラムに変えるコンピュータプログラムとの両方を包含する。
更に、コンピュータプログラム要素は、上述の方法の例示的な実施形態のプロシージャを満たすために必要なすべてのステップを提供することができてもよい。
本発明の更なる例示的な実施形態によれば、CD-ROMなどのコンピュータ可読媒体が、提示され、コンピュータ可読媒体は、それに記憶されたコンピュータプログラム要素を有し、そのコンピュータプログラム要素は、前述のセクションによって説明される。コンピュータプログラムは、他のハードウェアと共に又は他のハードウェアの一部として供給される光記憶媒体又はソリッドステート媒体などの適切な媒体上に記憶及び/又は配布され得るが、インターネット又は他の有線もしくは無線電気通信システムを介してなど、他の形態で配布されてもよい。
しかしながら、コンピュータプログラムは、ワールドワイドウェブのようなネットワーク上で提示されてもよく、そのようなネットワークからデータプロセッサのワーキングメモリにダウンロードされることができる。本発明の更なる例示的な実施形態によれば、コンピュータプログラム要素をダウンロードのために利用可能にするための媒体が、提供され、このコンピュータプログラム要素は、本発明の前述の実施形態のうちの1つによる方法を実行するように構成される。
本発明の実施形態は、それぞれ異なる主題を参照して説明されることに留意されたい。特に、いくつかの実施形態は、方法タイプの請求項を参照して説明され、他の実施形態は、装置タイプの請求項を参照して説明される。しかしながら、当業者は、上記及び下記の説明から、別段の通知がない限り、1つのタイプの主題に属する特徴の任意の組み合わせに加えて、異なる主題に関する特徴間の任意の組み合わせも、本出願で開示されると見なされることを理解するのであろう。しかしながら、全ての特徴は、特徴の単純な合計よりも高い相乗効果を提供するように組み合わせられることができる。
本発明は、図面及び前述の説明において図示され、詳細に説明されているが、そのような図示及び説明は、説明的又は例示的であり、限定的ではないと見なされるべきである。本発明は、開示された実施形態に限定されない。開示された実施形態に対する他の変形は、図面、開示及び従属請求項の検討から、請求項に記載の発明を実施する際に当業者によって理解及び達成されることができる。
請求項において、単語「有する」は、他の要素又はステップを排除するものではなく、不定冠詞「a」又は「an」は、複数性を排除するものではない。単一のプロセッサ又は他のユニットは、特許請求の範囲に記載されたいくつかのアイテムの機能を満たしてもよい。特定の手段が相互に異なる従属請求項において言及されているという単なる事実は、これらの手段の組み合わせが有利に使用されることができないことを示すものではない。請求項におけるいかなる参照符号も、範囲を限定するものとして解釈されるべきではない。
Claims (15)
- 画像データ入力部と、
データプロセッサと、
出力インタフェースと、
を有する、デジタルサブトラクション撮像のための装置において、
前記画像データ入力部は、対象の関心領域の3D画像データを受信するように構成され、前記3D画像データは、前記関心領域の複数の2D投影画像からの再構成に基づく前記対象の3D表現を含み、前記画像データ入力部は、前記関心領域の2DライブX線画像を受信するように構成され、
前記データプロセッサは、前記3D画像データと前記2DライブX線画像とを位置合わせするように構成され、前記2DライブX線画像に対応する前記3D画像データのマッチする姿勢が、決定され、前記データプロセッサは、前記決定されたマッチする姿勢に基づいて前記3D画像データからデジタル再構成放射線撮影を計算して、2Dマスク画像を生成し、前記2Dマスク画像の前記生成のために、適合2Dマスク画像を得るために前記2DライブX線画像に関連する現在のデータを使用するように構成され、前記2DライブX線画像に関連する前記データが、2Dライブ取得パラメータ及び前記2DライブX線画像のデータのグループのうちの少なくとも1つを有し、前記データプロセッサは、前記生成された適合2Dマスク画像を前記2DライブX線画像から減算するように構成され、
前記出力インタフェースは、前記関心領域における変化を強調するデジタル画像を提供するように構成される、
装置。 - 前記画像データ入力部は、前記関心領域の血管造影剤注入2DライブX線画像を前記2DライブX線画像として受信するように構成され、
前記データプロセッサは、前記生成された適合2Dマスク画像を、前記血管造影剤注入2DライブX線画像から減算するように構成され、
前記出力インタフェースは、前記関心領域内の血管構造を表すデジタルサブトラクション血管造影を提供するように構成される、
請求項1に記載の装置。 - 前記適合2Dマスク画像を得るために、
前記画像データ入力部は、前記複数の2D投影画像を有する前記関心領域の前記投影データを受信し、前記2DライブX線画像の取得から2Dデータ特性を受信するように構成され、
前記データプロセッサは、前記2Dデータ特性に基づいて、前記関心領域の前記投影データの前記再構成を適合させて、前記対象の適合された3D表現を有する適合された3D画像データを生成するように構成され、
前記データプロセッサは、前記デジタル再構成放射線撮影の計算のために、前記適合された3D画像データを使用するように構成され、
前記データプロセッサは、前記適合された3D画像データを前記位置合わせのためにも使用するように構成される、
請求項1又は2に記載の装置。 - 前記適合2Dマスク画像を得るために、
前記画像データ入力部は、前記2DライブX線画像の取得から2Dデータ特性を受信するように構成され、
前記データプロセッサは、前記2Dデータ特性に基づいて、前記3D画像データからの前記デジタル再構成放射線撮影の前記計算を適合させて、適合2Dマスク画像を生成するように構成される、
請求項1乃至3のいずれか一項に記載の装置。 - 前記2Dデータ特性は、
i)スペクトル、解像度、X線撮像エネルギのグループの少なくとも1つを有する2Dライブ取得パラメータ、及び/又は
ii)解像度及びコントラストのグループのうちの少なくとも1つを有する前記2DライブX線画像の画像データ、
に基づく、請求項1乃至4のいずれか一項に記載の装置。 - 前記3D画像データについて、前記適合された再構成は、前記2DライブX線画像に使用されると予想される所定の平均光子エネルギに対応する再構成された減衰値を有し、
前記データプロセッサは、各組織タイプが前記予想される光子エネルギに対して個別に調整された減衰値を有するように、骨及び軟組織タイプの領域にセグメント化された前記3D画像をセグメント化するように構成され、
前記データプロセッサは、前記デジタル再構成X線放射線撮影のために異なるスペクトル関連因子を適用するように構成され、
前記3D画像データは、血管なしである、
請求項1乃至5のいずれか一項に記載の装置。 - 前記データプロセッサは、前記2DライブX線画像内の前記対象の姿勢を推定するように構成され、
i)前記3D画像データ及び前記2DライブX線画像の位置合わせのために、前記データプロセッサは、前記3D画像データの前記対象の前記3D表現を前記推定された姿勢にアラインするように構成され、及び/又は
ii)前記関心領域の前記投影データの前記再構成のために、前記データプロセッサは、前記再構成に対して前記推定された姿勢を考慮するように構成される、
請求項1乃至6のいずれか一項に記載の装置。 - マッチする姿勢を見つけるために、前記データプロセッサは、
i)システムベースの幾何学的位置合わせに基づく姿勢推定と、
ii)セグメント化ベースの姿勢推定と、
iii)造影剤注入血管構造の存在が前記2DライブX線画像及び/又は前記3D画像データにおいて無視されるという制約と、
iv)前記3D画像データの異なる投影方向に対する複数のデジタル再構成放射線画像の生成、及び前記2Dライブ画像との可能なマッチングの識別のための前記複数のデジタル再構成放射線画像の使用と、
のグループのうちの少なくとも1つを提供するように構成される、
請求項1乃至7のいずれか一項に記載の装置。 - 前記マスク生成のために、前記データプロセッサは、前記3D再構成プラス投影プロセスの伝達関数を反転させるように構成される、請求項1乃至8のいずれか一項に記載の装置。
- 前記データプロセッサは、
前記2DライブX線画像を入力として用いた2D散乱推定、又は
3D散乱推定、
を提供するように構成され、
前記データプロセッサは、前記計算されたデジタル再構成X線画像に散乱信号を加えるように構成される、
請求項1乃至9のいずれか一項に記載の装置。 - 前記3D画像データの再構成のための補正として、前記データプロセッサは、
i)仮想ヒール効果補正、及び/又は
ii)仮想逆検出器ゲイン補正、
を提供するように構成される、
請求項1乃至10のいずれか一項に記載の装置。 - 前記装置が、
ディスプレイ、
を更に有し、
前記ディスプレイは、前記関心領域内の前記血管構造を表すデジタルを提示するために提供される、
請求項1乃至11のいずれか一項に記載の装置。 - X線源及びX線検出器を有するX線撮像装置と、
請求項1乃至12のいずれか一項に記載のデジタルサブトラクション撮像のための装置と、
有する、デジタルサブトラクション撮像のための撮像システムにおいて、
前記X線撮像装置は、前記2DライブX線画像を生成するように構成される、
撮像システム。 - デジタルサブトラクション撮像のための方法において、
対象の関心領域の3D画像データを受信するステップであって、前記3D画像データは、前記関心領域の複数の2D投影画像からの再構成に基づく前記対象の3D表現を有する、ステップと、
前記関心領域の2DライブX線画像を受信するステップと、
前記3D画像データと前記2DライブX線画像とを位置合わせするステップであって、前記2DライブX線画像に対応する前記3D画像データのマッチする姿勢が決定される、ステップと、
2Dマスク画像を生成するように、前記決定されたマッチする姿勢に基づいて、前記3D画像データからデジタル再構成放射線撮影を計算するステップであって、前記2Dマスク画像の生成のために、前記2DライブX線画像に関連する現在のデータが、適合2Dマスク画像を得るために使用され、前記2DライブX線画像に関連するデータは、2Dライブ取得パラメータ及び前記2DライブX線画像のデータのグループのうちの少なくとも1つを有する、ステップと、
前記生成された適合2Dマスク画像を前記2DライブX線画像から減算するステップと、
前記関心領域の変化を強調するデジタル画像を提供するステップと、
を有する方法。 - プロセッサによって実行されると、前記プロセッサが請求項14に記載の方法を実行することを可能にするコンピュータ可読コード又は命令を有する、コンピュータプログラム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP21290018.7A EP4062835A1 (en) | 2021-03-23 | 2021-03-23 | Subtraction imaging |
EP21290018.7 | 2021-03-23 | ||
PCT/EP2022/056447 WO2022200089A1 (en) | 2021-03-23 | 2022-03-14 | Subtraction imaging |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024513735A true JP2024513735A (ja) | 2024-03-27 |
Family
ID=75639831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023557761A Pending JP2024513735A (ja) | 2021-03-23 | 2022-03-14 | サブトラクション撮像 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20240164733A1 (ja) |
EP (2) | EP4062835A1 (ja) |
JP (1) | JP2024513735A (ja) |
CN (1) | CN117083021A (ja) |
WO (1) | WO2022200089A1 (ja) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7432924B2 (en) * | 2003-08-28 | 2008-10-07 | Kabushiki Kaisha Toshiba | 3D digital subtraction angiography image processing apparatus |
WO2015075047A1 (en) | 2013-11-20 | 2015-05-28 | Koninklijke Philips N.V. | Medical image viewing device for navigation in x-ray imaging, medical imaging system and method for providing improved x-ray image navigation information |
CN107106102B (zh) * | 2015-01-05 | 2020-11-17 | 皇家飞利浦有限公司 | 数字减影血管造影 |
-
2021
- 2021-03-23 EP EP21290018.7A patent/EP4062835A1/en not_active Withdrawn
-
2022
- 2022-03-14 US US18/282,595 patent/US20240164733A1/en active Pending
- 2022-03-14 WO PCT/EP2022/056447 patent/WO2022200089A1/en active Application Filing
- 2022-03-14 JP JP2023557761A patent/JP2024513735A/ja active Pending
- 2022-03-14 CN CN202280023744.8A patent/CN117083021A/zh active Pending
- 2022-03-14 EP EP22714186.8A patent/EP4312772A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
US20240164733A1 (en) | 2024-05-23 |
WO2022200089A1 (en) | 2022-09-29 |
EP4312772A1 (en) | 2024-02-07 |
EP4062835A1 (en) | 2022-09-28 |
CN117083021A (zh) | 2023-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10682112B2 (en) | Suppression of independent movements in a series of 2D X-ray fluoroscopy images using a 3D pre-operative volume | |
US7433507B2 (en) | Imaging chain for digital tomosynthesis on a flat panel detector | |
US11839501B2 (en) | Image creation device | |
Prell et al. | A novel forward projection-based metal artifact reduction method for flat-detector computed tomography | |
CN107019522B (zh) | 对操作者提供图像支持的x射线装置和数据载体 | |
US10083511B2 (en) | Angiographic roadmapping mask | |
JP7337556B2 (ja) | 医用画像処理装置、x線診断装置及び医用画像処理方法 | |
JP2009022754A (ja) | 放射線画像の位置揃えを補正する方法 | |
JP6704374B2 (ja) | 体脂肪率測定装置、方法およびプログラム | |
JP2004320771A (ja) | ディジタルサブトラクション血管造影法を実施するための方法 | |
CN109381205B (zh) | 用于执行数字减影血管造影的方法、混合成像装置 | |
JP2014014673A (ja) | 血管造影画像の取得方法 | |
US20230097849A1 (en) | Creation method of trained model, image generation method, and image processing device | |
CN106456080B (zh) | 用于修改x射线数据中的对tee探头的成像的设备 | |
JP7320116B2 (ja) | 対象を撮像する装置 | |
JP7167564B2 (ja) | X線撮影装置およびx線撮影装置の作動方法 | |
EP4287120A1 (en) | Guidance during medical procedures | |
US20240164733A1 (en) | Subtraction imaging | |
Close et al. | Accuracy assessment of layer decomposition using simulated angiographic image sequences | |
CN110678123A (zh) | 改进的x射线图像中的几何测量 | |
EP4162881A1 (en) | Enhancing angiograms | |
WO2022106305A1 (en) | Guidance for positioning a patient in medical imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230920 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20230920 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240402 |