JP2000298076A - 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 - Google Patents
鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置Info
- Publication number
- JP2000298076A JP2000298076A JP10766499A JP10766499A JP2000298076A JP 2000298076 A JP2000298076 A JP 2000298076A JP 10766499 A JP10766499 A JP 10766499A JP 10766499 A JP10766499 A JP 10766499A JP 2000298076 A JP2000298076 A JP 2000298076A
- Authority
- JP
- Japan
- Prior art keywords
- ctf
- sharpness
- rectangular wave
- data
- wave chart
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title abstract description 19
- 238000011156 evaluation Methods 0.000 claims description 51
- 238000012545 processing Methods 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012854 evaluation process Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 15
- 238000002601 radiography Methods 0.000 description 17
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 12
- 238000005070 sampling Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000013500 data storage Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
Landscapes
- Measurement Of Radiation (AREA)
Abstract
(57)【要約】
【課題】チャート法を用いるＣＴＦ測定においても、高
精度なＣＴＦの測定を可能とした、鮮鋭度評価方法およ
びこの方法を用いる鮮鋭度評価装置を提供すること。 【解決手段】矩形波チャートを露光して得られたデジタ
ルデータからＣＴＦを計算し、これに基づいて画像の鮮
鋭度の評価を行う方法であって、前記矩形波チャートを
露光して得られたパターンの振幅を計算する際、ピーク
値の中央値または該中央値±ｎ（ここで、ｎは同一周波
数での繰り返し回数をＮとしたとき、ｎ＜Ｎ／２−１で
表わされる数である）の平均値、あるいは、最大，最小
値からｍ個を除いた残りのデータの平均値を用いること
を特徴とする鮮鋭度評価方法および鮮鋭度評価装置。
精度なＣＴＦの測定を可能とした、鮮鋭度評価方法およ
びこの方法を用いる鮮鋭度評価装置を提供すること。 【解決手段】矩形波チャートを露光して得られたデジタ
ルデータからＣＴＦを計算し、これに基づいて画像の鮮
鋭度の評価を行う方法であって、前記矩形波チャートを
露光して得られたパターンの振幅を計算する際、ピーク
値の中央値または該中央値±ｎ（ここで、ｎは同一周波
数での繰り返し回数をＮとしたとき、ｎ＜Ｎ／２−１で
表わされる数である）の平均値、あるいは、最大，最小
値からｍ個を除いた残りのデータの平均値を用いること
を特徴とする鮮鋭度評価方法および鮮鋭度評価装置。
Description
【０００１】
【発明の属する技術分野】本発明は、鮮鋭度評価方法お
よびこれを用いる鮮鋭度評価装置に関し、より詳細に
は、ＤＲ（ＤigitalＲadiography），ＣＲ（Ｃomputed
Ｒadiography）をはじめとする各種ラジオグラフィー用
の装置において取り扱われる画像（以下、デジタルラジ
オグラフィー画像という）の鮮鋭度評価方法およびこの
方法を利用する鮮鋭度評価装置に関する。
よびこれを用いる鮮鋭度評価装置に関し、より詳細に
は、ＤＲ（ＤigitalＲadiography），ＣＲ（Ｃomputed
Ｒadiography）をはじめとする各種ラジオグラフィー用
の装置において取り扱われる画像（以下、デジタルラジ
オグラフィー画像という）の鮮鋭度評価方法およびこの
方法を利用する鮮鋭度評価装置に関する。
【０００２】
【従来の技術】従来から、増感紙／フィルム系の解像特
性の評価には、ＭＴＦ（ＭodulationＴransfer Ｆuncti
on ）による方法が、有力な手段として頻繁に用いられ
ているが、このアナログ画像系の手法は、デジタル画像
系にも適用可能である。ただし、デジタル画像系では、
離散的なデータ・サンプリングに起因するエリアシング
（折り返し誤差）のために、正確なＭＴＦが計算される
とは限らず、ＭＴＦの計算結果の解釈には、充分な注意
が必要である。
性の評価には、ＭＴＦ（ＭodulationＴransfer Ｆuncti
on ）による方法が、有力な手段として頻繁に用いられ
ているが、このアナログ画像系の手法は、デジタル画像
系にも適用可能である。ただし、デジタル画像系では、
離散的なデータ・サンプリングに起因するエリアシング
（折り返し誤差）のために、正確なＭＴＦが計算される
とは限らず、ＭＴＦの計算結果の解釈には、充分な注意
が必要である。
【０００３】これに関連して、藤田等は、ＩＩ（Ｉmage
Ｉntensifier）／ＴＶデジタル撮像系の解像特性の研究
において、上述のエリアシングの影響を含まない「プリ
サンプリング」ＭＴＦを測定している。この「プリサン
プリング」ＭＴＦは、幾何学的不鋭（アンシャープネ
ス），検出器の不鋭およびサンプリング・アパーチャの
不鋭を含んだ「デジタル系に固有」なＭＴＦである（医
用画像情報学会雑誌、Ｖol.6，Ｎo.1(1989) 参照）。
Ｉntensifier）／ＴＶデジタル撮像系の解像特性の研究
において、上述のエリアシングの影響を含まない「プリ
サンプリング」ＭＴＦを測定している。この「プリサン
プリング」ＭＴＦは、幾何学的不鋭（アンシャープネ
ス），検出器の不鋭およびサンプリング・アパーチャの
不鋭を含んだ「デジタル系に固有」なＭＴＦである（医
用画像情報学会雑誌、Ｖol.6，Ｎo.1(1989) 参照）。
【０００４】前述のように、ＭＴＦは鮮鋭度などの解像
特性の評価によく利用されるが、その測定には、正弦波
チャートが必要である。この正弦波チャートは作成が困
難なために、これに代わって、矩形波チャートを用いる
ＣＴＦ（Ｃontrast Ｔransfer Ｆunction ）の測定が行
われることが多い。測定されたＣＴＦは、Ｃoltmanの換
算式を用いてＭＴＦに変換することができる。
特性の評価によく利用されるが、その測定には、正弦波
チャートが必要である。この正弦波チャートは作成が困
難なために、これに代わって、矩形波チャートを用いる
ＣＴＦ（Ｃontrast Ｔransfer Ｆunction ）の測定が行
われることが多い。測定されたＣＴＦは、Ｃoltmanの換
算式を用いてＭＴＦに変換することができる。
【０００５】
【発明が解決しようとする課題】この矩形波チャートを
用いるＣＴＦの測定において問題となるのは、測定対象
である矩形波チャートを露光した蓄積性蛍光体シート
（以下、単に蛍光体シートという）の測定値にはノイズ
が入りやすく、このノイズの影響を受けやすいことであ
る。すなわち、このようなノイズをそのままにすると、
蛍光体シートを読み取る際に、このノイズの影響を受け
てデータが変化してしまい、ＣＴＦの測定精度が著しく
低下してしまうことである。
用いるＣＴＦの測定において問題となるのは、測定対象
である矩形波チャートを露光した蓄積性蛍光体シート
（以下、単に蛍光体シートという）の測定値にはノイズ
が入りやすく、このノイズの影響を受けやすいことであ
る。すなわち、このようなノイズをそのままにすると、
蛍光体シートを読み取る際に、このノイズの影響を受け
てデータが変化してしまい、ＣＴＦの測定精度が著しく
低下してしまうことである。
【０００６】上述の問題は、ＣＴＦの測定を行う対象
が、診断用途などの、厳密な表示性能が要求される用途
に使用されるデジタルラジオグラフィー装置である場合
には、極めて重大な問題ということになる。
が、診断用途などの、厳密な表示性能が要求される用途
に使用されるデジタルラジオグラフィー装置である場合
には、極めて重大な問題ということになる。
【０００７】本発明は、上記事情に鑑みてなされたもの
で、その目的とするところは、従来の技術における問題
を解消し、チャート法を用いるＣＴＦ測定においても、
高精度なＣＴＦの測定を可能とした、鮮鋭度評価方法お
よびこの方法を用いる鮮鋭度評価装置を提供することに
ある。
で、その目的とするところは、従来の技術における問題
を解消し、チャート法を用いるＣＴＦ測定においても、
高精度なＣＴＦの測定を可能とした、鮮鋭度評価方法お
よびこの方法を用いる鮮鋭度評価装置を提供することに
ある。
【０００８】
【課題を解決するための手段】上記目的を達成するた
め、本発明に係る鮮鋭度評価方法は、矩形波チャートを
露光して得られたデジタルデータからＣＴＦを計算し、
これに基づいて画像の鮮鋭度の評価を行う際に、前記矩
形波チャートを露光して得られたパターンの振幅の計算
において、ピーク値の中央値または該中央値±ｎ（ここ
で、ｎは同一周波数での繰り返し回数をＮとしたとき、
ｎ＜Ｎ／２−１で表わされる数である）の平均値、ある
いは、最大，最小値からｍ個を除いた残りのデータの平
均値を用いることを特徴とするものである。
め、本発明に係る鮮鋭度評価方法は、矩形波チャートを
露光して得られたデジタルデータからＣＴＦを計算し、
これに基づいて画像の鮮鋭度の評価を行う際に、前記矩
形波チャートを露光して得られたパターンの振幅の計算
において、ピーク値の中央値または該中央値±ｎ（ここ
で、ｎは同一周波数での繰り返し回数をＮとしたとき、
ｎ＜Ｎ／２−１で表わされる数である）の平均値、ある
いは、最大，最小値からｍ個を除いた残りのデータの平
均値を用いることを特徴とするものである。
【０００９】また、本発明に係る鮮鋭度評価方法では、
前記方法により計算したＣＴＦの値の中央値または該中
央値±ｎ’（ここで、ｎ’は同一周波数での繰り返し回
数をＮ’としたとき、ｎ’＜Ｎ’／２−１で表わされる
数である）の平均値、あるいは、最大，最小値からｍ個
を除いた残りのデータの平均値を採用することを特徴と
するものである。
前記方法により計算したＣＴＦの値の中央値または該中
央値±ｎ’（ここで、ｎ’は同一周波数での繰り返し回
数をＮ’としたとき、ｎ’＜Ｎ’／２−１で表わされる
数である）の平均値、あるいは、最大，最小値からｍ個
を除いた残りのデータの平均値を採用することを特徴と
するものである。
【００１０】さらに、本発明に係る鮮鋭度評価方法にお
いては、前記デジタルデータとしては、対数変換された
ものを用いることを特徴とする。またさらに、本発明に
係る鮮鋭度評価方法においては、算出されたＣＴＦから
Ｃoltmanの換算式を用いてＭＴＦを計算することを特徴
とするものである。
いては、前記デジタルデータとしては、対数変換された
ものを用いることを特徴とする。またさらに、本発明に
係る鮮鋭度評価方法においては、算出されたＣＴＦから
Ｃoltmanの換算式を用いてＭＴＦを計算することを特徴
とするものである。
【００１１】一方、上述の各鮮鋭度評価方法は、請求項
５〜８に示した通り、これを、コンピュータを用いた画
像の鮮鋭度評価装置として、具体化することも可能なも
のである。
５〜８に示した通り、これを、コンピュータを用いた画
像の鮮鋭度評価装置として、具体化することも可能なも
のである。
【００１２】すなわち、本発明は、矩形波チャートを露
光して得られたデジタルデータからＣＴＦを計算し、こ
れに基づいて画像の鮮鋭度の評価を行うものであり、前
記矩形波チャートを露光して得られたパターンの振幅を
計算する際に、ピーク値の中央値または該中央値±ｎ
（ここで、ｎは同一周波数での繰り返し回数をＮとした
とき、ｎ＜Ｎ／２−１）の平均値、あるいは、最大・最
小値からｍ個を除いた残りのデータの平均値を用いる評
価処理手段を有することを特徴とする鮮鋭度評価装置と
して、具体化することができる。
光して得られたデジタルデータからＣＴＦを計算し、こ
れに基づいて画像の鮮鋭度の評価を行うものであり、前
記矩形波チャートを露光して得られたパターンの振幅を
計算する際に、ピーク値の中央値または該中央値±ｎ
（ここで、ｎは同一周波数での繰り返し回数をＮとした
とき、ｎ＜Ｎ／２−１）の平均値、あるいは、最大・最
小値からｍ個を除いた残りのデータの平均値を用いる評
価処理手段を有することを特徴とする鮮鋭度評価装置と
して、具体化することができる。
【００１３】また、本発明に係る画像の鮮鋭度評価装置
は、矩形波チャートを露光して得られたデジタルデータ
に対して、請求項１に記載の方法により計算したＣＴＦ
の値の中央値または該中央値±ｎ’（ここで、ｎ’は同
一周波数での繰り返し回数をＮ’としたとき、ｎ’＜
Ｎ’／２−１で表わされる数である）の平均値、あるい
は、最大，最小値からｍ個を除いた残りのデータの平均
値を採用し、これに基づいて画像の鮮鋭度評価を行う評
価処理手段を有することを特徴とする。
は、矩形波チャートを露光して得られたデジタルデータ
に対して、請求項１に記載の方法により計算したＣＴＦ
の値の中央値または該中央値±ｎ’（ここで、ｎ’は同
一周波数での繰り返し回数をＮ’としたとき、ｎ’＜
Ｎ’／２−１で表わされる数である）の平均値、あるい
は、最大，最小値からｍ個を除いた残りのデータの平均
値を採用し、これに基づいて画像の鮮鋭度評価を行う評
価処理手段を有することを特徴とする。
【００１４】なお、本発明に係る画像の鮮鋭度評価装置
においても、前記デジタルデータとしては、対数変換さ
れたものを用いることが好ましい。またさらに、本発明
に係る画像の鮮鋭度評価装置は、上で算出されたＣＴＦ
からＣoltmanの換算式を用いてＭＴＦを計算する処理手
段を有することを特徴とするものである。
においても、前記デジタルデータとしては、対数変換さ
れたものを用いることが好ましい。またさらに、本発明
に係る画像の鮮鋭度評価装置は、上で算出されたＣＴＦ
からＣoltmanの換算式を用いてＭＴＦを計算する処理手
段を有することを特徴とするものである。
【００１５】上述のように構成したことにより、本発明
によれば、画像の鮮鋭度評価を行う際に、チャート法を
用いるＣＴＦ測定においても、正確なＣＴＦの測定が可
能になるとともに、簡単な操作で、容易かつ確実に、画
像の鮮鋭度評価を行うために好適に用い得る画像の鮮鋭
度評価装置を提供することができる。
によれば、画像の鮮鋭度評価を行う際に、チャート法を
用いるＣＴＦ測定においても、正確なＣＴＦの測定が可
能になるとともに、簡単な操作で、容易かつ確実に、画
像の鮮鋭度評価を行うために好適に用い得る画像の鮮鋭
度評価装置を提供することができる。
【００１６】
【発明の実施の形態】以下、添付の図面に示す好適実施
例に基づいて、本発明の実施の形態を詳細に説明する。
例に基づいて、本発明の実施の形態を詳細に説明する。
【００１７】図１は、本発明の一実施例に係るデジタル
ラジオグラフィー装置の鮮鋭度評価装置（以下、単に評
価装置という）の概略構成を示すブロック図である。本
実施例に係る評価装置は、図示されていないデジタルラ
ジオグラフィー装置において前述の矩形波チャートを露
光した蓄積性蛍光体シート（以下、単に蛍光体シートと
いう）２０を受け入れ、この蛍光体シート２０上の発光
強度データから、ＣＴＦを算出するという装置である。
ラジオグラフィー装置の鮮鋭度評価装置（以下、単に評
価装置という）の概略構成を示すブロック図である。本
実施例に係る評価装置は、図示されていないデジタルラ
ジオグラフィー装置において前述の矩形波チャートを露
光した蓄積性蛍光体シート（以下、単に蛍光体シートと
いう）２０を受け入れ、この蛍光体シート２０上の発光
強度データから、ＣＴＦを算出するという装置である。
【００１８】なお、ここでは、本実施例に係る評価装置
の評価対象となるデジタルラジオグラフィー装置におけ
るサンプリング周波数は１０ｐｉｘ／ｍｍであり、上述
の矩形波チャートとしては、０．０５，０．５，１．
０，１．５，２．０，２．５，３．０，４．０，５．
０，６．０，８．０，１０．０ｃ／ｍｍの空間周波数を
有する矩形波チャート（化成オプトニクス（株）製「Ｔ
ＹＰＥ１」）を用いるものとする。
の評価対象となるデジタルラジオグラフィー装置におけ
るサンプリング周波数は１０ｐｉｘ／ｍｍであり、上述
の矩形波チャートとしては、０．０５，０．５，１．
０，１．５，２．０，２．５，３．０，４．０，５．
０，６．０，８．０，１０．０ｃ／ｍｍの空間周波数を
有する矩形波チャート（化成オプトニクス（株）製「Ｔ
ＹＰＥ１」）を用いるものとする。
【００１９】なお、上述のサンプリング条件において
は、ナイキスト周波数が５ｃ／ｍｍとなることから、こ
の条件下で使用される本実施例に係る評価装置では、上
述の矩形波チャート上のパターンを読み取ったデジタル
データ（以下、単にデータという）からＣＴＦを算出す
る際に、次に示す各空間周波数の信号に基づくデータは
用いないようにしたものである。
は、ナイキスト周波数が５ｃ／ｍｍとなることから、こ
の条件下で使用される本実施例に係る評価装置では、上
述の矩形波チャート上のパターンを読み取ったデジタル
データ（以下、単にデータという）からＣＴＦを算出す
る際に、次に示す各空間周波数の信号に基づくデータは
用いないようにしたものである。
【００２０】すなわち、まず、上述の矩形波チャート上
のパターンのナイキスト周波数（５ｃ／ｍｍ）を越える
周波数に対応する値（６．０，８．０，１０．０ｃ／ｍ
ｍ）を用いないようにする他、ナイキスト周波数の整数
分の１の周波数に該当する５．０および２．５ｃ／ｍｍ
の空間周波数に対応する値（データ）をも用いないよう
にしている。
のパターンのナイキスト周波数（５ｃ／ｍｍ）を越える
周波数に対応する値（６．０，８．０，１０．０ｃ／ｍ
ｍ）を用いないようにする他、ナイキスト周波数の整数
分の１の周波数に該当する５．０および２．５ｃ／ｍｍ
の空間周波数に対応する値（データ）をも用いないよう
にしている。
【００２１】これは、周知のように、矩形波チャートを
用いるＣＴＦの測定では、サンプリング周波数とチャー
ト上のパターンの周波数が整数倍の関係にある場合に、
フェーズがちょうど合ってしまってレスポンスが実際よ
り高く出るとか、あるいは逆に、位相がちょうど半分ず
れている場合に、レスポンスが実際より低く出てしまう
といった現象（in phase, out phase などと呼ばれる）
が起き、測定値が安定しなくなることに対処するためで
ある。
用いるＣＴＦの測定では、サンプリング周波数とチャー
ト上のパターンの周波数が整数倍の関係にある場合に、
フェーズがちょうど合ってしまってレスポンスが実際よ
り高く出るとか、あるいは逆に、位相がちょうど半分ず
れている場合に、レスポンスが実際より低く出てしまう
といった現象（in phase, out phase などと呼ばれる）
が起き、測定値が安定しなくなることに対処するためで
ある。
【００２２】図１に戻って、図中、１１はデータ取得部
であり、上述の蛍光体シート２０上に形成されている矩
形波チャート上の各パターンに対応する部分からの発光
強度を読み取る機能を有するものである。１２はキーボ
ードなどの操作部、１３は表示部、１４は上記データ取
得部１１で取得した発光強度データと、使用した矩形波
チャートのタイプ，固有番号などを記憶するデータ記憶
部である。
であり、上述の蛍光体シート２０上に形成されている矩
形波チャート上の各パターンに対応する部分からの発光
強度を読み取る機能を有するものである。１２はキーボ
ードなどの操作部、１３は表示部、１４は上記データ取
得部１１で取得した発光強度データと、使用した矩形波
チャートのタイプ，固有番号などを記憶するデータ記憶
部である。
【００２３】１５は上述のデータ記憶部１４に記憶され
ている読み取り（発光強度）データとこれに対応する矩
形波チャートのタイプ，固有番号などに基づいて、後述
する手順に従ってＣＴＦを算出する評価処理部を示して
いる。また、１６はこの評価処理部１５で算出されたＣ
ＴＦから、後述するＣoltmanの換算式を用いてＭＴＦを
計算する演算部を示している。
ている読み取り（発光強度）データとこれに対応する矩
形波チャートのタイプ，固有番号などに基づいて、後述
する手順に従ってＣＴＦを算出する評価処理部を示して
いる。また、１６はこの評価処理部１５で算出されたＣ
ＴＦから、後述するＣoltmanの換算式を用いてＭＴＦを
計算する演算部を示している。
【００２４】以下、上述のように構成された本実施例の
動作の概要を説明する。
動作の概要を説明する。
【００２５】デジタルラジオグラフィー装置の評価を行
うためには、まず、当該デジタルラジオグラフィー装置
の蛍光体シート装填部に上述の蛍光体シート２０をセッ
トし、キーボードなどの操作部から指示を行って、前述
の所定の矩形波チャートを露光する。ここでは、前述の
ような周波数分布を有する矩形波チャートを用いるが、
矩形波チャートとしては、評価の都度、所望の矩形波チ
ャートをセットすれば良く、また、予め種々のタイプの
ものをセットしておいてその中から選択するようにして
も良い。
うためには、まず、当該デジタルラジオグラフィー装置
の蛍光体シート装填部に上述の蛍光体シート２０をセッ
トし、キーボードなどの操作部から指示を行って、前述
の所定の矩形波チャートを露光する。ここでは、前述の
ような周波数分布を有する矩形波チャートを用いるが、
矩形波チャートとしては、評価の都度、所望の矩形波チ
ャートをセットすれば良く、また、予め種々のタイプの
ものをセットしておいてその中から選択するようにして
も良い。
【００２６】前述のように、矩形波チャートには複数の
空間周波数に対応するパターンが含まれており、データ
取得部１１では、この矩形波チャートを露光した蛍光体
シート２０から、上述の複数のパターンに対応する発光
強度のデータを取得する。この発光強度のデータは、前
述のデータ記憶部１４に保存される。なお、ここで保存
されるデータは、評価対象となるデジタルラジオグラフ
ィー装置のＸ線照射線量を一定にして露光されたもので
あることが必要である。
空間周波数に対応するパターンが含まれており、データ
取得部１１では、この矩形波チャートを露光した蛍光体
シート２０から、上述の複数のパターンに対応する発光
強度のデータを取得する。この発光強度のデータは、前
述のデータ記憶部１４に保存される。なお、ここで保存
されるデータは、評価対象となるデジタルラジオグラフ
ィー装置のＸ線照射線量を一定にして露光されたもので
あることが必要である。
【００２７】なお、デジタルラジオグラフィー装置にお
いて取り扱うデータは、濃度情報である画像との対応付
けをよくするために、適宜の段階、つまり、読み取り
後、Ａ／Ｄ変換器によりアナログデータからデジタルデ
ータに変換する前、または変換後に、対数変換器により
対数値に変換しておくことが好ましい。
いて取り扱うデータは、濃度情報である画像との対応付
けをよくするために、適宜の段階、つまり、読み取り
後、Ａ／Ｄ変換器によりアナログデータからデジタルデ
ータに変換する前、または変換後に、対数変換器により
対数値に変換しておくことが好ましい。
【００２８】評価処理部１５は、データ記憶部１４に保
存されている蛍光体シート２０からの読み取りデータを
読み込み、ＣＴＦの算出を行う。この際に、矩形波チャ
ート上の空間周波数０．０５に対応するパターンから得
られた値を空間周波数０に相当するものとして、これを
用いて正規化を行う。図２（ａ），（ｂ）により、この
点を詳細に説明する。
存されている蛍光体シート２０からの読み取りデータを
読み込み、ＣＴＦの算出を行う。この際に、矩形波チャ
ート上の空間周波数０．０５に対応するパターンから得
られた値を空間周波数０に相当するものとして、これを
用いて正規化を行う。図２（ａ），（ｂ）により、この
点を詳細に説明する。
【００２９】図２（ａ）は、矩形波チャートのパターン
を示す模式図、図２（ｂ）は、これを読み取ったデータ
を示す模式図である。ここで、図２（ｂ）中のＣ0outは
基準値、Ｃx は任意の空間周波数に対応する値である。
これから下の式（１）に従ってＣＴＦの値Ｒx を求め
る。
を示す模式図、図２（ｂ）は、これを読み取ったデータ
を示す模式図である。ここで、図２（ｂ）中のＣ0outは
基準値、Ｃx は任意の空間周波数に対応する値である。
これから下の式（１）に従ってＣＴＦの値Ｒx を求め
る。
【数１】
【００３０】上述の処理において、評価処理部１５は、
この系のナイキスト周波数５ｃ／ｍｍに関連して、前述
の使用しない周波数のうち、２．５、５．０ｃ／ｍｍの
データは、これをそのまま用いることなく、例えば２．
５ｃ／ｍｍであれば２．０，３．０ｃ／ｍｍの２点のデ
ータを用いて、例えば２次関数近似などにより内挿して
得た値を用いるものとする。
この系のナイキスト周波数５ｃ／ｍｍに関連して、前述
の使用しない周波数のうち、２．５、５．０ｃ／ｍｍの
データは、これをそのまま用いることなく、例えば２．
５ｃ／ｍｍであれば２．０，３．０ｃ／ｍｍの２点のデ
ータを用いて、例えば２次関数近似などにより内挿して
得た値を用いるものとする。
【００３１】また、評価処理部１５は、ＣＴＦを算出す
る際に使用する周波数のデータについて、前述のよう
に、前記矩形波チャートを露光して得られたパターンの
振幅を計算する際、ピーク値の中央値を用いる。すなわ
ち、矩形波チャートを露光して得られたパターンの読み
取りを複数ライン（例えば、１００ライン）にわたって
行い、これらの測定値の単純な平均値ではなく、中央値
を使うようにするというものである。
る際に使用する周波数のデータについて、前述のよう
に、前記矩形波チャートを露光して得られたパターンの
振幅を計算する際、ピーク値の中央値を用いる。すなわ
ち、矩形波チャートを露光して得られたパターンの読み
取りを複数ライン（例えば、１００ライン）にわたって
行い、これらの測定値の単純な平均値ではなく、中央値
を使うようにするというものである。
【００３２】なお、ＣＴＦ算出時の処理としては、上述
のピーク値の中央値を用いる方式の他にも、ピーク値の
中央値±ｎ（ここで、ｎは同一周波数での繰り返し回数
をＮとしたとき、ｎ＜Ｎ／２−１で表わされる数であ
る）の平均値を用いる方式、あるいは、最大，最小値か
らｍ個を除いた残りのデータの平均値を用いる方式も、
好適に利用することが可能である。
のピーク値の中央値を用いる方式の他にも、ピーク値の
中央値±ｎ（ここで、ｎは同一周波数での繰り返し回数
をＮとしたとき、ｎ＜Ｎ／２−１で表わされる数であ
る）の平均値を用いる方式、あるいは、最大，最小値か
らｍ個を除いた残りのデータの平均値を用いる方式も、
好適に利用することが可能である。
【００３３】ＣＴＦを算出する際にこのような配慮を行
うことにより、蛍光体シート２０からのデータ読み取り
時におけるノイズの影響を排除し、高精度にＣＴＦを算
出することが可能になる。
うことにより、蛍光体シート２０からのデータ読み取り
時におけるノイズの影響を排除し、高精度にＣＴＦを算
出することが可能になる。
【００３４】演算部１６では、評価処理部１５によって
算出されたＣＴＦの値に基づいて，前述のＣoltmanの換
算式を用いてＭＴＦを計算する。Ｃoltmanの換算式は、
よく知られているように、
算出されたＣＴＦの値に基づいて，前述のＣoltmanの換
算式を用いてＭＴＦを計算する。Ｃoltmanの換算式は、
よく知られているように、
【数２】 ここで、Ｒsq（ω）はＣＴＦの値、Ｒ（ω）波ＭＴＦの
値を示す。であり、これにより、ＣＴＦとＭＴＦとは、
１対１対応に変換される。
値を示す。であり、これにより、ＣＴＦとＭＴＦとは、
１対１対応に変換される。
【００３５】上記実施例によれば、矩形波チャートを読
み取ったデータについて、測定時のノイズの影響を受け
る可能性のある部分のデータを排除してＣＴＦを算出す
るようにしたので、高精度のＣＴＦの算出が可能になる
という効果が得られる。なお、このようなノイズの影響
を受ける可能性のある部分のデータを排除する方式につ
いては、各種の方式が利用可能であることは、前述の通
りである。
み取ったデータについて、測定時のノイズの影響を受け
る可能性のある部分のデータを排除してＣＴＦを算出す
るようにしたので、高精度のＣＴＦの算出が可能になる
という効果が得られる。なお、このようなノイズの影響
を受ける可能性のある部分のデータを排除する方式につ
いては、各種の方式が利用可能であることは、前述の通
りである。
【００３６】また、演算部１６における、ＣＴＦからＭ
ＴＦを計算する過程においても、以下に説明するよう
に、用いるＣＴＦの値を選択することが有効である。す
なわち、例えば、前述のように、評価処理部１５で算出
されたＣＴＦの値の中央値またはこの中央値±ｎ’（こ
こで、ｎ’は同一周波数での繰り返し回数をＮ’とした
とき、ｎ’＜Ｎ’／２−１で表わされる数である）の平
均値、あるいは、最大，最小値からｍ個を除いた残りの
データの平均値を採用する方式などが好適に用い得る。
ＴＦを計算する過程においても、以下に説明するよう
に、用いるＣＴＦの値を選択することが有効である。す
なわち、例えば、前述のように、評価処理部１５で算出
されたＣＴＦの値の中央値またはこの中央値±ｎ’（こ
こで、ｎ’は同一周波数での繰り返し回数をＮ’とした
とき、ｎ’＜Ｎ’／２−１で表わされる数である）の平
均値、あるいは、最大，最小値からｍ個を除いた残りの
データの平均値を採用する方式などが好適に用い得る。
【００３７】以上説明した実施例においては、基準にな
る矩形波チャートとして所定のものをセットする場合を
例に挙げたが、前述の通り、矩形波チャートは、予め複
数のタイプのものをデジタルラジオグラフィー装置にセ
ットしておいて、この中から適切なものを選択するよう
にしてもよい。この場合、デジタルラジオグラフィー装
置側で、選択された矩形波チャートのタイプ、あるいは
固有番号などを識別して、評価装置側に伝達することが
好ましい。
る矩形波チャートとして所定のものをセットする場合を
例に挙げたが、前述の通り、矩形波チャートは、予め複
数のタイプのものをデジタルラジオグラフィー装置にセ
ットしておいて、この中から適切なものを選択するよう
にしてもよい。この場合、デジタルラジオグラフィー装
置側で、選択された矩形波チャートのタイプ、あるいは
固有番号などを識別して、評価装置側に伝達することが
好ましい。
【００３８】また、上記実施例においては、評価対象と
なるデジタルラジオグラフィー装置とこれを評価する鮮
鋭度評価装置とを別々の装置とした形で説明したが、こ
れらを一体化した、評価装置を内蔵したデジタルラジオ
グラフィー装置とすることも有効な態様である。この他
にも、１台の評価装置と複数のデジタルラジオグラフィ
ー装置とをオンラインで結ぶなど、異なる態様が採用可
能である。
なるデジタルラジオグラフィー装置とこれを評価する鮮
鋭度評価装置とを別々の装置とした形で説明したが、こ
れらを一体化した、評価装置を内蔵したデジタルラジオ
グラフィー装置とすることも有効な態様である。この他
にも、１台の評価装置と複数のデジタルラジオグラフィ
ー装置とをオンラインで結ぶなど、異なる態様が採用可
能である。
【００３９】
【発明の効果】以上、詳細に説明したように、本発明に
よれば、チャート法を用いるＣＴＦ測定においても、高
精度なＣＴＦの測定を可能とした、鮮鋭度評価方法およ
びこの方法を用いる鮮鋭度評価装置を実現することが可
能である。
よれば、チャート法を用いるＣＴＦ測定においても、高
精度なＣＴＦの測定を可能とした、鮮鋭度評価方法およ
びこの方法を用いる鮮鋭度評価装置を実現することが可
能である。
【図１】 本発明の一実施例に係るデジタルラジオグラ
フィー装置の評価装置の概略構成を示すブロック図であ
る。
フィー装置の評価装置の概略構成を示すブロック図であ
る。
【図２】 （ａ）は矩形波チャートのパターンを示す模
式図、（ｂ）はこれを読み取ったデータを示す模式図で
ある。
式図、（ｂ）はこれを読み取ったデータを示す模式図で
ある。
１１ データ取得手段 １２ 操作部 １３ 表示部 １４ データ記憶部 １５ 評価処理部 １６ 演算部 ２０ 蛍光体シート
Claims (8)
- 【請求項１】矩形波チャートを露光して得られたデジタ
ルデータからＣＴＦを計算し、これに基づいて画像の鮮
鋭度の評価を行う鮮鋭度評価方法であって、前記矩形波
チャートを露光して得られたパターンの振幅を計算する
際、ピーク値の中央値または該中央値±ｎ（ここで、ｎ
は同一周波数での繰り返し回数をＮとしたとき、ｎ＜Ｎ
／２−１で表わされる数である）の平均値、あるいは、
最大，最小値からｍ個を除いた残りのデータの平均値を
用いることを特徴とする鮮鋭度評価方法。 - 【請求項２】矩形波チャートを露光して得られたデジタ
ルデータからＣＴＦを計算し、これに基づいて画像の鮮
鋭度の評価を行う鮮鋭度評価方法であって、請求項１に
記載の方法により計算したＣＴＦの値の中央値または該
中央値±ｎ’（ここで、ｎ’は同一周波数での繰り返し
回数をＮ’としたとき、ｎ’＜Ｎ’／２−１で表わされ
る数である）の平均値、あるいは、最大，最小値からｍ
個を除いた残りのデータの平均値を採用することを特徴
とする鮮鋭度評価方法。 - 【請求項３】前記デジタルデータとしては、対数変換さ
れたものを用いる請求項１または２に記載の鮮鋭度評価
方法。 - 【請求項４】前記算出したＣＴＦからＣoltmanの換算式
を用いてＭＴＦを計算する請求項１〜３のいずれか１項
に記載の鮮鋭度評価方法。 - 【請求項５】矩形波チャートを露光して得られたデジタ
ルデータに対して、請求項１記載の計算機能を適用して
ＣＴＦを算出し、これに基づいて画像の鮮鋭度評価を行
う評価処理手段を有することを特徴とする鮮鋭度評価装
置。 - 【請求項６】矩形波チャートを露光して得られたデジタ
ルデータに対して、請求項２記載の計算機能を適用して
ＣＴＦを算出し、これに基づいて画像の鮮鋭度評価を行
う評価処理手段を有することを特徴とする鮮鋭度評価装
置。 - 【請求項７】前記デジタルデータとしては、対数変換さ
れたものを用いる請求項５または６に記載の鮮鋭度評価
装置。 - 【請求項８】前記評価処理手段に加えて、算出したＣＴ
ＦからＣoltmanの換算式を用いてＭＴＦを計算する処理
手段を有する請求項５〜７のいずれか１項に記載の鮮鋭
度評価装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10766499A JP2000298076A (ja) | 1999-04-15 | 1999-04-15 | 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10766499A JP2000298076A (ja) | 1999-04-15 | 1999-04-15 | 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2000298076A true JP2000298076A (ja) | 2000-10-24 |
Family
ID=14464894
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10766499A Withdrawn JP2000298076A (ja) | 1999-04-15 | 1999-04-15 | 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2000298076A (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007024889A (ja) * | 2005-07-12 | 2007-02-01 | Canon Inc | Ｏｔｆ測定システム |
| JP2007158846A (ja) * | 2005-12-06 | 2007-06-21 | Sharp Corp | 画像読取装置および方法 |
| JP2019075697A (ja) * | 2017-10-16 | 2019-05-16 | 日本放送協会 | 解像度特性測定装置およびそのプログラム |
-
1999
- 1999-04-15 JP JP10766499A patent/JP2000298076A/ja not_active Withdrawn
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007024889A (ja) * | 2005-07-12 | 2007-02-01 | Canon Inc | Ｏｔｆ測定システム |
| JP2007158846A (ja) * | 2005-12-06 | 2007-06-21 | Sharp Corp | 画像読取装置および方法 |
| JP2019075697A (ja) * | 2017-10-16 | 2019-05-16 | 日本放送協会 | 解像度特性測定装置およびそのプログラム |
| JP7058102B2 (ja) | 2017-10-16 | 2022-04-21 | 日本放送協会 | 解像度特性測定装置およびそのプログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3188491B2 (ja) | Ｘ線記録のダイナミック圧縮方法及びその装置 | |
| Fujita et al. | A simple method for determining the modulation transfer function in digital radiography | |
| Samei et al. | Intercomparison of methods for image quality characterization. i. modulation transfer function a | |
| EP0032237B1 (en) | Method of and apparatus for processing a radiographic image | |
| JPH0445689Y2 (ja) | ||
| CN101129266B (zh) | X线图像处理系统 | |
| US5550888A (en) | Method of displaying notably a digital X-ray image as a visible image, and device for carrying out the method | |
| US6934590B2 (en) | Quality control system for medical diagnostic apparatus | |
| US9996910B2 (en) | Radiographic image processing device, method, and recording medium | |
| US5297036A (en) | Method for the correction of the measurements of optical density made on a radiographic film | |
| JP4574181B2 (ja) | 画像処理方法及び装置 | |
| JP3459745B2 (ja) | 画像処理装置、放射線撮影装置及び画像処理方法 | |
| Boone | Scatter correction algorithm for digitally acquired radiographs: Theory and results | |
| Schaetzing | Management of pediatric radiation dose using Agfa computed radiography | |
| JP2005005792A (ja) | 画像処理方法および装置並びにｘ線撮影装置 | |
| Mah et al. | Evaluation of a quality control phantom for digital chest radiography | |
| US7796792B2 (en) | Method of identifying disturbing frequencies originating from the presence of an anti-scatter grid during acquisition of a radiation image | |
| JP2008073515A (ja) | Ｘ線画像システム | |
| JPH0614911A (ja) | Ｘ線診断方法及びその装置 | |
| JP2000298076A (ja) | 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 | |
| US9978132B2 (en) | Radiation image processing device, method, and program | |
| JP2003310593A (ja) | 位相情報復元方法及び位相情報復元装置、並びに、位相情報復元プログラム | |
| Higashide et al. | Application of a variable filter for presampled modulation transfer function analysis with the edge method | |
| JP2004275769A (ja) | Ｘ線露光値を求めて記録する方法およびｘ線診断装置 | |
| JP5534699B2 (ja) | Ｘ線診断装置及び画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20060704 |