JP5878119B2 - X線ct装置、及びその制御方法 - Google Patents
X線ct装置、及びその制御方法 Download PDFInfo
- Publication number
- JP5878119B2 JP5878119B2 JP2012520486A JP2012520486A JP5878119B2 JP 5878119 B2 JP5878119 B2 JP 5878119B2 JP 2012520486 A JP2012520486 A JP 2012520486A JP 2012520486 A JP2012520486 A JP 2012520486A JP 5878119 B2 JP5878119 B2 JP 5878119B2
- Authority
- JP
- Japan
- Prior art keywords
- filter
- ray
- projection data
- convolution
- value
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating 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/02—Investigating 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/04—Investigating 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/046—Investigating 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]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Description
まず、図1-1及び図1-2に基づいて、本発明が適用されるコーンビームX線CT装置の概略的な構成について説明する。図1-1は、本発明が適用されるコーンビームX線CT装置(Cアーム方式)1を示す概略構成図である。図1-2は、本発明が適用される移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置1aを示す概略構成図である。
撮影部10は、寝台17と、該寝台17に横臥された被検体2にX線を照射するX線源11と、該X線源11に対向して設置され被検体2を透過したX線を検出することによりX線透過像111を出力する2次元X線検出器12と、X線源11及び2次元X線検出器12を機械的に接続するC型アーム13と、該C型アーム13を保持するC型アーム保持体14と、該C型アーム保持体14を天井に取り付ける天井支持体15と、該天井支持体15を図示の状態で前後左右の2次元方向に移動可能に支持する天井レール16と、被検体2に造影剤を注入するインジェクタ18と、を備える。
制御演算部20は、撮影部10を制御する撮影部制御手段100と、撮影部10が出力したX線透過像111を収集して格納する画像収集手段110と、収集されたX線透過像111に基づいて3次元CT像を再構成する再構成手段200と、再構成手段200が生成した3次元CT像を表示する画像表示手段280と、再構成手段200がコンボリューションフィルタを生成するために使用する生成条件を入力するコンボリューションフィルタ入力手段320と、再構成手段200がフーリエ変換(以下Fast Fourier Transformを略してFFTと記載する)フィルタを生成するために使用する生成条件を入力するFFTフィルタ入力手段340と、を備える。また、前記コンボリューションフィルタとは、画像空間上でコンボリューション演算を用いてスムージングや、鮮鋭化等の画像処理を行なう際に、当該画素値及びその周辺の画素値に対し重畳する係数である。また、FFTフィルタとは、2次元の画像空間に対し1ラインごと(1次元)にFFT変換し、該FFT変換により生成したラインデータにおいて、周波数ごとに重畳する係数である。また、前記コンボリューションフィルタを生成するために使用する生成条件及び、前記FFTフィルタを生成するために使用する生成条件の詳細については、図5、図6などを用い後述する。
撮影部制御手段100は、C型アーム13の、回転中心軸4の回りの回転移動を制御する撮影系回転制御手段101と、天井支持体15の天井レール16上での位置を制御してC型アーム13の被検体2に対する位置を2次元的に制御する撮影系位置制御手段102と、X線管11tに流す管電流のON、OFF等を制御するX線照射制御手段103と、インジェクタ18が被検体2に注入する造影剤の注入量及び注入タイミングを制御するインジェクタ制御手段104と、寝台17の位置を制御して被検体2の位置を調整するための寝台制御手段105と、2次元X線検出器12によるX線透過像111の撮影を制御する検出系制御手段107と、を備える。なお、C型アーム13の回転方向は、前述したように、紙面に平行な方向に回転中心軸4が有りX線源11と2次元X線検出器12が回転中心軸4を中心に旋回するようになっていても(図1-1)、回転中心軸4が紙面と垂直な方向に存在しX線源11と2次元X線検出器12が紙面と平行な面内をスライド回転するようになっていても(図1-2)、あるいはその両方の回転動作を備えるものであってもよい。
再構成手段200は、前処理手段210と、フィルタ変換情報生成手段220と、平滑化手段230と、フィルタリング手段250と、逆投影手段260と、を備える。
コンボリューションフィルタ入力手段320は、コンボリューションフィルタ生成手段233が生成する2次元のコンボリューションフィルタの生成条件を設定する。以下、図5に基づいて、コンボリューションフィルタ入力手段320が使用するGUIの一例を説明する。図5は、コンボリューションフィルタ設定画面30の一例を示す模式図である。
FFTフィルタ入力手段340は、FFTフィルタ生成手段253が生成するFFTフィルタの生成条件を設定する。以下、図6に基づいて、FFTフィルタ入力手段340が使用するGUIの一例を説明する。図6は、FFTフィルタ設定画面40の一例を示す模式図である。図5と同様に、タブ41〜44は撮影部位選択タブ、ボタン45は撮影部位追加ボタンである。図6は胸部条件設定タブ42が選択された場合を表している。
リストボックス51は高空間分解能領域用フィルタ関数を選択するリストボックス、リストボックス52は低コントラスト領域用フィルタ関数を選択するリストボックスであり、FFTフィルタ生成手段253が使用する生成条件を選択する。ポイント53はFFTフィルタ関数閾値設定ポイント、ポイント54はFFTフィルタ関数変化量設定ポイントであり、左右にドラッグしてFFTフィルタ関数閾値μF、あるいは閾値境界線の傾きを変化させることによりFFTフィルタ関数変化量βFを変化できるようになっている。ここでいう「FFTフィルタ関数閾値μF」とは、FFTフィルタを用いた平滑化処理の大小を規定する値である。また、「FFTフィルタ関数変化量βF」とは、FFTフィルタ処理に用いられる二つのフィルタ関数(図6におけるフィルタ関数1とフィルタ関数2)との変化の割合を規定する値である。上記のフィルタ関数閾値μF、フィルタ関数変化量βF、及び選択した高空間分解能領域用フィルタ関数、低コントラスト領域用フィルタ関数を用いて、FFTフィルタ生成手段253がFFTフィルタを生成する処理の詳細については後述する。
次に、コーンビームX線CT装置1による撮影における動作の概要について説明する。
次に、再構成手段200により実行される再構成処理(S200)の流れを、図7に基づいて説明する。図7は、再構成手段200により実行される再構成処理(S200)の流れを示すフローチャートである。以下、図7のステップ順に沿って説明する。
前処理手段210は、画像収集手段110が収集した、被検体2と空気のX線透過像111に自然対数変換演算を施し、投影データ211に変換する(S210)。
フィルタ変換情報生成手段220は、投影データ211の各点近傍の特徴量(例えば、各点近傍の画素の画素値の平均値Acや標準偏差σc)を計算し、ステップS230でコンボリューションフィルタを生成するための、コンボリューションフィルタ変換パラメータを生成する(S220)。
平滑化手段230は、ステップS220が生成したコンボリューションフィルタ変換パラメータと、予めユーザがコンボリューションフィルタ入力手段320を用いて設定したコンボリューションフィルタの生成条件と、を用いて、投影データ211の各点ごとにコンボリューションフィルタを生成し、各点ごとに2次元のコンボリューション演算を施す(S230)。
フィルタ変換情報生成手段220は、コンボリューション演算処理を施した投影データ212の各横ラインデータの画素値の特徴量(平均値AF、標準偏差σF)を計算し、ステップS250でFFTフィルタを生成するための、FFTフィルタ変換パラメータを生成する(S240)。
フィルタリング手段250は、ステップS240が生成したFFTフィルタ変換パラメータと、予めユーザがFFTフィルタ入力手段340を用いて設定したFFTフィルタの生成条件と、を用いて、コンボリューション演算処理を施した投影データ212の各横ラインに対応したFFTフィルタを生成し、FFTフィルタ処理を施す(S250)。
逆投影手段260は、ステップS250によるFFTフィルタ処理後の投影データを用いて、逆投影演算を行う(S260)。
すべての投影データについて、ステップS210からステップS260までの処理を実行したかを判別する。すべての投影データを処理していない(NOの)場合は、ステップS210に戻り、次の投影データについてステップS210からステップS260までの処理を実行する。YESの場合は、再構成処理(S200)を終え、被検体2の3次元CT像を出力する(S270)。
投影データ読込み手段221は、ステップS210で生成された投影データ211を読込む(S221)。
ROI設定手段222は、図9に示すように、投影データ211の特徴量を計算するための、長方形あるいは正方形のROI(計算領域239)の大きさを設定する。ROIサイズは、座標点238を中心に、例えば、15×15〜25×25画素とする(S222)。
特徴量算出手段223は、投影データ211上の座標点を走査し、ステップS222で指定されたROI領域(計算領域239)内の画素値の平均値Ac及び標準偏差σcを計算することにより、投影データ211の各点の特徴量を求める。このとき、平均値から大きく外れたデータにより、標準偏差の計算値が影響を受けるのを抑えるために、一度ROI毎に、そのROI内の画素値の平均値、最大値、最小値を横軸とするヒストグラムを作成する。そして、ROI内の画素値の平均値から例えば±1/10の度数分布内にあるデータのみを使用して、標準偏差値を求めるようにしてもよい(S223)。なお、本実施形態では、投影データ211の各点について特徴量を求めたが、全ての点について特徴量を求めることなく、適宜間引きした点についてのみ特徴量を求めてもよい。その場合、間引きされた投影データ211の点の特徴量は、その点を含む計算領域から得られた画素値の特徴量を流用してもよい。
特徴量フィティング手段224は、ステップS223が計算した投影データ211各点の特徴量を、投影データの座標値の関数としてフィティング(すなわち、投影データ211の各点の横方向及び縦方向の位置情報(座標値)と、その位置における画素値の特徴量と、の対応付け)する(S224)。なお、ステップS223及びS224を省略することも可能である。この場合、例えば、撮影部位や撮影条件を基に、ユーザが画素値の特徴量(または特徴量の範囲)の予測値を入力し、その入力された画素値の特徴量を用いて、下記のフィルタ変換情報を算出してもよい。
フィルタ変換情報算出手段225は、ステップS224のフィティング結果をコンボリューションフィルタのパラメータに変換、すなわち、フィティング結果に含まれる画素値の特徴量を用いてコンボリューションフィルタのパラメータを生成する(S225)。このコンボリューションフィルタのパラメータ及びその処理の詳細について、以下図12に基づいて説明する。図12は、投影データ211の標準偏差σcとパラメータWaの関数が示す曲線を示す説明図である。
フィルタ変換情報保存手段226は、ステップS225が求めたフィルタ変換情報(コンボリューションフィルタ変換パラメータを投影データの座標値の関数とする情報)を保存する(S226)。
以下図13のステップ順に沿って説明する。
コンボリューションフィルタ変換情報読込み手段231は、コンボリューションフィルタ変換情報生成処理(S220)が生成したコンボリューションフィルタ変換情報(具体的には、ステップS226により投影データの座標値の関数として保存された、平滑化パラメータWa及びWb)を読込む(S231)。
画素スキャン手段232は、投影データ211の各点(以下、投影画素という)の座標値(画素の座標と画素値とに相当する)を走査する(S232)。
コンボリューションフィルタ生成手段233は、コンボリューションフィルタ変換情報を用いて、各投影画素の座標値に対応する平滑化パラメータWa、Wbを読み出し、これらを用いて各投影画素に対応するコンボリューションフィルタを生成する(S233)。
周辺画素読込み手段234は、コンボリューション演算の対象となる投影画素の周辺に位置する投影画素(周辺画素ともいう)の画素値を読み取る(S234)。周辺に位置する投影画素は、以下のコンボリューション演算に用いられるものである。
コンボリューション手段235は、ステップS233で生成されたコンボリューションフィルタと、ステップS232,S234で読み込まれた投影画素及び周辺画素の画素値を用いて、コンボリューション演算を行う(S235)。コンボリューション手段235は、例えば、図5のコンボリューションフィルタ設定画面において3×3のフィルタサイズが設定されると、式(3)に基づいて、コンボリューション演算を行う。
コンボリューション演算処理後投影データ保存手段236は、ステップS236のコンボリューション演算結果を保存する(S236)。
すべての投影画素について、ステップS232からステップS236までの処理を実行したかが判別される。すべての投影画素を処理していない(NOの)場合は、ステップS232に戻り、次の投影画素についてステップS232からステップS236までの処理を実行する。YESの場合は、平滑化処理(S230)を終え、FFTフィルタ変換情報生成処理(S240)へ進む。
投影データ読込み手段221は、ステップS230で生成されたコンボリューション演算処理後投影データ212を読込む(S241)。
ROI設定手段222は、図15に示すような、横長の長方形ROI(計算領域342)の大きさを設定する(S242)。
特徴量算出手段223は、長方形ROI(計算領域342)を、コンボリューション演算処理後投影データ212上で上下方向(又は縦方向ともいう)に走査し、特徴量(例えばROI内の画素の画素値の平均値AF、標準偏差σF)を計算する。このとき、ステップS223の計算のときと同様に、平均値から大きく外れたデータにより、標準偏差の計算値が影響を受けるのを抑えるために、一度、ROI内の画素値の平均値、最大値、最小値を横軸とするヒストグラムを作成し、ROI内の画素値の平均値から例えば±1/10の度数分布内にあるデータのみを使用して、標準偏差値を求めるようにしてもよい(S243)。
特徴量フィティング手段224は、ステップS243が計算した特徴量を、投影データの上下方向に関する座標値の関数としてフィティング(すなわち、投影データの上下方向の位置情報と、その位置における特徴量と、の対応付け)する(S244)。
フィルタ変換情報算出手段225は、ステップS244のフィティング結果をFFTフィルタのパラメータに変換する。このFFTフィルタのパラメータ及びその処理の詳細については、以下に説明する(S245)。
フィルタ変換情報保存手段226は、ステップS245で求められたFFTフィルタ変換情報を保存する(S246)。
FFTフィルタ変換情報読込み手段251は、FFTフィルタ変換情報生成処理(S240)が生成したFFTフィルタ変換情報(具体的には、ステップS246により投影データの上下方向に関する座標値の関数として保存された、高空間分解能領域用FFTフィルタ関数含有率)を読込む(S251)。
ラインデータ読込み手段252は、コンボリューション演算処理後投影データ212から一度にフィルタリング処理する横ラインデータ、例えば、図15における横ラインデータ352を読込む(S252)。
FFTフィルタ生成手段253は、各ラインデータに対し、そのラインデータの上下方向の座標値と、ステップS251で読み込んだFFTフィルタ変換情報と、を用いてラインデータの上下座標に対応する
FFTフィルタを生成する(S253)。FFTフィルタ生成手段253は、FFTフィルタ変換情報に基づいてラインデータの上下方向の座標値に対応する高空間分解能領域用FFTフィルタ関数含有率を算出する。そして、図6のGUI画面によりユーザが設定入力した高空間分解能領域用FFTフィルタ関数を、算出された含有率に従って含み、同じく図6のGUI画面によりユーザが設定入力した低コントラスト領域用のFFTフィルタ関数を、(1−上記含有率)に従って含むFFTフィルタを生成する(S253)。
FFT手段254は、横ラインデータ352を周波数データに変換する(S254)。
FFTフィルタ積算手段255は、周波数データにステップS253が生成したFFTフィルタを積算する(S255)。
逆FFT手段256は、周波数データを実空間のフィルタ処理後投影データに変換する(S256)。
フィルタ処理後投影データ保存手段257は、ステップS256が生成したフィルタ処理後投影データを保存する(S257)。
すべての横ラインデータについて、ステップS252からステップS257までの処理を実行したかを判別する。すべてのラインデータを処理していない(NOの)場合は、ステップS252に戻り、次のラインデータについてステップS252からステップS257までの処理を実行する。YESの場合は、フィルタリング処理(S250)を終え、逆投影処理(S260)へ進む。
Claims (10)
- X線を発生させるX線源と、
前記X線源に対向して配置され、被検体を透過した前記X線を検出して前記被検体の投影データを出力するX線検出器と、
前記X線源と前記X線検出器とを対向配置したまま回転させる回転手段と、
前記投影データに含まれる画素値の特徴量に応じて変化する画像処理フィルタを生成するフィルタ生成手段と、
前記投影データに対し、前記生成された画像処理フィルタを用いて再構成演算を行い、前記被検体のX線CT像を生成する再構成手段と、
前記X線CT像を表示する画像表示手段と、を備え、
前記画素値の特徴量は、前記投影データの各点の画素値と、前記各点の近傍の画素値との平均値と標準偏差から計算される、
ことを特徴とするX線CT装置。 - 前記フィルタ生成手段は、前記投影データの各点に対し、当該点の画素値の特徴量に対応した前記画像処理フィルタとしてコンボリューションフィルタを生成し、
前記再構成手段は、前記各点毎に、前記コンボリューションフィルタを重畳するコンボリューション手段を備える、
ことを特徴とする請求項1に記載のX線CT装置。 - 前記X線CT装置は、前記コンボリューションフィルタの生成条件を入力する第一入力手段と、
前記投影データの画素値の特徴量に対応して変化するパラメータを含む第一フィルタ変換情報を生成する第一フィルタ変換情報生成手段と、を更に備え、
前記フィルタ生成手段は、前記入力された生成条件及び前記第一フィルタ変換情報を用い、前記コンボリューションフィルタを生成する、
ことを特徴とする請求項2に記載のX線CT装置。 - 前記第一入力手段は、前記生成条件として、前記コンボリューションフィルタの横方向コンボリューションサイズ、横方向の平滑化処理の大小を規定する横方向フィルタ関数閾値、及び横方向の平滑化フィルタ関数の変化量を規定する横方向フィルタ関数変化量と、前記コンボリューションフィルタの縦方向コンボリューションサイズ、縦方向の平滑化処理の大小を規定する縦方向フィルタ関数閾値、及び縦方向の平滑化フィルタ関数の変化量を規定する縦方向フィルタ関数変化量と、の入力を受付け、
前記第一フィルタ変換情報生成手段は、前記投影データの画素値の特徴量と、前記横方向フィルタ関数閾値及び横方向フィルタ関数変化量と、を用いて横方向平滑化パラメータを算出するとともに、前記投影データの画素値の特徴量と、前記縦方向フィルタ関数閾値及び縦方向フィルタ関数変化量と、を用いて縦方向平滑化パラメータを算出し、
前記フィルタ生成手段は、前記投影データの各点に対し、当該点の画素値の特徴量と、その特徴量に応じた前記横方向平滑化パラメータ及び前記縦方向平滑化パラメータと、を用い、前記入力された横方向コンボリューションサイズ及び縦方向コンボリューションサイズからなる前記コンボリューションフィルタを生成する、
ことを特徴とする請求項3に記載のX線CT装置。 - 前記第一フィルタ変換情報生成手段は、前記投影データの各点を含む第一計算領域を設定する第一領域設定手段と、前記第一計算領域に含まれる画素の画素値の特徴量を算出する第一特徴量算出手段と、前記各点の座標値と前記画素値の特徴量とを対応付ける第一特徴量フィティング手段と、を備え、前記座標値に対応付けられた画素値の特徴量を用いて算出された前記横方向平滑化パラメータ及び前記縦方向平滑化パラメータと、前記座標値と、が対応付けられた前記第一フィルタ変換情報を生成し、
前記フィルタ生成手段は、前記投影データの各点の座標値と、前記第一フィルタ変換情報と、を用いて前記コンボリューションフィルタを生成する、
ことを特徴とする請求項4に記載のX線CT装置。 - 前記フィルタ生成手段は、前記投影データ、又は前記コンボリューションフィルタが重畳されたコンボリューション演算処理後投影データの、1次元フーリエ変換処理の対象となる各ラインデータに対し、当該ラインデータの画素値の特徴量に基づいて前記画像処理フィルタとして周波数に対応したFFTフィルタを生成し、
前記再構成手段は、前記投影データの各ラインデータ、又は前記コンボリューション演算処理後投影データの各ラインデータ毎に、前記FFTフィルタを用いたフィルタリング処理を行うフィルタリング手段を備える、
ことを特徴とする請求項2に記載のX線CT装置。 - 前記X線CT装置は、前記FFTフィルタの生成条件を入力する第二入力手段と、
前記投影データのラインデータの画素値の特徴量に対応して変化するパラメータを含む第二フィルタ変換情報を生成する第二フィルタ変換情報生成手段と、を更に備え、
前記フィルタ生成手段は、前記入力された生成条件及び前記第二フィルタ変換情報を用い、前記投影データの各ラインデータ、又は前記コンボリューション演算処理後投影データの各ラインデータに対し、前記FFTフィルタを生成する、
ことを特徴とする請求項6に記載のX線CT装置。 - 前記第二フィルタ変換情報生成手段は、前記投影データ又は前記コンボリューション演算処理後投影データのラインデータからなる第二計算領域を設定する第二領域設定手段と、前記第二計算領域に含まれる画素の画素値の特徴量を計算する第二特徴量算出手段と、前記ラインデータの長手方向に直交する方向における座標値と、前記画素値の特徴量と、を対応付ける第二特徴量フィティング手段と、を備え、前記座標値に対応付けられた画素値の特徴量を用いて算出されたパラメータと、前記座標値と、が対応付けられた第二フィルタ変換情報を生成し、
前記フィルタ生成手段は、前記各ラインデータの座標値と、前記第二フィルタ変換情報
と、を用いて前記FFTフィルタを生成する、
ことを特徴とする請求項7に記載のX線CT装置。 - 前記第二入力手段は、前記生成条件として、相対的に高空間分解能な領域に用いられる第一フィルタ関数と、相対的に低コントラスト領域に用いられる第二フィルタ関数と、前記第一フィルタ関数又は第二フィルタ関数の含有率の大小を規定するフィルタ関数閾値と、前記含有率の変化量を示すフィルタ関数変化量と、の入力を受け付け、前記第二フィルタ変換情報生成手段は、前記第二計算領域の画素値の特徴量に応じた前記含有率からなる前記パラメータを算出し、前記フィルタ生成手段は、前記第二フィルタ変換情報に基づいて、前記各ラインデータの前記第一フィルタ関数または前記第二フィルタ関数の含有率を算出し、その含有率に従って前記第一フィルタ関数と前記第二フィルタ関数とを積算することにより、前記FFTフィルタを生成する、ことを特徴とする請求項8に記載のX線CT装置。
- X線を発生させるX線源と、前記X線源に対向して配置され、被検体を透過した前記X線を検出して前記被検体の投影データを出力するX線検出器と、
前記X線源と前記X線検出器とを対向配置したまま回転させる回転手段と、を備えたX線CT装置の制御方法において、
前記投影データに含まれる画素値の特徴量に応じて変化する画像処理フィルタを生成するステップと、前記投影データに対し、前記生成された画像処理フィルタを用いて再構成演算を行い、前記被検体のX線CT像を生成するステップと、前記X線CT像を表示するステップと、を備え、
前記画素値の特徴量は、前記投影データの各点の画素値と、前記各点の近傍の画素値との平均値と標準偏差から計算される、
ことを特徴とするX線CT装置の制御方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012520486A JP5878119B2 (ja) | 2010-06-17 | 2011-06-16 | X線ct装置、及びその制御方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010138576 | 2010-06-17 | ||
JP2010138576 | 2010-06-17 | ||
JP2012520486A JP5878119B2 (ja) | 2010-06-17 | 2011-06-16 | X線ct装置、及びその制御方法 |
PCT/JP2011/063766 WO2011158893A1 (ja) | 2010-06-17 | 2011-06-16 | X線ct装置、及びその制御方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2011158893A1 JPWO2011158893A1 (ja) | 2013-08-19 |
JP5878119B2 true JP5878119B2 (ja) | 2016-03-08 |
Family
ID=45348286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012520486A Expired - Fee Related JP5878119B2 (ja) | 2010-06-17 | 2011-06-16 | X線ct装置、及びその制御方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9025848B2 (ja) |
JP (1) | JP5878119B2 (ja) |
CN (1) | CN102946807B (ja) |
WO (1) | WO2011158893A1 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5878119B2 (ja) | 2010-06-17 | 2016-03-08 | 株式会社日立メディコ | X線ct装置、及びその制御方法 |
US9177397B2 (en) * | 2011-06-21 | 2015-11-03 | Koninklijke Philips N.V. | Imaging apparatus |
JP6122269B2 (ja) * | 2011-12-16 | 2017-04-26 | キヤノン株式会社 | 画像処理装置、画像処理方法、及びプログラム |
JP6139821B2 (ja) * | 2012-03-22 | 2017-05-31 | 東芝メディカルシステムズ株式会社 | X線ct装置 |
JP5805689B2 (ja) * | 2013-03-08 | 2015-11-04 | 株式会社モリタ製作所 | X線ct撮影装置及びx線ct撮影方法 |
JP6165511B2 (ja) * | 2013-06-12 | 2017-07-19 | 東芝メディカルシステムズ株式会社 | X線コンピュータ断層撮影装置 |
US9224216B2 (en) * | 2013-07-31 | 2015-12-29 | Kabushiki Kaisha Toshiba | High density forward projector for spatial resolution improvement for medical imaging systems including computed tomography |
JP6283875B2 (ja) * | 2013-09-05 | 2018-02-28 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置、x線診断装置およびx線コンピュータ断層撮影装置 |
CN104142353B (zh) * | 2014-06-06 | 2017-01-11 | 清华大学 | 可见光成像辅助定位局部ct扫描方法及系统 |
JP6386981B2 (ja) * | 2015-08-31 | 2018-09-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像処理方法、画像処理装置及び放射線断層撮影装置並びにプログラム |
CN111462266B (zh) * | 2020-03-20 | 2023-08-04 | 北京东软医疗设备有限公司 | 图像重建方法、装置、ct设备及ct系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07136157A (ja) * | 1993-07-02 | 1995-05-30 | General Electric Co <Ge> | コンピュータ断層撮影イメージング装置 |
JP2003159245A (ja) * | 2001-09-05 | 2003-06-03 | Siemens Ag | 適応フィルタリングのための方法及び装置 |
JP2004000585A (ja) * | 2002-04-29 | 2004-01-08 | Ge Medical Systems Global Technology Co Llc | ノイズ軽減のための適応投影フィルタ処理スキーム |
JP2004070793A (ja) * | 2002-08-08 | 2004-03-04 | Ge Medical Systems Global Technology Co Llc | 3次元空間フィルタ装置および方法 |
JP2005160544A (ja) * | 2003-11-28 | 2005-06-23 | Toshiba Corp | 画像処理装置 |
WO2005110232A1 (ja) * | 2004-05-19 | 2005-11-24 | Hitachi Medical Corporation | 画像処理装置及びその方法 |
JP2006034785A (ja) * | 2004-07-29 | 2006-02-09 | Ge Medical Systems Global Technology Co Llc | X線ct画像処理方法およびx線ct装置 |
JP2007190182A (ja) * | 2006-01-19 | 2007-08-02 | Ge Medical Systems Global Technology Co Llc | 画像表示装置およびx線ct装置 |
JP2008148970A (ja) * | 2006-12-19 | 2008-07-03 | Ge Medical Systems Global Technology Co Llc | X線ctシステム |
JP4348989B2 (ja) * | 2003-04-15 | 2009-10-21 | 株式会社島津製作所 | 断層再構成装置およびそれを用いた断層撮影装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1223553A3 (en) * | 2000-10-17 | 2003-09-24 | Fuji Photo Film Co., Ltd. | Apparatus for suppressing noise by adapting filter characteristics to input image signal based on characteristics of input image signal |
JP4360817B2 (ja) * | 2002-10-18 | 2009-11-11 | 株式会社日立メディコ | 放射線断層撮影装置 |
WO2005077278A1 (ja) * | 2004-02-16 | 2005-08-25 | Hitachi Medical Corporation | 断層撮影像の再構成方法及び断層撮影装置 |
US8045776B2 (en) * | 2007-03-06 | 2011-10-25 | General Electric Company | Geometry-dependent filtering in CT method and apparatus |
DE102009007680A1 (de) * | 2008-08-13 | 2010-07-01 | CT Imaging GmbH | Verfahren und Vorrichtung zur interaktiven CT-Rekonstruktion |
US8244013B2 (en) * | 2008-09-10 | 2012-08-14 | Siemens Medical Solutions Usa, Inc. | Medical image data processing and interventional instrument identification system |
JP5878119B2 (ja) | 2010-06-17 | 2016-03-08 | 株式会社日立メディコ | X線ct装置、及びその制御方法 |
-
2011
- 2011-06-16 JP JP2012520486A patent/JP5878119B2/ja not_active Expired - Fee Related
- 2011-06-16 US US13/703,448 patent/US9025848B2/en active Active
- 2011-06-16 CN CN201180029701.2A patent/CN102946807B/zh not_active Expired - Fee Related
- 2011-06-16 WO PCT/JP2011/063766 patent/WO2011158893A1/ja active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07136157A (ja) * | 1993-07-02 | 1995-05-30 | General Electric Co <Ge> | コンピュータ断層撮影イメージング装置 |
JP2003159245A (ja) * | 2001-09-05 | 2003-06-03 | Siemens Ag | 適応フィルタリングのための方法及び装置 |
JP2004000585A (ja) * | 2002-04-29 | 2004-01-08 | Ge Medical Systems Global Technology Co Llc | ノイズ軽減のための適応投影フィルタ処理スキーム |
JP2004070793A (ja) * | 2002-08-08 | 2004-03-04 | Ge Medical Systems Global Technology Co Llc | 3次元空間フィルタ装置および方法 |
JP4348989B2 (ja) * | 2003-04-15 | 2009-10-21 | 株式会社島津製作所 | 断層再構成装置およびそれを用いた断層撮影装置 |
JP2005160544A (ja) * | 2003-11-28 | 2005-06-23 | Toshiba Corp | 画像処理装置 |
WO2005110232A1 (ja) * | 2004-05-19 | 2005-11-24 | Hitachi Medical Corporation | 画像処理装置及びその方法 |
JP2006034785A (ja) * | 2004-07-29 | 2006-02-09 | Ge Medical Systems Global Technology Co Llc | X線ct画像処理方法およびx線ct装置 |
JP2007190182A (ja) * | 2006-01-19 | 2007-08-02 | Ge Medical Systems Global Technology Co Llc | 画像表示装置およびx線ct装置 |
JP2008148970A (ja) * | 2006-12-19 | 2008-07-03 | Ge Medical Systems Global Technology Co Llc | X線ctシステム |
Also Published As
Publication number | Publication date |
---|---|
US9025848B2 (en) | 2015-05-05 |
CN102946807A (zh) | 2013-02-27 |
US20130094739A1 (en) | 2013-04-18 |
CN102946807B (zh) | 2015-10-14 |
JPWO2011158893A1 (ja) | 2013-08-19 |
WO2011158893A1 (ja) | 2011-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5878119B2 (ja) | X線ct装置、及びその制御方法 | |
US8965078B2 (en) | Projection-space denoising with bilateral filtering in computed tomography | |
US8280135B2 (en) | System and method for highly attenuating material artifact reduction in x-ray computed tomography | |
JP4537129B2 (ja) | トモシンセシス用途における対象物を走査するためのシステム | |
US9532755B2 (en) | X-ray computed tomographic imaging apparatus and method for same | |
JP2006043431A (ja) | ヘリカルマルチスライスctのための回復ノイズを伴うヘリカルウィンドミルアーチファクトを低減する方法 | |
JP2007181623A (ja) | X線ct装置 | |
JP5324883B2 (ja) | Ct装置および金属形状抽出方法 | |
JP3897925B2 (ja) | コーンビームct装置 | |
JP5669799B2 (ja) | 画像処理装置、放射線画像撮影システム、画像処理プログラム、及び画像処理方法 | |
JP3540916B2 (ja) | 3次元x線ct装置 | |
JP3583554B2 (ja) | コーンビームx線断層撮影装置 | |
JP5214110B2 (ja) | X線ct装置 | |
US11288775B2 (en) | Methods and systems for parametric noise modulation in x-ray imaging | |
JP2007175154A (ja) | X線ct装置 | |
JP7345292B2 (ja) | X線トモシンセシス装置、画像処理装置、および、プログラム | |
JP2009279301A (ja) | コーンビームx線ct装置 | |
JP5452841B2 (ja) | X線ct装置 | |
US20230145920A1 (en) | Systems and methods for motion detection in medical images | |
JP2001212135A (ja) | 放射線断層撮影装置および記憶媒体 | |
JP2003033348A (ja) | 3次元x線ct装置 | |
JP6199118B2 (ja) | X線ct装置及び医用画像処理方法 | |
JP2002204796A (ja) | 3次元x線ct装置 | |
JP2007275125A (ja) | X線ct装置およびそのx線ct画像再構成方法、x線ct画像撮影方法 | |
WO2018116791A1 (ja) | 医用画像処理装置及びそれを備えたx線ct装置、医用画像処理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140519 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150526 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150724 |
|
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: 20160105 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160127 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5878119 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |