JP6169558B2 - コントラスト依存の解像度をもつ画像 - Google Patents

コントラスト依存の解像度をもつ画像 Download PDF

Info

Publication number
JP6169558B2
JP6169558B2 JP2014501744A JP2014501744A JP6169558B2 JP 6169558 B2 JP6169558 B2 JP 6169558B2 JP 2014501744 A JP2014501744 A JP 2014501744A JP 2014501744 A JP2014501744 A JP 2014501744A JP 6169558 B2 JP6169558 B2 JP 6169558B2
Authority
JP
Japan
Prior art keywords
image
contrast
resolution
gradient
reconstructed
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.)
Active
Application number
JP2014501744A
Other languages
English (en)
Other versions
JP2014514048A (ja
JP2014514048A5 (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 JP2014514048A publication Critical patent/JP2014514048A/ja
Publication of JP2014514048A5 publication Critical patent/JP2014514048A5/ja
Application granted granted Critical
Publication of JP6169558B2 publication Critical patent/JP6169558B2/ja
Active 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Description

本発明は、概して、再構成画像から、コントラストに依存する解像度をもつ画像を生成することに関し、コンピュータトモグラフィ(CT)の特定の用途に関して記述される。しかしながら、本発明は、例えばハイブリッドPET/CT及び/又はMRI/CTシステムのような他のイメージングモダリティ、デジタルX線システム、及び/又は他のイメージングモダリティにも適応する。
コンピュータトモグラフィ(CT)スキャナは、縦軸又はz軸を中心に検査領域の周りを回転する回転可能なガントリに取り付けられるX線管を有する。X線管は、検査領域及び検査領域内の被検体又は対象を横切る放射線を放出する。検出器アレイは、検査領域をはさんでX線管と反対側で円弧をなす。検出器アレイは、検査領域(及び検査領域内の被検体又は対象)を横切り、それらを示す投影データを生成する放射線を検出する。再構成器は、投影データを再構成し、投影データを示すボリュメトリック画像データを生成する。
残念ながら、CTスキャナは、イオン化放射線を放出し、ゆえに、患者をイオン化放射線に曝露させ、これは、癌のリスクを増大させうる。概して、患者に与えられる放射線量は、複数のファクタに依存し、かかるファクタは、非制限的な例として、管電流(mA)、管電圧(kVp)、曝露時間(ヘリカルスキャンの場合)、並びにスライスの厚さ及び間隔(軸方向スキャンの場合)を含む。従って、与えられる線量は、上述したものの1又は複数を調整することによって低減されることができる。しかしながら、画像ノイズは、放射線量に反比例し、従って、放射線量の低減は、患者に与えられる線量を低減するだけでなく、取得されるデータにおける画像ノイズを増大させる。画像ノイズは、再構成の間に画像に伝わるので、画像品質を低下させ(すなわち、よりノイズの多い画像)、プロシージャの診断値を低下させることがある。
増大されるノイズを低減するための1つのアプローチは、例えば画像内のノイズに関する強いペナルティを用いる最尤法(ML)アプローチ(例えば「ペナルティ付き尤度」又は「「正則化」アプローチ)のような、統計的な反復再構成技法を使用することである。しかしながら、例えばコントラストレベルがノイズレベルに近い低コントラストの対象の場合、正則化は、単一ノイズ実現において、現れうるエッジをシャープにするが、実際、ノイズによる大きい量の不確実さを含む。これは、図1に見ることができ、プロファイル102は、スキャン対象の低コントラスト領域のシャープなFBP再構成を表し、プロファイル104及び106は、スキャン対象の低コントラスト領域の正則化ML反復再構成を表す。分かるように、MLプロファイル104及び106は、FBPプロファイル102より低い解像度を有し、MLプロファイル104及び106の中間値は、FBPプロファイル102ほどシャープでない。
複数ノイズ実現の場合、重く正則化された画像は、コントラスト依存の解像度を示し、平均画像では、正則化された画像は、より低い解像度を有するが、従来のフィルタード逆投影(FBP)画像の解像度により近く、高コントラスト領域で一層シャープな解像度を維持しつつ、低コントラスト領域で非常に滑らかなフィルタをもつ。残念なことに、正則化尤度反復再構成の前方投影及び逆投影の演算は、グラフィック処理ユニット(GPU)のような並列ハードウェアの加速があるとしても、計算コストが高い。結果として、反復正則化尤度に基づく再構成は、大きな計算コストを伴い、日常的な医療診療におけるそれらの使用を妨げ又は禁じうる。
本発明の見地は、上述の問題及びその他に対処する。
1つの見地によれば、1又は複数の再構成画像を処理する方法であって、第1の画像解像度を有する第1の再構成画像を入力するステップと、第2の画像解像度を有する第2の再構成画像を入力するステップと、含む。第1の解像度は、第2の解像度より大きい。方法は更に、第1及び第2の再構成画像を組み合わせて単一画像にすることによって、コントラストに依存する解像度(以下、コントラスト依存解像度)の画像を生成するステップを含み、ここで、第1の再構成画像が、コントラスト依存解像度の画像のより高コントラストの領域に関して、コントラスト依存解像度の画像により大きく寄与し、第2の再構成画像が、コントラスト依存解像度の画像のより低コントラスト領域に関して、コントラスト依存解像度の画像により大きく寄与する。
別の見地によれば、システムは、コントラスト依存解像度の重み関数に基づいて、第1の解像度を有する第1の再構成画像及び第2の解像度を有する第2の再構成画像を選択的にブレンドすることによって、コントラスト依存解像度の画像を生成するコントラスト依存解像度画像生成器を有する。
別の見地によれば、コンピューティングシステムの1又は複数のプロセッサによって実行されるとき、コンピューティングシステムに、単一の反復において、より高い及びより低いコントラスト領域を有するより高い及びより低解像度の再構成画像を選択的に組み合わせることによって、正則化反復再構成の複数の反復を使用して生成される画像と同じコントラスト依存解像度特性をもつコントラスト依存解像度の画像を生成させるためのコンピュータ可読命令コンピュータ可読命令によって符号化されるコンピューティング可読記憶媒体が提供される。
本発明は、さまざまなコンポーネント及びコンポーネントの取り合わせ並びにさまざまなステップ及びステップの取り合わせの形を取ることができる。図面は、好適な実施形態を示すことのみを目的とし、本発明を制限するものとして解釈されるべきでない。
スキャン対象の低いコントラスト領域に関する、フィルタード逆投影(FBP)及びペナルティ付き/正則化最尤(ML)反復再構成のプロファイルを示す図。 画像データプロセッサ及びコントラスト依存解像度の画像生成器に関連して例示のイメージングシステムを概略的に示す図。 画像データプロセッサ、コントラスト依存解像度の画像生成器及び画像コンバイナの例を概略的に示す図。 コントラスト依存解像度の画像生成器の重み付けマトリクス生成器の例を示す図。 イメージングされたデータを処理し、コントラスト依存解像度の画像を生成する方法の例を示す図。 スキャン対象の低いコントラスト領域に関して、フィルタード逆投影(FBP)及び複数反復最尤(ML)再構成のプロファイルと比べて、本願明細書に記述されるコントラスト依存解像度の再構成のプロファイルをグラフで示す図。
図2は、コンピュータトモグラフィ(CT)スキャナのようなイメージングシステム200を示す。イメージングシステム200は、概して静止しているガントリ202及び回転ガントリ204を有する。回転ガントリ204は、静止ガントリ202によって回転可能に支持され、縦軸又はz軸208を中心に検査領域206の周りを回転する。
X線管のような放射線源210は、回転ガントリ204によって回転可能に支持される。放射線源210は、回転ガントリ204と共に回転し、検査領域206を横切る放射線を放出する。X線源コリメータは、概して円錐、扇形、くさび形又他の形状の放射線ビームを形成するように、放射線をコリメートするコリメーション部材を有する。
1次元又は2次元の放射線感受性検出器アレイ212は、検査領域206をはさんで放射線源210と反対側で円弧をなす。検出器アレイ212は、z軸208方向に沿って延在する検出器の複数の行を有する。検出器アレイ212は、検査領域206を横切る放射線を検出し、検出された放射線を示す投影データを生成する。
寝台のような患者支持体214が、例えば検査領域206内の人間患者のような対象又は被検体を支持する。支持体214は、対象又は被検体をロードし、スキャンし、及び/又はアンロードするために、対象又は被検体を移動させるように構成される。
汎用コンピュータシステム又はコンピュータが、オペレータコンソール216として役立つ。コンソール216のプロセッサは、コンピュータ可読記憶媒体(例えば、物理メモリ)に記憶されるコンピュータ可読命令を実行し、コンピュータ可読命令は、オペレータが、十分な線量又は低線量のスキャンプロトコルを選択し、コントラスト依存解像度の画像を生成し、スキャンを開始する等のシステム200の操作を制御することを可能にする。
再構成器218は、投影データを再構成し、それを示すボリュメトリック画像データを生成する。再構成器218は、従来のフィルタード逆投影再構成、コーンビームアルゴリズム、反復アルゴリズム、及び/又は他のアルゴリズムを用いることができる。更に、再構成器218は、より高い、より低い及び/又は他の解像度を有する画像を生成するために利用されることができる。
画像データプロセッサ220は、画像データを処理する。図示される画像データプロセッサ220は、少なくともコントラスト依存解像度画像生成器222を含む。一例において、コントラスト依存解像度画像生成器222は、コントラスト依存解像度の重み関数を使用して、解像度の異なる複数の再構成画像を組み合わせて又は他のやり方でブレンドして、単一画像にする。再構成画像は、より高解像度の画像及びより高い解像度画像から生成されるより低解像度の画像、又は同じ投影データから再構成されるより高い及びより低解像度の画像再構成を含む。更に、画像は、より高い及びより低いコントラストの領域及び/又は対象を含むことができる。
後で詳しく述べるように、上述した重み関数は、画像のより高コントラストの領域に関して、より高解像度の画像により高い重みを(及びゆえにより低解像度の画像により低い重み)を適用し、画像のより低いコントラスト領域に関しては、より低解像度の画像により高い重みを(及びゆえに、より高解像度の画像により低い重み)適用する。一例において、結果として得られる画像は、計算コストがずっと少ないながらも、正則化反復再構成と同じコントラスト依存解像度の特性を有する。更に、単一の再構成画像のエッジは、従来の正則化反復画像によって与えられる信頼できるエッジの幾分間違った印象に代わって、エッジが知られている真の正確さのより正確な反映を示す。更に、バックグラウンド領域のノイズパターンは、一般的なFBPノイズパターンに比較的慣れている顧客にとって不自然さがより低くみえる。
図示される実施形態において、画像データプロセッサ220は、システム200の一部として示される。この例において、画像データプロセッサ220は、別個のコンポーネント(図示される)、再構成器218の一部、コンソール216の一部、及び/又は別のコンポーネントの一部でありうる。代替の実施形態において、画像データプロセッサ220は、例えば別個のコンピューティングシステムの一部として、システム200から遠隔していてもよく、及び/又はコンピューティングシステムにわたって分散されてもよい。更に、画像データプロセッサ220は、コンピュータ可読記憶媒体(例えば、物理メモリ)に埋め込まれ又は符号化され及び/又は搬送波又は信号によって運ばれる、1又は複数のコンピュータ可読命令を実行する1又は複数のプロセッサを通じて実現されることができる。
図3は、画像データプロセッサ(IDP)220及びコントラスト依存解像度画像生成器(CDRIG)222の非限定的な例を概略的に示す。他の実現例も本願明細書において企図される。
図示される画像データプロセッサ220は、より高解像度の画像、又は両方が同じ投影データから再構成されるより高解像度の画像及びより低解像度の画像の両方を入力として受け取る。後者の例が、図3に示される。この例において、より高解像度の画像及びより低解像度の画像は、コントラスト依存解像度画像生成器222に提供され、画像平滑化器302は省かれることができる。
しかしながら、より高解像度の画像のみが、画像データプロセッサ220に提供される場合、画像平滑化器302は、より高解像度の画像を平滑化して、より高解像度の画像からより低解像度の画像を生成し、より高解像度の画像及び生成されたより低解像度の画像の両方が、コントラスト依存解像度画像生成器222に提供される。
図示されるコントラスト依存解像度画像生成器222は、コントラスト及び勾配レベル検出器304を有する。コントラスト及び勾配レベル検出器304は、より低解像度の画像のコントラスト及び勾配レベルを決定する。コントラスト依存解像度画像生成器222は、より低解像度の画像のコントラスト及び勾配レベルを表す勾配画像を生成し、出力する。
重み関数生成器306は、如何にして、勾配画像に基づいて、より低解像度の画像及びより高解像度の画像を組み合わせて、コントラスト依存解像度をもつ単一画像を生成するかを規定する重み関数を生成する。生成された重み関数は、より高コントラストの領域では、高解像度の画像がより大きく寄与し、低コントラストの領域では、より低解像度の画像がより大きく寄与するように、画像のそれぞれ異なる領域を重み付けるという点で位置依存である。重み関数生成器306は、重み関数を表す信号を生成し出力する。
画像コンバイナ308は、重み関数生成器306によって決定された重み関数に基づいて、より高解像度の画像及びより低解像度の画像を組み合わせ、コントラスト依存解像度の画像を生成し、出力する。ここに記述されるように、画像コンバイナ308は、より高解像度の画像及びより低解像度の画像をブレンドするために重み関数を用い、より高解像度の画像が、より高コントラストの領域に関して、コントラスト依存解像度の画像により大きく寄与し、より低解像度の画像が、より低コントラストの領域に関して、コントラスト依存解像度の画像により大きく寄与するようにする。
図示される実施形態において、画像データプロセッサは、任意には、ノイズ低減器310を有し、ノイズ低減器310は、画像コンバイナ308を通じてより低解像度の画像と組み合わせられるより高解像度の画像のノイズを低減するために使用される。ノイズ低減器310は、より高解像度の画像における多くのノイズを除去するように構成され、それにより、かかるノイズが、組み合わされた画像において不連続性をもたらさなくなる。ノイズ低減器310は、エッジを維持し、ノイズを除去するように構成される任意のノイズ低減アルゴリズムを用いることができ、例えば双方向フィルタ、全変動最小化に関するROFアルゴリズム等を用いることができる。
図3において、ノイズ低減器310及び画像平滑化器302は、画像データプロセッサ220の一部として示されている。別の実施形態では、ノイズ低減器310又は画像平滑化器302の少なくとも一方が、コントラスト依存解像度画像生成器222の一部である。別の実施形態では、ノイズ低減器310又は画像平滑化器302の少なくとも一方が、例えば再構成器218、コンソール216及び/又は別のコンピューティングシステムの一部として、画像データプロセッサ220の外部に存在する。
図4は、コントラスト及び勾配レベル検出器304、重み関数生成器(WFG)306及び画像コンバイナ308の非限定的な例を概略的に示す。他の実現例も本願明細書において企図される。
ここで述べられるように、コントラスト及び勾配レベル検出器304は、より低解像度の画像のコントラスト及び勾配レベルを決定し、勾配画像を出力する。この例において、コントラスト及び勾配レベル検出器304は、式1に示されるように表現される勾配画像を生成する:
Figure 0006169558
ここで、gは、より低解像度の画像u(x,y)の勾配画像を表す。
勾配カットオフ設定器402は、勾配画像に、最小勾配大きさレベル(gmin)及び最大勾配大きさレベル(gmax)を設定し、各レベルは、より低解像度の画像のみ又はより高解像度の画像のみが、コントラスト依存解像度の画像を生成する際に使用される勾配カットオフをそれぞれ決定する。これらの値は、省略時設定されることができ、ユーザにより決定されることができ、及び/又は調整可能であり、及び概して所望の視覚的プレゼンテーションを提供するように設定されることができる。
勾配マップ生成器404が、式2に示されるように、0と1の値の間の勾配をスケールする勾配マップを生成する。
Figure 0006169558
勾配マップ変換コンポーネント406は、重みマップ上の滑らかな遷移を確実にするように、関数(例えばシヌソイド関数(図示される)又は他の関数)を勾配マップに適用し、式3に示される重み関数を与えるために、使用されることができる。
Figure 0006169558
画像コンバイナ308は、重み関数、より高解像度の画像、及びより低解像度の画像を入力し、式4に基づいて、コントラスト依存解像度の画像を生成する。
Figure 0006169558
ここで、コントラスト依存の解像度は、より高解像度の画像及びより低解像度の画像のブレンドされた又は重み付けされた合計である。
図5は、画像データを処理し、処理された画像データから、コントラスト依存解像度の画像を生成する例示の方法を示す。
ここに記述される方法のステップの順序は制限的なものではない。従って、他の順序も本願明細書において企図される。更に、1又は複数のステップが省かれることができ、及び/又は1又は複数の付加のステップが含められることができる。
ステップ502において、対象又は被検体が、イメージングシステム200又は他のイメージングシステムによってスキャンされる。
ステップ504において、スキャンから投影データが取得される。
ステップ506において、投影データが再構成され、少なくともより高解像度の画像及び任意には対応するより低解像度の画像が生成される。
ステップ508において、対応するより低解像度の画像が投影データから生成されない場合に、対応するより低解像度の画像が、より高解像度の画像から生成されることができる。
ステップ510において、より低解像度の画像のコントラスト及び勾配レベルが検出される。ここに記述されるように、この情報は、勾配画像において表現されることができる。
ステップ512において、より高解像度の画像及びより低解像度の画像を組み合わせるための重み関数が、ここに記述されるようにコントラスト及び勾配レベルに基づいて又は他のやり方で、生成される。
ステップ514において、任意には、より高解像度の画像が、エッジを維持しつつノイズを低減するアルゴリズムによって処理される。
ステップ516において、フィルタリングされたより高解像度の画像及びより低解像度の画像が、ここに記述されるように重み関数に基づいて組み合わせられ、コントラスト依存解像度の画像を生成する。
上述したものは、コンピュータプロセッサによって実行されるとき、上述したステップをプロセッサに実施させるためのコンピュータ可読命令によって実現されることができる。そのような場合、命令は、当該コンピュータに関連付けられ又は当該コンピュータにアクセス可能であるコンピュータ可読記憶媒体に記憶される。付加的に又は代替として、命令の1又は複数が、搬送波又は信号によって運ばれることができる。
図6は、シャープなFBP再構成のプロファイル602、複数反復ML再構成のプロファイル604、及び本願明細書に記述されるコントラスト依存解像度再構成のプロファイル606の間の解像度の比較を示す。図示されるように、コントラスト依存再構成の解像度及び複数反復ML再構成の解像度はほぼ同じである。しかしながら、コントラスト依存再構成は、複数反復ML再構成より計算コストが低い、シングルパスの再構成である。
本発明は、さまざまな実施形態に関してここに記述された。変形及び変更が、本願明細書の説明を読むことにより、当業者に見出されることができる。このような変形及び変更が、添付の請求項又はそれと等価なものの範囲にある限り、本発明は、すべてのこのような変更及び変更を含むものとして解釈されることが意図される。

Claims (13)

  1. 1又は複数の再構成画像を処理する方法であって、
    第1の画像解像度を有する第1の再構成画像を入力するステップと、
    第2の画像解像度を有する第2の再構成画像を入力するステップであって、前記第1の画像解像度が、前記第2の画像解像度より大きい、ステップと、
    前記第1の再構成画像及び前記第2の再構成画像を組み合わせて単一画像にすることにより、コントラスト依存解像度の画像を生成するステップであって、前記コントラスト依存解像度の画像のより高コントラストの領域をもたらすように前記第1の再構成画像が前記コントラスト依存解像度の画像により大きく寄与し、前記コントラスト依存解像度の画像のより低コントラストの領域をもたらすように前記第2の再構成画像が前記コントラスト依存解像度の画像により大きく寄与するようにされる、ステップと、
    を含む方法。
  2. 1又は複数の再構成画像を処理する方法であって、
    第1の画像解像度を有する第1の再構成画像を入力するステップと、
    第2の画像解像度を有する第2の再構成画像を入力するステップであって、前記第1の画像解像度が、前記第2の画像解像度より大きい、ステップと、
    前記第1の再構成画像及び前記第2の再構成画像を組み合わせて単一画像にすることにより、コントラスト依存解像度の画像を生成するステップであって、前記コントラスト依存解像度の画像のより高コントラストの領域について、前記第1の再構成画像がより大きく寄与し、前記コントラスト依存解像度の画像のより低コントラストの領域について、前記第2の再構成画像がより大きく寄与するようにされる、ステップと、
    を含み、
    前記コントラスト依存解像度の画像を生成する前記ステップは、コントラスト依存解像度の重み関数に基づいて、前記第1及び前記第2の再構成画像を組み合わせることを含み、
    前記重み関数は、前記第2の再構成画像のコントラスト及び勾配レベルに基づいて、如何にして、より低解像度の再構成画像及びより高解像度の再構成画像を組み合わせてコントラスト依存解像度の画像を生成するか、を規定し、
    前記第2の再構成画像のコントラスト及び勾配レベルを検出するステップと、
    前記検出されたコントラスト及び勾配レベルに基づいて勾配画像を生成するステップと、
    前記勾配画像に基づいて前記重み関数を生成するステップと、
    を更に含む、方法。
  3. 前記コントラスト依存解像度の画像を生成する前記ステップは、コントラスト依存解像度の重み関数に基づいて、前記第1及び前記第2の再構成画像を組み合わせることを含み、
    前記重み関数は、前記第2の再構成画像のコントラスト及び勾配レベルに基づいて、如何にして、より低解像度の再構成画像及びより高解像度の再構成画像を組み合わせてコントラスト依存解像度の画像を生成するか、を規定し、
    前記第2の再構成画像のコントラスト及び勾配レベルを検出するステップと、
    前記検出されたコントラスト及び勾配レベルに基づいて勾配画像を生成するステップと、
    前記勾配画像に基づいて前記重み関数を生成するステップと、
    を更に含む、
    請求項1に記載の方法。
  4. 前記勾配画像に、予め決められた最小勾配大きさレベル及び予め決められた最大勾配大きさレベルを設定するステップであって、前記予め決められた最小及び最大勾配大きさレベルは、それぞれ、前記第2の再構成画像のみが前記コントラスト依存解像度の画像のために使用される勾配レベル、及び前記第1の再構成画像のみが前記コントラスト依存解像度画像のために使用される勾配レベルを示す、ステップと、
    前記最小及び最大勾配大きさレベルを有する前記勾配画像に基づいて、前記重み関数を生成するステップと、
    を更に含む、請求項2又は3に記載の方法。
  5. 前記最小及び前記最大勾配大きさレベルを有する前記勾配画像に基づいて、0から1までの勾配をスケールする勾配マップを生成するステップと、
    前記勾配マップに基づいて、前記重み関数を生成するステップと、
    を更に含む、請求項4に記載の方法。
  6. 滑らかな重み遷移を確実にする関数により前記勾配マップを変換するステップと、
    前記変換された勾配マップに基づいて、前記重み関数を生成するステップと、
    を更に含む、請求項5に記載の方法。
  7. 前記第1の再構成画像のノイズを低減するステップと、
    前記ノイズが低減された第1の再構成画像及び前記第2の再構成画像を組み合わせることによって、前記コントラスト依存解像度の画像を生成するステップと、
    を更に含む、請求項1乃至6のいずれか1項に記載の方法。
  8. 前記第1の再構成画像のノイズを低減する前記ステップは、前記第1の再構成画像のエッジを維持しつつ、前記第1の再構成画像のノイズを低減するノイズ低減アルゴリズムを適用することを含む、請求項7に記載の方法。
  9. 前記第1の再構成画像は、イメージングシステムによる対象又は被検体のスキャンから投影データを再構成することによって生成され、前記第2の再構成画像は、前記投影データを再構成することによって又は前記第1の再構成画像から、生成される、請求項1乃至8のいずれか1項に記載の方法。
  10. コントラスト依存解像度の重み関数に基づいて、第1のより高い解像度を有する第1の再構成画像及び第2のより低い解像度を有する第2の再構成画像を選択的に組み合わせることによって、コントラスト依存解像度の画像を生成するコントラスト依存解像度画像生成器を有するシステムであって、
    前記コントラスト依存解像度画像生成器は、
    より低解像度の画像のコントラスト及び勾配レベルを検出し、それを示す勾配画像を生成するコントラスト及び勾配レベル検出器と、
    前記勾配画像に基づいて、コントラスト依存解像度の重み関数を生成する重み関数生成器と、
    前記コントラスト依存解像度の重み関数に基づいて、第1の再構成画像及び第2の再構成画像を組み合わせる画像コンバイナと、
    を有する、
    システム。
  11. コントラスト依存解像度の重み関数に基づいて、第1のより高い解像度を有する第1の再構成画像及び第2のより低い解像度を有する第2の再構成画像を選択的に組み合わせることによって、コントラスト依存解像度の画像を生成するコントラスト依存解像度画像生成器を有するシステムであって、
    前記重み関数は、コントラスト依存解像度の画像のより高コントラストの領域について、より高解像度の画像をより高く重み付けし、コントラスト依存解像度の画像のより低コントラストの領域について、より低解像度の画像をより高く重み付けし、
    前記コントラスト依存解像度画像生成器は、
    より低解像度の画像のコントラスト及び勾配レベルを検出し、それを示す勾配画像を生成するコントラスト及び勾配レベル検出器と、
    前記勾配画像に基づいて、コントラスト依存解像度の重み関数を生成する重み関数生成器と、
    前記コントラスト依存解像度の重み関数に基づいて、第1の再構成画像及び第2の再構成画像を組み合わせる画像コンバイナと、
    を有する、システム。
  12. 前記重み関数は、コントラスト依存解像度の画像のより高コントラストの領域をもたらすようにより高解像度の画像をより高く重み付けし、コントラスト依存解像度の画像のより低コントラストの領域をもたらすようにより低解像度の画像をより高く重み付けする、請求項10に記載のシステム。
  13. 前記重み関数生成器が、
    最小及び最大勾配レベルを前記勾配画像に設定する勾配カットオフ設定器と、
    前記勾配画像及び前記最小及び最大勾配レベルに基づいて勾配マップを生成する勾配マップ生成器と、
    前記勾配マップを重み関数に変換する勾配マップ変換コンポーネントと、
    を有する、請求項10乃至12の何れか一項に記載のシステム
JP2014501744A 2011-03-28 2012-03-16 コントラスト依存の解像度をもつ画像 Active JP6169558B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161468126P 2011-03-28 2011-03-28
US61/468,126 2011-03-28
PCT/IB2012/051266 WO2012131520A2 (en) 2011-03-28 2012-03-16 Contrast-dependent resolution image

Publications (3)

Publication Number Publication Date
JP2014514048A JP2014514048A (ja) 2014-06-19
JP2014514048A5 JP2014514048A5 (ja) 2015-04-16
JP6169558B2 true JP6169558B2 (ja) 2017-07-26

Family

ID=45937482

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014501744A Active JP6169558B2 (ja) 2011-03-28 2012-03-16 コントラスト依存の解像度をもつ画像

Country Status (7)

Country Link
US (1) US9336612B2 (ja)
EP (1) EP2691932B1 (ja)
JP (1) JP6169558B2 (ja)
CN (1) CN103608839B (ja)
BR (1) BR112013024733A2 (ja)
RU (1) RU2585790C2 (ja)
WO (1) WO2012131520A2 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6312401B2 (ja) * 2012-11-30 2018-04-18 キヤノン株式会社 画像処理装置、画像処理方法、及びプログラム
EP2984632B1 (en) 2013-04-10 2020-10-07 Koninklijke Philips N.V. Image quality index and/or imaging parameter recommendation based thereon
WO2014167450A1 (en) 2013-04-10 2014-10-16 Koninklijke Philips N.V. Reconstructed image data visualization
WO2014192831A1 (ja) * 2013-05-28 2014-12-04 株式会社東芝 医用画像診断装置及び制御方法
WO2015014678A1 (en) * 2013-07-30 2015-02-05 Koninklijke Philips N.V. Combined mri pet imaging
CN103955610B (zh) * 2014-04-22 2017-04-26 青岛大学附属医院 一种医学影像计算机辅助分析方法
CN106233327B (zh) 2014-04-23 2019-11-19 皇家飞利浦有限公司 去噪图像数据中对低对比度结构的恢复
WO2016042466A2 (en) * 2014-09-15 2016-03-24 Koninklijke Philips N.V. Iterative image reconstruction with a sharpness driven regularization parameter
CN104660899A (zh) * 2014-09-29 2015-05-27 上海华勤通讯技术有限公司 照片的处理方法
US10078057B2 (en) * 2015-09-02 2018-09-18 Fei Company Data processing in a tomographic imaging apparatus
US10255696B2 (en) 2015-12-11 2019-04-09 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
EP3234919B1 (en) 2015-12-11 2021-03-17 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
US20170294033A1 (en) * 2016-04-06 2017-10-12 Varex Imaging Corporation Dose efficient x-ray detector and method
WO2019036446A2 (en) * 2017-08-17 2019-02-21 The United States Of America, As Represented By The Secretary, Department Of Health And Human Services HYBRID TOMODENSITOMETRY SYSTEM WITH ADDITIONAL DETECTORS IMMEDIATELY NEAR THE BODY
US10719934B2 (en) * 2018-10-17 2020-07-21 Canon Medical Systems Corporation Devices, systems, and methods for medical imaging
CN109587462A (zh) * 2018-11-30 2019-04-05 中山大学 高对比度投影画面调整方法、装置、投影装置及投影系统

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4463375A (en) 1982-09-07 1984-07-31 The Board Of Trustees Of The Leland Standford Junior University Multiple-measurement noise-reducing system
US5754697A (en) * 1994-12-02 1998-05-19 Fu; Chi-Yung Selective document image data compression technique
US6208765B1 (en) * 1998-06-19 2001-03-27 Sarnoff Corporation Method and apparatus for improving image resolution
US6724942B1 (en) * 1999-05-24 2004-04-20 Fuji Photo Film Co., Ltd. Image processing method and system
JP4053185B2 (ja) * 1999-06-22 2008-02-27 富士フイルム株式会社 画像処理方法および装置
KR100344807B1 (ko) * 2000-01-12 2002-07-20 엘지전자주식회사 영상 신호 보정 장치 및 방법
JP4556276B2 (ja) * 2000-03-23 2010-10-06 ソニー株式会社 画像処理回路及び画像処理方法
US6778692B1 (en) * 2000-08-11 2004-08-17 General Electric Company Image processing method and apparatus including image improving circuit
WO2002048964A1 (en) * 2000-12-14 2002-06-20 Matsushita Electric Works, Ltd. Image processor and pattern recognition apparatus using the image processor
US6816627B2 (en) * 2001-04-12 2004-11-09 Lockheed Martin Corporation System for morphological image fusion and change detection
DE10150428A1 (de) 2001-10-11 2003-04-30 Siemens Ag Verfahren zur Erzeugung dreidimensionaler, mehrfachaufgelöster Volumenbilder eines Untersuchungsobjekts
JP4130114B2 (ja) * 2002-10-09 2008-08-06 株式会社日立メディコ 超音波イメージング装置及び超音波信号処理方法
JP4413504B2 (ja) * 2003-02-13 2010-02-10 株式会社東芝 医用画像処理装置、医用画像処理方法および医用画像処理プログラム
US7492967B2 (en) 2003-09-24 2009-02-17 Kabushiki Kaisha Toshiba Super-resolution processor and medical diagnostic imaging apparatus
JP2007508092A (ja) * 2003-10-13 2007-04-05 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁気共鳴画像化システム及び方法
CN100502786C (zh) * 2004-05-19 2009-06-24 株式会社日立医药 图像处理装置及其方法
JP4602703B2 (ja) * 2004-06-14 2010-12-22 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 放射線撮影装置および画像生成装置
US7623691B2 (en) * 2004-08-06 2009-11-24 Kabushiki Kaisha Toshiba Method for helical windmill artifact reduction with noise restoration for helical multislice CT
US20070014468A1 (en) * 2005-07-12 2007-01-18 Gines David L System and method for confidence measures for mult-resolution auto-focused tomosynthesis
WO2007012997A2 (en) * 2005-07-26 2007-02-01 Koninklijke Philips Electronics, N.V. Cardiac region detection from motion analysis of small scale reconstruction
CN101305397B (zh) * 2005-10-12 2012-09-19 有源光学有限公司 基于多个图像帧形成图像的方法,图像处理系统和数字照相机
US7983462B2 (en) 2005-11-22 2011-07-19 Purdue Research Foundation Methods and systems for improving quality of an image
JP2009533163A (ja) * 2006-04-13 2009-09-17 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コイル感度の波動伝搬に基づく評価
JP2009540942A (ja) * 2006-06-22 2009-11-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ デュアルx線管ゲーティング
JP2009238009A (ja) * 2008-03-27 2009-10-15 Toshiba Corp 画像処理装置、及び方法
DE102008021639A1 (de) 2008-04-30 2009-12-03 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum Erzeugen einer CT-Rekonstruktion eines Objekts mit einem hochaufgelösten interessierenden Objektbereich
TWI457853B (zh) * 2009-03-24 2014-10-21 Ind Tech Res Inst 提供深度資訊之影像處理方法及其影像處理系統
CN101915901A (zh) * 2010-08-17 2010-12-15 中国科学院深圳先进技术研究院 磁共振成像方法及装置

Also Published As

Publication number Publication date
CN103608839B (zh) 2017-04-12
US20140016849A1 (en) 2014-01-16
BR112013024733A2 (pt) 2016-12-27
CN103608839A (zh) 2014-02-26
JP2014514048A (ja) 2014-06-19
WO2012131520A2 (en) 2012-10-04
US9336612B2 (en) 2016-05-10
RU2013147816A (ru) 2015-05-10
RU2585790C2 (ru) 2016-06-10
EP2691932B1 (en) 2019-07-24
WO2012131520A3 (en) 2012-11-22
EP2691932A2 (en) 2014-02-05

Similar Documents

Publication Publication Date Title
JP6169558B2 (ja) コントラスト依存の解像度をもつ画像
JP6100772B2 (ja) 画像処理方法及びコンピューティング装置
JP5080986B2 (ja) 画像診断システム及び画像診断方法
JP5937093B2 (ja) 低線量ctノイズを除去するためのシステム及び方法
CN102270350B (zh) 结合四维噪声滤波器的迭代ct图像重建
US9189832B2 (en) Method and system for noise reduction in low dose computed tomography
JP6275826B2 (ja) ノイズ除去再構成画像データエッジ改善
US10789738B2 (en) Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator
US10878544B2 (en) Image data processing
US8938108B2 (en) Method for artifact reduction in cone-beam CT images
KR20120116940A (ko) 향상된 화상 데이터/선량 감소
Chang et al. Modeling and pre-treatment of photon-starved CT data for iterative reconstruction
JP6158910B2 (ja) 正則化による反復画像再構成
CN102270349A (zh) 在没有调节项的情况下对ct 图像的迭代重建
CN110073412B (zh) 使用交替否定的图像噪声估计
CN107341836B (zh) 一种ct螺旋扫描图像重建方法及装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150225

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150225

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20151211

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151217

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20160317

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160617

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160913

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20161213

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20170214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170313

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170628

R150 Certificate of patent or registration of utility model

Ref document number: 6169558

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250