JP7123927B2 - 交互打ち消しを使用する画像ノイズ推定 - Google Patents

交互打ち消しを使用する画像ノイズ推定 Download PDF

Info

Publication number
JP7123927B2
JP7123927B2 JP2019530016A JP2019530016A JP7123927B2 JP 7123927 B2 JP7123927 B2 JP 7123927B2 JP 2019530016 A JP2019530016 A JP 2019530016A JP 2019530016 A JP2019530016 A JP 2019530016A JP 7123927 B2 JP7123927 B2 JP 7123927B2
Authority
JP
Japan
Prior art keywords
projection data
image
noise
data
views
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
JP2019530016A
Other languages
English (en)
Other versions
JP2019535451A5 (ja
JP2019535451A (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 JP2019535451A publication Critical patent/JP2019535451A/ja
Publication of JP2019535451A5 publication Critical patent/JP2019535451A5/ja
Application granted granted Critical
Publication of JP7123927B2 publication Critical patent/JP7123927B2/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
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/401Imaging image processing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/419Imaging computed tomograph
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical

Landscapes

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

Description

以下は、広くは、体積画像データをノイズ除去するために画像ノイズを推定することに関し、特に、交互打ち消し(alternating negation)アプローチを使用して画像ノイズを推定することに関し、コンピュータ断層撮影(CT)及び/又は他のX線撮像のようなX線撮像に対する特定の応用とともに説明される。
CTスキャナは、z軸について検査領域の周りを回転する回転可能ガントリ上に取り付けられたX線管を含む。検出器アレイは、X線管から検査領域の反対側の角度的弧に対する。X線管は、検査領域並びにその中の対象及び/又はオブジェクトの一部を横切る放射線を放射する。検出器アレイは、検査領域を横切る放射線を検出し、それを示す投影データを生成する。再構成器は、投影データを再構成し、検査領域並びにその中の対象及び/又はオブジェクトの一部を示す体積(三次元)画像データを生成する。以下は、ここで使用される投影データを説明する。
図1において、焦点スポット102は、X線放射線104が放射される放射線源の領域を表す。コリメータ106のようなビーム成形装置は、検査領域108を横切るビーム108(例えば図示された例においてファンビーム)を規定する。検出器アレイ112は、検査領域110を横切って、焦点スポット102の反対側に配置される。検出器アレイ112は、少なくとも1行の検出器素子1141、...、114N(114)を含む。各検出器素子114は、そこに衝突するファンビームからX線放射線を検出し、焦点スポット102から検出器素子114までの経路に沿った合計X線減衰(線積分)を示す電気信号を生成する。
第三世代CTスキャナ技術を用いて、X線管(したがって焦点スポット102)及び検出器アレイ112は、検査領域110の周りの円形経路116に沿って協調して回転する。図1は、経路116上で1つの角度113(例えば、θ1)において焦点スポット102、ファンビーム108及び検出器アレイ112を示し、図2は、経路116上の異なる角度115(例えば、θN)において示す。任意の1つの角度における検出器素子114に対する線積分のセットが、ここでビューと称される。角度113における線積分のセットは、1つのビューを表し、角度115における線積分のセットは、他のビューを表す。スキャンに対するビューのセットは、ここで投影データと称される。
ノイズ除去アルゴリズムは、投影領域(線積分)において及び/又は画像領域(体積画像データ)においてノイズ除去するアルゴリズムを含む。反復的画像領域ノイズ除去に対して、各画素における入力画像ノイズレベルのマップを提供するノイズ画像が、空間的に変化するノイズレベルに対するノイズ除去強度を調整するのに使用される。US7706497B2は、画像領域ノイズ除去に対するノイズ画像を決定する方法を記載している。この方法は、投影データを奇数ビュー(例えば、θ1、θ3、θ5、...)及び偶数ビュー(例えば、θ2、θ4、θ6、...)(又は管の焦点がビューごとに往復して切り替えられる場合に異なる焦点)に分割することを含む。奇数ビューのセットは、奇数ビュー画像を生成するように1つの処理チェーンにより再構成され、偶数セットのビューは、偶数ビュー画像を生成するように異なる処理チェーンにより再構成される。奇数及び偶数ビュー画像は、ノイズ画像を作成するように減算される。このアイデアは、2つのセットのビューが同じ構造のスキャンからのものであることである。このようにして、両方の画像は、本質的に同じ構造を表し、2つの画像を減算することは、構造を相殺し、無相関のノイズのみを残す。
図3は、このプロセスをフロー図により示す。スキャンから生成された投影データは、302において受信され、304において奇数ビュー及び偶数ビューに分割される。奇数ビュー処理チェーン306及び308は、ファンビーム幾何構成から平行ビーム幾何構成に奇数ビューをリビニングし、リビニングされた奇数ビューを再構成し、奇数ビュー画像を生成する。偶数ビュー処理チェーン310及び312は、ファンビーム幾何構成から平行ビーム幾何構成に偶数ビューをリビニングし、リビニングされた偶数ビューを再構成し、偶数ビュー画像を生成する。314において、奇数及び偶数ビュー画像は、減算され、ノイズ画像を生成する。処理チェーン316及び318は、ファンビーム幾何構成から平行ビーム幾何構成に全てのビューをリビニングし、リビニングされたビューを再構成し、構造画像を生成する。320において、局所ノイズ分散推定が、ノイズ画像から推定され、構造画像に対するノイズマップを生成し、構造画像は、ノイズマップを使用してノイズ除去される。
不幸なことに、ノイズ画像を生成するこのアプローチは、ビューが半分に分割されるので角度方向又は半径方向の解像度を低減し、したがって、リビニング及び再構成チェーンにおいて使用される解像度に負の影響を持つ。リビニングプロセスにおいて、奇数及び偶数ビューは、異なるグリッドに対する補間により再サンプリングされる。リビニング後の奇数及び偶数データセットは、したがって、最悪の解像度を持ち、高周波ノイズ成分は、正確に再構成されることができない。結果として、最終的なノイズ推定は、再構成に対して使用されるフィルタカーネルに依存して不正確であり得、これは、最適未満のノイズ除去及びノイズ除去画像を生じうる。結果として、ノイズ画像を生成する他のアプローチに対する解決されていない要望が存在する。
ここに記載される態様は、上記の問題等に対処する。
一態様において、撮像システムは、X線放射線を放射するように構成された放射線源と、X線放射線を検出し、これを示す投影データを生成する検出器アレイであって、前記投影データは、複数のビューを含む、検出器アレイと、前記投影データを再構成し、ノイズのみ画像(noise only image)を生成するように構成された第1の処理チェーンとを含む。
他の態様において、方法は、撮像システムにより生成された投影データを受信するステップと、前記投影データを再構成し、ノイズのみ画像を生成するように構成された第1の処理チェーンを用いて前記投影データを処理するステップとを含む。
他の態様において、コンピュータ可読記憶媒体は、コンピュータ可読命令をエンコードされる。前記コンピュータ可読命令は、プロセッサにより実行される場合に、前記プロセッサに、X線撮像システムを用いてオブジェクト又は対象をスキャンさせ、投影データを生成させ、前記投影データを再構成してノイズのみ画像を生成するように構成された第1の処理チェーンを用いて前記投影データを処理させ、前記投影データを再構成して構造画像を生成するように構成された第2の処理チェーンを用いて前記投影データを処理させ、前記ノイズのみ画像に基づいて前記構造画像をノイズ除去させる。
本発明は、様々なコンポーネント及びコンポーネントの構成、並びに様々なステップ及びステップの構成の形を取りうる。図面は、好適な実施例を説明する目的のみであり、本発明を限定するように解釈されるべきではない。
投影データの複数のビューの1つを描写する。 複数のビューの異なるビューを描写する。 ビューを奇数ビュー及び偶数ビューに分割することにより画像領域ノイズ除去に対するノイズ画像を生成する従来技術の方法のフロー図である。 投影データプロセッサを含み、ノイズのみ画像を生成する処理チェーンを持つ一例の撮像システムを概略的に示す。 投影データプロセッサの一例を概略的に示す。 一実施例による一例の方法を示す。 ここに記載されるアプローチで生成されたノイズ画像と従来技術のアプローチで生成されたノイズ画像との比較を示す。
図4は、コンピュータ断層撮影(CT)スキャナのような撮像システム400を示す。撮像システム400は、概して静止しているガントリ402及び回転ガントリ404を含む。回転ガントリ404は、ベアリング等により静止ガントリ402により回転可能に支持され、長手軸又はz軸について検査領域406の周りを回転する。
X線管のような放射線源408は、回転ガントリ404により回転可能に支持される。放射線源408は、回転ガントリ404とともに回転し、検査領域406を横切る放射線を放射する。線源コリメータは、概してファン形状の放射線ビームを形成するように前記放射線をコリメートするコリメーション部材を含む。他の実施例において、前記ビームは、概してコーン、楔又は他の形状の放射線ビームを形成するように成形される。
一又は二次元放射線感知検出器アレイ410は、検査領域406を横切って放射線408の反対側の角度的弧に対する。検出器アレイ410は、z軸方向に沿って延在する検出の1以上の行を含む。検出器アレイ410は、検査領域406を横切る放射線を検出し、それを示す投影データを生成する。
ここに記載されるように、前記投影データは、複数のビューを含み、各ビューは、検査領域406の周りの放射線源408及び検出器アレイ410として、データが取得される回転の異なる角度に対応する。各ビューは、検出器アレイ410の検出器により出力された複数の値(線積分)を含み、各値は、放射線408から検出器アレイ410の当該値を生成する検出器までの経路に沿ったX線減衰を示す。
第1の(従来の)処理チェーン412は、リビニングプロセッサ414及び再構成プロセッサ416を含む。第2の処理チェーン418は、投影データプロセッサ420、リビニングプロセッサ414及び再構成プロセッサ416を含む。この例において、処理チェーン412及び418は、リビニングプロセッサ及び再構成プロセッサ416を共有する。変形例において、処理チェーン412及び418の各々は、独自のリビニングプロセッサ及び/又は再構成プロセッサを持つ。
第1の処理チェーン412に対して、リビニングプロセッサ414は、ファンビーム幾何構成から平行ビーム幾何構成に前記ビューをリビニングする。再構成プロセッサ416は、前記リビニングされたビューを再構成し、(「構造画像」と称される)体積画像データを生成する。第2の処理チェーン418に対して、投影データプロセッサ420は、リビニングプロセッサ414によるリビニング及び再構成プロセッサ416による再構成の前に前記ビューを前処理し、前記投影データにおいて表された構造は、リビニング及び再構成動作中に相殺し、ノイズのみ画像を生成する。
投影データプロセッサ420の一例は、図5に示される。乗算器502は、前記投影データのビューに所定の乗算係数504を乗算する。一例において、これは、第1のビューにプラス1の値を、第2のビューにマイナス1の値を、第3のビューにプラス1の値を、第4のビューにマイナス1の値を、...(n-1)番目のビューにプラス1の値を、n番目のビューにマイナス1の値を乗算することを含み、逆も同様であり、その場合、第1のビューがマイナス1の値を乗算され、第2のビューがマイナス1の値を乗算され、...以下同様である。このアプローチは、ここで交互打ち消しアプローチと称される。
一般に、所定の乗算係数504は、ビューごとにプラス1とマイナス1との間で交代する。変形例において、乗算係数504は、ビューごとに交代する必要はない。例えば、他の例において、乗算係数504は、代わりに、2又はそれ以上のビューごとに交代することができる。例えば、これは、第1のビューにプラス1の値を、第2のビューにプラス1の値を、第3のビューにマイナス1の値を、第4のビューにマイナス1の値を、...乗算することを含む。他の変形例において、所定の乗算係数504は、プラス及びマイナス1以外の値を含むことができる。
図4に戻って、第1の処理チェーン412の出力は、ノイズを持つスキャンされた構造を表す画素/ボクセルを持つ従来の体積画像データであり、第2の処理チェーン418の出力は、ノイズのみ画像である。この例において、両方の処理チェーン412及び418を用いて、リビニングプロセッサ414は、従来の又は他のファンから平行への(fan-to-parallel)ビームリビニングアルゴリズムを実施し、再構成プロセッサ416は、従来の又は他のフィルタ逆投影(FBP)再構成アルゴリズム及び/又は反復的再構成のような他の再構成アルゴリズム及び/又は他の再構成アルゴリズムを実施する。
ファンから平行へのビームリビニングの非限定的な例は、1997年11月26日に出願され、"Real time continuous CT imaging"と題されたUS6078639Aに記載されており、フィルタ逆投影再構成の非限定的な例は、2004年2月9日に出願され、"System and method for helical cone-beam computed tomography with exact reconstruction"と題されたUS2006/0140335A1に記載されており、両方の文献が、参照により全体的にここに組み込まれる。他のアプローチも、ここで考えられる。
ノイズ除去プロセッサ422は、前記再構成画像(3D画像データ)及び前記ノイズ画像の両方を受信し、前記ノイズ画像を使用して前記再構成画像をノイズ除去する。例えば、ノイズ除去プロセッサ422は、前記ノイズ画像からノイズ分散及び/又はノイズ標準偏差の推定値又はノイズマップを生成する。一例において、これは、小さな関心領域(ROI)の局所的分散又は標準偏差を取り、これを画像中で移動させることにより達成される。ノイズ標準偏差を使用する一例の画像領域ノイズ除去アルゴリズムは、2012年11月12日に出願され、"Image domain de-noising"と題されたUS9159122B2、2014年6月26日に出願され、"Methods of utilizing image noise information"と題されたUS2016/0140725A1、及び2015年10月29日に出願され、"Enhanced image data/dose reduction"と題されたUS8938110B2に記載されており、これらは、参照により全体的にここに組み込まれる。
カウチのような対象支持台424は、検査領域406において人間の患者のようなオブジェクト又は対象を支持する。対象支持台424は、前記オブジェクト又は対象を乗せる、スキャンする及び/又は降ろすために前記オブジェクト又は対象を移動するように構成される。汎用計算システム又はコンピュータは、オペレータコンソール426として機能する。コンソール426は、モニタのような人間可読出力装置と、キーボード、マウス等のような入力装置とを含む。コンソール426上にあるソフトウェアは、オペレータがグラフィカルユーザインタフェース(GUI)により又は他の方法で撮像システム400とインタラクトする及び/又は操作することを可能にする。これは、前記再構成画像をノイズ除去するのにノイズ除去プロセッサ422を使用するもののような撮像プロトコルを選択することを含む。
上記のことは、撮像システム400から直接的に来る投影データを示すが、変形例において、前記投影データは、画像保管通信システム(PACS)、放射線医学情報システム(RIS)、病院情報システム(HIS)及び/又は他のデータリポジトリのようなデータリポジトリから受信される及び/又は取り出されることができる。加えて又は代わりに、前記再構成画像、前記ノイズ画像及び/又は前記ノイズ除去された再構成画像は、前記データリポジトリに記憶されることができる。
加えて又は代わりに、投影データプロセッサ420、リビニングプロセッサ414、再構成プロセッサ416、及びノイズ除去プロセッサ422の少なくとも1つは、コンソール426及び/又は撮像システム400から離れた外部の計算装置の一部である。例えば、これらのコンポーネントの1以上が、「クラウド」ベースのサービスの一部であることができる、複数の装置にわたって分散される、他の撮像システムの一部である等であることができる。
図6は、一実施例による一例の方法をここに示す。
前記方法の動作の順序が限定的ではないと理解されるべきである。このようなものとして、他の順序が、ここで考えられる。加えて、1以上の動作が、省略されてもよく、及び/又は1以上の追加の動作が、含められてもよい。
602において、対象又はオブジェクトのスキャンからの投影データが、得られる。前記投影データは、スキャン中に撮像システム400及び/又は他の撮像システムにより生成され、そこから及び/又はデータリポジトリから得られることができる。
604及び606において、前記投影データは、ファンビームから平行ビームフォーマットに前記投影データをリビニングし、フィルタ逆投影により前記平行ビームデータを再構成し、前記スキャンされた対象又はオブジェクトの構造画像(体積画像データ)を生成することを含む従来の処理チェーンによって処理される。
608において、前記投影データは、また、ここに記載されるように、ビューにプラス1及びマイナス1を又は逆に乗算することを交互に行うことにより前処理され、前処理された投影データを生成する。
610及び612において、前記前処理された投影データは、ファンビームから平行ビームフォーマットに前記前処理された投影データをリビニングし、フィルタ逆投影により前記平行ビームデータを再構成することを含む従来の処理チェーンによって処理され、ノイズのみ画像を生成する。
614において、前記構造画像は、ここに記載されるように、前記ノイズ画像を使用してノイズ除去され、ノイズ除去された構造画像を生成する。
上記の方法は、コンピュータ可読記憶媒体上にエンコード又は埋め込まれたコンピュータ可読命令を用いて実施されてもよく、前記コンピュータ可読命令は、コンピュータプロセッサにより実行される場合に、前記プロセッサに、記載された動作を実行させる。加えて又は代わりに、前記コンピュータ可読命令の少なくとも1つは、信号、搬送波又は他の一時的媒体により運ばれる。
比較として、図7は、構造画像70、ここに記載されたアプローチを使用するノイズ画像704、及び図3に記載された従来のアプローチを使用するノイズ画像706を示す。ノイズ画像704及び706は、構造画像702の各画素/ボクセルに対するノイズ値を示す。ノイズ標準偏差は、構造画像702の関心領域(ROI)708の画素/ボクセル値に対して計算される。一般に、ROI708は、同じ又は同様のCT数(フラット)を持つとみられる同じ解剖学的組織の領域上に配置される。この例において、ROI708に対するノイズ標準偏差は、73である。
ノイズ画像704において、ROI710は、構造画像702内のROI708と同じ場所に配置される。従来技術のノイズ画像706において、ROI712も、構造画像702内のROI708と同じ場所に配置される。示されるように、ここに記載されたアプローチを使用することは、構造画像702のROI708に対するノイズ標準偏差に近く、従来のノイズ画像706のROI712に対するノイズ標準偏差と比較すると構造画像702のROI708に対するノイズ標準偏差により近いノイズ標準偏差を生成する。
本発明は、好適な実施例を参照して説明されている。修正例及び変更例は、先行する詳細な説明を読み、理解すると、他者が思いつきうる。本発明は、添付の請求項又はその同等物の範囲に入る限り、全てのこのような修正例及び変更例を含むと解釈されることを意図される。
以下、本発明の各種形態を付記する。
(付記1)
X線放射線を放射する放射線源と、
X線放射線を検出し、それを示す投影データを生成する検出器アレイであって、前記投影データが、複数のビューを含む、当該検出器アレイと、
前記投影データを再構成し、ノイズのみ画像を生成する第1の処理チェーンと、
を有する撮像システム。
(付記2)
前記第1の処理チェーンが、前記ビューに所定の乗算係数を乗算し、前処理された投影データを生成する投影データプロセッサを含む、付記1に記載の撮像システム。
(付記3)
前記所定の乗算係数が、前記前処理された投影データの1つおきのビューにマイナス1を乗算する、付記2に記載の撮像システム。
(付記4)
前記所定の乗算係数が、前記前処理された投影データの残りのビューにプラス1を乗算する、付記3に記載の撮像システム。
(付記5)
前記第1の処理チェーンが、前記前処理された投影データを平行ビームフォーマットにリビニングするリビニングプロセッサを含む、付記2乃至4のいずれか一項に記載の撮像システム。
(付記6)
前記第1の処理チェーンが、前記リビニングされた前処理された投影データを再構成し、前記ノイズのみ画像を生成する再構成プロセッサを含む、付記5に記載の撮像システム。
(付記7)
前記投影データを再構成し、構造及びノイズを含む体積画像データを生成する第2の処理チェーン、
を有する、付記1乃至6のいずれか一項に記載の撮像システム。
(付記8)
前記ノイズ画像を使用して前記体積画像データをノイズ除去し、ノイズ除去された体積画像データを生成するノイズ除去プロセッサ、
を有する、付記7に記載の撮像システム。
(付記9)
前記ノイズ除去プロセッサが、前記ノイズのみ画像からノイズ分散又はノイズ標準偏差の少なくとも1つを推定し、前記ノイズ分散又はノイズ標準偏差が、前記体積画像データをノイズ除去するのに使用される、付記8に記載の撮像システム。
(付記10)
撮像システムにより生成された投影データを受信するステップであって、前記投影データが複数のビューを含む、ステップと、
前記投影データを再構成してノイズのみ画像を生成する第1の処理チェーンを用いて前記投影データを処理するステップと、
を有する方法。
(付記11)
前記投影データを処理するステップが、
前処理された投影データを生成するように前記投影データのビューに所定の乗算係数を乗算するステップ、
を含む、付記10に記載の方法。
(付記12)
前記乗算するステップが、
偶数ビューにマイナス1又はプラス1の一方を乗算し、奇数ビューにマイナス1又はプラス1の他方を乗算するステップ、
を含む、付記11に記載の方法。
(付記13)
前記前処理された投影データを平行ビームフォーマットにリビニングするステップと、
前記リビニングされたデータを再構成し、前記ノイズのみ画像を生成するステップと、
を有する、付記11乃至12のいずれか一項に記載の方法。
(付記14)
前記投影データを再構成して構造の画像を生成する第2の処理チェーンを用いて前記投影データを処理するステップ、
を有する、付記10乃至12のいずれか一項に記載の方法。
(付記15)
前記投影データを平行ビームフォーマットにリビニングするステップと、
前記リビニングされたデータを再構成し、前記構造の画像を生成するステップと、
を含む、付記14に記載の方法。
(付記16)
関心領域の局所的ノイズ分散又はノイズ標準偏差を取り、これを画像中で移動させることにより前記ノイズのみ画像に基づいてノイズ分散又はノイズ標準偏差の少なくとも1つを推定するステップ、
を有する、付記13に記載の方法。
(付記17)
前記ノイズのみ画像に基づいて前記構造の画像をノイズ除去し、ノイズ除去された構造の画像を生成するステップ、
を有する、付記15乃至16のいずれか一項に記載の方法。
(付記18)
計算システムのプロセッサにより実行される場合に、前記プロセッサに、
X線撮像システムを用いてオブジェクト又は対象をスキャンさせ、複数のビューを含む投影データを生成させ、
前記投影データを再構成してノイズのみ画像を生成する第1の処理チェーンを用いて前記投影データを処理させ、
前記投影データを再構成して構造画像を生成する第2の処理チェーンを用いて前記投影データを処理させ、
前記ノイズのみ画像に基づいて前記構造画像をノイズ除去させる、
1以上のコンピュータ実行可能命令をエンコードされたコンピュータ可読記憶媒体。
(付記19)
前記プロセッサが、
前処理された投影データを生成するように前記投影データのビューに所定の乗算係数を乗算し、
前記前処理された投影データを平行ビームフォーマットにリビニングし、
前記ノイズのみ画像を生成するように前記リビニングされたデータを再構成する、
付記18に記載のコンピュータ可読記憶媒体。
(付記20)
前記プロセッサが、
前処理された投影データを生成するように前記投影データに交互打ち消しを使用する、
付記18に記載のコンピュータ可読記憶媒体。

Claims (11)

  1. X線放射線を放射する放射線源と、
    X線放射線を検出し、それを示す投影データを生成する検出器アレイであって、前記投影データが、複数のビューを含む、当該検出器アレイと、
    前記投影データを再構成し、ノイズのみ画像を生成する第1の処理チェーンと、
    を有する撮像システムであって、
    前記第1の処理チェーンが、前記ビューに所定の乗算係数を乗算し、前処理された投影データを生成する投影データプロセッサを含み、
    前記第1の処理チェーンが、前記前処理された投影データを平行ビームフォーマットにリビニングするリビニングプロセッサを含み、
    前記第1の処理チェーンが、前記リビニングされた前処理された投影データを再構成し、前記ノイズのみ画像を生成する再構成プロセッサを含む、
    撮像システム
  2. 前記第1の処理チェーンが、前記投影データの1つおきのビューにマイナス1を乗算し、前記投影データの残りのビューにプラス1を乗算する、請求項1に記載の撮像システム。
  3. 前記投影データを再構成し、構造及びノイズを含む体積画像データを生成する第2の処理チェーン、
    を有する、請求項1乃至2のいずれか一項に記載の撮像システム。
  4. 前記ノイズ画像を使用して前記体積画像データをノイズ除去し、ノイズ除去された体積画像データを生成するノイズ除去プロセッサ、
    を有する、請求項3に記載の撮像システム。
  5. 前記ノイズ除去プロセッサが、前記ノイズのみ画像からノイズ分散又はノイズ標準偏差の少なくとも1つを推定し、前記ノイズ分散又はノイズ標準偏差が、前記体積画像データをノイズ除去するのに使用される、請求項4に記載の撮像システム。
  6. 撮像システムにより生成された投影データを受信するステップであって、前記投影データが複数のビューを含む、ステップと、
    前記投影データの偶数ビューにマイナス1又はプラス1の一方を乗算し、前記投影データの奇数ビューに前記マイナス1又はプラス1の他方を乗算することにより、第1の処理チェーンを用いて前記投影データを処理するステップと、
    前記処理された投影データを再構成してノイズのみ画像を生成するステップと、
    前記前処理された投影データを平行ビームフォーマットにリビニングするステップと、
    前記リビニングされたデータを再構成し、前記ノイズのみ画像を生成するステップと、
    を有する方法。
  7. 前記投影データを再構成して構造の画像を生成する第2の処理チェーンを用いて前記投影データを処理するステップ、
    を有する、請求項6又は7に記載の方法。
  8. 前記投影データを平行ビームフォーマットにリビニングするステップと、
    前記リビニングされたデータを再構成し、前記構造の画像を生成するステップと、
    を含む、請求項7に記載の方法。
  9. 関心領域の局所的ノイズ分散又はノイズ標準偏差を取り、これを画像中で移動させることにより前記ノイズのみ画像に基づいてノイズ分散又はノイズ標準偏差の少なくとも1つを推定するステップ、
    を有する、請求項6に記載の方法。
  10. 前記ノイズのみ画像に基づいて前記構造の画像をノイズ除去し、ノイズ除去された構造の画像を生成するステップ、
    を有する、請求項8又は9に記載の方法。
  11. 計算システムのプロセッサにより実行される場合に、前記プロセッサに、
    撮像システムによって生成された、複数のビューを含む投影データを受け取るステップと、
    前記投影データの偶数ビューにマイナス1又はプラス1の一方を乗算し、前記投影データの奇数ビューに前記マイナス1又はプラス1の他方を乗算するよう構成された第1の処理チェーンを用いて前記投影データを処理するステップと、
    前記処理された投影データを再構成してノイズのみ画像を生成するステップと、
    前記前処理された投影データを平行ビームフォーマットにリビニングするステップと、
    前記リビニングされたデータを再構成し、前記ノイズのみ画像を生成するステップと、
    を実行させる、1以上のコンピュータ実行可能命令をエンコードされたコンピュータ可読記憶媒体。
JP2019530016A 2016-12-06 2017-12-06 交互打ち消しを使用する画像ノイズ推定 Active JP7123927B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201662430424P 2016-12-06 2016-12-06
US62/430,424 2016-12-06
PCT/EP2017/081593 WO2018104349A1 (en) 2016-12-06 2017-12-06 Image noise estimation using alternating negation

Publications (3)

Publication Number Publication Date
JP2019535451A JP2019535451A (ja) 2019-12-12
JP2019535451A5 JP2019535451A5 (ja) 2021-01-21
JP7123927B2 true JP7123927B2 (ja) 2022-08-23

Family

ID=60923439

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019530016A Active JP7123927B2 (ja) 2016-12-06 2017-12-06 交互打ち消しを使用する画像ノイズ推定

Country Status (5)

Country Link
US (1) US10984564B2 (ja)
EP (1) EP3552181B1 (ja)
JP (1) JP7123927B2 (ja)
CN (1) CN110073412B (ja)
WO (1) WO2018104349A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7395523B2 (ja) * 2021-02-02 2023-12-11 富士フイルムヘルスケア株式会社 医用画像処理装置および医用画像処理方法
TW202234049A (zh) * 2021-02-15 2022-09-01 日商濱松赫德尼古斯股份有限公司 放射線圖像獲取裝置、放射線圖像獲取系統及放射線圖像獲取方法
JP2024534861A (ja) 2021-08-31 2024-09-26 コーニンクレッカ フィリップス エヌ ヴェ 機械学習を使用した医用画像のノイズ除去

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090232269A1 (en) 2008-03-14 2009-09-17 Jiang Hsieh Methods and apparatus for noise estimation for multi-resolution anisotropic diffusion filtering
WO2015012323A1 (ja) 2013-07-25 2015-01-29 株式会社日立メディコ X線ct装置
JP2018000618A (ja) 2016-07-04 2018-01-11 株式会社日立製作所 画像再構成装置、x線ct装置、および、画像再構成方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6078639A (en) 1997-11-26 2000-06-20 Picker International, Inc. Real time continuous CT imaging
EP1599836B1 (en) 2003-02-14 2010-04-21 Koninklijke Philips Electronics N.V. System and method for helical cone-beam computed tomography with exact reconstruction
US7672421B2 (en) 2005-10-12 2010-03-02 Siemens Medical Solutions Usa, Inc. Reduction of streak artifacts in low dose CT imaging through multi image compounding
US7756312B2 (en) 2006-10-19 2010-07-13 General Electric Company Methods and apparatus for noise estimation
KR20120116940A (ko) 2009-11-25 2012-10-23 코닌클리케 필립스 일렉트로닉스 엔.브이. 향상된 화상 데이터/선량 감소
US8396273B2 (en) 2009-12-30 2013-03-12 General Electric Company Noise reduction method for dual-energy imaging
CN102947861B (zh) * 2010-06-21 2016-06-29 皇家飞利浦电子股份有限公司 用于在低剂量计算机断层摄影中降低噪声的方法和系统
US8761540B2 (en) 2011-06-14 2014-06-24 Kabushiki Kaisha Toshiba Method and system for estimating noise level
EP2783344B1 (en) 2011-11-23 2017-05-03 Koninklijke Philips N.V. Image domain de-noising
US9031297B2 (en) 2013-02-01 2015-05-12 Kabushiki Kaisha Toshiba Alternative noise map estimation methods for CT images
US9076237B2 (en) 2013-03-12 2015-07-07 Wisconsin Alumni Research Foundation System and method for estimating a statistical noise map in x-ray imaging applications
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
EP3014573B1 (en) 2013-06-28 2017-12-13 Koninklijke Philips N.V. Methods of utilizing image noise information
US9357976B2 (en) * 2013-10-24 2016-06-07 General Electric Company System and method of noise deletion in computed tomography
JPWO2016132880A1 (ja) * 2015-02-16 2017-11-30 株式会社日立製作所 演算装置、x線ct装置、及び画像再構成方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090232269A1 (en) 2008-03-14 2009-09-17 Jiang Hsieh Methods and apparatus for noise estimation for multi-resolution anisotropic diffusion filtering
WO2015012323A1 (ja) 2013-07-25 2015-01-29 株式会社日立メディコ X線ct装置
JP2018000618A (ja) 2016-07-04 2018-01-11 株式会社日立製作所 画像再構成装置、x線ct装置、および、画像再構成方法

Also Published As

Publication number Publication date
US20190385345A1 (en) 2019-12-19
CN110073412A (zh) 2019-07-30
WO2018104349A1 (en) 2018-06-14
EP3552181B1 (en) 2021-11-10
EP3552181A1 (en) 2019-10-16
CN110073412B (zh) 2023-11-03
US10984564B2 (en) 2021-04-20
JP2019535451A (ja) 2019-12-12

Similar Documents

Publication Publication Date Title
US9959640B2 (en) Iterative image reconstruction with a sharpness driven regularization parameter
CN103649990B (zh) 用于谱ct的图像处理
JP6169558B2 (ja) コントラスト依存の解像度をもつ画像
US9031299B2 (en) Low dose CT denoising
US7515676B2 (en) Method, apparatus, and computer program product for sinogram completion
US8938108B2 (en) Method for artifact reduction in cone-beam CT images
CN102270350A (zh) 结合四维噪声滤波器的迭代ct图像重建
WO2010128413A1 (en) System and method for generating a tomographic reconstruction filter
JP2011125700A (ja) X線画像の時間分解能を高めるシステム及び方法
US11238625B2 (en) Imaging system extended field-of-view
JP7123927B2 (ja) 交互打ち消しを使用する画像ノイズ推定
JP2021511875A (ja) スペクトルボリューム画像データを生成するように構成された非スペクトルコンピュータ断層撮影(ct)スキャナ
US9996951B2 (en) Computerized tomographic image exposure and reconstruction method
EP3513375A1 (en) Edge noise reduction
JP2018143574A (ja) X線ct装置及び画像処理方法
JP2021502842A (ja) 複数の異なるタイプの投影データのボクセル毎に1回のジオメトリ演算を行う単一のct逆投影器
WO2018050462A1 (en) Edge noise reduction

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201201

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20211013

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211019

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20211117

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220330

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220810

R150 Certificate of patent or registration of utility model

Ref document number: 7123927

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150