(57)【特許請求の範囲】
【請求項１】
コンピュータの情報処理を用いてデジタル画像の品質を評価する処理を行う検査画像品質評価システムであって、
検査対象物から得られた検査の対象となるデジタル画像のデータを入力する入力部と、
上記デジタル画像のデータから、デジタル化される前のアナログ像を推定して当該推定アナログ画像のデータを得る処理を行うアナログ像推定部と、
上記デジタル画像と上記推定アナログ画像とを周波数空間上で比較することにより、当該推定アナログ画像からデジタル画像への劣化の度合いを定量化した相対的な品質の評価値を計算する評価値計算部と、
上記評価値を含むデータ情報を人に対し出力する出力部と、を有すること、を特徴とする検査画像品質評価システム。
【請求項２】
請求項１記載の検査画像品質評価システムにおいて、
前記品質の評価値についての品質のレベルを計算するレベル計算部と、
前記レベルを含むデータ情報を人に対し出力する前記出力部と、を有し、
前記レベル計算部は、上記評価値を、判定基準情報と比較して上記レベルを判定する処理を行うこと、を特徴とする検査画像品質評価システム。
【請求項３】
請求項１記載の検査画像品質評価システムにおいて、
前記評価値計算部は、
前記デジタル画像と推定アナログ画像とを周波数空間で比較して比較値を得る比較部と、
前記デジタル画像から品質を表す１つ以上の特徴量を抽出する特徴抽出部と、
前記特徴量に重みをつけて加算することにより統合指標値を算出し、当該統合指標値と上記比較値とを統合することで前記品質の評価値として定量化する処理を行う評価値定量化部と、を有すること、を特徴とする検査画像品質評価システム。
【請求項４】
請求項２記載の検査画像品質評価システムにおいて、
前記判定基準情報を事前に設定する前処理部を有し、
前記前処理部は、前記対象の画像に関する人による目視の評価値を、前記判定基準情報として設定すること、を特徴とする検査画像品質評価システム。
【請求項５】
請求項２記載の検査画像品質評価システムにおいて、
前記判定基準情報を事前に設定する前処理部を有し、
前記前処理部は、前記対象の画像に関する人による目視の評価値を含む教師情報を設定し、
前記教師情報を入力として学習処理により前記判定基準情報を生成する学習部を有すること、を特徴とする検査画像品質評価システム。
【請求項６】
請求項５記載の検査画像品質評価システムにおいて、
前記教師情報として、対象のデジタル画像とアナログ画像の組と、それらに対する人の目視による品質のレベルとのセットを入力すること、を特徴とする検査画像品質評価システム。
【請求項７】
請求項６記載の検査画像品質評価システムにおいて、
前記学習部は、前記セットのデータ情報を入力して学習処理することにより、前記レベルの判定のための判定基準情報を生成し、
上記学習処理により、前記品質の評価値を統合指標値で算出する場合の特徴量の重みを決定すること、を特徴とする検査画像品質評価システム。
【請求項８】
請求項１記載の検査画像品質評価システムにおいて、
前記検査の対象となる画像は、放射線検出装置からの検出結果をデジタイザによりデジタル化したデジタル画像であり、
前記推定アナログ画像のデータを対象として欠陥部を検出する処理を行う検査部を有すること、を特徴とする検査画像品質評価システム。
【請求項９】
コンピュータの情報処理を用いて、デジタイザにより得られたデジタル画像の品質を評価することによりデジタイザの性能を保証する処理を行うデジタイザ保証システムであって、
前記デジタイザにより得られたデジタル画像のデータと、前記デジタイザの性能を評価するための評価パターン画像データと、を入力する入力部と、
上記デジタル画像のデータから、前記デジタイザによりデジタル化される前のアナログ像を推定して当該推定アナログ画像のデータを得る処理を行うアナログ像推定部と、
前記評価パターン画像データを用いた人の目視の評価値を入力して判定基準情報として設定する処理を行う前処理部と、
上記デジタル画像と上記推定アナログ画像とを周波数空間上で比較することにより、当該推定アナログ画像からデジタル画像への劣化の度合いを定量化した相対的な品質の評価値を計算する評価値計算部と、
前記品質の評価値を前記判定基準情報と比較することにより前記品質の評価値についてのレベルを判定するレベル計算部と、
上記レベルに基づき、前記デジタイザによる画質の劣化度合いを推定し、当該デジタイザの性能の良否を表す指標値を決定し、当該指標値を含む情報を人に対し出力する処理を行うデジタイザ保証部と、を有すること、を特徴とするデジタイザ保証システム。
【請求項１０】
コンピュータの情報処理を用いてデジタル画像の品質を評価する処理を行う検査画像品質評価方法であって、
検査対象物から得られた検査の対象となるデジタル画像のデータを入力する入力ステップと、
上記デジタル画像のデータから、デジタル化される前のアナログ像を推定して当該推定アナログ画像のデータを得る処理を行うアナログ像推定ステップと、
上記デジタル画像と上記推定アナログ画像とを周波数空間上で比較することにより、当該推定アナログ画像からデジタル画像への劣化の度合いを定量化した相対的な品質の評価値を計算する評価値計算ステップと、
上記品質の評価値を判定基準情報と比較することにより上記品質の評価値についての品質のレベルを判定するレベル計算ステップと、
上記評価値、レベルを含むデータ情報を人に対し出力する出力ステップと、を有すること、を特徴とする検査画像品質評価方法。
【請求項１１】
コンピュータの情報処理を用いてデジタル画像の品質を評価する処理を実行させる検査画像品質評価プログラムであって、
検査対象物から得られた検査の対象となるデジタル画像のデータを入力する処理を行うプログラムと、
上記デジタル画像のデータから、デジタル化される前のアナログ像を推定して当該推定アナログ画像のデータを得る処理を行うプログラムと、
上記デジタル画像と上記推定アナログ画像とを周波数空間上で比較することにより、当該推定アナログ画像からデジタル画像への劣化の度合いを定量化した相対的な品質の評価値を計算するプログラムと、
上記品質の評価値を判定基準情報と比較することにより上記品質の評価値についての品質のレベルを判定するプログラムと、
上記評価値、レベルを含むデータ情報を人に対し出力するプログラムと、を有すること、を特徴とする検査画像品質評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コンピュータ情報処理（デジタル画像データ処理）を用いて画像の品質を評価する技術に関する。特に、アナログ・デジタル変換手段であるデジタイザを用いて得られた対象物（被検査物）の画像の検査・診断等のために当該画像の品質（画質）を定量化・保証等する技術に関する。また特に、放射線検出装置を用いた放射線透過試験（非破壊検査）に使用する画像の品質を評価する技術に関する。
【背景技術】
【０００２】
例えば原子力プラントの製造現場や医療の現場では、放射線透過試験（非破壊検査）・診断等が行われている。その際、従来は、放射線検出装置を用いて対象物に対して放射線を透過して得られたフィルム像（アナログ像）を対象（検査画像）として検査員により目視確認することにより不良部（欠陥・疾患等）などを判定・検出していた。
【０００３】
また近年では、フィルムスキャナやイメージングプレート等に代表される種々のデジタイザが開発されており、フィルム像（アナログ像）を元にデジタイザにより得られたデジタル像（画像データ）を対象とした人による目視検査・診断等へと移行してきている。デジタル像の目視検査は、操作性や保管・検索のしやすさ等の観点で優位である。
【０００４】
更に今後は、上記デジタル像を対象としたデータ処理（画像解析処理等）により、人の目視に頼らない自動的な検査・診断等が加速するものと考えられている。
【０００５】
ここで、上記のようなフィルム像からデジタル像による検査への移行で問題となるのが、アナログからデジタルへの変換に伴う画質（品質）の劣化である。デジタル像を用いることにより、従来のフィルム像の目視判定では検出できていた微小な欠陥などが検出できなくなる可能性がある。
【０００６】
このため、フィルム像に対するデジタル像の相対的な画質（品質）を定量化することにより、当該デジタル像が検査・診断等に値すべき画質（品質）であるかどうか、判定や保証や作業支援等が実現できる方式などの必要性が生じている。
【０００７】
例えば原子力プラントの配管の放射線透過試験（非破壊検査）に使用する対象の画像（検査画像）についての品質を評価・保証するために好適な方式（検査画像品質評価方法など）が求められている。
【０００８】
上記デジタル像の画質の定量化の手法としては、欧州規格（EN Standards）  EN14096-1（非特許文献１）において記載されているように、フィルムスキャナのようなデジタイザの性能をＭＴＦ（Modulation Transfer Function）に基づいて評価する手法がある。また、EN14784-1（非特許文献２）において記載されているように、ＭＴＦとDuplex wire methodに基づいて、イメージングプレート（デジタイザ）とそれに付随する周辺機器・ソフトウェア等を包括したシステムを評価する手法がある。
【０００９】
しかしこれらの手法は、デジタル像自体の画質を絶対値で評価する手法であり、フィルム像に対するデジタル像の画質の相対値（相対的な画質）を定量化して出力することはできない。
【００１０】
また従来の手法では、デジタル像の画質（相対的な画質）の評価値について、実際に目視で検査を行う検査員の評価値（官能評価値）と整合する値（概略一致する値）を計算・出力することは困難である。
【００１１】
上記デジタル像の画質の検査・評価等に関する先行技術例としては、特開２０００−４６６５０号公報（特許文献１）、特開２００４−２９８６１７公報（特許文献２）などがある。
【００１２】
特許文献１（「画質検査装置」）では、人間の視覚と近い画質判定ができる画質検査を自動的に行うこと、人間の視覚による官能検査と相関の取れる画質検査を自動的に行うこと、等が記載されている（要約参照）。
【００１３】
特許文献２（「放射線撮影システム」等）では、放射線撮影システムの検査において、画像の定量評価及び目視評価の両方を簡便かつ低コストで行うと共に、放射線撮影システムの不変性評価の精度を上げることができるファントム等を提供する、等と記載されている（要約参照）。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開２０００−４６６５０号公報
【特許文献２】特開２００４−２９８６１７号公報
【非特許文献】
【００１５】
【非特許文献１】欧州規格  EN14096-1
【非特許文献２】欧州規格  EN14784-1
【非特許文献３】M.Mori,T.Kashiwa and Y.Aoki,2007,”Digital Image Evaluation method for Digital Radiography”ASME,ICONE18 Xi’an,China,Mat2010
【発明の概要】
【発明が解決しようとする課題】
【００１６】
上述の通り、アナログ像またはデジタル像の品質の評価の例である放射線検査・診断等は、人（目視検査員や医師等）による官能評価であることが多い。特に原子力プラントの配管の検査などにおいては、資格を持った検査員の判定（目視判定）が絶対である。
【００１７】
このため、前述のデジタイザによるデジタル像を用いた検査・診断等を行う場合、従来の対象であったフィルム像（アナログ像）に対してそのデジタル像がどの程度劣化しているかといったことを表す画質の相対値（相対的な画質の評価値）を自動的に計算（定量化）・出力することで、当該デジタル像が検査・診断等に値すべき品質であるかどうか、判定や保証や作業支援等が実現できることが求められている。
【００１８】
また特に、上記デジタル像の相対的な画質の評価値の自動的な計算・出力を行う場合、及び上記デジタル像を用いた自動的な検査・診断等を行う場合には、従来のフィルム像を対象とした人（検査員等）の目視判定の結果（官能評価値、目視感）と整合するような結果（評価値等）が得られることが好ましい。
【００１９】
以上を鑑み、本発明の主な目的は、上記デジタル画像を用いた品質評価・検査などの技術において、（１）アナログ像（フィルム像）に対するデジタル像の品質（画質）の相対的な評価値などを自動的に計算・出力することができる技術、（２）好ましくは、人の目視の評価値（目視感）と整合する評価値などを得ることができる技術、を提供することである。
【００２０】
これにより、当該デジタル像が検査・診断等に値すべき品質であるかどうか、判定や保証や作業支援等が実現できることである。関連して、上記デジタル画像（検査画像）を用いた自動的な検査・診断等を実現できること、特に目視検査の場合と整合するように不良部などを検出できること、デジタイザに関する性能の保証ができること、等である。
【課題を解決するための手段】
【００２１】
上記目的を達成するため、本発明のうち代表的な形態は、コンピュータの情報処理を用いてデジタル画像の品質を評価する処理を行う検査画像品質評価システム（対応する方法、プログラム等）であって、以下に示す構成を有することを特徴とする。
【００２２】
（０）本システムは、デジタイザなどで取得したデジタル画像（検査画像）の品質を評価する処理を行うシステムであって、アナログ像（フィルム像）に対するデジタル像の相対的な品質（画質）の評価値を、人（検査員）の目視の評価値と整合するように計算・出力する処理機能を有する。本システムは、対象となるデジタル画像（デジタル像）を入力する入力部と、前記デジタル画像から、前記デジタイザなどでデジタル化する前のアナログ画像（アナログ像）を推定し推定アナログ像データを得る処理を行うアナログ像推定部と、前記デジタル画像と推定アナログ画像から上記相対的な品質の評価値（Ｈ）を計算する処理部と、前記評価値（Ｈ）に関するレベル（Ｌ）を決定するための判定基準情報を、人による当該デジタル画像及び対応アナログ画像の目視の評価結果（評価値）を既知の評価値として設定または学習処理することにより生成する前処理部と、前記評価値（Ｈ）を前記判定基準情報と比較して前記レベル（Ｌ）を判定するレベル判定部と、上記検査画像（デジタル画像または推定アナログ画像）を用いて自動的に欠陥の検出などの検査処理を行う検査部と、上記画像データ、評価値（Ｈ）、レベル（Ｌ）、検査結果情報などを人（検査員）に対し出力処理する出力部と、を有する。
【００２３】
前記評価値計算部は、詳しくは、前記デジタル画像と推定アナログ画像とを周波数空間で比較して比較値を得る比較部と、前記デジタル画像から品質や特徴を表す１つ以上の特徴量を抽出する特徴抽出部と、前記特徴量と前記周波数空間上の比較値とを統合して、前記デジタル画像の推定アナログ画像に対する相対的な品質の評価値（Ｈ）を算出する評価値定量化部とを有する。上記構成により、人の目視感と整合する評価値（Ｈ）、レベル（Ｌ）を得る。
【００２４】
（１）第１方式：  本システムは、検査等の対象となるデジタル画像のデータを入力する手段と、上記デジタル画像のデータから、デジタル化される前のフィルム像（アナログ像）を推定して当該推定アナログ画像のデータを得る処理を行う手段と、上記デジタル画像と上記推定アナログ画像とを周波数空間上で比較（例えば差分演算）することにより、アナログ像からデジタル像への劣化の度合いを定量化した相対的な品質の評価値（Ｈ）を計算する手段と、上記評価値（Ｈ）のデータ情報を人に対し出力する手段と、を有する。上記検査画像は、例えば欠陥などの不良部を検出するために放射線検出装置及びデジタイザにより取得されるデジタル画像である。
【００２５】
（２）また本システムは、上記品質の評価値（Ｈ）についてのレベル（Ｌ）、即ち検査等にあたり人の目視感に基づく判定基準に照らしてどの程度信頼できる値であるかを表すレベル（Ｌ）を計算する手段と、当該レベル（Ｌ）のデータ情報を人に対し出力する手段と、を有する。上記レベル（Ｌ）を計算する手段は、上記評価値（Ｈ）を判定基準情報（人の目視感に基づく設定情報）と比較することによりレベル（Ｌ）を判定する。
【００２６】
（３）また特に、上記評価値（Ｈ）の計算手段は、前記デジタル画像から、例えば分解能，コントラスト，ノイズ周波数などの複数の種類の画質の特徴量を抽出する処理を行う手段と、上記複数の各々の特徴量に重み（係数）をつけて加算することにより１つの統合指標値を算出し、当該統合指標値と、前記デジタル画像と推定アナログ画像との周波数空間上の比較結果とを統合して、前記品質の評価値（Ｈ）として算出する処理を行う手段とを有する。
【００２７】
（４）また、人の目視の評価値と整合するような評価値（Ｈ）、レベル（Ｌ）を得る手段として、人の目視の評価値を反映した判定基準情報を得る前処理部を有する。前処理部は、例えば、人による目視の評価値を事前に判定基準情報として設定する。あるいは、前処理部は、人による目視の評価値を教師情報として学習処理により判定基準情報を生成（更新）する。特に、評価値（Ｈ）が人の目視の評価値と整合するような回帰式を事前に決定する学習処理を行うことにより判定基準情報を算出する学習部を有する。前処理部は、例えば、複数のデジタル画像及び対応する複数のアナログ画像と、それらに対する人の目視による評価値とのセットを事前に教師情報として入力する。学習部は、上記セット（教師情報）の入力から、レベル（Ｌ）の判定のための判定基準情報を生成する。また前処理部は、例えば上記セットを用いて、前記統合指標値を算出する場合の各々の特徴量の重み（係数）を事前に算出・設定する。
【００２８】
（５）第２方式：  本システムは、第２方式の検査画像品質評価処理を行う第２方式処理部を有する。第２方式処理部は、デジタル画像（検査画像）を入力する処理と、デジタル画像から１つ以上の特徴量を抽出して品質評価値（Ｈ）（特に統合指標値）を計算する処理と、上記品質評価値（Ｈ）及び判定基準情報からレベル（Ｌ）を決定する処理と、レベル（Ｌ）情報を出力する処理と、を行う。
【００２９】
（６）デジタイザ保証方式：  本システムは、デジタイザの性能（アナログからデジタルへの変換に係わる画質の劣化の度合い）を評価し保証する処理を行うデジタイザ保証部を有する。デジタイザ保証部は、デジタイザによるデジタル画像を入力する処理と、デジタイザの性能を評価するための評価パターン画像データを入力する処理と、上記デジタル画像からデジタイザによってデジタル化する前のアナログ像を推定する処理と、前記評価パターンデータを用いた人の目視の評価値を入力する処理と、上記デジタル画像と推定アナログ画像とを周波数空間上で比較することによりアナログ像に対するデジタル画像の相対的な品質の評価値（Ｈ）を計算する処理と、当該評価値（Ｈ）を判定基準情報と比較することにより評価値（Ｈ）のレベル（Ｌ）を計算し、当該計算した値（Ｈ，Ｌ）に基づき、デジタイザによる画質の劣化度合いを推定し、性能の良否を表す指標値を決定する処理と、上記性能の良否を表す指標値を含む結果情報を人に対し出力する処理とを行う。
【発明の効果】
【００３０】
本発明のうち代表的な形態によれば、上記デジタル画像を用いた品質評価・検査などの技術において、（１）アナログ像に対するデジタル像の品質の相対的な評価値などを自動的に計算・出力することができる。（２）好ましくは、人の目視の評価値（目視感）と整合する評価値などを得ることができる。
【００３１】
これにより、当該デジタル像が検査・診断等に値すべき品質であるかどうか、判定や保証や作業支援等が実現できる。関連して、上記デジタル画像（検査画像）を用いた自動的な検査・診断等を実現でき、特に目視検査の場合と整合するように不良部などを検出でき、デジタイザに関する性能の保証ができる。
【図面の簡単な説明】
【００３２】
【図１】本発明の一実施の形態のシステム（検査画像品質評価システムを含んで成るシステム）の全体構成を示す図である。
【図２】アナログ像に対するデジタル像の画質の定量化の概念を示す図である。
【図３】一般的な画像劣化モデルの一例を示す図である。
【図４】本発明の実施の形態１（第１方式）における画像品質評価処理を示す図である。
【図５】実施の形態１におけるノイズ除去処理の一例を示す図である。
【図６】点像広がり関数（ＰＳＦ）と画像劣化の関係を示す図である。
【図７】実施の形態１における点像広がり関数（ＰＳＦ）逆変換によるアナログ像推定処理の一例を示す図である。
【図８】実施の形態１におけるパワースペクトルの比較処理の一例を示す図である。
【図９】実施の形態１における目視感と整合する値（Ｈ，Ｌ）の計算処理の一例を示す図である。
【図１０】実施の形態１におけるデジタル画像のレベル（Ｌ）を決定する処理の一例を示す図である。
【図１１】実施の形態１における前処理部による判定基準の設定処理の一例を示す図である。
【図１２】実施の形態１における検査部による検査処理の一例を示す図である。
【図１３】本発明の実施の形態２（第２方式）における、品質レベル決定処理の例を示す図である。
【図１４】本発明の実施の形態３のシステムにおける、デジタイザ保証処理の一例を示す図である。
【図１５】従来技術における放射線検査（アナログ像の目視検査）の基本的な流れを示す図である。
【発明を実施するための形態】
【００３３】
以下、本発明の実施の形態（検査画像品質評価システムなど）を図面に基づいて詳細に説明する。なお、実施の形態を説明するための全図において、同一部には原則として同一符号を付し、その繰り返しの説明は省略する。以下の実施の形態では、対象として、放射線透過による欠陥検査及びそのデジタル画像を例に説明する。
【００３４】
実施の形態１（図１〜図１２）では、第１方式の検査画像評価処理を行う。第１方式では、デジタル像からのアナログ像の推定処理を含む。実施の形態１の処理では、デジタル像（検査画像）の入力、デジタル像の品質の評価値（Ｈ）及びそのレベル値（Ｌ）の自動的な計算、その結果情報の出力処理、検査画像からの欠陥検出を含む検査処理などを行う。
【００３５】
実施の形態２（図１３）では、第２方式（特に欧州の規格に対応した方式）の検査画像評価処理を行う。第２方式では、第１方式と共通する考え方の処理を含みつつ、デジタル像からのアナログ像の推定処理を含まずに、デジタル像の品質の評価値（Ｈ）及びそのレベル値（Ｌ）の自動的な計算などを行う。
【００３６】
実施の形態３（図１４）のデジタイザ保証方式（システム）では、実施の形態１と同様の考え方での評価値（Ｈ）やレベル値（Ｌ）の計算処理などを用いつつ、デジタイザの性能を保証するための処理を行う。なお実施の形態１では、推定や評価の対象は、アナログ信号源とデジタイザとを含んだ広く不特定な範囲となる。それに対し、実施の形態３では、推定や評価の対象は、実施の形態１よりも限定された、デジタイザの部分となる。
【００３７】
＜従来方式（アナログ像の目視判定）＞
図１５は、前提技術として、従来方式であるアナログ像の目視検査（目視判定）を行う場合の放射線検査の基本的な流れ（概念）を示している。放射線検出装置９０１と、現像部９０２と、目視検査部９０３とを有する。
【００３８】
放射線検出装置９０１では、放射線源９２１から放射線を対象（検査対象物）９２２に照射し、透過された結果（強度分布）を、検出部９２３（放射線専用のフィルム）で検出する。検出部９２３で得られた出力（透過像）を、現像部９０２でアナログ像（フィルム像）Ａ１として現像する。そして、現像後のアナログ像（フィルム像）Ａ１を対象として、目視検査部９０３で検査員（Ｕ）により目視判定（目視検査）を行う。これによりアナログ像（フィルム像）Ａ１から欠陥部などが検出される。
【００３９】
また、デジタル像を対象として目視判定する形態の場合は、例えば現像部９０２からのアナログ像（フィルム像）Ａ１をフィルムスキャナ等（デジタイザ）でデジタル化したデジタル像を得る。そして目視検査部９０３で検査員（Ｕ）により当該デジタル像を対象とした目視判定を行う。
【００４０】
近年、工業製品の非破壊検査、医療現場でのレントゲン診断など、放射線による検査を行う現場では、上記放射線検出装置９０１の検出部９２３は、イメージングプレートなどのデジタイザにとって代わり、デジタイザにより検査画像として直接的にデジタル画像を得ることができるようになってきた。このデジタル画像を対象に目視検査を行う場合、現像部９０２が不要となる。同時に今後は、目視検査部９０３での人による目視判定が、コンピュータの画像データ処理による自動検査・診断等に置き換わることが想定される。
【００４１】
一方、デジタイザを用いたデジタル画像による検査では、従来のフィルム像（アナログ像）に比べてデータ画像の画質が劣化する可能性がある。このため、例えばフィルム像による目視検査からデジタル画像による自動検査へ移行するためには、フィルム像に対するデジタル画像の相対的な画質を定量化し、デジタイザより得られるデジタル画像が、検査等を行う上で必要な品質を保っているかどうかを定量的に評価・保証する必要がある。
【００４２】
＜実施の形態１＞
図１〜図１２を用いて、本発明の実施の形態１の検査画像品質評価システム、及び対応する検査画像品質評価方法などについて説明する。
【００４３】
［システム］
図１で、実施の形態１のシステムの全体構成例、及び検査の基本的な流れなどを示している。なお図１のシステムでは、従来のアナログ像やデジタル像の目視検査の構成や、後述の実施の形態２，３のシステムの機能をも統合したシステム形態（即ち各種の方式の検査が可能な形態）を示しているが、これらは適宜分離した形態が可能である。
【００４４】
本システム全体は、放射線検出装置（放射線検査装置）１０１と、現像装置１０２と、デジタイザ（イメージングプレート）１０３と、フィルムスキャナ１０４と、目視検査部１０５と、ＰＣ（検査装置）１とが接続される構成である。
【００４５】
ＰＣ（検査装置）１は、検査ＳＷ（ソフトウェア）１０を有し、そのソフトウェアプログラム処理により、図１や図４以降に示すような、検査画像品質評価処理、及び自動検査処理などの機能を提供する。ＰＣ１は、プロセッサ、メモリ、ディスク、入力装置、出力装置、通信インタフェース装置、バスなどの公知の要素で構成される。例えばプロセッサはディスクから本実施の形態のプログラムをメモリにロードして実行することにより検査ＳＷ１０を実現する。
【００４６】
図１のシステムでは、検査方式として、（１）アナログ像Ａ１の目視検査（従来技術）と、（２）デジタル像Ｄ１，Ｄ２の目視検査（従来技術）と、（３）本検査装置１によるデジタル像Ｄ１，Ｄ２の自動検査と、が可能である。なおデジタイザとして特にイメージングプレート（１０３）の場合で説明するが、フィルムスキャナ１０４等の他の手段の場合でも同様である。
【００４７】
放射線検出装置１０１は、図１５の放射線検出装置９０１を、デジタル対応に置き換えた構成であり、デジタイザ１０３に接続される。現像装置１０２は、図１５の現像部９０２と同様であり、放射線検出結果からフィルム像（アナログ像）Ａ１を得る。デジタイザ（イメージングプレート）１０３は、放射線検出装置１０１に応じたイメージングプレートであり、放射線検出結果からアナログ・デジタル変換によりデジタル像（データ）Ｄ１を得る。なお１０３は放射線検出装置１０１内の検出部として構成されてもよい。フィルムスキャナ１０４は、現像装置１０２に対応したデジタイザの例であり、フィルム像（アナログ像）Ａ１を入力しアナログ・デジタル変換によりデジタル像（データ）Ｄ２を得る。目視検査部１０５は、図１５の目視検査部９０３と同様であり、アナログ像Ａ１やデジタル像Ｄ１，Ｄ２を対象（検査画像）として検査員（Ｕ）の目視判定による検査を行う場合に使用する。
【００４８】
検査装置１の検査ＳＷ１０は、図示する各処理部（プログラムモジュール等）により各データ情報を処理する。主な処理部として、入力部１１、アナログ像推定部１２、評価値（Ｈ）計算部１３、レベル（Ｌ）計算部１４、設定部１５、出力部（表示部）１６、検査部（欠陥検出部）１７、第２方式処理部１８（実施の形態２の機能）、デジタイザ保証部１９（実施の形態３の機能）などを備える。また主な処理部を構成する副次的な処理部として、ノイズ除去部２１、逆変換部２２、ＦＦＴ部２３、比較部２４、特徴抽出部２５、評価値（Ｈ）定量化部２６、レベル（Ｌ）判定部２７、前処理部２８、学習部２９などを有する。記憶部に管理される各データ情報として、デジタル画像データ（ＧＤ）３１、推定アナログ画像データ（ＧＡ）３２、周波数スペクトルデータ３３、品質評価値（Ｈ）３４、品質レベル値（Ｌ）３５、判定基準情報３６、評価パターンデータ３７、検出欠陥情報３８、教師情報３９などを有する。
【００４９】
入力部１１は、検査装置１の外部から、デジタル像Ｄ１，Ｄ２等に対応する、検査対象のデジタル画像データ（ＧＡ）３１を入力・取得する処理を行う。
【００５０】
アナログ像推定部１２は、デジタル画像データ（ＧＡ）３１から、元のアナログ画像を推定する処理を行い、推定アナログ画像データ（ＧＡ）３２を得る。
【００５１】
評価値（Ｈ）計算部１３は、デジタル画像データ（ＧＡ）３１と推定アナログ画像データ（ＧＡ）３２とを用いて、検査対象のアナログ画像からデジタル画像への変換に係わる相対的な品質の評価値（Ｈ）を計算する処理を行う。言い換えると、Ｈ計算部１３は、アナログ像（ＧＡ）に対するデジタル像（ＧＤ）の画質の劣化度合いを評価し定量化した値（Ｈ）を求める。Ｈ計算部１３は、詳しくは、比較部２４、特徴抽出部２５、Ｈ定量化部２６などを用いて実現される。
【００５２】
レベル（Ｌ）計算部１４は、Ｈ計算部１３により得た対象画像の品質評価値（Ｈ）に関する品質レベル（Ｌ）を計算する処理を行う。Ｌ計算部１４は、詳しくは、レベル（Ｌ）判定部２７、前処理部２８などを用いて実現される。
【００５３】
品質レベル（Ｌ）は、本システム（Ｈ計算部１３）で自動的に求めた品質評価値（Ｈ）に関するレベル判定値であり、検査員（Ｕ）の目視評価値（目視感）を反映した判定基準（３６）に照らしたとき、どの程度信頼できるか（当該デジタル像が検査に値すべき品質かどうか）を表す値である。なお評価値（Ｈ）やレベル（Ｌ）は説明上の呼称であるため適宜言い換え可能である。
【００５４】
設定部１５は、本システムでの設定情報を検査装置１の表示画面で検査員（Ｕ）により入力・設定可能とする処理を行う。特に、設定部１５は、前処理部２８での判定基準情報３６またはその生成のための教師情報３９などを入力・設定可能である。
【００５５】
出力部（表示部）１６は、本検査ＳＷ１０での各種データ情報や処理結果情報などを検査装置１の画面に表示して検査員（Ｕ）により確認可能とするといった出力処理を行う。
【００５６】
検査部（欠陥検出部）１７は、検査対象の画像データ（ＧＤまたはＧＡ）から、欠陥などの不良部を自動的に検出してその検出欠陥情報３８を得る処理を行う。特に、検査部１７は、推定アナログ画像データ（ＧＡ）を対象として検査処理を行うことができ、これにより検査の高感度化を可能とする。検出欠陥情報３８は、例えば対象画像データから公知の画像解析処理等により検出された欠陥領域・座標などの情報である。
【００５７】
第２方式処理部１８は、実施の形態２の機能に対応し、デジタル画像から元のアナログ画像を推定せずに、評価値（Ｈ）及びレベル（Ｌ）などを計算する、第２方式での処理を行う。
【００５８】
実施の形態２では、特に、複数の特徴量（特徴抽出部２５で算出した値）から統合指標値として品質評価値（Ｈ０）を求め、当該統合指標値（Ｈ０）及び判定基準情報３６から品質レベル（Ｌ）を決定する。また、第２方式処理部１８（または設定部１５）は、統合指標値（Ｈ０）の計算で用いる重み値などを設定する機能も含む。
【００５９】
デジタイザ保証部１９は、実施の形態３の機能に対応し、上記各処理部、及び評価パターンデータ３７を用いて、デジタイザ１０３の性能（アナログ・デジタル変換での劣化度合いを示す）を、前述の評価値（Ｈ）及びレベル（Ｌ）の相当値を用いて計算（評価）し、当該性能を示す値（情報）を検査員（Ｕ）に対し表示する。これにより、検査員（Ｕ）によりデジタイザ１０３の性能の良否を容易に判定・確認可能としてデジタイザ１０３の性能を保証する。
【００６０】
ノイズ除去部２１は、デジタル画像データ（ＧＤ）から、ノイズを除去する処理を行う（後述、図５等）。
【００６１】
逆変換部２２は、ノイズ除去後のデジタル画像データ（ＧＤ’）から、劣化モデル（図３）に基づく逆演算（点像広がり関数（ＰＳＦ）逆変換処理、デコンボリューション処理）により、デジタル画像の元のアナログ画像を推定して、復元画像に相当する推定アナログ画像データ（ＧＡ）を得る処理を行う。
【００６２】
ＦＦＴ部２３は、ＦＦＴ（高速フーリエ変換）処理を行い、周波数スペクトル（パワースペクトル）データ３３を得る処理を行う。
【００６３】
比較部２４は、評価値（Ｈ）の算出のために、デジタル画像データ（ＧＤ）の周波数スペクトルデータ（３３Ｄ）と推定アナログ画像データ（ＧＡ）の周波数スペクトルデータ（３３Ａ）とを比較（例えば差分演算）する処理を行い（図４，図８等）、その比較値を得る。例えば、差分演算として、最大周波数の差、差の割合の平均値、などを算出する。
【００６４】
特徴抽出部２５は、評価値（Ｈ）の算出のために、デジタル画像データ（ＧＤ）及び推定アナログ画像データ（ＧＡ）から、１つ以上の種類の特徴量を抽出する処理を行う（図図１０等）。特徴量の種類としては、例えば、パワースペクトルの高周波成分、分散値などや、分解能、コントラスト、ノイズ周波数などがある。
【００６５】
評価値（Ｈ）定量化部２６は、比較部２４の結果（比較値）や、特徴抽出部２５の結果（特徴量）を用いて、最終的な品質評価値（Ｈ）を定量化する処理を行い、品質評価値（Ｈ）３４のデータ情報として得る（図８等）。例えば、評価値（Ｈ）定量化部２６は、特徴量と周波数スペクトル空間上の比較値を統合した品質評価値（Ｈ）を算出する。例えば、評価値（Ｈ）定量化部２６は、複数の種類の特徴量（ないし品質評価値）とそれぞれの重み値を用いて統合指標値を算出して品質評価値（Ｈ）とする（図１０等）。
【００６６】
レベル（Ｌ）判定部２７は、品質評価値（Ｈ）ないし特徴量と、判定基準情報３６とを比較することで、品質レベル（Ｌ）を判定（決定）する処理を行い、品質レベル値（Ｌ）３５のデータ情報として得る（図９，図１０等）。
【００６７】
前処理部２８は、前処理（検査や評価よりも事前の処理）により、判定基準情報３６（品質レベル（Ｌ）の判定のための情報）を得る処理を行う（図１１等）。特に、前処理部２８は、学習部２９を用いて、検査員（Ｕ）の目視による評価結果（官能評価値）を既知の評価値として学習させることにより、判定基準情報３６を生成する。例えば前処理として、検査員（Ｕ）により判定基準情報３６を事前に入力・設定することができる。また例えば前処理として、学習部２９で教師あり学習処理により判定基準情報３６を生成するための情報（学習用の教師情報３９）を設定することができる。学習用の教師情報３９として、例えば、対応するアナログ画像とデジタル画像の組と、それらを対象とした検査員の目視判定による品質評価値・レベル値と、のセットのデータ情報を設定できる。
【００６８】
学習部２９は、品質評価値（Ｈ）と、上記セットのデータ情報（学習用の教師情報３９）とを用いて、教師あり学習処理を行い、判定基準情報３６を生成（更新）する処理を行う。学習処理は、例えば、回帰式（回帰分析）や、ＳＶＭ（Support Vector Machine）などの手法を適用可能である。本システムでは、上記セット（学習用の教師情報３９）をもとに生成／設定される判定基準情報３６として、検査員（Ｕ）の目視判定による品質評価値（官能評価値）を反映させる形としている。したがって、本システムで自動的に計算される評価値（Ｈ）及びレベル（Ｌ）について、人の目視感と整合した値を得ることができる。
【００６９】
［アナログ像に対するデジタル像の画質の定量化］
図２は、本発明（実施の形態）で用いる、アナログ像（フィルム像）に対するデジタル像の画質の定量化の概念について示す。２０１はアナログ像（フィルム像）である（例えば図１のＡ１対応）。２０２は対応するデジタル画像である（例えばＤ１対応）。フィルム像２０１に対し、デジタル化の後であるデジタル画像２０２は、解像度が低下するなどの劣化が生じていると仮定する。２０３，２０４は、アナログ像２０１、デジタル画像２０２をそれぞれフーリエ変換（ＦＦＴ）して得られたパワースペクトル（周波数スペクトル）である。上記劣化などが生じている場合、デジタル画像は高周波成分が落ちるため、両者のパワースペクトル２０３，２０４を比較（２０５）することで、アナログ側の２０３に比べてデジタル側の２０４の高周波成分がどれだけ減少したかによって定量化が可能となる。
【００７０】
しかし、現実には、フィルム像２０１に対してデジタル処理をしてパワースペクトル２０３を得ることはできない。これに対し、本発明（実施の形態）では、デジタル画像２０２（ＧＤ）から劣化（デジタル化）前のアナログ像２０１（ＧＡ）を推定する処理を行う（アナログ像推定部１２）。当該推定したアナログ像（ＧＡ）からは、デジタル処理をしてパワースペクトル２０３（３３Ａ）を得ることができる。そして、当該推定アナログ像（ＧＡ）とデジタル画像（ＧＤ）とを周波数空間（パワースペクトル）で比較（２０５）することにより（比較部２４）、アナログ像（ＧＡ）に対するデジタル画像（ＧＤ）の相対的な品質（画質）を定量化し（Ｈ定量化部２６）、品質評価値（Ｈ）を算出する（Ｈ計算部１３）。品質評価値（Ｈ）の例としては、解像度○％、輝度分解能○％、最小欠陥検出サイズ○％等化といった、アナログ像を１００％としたときの画質を表す特徴の比率や、検出可能な最小欠陥サイズの比率などがある。
【００７１】
［一般的な劣化モデル］
図３（ａ）は、画像の一般的な劣化モデルを示す。３００は、例として２つの楕円形のパターンが撮像されたアナログ像である。３０１は、アナログ像３００のａ−ｂ上の断面輝度波形である。３０２は、アナログ像３００がデジタル化される際の劣化の例（点像広がり関数変換）を示す。３０２において、アナログ像３００に対し点像広がり関数（ＰＳＦ）３１３が畳み込まれること（３１１）により、点像のボケが生じ、さらにノイズ３１４が加算（３１２）される。３０５は、劣化（３０２）の後のデジタル像であり、３０６がその断面輝度波形である。デジタル像３０５はアナログ像３００に比べて解像度が下がっていることを示す。
【００７２】
そこで、本発明（実施の形態）では、図３（ｂ）に示すように、上記デジタル画像（ＧＤ）について、上記劣化モデル（図３（ａ））の逆演算（ＰＳＦ逆変換（Ｓ１４））により、アナログ像（ＧＡ）を推定する（アナログ像推定部１２）。
【００７３】
図３（ａ）の演算モデル（劣化モデル）をまとめると、［アナログ画像］＊［ＰＳＦ］＋［ノイズ］→［デジタル画像］である。＊（３１１）は畳み込み（コンボリューション）を示す。ＰＳＦは点像広がり関数の略である。図３（ｂ）の逆演算モデル（アナログ像推定モデル）をまとめると、（［デジタル画像］−［ノイズ］）＊［逆ＰＳＦ］→［推定アナログ画像（ＧＡ）］である。
【００７４】
［第１方式：品質評価値（Ｈ）計算］
図４は、実施の形態１（第１方式）における評価値（Ｈ）計算部１３の処理として、図２の画質定量化モデル及び図３（ｂ）のアナログ像推定モデルに基づく、デジタル画像（ＧＤ）の相対的な画質の評価処理の手順を示す（Ｓ１１等は手順などに対応する）。
【００７５】
Ｓ１１は、入力部１１により得られたデジタル画像データ（ＧＤ）である。Ｓ１１のデジタル画像（ＧＤ）を入力として、まずＳ１２でノイズ除去部２１によるノイズ除去処理を行う。これによりＳ１３でノイズが除去された画像データ（ＧＤ’）を得る。次にＳ１４で、ノイズ除去後の画像（ＧＤ’）を入力として、点像広がり関数（ＰＳＦ）の逆変換処理（デコンボリューション処理）を逆変換部２２により行う。これによりＳ１５で、推定アナログ画像データ（ＧＡ）（劣化の無い復元画像）を得る。Ｓ１３〜Ｓ１５の処理はアナログ像推定部１２の処理（逆変換部２２の処理を含む）である。
【００７６】
Ｓ１６，Ｓ１７では、Ｓ１１のデジタル画像（ＧＤ）と、Ｓ１５の推定アナログ画像（ＧＡ）との各々をフーリエ変換（ＦＦＴ）処理する（ＦＦＴ部２３の処理）。Ｓ１８の３３Ｄはデジタル画像（ＧＤ）より得られたパワースペクトル、Ｓ１９の３３Ａは、推定アナログ画像（ＧＡ）より得られたパワースペクトルである（周波数スペクトルデータ３３対応）。そしてＳ２０では、これらのパワースペクトル３３Ｄ，３３Ａを比較部２４で比較処理することにより、デジタル画像（ＧＤ）の相対的な画質を定量化し、品質評価値（Ｈ）を得る（Ｈ定量化部２６の処理）。
【００７７】
Ｓ２０の比較（定量化）の一例としては、アナログ側のパワースペクトル３３Ａの高周波成分に対してデジタル側のパワースペクトル３３Ｄの高周波成分がどこまでカットされたかで評価する。または、それぞれのパワースペクトルの分散値（各周波数に対するパワーの分散値）の比で評価する。更には、その両者の重み付き線形和で評価してもよい。いずれにせよ、パワースペクトル（３３）から１つ以上の特徴量を算出し（特徴抽出部２５の処理）、それらを統合した値で評価することで、最終的な品質評価値（Ｈ）を得る。算出された品質評価値（Ｈ）を含む情報は出力部１６により出力処理される。
【００７８】
［ノイズ除去］
図５は、Ｓ１２のノイズ除去部２１によるノイズ除去処理の例を示す。ここではノイズとしては白色ノイズを仮定している。図５（ａ）で、５００はノイズの無い画像のパワースペクトル、５０１は白色ノイズのパワースペクトルである。白色ノイズのパワースペクトルは全周波数領域で一定となる。５０２は白色ノイズ（５０１）が付加（加算）された画像のパワースペクトルである。５０１のノイズ成分の分だけ、全周波数領域でオフセットが加算される。従って、白色ノイズを除去するためには、そのオフセット量を減算すればよい。
【００７９】
図５（ｂ）で、ノイズを除去したい対象の画像５０３（ＧＤ）に対し、Ｓ１２のノイズ除去の手順では、まずＳ５０１でフーリエ変換（ＦＦＴ）処理し（ＦＦＴ部２３）、パワースペクトルを得る。そしてＳ５０２で、全周波数領域から一定のパワー（設定値５１１）を減算し、ノイズの除去された画像５０４（ＧＤ’）を得る。ここで、５１１として減算するオフセット量は、ある程度、ノイズがモデル化可能であれば事前にセットしておく。不明な場合は、得られたパワースペクトルから下位Ｎ％の値をオフセット値とする等、統計的に求めてもよい。また、ノイズとしては上記の白色ノイズに限定する必要は無く、画像取得時の特性によりそのモデルが既知であれば、その部分を除去すればよい。
【００８０】
［アナログ像推定］
以上（図５）の通りにノイズを除去した後、本発明（実施の形態）では、画像の劣化度を推定し、画像（アナログ像）を復元（推定）する（アナログ像推定部２１）。ここでは、画像劣化（ボケ）の要因となる点像広がり関数（ＰＳＦ）をガウス関数と仮定する。画像の劣化は、下記の式１の通り、平均０、分散σ２の２次元ガウス関数が畳み込まれることによって生じる。
【００８１】
【数１】
0005624975.tif 000002
【００８２】
ただし、式１におけるf(x,y)は劣化前の画像、G(x,y)は２次元ガウス関数（下記の式２）である。
【００８３】
【数２】
0005624975.tif 000003
【００８４】
ここで、上記のσが大きいほど、入力画像f(x,y)の解像度劣化は激しくなる。また、劣化が激しくなるほど、そのパワースペクトルは高周波成分が落ちることになる。
【００８５】
図６は、その作用（ＰＳＦと画像劣化の関係）を示す。第１ケースにおいて、６０１は、入力画像である劣化前の画像（f(x,y)）、６０２は、分散σ１の点像広がり関数（ＰＳＦ）（式２）、６０３は、劣化後の画像(F(x,y)）、６０４は、劣化後の画像（６０３）のフーリエ変換（ＦＦＴ）後のパワースペクトル（分散Ｖ１）である。同様に、第２ケースにおいて、６０５は、分散σ２の点像広がり関数（ＰＳＦ）、６０６は、劣化後の画像(F(x,y)）、６０７は、劣化後の画像（６０６）のフーリエ変換（ＦＦＴ）後のパワースペクトル（分散Ｖ２）である。
【００８６】
σ１≪σ２の場合（σ１よりもσ２の方がずっと大きい）、入力画像６０１に対し、劣化は６０３よりも６０６の方が大きく、パワースペクトルの分散は、反対にＶ１≫Ｖ２となることを示している。
【００８７】
これらのこと（第１、第２ケース）から、６０８に示すような関係式が成り立つ。すなわち、横軸をＰＳＦの分散σｉ、縦軸をＰＳＦにより劣化した画像のパワースペクトルの分散Ｖｉとすると、Ｖｉはσｉが大きくなるにつれて減少する。この関係式（６０８）を算出することで劣化後画像（デジタル画像（ＧＤ））から劣化前画像（アナログ画像（ＧＡ））を復元（推定）することが可能となる。例えば、Ｖｉとσｉが線形の関係にあれば、下記の（式３）の係数ａ，ｂを決定する。
【００８８】
【数３】
0005624975.tif 000004
【００８９】
上記の係数ａ，ｂは、シミュレーションをベースに求めてもよいし、事例（実際に何点かのＰＳＦとそれにより得られる劣化後画像のパワースペクトルの関係から算出する等）により求めてもよい。
【００９０】
［ＰＳＦ逆変換によるアナログ像推定］
図７は、アナログ像推定部１２（逆変換部２２）の処理として、図４のＳ１４の点像広がり関数（ＰＳＦ）逆変換によるアナログ像推定（画像復元）処理の例を示す。上述の通り、点像広がり関数（ＰＳＦ）はガウス関数を仮定する。７００は、復元対象となるデジタル画像（ＧＤ’）である。アナログ像推定部１２（逆変換部２２）では、７００（ＧＤ’）に対し、まずＳ７０１でフーリエ変換（ＦＦＴ）処理し（ＦＦＴ部２３）、パワースペクトルを得る。次にＳ７０２でパワースペクトルの分散を抽出する。分散、即ち周波数成分の広がりを抽出することで、画像７００（ＧＤ’）の劣化の度合いを得る。分散の演算方法の例として、下記の式４のＶを求める。
【００９１】
【数４】
0005624975.tif 000005
【００９２】
なお式４で、ω：各周波数、μ：平均周波数、Sx(ω)：各周波数ωのパワースペクトルである。
【００９３】
１２（２２）は、上記分散Ｖを、ｘ方向、ｙ方向の周波数それぞれより算出し、その二乗和をパワースペクトルの分散とする。本演算は、周波数成分の広がりを表す要素の一例であって、他の演算方法で定義してもよい。
【００９４】
次にＳ７０３で、１２（２２）は、上記抽出された分散値に応じて、点像広がり関数（ＰＳＦ）を算出する。実際にはＰＳＦをガウス関数と仮定しているので、そのσを算出すればよい。
【００９５】
このＰＳＦの算出方法の一例としては、Ｓ７０５に示す事前のキャリブレーションがある。例えば、基準のゲージ（評価パターンデータ３７）などを対象として撮像した画像に対し、σの大きさをいくつか変えたσ１，σ２，……のガウス関数を畳み込み、そのパワースペクトルＶ１，Ｖ２，……を得る。これらを図６の６０８のように関数にあてはめて関係式を推定しておく。１２（２２）は、このような事前のキャリブレーション（Ｓ７０５）により推定した関係式（６０８）を用いて、点像広がり関数（ＰＳＦ）のσを算出する。そして、１２（２２）は、そのσを用いて、Ｓ７０４で、畳み込み（ＰＳＦ）の逆変換を行うことにより、復元画像である推定アナログ画像データ７０１（ＧＡ）を得る。
【００９６】
以上のように、白色ノイズを仮定したノイズの除去（図５）と、ガウス関数を仮定した点像広がり関数（ＰＳＦ）の推定（図７）により、アナログ像（ＧＡ）を推定する。なお上記のアナログ像の推定の手法に限らず、他の手法の組み合わせであっても構わない。例えば、基本的なデコンボリューション法として、ヤコビ法、フーリエ変換法などを用いてもよい。評価したいデジタル画像からアナログ像を推定するアプローチをとり、それらを比較することでデジタル画像の品質を定量化することが本発明の範囲である。
【００９７】
［比較処理］
次に、図８を用いて、比較部２４を用いたデジタル画像（ＧＤ）と推定アナログ像（ＧＡ）との比較処理（図１のＳ２０）の一例を説明する。図１のＳ２０では、両画像（ＧＤ，ＧＡ）のパワースペクトル（３３Ｄ，３３Ａ）の違いによって、アナログ像に対するデジタル像の劣化度合いを評価している。
【００９８】
図８で、８００は上述の通りデジタル画像（ＧＤ）から推定したアナログ像（ＧＡ）を示す。８０１はデジタル画像（ＧＤ）を示す。８０２（３３Ａ），８０３（３３Ｄ）は、８００（ＧＡ），８０１（ＧＤ）の画像をそれぞれフーリエ変換（ＦＦＴ）処理し（Ｓ１７，Ｓ１６）、得られたパワースペクトルである。これら（８０２，８０３）に対し、Ｓ８０１の比較部２４の処理において、まずＳ８０１の差分演算で差分（比較値）８０４を求める。そしてＳ８０２で、Ｈ定量化部２６の処理により、差分（比較値）８０４から、アナログ像（ＧＡ）に対するデジタル画像（ＧＤ）の品質評価値（Ｈ）（最終的には１つ。例えば後述の統合指標値Ｈ０）を算出する。なおここで、手法（機能）に応じて、評価値（Ｈ）は、１つに限らず、２つ以上、算出することができる（２つ以上を算出する場合、区別のため、特徴量と称する）。
【００９９】
評価値（Ｈ）の算出方法の一例として、下記の式５の通り、最大周波数の差である。
【０１００】
【数５】
0005624975.tif 000006
【０１０１】
評価値（Ｈ）の別の算出方法の例として、下記の式６の通り、アナログ像に対する、差の割合の平均値である。
【０１０２】
【数６】
0005624975.tif 000007
【０１０３】
式５，式６の評価値（Ｈ）の算出例のいずれかでもよいし、両者を統合して算出してもよい。また別の品質評価値を算出してもよい。
【０１０４】
ここで、これらの評価値（Ｈ）は、従来は目視検査の対象となるフィルム像に対するデジタル画像の品質の劣化度である（本実施の形態ではＧＡに対するＧＤの劣化度）。このため、この劣化度（評価値Ｈ）であれば当該画像の検査を行う上で問題無いといったことを確認するためのレベル（Ｌ）（その判定）については、用途などに応じて検査員（Ｕ）の目視判定による評価値・レベルと整合させる必要がある場合もある。この整合が必要な場合（用途）に対応して、本発明（実施の形態）では、以下に詳述するように、検査員（Ｕ）の目視感と整合する評価値（Ｈ）・レベル（Ｌ）を計算・出力することが可能な構成である。
【０１０５】
［目視感と整合する値（Ｈ，Ｌ）の計算処理］
図９は、更に、Ｈ計算部１３、Ｌ計算部１４などによる処理として、前述のデジタル画像（ＧＤ）から計算される品質評価値（Ｈ）についての品質レベル（Ｌ）の判定について、検査員（Ｕ）の目視感と整合するように値（Ｌ）を決定する処理の流れを示す。
【０１０６】
Ｓ９００で、Ｈ計算部１３（Ｈ定量化部２６または特徴抽出部２５）は、図４等の例に従い、入力のデジタル画像（ＧＤ）をもとに、１つ以上の品質評価値（Ｈ）ないし特徴量を算出する。例えば、まず評価対象のデジタル画像（ＧＤ）が入力されると、Ｈ定量化部２６または特徴抽出部２５では、アナログ像を推定するために、ノイズを除去し（前記Ｓ１２）、ノイズ除去後のデジタル画像（ＧＤ’）について、点像広がり関数（ＰＳＦ）の推定と逆変換の処理によりアナログ像（ＧＡ）を推定する（前記Ｓ１４）。次に、推定アナログ像（ＧＡ）に対するデジタル像（ＧＤ）の品質評価値（Ｈ）を、両者のフーリエ変換（ＦＦＴ）後のパワースペクトルの比較により算出する（前記Ｓ２０）。この評価値Ｈはデジタル画像（ＧＤ）の特徴量を示す。ここでは、複数の特徴量を１つの値に統合したものを評価値Ｈと称す。
【０１０７】
Ｓ９０１では、Ｌ計算部１４（Ｌ判定部２７）は、上記の評価値Ｈ（特徴量）を、事前に設定してある判定基準（３６）と比較し、Ｓ９０２で、当該判定基準（３６）に基づいて、デジタル画像（ＧＤ）の評価値ＨについてのレベルＬを決定（判定）する。検査員（Ｕ）は、出力部１６により出力処理される上記結果情報（Ｈ，Ｌ）を参照して、品質判定（当該デジタル画像（ＧＤ）を検査に用いて問題無いかといった判断・確認）をすることができる。評価値Ｈは上述の通り、画像の品質を示す特徴量（計算値）そのものであり、品質（Ｌ）は評価値をより粗く量子化した値であり、検査員の目視感との整合が可能な値、すなわち、人間の評価分解能と同程度の値となる。一例としては、ＯＫ・ＮＧのように２値であってもいいし、悪・普・良のように３値でもいいし、レベル１・レベル２・レベル３……といった多値であってもいい。
【０１０８】
［レベル判定処理］
図１０は、上述のレベル（Ｌ）判定等に関する具体的な概念・処理例を示す。なおここおでは評価値（Ｈ）として複数の特徴量（例：２種類：Ａ，Ｂ）を抽出して統合評価値（０）を計算してレベルＬを判定する場合で説明するが、１つの評価値（特徴量）を対象としてレベルＬを判定する形も可能である。
【０１０９】
図１０で、１０００は、ある撮像条件で撮像して得られたデジタル画像（ＧＤ）である。このような画像（ＧＤ）（１０００）がＨ計算部１３（Ｈ定量化部２６または特徴抽出部２５）に入力されると、まず当該画像の各々について品質評価値Ｈを算出する。ここでは、特徴抽出部２５により、２つの異なる種類の評価値（特徴量）Ａ，Ｂ（例えば、パワースペクトルの分散値やＭＴＦ値など）を算出する。
【０１１０】
一方、前処理部２８において、事前に、画像（ＧＤ）（複数枚）に関する検査員（Ｕ）の目視による評価値に基づく判定基準（３６）を決定（設定）しておく。
【０１１１】
１００２は、上記２つの評価値（特徴量Ａ，Ｂ）に基づいて決定した判定基準（３６）の例を示す。特徴量Ａ，Ｂを軸とする２次元の特徴空間において、各レベル（Ｌ）の取り得る範囲が判定基準情報３６（例えば閾値の直線）として設定されている。
【０１１２】
本例では、IF  「ThA1＜Ａ≦ThA2」  AND  「ThB1＜Ｂ≦ThB2」  THEN  Level１、といったように、閾値（ThA1等）でレベルＬ（例えばLevel１，Level２，……等）の範囲が決定可能となるように、各レベルＬの範囲を、各軸に平行になるように矩形で区切った判定基準（３６）を示している。そして、この各レベルＬの範囲が設定された特徴空間上で、デジタル画像（ＧＤ）１０００から算出した評価値（特徴量Ａ，Ｂ）がどこに位置するかによって、１０００のレベルＬを決定する（レベル判定部２７、図９のＳ９０２、図１０の１００３）。
【０１１３】
１００２では、デジタル画像（ＧＤ）１０００から算出した評価値（Ａ，Ｂ）が例えばLevel３の領域内に位置することを示しており、この場合、レベルＬがLevel３に決定される（１００３）。
【０１１４】
上記例では、判定基準（３６）を特徴軸に平行な直線で区切った矩形領域で示したが、特徴軸に対して傾きをもった直線で区切ることも可能である。また直線ではなく、特徴空間上で各評価値の領域を曲線で区切ることも可能である。また、領域の設定は、設定部１５の処理で、IF  「ThA1＜Ａ≦ThA2」  AND  「ThB1＜Ｂ≦ThB2」  THEN  Level１、といった論理における閾値（ThA1，ThA2，ThB1，ThB2等）を、ユーザ（Ｕ）の入力により受け付けて、判定基準情報３６として設定することもできる。
【０１１５】
本発明において品質評価値（Ｈ）のレベル（Ｌ）を決定するための判定基準（３６）は、上記例（図１０）の２次元の特徴空間での設定に限らず、より高次元の特徴空間においても設定可能である。例えば３種類の評価値（特徴量）を用いた３次元の特徴空間において設定してもよく、その場合、判定基準（３６）は、平面、曲面で設定可能である。さらにより高次元の特徴空間においては、判定基準（３６）は、超平面、超曲面となる。ただこの場合、ユーザ（Ｕ）が複数の特徴軸に対して閾値などを設定するのは困難である。この点、本発明（実施の形態）では、以下に説明するように、教師あり学習による判定基準（３６）の自動設定機能を有する。これにより判定基準（３６）の設定が容易である。
【０１１６】
なお本システムでのレベル（Ｌ）判定の機能は、上記構成例のように検査員（Ｕ）の目視判定の評価値（それに基づく判定基準情報３６）と整合させたレベルＬを判定する場合と、検査員（Ｕ）の目視判定の評価値とは独立で自動的に数値計算でレベルＬを判定する場合（判定基準情報３６に検査員（Ｕ）の目視判定の評価値を反映しない場合）と、のいずれも可能である（設定部１５で選択可能）。
【０１１７】
［前処理、学習処理］
図１１は、前処理部２８による前述のレベル（Ｌ）判定のための判定基準（３６）を設定する処理例を示す。特に、学習部２９を用いて、検査員（Ｕ）であるユーザの観点（目視判定の官能評価値）を反映するように判定基準（３６）を学習処理で自動的に生成する処理例を示す。
【０１１８】
１１００は、対象のアナログ像（フィルム像）（例えばＡ１）であり、１１０１は、それ（Ａ１）をデジタル化して得られたデジタル画像（Ｄ１）のデータ（ＧＤ）である。１１００と１１０１の組は複数を使用する必要がある。
【０１１９】
Ｓ１１２は、検査員（Ｕ）が目視検査部１０５でアナログ像１１００に対するデジタル像１１０１のレベル（Ｌ）を、目視判定により得る処理（手順）を示す。Ｓ１１２の当該デジタル画像１１０１とその目視判定によるレベル値とのセットのデータ情報を、学習部２９の教師あり学習処理での教師情報３９として使用する。
【０１２０】
まず、レベル（Ｌ）が既知である画像、即ち教師付きのデジタル画像（１１０１）から、前述の通り（図８のＳ２０等と同様）、Ｈ計算部１３（２６，２５）により１つ以上の品質評価値（特徴量）１１０３を算出する（Ｓ１１１）。この処理を、前述の各レベル（Ｌ）を網羅するような複数のセットで繰返し、複数かつ様々なレベルに対応する品質評価値（特徴量）を算出する。
【０１２１】
そして、当該デジタル画像から算出された各品質評価値（特徴量）１１０３と、それに対応するデジタル画像の目視判定によるレベル値とのセット（教師情報３９）を、学習部２９において学習させる処理（Ｓ１１４）により、判定基準情報３６を生成・出力する。この学習は、回帰式（回帰分析）や、パラメトリックな手法、ノンパラメトリックな手法など、各種の一般的な識別器を使用することができる。その一例としては、サポートベクターマシン（ＳＶＭ）がある。
【０１２２】
また、上記学習処理（Ｓ１１４）は、非学習の単純な設定とすることも可能である。例えば前処理部２８／設定部１５からユーザ（Ｕ）入力により判定基準情報３６を設定可能である。
【０１２３】
以上の通り、本システムでは、各デジタル画像（ＧＤ）のレベル（Ｌ）を、検査員（Ｕ）などのユーザの目視判定の評価結果を教師情報３９として学習処理し、判定基準（３６）を生成し設定する。これにより、本システムで自動的に計算されるレベル（Ｌ）について、実際の検査員（Ｕ）による目視評価結果と整合のとれた値を得ることができる。
【０１２４】
［検査処理］
図１２は、実施の形態１における検査部１７による欠陥検出処理を示す。検査対象画像としては、上記評価対象としたデジタル画像データ（ＧＤ）、及び推定アナログ画像データ（ＧＡ）のいずれを用いることも可能である。特に推定アナログ画像データ（ＧＡ）を対象として高感度化検査が可能である。検査員（Ｕ）は、本検査装置１の検査ＳＷ１０を用いて、表示画面で情報を参照しながら、検査を実行（命令等）することができ、例えば検査の結果の検出欠陥情報３８（欠陥領域など）を画面で確認することができる。
【０１２５】
まずＳ１２０で、入力部１１などにより検査対象の画像、ここではまずデジタル画像データ（ＧＤ）を取得する。この段階で当該画像（ＧＤ）は理想像（アナログ像）に対して劣化がある。このためＳ１２１で、前記Ｓ１２，Ｓ１４のように画像復元（アナログ像推定）を行って理想像（アナログ像）を推定し、推定アナログ画像データ（ＧＡ）を得る。そしてＳ１２２で、推定アナログ画像データ（ＧＡ）を対象として、検査部１７による検査処理（高感度化検査）を実行する。検査処理の内容としては、公知の画像解析処理などにより欠陥部に相当する領域や座標を検出し、結果を検出欠陥情報３８として記憶部に記憶する。またＳ１２３で、出力部１７により、対象画像とともに検出欠陥情報３８を検査装置１の画面に表示する。
【０１２６】
［効果等（１）］
以上、実施の形態１により、検査員（Ｕ）が対象物のフィルム像もしくは対象物そのものを直接的に目視で観察するのと同様の画像をアナログ像推定処理（１２）により取得し、かつ、その推定アナログ画像（ＧＡ）に対するデジタル画像（ＧＤ）の品質・劣化度合いについて、検査員（Ｕ）の目視評価結果を既知の評価値として反映した情報（３９）から生成した判定基準（３６）を用いて画像の品質レベル（Ｌ）を判定する。これにより、検査員（Ｕ）の官能評価値と整合のとれたレベル値（Ｌ）を得ることができる。また上記推定アナログ画像（ＧＡ）を用いることで、高感度な検査が可能である。また、目視感と整合させるための判定基準（３６）の設定手段を備えることにより、最終的な品質の事前予測が可能となる。
【０１２７】
上記第１方式として、デジタル画像の品質を、アナログ像（ＧＡ）の推定とそれに対する劣化度によって判定する構成例を示した。これに対し、以下に示すように、本発明では、実施の形態２（第２方式）として、アナログ像（ＧＡ）を推定せずに品質レベルを判定することを可能とする。
【０１２８】
＜実施の形態２＞
図１３を用いて、本発明の実施の形態２のシステム、方法等について説明する。基本構成は実施の形態１と同様（図１等）であり、異なる要素として、第２方式処理部１８での処理を行う。第２方式では、デジタル画像（ＧＤ）のレベル（Ｌ）判定の別の例を示すものである。
【０１２９】
［第２方式処理］
図１３は、第２方式処理部１８による品質評価値（Ｈ）及びレベル（Ｌ）決定処理を示す。１３０１は、評価対象となるデジタル画像（ＧＤ）である（入力部１１から入力処理）。Ｓ１３１では、第２方式処理部１８は、例えば特徴抽出部２５を用いて、デジタル画像（ＧＤ）１３０１から、品質を決めるための特徴量を抽出する。画質を表す特徴量の例としては、解像度（ＭＴＦ）、コントラスト、ＳＮ比、ウィナースペクトル（ノイズ周波数）などがあり、ここではそれらを用いるが、それ以外でもよい。Ｓ１３１で例えばそれら複数の種類の特徴量Ａ，Ｂ，Ｃ……を抽出するとする。
【０１３０】
Ｓ１３２では、Ｓ１３１で算出した１つ以上の特徴量に対し、判定基準（３６）に基づいてレベル（Ｌ）を決定する。前述（図１０）の通り、複数の特徴量を用いて多次元の特徴空間における判定基準（３６）に基づいてレベル（Ｌ）を決定することも可能であるが、ここ（Ｓ１３２）では、複数の特徴量Ａ，Ｂ，Ｃ……から１つの統合品質評価値（統合指標値）（Ｈ０）を算出して一次元でレベル判定する処理例を示す。この算出の一例としては、下記の（式７）の通り、各特徴量（f1,f2,・・・）にそれぞれ重み（a,b,・・・）を付けた線形和によって統合品質評価値（Ｈ０）を算出する。
【０１３１】
【数７】
0005624975.tif 000008
【０１３２】
式７で、重み（a,b,・・・）は、図１１で説明した通り、学習機能（学習部２９の学習処理、例えば回帰式）などを用いて事前に決定・設定する。もしくは、設定部１５を用いて、ユーザ（Ｕ）からの設定情報を受け付けることも可能である。
【０１３３】
そしてＳ１３３で、第２方式処理部１８は、事前に設定済みの判定基準（３６）に基づき、上記の統合品質評価値（Ｈ０）から品質レベル（Ｌ）を決定する。１３０２は、各デジタル画像（ＧＤ）について、横軸にＳ１３２の統合品質評価値（Ｈ０）、縦軸にＳ１３３の品質レベル（Ｌ）をとった散布図である。散布図中のステップ状の線が、事前に設定された判定基準（３６）である。
【０１３４】
上記判定基準（３６）は、例えば、事前に、デジタル画像（ＧＤ）からＳ１３１，１３２により算出される値とそれに対応する検査員（Ｕ）の目視判定の評価結果（デジタル画像に対するレベル）を教師情報（３９）として与えて、１３０２のような散布図を作り、ＳＶＭなどの識別器を用いた機械学習により算出する。
【０１３５】
上記処理例では、画質を表す特徴量として解像度（ＭＴＦ）などを用いる例を述べたが、図８で示したように、推定アナログ画像（ＧＡ）とのパワースペクトルの比較から算出される特徴（比較値）も加えて統合的に判定することも可能である。つまり、品質評価値（Ｈ）は、評価対象となるデジタル画像（ＧＡ）のみから算出するもの、推定アナログ画像（ＧＡ）との比較で算出するもの等、様々に可能であり、それら１つ以上の品質評価値（ないし特徴量）を用いてデジタル画像の品質を保証するためのレベル（Ｌ）を判定すること等が本発明の範囲である。
【０１３６】
［効果等（２）］
以上、実施の形態２により、実施の形態１と同様の効果が得られる。
【０１３７】
＜実施の形態３＞
図１４を用いて、本発明の実施の形態３のシステム、方法等について説明する。基本構成は実施の形態１と同様（図１等）であり、異なる要素として、デジタイザ保証部１９での処理を行う。
【０１３８】
［デジタイザ保証］
図１４は、前述（実施の形態１）の画像品質評価と同様の処理に基づいて、デジタイザ１０３の性能の保証を行う場合の処理の流れを示す。
【０１３９】
例えば各種のフィルム像（Ａ１）をフィルムスキャナ１０４などのデジタイザを用いてデジタル画像化する際、そのデジタイザの性能が問題ないかどうかを保証する必要がある場合がある。本発明（実施の形態３）では、特定の評価用のフィルム像（評価ゲージ）（対応する評価パターンデータ３７）を用いて、デジタイザの性能の保証のための情報を計算・出力する処理を行う。以下では、図１のデジタイザ１０３（イメージングプレート）を対象とした例で説明する。
【０１４０】
１４０は、評価パターンデータ３７（評価ゲージ）である画像データの一例である。この画像データは、例えば、モノクロで規則的に各角度の線が含まれている。これ以外でも様々な周波数を持つパターン（評価しやすい情報）が含まれていればよい。
【０１４１】
まずＳ１４１で、評価対象となるデジタイザ１０３による評価パターンデータ３７（１３０）の入力に対応したデジタル画像（ＧＤ）を取得する。ここから、Ｓ１４２では、図１のＳ１２，Ｓ１４と同様の処理により、画像復元（アナログ像推定）を行い、推定アナログ画像データ(ＧＡ)を得る。Ｓ１４３で、推定アナログ像（ＧＡ）とデジタル画像（ＧＤ）とを比較し、デジタイザ１０３による、アナログ像に対するデジタル画像の劣化度（前述のＨ，Ｌ相当）を推定・評価する。そしてＳ１４４のデジタイザ保証処理では、上記の劣化度から、デジタイザ１０３から取得されたデジタル画像（ＧＤ）が検査に耐えうる品質を保っているか否かをレベル判定する。または、解像度Ｘ％劣化といった定量値を算出することで性能を評価する。
【０１４２】
上記Ｓ１４３で劣化度を推定するための基準としては、Ｓ１４５で、評価パターンデータ３７のデジタル画像とそれを検査員（Ｕ）による目視で照合・判定することで設定した劣化度とのセットを事前に入手して機械学習などによって当該基準情報を設定する（前述の前処理に対応する）。
【０１４３】
［効果等（３）］
以上、実施の形態３では、更に、品質評価値（Ｈ）、レベル（Ｌ）に相当する性能指標値を用いて、デジタイザ１０３の性能を保証することが可能となる。これにより、デジタイザ１０３の管理、機差の調整などが実現できる。
【０１４４】
以上、本発明者によってなされた発明を実施の形態に基づき具体的に説明したが、本発明は前記実施の形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能であることは言うまでもない。例えば、上述した放射線透過によるフィルム像より取得したデジタル画像の品質評価の形態に限らず、液晶パネルや有機ＥＬパネルなどのフラットパネルの画質検査にも同様に適用可能である。
【０１４５】
また、第１方式のようにデジタル画像から推定したアナログ画像と、第２方式のようにアナログ像を推定しない場合のデータとを組み合わせて品質を判定する形態なども可能である。
【産業上の利用可能性】
【０１４６】
本発明は、放射線透過試験などを含む各種のデジタル画像の検査・診断等に利用可能である。
【符号の説明】
【０１４７】
１…ＰＣ（検査装置）、１０…検査ＳＷ（ソフトウェア）、１１…入力部、１２…アナログ像推定部、１３…評価値（Ｈ）計算部、１４…レベル（Ｌ）計算部、１５…設定部、１６…出力部（表示部）、１７…検査部（欠陥検出部）、１８…第２方式処理部、１９…デジタイザ保証部、２１…ノイズ除去部、２２…逆変換部、２３…ＦＦＴ部、２４…比較部、２５…特徴抽出部、２６…評価値（Ｈ）定量化部、２７…レベル（Ｌ）判定部、２８…前処理部、２９…学習部、３１…デジタル画像データ（ＧＤ）、３２…推定アナログ画像データ（ＧＡ）、３３…周波数スペクトルデータ、３４…品質評価値（Ｈ）、３５…品質レベル値（Ｌ）、３６…判定基準情報、３７…評価パターンデータ、３８…検出欠陥情報、３９…教師情報、１０１…放射線検査装置（放射線検出装置）、１０２…現像装置、１０３…デジタイザ（イメージングプレート）、１０４…フィルムスキャナ、１０５…目視検査部。
【図１】
0005624975.tif 000009
【図９】
0005624975.tif 000010
【図１５】
0005624975.tif 000011
【図２】
0005624975.tif 000012
【図３】
0005624975.tif 000013
【図４】
0005624975.tif 000014
【図５】
0005624975.tif 000015
【図６】
0005624975.tif 000016
【図７】
0005624975.tif 000017
【図８】
0005624975.tif 000018
【図１０】
0005624975.tif 000019
【図１１】
0005624975.tif 000020
【図１２】
0005624975.tif 000021
【図１３】
0005624975.tif 000022
【図１４】
0005624975.tif 000023