JP2010005373A - 放射線画像補正方法および装置ならびにプログラム - Google Patents
放射線画像補正方法および装置ならびにプログラム Download PDFInfo
- Publication number
- JP2010005373A JP2010005373A JP2009025677A JP2009025677A JP2010005373A JP 2010005373 A JP2010005373 A JP 2010005373A JP 2009025677 A JP2009025677 A JP 2009025677A JP 2009025677 A JP2009025677 A JP 2009025677A JP 2010005373 A JP2010005373 A JP 2010005373A
- Authority
- JP
- Japan
- Prior art keywords
- image
- radiation
- radiographic image
- correction
- region
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims description 29
- 238000003702 image correction Methods 0.000 title claims description 24
- 230000005855 radiation Effects 0.000 claims abstract description 137
- 238000012937 correction Methods 0.000 claims description 119
- 210000000038 chest Anatomy 0.000 claims description 36
- 230000001678 irradiating effect Effects 0.000 claims description 21
- 238000012886 linear function Methods 0.000 claims description 13
- 239000000284 extract Substances 0.000 claims description 4
- 230000000306 recurrent effect Effects 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 12
- 210000004072 lung Anatomy 0.000 description 11
- 210000001991 scapula Anatomy 0.000 description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 238000002601 radiography Methods 0.000 description 4
- 210000000746 body region Anatomy 0.000 description 2
- 230000037237 body shape Effects 0.000 description 2
- 239000000428 dust Substances 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 210000003109 clavicle Anatomy 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010422 painting Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Classifications
-
- G06T5/90—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/28—Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
Abstract
【課題】放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像における再現性のない濃度ムラを補正する。
【解決手段】放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正する。
【選択図】図2
【解決手段】放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正する。
【選択図】図2
Description
本発明は、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する放射線画像補正方法および装置ならびにプログラムに関するものである。
医療診断を目的とする放射線撮影において、放射線照射手段により被写体に向けて照射された放射線を検出して電気信号に変換することにより、被写体の放射線画像を取得する放射線画像検出器が知られている。放射線画像検出器としては、放射線を照射するとこの放射線エネルギーの一部が蓄積され、その後可視光等の励起光を照射すると蓄積されたエネルギーに応じて輝尽発光を示す蓄積性蛍光体(輝尽性蛍光体)を利用して、被写体の放射線像を一旦蓄積性蛍光体シートに記録し、この蓄積性蛍光体シートをレーザ光などの励起光で走査して輝尽発光光を生ぜしめ、その輝尽発光光を検出して上記放射線像を表す画像信号を得るCR(Computed Radiography)方式の放射線画像検出器と、放射線を照射するとこの放射線エネルギーに応じた電荷を発生する固体センサーを利用して、被写体の放射線像を電荷に変換して蓄積し、その蓄積した電荷を薄膜トランジスタまたは光の照射により電荷を発生する半導体材料を利用して読み取る方式の放射線画像検出器等がある。
ところで、放射線撮影においては、放射線照射手段における放射線照射の強度ムラ、放射線画像検出器の受光面における感度ムラ等による影響が被写体の正確な放射線画像情報を得る上で障害となる。
この問題を解決するため、特許文献1においては、事前に被写体を透過させることなく撮影して得られた放射線画像から補正データを作成し、この補正データを用いて、被写体を透過させて撮影して得られた放射線画像の画像信号を補正する方法が提案されている。
しかしながら、この特許文献1において提案されている方法では、毎回の撮影において同じように発生するムラを補正することは可能であるが、撮影の度に不規則的に発生するムラ(再現性のないムラ)を補正することはできない。
これに対して、特許文献2には、放射線画像の各画素の信号値を水平方向(あるいは垂直方向)に加算することによりプロファイルを作成し、そのプロファイルが急激に変化している箇所を抽出し、その箇所のプロファイルがなだらかに変化するように放射線画像の画像信号を補正することにより、光電変換するときに発生した振動、蓄積性蛍光体シートから画像信号を読取る際に用いられる光学系に付着したゴミ、そのゴミの移動等に起因する、再現性のないムラを補正する方法が提案されている。
また、非特許文献1では、放射線画像中の、放射線が直接入射した領域である素抜け領域の画素値を用いて放射線画像の濃度ムラを近似した補正データを作成し、その放射線画像から補正データを減算することにより、放射線画像における再現性のないムラを補正する方法が提案されている。
Bernardino Arcay Varela et al., "Preprocessing of Radiological Images: Comparison of the Application of Polynomic Algorithms and Artificial Neural Networks to the Elimination of Variations in Background Luminosity", Lect. Notes Comput. Sci.1607, 452-459, 1999
ところで、回診車などで回診撮影を行う場合には、X線撮影専用の撮影室と異なり、ベッドに横になった被検者の背中側にイメージングプレート(蓄積性蛍光体シート)、フラットパネルディテクター等の放射線画像検出器を設置して撮影を行うため、ベッドが柔らかいなどの撮影環境の影響を受けやすく、放射線照射手段による放射線の照射方向に対して放射線画像検出器の位置や配置方向が撮影の度に不規則的にずれてしまい、撮影して得られた放射線画像の全体にかけてなだらかに変化する、再現性のない濃度ムラが発生する場合がある。
しかしながら、上記特許文献2の方法では、作成したプロファイルが急激に変化している箇所を抽出し、その箇所のプロファイルがなだらかに変化するように画像信号を補正するので、放射線画像の全体にかけてなだらかに変化する、上記再現性のない濃度ムラは補正することができない。
また、上記非特許文献1の方法では、素抜け領域の画素値を用いて放射線画像全体の濃度ムラを近似した補正データを作成するので、素抜け領域が、放射線画像中に少なくともその補正データが求められる程度に広がっていることが必要となるが、被写体である患者の姿勢や撮影部位などによっては、放射線画像の中に素抜け領域が存在しなかったり、存在しても小さすぎたりする場合があり、このような場合には、放射線画像全体における再現性のない濃度ムラを補正することができない。
本発明は、上記事情に鑑み、放射線画像中における素抜け領域の存在の有無に拘わらず、濃度ムラ、特に再現性のない濃度ムラを補正することが可能な画像補正方法および装置ならびにプログラムを提供することを目的とするものである。
本発明の放射線画像補正方法は、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する放射線画像補正方法において、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度(画素値)になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正することを特徴とするものである。
本発明の放射線画像補正装置は、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する放射線画像補正装置において、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得する領域取得手段と、取得した領域の画素値を用いて放射線画像の画素値を補正する補正手段とを備えたことを特徴とするものである。
上記装置においては、放射線画像が人体の胸部正面放射線画像であり、領域取得手段が、胸部正面放射線画像中の、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域を略等濃度になる領域として取得するものであってもよい。
また、上記装置は、略等濃度になる領域を特定する情報を入力するための特定情報入力手段をさらに備え、領域取得手段が、特定情報入力手段により入力された情報に基づいて略等濃度になる領域を取得するものであってもよい。
また、補正手段は、略等濃度になる領域の画素値を用いて放射線画像の濃度ムラを近似した線形関数を求め、求められた線形関数を用いて放射線画像の画素値を補正するものであってもよい。
補正手段は、肋骨の重なりが撮影された領域の長手方向の少なくとも一方の長さが所定の閾値以下であるときは、肋骨の重なりが撮影された領域の画素値を用いて放射線画像の濃度ムラを近似した線形関数を求め、両方の長さが所定の閾値より大きいときは、肋骨の重なりが撮影された領域の画素値を用いて放射線画像の濃度ムラを近似した非線形関数を求め、求められた関数を用いて放射線画像の画素値を補正するものであってもよい。
また、上記装置は、濃度ムラを近似した関数のパラメータを表示する表示手段と、表示手段により表示されたパラメータの修正値を入力するための修正値入力手段とをさらに備え、補正手段が、修正値入力手段により入力された修正値を用いてパラメータを修正し、修正後のパラメータにより規定される関数を用いて放射線画像の画素値を補正するものであってもよい。
また、補正手段による補正に用いられた関数を、補正が施された放射線画像に関連付けて保存する保存手段をさらに備えたものであってもよい。
また、放射線画像が、人体の胸部正面放射線画像であり、補正手段が、胸部正面放射線画像から人体の体軸を抽出し、略等濃度になる領域の画素値を用いて、胸部正面放射線画像におけるその抽出された体軸に直交する方向に存在する濃度ムラを近似した関数を求め、求められた関数を用いて胸部正面放射線画像の画素値を補正するものであってもよい。
また、胸部正面放射線画像における、体軸方向に存在する濃度ムラおよび体軸に直交する方向に存在する濃度ムラと、体軸に直交する方向に存在する濃度ムラのうち、いずれの濃度ムラを補正手段による補正の対象とするかを特定する情報を入力するための対象特定情報入力手段をさらに備え、補正手段が、対象特定情報入力手段により入力された情報に基づいて、体軸方向に存在する濃度ムラおよび体軸に直交する方向に存在する濃度ムラを近似した関数と、体軸に直交する方向に存在する濃度ムラを近似した関数のうち、いずれかの関数を求め、求められた関数を用いて胸部正面放射線画像の画素値を補正するものであってもよい。
本発明のプログラムは、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正するための放射線画像補正プログラムであって、コンピュータに、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正することを実行させるためのものである。
略等濃度になる領域とは、その領域において検出された放射線量の最大値と最小値の比が2倍程度かそれ以下となる領域であって、放射線画像の濃度(画素値)は、この検出された放射線量の比が差として表現された空間へ対数変換されたものである。
なお、略等濃度になる領域としては、放射線画像の全体の濃度ムラを近似する関数が求められる程度に広がっている領域であることが望ましい。
また、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であるとは、もし被写体がなければ、放射線画像検出器の全領域に同一の照射量で放射線が照射されることを意味する。
また、放射線画像検出器の全領域とは、実質的な画像信号の検出領域を意味する。
また、放射線画像検出器とは、イメージングプレートや、フラットパネルディテクター等を意味する。
本発明の放射線画像補正方法および装置ならびにプログラムによれば、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する際、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正することにより、放射線画像中における素抜け領域の存在の有無に拘わらず、その放射線の照射が均一であれば略等濃度になるはずの領域の画素値に表れる照射ムラの影響を基に、放射線画像全体における再現性のない濃度ムラさえ補正することができる。
上記方法および装置ならびにプログラムにおいては、人体の胸部正面放射線画像中の、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域を略等濃度になる領域として取得することができる。この領域は、放射線撮影において、放射線照射手段により被写体に向けて照射される放射線の強度を放射線画像検出器の全領域に対して均一にした場合、略等濃度になる放射線画像中の領域である。
以下、図面を参照して本発明の放射線画像補正装置の実施の形態について説明する。なお、図1に示す本発明の一実施形態である放射線画像補正装置1は、補助記憶装置に読み込まれた放射線画像補正プログラムをコンピュータ(たとえばパーソナルコンピュータ等)上で実行することにより実現される。このとき、この放射線画像補正プログラムは、CD‐ROM等の情報記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。
図1に示すように、放射線画像補正装置1は、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正するものであって、領域取得部10、補正部20、入力部30、表示部40、保存部50などを備えている。
領域取得部10は、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得するものであり、たとえば図2に示すように、放射線画像が人体の胸部正面放射線画像である場合、胸部正面放射線画像中の、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域A(図中白抜きでの実線で囲まれた領域)を略等濃度になる領域として取得する。
領域取得部10は、このように略等濃度になる領域を自動検出して取得するものであってもよいし、ユーザーによる略等濃度になる領域を特定する情報を入力部30から受け付け、その情報に基づいて略等濃度になる領域を取得するものであってもよい。ここで、略等濃度になる領域を特定する情報は、その情報から直接的に等濃度になる領域が定まるような情報であり、たとえば、等濃度になる領域として特定したい領域が一意的に決まるよう、その領域を囲むように指定された幾つかの点の位置情報、塗りつぶしで等濃度になる領域として指定された領域の情報などを含む。
略等濃度になる領域として肋骨の重なりが撮影された領域を自動検出する場合には、たとえば特開2003‐6661号公報に記載されている技術を用いて、右肺外側輪郭、右肺内側輪郭、右肺下側輪郭、左肺外側輪郭、左肺内側輪郭、左肺下側輪郭で囲まれた領域を決定し、決定した領域の左右の輪郭から外側に、肋骨の重なりが撮影された領域が一般的に有する程度の幅の領域を抽出することによりその肋骨の重なりが撮影された領域を取得する。また、ここで説明した取得方法に限らず、他の技術を用いて取得するようにしてもよい。
補正部20は、領域取得部10において取得された略等濃度になる領域の画素値を用いて放射線画像の濃度ムラを近似した補正データを作成する補正データ作成部21と、作成した補正データを用いて放射線画像の画素値を補正する補正処理部22と、胸部正面放射線画像から人体の体軸を抽出する体軸抽出部23を備える。
補正データ作成部21は、まず、たとえば下記の式(1)のように、略等濃度になる領域の画像信号を近似した線形モデルZ´(x, y)を作成する。ここで、定数a、b、cは、例えば最小二乗法により最も確からしい線形モデルZ´(x, y)を構成するように選ばれたパラメータである。具体的には、略等濃度になる領域内の各画素の画素値A(x, y) とZ´(x, y)との差の2乗和が最も小さくなるようなパラメータを探索して求める。
次に、その求められたパラメータa,bおよび放射線画像の中心の座標(xc, yc)を用いて、下記の式(2)のように、放射線画像の水平方向xおよび垂直方向yにおける濃度ムラを近似した線形モデルZ(x, y)を補正データとして作成し、作成された補正データを補正処理部22に出力する。
図2に示す画像Zは、胸部正面放射線画像である放射線画像P中の、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域A内の各画素の画素値を用いて、放射線画像Pの濃度ムラを線形関数に近似して取得した補正データの一例を示すものである。
ここでは、補正データ作成部21が、放射線画像の濃度ムラを線形関数に近似して補正データを作成するものである場合について説明したが、肋骨の重なりが撮影された領域Aの長手方向の長さの少なくとも一方が所定の閾値以下であるときは、放射線画像の濃度ムラを線形関数に近似して補正データを作成し、その長さがいずれも所定の閾値より大きいときは、放射線画像の濃度ムラを非線形関数に近似して補正データを作成するものであってもよい。例えば、被写体の体軸がその放射線画像の中心を通って該放射線画像の垂直方向yに伸びている場合に、肋骨の重なりが撮影された領域Aの長手方向の長さの少なくとも一方がその放射線画像の垂直方向yの長さの7割以下のときは、放射線画像の濃度ムラを線形関数に近似して補正データを作成し、いずれの長さも7割より大きいときは、放射線画像の濃度ムラを非線形モデルに近似して補正データを作成することができる。
以下、放射線画像の濃度ムラを近似した非線形モデルの生成方法について説明する。まず、たとえば下記の式(3)のように、略等濃度になる領域の画像信号を近似した非線形モデルZ´(x, y)を作成する。ここで、定数A、B、C、Dは、上記線形モデルの場合と同様に、最小二乗法等の手法により最も確からしい非線形モデルZ´(x, y)を構成するように選ばれたパラメータである。
次に、その求められたパラメータA、B、Cおよび放射線画像の中心の座標(xc, yc)を用いて、下記の式(4)のように、放射線画像の水平方向xおよび垂直方向yにおける濃度ムラを近似した非線形モデルZ(x, y)を補正データとして作成し、作成された補正データを補正処理部22に出力する。
なお、肋骨の重なりが撮影された領域Aは、被写体の体型によっては体軸方向に濃度変化を伴う場合があるので、その領域Aの画素値を用いて補正データを作成する場合には、体軸に直交する方向に存在する濃度ムラのみを補正の対象とし、その体型に起因して現れる濃度ムラを補正の対象から除外することもできる。
以下に、図3を参照して、体軸に直交する方向に存在する濃度ムラのみを近似した補正データを取得する方法について説明する。
まず、体軸抽出部23が、胸部正面放射線画像から人体の体軸26を抽出する。体軸抽出部23は、例えば特願2008-037145号明細書に記載されているように、ガボールフィルタを用いて胸部画像の各画素において、椎体の太さに対応するエッジ方向値および椎体の太さに対応するエッジ強度値を含むエッジ成分値を抽出し、次に、胸部画像中の少なくとも鎖骨と肋骨が重なる領域である左右の測端部を含まない領域として設定された関心領域中の、各画素における最も高いエッジ強度値に対応するエッジ方向値を、その画素における最も高いエッジ強度値により重み付け平均して得られた方向を椎体方向と推定する。次に、推定された椎体方向に対して略垂直に胸部画像を走査し、所定の画素値以下である各画素を椎体領域として抽出し、抽出された椎体領域の正中線を人体の正中線として検出する。
次に、補正データ作成部21が、画像のy軸を基準とした被写体の体軸26の傾斜角φを取得し、下記の式(5)のように、略等濃度になる領域の画像信号を近似した線形モデルZ´(x, y)を作成する。ここで、定数a、cは、上述した他の補正データ作成方法と同様に、最小二乗法等の手法により最も確からしい線形モデルZ´(x, y)を構成するように選ばれたパラメータである。
次に、その求められたパラメータaを用いて、下記の式(6)のように、体軸に直交する方向に存在する濃度ムラを近似した線形モデルZ(x, y)を補正データとして作成する。
なお、肋骨の重なりが撮影された領域Aの画素値を用いて胸部正面放射線画像の画素値を補正する場合に、体軸方向に存在する濃度ムラおよび体軸に直交する方向に存在する濃度ムラと、体軸に直交する方向に存在する濃度ムラのうち、いずれの濃度ムラを補正手段による補正の対象とするかは、領域Aの体軸方向における濃度勾配の程度に基づいて自動決定されるようにしてもよいし、ユーザーが入力部30を介して入力した、そのいずれの濃度ムラを補正の対象とするかを特定する情報に基づいて決定されるようにしてもよい。ここで、そのいずれの濃度ムラを補正の対象とするかを特定する情報は、たとえば表示部40の表示画面上に、補正の対象としてそのいずれの濃度ムラを選択して設定できるように表示された設定画面において、ユーザーがそのいずれの濃度ムラを選択して設定した設定情報である。
なお、以上のように濃度ムラを近似した関数(補正データ)のパラメータは表示部40により図示していないモニターの画面上に表示され、ユーザーが画面上に表示されたパラメータに対して、入力部30により修正値を入力すると、補正データ作成部21が、ユーザーによる修正値を入力部30から受け付け、その入力された修正値を用いてパラメータを修正する。そして、修正後のパラメータにより規定される関数を補正データとして補正処理部22に出力することもできる。
図4は、表示部40により濃度ムラを近似した線形関数のパラメータa、bをモニターの画面上に表示する態様の一例を示すものであり、グラフ41上でパラメータa、bの値により規定される点Qをユーザーがカーソルで移動させたり、新たな位置を指定する等の操作により、パラメータa、bを所望の値に修正することができる。
補正処理部22は、補正データ作成部21において作成された補正データを用いて放射線画像の画素値を補正するものであり、図2に示すように、放射線画像Pの各画素P(x, y)から、補正データ作成部21により作成された補正データZ(x, y)を減算することにより、濃度ムラを補正した各画素値S(x, y)を算出してムラ補正済み画像Sを作成する。
以上のように作成された補正済み画像Sは、表示部40によりモニターの画面上に表示される。そのとき、補正処理部22による補正に用いられた補正データZ(x, y)のパラメータや、補正データZ(x, y)の作成に用いられた略等濃度になる領域の輪郭等の情報も画面上に表示させることにより、どのような補正処理が行われたかをユーザーに分かりやすく提示することができる。
なお、ユーザーが画面上に表示されたパラメータに対して、入力部30により修正値を入力すると、補正データ作成部21が、ユーザーによる修正値を入力部30から受け付け、その入力された修正値を用いてパラメータを修正し、補正処理部22が、修正後のパラメータにより規定される関数を用いて放射線画像の画素値を改めて補正し、ムラ補正済み画像Sを作成するようにしてもよい。
また、記録部50(保存手段)は、上記補正が施された放射線画像Sと、その補正に用いられた補正データZ(x, y)を互いに関連付けて図示しない記録媒体に保存するものである。
上記構成により、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する際、まず、領域取得部10が、放射線画像Pの被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得する。次に、補正データ作成部21が、領域取得部10において取得した略等濃度になる領域の画素値を用いて放射線画像の濃度ムラを近似した補正データZ(x, y)を作成し、補正処理部22が、補正データ作成部21において作成された補正データを用いて放射線画像Pの画素値を補正し、補正済み画像Sを作成する。次に、表示部40が、補正処理部22において作成された補正済み画像S、その補正に用いられた補正データZ(x, y)、補正データZ(x, y)の作成に用いられた略等濃度になる領域の輪郭等をモニターの画面上に表示させる。
ユーザーが画面上に表示されたパラメータに対して、入力部30により修正値を入力すると、補正データ作成部21が、ユーザーによる修正値を入力部30から受け付け、その入力された修正値を用いてパラメータを修正し、補正処理部22が、修正後のパラメータにより規定される関数を用いて放射線画像の画素値を改めて補正し、ムラ補正済み画像Sを作成する。そして、記録部50が、最終的に作成されたムラ補正済み画像Sおよびその補正に用いられた補正データZ(x, y)を互いに関連付けて記録媒体に保存する。
なお、ここでは、補正データ作成部21において作成した補正データZ(x, y)をそのまま用いて補正を行うことにより補正済み画像Sを作成した後、その結果をユーザーに提示して、必要に応じて補正データに修正を加え、修正後の補正データを用いて改めて補正を行う場合について説明したが、最初の補正を行なう前に、補正データ作成部21において作成した補正データZ(x, y)や、その補正データZ(x, y)の作成に用いられた略等濃度になる領域の輪郭等を一旦モニターの画面上に表示させ、必要に応じて補正データにユーザーによる修正を受け付け、修正後の補正データを用いて補正を行うようにしてもよい。
上記実施の形態によれば、放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した被写体の放射線画像の画素値を補正する際、放射線画像の被写体が撮影された領域中の、放射線照射手段による放射線の照射が放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、取得した領域の画素値を用いて放射線画像の画素値を補正することにより、放射線画像中における素抜け領域の存在の有無に拘わらず、その放射線画像中の放射線の照射が均一であれば略等濃度になるはずの領域の画素値に表れる照射ムラの影響を基に、放射線画像全体における再現性のない濃度ムラさえを補正することができる。
また、略等濃度になる領域として取得された肋骨の重なりが撮影された領域の長手方向の長さが所定の閾値より大きいときは、非線形の関数を用いることにより放射線画像の濃度ムラをより正確に近似することができる。なお、少なくとも一方の長さが所定の閾値以下であるときには、非線形関数で濃度ムラを近似すると、不適切な補正データが生成される虞があるので、線形の関数を用いて補正データを生成することがより好ましい。
なお、上記実施の形態では、放射線画像が人体の胸部正面放射線画像であるときに、領域取得部10が、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域を略等濃度になる領域として取得する場合を例示して説明しているが、それに限らず、たとえば、図5に示すような胸郭の左右の上部外側にある肩甲骨の領域B(図中白抜きでの実線で囲まれた領域)を略等濃度になる領域として取得し、その領域Bの画素値を用いて放射線画像の体軸に直交する方向に存在する濃度ムラを近似した補正データを作成し、放射線画像から、その作成した補正データを減算することにより、放射線画像全体における体軸に直交する方向の、再現性のない濃度ムラを補正するようにしてもよい。この領域Bは、体軸に対して左右対称で、かつ胸郭の外側(左右差があり得る肺野の影響を受けない位置)に存在するので略等濃度とみなせるとともに、胸部正面放射線画像に高確率で写り込んでいるので実用的である。
肩甲骨の領域Bの取得は、領域取得部10に肩甲骨の領域Bを自動決定するためのプログラムを設け、そのプログラムにより、例えば特開2003‐6661号公報に記載されている技術を用いて左右両肺野の外側の輪郭線を抽出し、両肺野の上部外側に図5に示すような三角形の領域Bを設定することをコンピュータに実行させて行う。ここで、三角形の領域Bの設定は、まず、体軸に直交する方向に伸びる上辺の高さが肺野の上端の高さに一致するように設定するとともに、内側の辺は肺上部の胸郭(肋骨重なりによる輪郭線)に接するように設定する。次いで、上辺の長さを、所定の値(例えば、体軸から肺野の外側の輪郭線までの最大距離の1/2)に設定し、体軸に平行になるように外側の辺を設定する。
なお、領域取得部10は、このように略等濃度になる領域である肩甲骨の領域Bを自動で取得するものであってもよいし、ユーザーが、肩甲骨の領域Bが一意的に決まるよう、その領域を囲むような幾つかの点の位置情報を入力し、あるいは塗りつぶし等で肩甲骨の領域Bの具体的な位置を特定する情報を入力することにより、その情報に基づいて肩甲骨の領域Bを取得するものであってもよい。
1 放射線画像補正装置
10 領域取得部
20 補正部
30 入力部
40 表示部
50 保存部
A、B 略等濃度になる領域
Z 補正データ
P 放射線画像
10 領域取得部
20 補正部
30 入力部
40 表示部
50 保存部
A、B 略等濃度になる領域
Z 補正データ
P 放射線画像
Claims (11)
- 放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した前記被写体の放射線画像の画素値を補正する放射線画像補正方法において、
前記放射線画像の前記被写体が撮影された領域中の、前記放射線照射手段による放射線の照射が前記放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、
該取得した領域の画素値を用いて前記放射線画像の画素値を補正する
ことを特徴とする放射線画像補正方法。 - 放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した前記被写体の放射線画像の画素値を補正する放射線画像補正装置において、
前記放射線画像の前記被写体が撮影された領域中の、前記放射線照射手段による放射線の照射が前記放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得する領域取得手段と、
該取得した領域の画素値を用いて前記放射線画像の画素値を補正する補正手段と
を備えたことを特徴とする放射線画像補正装置。 - 前記放射線画像が、人体の胸部正面放射線画像であり、
前記領域取得手段が、前記胸部正面放射線画像中の、胸郭の左右の外輪郭部に存在する肋骨の重なりが撮影された領域を前記略等濃度になる領域として取得するものであることを特徴とする請求項2記載の放射線画像補正装置。 - 前記略等濃度になる領域を特定する情報を入力するための領域特定情報入力手段をさらに備え、
前記領域取得手段が、前記領域特定情報入力手段により入力された前記情報に基づいて前記略等濃度になる領域を取得するものであることを特徴とする請求項2または3記載の放射線画像補正装置。 - 前記補正手段が、前記略等濃度になる領域の画素値を用いて前記放射線画像の濃度ムラを近似した線形関数を求め、該求められた線形関数を用いて前記放射線画像の画素値を補正するものであることを特徴とする請求項2から4のいずれか1項記載の放射線画像補正装置。
- 前記補正手段が、前記肋骨の重なりが撮影された領域の長手方向の長さの少なくとも一方が所定の閾値以下であるときは、前記肋骨の重なりが撮影された領域の画素値を用いて前記放射線画像の濃度ムラを近似した線形関数を求め、前記肋骨の重なりが撮影された領域の長手方向の長さがいずれも前記所定の閾値より大きいときは、前記肋骨の重なりが撮影された領域の画素値を用いて前記放射線画像の濃度ムラを近似した非線形関数を求め、該求められた関数を用いて前記放射線画像の画素値を補正するものであることを特徴とする請求項3記載の放射線画像補正装置。
- 前記濃度ムラを近似した関数のパラメータを表示する表示手段と、
該表示手段により表示されたパラメータの修正値を入力するための修正値入力手段とをさらに備え、
前記補正手段が、前記修正値入力手段により入力された修正値を用いて前記パラメータを修正し、該修正後のパラメータにより規定される関数を用いて前記放射線画像の画素値を補正するものであることを特徴とする請求項5または6記載の放射線画像補正装置。 - 前記補正手段による補正に用いられた関数を、前記補正が施された放射線画像に関連付けて保存する保存手段をさらに備えたことを特徴とする請求項5から7のいずれか1項記載の放射線画像補正装置。
- 前記放射線画像が、人体の胸部正面放射線画像であり、
前記補正手段が、前記胸部正面放射線画像から前記人体の体軸を抽出し、前記略等濃度になる領域の画素値を用いて、前記胸部正面放射線画像における前記抽出された体軸に直交する方向に存在する濃度ムラを近似した関数を求め、該求められた関数を用いて前記胸部正面放射線画像の画素値を補正するものであることを特徴とする請求項2記載の放射線画像補正装置。 - 前記胸部正面放射線画像における、前記体軸方向に存在する濃度ムラおよび前記体軸に直交する方向に存在する濃度ムラと、前記体軸に直交する方向に存在する濃度ムラのうち、いずれの濃度ムラを前記補正手段による補正の対象とするかを特定する情報を入力するための対象特定情報入力手段をさらに備え、
前記補正手段が、前記対象特定情報入力手段により入力された前記情報に基づいて、前記体軸方向に存在する濃度ムラおよび前記体軸に直交する方向に存在する濃度ムラを近似した関数と、前記体軸に直交する方向に存在する濃度ムラを近似した関数のうち、いずれかの関数を求め、該求められた関数を用いて前記胸部正面放射線画像の画素値を補正するものであることを特徴とする請求項9記載の放射線画像補正装置。 - 放射線照射手段により被写体に向けて照射された放射線を放射線画像検出器で検出することにより取得した前記被写体の放射線画像の画素値を補正するための放射線画像補正プログラムであって、
コンピュータに、
前記放射線画像の前記被写体が撮影された領域中の、前記放射線照射手段による放射線の照射が前記放射線画像検出器の全領域に対して均一であれば略等濃度になる領域を取得し、
該取得した領域の画素値を用いて前記放射線画像の画素値を補正することを実行させるための放射線画像補正プログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009025677A JP2010005373A (ja) | 2008-02-14 | 2009-02-06 | 放射線画像補正方法および装置ならびにプログラム |
US12/371,279 US20090208087A1 (en) | 2008-02-14 | 2009-02-13 | Radiographic image correction method, apparatus and recording-medium stored therein program |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008032784 | 2008-02-14 | ||
JP2008140959 | 2008-05-29 | ||
JP2009025677A JP2010005373A (ja) | 2008-02-14 | 2009-02-06 | 放射線画像補正方法および装置ならびにプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010005373A true JP2010005373A (ja) | 2010-01-14 |
Family
ID=40955169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009025677A Abandoned JP2010005373A (ja) | 2008-02-14 | 2009-02-06 | 放射線画像補正方法および装置ならびにプログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US20090208087A1 (ja) |
JP (1) | JP2010005373A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012102052A1 (ja) * | 2011-01-27 | 2012-08-02 | 富士フイルム株式会社 | 放射線画像撮影装置および方法 |
JP2018143699A (ja) * | 2017-03-09 | 2018-09-20 | 富士フイルム株式会社 | 医用画像撮影制御装置、方法およびプログラム |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009285145A (ja) * | 2008-05-29 | 2009-12-10 | Fujifilm Corp | 放射線画像補正装置および方法ならびにプログラム |
JP2018157884A (ja) * | 2017-03-22 | 2018-10-11 | コニカミノルタ株式会社 | X線動画像処理装置 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01125675A (ja) * | 1987-08-03 | 1989-05-18 | Toshiba Corp | ディジタルx線画像における自動化分析のための方法および装置 |
JPH0449950A (ja) * | 1990-06-15 | 1992-02-19 | Fuji Photo Film Co Ltd | 画像濃度調整装置 |
JPH1063838A (ja) * | 1996-08-19 | 1998-03-06 | Fuji Photo Film Co Ltd | 画像処理方法および装置 |
JPH10105701A (ja) * | 1996-09-25 | 1998-04-24 | Fuji Photo Film Co Ltd | 放射線画像強調処理方法および装置 |
JP2002269537A (ja) * | 2001-03-14 | 2002-09-20 | Canon Inc | 画像処理装置、方法およびプログラム |
JP2002366961A (ja) * | 2001-04-03 | 2002-12-20 | Fuji Photo Film Co Ltd | 画像輪郭検出方法および装置 |
JP2003000579A (ja) * | 2001-05-10 | 2003-01-07 | Agfa Gevaert Nv | 放射線写真の不均一性の遡及的修正方法 |
JP2005211439A (ja) * | 2004-01-30 | 2005-08-11 | Fuji Photo Film Co Ltd | 異常陰影表示装置およびそのプログラム |
JP2006000223A (ja) * | 2004-06-15 | 2006-01-05 | Canon Inc | X線ct装置 |
JP2006130212A (ja) * | 2004-11-09 | 2006-05-25 | Fuji Photo Film Co Ltd | 異常陰影候補検出方法、および装置、ならびにプログラム |
JP2006149960A (ja) * | 2004-12-01 | 2006-06-15 | Matsushita Electric Ind Co Ltd | X線撮影装置 |
JP2007068715A (ja) * | 2005-09-06 | 2007-03-22 | Canon Inc | 医用画像処理装置及び方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60234643A (ja) * | 1984-05-09 | 1985-11-21 | 富士写真フイルム株式会社 | 放射線画像情報再生方法及び装置 |
US5471987A (en) * | 1993-03-30 | 1995-12-05 | Konica Corporation | Method of compressing a dynamic range for a radiation image |
US6314198B1 (en) * | 1996-09-25 | 2001-11-06 | Canon Kabushiki Kaisha | Radiographic, digital image processing system |
JP2003190125A (ja) * | 2001-12-27 | 2003-07-08 | Konica Corp | 画像処理装置、画像処理方法、プログラム、及び記憶媒体 |
JP4908289B2 (ja) * | 2007-03-30 | 2012-04-04 | 富士フイルム株式会社 | 画像処理装置、方法及びプログラム |
US20090060372A1 (en) * | 2007-08-27 | 2009-03-05 | Riverain Medical Group, Llc | Object removal from images |
-
2009
- 2009-02-06 JP JP2009025677A patent/JP2010005373A/ja not_active Abandoned
- 2009-02-13 US US12/371,279 patent/US20090208087A1/en not_active Abandoned
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01125675A (ja) * | 1987-08-03 | 1989-05-18 | Toshiba Corp | ディジタルx線画像における自動化分析のための方法および装置 |
JPH0449950A (ja) * | 1990-06-15 | 1992-02-19 | Fuji Photo Film Co Ltd | 画像濃度調整装置 |
JPH1063838A (ja) * | 1996-08-19 | 1998-03-06 | Fuji Photo Film Co Ltd | 画像処理方法および装置 |
JPH10105701A (ja) * | 1996-09-25 | 1998-04-24 | Fuji Photo Film Co Ltd | 放射線画像強調処理方法および装置 |
JP2002269537A (ja) * | 2001-03-14 | 2002-09-20 | Canon Inc | 画像処理装置、方法およびプログラム |
JP2002366961A (ja) * | 2001-04-03 | 2002-12-20 | Fuji Photo Film Co Ltd | 画像輪郭検出方法および装置 |
JP2003000579A (ja) * | 2001-05-10 | 2003-01-07 | Agfa Gevaert Nv | 放射線写真の不均一性の遡及的修正方法 |
JP2005211439A (ja) * | 2004-01-30 | 2005-08-11 | Fuji Photo Film Co Ltd | 異常陰影表示装置およびそのプログラム |
JP2006000223A (ja) * | 2004-06-15 | 2006-01-05 | Canon Inc | X線ct装置 |
JP2006130212A (ja) * | 2004-11-09 | 2006-05-25 | Fuji Photo Film Co Ltd | 異常陰影候補検出方法、および装置、ならびにプログラム |
JP2006149960A (ja) * | 2004-12-01 | 2006-06-15 | Matsushita Electric Ind Co Ltd | X線撮影装置 |
JP2007068715A (ja) * | 2005-09-06 | 2007-03-22 | Canon Inc | 医用画像処理装置及び方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012102052A1 (ja) * | 2011-01-27 | 2012-08-02 | 富士フイルム株式会社 | 放射線画像撮影装置および方法 |
JP2018143699A (ja) * | 2017-03-09 | 2018-09-20 | 富士フイルム株式会社 | 医用画像撮影制御装置、方法およびプログラム |
US11185312B2 (en) | 2017-03-09 | 2021-11-30 | Fujifilm Corporation | Medical image capturing control device, method, and program |
Also Published As
Publication number | Publication date |
---|---|
US20090208087A1 (en) | 2009-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1788525B1 (fr) | Procédé d'imagerie radiographique pour la reconstruction tridimensionnelle, dispositif et programme d'ordinateur pour mettre en oeuvre ce procédé | |
JP5124363B2 (ja) | エネルギーサブトラクション処理装置、方法、およびプログラム | |
JP5224898B2 (ja) | 放射線画像撮影装置及びその駆動方法 | |
JP2010263961A (ja) | X線画像撮影装置およびx線画像撮影装置の制御方法 | |
JP6363573B2 (ja) | 線源画像面間距離取得装置、方法およびプログラム、並びに放射線画像処理装置、方法およびプログラム | |
WO2015133123A1 (ja) | 放射線画像処理装置および方法並びにプログラム | |
US11860111B2 (en) | Image reconstruction method for X-ray measuring device, structure manufacturing method, image reconstruction program for X-ray measuring device, and X-ray measuring device | |
JP3888046B2 (ja) | 放射線画像処理方法および放射線画像処理装置 | |
EP1903499A2 (en) | Radiological image capturing system and radiological image capturing method | |
JP2002085392A (ja) | 放射線画像処理方法および放射線画像処理装置 | |
JP4159701B2 (ja) | デジタル放射線画像の評価方法および評価装置 | |
WO2014050045A1 (ja) | 体動検出装置および方法 | |
JP4129598B2 (ja) | 画像処理装置及び画像処理方法 | |
JP2010005373A (ja) | 放射線画像補正方法および装置ならびにプログラム | |
US7260254B2 (en) | Comparing images | |
JP6105903B2 (ja) | 画像処理装置、画像処理方法、放射線撮影システム及びプログラム | |
JP2001238868A (ja) | 画像処理方法及び画像処理装置 | |
JP2005109867A (ja) | 画像処理方法および画像処理装置ならびに画像処理プログラム | |
JP2001120524A (ja) | 放射線画像処理装置 | |
JP3731400B2 (ja) | 画像処理方法および画像処理装置 | |
JP2009285145A (ja) | 放射線画像補正装置および方法ならびにプログラム | |
JP3661797B2 (ja) | 放射線画像の位置合せ用対応点決定方法 | |
JP2008131981A (ja) | 画質評価演算方法および装置ならびにプログラム | |
JP2001245140A (ja) | 画像処理方法および画像処理装置 | |
KR101709406B1 (ko) | Cr 시스템에서 방사선량을 산출하는 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110708 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130221 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130312 |
|
A762 | Written abandonment of application |
Free format text: JAPANESE INTERMEDIATE CODE: A762 Effective date: 20130508 |