JP2003141535A - 画像評価方法および装置並びにプログラム - Google Patents
画像評価方法および装置並びにプログラムInfo
- Publication number
- JP2003141535A JP2003141535A JP2001337418A JP2001337418A JP2003141535A JP 2003141535 A JP2003141535 A JP 2003141535A JP 2001337418 A JP2001337418 A JP 2001337418A JP 2001337418 A JP2001337418 A JP 2001337418A JP 2003141535 A JP2003141535 A JP 2003141535A
- Authority
- JP
- Japan
- Prior art keywords
- area
- evaluation
- image
- signal
- target image
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 18
- 238000011156 evaluation Methods 0.000 claims abstract description 270
- 238000001514 detection method Methods 0.000 claims abstract description 73
- 238000012360 testing method Methods 0.000 claims description 10
- 238000007689 inspection Methods 0.000 claims description 8
- 238000012795 verification Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 abstract description 21
- 238000002601 radiography Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 20
- 230000005855 radiation Effects 0.000 description 10
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 7
- 235000015220 hamburgers Nutrition 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 239000000428 dust Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 229920006397 acrylic thermoplastic Polymers 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006996 mental state Effects 0.000 description 1
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- ISXSCDLOGDJUNJ-UHFFFAOYSA-N tert-butyl prop-2-enoate Chemical compound CC(C)(C)OC(=O)C=C ISXSCDLOGDJUNJ-UHFFFAOYSA-N 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/001—Industrial image inspection using an image reference approach
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】
【課題】 画像データを生成する放射線画像撮影装置や
画像処理を行う画像処理装置において得られる画像の評
価を、ＣＤＲＡＤファントム等を用いて効率よく行う。 【解決手段】 ＣＤＲＡＤファントムの放射線画像を撮
影して評価対象画像データＳ０を得、これを用いて画像
の評価を行う。ＣＤＲＡＤファントムと同様のパターン
の画像を表す参照画像データＳ１と評価対象画像データ
Ｓ０との位置合わせを行い、位置合わせ後、評価対象画
像および参照画像から格子線により囲まれる格子領域を
検出する。そして、評価対象画像の全ての格子領域に存
在する信号領域の位置を検出し、正解テーブルと検出結
果を照合して格子領域毎の正誤結果を得、この正誤結果
に基づいて評価対象画像の評価を行う。
画像処理を行う画像処理装置において得られる画像の評
価を、ＣＤＲＡＤファントム等を用いて効率よく行う。 【解決手段】 ＣＤＲＡＤファントムの放射線画像を撮
影して評価対象画像データＳ０を得、これを用いて画像
の評価を行う。ＣＤＲＡＤファントムと同様のパターン
の画像を表す参照画像データＳ１と評価対象画像データ
Ｓ０との位置合わせを行い、位置合わせ後、評価対象画
像および参照画像から格子線により囲まれる格子領域を
検出する。そして、評価対象画像の全ての格子領域に存
在する信号領域の位置を検出し、正解テーブルと検出結
果を照合して格子領域毎の正誤結果を得、この正誤結果
に基づいて評価対象画像の評価を行う。
Description
【０００１】
【発明の属する技術分野】本発明は、画像撮影装置や画
像処理装置等において得られる画像を評価する画像評価
方法および装置並びに画像評価方法をコンピュータに実
行させるためのプログラムに関するものである。
像処理装置等において得られる画像を評価する画像評価
方法および装置並びに画像評価方法をコンピュータに実
行させるためのプログラムに関するものである。
【０００２】
【従来の技術】周波数強調処理、平滑化処理、ノイズ除
去処理のような画像処理アルゴリズムを開発するに当た
り、画像処理の効果をどのように評価するかは、開発者
にとって大きな問題である。
去処理のような画像処理アルゴリズムを開発するに当た
り、画像処理の効果をどのように評価するかは、開発者
にとって大きな問題である。
【０００３】画像の画質を表す指標としてＤＱＥ値（De
tective Quantum Efficiency）が用いられている。ＤＱ
Ｅ値とは、下記式により算出される画質を表す値であっ
て、ＤＱＥ値が大きくなるにしたがって画質は向上す
る。
tective Quantum Efficiency）が用いられている。ＤＱ
Ｅ値とは、下記式により算出される画質を表す値であっ
て、ＤＱＥ値が大きくなるにしたがって画質は向上す
る。
【０００４】ＤＱＥ（ｕ）＝｛(log１０ｅ)２×γ２
×ＭＴＦ２(u)｝／｛ｑ×ＷＳ(u)｝ 但し、ＭＴＦ（Modulation Transfer Function）は、Ｃ
ＴＦチャート（ContrastTransfer Function Chart）を
撮影することにより得られるものであり、各周波数帯域
毎の画像信号の解像度の大きさ、すなわち鮮鋭性を表す
指標である。また、ＷＳ(u)はウィナースペクトルであ
り、粒状性を表す指標である。なお、ｕは空間周波数、
γは画像を得る際に用いたフイルムの階調特性を示すγ
値、ｑは量子数を表す。
×ＭＴＦ２(u)｝／｛ｑ×ＷＳ(u)｝ 但し、ＭＴＦ（Modulation Transfer Function）は、Ｃ
ＴＦチャート（ContrastTransfer Function Chart）を
撮影することにより得られるものであり、各周波数帯域
毎の画像信号の解像度の大きさ、すなわち鮮鋭性を表す
指標である。また、ＷＳ(u)はウィナースペクトルであ
り、粒状性を表す指標である。なお、ｕは空間周波数、
γは画像を得る際に用いたフイルムの階調特性を示すγ
値、ｑは量子数を表す。
【０００５】ここで、線形的な画像処理の場合、鮮鋭性
と粒状性とが同じように改善されるため、ＤＱＥ自体は
変化しない。一方、非線形な画像処理の場合、鮮鋭性と
粒状性とを独立して処理できるため、鮮鋭性および粒状
性についての画像処理の有無により、ＤＱＥを大きく改
善することができる。
と粒状性とが同じように改善されるため、ＤＱＥ自体は
変化しない。一方、非線形な画像処理の場合、鮮鋭性と
粒状性とを独立して処理できるため、鮮鋭性および粒状
性についての画像処理の有無により、ＤＱＥを大きく改
善することができる。
【０００６】一方、放射線画像の分野での画像の画質を
評価する代表的な方法として、バーガーファントムを用
いた視覚評価方法が知られている。バーガーファントム
とは、直径および厚さが段階的に変化する複数の円柱状
のアクリルを並べて配列したファントムである。このバ
ーガーファントムを撮影することにより得られる画像を
観察することにより、システムにおいて得られる画像の
空間分解能およびコントラスト分解能の双方を評価する
ことができる。なお、円柱状アクリルの径が空間分解能
を評価する指標であり、その厚さがコントラスト分解能
を評価する指標となっている。ここで、画質の良好な画
像を得ることができるシステムにおいて得られる画像ほ
ど、より径が小さくかつより厚さが小さい円柱状アクリ
ルの像を視認することができる。したがって、画像の評
価を行うシステムにおいて、バーガーファントムの画像
を得、その画像に含まれる円柱状アクリルに対応する部
分が視認できる限界を認識することにより、そのシステ
ムにおいて得られる画像の評価を行うことができる。
評価する代表的な方法として、バーガーファントムを用
いた視覚評価方法が知られている。バーガーファントム
とは、直径および厚さが段階的に変化する複数の円柱状
のアクリルを並べて配列したファントムである。このバ
ーガーファントムを撮影することにより得られる画像を
観察することにより、システムにおいて得られる画像の
空間分解能およびコントラスト分解能の双方を評価する
ことができる。なお、円柱状アクリルの径が空間分解能
を評価する指標であり、その厚さがコントラスト分解能
を評価する指標となっている。ここで、画質の良好な画
像を得ることができるシステムにおいて得られる画像ほ
ど、より径が小さくかつより厚さが小さい円柱状アクリ
ルの像を視認することができる。したがって、画像の評
価を行うシステムにおいて、バーガーファントムの画像
を得、その画像に含まれる円柱状アクリルに対応する部
分が視認できる限界を認識することにより、そのシステ
ムにおいて得られる画像の評価を行うことができる。
【０００７】しかしながら、バーガーファントムは、予
め決められた場所に円柱状アクリルが存在する。したが
って、観察者はその場所に円柱状アクリルの像が存在す
るということを認識しているため、円柱状アクリルの像
が見えないのに見えてしまうという観察者の心理的な影
響により評価結果が左右されるという問題がある。
め決められた場所に円柱状アクリルが存在する。したが
って、観察者はその場所に円柱状アクリルの像が存在す
るということを認識しているため、円柱状アクリルの像
が見えないのに見えてしまうという観察者の心理的な影
響により評価結果が左右されるという問題がある。
【０００８】このため、バーガーファントムに代えてＣ
ＤＲＡＤファントム（Contrast Detail Digital/Conven
tional Radiography Phantom, Northwest X-ray Inc.
製）を用いた評価が行われている。
ＤＲＡＤファントム（Contrast Detail Digital/Conven
tional Radiography Phantom, Northwest X-ray Inc.
製）を用いた評価が行われている。
【０００９】ＣＤＲＡＤファントムは、１５×１５の格
子状に領域が分割されたアクリル板からなり、各格子領
域内に直径および深さを段階的に変更した穴を形成して
なるものである。具体的には、上下方向に直径が、左右
方向に深さが変化するように穴が形成されている。ま
た、上段の３行目までは１つの格子領域に１つの穴が形
成されているが、それ以外の領域は、格子領域の中心お
よび格子領域の４隅のいずれか１箇所に中心と同一の直
径および深さを有する穴が形成されている。そして、こ
のＣＤＲＡＤファントムの放射線画像を取得してＣＤＲ
ＡＤファントム画像を得、このＣＤＲＡＤファントム画
像の評価を行うことにより、ＣＤＲＡＤファントムの放
射線画像を得たシステムにおいて得られる画像の評価を
行うことができる。また、ＣＤＲＡＤファントム画像に
対して画像処理を行って処理済みのＣＤＲＡＤファント
ム画像を得ることにより、画像処理を行うシステムにお
いて得られる画像の評価も行うことができる。以下、画
像の評価について説明する。
子状に領域が分割されたアクリル板からなり、各格子領
域内に直径および深さを段階的に変更した穴を形成して
なるものである。具体的には、上下方向に直径が、左右
方向に深さが変化するように穴が形成されている。ま
た、上段の３行目までは１つの格子領域に１つの穴が形
成されているが、それ以外の領域は、格子領域の中心お
よび格子領域の４隅のいずれか１箇所に中心と同一の直
径および深さを有する穴が形成されている。そして、こ
のＣＤＲＡＤファントムの放射線画像を取得してＣＤＲ
ＡＤファントム画像を得、このＣＤＲＡＤファントム画
像の評価を行うことにより、ＣＤＲＡＤファントムの放
射線画像を得たシステムにおいて得られる画像の評価を
行うことができる。また、ＣＤＲＡＤファントム画像に
対して画像処理を行って処理済みのＣＤＲＡＤファント
ム画像を得ることにより、画像処理を行うシステムにお
いて得られる画像の評価も行うことができる。以下、画
像の評価について説明する。
【００１０】図２は、ＣＤＲＡＤファントムの放射線画
像（以下ＣＤＲＡＤファントム画像とする）を示す図で
ある。穴の直径は下側ほど小さく、穴の深さは左側ほど
小さくなっているため、図２に示すように左下の格子領
域ほど穴の像が見えにくくなっている。なお、ＣＤＲＡ
Ｄファントム画像において格子領域の左側の数字は穴の
直径を、下側の数字は穴の深さを示す。
像（以下ＣＤＲＡＤファントム画像とする）を示す図で
ある。穴の直径は下側ほど小さく、穴の深さは左側ほど
小さくなっているため、図２に示すように左下の格子領
域ほど穴の像が見えにくくなっている。なお、ＣＤＲＡ
Ｄファントム画像において格子領域の左側の数字は穴の
直径を、下側の数字は穴の深さを示す。
【００１１】まず、上側３段の格子領域については、穴
の像が見えなくなった限界を指示し、下側１２段の格子
領域については、穴の位置を指示する。なお、穴は格子
領域の４隅のいずれか１箇所に形成されているため、右
上、右下、左上および左下というように指示を行う。
の像が見えなくなった限界を指示し、下側１２段の格子
領域については、穴の位置を指示する。なお、穴は格子
領域の４隅のいずれか１箇所に形成されているため、右
上、右下、左上および左下というように指示を行う。
【００１２】指示の後、全ての格子領域（１５×１５＝
２２５）について、指示した結果と正解位置とを比較
し、格子領域毎に正誤結果を求める。
２２５）について、指示した結果と正解位置とを比較
し、格子領域毎に正誤結果を求める。
【００１３】なお、１／４の確率により偶然に正解とな
ってしまう場合がある。これを排除するために、以下の
（１）から（４）のルールにより正誤結果を補正する。
ってしまう場合がある。これを排除するために、以下の
（１）から（４）のルールにより正誤結果を補正する。
【００１４】（１）正誤を判定する格子領域（以下注目
格子領域とする）が正解の場合、上下左右の最近傍の４
つの格子領域のうち２以上の格子領域が正解であれば、
その注目領域も正解とする。なお、注目格子領域がＣＤ
ＲＡＤファントム画像の辺に接している場合には、最近
傍の３つの格子領域のうち２以上の格子領域が正解であ
れば、その注目領域も正解とする。
格子領域とする）が正解の場合、上下左右の最近傍の４
つの格子領域のうち２以上の格子領域が正解であれば、
その注目領域も正解とする。なお、注目格子領域がＣＤ
ＲＡＤファントム画像の辺に接している場合には、最近
傍の３つの格子領域のうち２以上の格子領域が正解であ
れば、その注目領域も正解とする。
【００１５】（２）最近傍の格子領域が２つしかない、
すなわち注目格子領域がＣＤＲＡＤファントム画像の４
隅にある場合については、注目格子領域が正解の場合、
最近傍の格子領域の１つが正解であれば、注目格子領域
も正解とする。
すなわち注目格子領域がＣＤＲＡＤファントム画像の４
隅にある場合については、注目格子領域が正解の場合、
最近傍の格子領域の１つが正解であれば、注目格子領域
も正解とする。
【００１６】（３）注目格子領域が不正解の場合、最近
傍の４つまたは３つの格子領域のうち３以上の格子領域
が正解であれば、注目格子領域を正解とする。
傍の４つまたは３つの格子領域のうち３以上の格子領域
が正解であれば、注目格子領域を正解とする。
【００１７】（４）注目格子領域がＣＤＲＡＤファント
ム画像の４隅にある場合において、注目格子領域が不正
解の場合、最近傍の格子領域の全てが正解であれば、注
目格子領域も正解とする。
ム画像の４隅にある場合において、注目格子領域が不正
解の場合、最近傍の格子領域の全てが正解であれば、注
目格子領域も正解とする。
【００１８】そして、補正後の正誤結果を参照し、視認
可能であった穴の直径および穴の深さの限界値を、穴の
深さの対数値を横軸に、穴の直径の対数値を縦軸とした
両対数グラフにプロットすることにより、ＣＤ曲線を得
る。ここで、ＣはＣｏｎｔｒａｓｔを、ＤはＤｅｔａｉ
ｌを表す。ＣＤ曲線においては、より小さい直径かつよ
り深さが小さい穴が見えるほど、ＣＤＲＡＤファントム
画像を得たシステムが良好なものであるということが分
かる。
可能であった穴の直径および穴の深さの限界値を、穴の
深さの対数値を横軸に、穴の直径の対数値を縦軸とした
両対数グラフにプロットすることにより、ＣＤ曲線を得
る。ここで、ＣはＣｏｎｔｒａｓｔを、ＤはＤｅｔａｉ
ｌを表す。ＣＤ曲線においては、より小さい直径かつよ
り深さが小さい穴が見えるほど、ＣＤＲＡＤファントム
画像を得たシステムが良好なものであるということが分
かる。
【００１９】図２に示すように、ＣＤＲＡＤファントム
においては、格子領域の４隅のいずれかに存在する穴の
位置は予め分かっていないため、観察者の主観が入りに
くくなり、その結果、バーガーファントムを用いる場合
よりも客観的な評価結果を得ることができる。
においては、格子領域の４隅のいずれかに存在する穴の
位置は予め分かっていないため、観察者の主観が入りに
くくなり、その結果、バーガーファントムを用いる場合
よりも客観的な評価結果を得ることができる。
【００２０】
【発明が解決しようとする課題】上述したように、ＣＤ
ＲＡＤファントムを用いた画像の評価は、１つの画像に
ついて２２５の格子領域の正誤を判定して正誤結果を
得、さらにその正誤結果を補正してＣＤ曲線を作成する
という手順が必要であり、１画像当たりの評価時間は１
５分程度を要する。また、ＣＤＲＡＤファントムを用い
た画像の評価においても、バーガーファントムを用いた
画像の評価と同様に、評価を行う観察者の主観が入る場
合があるため、観察者による統計的な差異を少なくする
ために、３人以上の観察者が必要となる。また、撮影に
よる条件の差異を統計的に少なくするために、３以上の
ＣＤＲＡＤファントム画像を撮影して評価を行う必要が
ある。さらに、画像処理装置において画像データに施す
画像処理のパラメータを変化させて画像の評価を行う場
合には、変化させたパラメータ毎に３以上のＣＤＲＡＤ
ファントム画像を用意する必要がある。したがって、Ｃ
ＤＲＡＤファントム画像を用いて画像の評価を行うには
非常に手間がかかり、評価のための人件費等のコストが
増大するという問題がある。
ＲＡＤファントムを用いた画像の評価は、１つの画像に
ついて２２５の格子領域の正誤を判定して正誤結果を
得、さらにその正誤結果を補正してＣＤ曲線を作成する
という手順が必要であり、１画像当たりの評価時間は１
５分程度を要する。また、ＣＤＲＡＤファントムを用い
た画像の評価においても、バーガーファントムを用いた
画像の評価と同様に、評価を行う観察者の主観が入る場
合があるため、観察者による統計的な差異を少なくする
ために、３人以上の観察者が必要となる。また、撮影に
よる条件の差異を統計的に少なくするために、３以上の
ＣＤＲＡＤファントム画像を撮影して評価を行う必要が
ある。さらに、画像処理装置において画像データに施す
画像処理のパラメータを変化させて画像の評価を行う場
合には、変化させたパラメータ毎に３以上のＣＤＲＡＤ
ファントム画像を用意する必要がある。したがって、Ｃ
ＤＲＡＤファントム画像を用いて画像の評価を行うには
非常に手間がかかり、評価のための人件費等のコストが
増大するという問題がある。
【００２１】さらに、ＣＤＲＡＤファントム画像を用い
た画像の評価は、人間による視覚評価であるため、ＣＤ
ＲＡＤファントム画像を観察する環境、すなわち、評価
を行う場所の明るさ、評価を行った時刻、評価を行った
順番、観察者の体調や精神状態等に依存して評価結果が
変化する場合があり、安定した評価結果を得ることがで
きないという問題もある。
た画像の評価は、人間による視覚評価であるため、ＣＤ
ＲＡＤファントム画像を観察する環境、すなわち、評価
を行う場所の明るさ、評価を行った時刻、評価を行った
順番、観察者の体調や精神状態等に依存して評価結果が
変化する場合があり、安定した評価結果を得ることがで
きないという問題もある。
【００２２】本発明は上記事情に鑑みなされたものであ
り、ＣＤＲＡＤファントム等を用いた画像の評価を効率
よく行うことを目的とする。
り、ＣＤＲＡＤファントム等を用いた画像の評価を効率
よく行うことを目的とする。
【００２３】
【課題を解決するための手段】本発明による画像評価方
法は、少なくとも１つの信号領域が所定位置に配設され
た複数の評価領域からなるパターンであって、前記信号
領域のサイズおよびコントラストが該各評価領域毎に段
階的に異なるパターンについての評価対象画像を表す評
価対象画像データの入力を受け付け、前記パターンにつ
いての参照画像を表す参照画像データおよび前記評価対
象画像データに基づいて、前記評価対象画像における全
ての前記評価領域に存在する前記信号領域の位置を検出
し、該検出結果に基づいて、前記評価対象画像の評価を
行うことを特徴とするものである。
法は、少なくとも１つの信号領域が所定位置に配設され
た複数の評価領域からなるパターンであって、前記信号
領域のサイズおよびコントラストが該各評価領域毎に段
階的に異なるパターンについての評価対象画像を表す評
価対象画像データの入力を受け付け、前記パターンにつ
いての参照画像を表す参照画像データおよび前記評価対
象画像データに基づいて、前記評価対象画像における全
ての前記評価領域に存在する前記信号領域の位置を検出
し、該検出結果に基づいて、前記評価対象画像の評価を
行うことを特徴とするものである。
【００２４】「所定位置」とは、ＣＤＲＡＤファントム
画像のように、評価領域としての格子領域の中心位置お
よび４隅のいずれか１箇所の位置のように予め定められ
た位置のみならず、少なくとも１つの信号領域を評価領
域内のランダムな位置に配設したものであってもよい。
とくに、信号領域が１つのみしか存在しない場合には、
ランダムな位置に配設することが好ましい。
画像のように、評価領域としての格子領域の中心位置お
よび４隅のいずれか１箇所の位置のように予め定められ
た位置のみならず、少なくとも１つの信号領域を評価領
域内のランダムな位置に配設したものであってもよい。
とくに、信号領域が１つのみしか存在しない場合には、
ランダムな位置に配設することが好ましい。
【００２５】「コントラスト」とは、評価領域内におけ
る信号領域と信号領域以外の領域とにおける濃度の差異
を表すものであり、濃度差が大きいほどコントラストは
大きく、濃度差が小さいほどコントラストは小さいもの
となる。
る信号領域と信号領域以外の領域とにおける濃度の差異
を表すものであり、濃度差が大きいほどコントラストは
大きく、濃度差が小さいほどコントラストは小さいもの
となる。
【００２６】「参照画像データ」とは、評価対象画像と
同様のパターンを表す画像データであるが、評価対象画
像データよりも、より小さいサイズおよび小さいコント
ラストの信号領域まで認識可能な信号値を有する画像デ
ータをいう。したがって、参照画像データを再生するこ
とにより得られる参照画像においては、評価対象画像よ
りも小さいサイズおよび小さいコントラストの信号領域
まで認識することができるものである。もちろん、全て
のサイズおよびコントラストの信号領域を認識できるも
のであってもよい。
同様のパターンを表す画像データであるが、評価対象画
像データよりも、より小さいサイズおよび小さいコント
ラストの信号領域まで認識可能な信号値を有する画像デ
ータをいう。したがって、参照画像データを再生するこ
とにより得られる参照画像においては、評価対象画像よ
りも小さいサイズおよび小さいコントラストの信号領域
まで認識することができるものである。もちろん、全て
のサイズおよびコントラストの信号領域を認識できるも
のであってもよい。
【００２７】なお、本発明による画像評価方法において
は、前記参照画像および前記評価対象画像の位置合わせ
を行った後に、前記信号領域の位置を検出することが好
ましい。
は、前記参照画像および前記評価対象画像の位置合わせ
を行った後に、前記信号領域の位置を検出することが好
ましい。
【００２８】また、本発明による画像評価方法において
は、前記評価対象画像および前記参照画像のそれぞれの
前記評価領域を分割して、前記評価対象画像および前記
参照画像のそれぞれの前記評価領域において相対応する
複数の分割領域を得、前記評価対象画像および前記参照
画像の相対応する前記分割領域内の信号値に基づいて、
前記信号領域の位置を検出するようにしてもよい。
は、前記評価対象画像および前記参照画像のそれぞれの
前記評価領域を分割して、前記評価対象画像および前記
参照画像のそれぞれの前記評価領域において相対応する
複数の分割領域を得、前記評価対象画像および前記参照
画像の相対応する前記分割領域内の信号値に基づいて、
前記信号領域の位置を検出するようにしてもよい。
【００２９】この場合、前記信号領域のサイズに応じ
て，前記分割領域のサイズを変更することが好ましい。
て，前記分割領域のサイズを変更することが好ましい。
【００３０】さらにこの場合、前記評価対象画像および
前記参照画像の相対応する前記分割領域毎の複数の相関
値を算出し、該複数の相関値のうち最大相関値を得た前
記分割領域の位置を前記信号領域の位置として検出する
ことが好ましい。
前記参照画像の相対応する前記分割領域毎の複数の相関
値を算出し、該複数の相関値のうち最大相関値を得た前
記分割領域の位置を前記信号領域の位置として検出する
ことが好ましい。
【００３１】また、本発明による画像評価方法において
は、前記評価対象領域画像の前記評価領域を１つの領域
に１つの前記信号領域が含まれるように分割して略同一
サイズの複数の分割領域を得、前記参照画像の前記評価
領域から、該評価領域に対応する前記評価対象画像の前
記評価領域における前記分割領域に対応するサイズを有
し、１つの前記信号領域を含む検定領域を得、前記評価
対象画像および前記参照画像の相対応する前記評価領域
における前記複数の分割領域および前記検定領域の信号
値に基づいて、前記信号領域の位置を検出するようにし
てもよい。
は、前記評価対象領域画像の前記評価領域を１つの領域
に１つの前記信号領域が含まれるように分割して略同一
サイズの複数の分割領域を得、前記参照画像の前記評価
領域から、該評価領域に対応する前記評価対象画像の前
記評価領域における前記分割領域に対応するサイズを有
し、１つの前記信号領域を含む検定領域を得、前記評価
対象画像および前記参照画像の相対応する前記評価領域
における前記複数の分割領域および前記検定領域の信号
値に基づいて、前記信号領域の位置を検出するようにし
てもよい。
【００３２】この場合、前記信号領域のサイズに応じ
て，前記分割領域のサイズおよび前記検定領域のサイズ
を変更することが好ましい。
て，前記分割領域のサイズおよび前記検定領域のサイズ
を変更することが好ましい。
【００３３】また、この場合、前記評価対象画像および
前記参照画像の相対応する前記評価領域における前記複
数の分割領域および前記検定領域との複数の相関値を算
出し、該複数の相関値のうち最大相関値を得た前記分割
領域の位置を前記信号領域の位置として検出するように
してもよい。
前記参照画像の相対応する前記評価領域における前記複
数の分割領域および前記検定領域との複数の相関値を算
出し、該複数の相関値のうち最大相関値を得た前記分割
領域の位置を前記信号領域の位置として検出するように
してもよい。
【００３４】本発明による画像評価装置は、少なくとも
１つの信号領域が所定位置に配設された複数の評価領域
からなるパターンであって、前記信号領域のサイズおよ
びコントラストが該各評価領域毎に段階的に異なるパタ
ーンについての評価対象画像を表す評価対象画像データ
の入力を受け付ける入力手段と、前記パターンについて
の参照画像を表す参照画像データおよび前記評価対象画
像データに基づいて、前記評価対象画像における全ての
前記評価領域に存在する前記信号領域の位置を検出する
検出手段と、該検出結果に基づいて、前記評価対象画像
の評価を行う評価手段とを備えたことを特徴とするもの
である。
１つの信号領域が所定位置に配設された複数の評価領域
からなるパターンであって、前記信号領域のサイズおよ
びコントラストが該各評価領域毎に段階的に異なるパタ
ーンについての評価対象画像を表す評価対象画像データ
の入力を受け付ける入力手段と、前記パターンについて
の参照画像を表す参照画像データおよび前記評価対象画
像データに基づいて、前記評価対象画像における全ての
前記評価領域に存在する前記信号領域の位置を検出する
検出手段と、該検出結果に基づいて、前記評価対象画像
の評価を行う評価手段とを備えたことを特徴とするもの
である。
【００３５】なお、本発明による画像評価装置において
は、前記参照画像および前記評価対象画像の位置合わせ
を行う位置合わせ手段をさらに備えるものとし、前記検
出手段を、該位置合わせ手段による位置合わせを行った
後に、前記信号領域の位置を検出する手段とすることが
好ましい。
は、前記参照画像および前記評価対象画像の位置合わせ
を行う位置合わせ手段をさらに備えるものとし、前記検
出手段を、該位置合わせ手段による位置合わせを行った
後に、前記信号領域の位置を検出する手段とすることが
好ましい。
【００３６】また、本発明による画像評価装置において
は、前記検出手段を、前記評価対象画像および前記参照
画像のそれぞれの前記評価領域を分割して、前記評価対
象画像および前記参照画像のそれぞれの前記評価領域に
おいて相対応する複数の分割領域を得、前記評価対象画
像および前記参照画像の相対応する前記分割領域内の信
号値に基づいて、前記信号領域の位置を検出する手段と
してもよい。
は、前記検出手段を、前記評価対象画像および前記参照
画像のそれぞれの前記評価領域を分割して、前記評価対
象画像および前記参照画像のそれぞれの前記評価領域に
おいて相対応する複数の分割領域を得、前記評価対象画
像および前記参照画像の相対応する前記分割領域内の信
号値に基づいて、前記信号領域の位置を検出する手段と
してもよい。
【００３７】この場合、前記検出手段を、前記信号領域
のサイズに応じて，前記分割領域のサイズを変更する手
段としてもよい。
のサイズに応じて，前記分割領域のサイズを変更する手
段としてもよい。
【００３８】またこの場合、前記検出手段を、前記評価
対象画像および前記参照画像の相対応する前記分割領域
毎の複数の相関値を算出し、該複数の相関値のうち最大
相関値を得た前記分割領域の位置を前記信号領域の位置
として検出する手段としてもよい。
対象画像および前記参照画像の相対応する前記分割領域
毎の複数の相関値を算出し、該複数の相関値のうち最大
相関値を得た前記分割領域の位置を前記信号領域の位置
として検出する手段としてもよい。
【００３９】さらに、本発明による画像評価装置におい
ては、前記検出手段を、前記評価対象領域画像の前記評
価領域を１つの領域に１つの前記信号領域が含まれるよ
うに分割して略同一サイズの複数の分割領域を得、前記
参照画像の前記評価領域から、該評価領域に対応する前
記評価対象画像の前記評価領域における前記分割領域に
対応するサイズを有し、１つの前記信号領域を含む検定
領域を得、前記評価対象画像および前記参照画像の相対
応する前記評価領域における前記複数の分割領域および
前記検定領域の信号値に基づいて、前記信号領域の位置
を検出する手段としてもよい。
ては、前記検出手段を、前記評価対象領域画像の前記評
価領域を１つの領域に１つの前記信号領域が含まれるよ
うに分割して略同一サイズの複数の分割領域を得、前記
参照画像の前記評価領域から、該評価領域に対応する前
記評価対象画像の前記評価領域における前記分割領域に
対応するサイズを有し、１つの前記信号領域を含む検定
領域を得、前記評価対象画像および前記参照画像の相対
応する前記評価領域における前記複数の分割領域および
前記検定領域の信号値に基づいて、前記信号領域の位置
を検出する手段としてもよい。
【００４０】この場合、前記検出手段を、前記信号領域
のサイズに応じて，前記分割領域のサイズおよび前記検
定領域のサイズを変更する手段としてもよい。
のサイズに応じて，前記分割領域のサイズおよび前記検
定領域のサイズを変更する手段としてもよい。
【００４１】またこの場合、前記検出手段を、前記評価
対象画像および前記参照画像の相対応する前記評価領域
における前記複数の分割領域および前記検定領域との複
数の相関値を算出し、該複数の相関値のうち最大相関値
を得た前記分割領域の位置を前記信号領域の位置として
検出する手段としてもよい。
対象画像および前記参照画像の相対応する前記評価領域
における前記複数の分割領域および前記検定領域との複
数の相関値を算出し、該複数の相関値のうち最大相関値
を得た前記分割領域の位置を前記信号領域の位置として
検出する手段としてもよい。
【００４２】なお、本発明による画像評価方法をコンピ
ュータに実行させるためのプログラムとして提供しても
よい。
ュータに実行させるためのプログラムとして提供しても
よい。
【００４３】
【発明の効果】本発明によれば、評価対象画像データの
入力が受け付けられ、参照画像データおよび評価対象画
像データに基づいて、評価対象画像における全ての評価
領域に存在する信号領域の位置が検出され、検出結果に
基づいて評価対象画像の評価が行われる。このため、評
価対象画像の評価を行う観察者の主観的な要素や評価環
境が評価結果に反映されることがなくなり、その結果、
より客観的かつ安定した評価結果を得ることができる。
また、評価対象画像データの入力を行うのみで客観性の
高い評価結果を得ることができるため、多数の評価対象
画像を多数の評価者によって評価する必要がなくなり、
その結果、得られる画像の画質の評価を効率よく行うこ
とができる。
入力が受け付けられ、参照画像データおよび評価対象画
像データに基づいて、評価対象画像における全ての評価
領域に存在する信号領域の位置が検出され、検出結果に
基づいて評価対象画像の評価が行われる。このため、評
価対象画像の評価を行う観察者の主観的な要素や評価環
境が評価結果に反映されることがなくなり、その結果、
より客観的かつ安定した評価結果を得ることができる。
また、評価対象画像データの入力を行うのみで客観性の
高い評価結果を得ることができるため、多数の評価対象
画像を多数の評価者によって評価する必要がなくなり、
その結果、得られる画像の画質の評価を効率よく行うこ
とができる。
【００４４】また、参照画像と評価対象画像との位置合
わせを行った後に、信号領域の位置を検出することによ
り、信号領域の位置検出の精度を向上させることができ
る。
わせを行った後に、信号領域の位置を検出することによ
り、信号領域の位置検出の精度を向上させることができ
る。
【００４５】さらに、評価対象画像および参照画像のそ
れぞれの評価領域を分割して複数の分割領域を得、評価
対象画像および参照画像の相対応する分割領域内の信号
値に基づいて信号領域の位置を検出することにより、よ
り簡易に信号領域の位置検出を行うことができる。
れぞれの評価領域を分割して複数の分割領域を得、評価
対象画像および参照画像の相対応する分割領域内の信号
値に基づいて信号領域の位置を検出することにより、よ
り簡易に信号領域の位置検出を行うことができる。
【００４６】この場合、信号領域のサイズに応じて分割
領域のサイズを変更することにより、信号領域検出の精
度をより向上させることができる。
領域のサイズを変更することにより、信号領域検出の精
度をより向上させることができる。
【００４７】またこの場合、評価対象画像および参照画
像の相対応する分割領域毎の複数の相関値を算出し、複
数の相関値のうち最大相関値を得た分割領域の位置を信
号領域の位置として検出することにより、より簡易な演
算により信号領域の位置検出を行うことができる。
像の相対応する分割領域毎の複数の相関値を算出し、複
数の相関値のうち最大相関値を得た分割領域の位置を信
号領域の位置として検出することにより、より簡易な演
算により信号領域の位置検出を行うことができる。
【００４８】また、評価対象領域画像の評価領域を１つ
の領域に１つの信号領域が含まれるように分割して略同
一サイズの複数の分割領域を得、参照画像の評価領域か
ら分割領域に対応するサイズを有し、１つの信号領域を
含む検定領域を生成し、評価対象画像および参照画像の
相対応する評価領域における複数の分割領域および検定
領域の信号値に基づいて信号領域の位置を検出すること
により、より簡易に信号領域の位置検出を行うことがで
きる。
の領域に１つの信号領域が含まれるように分割して略同
一サイズの複数の分割領域を得、参照画像の評価領域か
ら分割領域に対応するサイズを有し、１つの信号領域を
含む検定領域を生成し、評価対象画像および参照画像の
相対応する評価領域における複数の分割領域および検定
領域の信号値に基づいて信号領域の位置を検出すること
により、より簡易に信号領域の位置検出を行うことがで
きる。
【００４９】この場合、信号領域のサイズに応じて分割
領域および検定領域のサイズを変更することにより、信
号領域検出の精度をより向上させることができる。
領域および検定領域のサイズを変更することにより、信
号領域検出の精度をより向上させることができる。
【００５０】またこの場合、評価対象画像および参照画
像の相対応する評価領域における複数の分割領域および
検定領域との複数の相関値を算出し、複数の相関値のう
ち最大相関値を得た分割領域の位置を信号領域の位置と
して検出することにより、より簡易な演算により信号領
域の位置検出を行うことができる。
像の相対応する評価領域における複数の分割領域および
検定領域との複数の相関値を算出し、複数の相関値のう
ち最大相関値を得た分割領域の位置を信号領域の位置と
して検出することにより、より簡易な演算により信号領
域の位置検出を行うことができる。
【００５１】
【発明の実施の形態】以下図面を参照して本発明の実施
形態について説明する。図１は本発明の実施形態による
画像評価装置の構成を示す概略ブロック図である。な
お、本実施形態による画像評価装置は、被写体を撮影し
てこの被写体の放射線画像を表す画像データを得る放射
線画像撮影装置あるいは放射線画像撮影装置において得
られた画像データに対して画像処理を施して処理済み画
像データを得る画像処理装置において得られる画像の評
価を行うためのものである。また、本実施形態において
は、ＣＤＲＡＤファントムを撮影することにより得られ
たＣＤＲＡＤファントム画像を表す画像データを評価の
対象とするものである。
形態について説明する。図１は本発明の実施形態による
画像評価装置の構成を示す概略ブロック図である。な
お、本実施形態による画像評価装置は、被写体を撮影し
てこの被写体の放射線画像を表す画像データを得る放射
線画像撮影装置あるいは放射線画像撮影装置において得
られた画像データに対して画像処理を施して処理済み画
像データを得る画像処理装置において得られる画像の評
価を行うためのものである。また、本実施形態において
は、ＣＤＲＡＤファントムを撮影することにより得られ
たＣＤＲＡＤファントム画像を表す画像データを評価の
対象とするものである。
【００５２】図１に示すように、本実施形態による画像
評価装置は、評価の対象となる評価対象画像を表す評価
対象画像データＳ０の入力を受け付ける入力手段１、評
価を行うための参照画像を表す参照画像データＳ１を記
憶した記憶手段２、評価対象画像と参照画像との位置合
わせを行う位置合わせ手段３、位置合わせ後の評価対象
画像および参照画像について、ＣＤＲＡＤファントムの
格子線により区画される複数の格子領域を検出する格子
領域検出手段４、各格子領域について信号領域の位置を
検出して検出結果を取得する検出手段５、検出手段５に
おける信号領域位置の検出結果から正誤結果を求めて出
力する正誤結果出力手段６、および正誤結果に基づいて
評価結果を生成してこれを出力する評価結果出力手段７
を備える。
評価装置は、評価の対象となる評価対象画像を表す評価
対象画像データＳ０の入力を受け付ける入力手段１、評
価を行うための参照画像を表す参照画像データＳ１を記
憶した記憶手段２、評価対象画像と参照画像との位置合
わせを行う位置合わせ手段３、位置合わせ後の評価対象
画像および参照画像について、ＣＤＲＡＤファントムの
格子線により区画される複数の格子領域を検出する格子
領域検出手段４、各格子領域について信号領域の位置を
検出して検出結果を取得する検出手段５、検出手段５に
おける信号領域位置の検出結果から正誤結果を求めて出
力する正誤結果出力手段６、および正誤結果に基づいて
評価結果を生成してこれを出力する評価結果出力手段７
を備える。
【００５３】入力手段１は、評価対象画像データＳ０の
入力を受け付けるものである。具体的には、評価対象画
像データＳ０が記憶されたメディアから評価対象画像デ
ータＳ０を読み出すメディアドライブ、ネットワーク経
由で転送された評価対象画像データＳ０を受け付けるイ
ンターフェース等が用いられる。
入力を受け付けるものである。具体的には、評価対象画
像データＳ０が記憶されたメディアから評価対象画像デ
ータＳ０を読み出すメディアドライブ、ネットワーク経
由で転送された評価対象画像データＳ０を受け付けるイ
ンターフェース等が用いられる。
【００５４】評価対象画像データＳ０は、ＣＤＲＡＤフ
ァントムの放射線画像を撮影することにより得られるも
のであり、ＣＤＲＡＤファントム画像を表すものであ
る。ＣＤＲＡＤファントム画像を図２に示す。ＣＤＲＡ
Ｄファントムに形成されている穴の直径は下側ほど小さ
く、穴の深さは左側ほど小さいため、図２に示すように
左下の格子領域ほど穴の像が見えにくくなっている。な
お、ＣＤＲＡＤファントム画像において格子領域の左側
の数字は穴の直径を、下側の数字は穴の深さを示す。ま
た、各格子領域においてＣＤＲＡＤファントムの穴の位
置に対応する位置が信号領域となっている。
ァントムの放射線画像を撮影することにより得られるも
のであり、ＣＤＲＡＤファントム画像を表すものであ
る。ＣＤＲＡＤファントム画像を図２に示す。ＣＤＲＡ
Ｄファントムに形成されている穴の直径は下側ほど小さ
く、穴の深さは左側ほど小さいため、図２に示すように
左下の格子領域ほど穴の像が見えにくくなっている。な
お、ＣＤＲＡＤファントム画像において格子領域の左側
の数字は穴の直径を、下側の数字は穴の深さを示す。ま
た、各格子領域においてＣＤＲＡＤファントムの穴の位
置に対応する位置が信号領域となっている。
【００５５】記憶手段２は参照画像データＳ１を記憶す
るものであり、具体的にはコンピュータに装備されたハ
ードディスクが用いられる。参照画像データＳ１は評価
対象画像データＳ０と同様に、ＣＤＲＡＤファントムの
放射線画像を撮影することにより取得されるが、評価対
象画像データＳ０を取得する場合よりも１０倍以上の高
い線量により撮影を行う等して、評価対象画像よりもＳ
／Ｎが高いものとなっている。したがって、参照画像は
評価対象画像よりもより小さい穴の像およびより深さが
小さい穴の像を視認することができる。
るものであり、具体的にはコンピュータに装備されたハ
ードディスクが用いられる。参照画像データＳ１は評価
対象画像データＳ０と同様に、ＣＤＲＡＤファントムの
放射線画像を撮影することにより取得されるが、評価対
象画像データＳ０を取得する場合よりも１０倍以上の高
い線量により撮影を行う等して、評価対象画像よりもＳ
／Ｎが高いものとなっている。したがって、参照画像は
評価対象画像よりもより小さい穴の像およびより深さが
小さい穴の像を視認することができる。
【００５６】なお、参照画像データＳ１は、全ての穴の
像を視認可能なように、ＣＤＲＡＤファントムを参照し
て人為的に作成してもよい。
像を視認可能なように、ＣＤＲＡＤファントムを参照し
て人為的に作成してもよい。
【００５７】位置合わせ手段３は、評価対象画像を基準
として参照画像と評価対象画像との位置合わせを行う。
具体的には、線が交差する位置を検出するクロス型フィ
ルタにより評価対象画像をフィルタリング処理して評価
対象画像上の格子点を検出する。そして格子点を基準と
したテンプレート領域を設定し、このテンプレート領域
を用いて参照画像上においてテンプレートマッチングを
行い、参照画像上に評価対象画像の格子点に対応する対
応格子点を設定する。そして、アフィン変換により参照
画像上における対応格子点の位置を評価対象画像上にお
ける格子点の位置に合わせることにより、参照画像およ
び評価対象画像の位置合わせを行う。なお、クロス型フ
ィルタを用いた位置合わせ方法は、特開平８−８３３３
６号公報に記載されている。
として参照画像と評価対象画像との位置合わせを行う。
具体的には、線が交差する位置を検出するクロス型フィ
ルタにより評価対象画像をフィルタリング処理して評価
対象画像上の格子点を検出する。そして格子点を基準と
したテンプレート領域を設定し、このテンプレート領域
を用いて参照画像上においてテンプレートマッチングを
行い、参照画像上に評価対象画像の格子点に対応する対
応格子点を設定する。そして、アフィン変換により参照
画像上における対応格子点の位置を評価対象画像上にお
ける格子点の位置に合わせることにより、参照画像およ
び評価対象画像の位置合わせを行う。なお、クロス型フ
ィルタを用いた位置合わせ方法は、特開平８−８３３３
６号公報に記載されている。
【００５８】なお、位置合わせは参照画像を基準とし、
評価対象画像を参照画像に合わせるように行ってもよ
い。
評価対象画像を参照画像に合わせるように行ってもよ
い。
【００５９】格子領域検出手段４は、以下のようにして
格子領域を検出する。まず、図３に示すように、位置合
わせされた参照画像を対象として、参照画像の縦方向お
よび横方向における中心付近３行分程度についての、参
照画像の平均プロファイルを求める。ＣＤＲＡＤファン
トムの外側の領域は直接放射線が照射される領域である
ためプロファイルの信号値は大きくなり、格子線におい
てはプロファイルの信号値は小さくなる。また、平均プ
ロファイルを求めることにより、ＣＤＲＡＤファントム
の穴の像および画像上のノイズのような局所的な信号は
平滑化されるため、ＣＤＲＡＤファントムの表面につい
ての平均プロファイルは、格子線の信号値よりも大きく
略一定の信号値を有するものとなる。なお、平均プロフ
ァイルは参照画像の全領域について求めてもよいが、演
算時間短縮のために中心付近の３行程度の領域における
平均プロファイルを求めているものである。
格子領域を検出する。まず、図３に示すように、位置合
わせされた参照画像を対象として、参照画像の縦方向お
よび横方向における中心付近３行分程度についての、参
照画像の平均プロファイルを求める。ＣＤＲＡＤファン
トムの外側の領域は直接放射線が照射される領域である
ためプロファイルの信号値は大きくなり、格子線におい
てはプロファイルの信号値は小さくなる。また、平均プ
ロファイルを求めることにより、ＣＤＲＡＤファントム
の穴の像および画像上のノイズのような局所的な信号は
平滑化されるため、ＣＤＲＡＤファントムの表面につい
ての平均プロファイルは、格子線の信号値よりも大きく
略一定の信号値を有するものとなる。なお、平均プロフ
ァイルは参照画像の全領域について求めてもよいが、演
算時間短縮のために中心付近の３行程度の領域における
平均プロファイルを求めているものである。
【００６０】次いで、しきい値Ｔｈ０を用いた２値化処
理により、平均プロファイルから格子線の位置を検出す
る。ＣＤＲＡＤファントムは格子領域が１５×１５存在
するため、１６本の格子線が参照画像の縦横においてそ
れぞれ検出される。しきい値Ｔｈ０は、平均プロファイ
ルの最小値（格子線の最小値）とＣＤＲＡＤファントム
表面の平均信号値との中間値を用いる。なお、ＣＤＲＡ
Ｄファントムに含まれる文字の部分を誤検出する場合が
あるため、しきい値Ｔｈ０を設定後、１６本の格子線が
検出されるまでしきい値Ｔｈ０の値を段階的に小さく変
化させるようにしてもよい。そして、縦横それぞれ１６
本の格子線を検出した後、格子線により区画される２２
５の領域を格子領域として検出する。なお、格子領域の
４隅は、参照画像上に設定される座標値により表すこと
ができる。
理により、平均プロファイルから格子線の位置を検出す
る。ＣＤＲＡＤファントムは格子領域が１５×１５存在
するため、１６本の格子線が参照画像の縦横においてそ
れぞれ検出される。しきい値Ｔｈ０は、平均プロファイ
ルの最小値（格子線の最小値）とＣＤＲＡＤファントム
表面の平均信号値との中間値を用いる。なお、ＣＤＲＡ
Ｄファントムに含まれる文字の部分を誤検出する場合が
あるため、しきい値Ｔｈ０を設定後、１６本の格子線が
検出されるまでしきい値Ｔｈ０の値を段階的に小さく変
化させるようにしてもよい。そして、縦横それぞれ１６
本の格子線を検出した後、格子線により区画される２２
５の領域を格子領域として検出する。なお、格子領域の
４隅は、参照画像上に設定される座標値により表すこと
ができる。
【００６１】なお、評価対象画像についても参照画像と
同様に平均プロファイルの算出および格子線の検出を行
うことにより格子領域を検出する。
同様に平均プロファイルの算出および格子線の検出を行
うことにより格子領域を検出する。
【００６２】検出手段５は下記のようにして信号領域の
位置を検出する。なお、信号領域の位置の検出は、ＣＤ
ＲＡＤファントム画像における下側の１２行、すなわち
穴の像を２つ含む格子領域について行われる。まず、評
価対象画像および参照画像の相対応する格子領域につい
て、穴の像を表す信号が存在する信号領域を検出する。
図４は格子領域を示す図である。評価対象画像および参
照画像の格子領域内において信号領域が存在する可能性
がある領域は、中心位置を除いては４隅のいずれか１箇
所であるため、図４に破線で示すように、評価対象画像
および参照画像の格子領域の４隅の近傍に、信号領域検
出のための所定範囲を有する検出領域１０Ａ〜１０Ｄ、
１１Ａ〜１１Ｄがそれぞれ設定される。なお、検出領域
１０Ａ〜１０Ｄ、１１Ａ〜１１Ｄは、格子領域の４隅の
座標値を基準として設定される。
位置を検出する。なお、信号領域の位置の検出は、ＣＤ
ＲＡＤファントム画像における下側の１２行、すなわち
穴の像を２つ含む格子領域について行われる。まず、評
価対象画像および参照画像の相対応する格子領域につい
て、穴の像を表す信号が存在する信号領域を検出する。
図４は格子領域を示す図である。評価対象画像および参
照画像の格子領域内において信号領域が存在する可能性
がある領域は、中心位置を除いては４隅のいずれか１箇
所であるため、図４に破線で示すように、評価対象画像
および参照画像の格子領域の４隅の近傍に、信号領域検
出のための所定範囲を有する検出領域１０Ａ〜１０Ｄ、
１１Ａ〜１１Ｄがそれぞれ設定される。なお、検出領域
１０Ａ〜１０Ｄ、１１Ａ〜１１Ｄは、格子領域の４隅の
座標値を基準として設定される。
【００６３】そして、図５に示すように、相対応する格
子領域における相対応する検出領域１０Ａ〜１０Ｄ、１
１Ａ〜１１Ｄ毎に、検出領域１０Ａ〜１０Ｄ、１１Ａ〜
１１Ｄ内の信号値の正規化相互相関値Ｃを下記の式
（１）により算出し、相関値Ｃが最大となる領域を、そ
の格子領域において穴の像が存在する信号領域として検
出する。
子領域における相対応する検出領域１０Ａ〜１０Ｄ、１
１Ａ〜１１Ｄ毎に、検出領域１０Ａ〜１０Ｄ、１１Ａ〜
１１Ｄ内の信号値の正規化相互相関値Ｃを下記の式
（１）により算出し、相関値Ｃが最大となる領域を、そ
の格子領域において穴の像が存在する信号領域として検
出する。
【００６４】
【数１】
但し、ｆ(x,y)：評価対象画像における検出領域内の信
号値 ｆave：評価対象画像における検出領域内の信号値の平
均値 ｇ(x,y)：参照画像における検出領域内の信号値 ｇave：参照画像における検出領域内の信号値の平均値
号値 ｆave：評価対象画像における検出領域内の信号値の平
均値 ｇ(x,y)：参照画像における検出領域内の信号値 ｇave：参照画像における検出領域内の信号値の平均値
【００６５】なお、図５に示す格子領域の場合、検出領
域１０Ｂについては穴の像が存在するため、検出領域１
１Ｂとの相関値Ｃは１に近い値となる。一方、他の検出
領域１０Ａ，１０Ｃ，１０Ｄについては、画像中に存在
するノイズの影響により、検出領域１１Ａ，１１Ｃ，１
１Ｄとの相関値Ｃは０に近い値となる。したがって、図
５に示す格子領域の場合、評価対象画像の格子領域にお
ける検出領域１０Ｂが信号領域として検出される。
域１０Ｂについては穴の像が存在するため、検出領域１
１Ｂとの相関値Ｃは１に近い値となる。一方、他の検出
領域１０Ａ，１０Ｃ，１０Ｄについては、画像中に存在
するノイズの影響により、検出領域１１Ａ，１１Ｃ，１
１Ｄとの相関値Ｃは０に近い値となる。したがって、図
５に示す格子領域の場合、評価対象画像の格子領域にお
ける検出領域１０Ｂが信号領域として検出される。
【００６６】ここで、図６に示すように、評価対象画像
の格子領域にゴミが存在する場合、目視評価では検出領
域１０Ｂ，１０Ｃのどちらが信号領域であるかが分から
ない。これに対して、参照画像の対応する格子領域内の
検出領域との正規化相互相関値Ｃを求めた場合、検出領
域１０Ｂと検出領域１１Ｂとの相関値Ｃは１に近い値と
なるが、検出領域１０Ｃと検出領域１１Ｃとの相関値は
０に近い値となる。したがって、目視により画質を評価
する場合よりも精度よく信号領域の位置を検出すること
ができる。
の格子領域にゴミが存在する場合、目視評価では検出領
域１０Ｂ，１０Ｃのどちらが信号領域であるかが分から
ない。これに対して、参照画像の対応する格子領域内の
検出領域との正規化相互相関値Ｃを求めた場合、検出領
域１０Ｂと検出領域１１Ｂとの相関値Ｃは１に近い値と
なるが、検出領域１０Ｃと検出領域１１Ｃとの相関値は
０に近い値となる。したがって、目視により画質を評価
する場合よりも精度よく信号領域の位置を検出すること
ができる。
【００６７】なお、上から３段目までの格子領域につい
ては、評価対象画像および参照画像の格子領域の中心位
置に検出領域を設定し、相対応する格子領域の検出領域
の相関値Ｃを算出する。
ては、評価対象画像および参照画像の格子領域の中心位
置に検出領域を設定し、相対応する格子領域の検出領域
の相関値Ｃを算出する。
【００６８】一方、信号領域の検出を下記のようにして
行ってもよい。以下説明する方法を方法２とし、上述し
た方法を方法１とする。まず、図７（ａ）に示すよう
に、評価対象画像の格子領域に検出領域１０Ａ〜１０Ｄ
を設定する。一方、参照画像の格子領域については、図
７（ｂ）に示すように、格子領域の中心位置にある穴の
像を囲む検出領域１２を設定する。なお、検出領域１０
Ａ〜１０Ｄと検出領域１２とは同一の大きさを有する。
そして、相対応する格子領域毎に検出領域１０Ａ〜１０
Ｄ内の信号値と検出領域１２内の信号値との正規化相互
相関値Ｃを上記の式（１）により算出し、相関値Ｃが最
大となる検出領域を、その格子領域において穴の像が存
在する信号領域として検出する。図７に示す格子領域の
場合、評価対象画像の格子領域における検出領域１０Ｂ
が信号領域として検出される。
行ってもよい。以下説明する方法を方法２とし、上述し
た方法を方法１とする。まず、図７（ａ）に示すよう
に、評価対象画像の格子領域に検出領域１０Ａ〜１０Ｄ
を設定する。一方、参照画像の格子領域については、図
７（ｂ）に示すように、格子領域の中心位置にある穴の
像を囲む検出領域１２を設定する。なお、検出領域１０
Ａ〜１０Ｄと検出領域１２とは同一の大きさを有する。
そして、相対応する格子領域毎に検出領域１０Ａ〜１０
Ｄ内の信号値と検出領域１２内の信号値との正規化相互
相関値Ｃを上記の式（１）により算出し、相関値Ｃが最
大となる検出領域を、その格子領域において穴の像が存
在する信号領域として検出する。図７に示す格子領域の
場合、評価対象画像の格子領域における検出領域１０Ｂ
が信号領域として検出される。
【００６９】ここで、図６に示すように、評価対象画像
の格子領域にゴミが存在する場合、検出領域１０Ｂのみ
ならず検出領域１０Ｃと検出領域１２との相関値Ｃが大
きくなるため、検出領域１０Ｃを誤って信号領域として
検出するおそれがある。したがって、方法２は方法１と
比較すると検出精度は落ちる。
の格子領域にゴミが存在する場合、検出領域１０Ｂのみ
ならず検出領域１０Ｃと検出領域１２との相関値Ｃが大
きくなるため、検出領域１０Ｃを誤って信号領域として
検出するおそれがある。したがって、方法２は方法１と
比較すると検出精度は落ちる。
【００７０】また、図８（ａ）に示すように、穴の径が
小さくなると検出領域１０Ａ〜１０Ｄに対する信号領域
の割合が小さくなるため、算出される相関値Ｃが小さく
なり、信号領域を精度よく検出することができなくな
る。このため、図８（ｂ）に示すように、検出領域１０
Ａ〜１０Ｄのサイズを穴の像の直径すなわち信号領域の
サイズに応じて変化させることにより、相関値Ｃを精度
よく算出でき、その結果信号領域をより精度よく検出す
ることができる。
小さくなると検出領域１０Ａ〜１０Ｄに対する信号領域
の割合が小さくなるため、算出される相関値Ｃが小さく
なり、信号領域を精度よく検出することができなくな
る。このため、図８（ｂ）に示すように、検出領域１０
Ａ〜１０Ｄのサイズを穴の像の直径すなわち信号領域の
サイズに応じて変化させることにより、相関値Ｃを精度
よく算出でき、その結果信号領域をより精度よく検出す
ることができる。
【００７１】そして、全ての格子領域について信号領域
を検出し、検出した信号領域の位置を検出結果テーブル
として出力する。図９は検出結果テーブルを示す図であ
る。検出結果テーブルにおいて、各格子領域に記載され
ている「０」，「１」，「２」，「３」は、検出した信
号領域の位置が格子領域の４隅の左上、右上、右下、左
下である場合をそれぞれ示す。
を検出し、検出した信号領域の位置を検出結果テーブル
として出力する。図９は検出結果テーブルを示す図であ
る。検出結果テーブルにおいて、各格子領域に記載され
ている「０」，「１」，「２」，「３」は、検出した信
号領域の位置が格子領域の４隅の左上、右上、右下、左
下である場合をそれぞれ示す。
【００７２】ここで、格子領域によっては、複数の検出
領域１０Ａ〜１０Ｄにおいて相関値Ｃが同一の値となっ
たり、ノイズ等の影響により相関値Ｃが１よりも大きく
なってしまう場合がある。このような場合は、その格子
領域については信号領域の位置が検出できなかったもの
として、「９」を出力する。なお、上から３段目までの
格子領域については、ノイズ等の影響により相関値Ｃが
１よりも大きくなってしまった場合に「９」を出力し、
それ以外の場合には「０」を出力する。
領域１０Ａ〜１０Ｄにおいて相関値Ｃが同一の値となっ
たり、ノイズ等の影響により相関値Ｃが１よりも大きく
なってしまう場合がある。このような場合は、その格子
領域については信号領域の位置が検出できなかったもの
として、「９」を出力する。なお、上から３段目までの
格子領域については、ノイズ等の影響により相関値Ｃが
１よりも大きくなってしまった場合に「９」を出力し、
それ以外の場合には「０」を出力する。
【００７３】正誤結果出力手段６は、検出手段５から出
力された検出結果テーブルと正解テーブルＴtとを用い
て信号検出位置の正誤結果を求める。
力された検出結果テーブルと正解テーブルＴtとを用い
て信号検出位置の正誤結果を求める。
【００７４】図１０は、信号領域の位置の正解位置を表
す正解テーブルＴtを示す図である。なお、正解テーブ
ルＴtは、ＣＤＲＡＤファントムを観察することによ
り、予め人為的に作成されてなり、記憶手段２に記憶さ
れているものである。そして、検出結果テーブルと正解
テーブルＴtとを比較し、各格子領域についての正誤結
果を図１１に示すような正誤テーブルＴtfとして生成す
る。図１１においてＴが正解を、Ｆが不正解を示す。
す正解テーブルＴtを示す図である。なお、正解テーブ
ルＴtは、ＣＤＲＡＤファントムを観察することによ
り、予め人為的に作成されてなり、記憶手段２に記憶さ
れているものである。そして、検出結果テーブルと正解
テーブルＴtとを比較し、各格子領域についての正誤結
果を図１１に示すような正誤テーブルＴtfとして生成す
る。図１１においてＴが正解を、Ｆが不正解を示す。
【００７５】なお、正誤結果においては、１／４の確率
により偶然に正解となってしまう場合がある。これを排
除するために、以下の（１）から（４）のルールにより
正誤結果を補正する。
により偶然に正解となってしまう場合がある。これを排
除するために、以下の（１）から（４）のルールにより
正誤結果を補正する。
【００７６】（１）正誤を判定する格子領域（以下注目
格子領域とする）が正解の場合、上下左右の最近傍の４
つの格子領域のうち２以上の格子領域が正解であれば、
その注目領域も正解とする。なお、注目格子領域がＣＤ
ＲＡＤファントム画像の辺に接している場合には、最近
傍の３つの格子領域のうち２以上の格子領域が正解であ
れば、その注目領域も正解とする。
格子領域とする）が正解の場合、上下左右の最近傍の４
つの格子領域のうち２以上の格子領域が正解であれば、
その注目領域も正解とする。なお、注目格子領域がＣＤ
ＲＡＤファントム画像の辺に接している場合には、最近
傍の３つの格子領域のうち２以上の格子領域が正解であ
れば、その注目領域も正解とする。
【００７７】（２）最近傍の格子領域が２つしかない、
すなわち注目格子領域がＣＤＲＡＤファントム画像の４
隅にある場合については、注目格子領域が正解の場合、
最近傍の格子領域の１つが正解であれば、注目格子領域
も正解とする。
すなわち注目格子領域がＣＤＲＡＤファントム画像の４
隅にある場合については、注目格子領域が正解の場合、
最近傍の格子領域の１つが正解であれば、注目格子領域
も正解とする。
【００７８】（３）注目格子領域が不正解の場合、最近
傍の４つまたは３つの格子領域のうち３以上の格子領域
が正解であれば、注目格子領域を正解とする。
傍の４つまたは３つの格子領域のうち３以上の格子領域
が正解であれば、注目格子領域を正解とする。
【００７９】（４）注目格子領域がＣＤＲＡＤファント
ム画像の４隅にある場合において、注目格子領域が不正
解の場合、最近傍の格子領域の全てが正解であれば注目
格子領域も正解とする。
ム画像の４隅にある場合において、注目格子領域が不正
解の場合、最近傍の格子領域の全てが正解であれば注目
格子領域も正解とする。
【００８０】図１２は補正後の正誤テーブルＴtfを示す
図である。図１２において太線で囲んだ格子領域が補正
された格子領域である。補正された正誤テーブルＴtfは
評価結果出力手段７に出力される。
図である。図１２において太線で囲んだ格子領域が補正
された格子領域である。補正された正誤テーブルＴtfは
評価結果出力手段７に出力される。
【００８１】評価結果出力手段７においては、入力され
た正誤テーブルＴtfに基づいて、正解であった穴の直径
および穴の深さの限界値を両対数グラフにプロットした
ＣＤ曲線を評価結果として出力する。具体的には、全て
の穴の深さについて、正解となった格子領域のうち、最
も穴の径が小さい格子領域についてその穴の直径の値
（以下限界値とする）を求め、穴の深さの対数値を横軸
に、穴の直径の対数値を縦軸にとった両対数グラフ上
に、各穴の直径の限界値をプロットすることにより、Ｃ
Ｄ曲線を作成することができる。
た正誤テーブルＴtfに基づいて、正解であった穴の直径
および穴の深さの限界値を両対数グラフにプロットした
ＣＤ曲線を評価結果として出力する。具体的には、全て
の穴の深さについて、正解となった格子領域のうち、最
も穴の径が小さい格子領域についてその穴の直径の値
（以下限界値とする）を求め、穴の深さの対数値を横軸
に、穴の直径の対数値を縦軸にとった両対数グラフ上
に、各穴の直径の限界値をプロットすることにより、Ｃ
Ｄ曲線を作成することができる。
【００８２】図１３は、図１２に示す正誤テーブルＴtf
から求めたＣＤ曲線を示す図である。ＣＤ曲線は、より
小さい直径かつより深さが小さい穴が見えるほど、グラ
フの左下に接近する。したがって、ＣＤ曲線がグラフの
左下に接近するほど、評価対象画像データＳ０を得たシ
ステムが高画質の画像を得ることができるものであると
いうことが分かる。
から求めたＣＤ曲線を示す図である。ＣＤ曲線は、より
小さい直径かつより深さが小さい穴が見えるほど、グラ
フの左下に接近する。したがって、ＣＤ曲線がグラフの
左下に接近するほど、評価対象画像データＳ０を得たシ
ステムが高画質の画像を得ることができるものであると
いうことが分かる。
【００８３】なお、ＣＤ曲線とともに、信号領域の位置
の正解率および確信度を評価結果に含めてもよい。ここ
で、確信度とは、２２５の全ての格子領域において求め
た相関値Ｃのうち、正解であった格子領域における相関
値の総和、または｛（正解であった格子領域における相
関値の総和）／（正解であった格子領域の数）｝×１０
０として算出される。
の正解率および確信度を評価結果に含めてもよい。ここ
で、確信度とは、２２５の全ての格子領域において求め
た相関値Ｃのうち、正解であった格子領域における相関
値の総和、または｛（正解であった格子領域における相
関値の総和）／（正解であった格子領域の数）｝×１０
０として算出される。
【００８４】次いで、本実施形態の動作について説明す
る。図１４は本実施形態の動作を示すフローチャートで
ある。まず、評価対象画像データＳ０の入力を受け付け
（ステップＳ１）、記憶手段２に記憶された参照画像デ
ータＳ１を用いて、位置合わせ手段３において、評価対
象画像および参照画像の位置合わせを行う（ステップＳ
２）。位置合わせ後、格子領域検出手段４において、評
価対象画像および参照画像から格子領域を検出する（ス
テップＳ３）。次いで、検出手段５において、評価対象
画像の各格子領域から信号領域の位置を検出する（ステ
ップＳ４）。そして、正誤結果出力手段６において、検
出結果に基づいて信号領域の位置の正誤結果を出力し
（ステップＳ５）、評価結果出力手段７において正誤結
果に基づいて評価結果を出力し（ステップＳ６）、処理
を終了する。
る。図１４は本実施形態の動作を示すフローチャートで
ある。まず、評価対象画像データＳ０の入力を受け付け
（ステップＳ１）、記憶手段２に記憶された参照画像デ
ータＳ１を用いて、位置合わせ手段３において、評価対
象画像および参照画像の位置合わせを行う（ステップＳ
２）。位置合わせ後、格子領域検出手段４において、評
価対象画像および参照画像から格子領域を検出する（ス
テップＳ３）。次いで、検出手段５において、評価対象
画像の各格子領域から信号領域の位置を検出する（ステ
ップＳ４）。そして、正誤結果出力手段６において、検
出結果に基づいて信号領域の位置の正誤結果を出力し
（ステップＳ５）、評価結果出力手段７において正誤結
果に基づいて評価結果を出力し（ステップＳ６）、処理
を終了する。
【００８５】このように、本実施形態によれば、評価対
象画像データＳ０および参照画像データＳ１に基づいて
放射線画像撮影装置や画像処理装置において得られる画
像の評価を行うようにしたため、評価対象画像の評価を
行う観察者の主観的な要素や評価環境が評価結果に反映
されることがなくなり、その結果、より客観的かつ安定
した評価結果を得ることができる。また、評価対象画像
データＳ０の入力を行うのみで客観性の高い評価結果を
得ることができるため、多数の評価対象画像を多数の評
価者によって評価する必要がなくなり、その結果、画像
の評価を効率よく行うことができる。
象画像データＳ０および参照画像データＳ１に基づいて
放射線画像撮影装置や画像処理装置において得られる画
像の評価を行うようにしたため、評価対象画像の評価を
行う観察者の主観的な要素や評価環境が評価結果に反映
されることがなくなり、その結果、より客観的かつ安定
した評価結果を得ることができる。また、評価対象画像
データＳ０の入力を行うのみで客観性の高い評価結果を
得ることができるため、多数の評価対象画像を多数の評
価者によって評価する必要がなくなり、その結果、画像
の評価を効率よく行うことができる。
【００８６】なお、上記実施形態においてはＣＤＲＡＤ
ファントムを撮影することにより、評価対象画像データ
Ｓ０を得ているが、画像処理装置において得られる画像
の評価を行う場合においては、ＣＤＲＡＤファントム画
像を表すＣＤＲＡＤファントム画像データを人為的に生
成し、このＣＤＲＡＤファントム画像データに対して、
評価を行う画像処理装置において画像処理を施すことに
より得られた処理済み画像データを評価対象画像データ
Ｓ０として用いることもできる。この場合、参照画像デ
ータＳ１としては、画像処理前のＣＤＲＡＤファントム
画像データをそのまま用いることができる。
ファントムを撮影することにより、評価対象画像データ
Ｓ０を得ているが、画像処理装置において得られる画像
の評価を行う場合においては、ＣＤＲＡＤファントム画
像を表すＣＤＲＡＤファントム画像データを人為的に生
成し、このＣＤＲＡＤファントム画像データに対して、
評価を行う画像処理装置において画像処理を施すことに
より得られた処理済み画像データを評価対象画像データ
Ｓ０として用いることもできる。この場合、参照画像デ
ータＳ１としては、画像処理前のＣＤＲＡＤファントム
画像データをそのまま用いることができる。
【００８７】また、上記実施形態においては、ＣＤＲＡ
Ｄファントムを撮影することにより評価対象画像データ
Ｓ０を得ているが、ＣＤＲＡＤファントムのみならず、
例えば、格子領域の４隅のいずれかの位置にのみ穴が形
成されたファントムを撮影することにより評価対象画像
データＳ０を得てもよい。このように、本実施形態にお
いては、評価対象画像および参照画像としては、少なく
とも１つの信号領域が所定位置に配設された複数の評価
領域からなるパターンであって、信号領域のサイズおよ
びコントラストが各評価領域毎に段階的に異なるパター
ンを有する画像であれば、いかなる画像をも用いること
ができる。
Ｄファントムを撮影することにより評価対象画像データ
Ｓ０を得ているが、ＣＤＲＡＤファントムのみならず、
例えば、格子領域の４隅のいずれかの位置にのみ穴が形
成されたファントムを撮影することにより評価対象画像
データＳ０を得てもよい。このように、本実施形態にお
いては、評価対象画像および参照画像としては、少なく
とも１つの信号領域が所定位置に配設された複数の評価
領域からなるパターンであって、信号領域のサイズおよ
びコントラストが各評価領域毎に段階的に異なるパター
ンを有する画像であれば、いかなる画像をも用いること
ができる。
【図１】本発明の実施形態による画像評価装置の構成を
示す概略ブロック図
示す概略ブロック図
【図２】ＣＤＲＡＤファントム画像を示す図
【図３】平均プロファイルの算出を説明するための図
【図４】格子領域を示す図
【図５】相関値の算出を説明するための図
【図６】ゴミが存在する格子領域の例を示す図
【図７】（ａ）は評価対象画像の格子領域に設定された
検出領域を示す図、（ｂ）は参照画像の格子領域に設定
された検出領域を示す図
検出領域を示す図、（ｂ）は参照画像の格子領域に設定
された検出領域を示す図
【図８】信号領域のサイズに応じた検出領域のサイズの
変更を説明するための図
変更を説明するための図
【図９】検出結果テーブルを示す図
【図１０】正解テーブルを示す図
【図１１】正誤テーブルを示す図
【図１２】補正された正誤テーブルを示す図
【図１３】ＣＤ曲線を示す図
【図１４】本実施形態の動作を示すフローチャート
１ 入力手段
２ 記憶手段
３ 位置合わせ手段
４ 格子領域検出手段
５ 検出手段
６ 正誤結果出力手段
７ 評価結果出力手段
１０Ａ〜１０Ｄ，１１Ａ〜１１Ｄ，１２ 検出領域
Claims (24)
- 【請求項１】 少なくとも１つの信号領域が所定位置
に配設された複数の評価領域からなるパターンであっ
て、前記信号領域のサイズおよびコントラストが該各評
価領域毎に段階的に異なるパターンについての評価対象
画像を表す評価対象画像データの入力を受け付け、 前記パターンについての参照画像を表す参照画像データ
および前記評価対象画像データに基づいて、前記評価対
象画像における全ての前記評価領域に存在する前記信号
領域の位置を検出し、 該検出結果に基づいて、前記評価対象画像の評価を行う
ことを特徴とする画像評価方法。 - 【請求項２】 前記参照画像および前記評価対象画像
の位置合わせを行った後に、前記信号領域の位置を検出
することを特徴とする請求項１記載の画像評価方法。 - 【請求項３】 前記評価対象画像および前記参照画像
のそれぞれの前記評価領域を分割して、前記評価対象画
像および前記参照画像のそれぞれの前記評価領域におい
て相対応する複数の分割領域を得、 前記評価対象画像および前記参照画像の相対応する前記
分割領域内の信号値に基づいて、前記信号領域の位置を
検出することを特徴とする請求項１または２記載の画像
評価方法。 - 【請求項４】 前記信号領域のサイズに応じて，前記
分割領域のサイズを変更することを特徴とする請求項３
記載の画像評価方法。 - 【請求項５】 前記評価対象画像および前記参照画像
の相対応する前記分割領域毎の複数の相関値を算出し、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出することを特徴と
する請求項３または４記載の画像評価方法。 - 【請求項６】 前記評価対象領域画像の前記評価領域
を１つの領域に１つの前記信号領域が含まれるように分
割して略同一サイズの複数の分割領域を得、前記参照画
像の前記評価領域から、該評価領域に対応する前記評価
対象画像の前記評価領域における前記分割領域に対応す
るサイズを有し、１つの前記信号領域を含む検定領域を
得、 前記評価対象画像および前記参照画像の相対応する前記
評価領域における前記複数の分割領域および前記検定領
域の信号値に基づいて、前記信号領域の位置を検出する
ことを特徴とする請求項１または２記載の画像評価方
法。 - 【請求項７】 前記信号領域のサイズに応じて，前記
分割領域のサイズおよび前記検定領域のサイズを変更す
ることを特徴とする請求項６記載の画像評価方法。 - 【請求項８】 前記評価対象画像および前記参照画像
の相対応する前記評価領域における前記複数の分割領域
および前記検定領域との複数の相関値を算出し、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出することを特徴と
する請求項６または７記載の画像評価方法。 - 【請求項９】 少なくとも１つの信号領域が所定位置
に配設された複数の評価領域からなるパターンであっ
て、前記信号領域のサイズおよびコントラストが該各評
価領域毎に段階的に異なるパターンについての評価対象
画像を表す評価対象画像データの入力を受け付ける入力
手段と、 前記パターンについての参照画像を表す参照画像データ
および前記評価対象画像データに基づいて、前記評価対
象画像における全ての前記評価領域に存在する前記信号
領域の位置を検出する検出手段と、 該検出結果に基づいて、前記評価対象画像の評価を行う
評価手段とを備えたことを特徴とする画像評価装置。 - 【請求項１０】 前記参照画像および前記評価対象画
像の位置合わせを行う位置合わせ手段をさらに備え、 前記検出手段は、該位置合わせ手段による位置合わせを
行った後に、前記信号領域の位置を検出する手段である
ことを特徴とする請求項９記載の画像評価装置。 - 【請求項１１】 前記検出手段は、前記評価対象画像
および前記参照画像のそれぞれの前記評価領域を分割し
て、前記評価対象画像および前記参照画像のそれぞれの
前記評価領域において相対応する複数の分割領域を得、 前記評価対象画像および前記参照画像の相対応する前記
分割領域内の信号値に基づいて、前記信号領域の位置を
検出する手段であることを特徴とする請求項９または１
０記載の画像評価装置。 - 【請求項１２】 前記検出手段は、前記信号領域のサ
イズに応じて，前記分割領域のサイズを変更する手段で
あることを特徴とする請求項１１記載の画像評価装置。 - 【請求項１３】 前記検出手段は、前記評価対象画像
および前記参照画像の相対応する前記分割領域毎の複数
の相関値を算出し、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出する手段であるこ
とを特徴とする請求項１１または１２記載の画像評価装
置。 - 【請求項１４】 前記検出手段は、前記評価対象領域
画像の前記評価領域を１つの領域に１つの前記信号領域
が含まれるように分割して略同一サイズの複数の分割領
域を得、 前記参照画像の前記評価領域から、該評価領域に対応す
る前記評価対象画像の前記評価領域における前記分割領
域に対応するサイズを有し、１つの前記信号領域を含む
検定領域を得、 前記評価対象画像および前記参照画像の相対応する前記
評価領域における前記複数の分割領域および前記検定領
域の信号値に基づいて、前記信号領域の位置を検出する
手段であることを特徴とする請求項９または１０記載の
画像評価装置。 - 【請求項１５】 前記検出手段は、前記信号領域のサ
イズに応じて，前記分割領域のサイズおよび前記検定領
域のサイズを変更する手段であることを特徴とする請求
項１４記載の画像評価装置。 - 【請求項１６】 前記検出手段は、前記評価対象画像
および前記参照画像の相対応する前記評価領域における
前記複数の分割領域および前記検定領域との複数の相関
値を算出し、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出する手段であるこ
とを特徴とする請求項１４または１５記載の画像評価装
置。 - 【請求項１７】 少なくとも１つの信号領域が所定位
置に配設された複数の評価領域からなるパターンであっ
て、前記信号領域のサイズおよびコントラストが該各評
価領域毎に段階的に異なるパターンについての評価対象
画像を表す評価対象画像データの入力を受け付ける手順
と、 前記パターンについての参照画像を表す参照画像データ
および前記評価対象画像データに基づいて、前記評価対
象画像における全ての前記評価領域に存在する前記信号
領域の位置を検出する手順と、 該検出結果に基づいて、前記評価対象画像の評価を行う
手順とを有する画像評価方法をコンピュータに実行させ
るためのプログラム。 - 【請求項１８】 前記参照画像および前記評価対象画
像の位置合わせを行う手順をさらに有し、 前記検出する手順は、該位置合わせを行った後に、前記
信号領域の位置を検出する手順である請求項１７記載の
プログラム。 - 【請求項１９】 前記検出する手順は、前記評価対象
画像および前記参照画像のそれぞれの前記評価領域を分
割して、前記評価対象画像および前記参照画像のそれぞ
れの前記評価領域において相対応する複数の分割領域を
得る手順と、 前記評価対象画像および前記参照画像の相対応する前記
分割領域内の信号値に基づいて、前記信号領域の位置を
検出する手順とを有する請求項１７または１８記載のプ
ログラム。 - 【請求項２０】 前記分割領域を得る手順は、前記信
号領域のサイズに応じて，前記分割領域のサイズを変更
する手順である請求項１９記載のプログラム。 - 【請求項２１】 前記検出する手順は、前記評価対象
画像および前記参照画像の相対応する前記分割領域毎の
複数の相関値を算出する手順と、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出する手順とを有す
る請求項１９または２０記載のプログラム。 - 【請求項２２】 前記検出する手順は、前記評価対象
領域画像の前記評価領域を１つの領域に１つの前記信号
領域が含まれるように分割して略同一サイズの複数の分
割領域を得る手順と、 前記参照画像の前記評価領域から、該評価領域に対応す
る前記評価対象画像の前記評価領域における前記分割領
域に対応するサイズを有し、１つの前記信号領域を含む
検定領域を得る手順と、 前記評価対象画像および前記参照画像の相対応する前記
評価領域における前記複数の分割領域および前記検定領
域の信号値に基づいて、前記信号領域の位置を検出する
手順とを有する請求項１７または１８記載のプログラ
ム。 - 【請求項２３】 前記分割領域を得る手順は、前記信
号領域のサイズに応じて，前記分割領域のサイズおよび
前記検定領域のサイズを変更する手順である請求項２２
記載のプログラム。 - 【請求項２４】 前記検出する手順は、前記評価対象
画像および前記参照画像の相対応する前記評価領域にお
ける前記複数の分割領域および前記検定領域との複数の
相関値を算出する手順と、 該複数の相関値のうち最大相関値を得た前記分割領域の
位置を前記信号領域の位置として検出する手順とを有す
る請求項２２または２３記載のプログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001337418A JP3945622B2 (ja) | 2001-11-02 | 2001-11-02 | 画像評価方法および装置並びにプログラム |
| US10/285,492 US7158691B2 (en) | 2001-11-02 | 2002-11-01 | Image evaluating method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001337418A JP3945622B2 (ja) | 2001-11-02 | 2001-11-02 | 画像評価方法および装置並びにプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003141535A true JP2003141535A (ja) | 2003-05-16 |
| JP3945622B2 JP3945622B2 (ja) | 2007-07-18 |
Family
ID=19152061
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001337418A Expired - Fee Related JP3945622B2 (ja) | 2001-11-02 | 2001-11-02 | 画像評価方法および装置並びにプログラム |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US7158691B2 (ja) |
| JP (1) | JP3945622B2 (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006527947A (ja) * | 2003-06-18 | 2006-12-07 | ブリティッシュ・テレコミュニケーションズ・パブリック・リミテッド・カンパニー | 映像品質評価方法及びシステム |
| JP2008043430A (ja) * | 2006-08-11 | 2008-02-28 | Shimadzu Corp | Ｘ線画像評価方法およびその装置 |
| JP2011524267A (ja) * | 2008-05-28 | 2011-09-01 | ペッパール ウント フュフス ゲゼルシャフト ミット ベシュレンクテル ハフツング | 印刷物を点検するための方法および機構とコンピュータプログラムとコンピュータプログラム製品 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004357868A (ja) * | 2003-06-03 | 2004-12-24 | Fuji Photo Film Co Ltd | Ｑｃファントム |
| US20090268953A1 (en) * | 2008-04-24 | 2009-10-29 | Apteryx, Inc. | Method for the automatic adjustment of image parameter settings in an imaging system |
| AT506865B1 (de) * | 2008-05-20 | 2010-02-15 | Siemens Vai Metals Tech Gmbh | Vorrichtung zum verbessern von genauigkeitseigenschaften von handhabungsgeräten |
| US9076195B2 (en) * | 2013-08-29 | 2015-07-07 | The Boeing Company | Methods and apparatus to identify components from images of the components |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4331021A (en) * | 1980-09-11 | 1982-05-25 | The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services | Contrast resolution tissue equivalent ultrasound test object |
| US5526442A (en) * | 1993-10-04 | 1996-06-11 | Hitachi Medical Corporation | X-ray radiography method and system |
| US5651046A (en) * | 1995-06-22 | 1997-07-22 | Duke University | Anatomic phantom for evaluation of projection radiographic imaging systems |
| US6694047B1 (en) * | 1999-07-15 | 2004-02-17 | General Electric Company | Method and apparatus for automated image quality evaluation of X-ray systems using any of multiple phantoms |
| US6674834B1 (en) * | 2000-03-31 | 2004-01-06 | Ge Medical Systems Global Technology Company, Llc | Phantom and method for evaluating calcium scoring |
-
2001
- 2001-11-02 JP JP2001337418A patent/JP3945622B2/ja not_active Expired - Fee Related
-
2002
- 2002-11-01 US US10/285,492 patent/US7158691B2/en not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006527947A (ja) * | 2003-06-18 | 2006-12-07 | ブリティッシュ・テレコミュニケーションズ・パブリック・リミテッド・カンパニー | 映像品質評価方法及びシステム |
| JP4741483B2 (ja) * | 2003-06-18 | 2011-08-03 | ブリティッシュ・テレコミュニケーションズ・パブリック・リミテッド・カンパニー | 映像品質評価方法及びシステム |
| US8508597B2 (en) | 2003-06-18 | 2013-08-13 | British Telecommunications Public Limited Company | Method and system for video quality assessment |
| JP2008043430A (ja) * | 2006-08-11 | 2008-02-28 | Shimadzu Corp | Ｘ線画像評価方法およびその装置 |
| JP2011524267A (ja) * | 2008-05-28 | 2011-09-01 | ペッパール ウント フュフス ゲゼルシャフト ミット ベシュレンクテル ハフツング | 印刷物を点検するための方法および機構とコンピュータプログラムとコンピュータプログラム製品 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3945622B2 (ja) | 2007-07-18 |
| US20030086626A1 (en) | 2003-05-08 |
| US7158691B2 (en) | 2007-01-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3326070B2 (ja) | 画像診断支援装置の画像処理方法 | |
| EP0671707B1 (en) | Method for adjusting positions of radiation images | |
| JP4601134B2 (ja) | 形状特徴に基づく欠陥検出の方法及び装置 | |
| US20020059520A1 (en) | Informarion processing apparatus, method for controlling the same, and storage medium | |
| EP1557792A2 (en) | Medical image processing apparatus | |
| EP1958158B1 (en) | Method for detecting streaks in digital images | |
| US20080012967A1 (en) | Defective-area correction apparatus, method and program and radiation detection apparatus | |
| JP2005296605A (ja) | 放射線写真画像を診断関連領域と診断非関連領域とにセグメント化する方法 | |
| KR20060125525A (ko) | 조사 영역 추출 방법 및 화상 처리 장치 | |
| KR102307995B1 (ko) | 딥러닝을 이용한 갑상선 암의 림프절 전이 진단 시스템 및 이의 동작 방법 | |
| JP4767721B2 (ja) | ディジタル医用検査画像の画像処理方法 | |
| US20200305706A1 (en) | Image processing method for glaucoma detection and computer program products thereof | |
| JP3945622B2 (ja) | 画像評価方法および装置並びにプログラム | |
| JP2005095501A (ja) | 医用画像処理装置 | |
| US7248729B2 (en) | Image processing apparatus and image processing method | |
| JP2009039411A (ja) | 放射線画像の品質評価方法、品質評価装置及びプログラム | |
| JP2001202507A (ja) | 画像の位置合わせ処理方法および位置合わせ処理装置 | |
| JP3239186B2 (ja) | 放射線画像の照射野領域抽出装置 | |
| EP0702321B1 (en) | Method and apparatus for classifying picture elements in radiation images | |
| JPH03102477A (ja) | 放射線画像処理装置 | |
| JP2007202811A (ja) | 照射野認識装置、照射野認識方法およびそのプログラム | |
| EP2856429B1 (en) | Methods and apparatus for image processing, and laser scanning ophthalmoscope having an image processing apparatus | |
| JP2004302581A (ja) | 画像処理方法並びに画像処理装置 | |
| JP2006293522A (ja) | 直線検出装置、直線検出方法およびそのプログラム | |
| JP2006288467A (ja) | 照射野判定装置、照射野判定方法およびそのプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040305 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061019 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061031 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061204 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061226 |
|
| 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: 20070403 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070404 |
|
| 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: 20110420 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120420 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130420 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |