JP5848759B2 - 低線量ctイメージングを実行するための方法及びシステム - Google Patents

低線量ctイメージングを実行するための方法及びシステム Download PDF

Info

Publication number
JP5848759B2
JP5848759B2 JP2013515994A JP2013515994A JP5848759B2 JP 5848759 B2 JP5848759 B2 JP 5848759B2 JP 2013515994 A JP2013515994 A JP 2013515994A JP 2013515994 A JP2013515994 A JP 2013515994A JP 5848759 B2 JP5848759 B2 JP 5848759B2
Authority
JP
Japan
Prior art keywords
image data
high resolution
resolution image
projection data
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.)
Expired - Fee Related
Application number
JP2013515994A
Other languages
English (en)
Other versions
JP2013529491A5 (ja
JP2013529491A (ja
Inventor
ラズ カルミ
ラズ カルミ
アミール リヴネ
アミール リヴネ
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2013529491A publication Critical patent/JP2013529491A/ja
Publication of JP2013529491A5 publication Critical patent/JP2013529491A5/ja
Application granted granted Critical
Publication of JP5848759B2 publication Critical patent/JP5848759B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • 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/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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/58Testing, adjusting or calibrating thereof
    • A61B6/582Calibration
    • A61B6/583Calibration using calibration phantoms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/436Limited angle

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Description

以下は一般にCTデータ収集及び再構成に、より具体的には高分解能再構成を伴う低線量CTに関する。
CTスキャナは生体組織に損傷を与え得る電離放射線を放出し、典型的な線量において癌、腫瘍及び遺伝子損傷のリスクの増加をもたらし、高線量において皮膚火傷及び脱毛を生じ得る。従って、CTスキャン中の電離放射線への患者暴露を軽減する(すなわち患者線量を減らす)様々な方法が提案されている。
文献で提案されている1つの方法は、圧縮センシング原理を用いるものであった。その目的は欠けている投影を事前画像などの追加情報で補い、一般的スパース制約(sparsity constraint)を導入することによって著しくアンダーサンプリングされたデータからアーチファクトのない断層画像を再構成することである。しかしながら、ほとんどの臨床例において、有益な情報は画像ドメインとサイノグラムドメインの両方に広範囲に広がっているため、CT画像は顕著なスパース特性を持たない。
結果として、圧縮センシング法を利用するために、標的画像と同様の特徴を持つ事前画像が必要である。その場合、2つの画像間の差分はさらに専用再構成中に利用されることができるスパース特性を持ち得る。事前画像は、例えばCT灌流などにおける標的スキャンの少し前(若しくは後)にとられるCTスキャンであり、又は心臓CTにおける全角度サンプリング低時間分解能画像であり得る。PICCS及びHYPRなどの技術はこうした事前スキャンに基づく。
あいにく、放射線量減少は必然的に画像ノイズに影響を及ぼし、これは主に検出器に到達するx線光子の固有ポアソン(すなわち"量子")ノイズによって支配される。加えて、一般的なCTシステムにおいて非常に低い線量で機能しようとする試みは、かなり過剰な画像ノイズとアーチファクトを作り出す。これは検出器素子によって作られる電子信号が電子ノイズのレベルに近い場合に発生する。
現在の臨床診療において、CTスキャナは要件が著しく異なり得る多くの異なる用途のために使用される。例えば、心臓スキャンは通常比較的短期間に高いx線束密度を要するが(高い管電流によって得られる)、肺スキャンは非常に低い管電流で行われ得る。このため、放射線検出器が非常に低い及び高いx線束密度の両方において信頼できるデータを与えることが重要である。
シンチレータピクセルに結合した電流積分フォトダイオードに基づく従来の積分検出器は、低信号を検出し同時に広ダイナミックレンジを持つ能力が限られている。通常はその場合、フォトダイオード暗電流と電子ノイズの両方によって影響されるノイズレベルは約10‐50平均x線量子に相当する。正確な数は特定設計と動作条件に依存する。信頼できる検出は測定値が約2倍以上ノイズよりも著しく大きい場合になされ得るので、ノイズレベルは最低検出限界を規定する。
従来の積分検出器は良好な直線性が通常非常に大きく1:100,000を超え得るフルダイナミックレンジを提供するが、実際の問題は主に単一読み取りあたり少数のx線量子、すなわち約1‐100程度のx線量子の信頼できる検出である。この検出範囲は、高減衰オブジェクトを通って横断する多くのビューがこのような低い値に達し得るため、低x線量において機能するために極めて重要である。低信号問題は特に小さいピクセルを持つ検出器アレイが高分解能スキャナを実現するために考慮される場合さらに頻発し得る。同様の制限は放射束が二重検出チャネル間で分割される二重エネルギー応用のために作られる二重層検出器において起こり得る。
本願の態様は上述の問題などに対処する。
一態様によれば、方法はアンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに基づいて高分解能画像データを生成するステップを含む。アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データは同じスキャンの異なる収集間隔中に収集される。
別の態様によれば、システムはスキャンの異なる積分期間中に高流束と低流束の間で放出放射束を交互に変調するように構成される線源と、流束の変調と協調して高分解能と低分解能の間で検出器ピクセル多重化を交互にスイッチするように構成される検出器アレイと、アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに対応する投影データに基づいて高分解能画像データを再構成するように構成される再構成器とを含む。
別の態様によれば、コンピュータ可読記憶媒体は、コンピュータのプロセッサによって実行されるときに、同じスキャンから得られるアンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに基づいて完全な高分解能画像データを再構成するために圧縮センシング再構成アルゴリズムをプロセッサに利用させる命令でエンコードされる。
本発明は様々な構成要素と構成要素の配置において、及び様々なステップとステップの配置において具体化し得る。図面は好適な実施形態を例示する目的に過ぎず、本発明を限定するものと解釈されない。
イメージングシステム例を図示する。 例示的な放射束変調/検出器ピクセル多重化レベル/グループ及びパターンを図示する。 例示的な放射束変調/検出器ピクセル多重化レベル/グループ及びパターンを図示する。 例示的な放射束変調/検出器ピクセル多重化レベル/グループ及びパターンを図示する。 例示的な放射束変調/検出器ピクセル多重化レベル/グループ及びパターンを図示する。 例示的な放射束変調/検出器ピクセル多重化レベル/グループ及びパターンを図示する。 アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データから完全な高分解能画像データを再構成するための方法例を図示する。 図7における完全な高分解能画像データを再構成するための方法例を図示する。
図1はコンピュータ断層撮影(CT)スキャナなどのイメージングシステム100を図示する。イメージングシステム100は固定ガントリ102と固定ガントリ102によって回転可能に支持される回転ガントリ104を含む。回転ガントリ104は長手方向すなわちz軸まわりに検査領域106の周囲を回転する。長椅子などの支持台118は検査領域106内で対象を支持する。支持台118はスキャンの前、最中及び/又は後にx、y、及び/又はz軸について対象を様々に位置づけるために使用されることができる。
x線管などの線源108は回転ガントリ104によって支持され、回転ガントリ104とともに回転し、放射線を放出する。線源コントローラ110は線源108を制御する。以下でより詳細に記載する通り、一実施形態において線源コントローラ110はスキャンの異なる積分間隔/期間中に少なくとも第1及び第2の異なる流束間で放出放射線の流束を変調するように線源108を制御することができる。第1(若しくは第2)の流束が第2(若しくは第1)の流束よりも大きい場合、スキャン中の第1及び第2の流束間の変調は線源108が高い流束を持つ放射線しか放出しない同じスキャンと比較して患者線量を減らす。
検出器ピクセルの単一の若しくは複数の列を持つ放射線検出器アレイ112が線源108の反対側に位置し、検査領域106を横断する放射線を検出し、それを示す投影データを生成する。検出器コントローラ114は放射線検出器アレイ112を制御する。以下でより詳細に記載する通り、一実施形態においてコントローラ110は個々のピクセル若しくはピクセルのより大きなグループが投影を検出するために使用されるように検出器ピクセル多重化を選択的に変更する。一般に、個々のピクセルはピクセルのより大きなグループと比較して比較的高分解能を提供する。
再構成器124は投影データを再構成し検査領域106を示すボリュメトリック画像データを生成する。再構成器124は例えば再構成アルゴリズムバンク116内のアルゴリズム及び/又は他のアルゴリズムなど、様々な再構成アルゴリズムを利用することができる。以下でより詳細に記載する通り、再構成器124はアンダーサンプリングされた高分解能投影データと低分解能再構成画像データが完全な高分解能ボリュメトリック画像データを再構成するために使用されるアルゴリズムを利用することができる。かかるアルゴリズムは患者線量を減らし完全な高分解能画像データを生成しながら、低放射束に伴う検出限界を緩和することを可能にする。
汎用コンピュータシステムがオペレータコンソール120として機能し、これはディスプレイ及び/又はプリンタなどの人が読み取り可能な出力装置とキーボード及び/又はマウスなどの入力装置を含む。コンソール120内のソフトウェアはオペレータがシステム100の操作を制御することを可能にし、例えばユーザが、放射線放出流束が変調され検出器ピクセル多重化がそれと協調して変化する(空間と時間においてレジストレーションされた高分解能及び低分解能投影データをもたらす)スキャン技術を選択し、得られる投影データから完全な高分解能画像データを再構成するための再構成アルゴリズムを選択することを可能にする。
上記で簡単に述べた通り、線源コントローラ110は放射束を変調するように線源108を制御することができ、検出器コントローラ114はピクセル多重化を変化させるように検出器アレイ112を制御することができる。これを行うために様々な方法が利用されることができ、方法は特定臨床応用、最適化、画像品質の妥協(例えば分解能、ノイズ、アーチファクトなどについて)、患者放射線量、システム能力及び性能、及び/又は他の要因など、様々な要因に基づき得る。
一例として、ピクセル多重化はより大きな有効ピクセル(通常は2若しくは4ピクセルのグループが使用される)への数検出器ピクセルのアナログ多重化を通して実現され得る。この場合、より大きなピクセルグループは小さな基本ピクセルとおおよそ同じ絶対電子ノイズレベルを持つが、同時に大きなピクセルグループに衝突する平均x線束は面積比に等しい倍率だけ大きい。従って、有効ピクセル面積の増加に対して信号対ノイズ比が改善される。空間分解能はより大きなピクセルを用いると減少する。
放射束変調は熱陰極などの電子放出体の温度を変えることによって、x線管の電子源と陽極の間の電場に影響を及ぼすようにパルス高電圧源でx線管に電力供給することによって、電子放出体の正面の電場を変化させることによって、x線管の陽極の表面に衝突する電子ビームの電気及び/又は磁気偏向を適用することによって、回転陽極の特殊な幾何学的構造を用いることによって、又は異なる材料から陽極を構成することなどによって、実現され得る。ある時間間隔において所望の平均放射束を実現するための方法は、'オン'及び'オフ'状態の間で放射線の非常に高速な連続スイッチングを用いることである。
図2、3、4及び5はそれぞれ互いに協調する変調放射束と多重化検出器ピクセルの限定されない実施例を図示する。図2(a)、3(a)、4(a)及び5(a)で、y軸は相対強度すなわち流束をあらわし、全図面でx軸は時間をあらわす。
図2(a)において流束は2つのレベル202と204の間を行ったり来たりし、低い方のレベル202は高い方のレベル204の1/4である。変調パターン206は流束が1積分期間(収集間隔、ビューなど)は高レベル204にあり、次の2積分期間は低レベル202にあるように流束を変調する。このパターンは経時的に繰り返される。図2(b)は対応する検出器多重化パターン208を示し、単一の小さな検出器ピクセル210が高レベル204の間に放射線を検出し、単一の小さな検出器ピクセルのグループ212が低レベル202の間に放射線を検出する。この実施例において、グループサイズは4検出器ピクセルであり、グループ形状はx軸(すなわちスキャナの角度方向)に沿った2つの検出器ピクセルとz軸に沿った2列の検出器ピクセルに及ぶ2次元アレイ(若しくはマトリクス)である。
図3(a)及び3(b)では、流束変調レベル202及び204並びに変調パターン206は図2(a)と同じである。しかしながら、多重化パターン302は高レベル204の間に放射線を検出するために単一検出器ピクセル210を使用し、低レベル202の間にz軸方向に沿った4検出器ピクセルのグループ304を使用することを含む。図4(a)及び4(b)では、低流束レベル402は高流束レベル204の半分であり、変調パターン404は1積分期間は高レベル204において、次の5積分期間は低レベル402において流束を変調し、多重化パターン406は高流束レベル204の間は単一検出器ピクセル210を、低流束レベル402の間はz軸方向に沿った2ピクセルのグループ408を使用することを含む。
図5(a)及び5(b)では、流束レベルは図1(a)及び2(a)と同じであり、流束変調パターンは図3(a)と同じであり、ピクセルグループは図2(b)と同じであり、ピクセル多重化パターン406は図4(b)と同じである。図2‐5において、総放射線量は高流束と単一ピクセルが各積分期間に使用されるスキャンの100%線量と比較して、それぞれ50.0%、50.0%、58.33%、及び37.5%に削減される。
上記実施例において、検出信号対ノイズ比は単一ピクセルが全ビュー/積分期間に使用される標準スキャンのものと等しく、相対放射線レベルは全ビュー/積分期間で1であることに留意されたい。他の実施形態において、多重化は信号対ノイズ比が少なくとも2つの積分期間で異なるように実行され得る。さらに、異なる及び/又はより大きなピクセルのグループ(例えば6、8、10、16など)が利用され得る。さらに、2つよりも多くの異なる放射強度レベル及び/又は2つよりも多くの異なるピクセルグループが利用されることができる。
さらに、2つの高分解能投影間の時間差はスキャン持続期間内に変化し得る。加えて、x線スペクトルはスキャン中に変化してもしなくてもよい。さらに、放射線の変調位相(若しくは基準時点に対する全シーケンスのシフト)が時間において調節可能であり得る。また対象形状において積分期間あたりのx線強度変調はステップ関数であり、例えばおおよそ10‐50マイクロ秒程度のいくらか遅い反応(図6の曲線602)も本明細書で考慮されることに留意されたい。
二重エネルギー応用のために作られる二重層検出器において、多重化期間は上方及び下方ピクセルを低ノイズで有効な従来の単層検出器ピクセルに組み合わせるために使用され得る。異なるスペクトル画像の再構成はアンダーサンプリングされた二重エネルギー投影と不完全な全スペクトル投影を利用する。
一例において、ピクセル多重化は相補型Nチャネル及びPチャネルCMOSトラジスタから作られる従来のCMOSスイッチに基づき得る。所要+Vc(制御電圧)をNチャネルゲートに、−VcをPチャネルゲートに印加することによって、スイッチがショート若しくは開接点として駆動され得る。N及びPトランジスタの構成を用いることでスイッチングシーケンス中に誘導される過剰電荷注入を削減することができる。結合する2つのトランジスタがよくマッチする場合、1fCbをはるかに下回るスイッチングシーケンス中の電荷注入が実現され得る。一部の検出器構成においてスイッチング中の注入電荷は無視できる可能性がある。
他の構成において、注入電荷は無視できず、従ってスイッチング後及び新たな読み取りが開始される前に直ちにこの電荷をリセットする特殊回路が組み入れられるべきである。これは例えばCT検出器電子機器において今日既に使用されている標準的技術によってなされることができる。スイッチング期間は数ナノ秒ほど低くなり得る。しかしながら、任意の追加リセット機構を含む総スイッチング時間は積分期間に従って設定され得る。例えばシステム100がおおよそ100‐300マイクロ秒の積分期間で構成される場合、スイッチング期間は数マイクロ秒程度に設定され得る。'積分期間'などという語は複数のイメージングビューの個々の持続時間を決定する対応する任意の一般的収集技術をあらわすために使用されることが留意されるべきである。
上記の通り、流束変調及びピクセル多重化が同期される。一例において、この同期は放射線変調若しくはピクセル多重化のいずれかの制御信号の時間位相を(例えばシステム較正スキーム中に)調節することによって制御されることができる。較正は臨床スキャンの前に、若しくは他の方法で、一回実行され得る。一例として、較正手順中、計画された交互スキャン構成が(空気中若しくはファントム上で)実行され、データが記録される。そして、相対的交代位相がわずかに変更され、測定が繰り返される。小ピクセルデータが全トライアル間で最高信号を実現し、多重化ピクセルグループのデータが最低信号を実現する位相設定を見つけるために反復シーケンスが実行され得る。
上記の通り、再構成器124はアンダーサンプリングされた高分解能投影データと、アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データから生成される低分解能再構成画像データに基づいて、完全な高分解能画像データを再構成することができる。一実施形態において、再構成アルゴリズムは式1と2を同時に解くことによって高分解能画像データを生成する。
Figure 0005848759
Figure 0005848759
式1は圧縮センシングのスパース性考察に関し、式2は断層画像再構成に関し、Ψはスパース化変換(sparsifying transform)であり、Xは高分解能画像データをあらわし、BはXの3D空間分解能をXのものに削減するブラーリング(blurring)変換であり、Xはよく再構成された基準低分解能画像データであり、Mは全関連スキャナ特性を含むシステム変換であり、Yはアンダーサンプリングされた高分解能投影データである。
一例において、式1は全変動法を用いてノルム最小化を介して処理されることができ、式2は最小二乗解という意味で反復再構成法(例えばART、MLEM)若しくはポアソンノイズモデルに基づく最適化によって処理されることができる。しかしながら、他の適切な数学的方法が代替的に使用され、本明細書で考慮され得る。変換Bは、フィルタ特性が高分解能及び低分解能モードの既知の変調伝達関数(MTF)から得られる、画像ボクセル空間内若しくはフーリエ変換空間内で操作される画像空間フィルタであり得る。適切なブラーリング変換は平滑化ローパスフィルタであり得る。ブラーリング変換Bは式1の反復解に組み込まれる。かかるプロセスに適切なスキームは式3に示される。
Figure 0005848759
tは反復シーケンスをあらわし、Iは更新された高分解能画像をあらわし、Iは基準画像をあらわし、Bはブラーリング変換をあらわし、αはプリセットパラメータをあらわし、TVは全変動作用素であり、'del'演算子(∇)は(各ボクセルについて)そのボクセル内の変化あたりの全変動の相対勾配を与える。
ブラーリング変換Bはシステム較正若しくはプリセットとして一度計算され得る。必要なパラメータが全てわかっているので、これは分析的に計算され得るか、又はスキャナのコンピュータモデルによってシミュレートされ得る。高分解能及び低分解能モードでスキャンされ得るファントムに対する測定較正手順を設計することもまた可能である。高分解能画像を低分解能画像へ変更する適切な変換を見つけるために画像特性が使用され得る。
図7はアンダーサンプリングされた高分解能投影データと不完全な低分解能投影データから高分解能画像データを再構成するための方法を例示する。
702においてスキャンが実行され、ここで放射束が変調され検出器ピクセルが協調して多重化される。限定されない実施例として、流束変調と検出器ピクセル多重化は図2‐5、その組み合わせ、及び/又はその他に関連して記載され得る。
704において、アンダーサンプリングされた高分解能投影データと低分解能投影データが組み合わされて低分解能投影データの完全なセットを生成する。一例において、複数の空間的に隣接した高分解能投影データが組み合わされて有効な低分解能投影データを生成し得る。
706において、低分解能投影データの完全なセットが再構成されて基準低分解能画像データを生成する。
708において、アンダーサンプリングされた高分解能投影データと基準低分解能画像データが利用されて完全な高分解能画像データを再構成する。本明細書で論じる通り、図8若しくはその他に記載の通り完全な高分解能画像データを再構成するために圧縮センシング再構成が利用され得る。
図8は図7の方法の動作708で使用され得る適切な圧縮センシングフロー図を図示する。
802において、アンダーサンプリングされた高分解能投影データと再構成された基準低分解能画像データが得られる。
804において、再構成された基準低分解能画像データが鮮明化される。例えば、一実施形態において再構成された基準低分解能画像データは画像再構成への初期推定を与えるデコンボリューション法を介して鮮明化される。
806において、アンダーサンプリングされた高分解能投影データと鮮明化された再構成された基準低分解能画像データが利用されて中間高分解能画像データを再構成する。再構成法は反復断層再構成であり得る。
808において、再構成された中間高分解能画像データがフィルタされる。例えば、一実施形態において、再構成された低分解能画像データが例えば数式1‐3に関連して本明細書で記載したブラーリング変換Bを用いて不鮮明化される。
810において、フィルタされた再構成された中間高分解能画像データと再構成された低分解能画像データの差分をとることによって差分画像データが生成される。
812において、差分画像データが所定基準を満たすかどうかが決定される。
差分画像データが所定基準を満たさない場合、814において中間高分解能画像データが最適化され、808におけるフィルタ後に再構成された低分解能画像により近くなる新たな中間高分解能画像を生成するために、動作808から812が繰り返される。最適化は再構成パラメータ、スパース性、全変動、正則化、及び/又は他の要因を考慮し得る。所定基準は所定数の反復でもあり得る。
差分画像データが所定基準を満たす場合、816において停止基準が満たされるかどうかが決定される。基準は所定数の反復、所定誤差閾値、反復結果間の差、及び/又は他の基準のうちの1つ以上を含み得る。
停止基準が満たされない場合、806の最初の動作に使用された鮮明化された再構成された基準低分解能画像データを交換する中間高分解能画像データを用いて動作806‐816が繰り返される。一般に、反復プロセスは継続し、各ステップにおいて高分解能画像データは所定解により近くなる。
停止基準が満たされる場合、818において高分解能画像データが出力される。
上記動作は、(複数の)コンピュータプロセッサによって実行されるときに(複数の)プロセッサに本明細書に記載の動作を実行させるコンピュータ可読命令を用いて実施され得る。かかる場合、命令は関連コンピュータに付随する及び/又は他の方法でアクセス可能なメモリなどのコンピュータ可読記憶媒体に記憶される。
"高い"、"より高い"、"低い"、及び"より低い"という語は本明細書において相対的レベルをあらわすために使用され、"高分解能"とは特定用途における目標分解能をあらわし、"低分解能"とはその用途における目標分解能よりも低い結果をあらわすことが留意される。
本発明は好適な実施形態に関して記載されている。修正及び変更は上記の詳細な説明を読んで理解することで想到され得る。本発明はかかる修正及び変更が添付の請求項若しくはその均等物の範囲内にある限り全て含むものと解釈されることが意図される。

Claims (20)

  1. アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに基づいて高分解能画像データを生成するステップを有する方法であって、
    前記アンダーサンプリングされた高分解能投影データと前記不完全な低分解能投影データは同じスキャンの異なる収集間隔中に収集される、方法。
  2. 前記不完全な低分解能投影データを前記アンダーサンプリングされた高分解能投影データで補完するステップをさらに有する、請求項1に記載の方法。
  3. 前記不完全な低分解能投影データに基づいて低分解能画像データを再構成するステップと、
    基準低分解能画像データと前記アンダーサンプリングされた高分解能投影データに基づいて前記高分解能画像データを再構成するステップとをさらに有する、請求項2に記載の方法。
  4. 前記基準低分解能画像データを鮮明化するステップと、
    前記鮮明化された基準低分解能画像データと前記アンダーサンプリングされた高分解能投影データに基づいて完全な高分解能画像データを再構成するステップとをさらに有する、請求項3に記載の方法。
  5. 前記鮮明化するステップが前記基準低分解能画像データをデコンボリューションするステップを含む、請求項4に記載の方法。
  6. 前記基準低分解能画像データと前記アンダーサンプリングされた高分解能投影データに基づいて中間高分解能画像データを再構成するステップと、
    前記高分解能画像データを生成するプロセスにおいて前記中間高分解能画像データをフィルタするステップをさらに有する、請求項3乃至5のいずれか一項に記載の方法。
  7. 前記中間高分解能画像データをフィルタするステップが前記中間高分解能画像データを不鮮明化するステップを含む、請求項6に記載の方法。
  8. 前記中間高分解能画像データが、決定された高分解能スキャンモードと低分解能スキャンモードに対応するMTFに基づいてフィルタされる、請求項6乃至7のいずれか一項に記載の方法。
  9. 前記フィルタされた中間高分解能画像データと前記基準低分解能画像データに基づいて差分画像データを決定するステップと、
    前記差分画像データが所定基準を満たすまで前記中間高分解能画像データを最適化するステップとをさらに有し、前記最適化された中間高分解能画像データが完全な高分解能画像データとして出力される、請求項6乃至8のいずれか一項に記載の方法。
  10. 前記アンダーサンプリングされた高分解能投影データと前記不完全な低分解能投影データが、高分解能データ収集と不完全な低分解能データ収集が交互に行われるイメージング手順中に収集される、請求項1乃至9のいずれか一項に記載の方法。
  11. 前記高分解能データ収集が、第1の流束を持つ放射線を放出するステップと、第1の領域を持つ検出器ピクセルを介して該放射線を検出するステップとを含み、前記低分解能データ収集が、第2の流束を持つ放射線を放出するステップと、第2の領域を持つように組み合わされた2つ以上の検出器ピクセルを介して該放射線を検出するステップとを含み、前記第1の流束は前記第2の流束よりも大きく、前記第1の領域は前記第2の領域よりも小さい、請求項10に記載の方法。
  12. 前記高分解能画像データを生成するために圧縮センシング再構成アルゴリズムを利用するステップをさらに有する、請求項1乃至11のいずれか一項に記載の方法。
  13. スキャンの異なる積分期間中に高流束と低流束の間で放出放射束を交互に変調する線源と、
    前記流束の変調と協調して高分解能と低分解能の間で検出器ピクセル多重化を交互にスイッチする検出器アレイと、
    アンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに対応する投影データに基づいて高分解能画像データを再構成する再構成器とを有する、システム。
  14. 前記再構成器が前記高分解能画像データを再構成するために圧縮センシング再構成アルゴリズムを利用する、請求項13に記載のシステム。
  15. 前記再構成器が前記低分解能投影データと前記アンダーサンプリングされた高分解能投影データに基づいて低分解能画像データを再構成し、前記低分解能画像データを鮮明化し、前記鮮明化された低分解能画像データに一部基づいて中間高分解能画像データを生成する、請求項14に記載のシステム。
  16. 前記再構成器が前記アンダーサンプリングされた高分解能投影データと前記鮮明化された低分解能画像データに基づいて中間高分解能画像データを再構成し、前記中間高分解能画像データを不鮮明化し、前記不鮮明化された中間高分解能画像データに一部基づいて前記高分解能画像データを生成する、請求項15に記載のシステム。
  17. 前記検出器アレイが検出器ピクセルの複数の列を含み、前記ピクセルの小さなグループが高分解能投影データを生成するために利用され、前記ピクセルの大きなグループが低分解能投影データを生成するために利用される、請求項13乃至16のいずれか一項に記載のシステム。
  18. 前記アンダーサンプリングされた高分解能投影データと前記不完全な低分解能投影データに基づいて低分解能投影データの完全なセットが生成され、前記高分解能画像データが該完全な低分解能投影データに一部基づいて生成される、請求項13乃至17のいずれか一項に記載のシステム。
  19. 高分解能及び低分解能収集の両方に対して所定の信号対ノイズ比が維持される、請求項13乃至18のいずれか一項に記載のシステム。
  20. コンピュータのプロセッサによって実行されるときに、該プロセッサに、
    同じスキャンから得られるアンダーサンプリングされた高分解能投影データと不完全な低分解能投影データに基づいて完全な高分解能画像データを再構成するために圧縮センシング再構成アルゴリズムを利用させる、コンピュータ実行可能命令でエンコードされるコンピュータ可読記憶媒体。
JP2013515994A 2010-06-21 2011-04-27 低線量ctイメージングを実行するための方法及びシステム Expired - Fee Related JP5848759B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US35682010P 2010-06-21 2010-06-21
US61/356,820 2010-06-21
PCT/IB2011/051850 WO2011161558A1 (en) 2010-06-21 2011-04-27 Method and system for performing low- dose ct imaging

Publications (3)

Publication Number Publication Date
JP2013529491A JP2013529491A (ja) 2013-07-22
JP2013529491A5 JP2013529491A5 (ja) 2015-07-09
JP5848759B2 true JP5848759B2 (ja) 2016-01-27

Family

ID=44120321

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013515994A Expired - Fee Related JP5848759B2 (ja) 2010-06-21 2011-04-27 低線量ctイメージングを実行するための方法及びシステム

Country Status (6)

Country Link
US (1) US9262845B2 (ja)
EP (1) EP2583250B1 (ja)
JP (1) JP5848759B2 (ja)
CN (1) CN102947864B (ja)
RU (1) RU2571564C2 (ja)
WO (1) WO2011161558A1 (ja)

Families Citing this family (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5608001B2 (ja) * 2010-07-13 2014-10-15 富士フイルム株式会社 放射線画像撮影方法および装置
JP5818588B2 (ja) * 2011-09-05 2015-11-18 株式会社東芝 放射線検出データ処理装置及び方法
US9147229B2 (en) * 2012-01-20 2015-09-29 Kabushiki Kaisha Toshiba Method and system for image denoising using discrete total variation (TV) minimization with one-direction condition
EP2915180B1 (en) * 2012-10-30 2018-12-05 California Institute of Technology Fourier ptychographic imaging systems, devices, and methods
US9864184B2 (en) 2012-10-30 2018-01-09 California Institute Of Technology Embedded pupil function recovery for fourier ptychographic imaging devices
US10652444B2 (en) 2012-10-30 2020-05-12 California Institute Of Technology Multiplexed Fourier ptychography imaging systems and methods
CN103413338B (zh) * 2013-05-29 2016-04-13 中国工程物理研究院流体物理研究所 一种基于广义变分最小化的少量投影ct图像重建方法
CN104240270B (zh) * 2013-06-14 2017-12-05 同方威视技术股份有限公司 Ct成像方法和系统
CN105659143B (zh) 2013-07-31 2019-03-22 加州理工学院 孔径扫描傅立叶重叠关联成像
KR20150016032A (ko) 2013-08-02 2015-02-11 삼성전자주식회사 영상 복원 모드 선택이 가능한 영상 복원 방법 및 그 장치
JP2016530567A (ja) 2013-08-22 2016-09-29 カリフォルニア インスティチュート オブ テクノロジー 可変照明フーリエタイコグラフィー撮像装置、システム、及び方法
JP6466057B2 (ja) * 2013-09-04 2019-02-06 キヤノンメディカルシステムズ株式会社 医用画像診断装置
DE102013217852B3 (de) * 2013-09-06 2014-10-30 Siemens Aktiengesellschaft Verfahren und Röntgensystem zur Zwei-Energiespektren-CT-Abtastung und Bildrekonstruktion
US11468557B2 (en) 2014-03-13 2022-10-11 California Institute Of Technology Free orientation fourier camera
JP6072723B2 (ja) 2014-04-21 2017-02-01 株式会社日立製作所 磁気共鳴イメージング装置、及び画像撮像方法
US10162161B2 (en) 2014-05-13 2018-12-25 California Institute Of Technology Ptychography imaging systems and methods with convex relaxation
WO2016007605A1 (en) * 2014-07-08 2016-01-14 The General Hospital Corporation System and method for motion-free computed tomography
EP3195265B1 (en) * 2014-09-15 2018-08-22 Koninklijke Philips N.V. Iterative image reconstruction with a sharpness driven regularization parameter
AU2015369663A1 (en) 2014-12-22 2017-05-11 California Institute Of Technology Epi-illumination fourier ptychographic imaging for thick samples
CN107209362B (zh) 2015-01-21 2020-11-06 加州理工学院 傅立叶重叠关联断层摄影
CN107209123B (zh) 2015-01-26 2020-08-11 加州理工学院 多孔傅立叶重叠关联和荧光成像
JP6675214B2 (ja) * 2015-03-12 2020-04-01 キヤノンメディカルシステムズ株式会社 X線ct装置及びデータ圧縮復元方法
JP2018509622A (ja) 2015-03-13 2018-04-05 カリフォルニア インスティチュート オブ テクノロジー フーリエタイコグラフィ手法を用いるインコヒーレント撮像システムにおける収差補正
US9993149B2 (en) 2015-03-25 2018-06-12 California Institute Of Technology Fourier ptychographic retinal imaging methods and systems
WO2016187591A1 (en) 2015-05-21 2016-11-24 California Institute Of Technology Laser-based fourier ptychographic imaging systems and methods
DE102015007939A1 (de) 2015-06-19 2016-12-22 Universität Stuttgart Verfahren und Computerprogrammprodukt zum Erzeugen eines hochaufgelösten 3-D-Voxeldatensatzes mit Hilfe eines Computertomographen
CN106530366B (zh) * 2015-09-09 2019-04-16 清华大学 能谱ct图像重建方法及能谱ct成像系统
CA3005439A1 (en) * 2015-11-20 2017-05-26 Integrated Dynamic Electron Solutions, Inc. Temporal compressive sensing systems
US10789737B2 (en) 2015-12-22 2020-09-29 Carestream Health, Inc. Tomographic image acquisition using asymmetrical pixel binning
EP3661334B1 (en) * 2016-05-31 2021-07-07 Koninklijke Philips N.V. Apparatus for generating x-rays
US10568507B2 (en) 2016-06-10 2020-02-25 California Institute Of Technology Pupil ptychography methods and systems
US11092795B2 (en) 2016-06-10 2021-08-17 California Institute Of Technology Systems and methods for coded-aperture-based correction of aberration obtained from Fourier ptychography
WO2018046454A1 (en) 2016-09-08 2018-03-15 Koninklijke Philips N.V. Radiation detector and x-ray imaging system
CN106683146B (zh) * 2017-01-11 2021-01-15 上海联影医疗科技股份有限公司 一种图像重建方法和图像重建算法的参数确定方法
CN106989835B (zh) * 2017-04-12 2023-07-11 东北大学 基于压缩感知的光子计数x射线能谱探测装置及成像系统
WO2019090149A1 (en) 2017-11-03 2019-05-09 California Institute Of Technology Parallel digital imaging acquisition and restoration methods and systems
CN109978809B (zh) * 2017-12-26 2022-02-22 同方威视技术股份有限公司 图像处理方法、装置及计算机可读存储介质
KR102128765B1 (ko) * 2018-05-16 2020-07-01 가천대학교 산학협력단 의료영상기기에서 압축센싱 활용을 위한 샘플링 패턴 산출장치 및 그 방법
US11408983B2 (en) * 2018-10-01 2022-08-09 Infineon Technologies Ag Lidar 2D receiver array architecture
US11039801B2 (en) 2019-07-02 2021-06-22 GE Precision Healthcare LLC Systems and methods for high-resolution spectral computed tomography imaging
CN112401912A (zh) * 2020-12-10 2021-02-26 杭州美诺瓦医疗科技股份有限公司 一种具有更低放射剂量的儿童骨龄成像方法及其成像装置
KR102540320B1 (ko) * 2021-07-20 2023-06-07 연세대학교 산학협력단 적응적 샘플링 기반 cbct 영상 획득 방법 및 장치

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6437338B1 (en) 1999-09-29 2002-08-20 General Electric Company Method and apparatus for scanning a detector array in an x-ray imaging system
WO2001039558A1 (en) 1999-11-23 2001-05-31 Koninklijke Philips Electronics N.V. X-ray examination apparatus with exposure control
US6389096B1 (en) * 2000-11-22 2002-05-14 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for providing additional computed tomography imaging modes
US7054406B2 (en) 2002-09-05 2006-05-30 Kabushiki Kaisha Toshiba X-ray CT apparatus and method of measuring CT values
JP4880587B2 (ja) * 2004-04-13 2012-02-22 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コンピュータ断層撮影のための動的線量制御
CN100536777C (zh) * 2004-04-13 2009-09-09 皇家飞利浦电子股份有限公司 用于计算机断层摄影的动态剂量控制
US20090310737A1 (en) * 2004-12-22 2009-12-17 Koninklijke Philips Electronics N.V. Method for computer tomography, and computer tomograph
EP1902424B1 (en) 2005-07-08 2009-11-11 Wisconsin Alumni Research Foundation Highly constrained image reconstruction method
EP1952342B1 (en) * 2005-11-17 2014-01-08 Philips Intellectual Property & Standards GmbH Method for displaying high resolution image data together with time-varying low resolution image data
CN100565336C (zh) * 2005-11-21 2009-12-02 清华大学 成像系统
US7317783B2 (en) 2006-04-21 2008-01-08 Pavel Dolgonos Reduced X-Ray exposure using power modulation
EP2027593A1 (en) 2006-05-22 2009-02-25 Philips Intellectual Property & Standards GmbH X-ray tube whose electron beam is manipulated synchronously with the rotational anode movement
US7881510B2 (en) * 2007-06-08 2011-02-01 Allegheny-Singer Research Institute Method and apparatus for forming an image with dynamic projective data
JP5220125B2 (ja) 2007-12-20 2013-06-26 ウイスコンシン アラムナイ リサーチ ファウンデーシヨン 先験的画像制限画像再構成法
US8194937B2 (en) 2007-12-20 2012-06-05 Wisconsin Alumni Research Foundation Method for dynamic prior image constrained image reconstruction
EP2240906B1 (en) 2008-01-14 2013-12-11 Wisconsin Alumni Research Foundation Method for prior image constrained progressive image reconstruction
US8135186B2 (en) * 2008-01-25 2012-03-13 Purdue Research Foundation Method and system for image reconstruction
US7697658B2 (en) 2008-02-01 2010-04-13 Virginia Tech Intellectual Properties, Inc. Interior tomography and instant tomography by reconstruction from truncated limited-angle projection data
US8472688B2 (en) 2008-04-17 2013-06-25 Wisconsin Alumni Research Foundation Method for image reconstruction employing sparsity-constrained iterative correction
WO2009155418A2 (en) * 2008-06-18 2009-12-23 Wright State University Computed tomography scanners, x-ray filters and methods thereof
US8952333B2 (en) * 2009-11-02 2015-02-10 Virginia Tech Intellectual Properties, Inc. Methods for improved single photon emission computed tomography using exact and stable region of interest reconstructions
US8705828B2 (en) * 2011-08-31 2014-04-22 Carestream Health, Inc. Methods and apparatus for super resolution scanning for CBCT system and cone-beam image reconstruction

Also Published As

Publication number Publication date
CN102947864A (zh) 2013-02-27
WO2011161558A1 (en) 2011-12-29
US9262845B2 (en) 2016-02-16
EP2583250B1 (en) 2014-12-17
US20130083886A1 (en) 2013-04-04
CN102947864B (zh) 2015-08-12
JP2013529491A (ja) 2013-07-22
RU2571564C2 (ru) 2015-12-20
EP2583250A1 (en) 2013-04-24
RU2013102504A (ru) 2014-07-27

Similar Documents

Publication Publication Date Title
JP5848759B2 (ja) 低線量ctイメージングを実行するための方法及びシステム
EP2049918B1 (en) Stereo tube computed tomography
JP5180181B2 (ja) コンピュータ断層撮影データ収集装置及び方法
US8731266B2 (en) Method and system for correcting artifacts in image reconstruction
JP2009545395A5 (ja)
US7396162B1 (en) Scatter correction for CT method and apparatus
US7995702B2 (en) System and method of data interpolation in fast kVp switching dual energy CT
JP2004188187A (ja) アーティファクト低減を容易にする方法及び装置
CN110179486B (zh) 多能量ct成像系统及其应用
EP2731504A2 (en) Imaging system detector calibration
US7362843B2 (en) System and method for reconstruction of cone beam tomographic projections with missing data
US8755488B2 (en) Voltage modulated x-ray tube
JP2013192951A (ja) 画像再構成装置、画像再構成方法およびx線コンピュータ断層撮影装置
JP7091262B2 (ja) コンピュータ断層撮影における空間分解能を向上させるためのシステムおよび方法
JP2010142478A (ja) X線ct装置
US9262825B2 (en) Image reconstruction in interleaved multi-energy imaging
Baek et al. Initial results with a multisource inverse-geometry CT system
CN106462983A (zh) 计算机化层析图像曝光和重构方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140424

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150226

A524 Written submission of copy of amendment under article 19 pct

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20150522

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: 20151105

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20151127

R150 Certificate of patent or registration of utility model

Ref document number: 5848759

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees