JP4159701B2 - デジタル放射線画像の評価方法および評価装置 - Google Patents
デジタル放射線画像の評価方法および評価装置 Download PDFInfo
- Publication number
- JP4159701B2 JP4159701B2 JP11223599A JP11223599A JP4159701B2 JP 4159701 B2 JP4159701 B2 JP 4159701B2 JP 11223599 A JP11223599 A JP 11223599A JP 11223599 A JP11223599 A JP 11223599A JP 4159701 B2 JP4159701 B2 JP 4159701B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- value
- digital
- evaluation
- radiation
- 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
- 238000011156 evaluation Methods 0.000 title claims description 103
- 238000012937 correction Methods 0.000 claims description 71
- 230000005855 radiation Effects 0.000 claims description 47
- 238000000034 method Methods 0.000 claims description 25
- 238000013441 quality evaluation Methods 0.000 claims description 19
- 238000001228 spectrum Methods 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 230000004936 stimulating effect Effects 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 34
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 5
- 238000002601 radiography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 229910001632 barium fluoride Inorganic materials 0.000 description 1
- -1 barium fluoride halide Chemical class 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or 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/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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
- Radiography Using Non-Light Waves (AREA)
- Measurement Of Radiation (AREA)
Description
【発明の属する技術分野】
本発明は、デジタル放射線画像等の粒状性および画質に関する評価を行なう評価技術に係り、特にデジタルラジオグラフィ装置によるX線を用いたデジタル画像の粒状性および画質を評価するデジタル放射線画像の評価方法および評価装置に関する。
【0002】
【従来の技術】
従来より医療分野においては、X線等を利用した種々の医療診断が行われている。これらは、被写体を透過したX線等を検出してアナログ画像を生成し、この画像を医療診断に供するものである。例えば、古くからX線胸部単純撮影画像が広く用いられている。
また、近年DR(デジタル・ラジオグラフィ)、CT(コンピュータ断層像撮影)、MR(核磁気共鳴像撮影)等の各種デジタル放射線画像が、例えばCRTディスプレイに表示されたり、プリンタ等によって出力され、医療現場において病気の診断に利用されている。
このようなデジタル放射線画像の中でも、特にX線を用いたCR(コンピューテッド・ラジオグラフィ)システムにおけるデジタル放射線画像(デジタルラジオグラフィ)が医療分野では広く普及している。
【0003】
このCRシステムでは、輝尽性蛍光体を塗布したX線検出器(蓄積蛍光体シート)は、遮光性の高いカセッテと呼ばれるケースに入れられた状態で被写体を透過したX線の照射を受ける。そして、カセッテから蓄積蛍光体シートを取り出し、レーザ光線を用いた光学系で蓄積蛍光体シートの撮影面全体がレーザ光線で走査される(スキャンされる)。その際、蓄積蛍光体シートの輝尽性蛍光体は、レーザ光線の照射を受けると、吸収したX線量に比例して輝尽発光するため、この発光した光を集光ガイドで集めて、フォトセンサーで検出し、その後デジタル変換してデジタル画像データを得ることができる。このデジタル画像データは、階調処理や空間周波数強調などの各処理が行われ、処理後の画像がCRTディスプレイに表示されたり、プリンタ等によって出力される。
【0004】
このようなCRシステムで得られるデジタル画像は、デジタルデータを基本とするため、上記階調処理や空間周波数強調等の様々なフィルタを設けて画像を自由に修正することが可能である。
例えば、医療用X線画像は、患者の被曝の点からX線の照射線量を多くすることは好ましくないため、少ない照射線量で撮影し、その後得られたデジタル画像について階調処理を行うことで、X線の照射線量に依存せずに同等の階調、濃度の画像を作成することが可能である。
【0005】
デジタル画像の粒状性については、X線のゆらぎに起因するもの、輝尽蛍光体の発光する光のゆらぎに起因するものおよび蓄積蛍光体シートの構造に起因するもの等があげられる。デジタル画像は、撮影する目的に応じてX線の照射線量を広い範囲で自由に選択できるが、デジタル画像の粒状性や画質は、X線の照射線量に大きく依存するため、デジタル画像の粒状性や画質を検査、評価するには、一定のX線照射条件のもとで比較する必要がある。しかし、毎回の検査、評価時に正確に一定のX線照射量を照射することは困難であり、その照射線量の差が検査、評価値の誤差の原因となる。
【0006】
そのため、このようなCRシステムでは、画質や画像の粒状性の程度を絶えず把握しまた管理しつつ、撮影画像によって適切かつきめ細かな診断が安定的に行えるように、CRシステムによる画像の粒状性や画質を所定のX線照射線量に対して知ることが必要である。
【0007】
【発明が解決しようとする課題】
従来は、CRシステムの画像の画質や粒状性について、X線を所定の照射線量を定めて画像を得ようとしても、ぴったりとX線照射線量を合わせることはできなかったため、その照射線量の誤差が粒状性や画質の誤差の原因になっていた。
従って、CRシステムの画像の粒状性や画質の評価値を正しく管理することを難しくしていた。このようなデジタル放射線画像の粒状性や画質の評価の問題は、X線に限られず、α線やβ線等の各種放射線による画像においても困難である。
【0008】
そこで、本発明は、デジタル放射線画像の粒状性や画質について所定の照射線量での評価値を精度よく得ることができないといった問題点を解消し、正確な所定の照射線量に対するデジタル画像の粒状性や画質の評価値を正確に得ることができる評価方法および評価装置を提供することを課題とする。
【0009】
【課題を解決するための手段】
上記目的を達成するために、本発明は、X線を含む放射線を用いて放射線検出器に記録される画像に所定の処理を施してデジタル画像を得、このデジタル画像の粒状性または画質に関する評価を行うに際し、前記放射線の照射線量を変えて記録される画像を2画像以上得、これらの画像に基づいたデジタル画像の粒状性または画質に関する評価値を画像ごとに2つ以上得、この2つ以上の評価値から、前記放射線の照射線量に対する評価値の補正曲線を定め、この補正曲線を用いて、所定の照射線量におけるデジタル画像の粒状性または画質の評価値を得ることを特徴とするデジタル放射線画像の評価方法を提供するものである。
【0010】
その際、前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を用いて記録された画像に所定の処理を施して得られるデジタル放射線画像であって、
前記評価値は、デジタル画像の粒状性を評価するためのRMS粒状値の自乗値あるいはウィーナースペクトルの各周波数の値であり、照射線量Xに対して前記補正値の補正曲線は、A/X+B(AおよびBは定数)で表され、この補正曲線の定数AおよびBを前記2つ以上の評価値から定めることで補正曲線を定めることが好ましく、また、前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を用いて記録された画像に所定の処理を施して得られるデジタル放射線画像であって、
前記評価値は、デジタル画像の画質を評価するNEQ値であり、照射線量Xに対して前記補正値の補正曲線は、1/(C/X+D)(CおよびDは定数)で表され、この補正曲線の定数CおよびDを前記2つ以上の評価値から定めることで補正曲線を定めることが好ましい。
さらに、前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を塗布した放射線検出器上に形成された画像に所定の処理を施して得られるデジタル放射線画像であって、前記放射線検出器を少なくとも2領域に分けて放射線の照射線量を変えて分割撮影を行い、放射線の照射線量の異なる2画像以上を得、上記記載のデジタル放射線画像の評価方法を行うことが好ましい。その際、前記デジタル放射線画像の画像データは、対数変換されたものを用いるのが好ましい。
【0011】
また、本発明は、X線を含む放射線の照射線量に応じて輝尽発光する蛍光材料の塗布された放射線検出器上に記録された画像に所定の処理を施して得られるデジタル放射線画像を評価する評価装置であって、照射線量を変えて少なくとも2つ以上の画像を同一の放射線検出器に分割撮影し、あるいは別個の放射線検出器に撮影し、この撮影された画像に所定の処理を施してデジタル放射線画像を得るデジタル画像取得手段と、このデジタル画像取得手段によって得られたデジタル画像の各々から粒状性を評価するためのRMS粒状値の自乗値あるいはウィーナースペクトルのスペクトル値を得、この得られた値から照射線量Xに対してA/X+B(AおよびBは定数)なる補正曲線の定数AおよびBを定めることによって補正曲線を決定し、この決定された補正曲線から所定の照射線量における粒状性の評価値を算出する粒状性評価手段、および、前記デジタル画像取得手段によって得られたデジタル画像の各々から画質を評価するためのNEQ値を得、この得られた値から照射線量Xに対して1/(C/X+D)(CおよびDは定数)なる補正曲線の定数CおよびDを定めることによって補正曲線を決定し、この決定された補正曲線から所定の照射線量における画質の評価値を算出する画質評価手段のうちの少なくとも一つとを備えたことを特徴とするデジタル放射線画像の評価装置を提供するものである。
【0012】
【発明の実施の形態】
以下、本発明のデジタル放射線画像の評価方法を実施する本発明のデジタル放射線画像の評価装置について、添付の図面に示される好適実施例を基に詳細に説明する。なお以下の説明では、X線を例とした放射線画像に関するものであるが、X線による画像に限られず、α線やβ線等の各種放射線による画像であってもよい。
【0013】
図1に、本発明のデジタル放射線画像評価装置の好適実施例であるデジタルX線画像評価装置10を示す。
デジタルX線画像評価装置10は、デジタルX線画像診断システムで得られるデジタル画像の粒状性や画質に関する評価を行なう装置であって、X線の照射されたX線検出器である蓄積蛍光体シート上に記録された画像からデジタル画像データを得るデジタル画像データ取得手段12と、画像の粒状性の評価をRMS粒状値の自乗値あるいはウィーナースペクトルのスペクトル値で行なう粒状性評価手段14と、デジタル画像の信号対雑音比を意味するNEQ値によって画質を評価する画質評価手段16と、得られたデジタル画像データからデジタル画像を表示し、粒状性評価値や画質評価値の算出を行なうグラフ表示を行なうモニタ18と、評価値を記録保存、管理する評価値記録部20とを主に有する。
【0014】
まず、デジタルX線画像評価装置10で処理されるX線検出器である蓄積蛍光体シートは、X線の照射線量に応じて輝尽発光する蛍光材料(例えば、フッ化ハロゲン化バリウム)が塗布されたシートであり、この蓄積蛍光体シートに記録される画像は、蓄積蛍光体シートの記録される領域を少なくとも2領域に分けてX線の照射線量を変えて分割撮影を行ったものであり、被写体に遮られることなく、X線が蓄積蛍光体シートに均一に照射された画像である。そのため、この画像には、X線のゆらぎや蓄積蛍光体シートの構造上の雑音成分等が含まれる。
なお本発明では、1つの蓄積蛍光体シート上に、X線の照射線量を変えて分割撮影するが、これに限られず、別個の蓄積蛍光体シート上に、X線の照射線量を変えた画像を別々に得てもよい。
【0015】
デジタル画像データ取得手段12は、上記蓄積蛍光体シート上にX線が照射され記録された画像を、所定のレーザビーム、例えば、He−Neレーザビーム等を走査して、X線の照射線量に応じて蛍光材料を輝尽発光させ、この発光した光をフォトセンサーを用いて信号化する。
【0016】
すなわち、レーザビームが光源から射出され、ビームエクスパンダによってビーム径を拡径され、ガルバノメータミラー等の光偏向器によって一次元的に(主走査方向に)走査される。主走査方向に走査されたレーザビームは、fθレンズによって、所定の読み取り位置に所定のビーム径で結像される。一方、放射線画像を蓄積記録された蛍光体シートは、ローラおよびエンドレスベルトからなる副走査搬送手段によって、その表面を読み取り位置に保たれつつ、主走査方向と直交する副走査方向に副走査搬送される。
前述のように、レーザビームは、主走査方向に走査されているので、副走査方向に搬送される蛍光体シートは、結果的に、レーザビームによって2次元的に全面を走査される。
【0017】
蛍光体シートのレーザービームによって走査された位置からは、蓄積記録された放射線画像に応じた輝尽発光光が射出される。この輝尽発光光は、光ガイドによって伝播されて、フィルタによって輝尽発光光以外の成分を除去され、フォトマルチプライヤによって、光電的に読み取られる。
このフォトマルチプライヤからの出力信号は、放射線画像の画像信号として出力され、Log変換を行い、その後AD変換を行なって、デジタル画像の画素単位の階調値であるデジタル画像データを得る。
ここで、蓄積蛍光体シートは上述したように記録領域が少なくとも2分割され、X線の照射線量を変えて、分割撮影されているが、蓄積蛍光体シート全体がレーザービームによって走査され読み取られ、すなわち、X線の照射線量を変えて分割撮影された画像が一度に読み取られる。
モニタ18には、蓄積蛍光体シート全体に渡って撮影されて得られたデジタル画像が表示され、画像分割され、粒状性評価値や画質評価値を各々得ることができる。
X線の照射線量を変えた画像を別個の蓄積蛍光体シートに得た場合は、別々にデジタル画像が表示される。
【0018】
粒状性評価手段14は、本発明の特徴とする部分であり、デジタルX線画像診断システムの画像の粒状性を評価する部分であり、その評価はあるX線照射線量におけるRMS粒状値の自乗値あるいはウィーナースペクトルの値(以降、WS値という)で行なう。
RMS粒状値の自乗値とは、デジタル画像データ取得手段12で得られたデジタル画像データのうち、画像の所定方向の1ラインの画像データ列または複数ラインの平均の画像データ列である画素単位の階調値について、上記所定方向の1ラインまたは複数ラインの階調値の平均値からの差分の自乗和をとった値、すなわちRMS変動値(root mean aquare variation)の自乗値をいう。このようなRMS粒状値の自乗値は、画像の方向それぞれに関して得ることができ、本実施例では、画像の方向をデジタル画像データ取得手段12でデータを得る際、レーザビームを走査する方向を主走査方向、および、それに直交する方向を副走査方向とし、その各々の方向に対してRMS粒状値の自乗値を別個に設けている。
また、WS値とは、画像データである階調値をフーリエ変換したパワースペクトルの所定の空間周波数でのスペクトル値であり、画像の雑音成分を表す評価値である。それゆえ、ウィーナースペクトルのスペクトル値は、注目する空間周波数によって変化する値である。
【0019】
粒状性評価手段14では、RMS粒状値の自乗値は、通常の公知の方法で計算され、WS値は、画像データである階調値をフーリエ変換して求めることで得られる。
この得られた評価値はX線の照射線量によって変化するため、X線の照射線量を変えて得られた2つ以上の評価値からX線の照射線量に対する補正曲線を算出して、あるX線照射線量における評価値を得ることができる。
【0020】
画質評価手段16は、デジタル画像の信号対雑音比を意味するNEQ(Noise Equivalent Quanta )値を用いて画質を評価する部分で、本発明の特徴とする部分である。
NEQ値とは、以下の式で表され、ノイズ透過量子数と呼ばれるものである。
NEQ値 = γ2 ×(loge 2)2 ×MTF2 /(WS値)2
すなわち、MTF(Modulation Transfer Function)の自乗値をWS値で割った値に比例し、この値は、雑音成分が低く、すなわちWS値が小さく、画像の先鋭度が高い、すなわちMTFの自乗値が高い程、NEQ値は高くなり、画像の画質として良好であることを示すものである。
【0021】
画質評価手段16では、MTFの自乗値は、チャート法等の公知の方法で得、WS値は、、画像データである階調値について画像全体の階調値の平均値からの差分を取り、この差分データを空間周波数にフーリエ変換し、この変換によって得られるパワースペクトルの所定の空間周波数でのスペクトル値を求めることで得られ、NEQ値を求める。
この求められたNEQ値はX線の照射線量によって変化するため、X線の照射線量を変えて得られた2つ以上のNEQ値からX線の照射線量に対する補正曲線を算出して、あるX線照射線量における評価値を得ることができる。
【0022】
モニタ18は、デジタル画像データ取得手段12で得られた画像を表示する部分であり、粒状性評価手段14や画質評価手段16で得られた評価値等を表示する部分である。また、蓄積蛍光体シートを分割して、異なるX線照射線量で画像を得た場合、蓄積蛍光体シート全体の画像を得た後画像が分割され、粒状性評価値や画質評価値を各々得ることができる。
評価値記録部20は、粒状性評価手段14や画質評価手段16で得られた評価値を保存記録する部分である。
なお、本発明のデジタル画像の評価装置には、デジタルX線画像評価装置10におけるモニタ18のような画像表示装置が必ずしも設けられる必要はなく、デジタル画像の評価装置の外部に接続されてもよい。
【0023】
デジタルX線画像評価装置10は、以上のように構成される。
つぎに、デジタルX線画像評価装置10の作用について説明する。
【0024】
まず、X線検出器である蓄積蛍光体シートは、蓄積蛍光体シートの記録される領域を少なくとも2領域に分けられ、X線の照射線量を変えて分割撮影を行う。分割撮影の仕方は特に限定されず、例えば、分割撮影する領域のうち撮影されない領域に、X線を透過しない金属板、たとえば鉛板等で覆い、順次この金属板の覆いの位置を変えて、X線の照射線量を変えて照射を行なう。
【0025】
デジタル画像データ取得手段12では、X線の照射線量に応じて輝尽発光する蛍光材料が塗布された蓄積蛍光体シート上にX線が照射されたものを、所定のレーザビームで走査する。蛍光体シートのレーザービームによって走査された位置からは、蓄積記録された放射線画像に応じた輝尽発光光が射出される。この輝尽発光光は、光ガイドによって伝播されて、フィルタによって輝尽発光光以外の成分を除去され、フォトマルチプライヤによって、光電的に読み取られ信号化される。その後Log変換を行いAD変換を行なってデジタル画像データを得る。得られたデジタル画像データは、モニタ18に送られ、画像表示される。
蓄積蛍光体シートを分割して、異なるX線照射線量で分割撮影した蓄積蛍光体シート全体の画像を得た後分割される。予め蓄積蛍光体シートを分割して、異なるX線照射線量で分割撮影する方が2つの画像を得るのに容易だからである。
【0026】
つぎに、得られた画像について、粒状性評価が行なわれる。粒状性評価は、粒状性評価手段14で行なわれ、所定のX線照射線量での粒状性の評価値を得る。
まず、得られたデジタルデータから、RMS粒状値(以降、RMS値という)の自乗値あるいは、ある空間周波数でのWS値を算出する。RMS値の自乗値は、画像の一方向の所定のライン上の画素単位の階調値のRMS値の自乗値を通常の公知の方法によって算出する。WS値については、画像データである階調値について画像全体の階調値の平均値からの差分を取り、この差分データを空間周波数にフーリエ変換し、この変換によって得られるパワースペクトルの所定の空間周波数でのスペクトル値を求めることで得られる。RMS値の自乗値およびWS値は、X線照射線量の異なる分割された複数の画像から得られるため、画像毎にRMS値の自乗値あるいはWS値のデータが得られ、合計2つ以上のRMS値の自乗値またはWS値を得る。
【0027】
つぎに、X線照射線量を横軸に、RMS値の自乗値あるいはWS値を縦軸にした対数グラフに上記RMS値の自乗値またはWS値のデータをプロットする。X線は照射線量を変えて少なくとも2画像以上得ているので、2つ以上のRMS値の自乗値またはWS値のデータをグラフにプロットする。また、RMS値の自乗値は、主走査方向および副走査方向のそれぞれに分けてプロットすることができる。
【0028】
つぎに、このプロットされた少なくとも2つのデータに基づいて、RMS値の自乗値あるいはWS値の補正曲線を決定する。RMS値の自乗値やWS値の補正曲線の決定は、補正曲線が一般にX線照射線量Xに対して、A/X+B(AおよびBは定数)で表されているので、補正曲線がプロットされた少なくとも2点のデータの上を通るように、定数AおよびBを定めることによって行なわれる。3点以上のプロットに対して補正曲線を定めるには、最小自乗によるカーブフィット方法等公知の方法が用いられる。また、その際データに重み付けをもちいてもよい。
特に、補正曲線がプロットされたグラフ上の2点のデータの上を通るように、補正曲線の定数AおよびBを定めた場合、プロットされたグラフ上の3点以上のデータから得られる補正曲線と大きな差はなく、従来プロットされた2点間を直線的に内挿し所定量の照射線量に対して評価値を精度よく得ることができなかった問題を同じ2点のデータに基づいて解決することができる。
【0029】
ところで、このような補正曲線が一定の関数の形式で表されることは、本出願の発明者らが見出したものであり、従来のアナログの放射線画像にはないデジタル放射線画像に特有のものである。
つまり、X線の場合、照射線量Xに対して、A/X+B(AおよびBは定数)の関数で表され、X線の照射線量Xが多い領域では、照射線量に係わらず一定値を有するが、照射線量Xが少ない領域では照射線量が小さくなるに連れて粒状性が急激に悪化する特徴を示す。X線の照射線量が小さいとX線の粒状性が悪化するのは、一般にデジタルX線画像診断システムでの画像の粒状性は、X線のゆらぎや輝尽蛍光材の発光する光のゆらぎや蓄積蛍光体シート自体の構造ノイズに起因するが、X線の照射線量が小さい場合、粒状性の評価にとってこのX線のゆらぎに起因するノイズが支配的となるからである。
なお、このような特徴はX線に限らず、α線やβ線等の各種放射線にも見られる。そのため、本発明は、α線やβ線等の各種放射線にも適用される。
【0030】
図2には、画像の粒状性を評価する評価方法を説明するグラフが示され、画像の主走査方向および副走査方向について、X線照射線量の異なる2つのRMS値の自乗値が対数変換されてプロットされ(主走査方向は記号「+」で、副走査方向は記号「口」でプロットされている)、さらに、未知数AおよびBを定めることによって定まった補正曲線A/X+Bを主走査方向および副走査方向の各々について重ね書きされている。また、補正曲線A/X+Bの構成要素A/Xの項およびBの項についても、主走査方向および副走査方向の各々について重ね書きされている。
得られた補正曲線A/X+Bの定数AおよびBは、図2のグラフの右上欄に示され、例えば主走査方向の定数Aは、0.62×10-4、定数Bは0.10×10-4であることが判る。
図2ではRMS値の自乗値の補正曲線を実施例として挙げたが、WS値の補正曲線を定めてもよいのは勿論である。
【0031】
表1の測定値の欄には、図2においてRMS値の自乗値の補正曲線に用いられた2つのデータが、主走査方向および副走査方向について、X線照射線量とともに示されている。一方、計算値の欄には、測定に使用されたX線照射線量でのRMS値の自乗値の計算値が補正曲線から算出されている。表1では、補正曲線A/X+Bの未知定数AおよびBは2つのデータから求められているので、RMS値の自乗値の実測値と計算値との間には差異は生じない。さらに計算値の下欄には、主走査方向および副走査方向のRMS値の自乗値が、X線の照射線量0.1(mR)、0.5(mR)、1.0(mR)、5.0(mR)、10(mR)、100(mR)の各々において示される。
この計算値から、所定のX線照射線量、例えば0.5(mR)や1.0(mR)でのRMS値の自乗値を得、所定のX線照射線量での評価値を得る。
上記例では、RMS値の自乗値であるが、WS値であってもよい。
【0032】
【表1】
【0033】
つぎに、得られた画像についての画質の評価が行なわれる。画質の評価は、画質評価手段16で行なわれ、所定のX線照射線量での画質の評価値であるNEQ値を得る。
まず、得られたデジタルデータから、画質の評価値であるNEQ値を算出する。その際、蓄積蛍光シートの領域を分割して、X線照射線量を変えて分割撮影された複数の画像から、NEQ値のデータを2つ以上得る。
NEQ値は、MTFの自乗値をWS値で割った値に比例し、MTF値およびWS値はそれぞれ空間周波数に依存するため、空間周波数ごとにNEQ値は定まる。MTFは、従来の公知の例による方法、例えばチャート法またはスリット法を用いる。WS値は、、画像データである階調値について画像全体の階調値の平均値からの差分を取り、この差分データを空間周波数にフーリエ変換し、この変換によって得られるパワースペクトルの所定の空間周波数でのスペクトル値を求めることで得られる。ここでウィーナースペクトルは、上述したように、X線照射線量Xに対して、補正曲線A/X+Bの形でスペクトル値を変えるため、ウィーナースペクトルのスペクトル値を分母に持つNEQ値では、X線の照射線量に対して、1/(C/X+D)(CおよびDは定数)の形で変化する。そこで、X線の照射線量の異なる画像のNEQ値を2つ以上得、これより補正曲線1/(C/X+D)の未知定数CおよびDを定め、補正曲線を算出する。
【0034】
補正曲線の未知数CおよびDの算出方法は、RMS値の自乗値について上述した場合と同じように、X線の照射線量を横軸に、NEQ値を縦軸にした対数変換されたグラフにデータをプロットする。X線は照射線量を変えて少なくとも2画像以上得ているので、2つ以上のNEQ値のデータをグラフにプロットする。その際、NEQ値は、空間周波数に依存するので、一定の間隔の空間周波数、例えば0.5サイクル/mm毎にNEQ値をプロットする。
つぎに、このプロットされた少なくとも2つのデータに基づいて、NEQ値の補正曲線を決定する。補正曲線の決定は、補正曲線1/(C/X+D)がプロットされた少なくとも2点のデータの上を通るように、定数CおよびDを定めることによって行なわれる。3点以上のプロットに対して補正曲線を定めるには、最小自乗によるカーブフィット方法等公知の方法が用いられる。また、その際データに重み付けをもちいてもよい。
特に、補正曲線がプロットされたグラフ上の2点のデータの上を通るように、補正曲線の定数CおよびDを定めた場合、プロットされたグラフ上の3点以上のデータから得られる補正曲線と大きな差はなく、従来プロットされた2点間を直線的に内挿し所定量の照射線量に対して評価値を精度よく得ることができなかった問題を同じ2点のデータに基づいて解決することができる。
【0035】
図3には、画像の画質を評価する評価方法を説明するグラフが示され、X線照射線量の異なる2つのNEQ値が空間周波数0.5(サイクル/mm)から3.0(サイクル/mm)まで0.5(サイクル/mm)毎にプロットされている。また、未知の定数CおよびDを定めることで定まったNEQ値の補正曲線1/(C/X+D)が重ね書きされている。得られた補正曲線1/(C/X+D)の定数CおよびDは、図3のグラフの右側の欄に示され、例えば定数Cは、空間周波数0.5(サイクル/mm)では1.480×10-5、定数Dは1.743×10-6であることが判る。このようにしてNEQ値の補正曲線を定めることができる。
【0036】
表2の測定値の欄には、図3で補正曲線に用いられたX線照射線量0.592(mR)および9.237(mR)でのNEQ値が空間周波数0.5(サイクル/mm)毎に示されている。一方、その下欄には、求められた補正曲線でのX線の照射線量0.1(mR)、0.5(mR)、1.0(mR)、5.0(mR)、10(mR)、100(mR)の各々におけるNEQ値の計算値が示される。この計算値から、所定のX線照射線量、例えば0.5(mR)や1.0(mR)でのNEQ値の評価値を得る。
【0037】
【表2】
【0038】
このようにして、粒状性評価値であるRMS値の自乗値あるいはWS値および画質評価値であるNEQ値をそれぞれ粒状性評価手段14および画質評価手段16で得る際、図2や図3に示されるグラフや表1や表2に示される一覧表が適宜モニタ18に表示される。
このモニタ18の表示結果を見ることによって、オペレータは、画像の粒状性や画質を評価し、管理することができる。また、その結果を評価値記憶部20に記憶させ、一定期間毎に評価を行なうことによって、画像の粒状性や画質についての経時変化を記録することができる。
なお、本実施例では、得られた画像データは対数変換されたものを用いることで、対数グラフ上で容易に補正曲線を定めることができるが、必ずしも対数変換される必要はなく、対数変換されない画像データより補正曲線を定めてもよい。
【0039】
以上、X線を用いたデジタル放射線画像の評価方法は、デジタル画像を得るシステムに特有の雑音成分を原因とする粒状性の補正曲線を定めることによって、所定の照射線量での粒状性および画質の評価を行なうものであるが、通常のカメラで撮影されたフィルムに対して、可視光線を照射してフィルムの撮影された画像情報を担持する投影光を得、この投影光をフォトセンサーを用いて信号化し、その後得られた信号をLog変換およびAD変換をしてデジタル画像データを得る撮影フィルムのデジタル画像においても、X線のデジタル画像と同様の評価方法を用いて、デジタル画像の粒状性や画質を評価することができる。
【0040】
その他、本発明は上記実施例に限定はされず、本発明の要旨を逸脱しない範囲において、各種の改良および変更を行ってもよいのはもちろんである。
【0041】
【発明の効果】
以上、詳細に説明したように、本発明において、放射線の照射線量を変えて得られる2つ以上の画像に基づいたデジタル画像の粒状性または画質に関する評価値を2つ以上得、この2つ以上の評価値から照射線量に対する評価値の補正曲線を定め、この補正曲線を用いて、所定の照射線量におけるデジタル画像の粒状性または画質の評価値を得るため、正確な所定の照射線量に対して粒状性や画質の評価値を正確に得ることができ、その結果、例えば医療診断において安定したきめの細かな診断を行なうことができる。
【図面の簡単な説明】
【図1】 本発明の画像評価方法を実施する本発明の画像評価装置の一実施例の概略ブロック図である。
【図2】 画像の粒状性を評価する評価方法を説明する説明図である。
【図3】 画像の画質の評価方法を説明する説明図である。
【符号の説明】
10 デジタルX線画像評価装置
12 デジタル画像データ取得手段
14 粒状性評価手段
16 画質評価手段
18 モニタ
20 評価値記録部
Claims (6)
- X線を含む放射線を用いて放射線検出器に記録される画像に所定の処理を施してデジタル画像を得、このデジタル画像の粒状性または画質に関する評価を行うに際し、
前記放射線の照射線量を変えて記録される画像を2画像以上得、これらの画像に基づいたデジタル画像の粒状性または画質に関する評価値を画像ごとに2つ以上得、
この2つ以上の評価値から、前記放射線の照射線量に対する評価値の補正曲線を定め、
この補正曲線を用いて、所定の照射線量におけるデジタル画像の粒状性または画質の評価値を得ることを特徴とするデジタル放射線画像の評価方法。 - 前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を用いて記録された画像に所定の処理を施して得られるデジタル放射線画像であって、
前記評価値は、デジタル画像の粒状性を評価するためのRMS粒状値の自乗値あるいはウィーナースペクトルのスペクトル値であり、照射線量Xに対して前記補正値の補正曲線は、A/X+B(AおよびBは定数)で表され、この補正曲線の定数AおよびBを前記2つ以上の評価値から定めることで補正曲線を定める請求項1に記載のデジタル放射線画像の評価方法。 - 前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を用いて記録された画像に所定の処理を施して得られるデジタル放射線画像であって、
前記評価値は、デジタル画像の画質を評価するNEQ値であり、照射線量Xに対して前記補正値の補正曲線は、1/(C/X+D)(CおよびDは定数)で表され、この補正曲線の定数CおよびDを前記2つ以上の評価値から定めることで補正曲線を定める請求項1に記載のデジタル放射線画像の評価方法。 - 前記デジタル画像は、放射線の照射線量に応じて輝尽発光する蛍光材料を塗布した放射線検出器上に形成された画像に所定の処理を施して得られるデジタル放射線画像であって、
前記放射線検出器を少なくとも2領域に分けて放射線の照射線量を変えて分割撮影を行い、放射線の照射線量の異なる2画像以上を得、請求項1〜3のいずれかに記載のデジタル放射線画像の評価方法を行うことを特徴とするデジタル放射線画像の評価方法。 - 前記デジタル放射線画像の画像データは、対数変換されたものを用いる請求項1〜4のいずれかに記載のデジタル放射線画像の評価方法。
- X線を含む放射線の照射線量に応じて輝尽発光する蛍光材料の塗布された放射線検出器上に記録された画像に所定の処理を施して得られるデジタル放射線画像を評価する評価装置であって、
照射線量を変えて少なくとも2つ以上の画像を同一の放射線検出器に分割撮影し、あるいは別個の放射線検出器に撮影し、この撮影された画像に所定の処理を施してデジタル放射線画像を得るデジタル画像取得手段と、
このデジタル画像取得手段によって得られたデジタル画像の各々から粒状性を評価するためのRMS粒状値の自乗値あるいはウィーナースペクトルのスペクトル値を得、この得られた値から照射線量Xに対してA/X+B(AおよびBは定数)なる補正曲線の定数AおよびBを定めることによって補正曲線を決定し、この決定された補正曲線から所定の照射線量における粒状性の評価値を算出する粒状性評価手段、
および、
前記デジタル画像取得手段によって得られたデジタル画像の各々から画質を評価するためのNEQ値を得、この得られた値から照射線量Xに対して1/(C/X+D)(CおよびDは定数)なる補正曲線の定数CおよびDを定めることによって補正曲線を決定し、この決定された補正曲線から所定の照射線量における画質の評価値を算出する画質評価手段のうちの少なくとも一つとを備えたことを特徴とするデジタル放射線画像の評価装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP11223599A JP4159701B2 (ja) | 1999-04-20 | 1999-04-20 | デジタル放射線画像の評価方法および評価装置 |
US09/553,179 US6801647B1 (en) | 1999-04-20 | 2000-04-20 | Method and apparatus for evaluating digital radiation image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP11223599A JP4159701B2 (ja) | 1999-04-20 | 1999-04-20 | デジタル放射線画像の評価方法および評価装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000305208A JP2000305208A (ja) | 2000-11-02 |
JP4159701B2 true JP4159701B2 (ja) | 2008-10-01 |
Family
ID=14581638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP11223599A Expired - Fee Related JP4159701B2 (ja) | 1999-04-20 | 1999-04-20 | デジタル放射線画像の評価方法および評価装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6801647B1 (ja) |
JP (1) | JP4159701B2 (ja) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7123759B1 (en) * | 1999-06-30 | 2006-10-17 | Fuji Photo Film Co., Ltd. | Image transfer and output method and system using the same |
KR100425379B1 (ko) * | 2000-11-13 | 2004-03-30 | 남상희 | 디지털 방사선 검출장치 및 방사선 검출방법 |
JP3945976B2 (ja) * | 2000-11-21 | 2007-07-18 | 富士フイルム株式会社 | 画像表示方法及び装置 |
JP2005028114A (ja) * | 2003-06-18 | 2005-02-03 | Canon Inc | 放射線撮影装置及び放射線撮影方法 |
JP4340889B2 (ja) * | 2004-07-15 | 2009-10-07 | 横河電機株式会社 | 検査装置 |
US7071690B2 (en) * | 2004-11-23 | 2006-07-04 | The Board Of Trustees Of The Leland Stanford Junior University | Reduction of blurring in view angle tilting MRI |
US7991242B2 (en) | 2005-05-11 | 2011-08-02 | Optosecurity Inc. | Apparatus, method and system for screening receptacles and persons, having image distortion correction functionality |
WO2006119603A1 (en) | 2005-05-11 | 2006-11-16 | Optosecurity Inc. | Method and system for screening luggage items, cargo containers or persons |
US7899232B2 (en) | 2006-05-11 | 2011-03-01 | Optosecurity Inc. | Method and apparatus for providing threat image projection (TIP) in a luggage screening system, and luggage screening system implementing same |
US8494210B2 (en) | 2007-03-30 | 2013-07-23 | Optosecurity Inc. | User interface for use in security screening providing image enhancement capabilities and apparatus for implementing same |
WO2009128213A1 (ja) * | 2008-04-14 | 2009-10-22 | 株式会社日立メディコ | 医用診断装置および医用診断装置の画質改善方法 |
CA3069211C (en) | 2011-09-07 | 2022-06-21 | Rapiscan Systems, Inc. | X-ray inspection system that integrates manifest data with imaging/detection processing |
JP5774447B2 (ja) | 2011-10-28 | 2015-09-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置および被曝線量計算方法並びにプログラム |
JP5647639B2 (ja) * | 2012-03-19 | 2015-01-07 | 富士フイルム株式会社 | 放射線撮影情報管理システムおよび方法、並びにプログラム |
CN109074889B (zh) | 2016-02-22 | 2022-10-18 | 拉皮斯坎系统股份有限公司 | 用于检测货物中的危险品和违禁品的系统和方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3051605B2 (ja) * | 1993-07-14 | 2000-06-12 | 富士写真フイルム株式会社 | 放射線画像形成方法とそれに用いるハロゲン化銀写真材料 |
US5420441A (en) * | 1993-11-23 | 1995-05-30 | Eastman Kodak Company | Automated technique for calibrating a storage phosphor reader |
US5574768A (en) * | 1994-05-26 | 1996-11-12 | The University Of Sydney | Screen film cassette |
-
1999
- 1999-04-20 JP JP11223599A patent/JP4159701B2/ja not_active Expired - Fee Related
-
2000
- 2000-04-20 US US09/553,179 patent/US6801647B1/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2000305208A (ja) | 2000-11-02 |
US6801647B1 (en) | 2004-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4159701B2 (ja) | デジタル放射線画像の評価方法および評価装置 | |
JPH03170129A (ja) | 画像処理装置 | |
JP2009285356A (ja) | 医療用撮影システム、画像処理装置、画像処理方法、およびプログラム | |
JP2004105437A (ja) | 医用画像処理装置及び医用画像撮影システム | |
US6819734B2 (en) | Method for removing rings and partial rings in computed tomography images | |
Schaetzing | Management of pediatric radiation dose using Agfa computed radiography | |
JPH0214378A (ja) | 放射線画像読取再生装置 | |
JP4307877B2 (ja) | 画像処理装置及び画像処理方法 | |
JP2005028037A (ja) | 医用画像処理装置及び医用画像処理方法 | |
JP3270153B2 (ja) | コンピュータトモグラフ | |
JP2002152593A (ja) | エネルギーサブトラクション方法および装置並びに記録媒体 | |
JP2002034961A (ja) | 放射線撮影装置及び放射線撮影方法 | |
JP4258855B2 (ja) | 放射線画像撮影表示方法および放射線画像撮影表示装置 | |
JP2952519B2 (ja) | 放射線画像の階調変換装置 | |
JPH05217689A (ja) | X線撮影装置およびx線撮影方法 | |
JP2009082498A (ja) | 乳房用の放射線撮影装置、放射線撮影システム、情報処理装置、及び入射線量導出方法 | |
US6891965B2 (en) | Image display method and apparatus | |
JP2001120524A (ja) | 放射線画像処理装置 | |
JP2983421B2 (ja) | 骨計測方法及び装置 | |
Theodorakou et al. | A novel method for producing x-ray test objects and phantoms | |
JP2001149359A (ja) | 撮像装置、画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 | |
Seeram et al. | Computed radiography: physics and technology | |
JPH1048341A (ja) | 画像生成方法と装置 | |
US7262426B2 (en) | Method, apparatus and program for image data correction | |
JP2000157518A (ja) | 放射線画像処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040830 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061215 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080219 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080410 |
|
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: 20080708 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080716 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110725 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110725 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120725 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120725 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130725 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |