JP5186833B2 - 画像品質管理支援方法 - Google Patents

Description

本発明は、画像品質管理支援方法に関する。

近年、医療分野において、患者の病状の診断のための画像を得るために放射線画像読取システムが用いられている。放射線画像読取システムでは、放射線撮影により被写体を透過した放射線エネルギーの一部を輝尽性蛍光体シートの内部に一旦蓄積し、放射線エネルギーが蓄積された輝尽性蛍光体シートを励起光で２次元的に走査して輝尽発光を生じさせ、この輝尽発光光を光電的に読み取って画像信号を得て、これをＡ／Ｄ変換してデジタルの画像データを取得する。

このような放射線画像読取システムにおいて得られた画像の画質が悪いと、適正な診断を行うことができなくなる。そのため、オペレータやサービスマンによって、画像読取装置やカセッテのＱＡ（Quality Assurance:品質保証）やＱＣ（Quality Control：品質管理）を行なうための検査方法がある。

この検査方法としては、放射線透過性の基板に放射線を遮蔽する部材や吸収する部材が配置されたＱＣファントムを撮影してＱＣファントムの放射線画像を得て、この放射線画像を検査対象となる画像読取装置で読み取り、読み取った放射線画像の画像データに品質評価結果を取得するための演算（以下、ＱＣ演算という）を施す検査が知られている。

検査精度の信頼性を確保するためには、ＱＣ演算時にＱＣファントムに内在された各種部材の位置を特定とする必要がある。従って、検査対象領域に対応する位置にＱＣファントムを設置しなくてはならず、また、ＱＣファントムの位置決めや位置ズレの補正を行なう技術が開示されている。

例えば、Ｘ線システムの評価を行なうシステムにおいて、Ｘ線ファントムを使用して、Ｘ線システム・パラメータの測定を行なう場所である関心領域（ＲＯＩ）を予測する場合、Ｘ線ファントムに設けられた位置決め用タブを用いて、Ｘ線ファントムの位置決めを行なう方法が開示されている（特許文献１参照）。

また、画質測定パターンの画像の位置を検出するために用いられるマーカが複数配置されたファントムを用いて撮影を行い、撮影で得られた画像データに基づいて、マーカ画像の位置を検出し、検出されたマーカ画像の位置と、予め記憶されているファントム画像におけるマーカ画像の基準位置とを比較することにより、ファントム画像の基準位置からの平行方向及び回転方向のズレ（位置ズレ）量を算出し、当該位置ズレ量に基づいてファントム画像の位置ズレ量を補正する技術が開示されている（特許文献２参照）。

更に、検査精度の信頼性を確保するために、画質劣化を早期に認識して、画質劣化の原因の対応を指示する必要がある。
そこで、端末より送信されたファントムを撮影して得られた画像及び当該画像の解析結果である画質パラメータを記憶する画質パラーメータデータベースに記憶し、画質パラメータが所定の基準値内に含まれるか否かの判定結果のレポートを端末に送信する画質管理システムが開示されている（特許文献３参照）。
特開２００１−５４５１４号公報 特開２００４−２９８６１７号公報 特開２００５−２０５１９３号公報

しかしながら、ＱＣファントムは、１施設に１つ以上設置されることが一般的であるが、多数のカセッテ型放射線撮影装置（ＣＲ）及び読取装置を設置している施設では、定期的に実施される画質評価において、放射線撮影装置と読取装置との組み合わせに相当する回数のＱＣファントムの撮影が必要となり、ＱＣファントムの使用回数が非常に多くなる傾向にある。このＱＣファントムの撮影における設置作業は、技師の手作業により行なわれており、当該設置作業においてＱＣファントムが破損したり、傾いて設置されたり、放射線源に対向する規定面に設置されていないと、位置合せが良好に行なわれなくなることが想定される。

特許文献１及び２のように、ＱＣファントムの撮影時に位置ズレが生じないように、位置決めタブが設けたり、撮影時の位置ズレを補正するために位置合せマークが付与されたりしたとしても、撮影後になっては、ＱＣファントムが破損していたり、不良な設置であったりする等のＱＣファントム自体や取り扱いに起因するものであるかの識別が困難である。特に、特許文献２のように、ＱＣファントムの画像に対して補正を行なう場合、補正が行なわれた画像で画質評価が行なわれることとなり、ＱＣファントムの取り扱いに起因する評価不良の識別が遅れる恐れがある。また、特許文献３のように、画質パラメータの評価結果に基づいて、画質の劣化の要因を早期に発見する方法では、ＱＣファントム自体や取り扱いに起因する劣化であるかの識別ができない。

本発明の課題は、上記問題に鑑みて、ＱＣファントムの取り扱いに起因する画質評価の不良を発見でき、ＱＣファントムの取り扱いに対する対策を講じさせ得ることである。

請求項１に記載の発明によれば、左右非対称及び上下非対称の形状を有する設置判別部材と、画質評価項目に対応した大きさ又は形状を有する評価用部材とが内在されたＱＣファントムが撮影された放射線画像を取得する画像取得工程と、前記ＱＣファントムを個別的に識別する識別情報を取得する識別情報取得工程と、前記画像取得工程により取得された放射線画像に基づいて前記設置判別部材の位置を特定する位置特定工程と、前記位置特定工程により特定された前記設置判別部材の位置に基づいて、前記ＱＣファントムの設置情報を算出する設置情報算出工程と、前記画像取得工程により取得された放射線画像に基づいて前記評価用部材に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する評価演算工程と、前記評価演算工程により生成された評価結果情報と、前記識別情報取得工程により取得された前記識別情報と、前記設置情報算出工程により算出された前記設置情報とが対応付けられた画像品質評価情報を記憶する記憶工程と、を含む画像品質管理支援方法が提供される。

また、前記記憶工程により記憶された複数の前記画像品質評価情報を読み出す読出工程と、読み出された前記複数の画像品質評価情報に含まれる前記設置情報に基づいて、前記ＱＣファントムに対する取り扱いの注意喚起情報を報知するか否かを判定する報知判定工程と、前記報知判定工程により前記注意喚起情報を報知すると判定された場合、当該注意喚起情報を報知する報知工程と、を含む。

さらに、前記設置情報は、前記ＱＣファントムの撮影時における前記ＱＣファントムの基準位置からの傾き量、又は、前記ＱＣファントムの撮影時において前記ＱＣファントムが照射線源に対向する規定面と対向する面に設置されて撮影されたことを示す非規定面情報を含む。

さらにまた、前記報知判定工程は、所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、且つ、前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、前記ＱＣファントムの撮影時の取り扱い手順の見直しを指示する前記注意喚起情報を報知すると判定する。

請求項２に記載の発明によれば、請求項１に記載の画像品質管理支援方法において、前記報知判定工程は、所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、且つ、前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数未満、及び、前記非規定面情報を含む画像品質評価情報の数が予め設定された数未満の場合、前記ＱＣファントムの破損確認を指示する前記注意喚起情報を報知すると判定する。
請求項３に記載の発明によれば、左右非対称及び上下非対称の形状を有する設置判別部材と、画質評価項目に対応した大きさ又は形状を有する評価用部材とが内在されたＱＣファントムが撮影された放射線画像を取得する画像取得工程と、前記ＱＣファントムを個別的に識別する識別情報を取得する識別情報取得工程と、前記画像取得工程により取得された放射線画像に基づいて前記設置判別部材の位置を特定する位置特定工程と、前記位置特定工程により特定された前記設置判別部材の位置に基づいて、前記ＱＣファントムの設置情報を算出する設置情報算出工程と、前記画像取得工程により取得された放射線画像に基づいて前記評価用部材に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する評価演算工程と、前記評価演算工程により生成された評価結果情報と、前記識別情報取得工程により取得された前記識別情報と、前記設置情報算出工程により算出された前記設置情報とが対応付けられた画像品質評価情報を記憶する記憶工程と、を含む画像品質管理支援方法が提供される。
また、前記記憶工程により記憶された複数の前記画像品質評価情報を読み出す読出工程と、読み出された前記複数の画像品質評価情報に含まれる前記設置情報に基づいて、前記ＱＣファントムに対する取り扱いの注意喚起情報を報知するか否かを判定する報知判定工程と、前記報知判定工程により前記注意喚起情報を報知すると判定された場合、当該注意喚起情報を報知する報知工程と、を含む。
さらに、前記設置情報は、前記ＱＣファントムの撮影時における前記ＱＣファントムの基準位置からの傾き量、又は、前記ＱＣファントムの撮影時において前記ＱＣファントムが照射線源に対向する規定面と対向する面に設置されて撮影されたことを示す非規定面情報を含む。
さらにまた、前記報知判定工程は、所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、且つ、前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数未満、及び、前記非規定面情報を含む画像品質評価情報の数が予め設定された数未満の場合、前記ＱＣファントムの破損確認を指示する前記注意喚起情報を報知すると判定する。

請求項４に記載の発明によれば、請求項１から３のいずれか一項に記載の画像品質管理支援方法において、前記ＱＣファントムは、自己の識別情報を示す識別部材を内在しており、前記識別情報取得工程は、画像取得工程により取得された放射線画像に基づいて前記識別部材が示す識別情報を取得する。

本発明によれば、放射線画像の被写体であるＱＣファントムを識別する識別情報と当該ＱＣファントムの設置情報とを当該放射線画像の評価結果情報に対応付けて画像品質評価情報として記憶することができるため、各ＱＣファントムの撮影時における設置傾向を把握することが可能となり、ＱＣファントムの取り扱いに起因する画質評価の不良において、ＱＣファントムの取り扱いに対する対策を講じさせ得ることができる。

以下、図を参照して本発明の実施の形態を詳細に説明する。
まず、構成を説明する。
図１に、放射線画像読取システムＮＳの全体構成図を示す。
図１に示すように、放射線画像読取システムＮＳは、撮影装置１０と、読取装置２０と、コンソール３０と、イメージャ４０と、サービスセンタ５０と、を備えて構成される。読取装置２０、コンソール３０及びイメージャ４０は、通信ネットワークＮ１を介して互いにデータ送受信可能に接続され、また、コンソール３０とサービスセンタ５０とは、通信ネットワークＮ２を介して互いにデータ送受信可能に接続されている。なお、各装置の台数は、特に限定されない。

通信ネットワークＮ１、Ｎ２は、ＬＡＮ（Local Area Network）やＷＡＮ（Wide Area Network）、インターネット等の様々な回線形態を適用可能である。なお、病院等の医療機関内で許可されるのであれば、無線通信や赤外線通信であってもよいが、重要な患者情報を含むため、送受信される情報は暗号化することが好ましい。また、病院内の通信方式としては、一般的に、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格が用いられており、上述した通信ネットワークＮ１、Ｎ２上の各装置間の通信では、ＤＩＣＯＭ ＭＷＭ(Modality Worklist Management)やＤＩＣＯＭ ＭＰＰＳ（Modality Performed Procedure Step）が用いられる。

以下、放射線画像読取システムＮＳの各構成装置について説明する。
撮影装置１０は、放射線照射スイッチ（図示せず）や、Ｘ線源などの放射線源１１等を有し、被写体Ｍに放射線を照射して撮影を行う。撮影装置１０において撮影を行う際には、被写体Ｍの放射線源１１と反対方向に、輝尽性蛍光体シートを収容したカセッテ１２が配置される。

カセッテ１２は、放射線エネルギーを蓄積する輝尽性蛍光体シートを備える放射線画像変換シートが収容された持ち運び可能な放射線変換媒体であり、撮影装置１０の放射線源１１からの照射放射線量に対する被写体Ｍの放射線透過率分布に従った放射線量を内設された輝尽性蛍光体シートの輝尽性蛍光体層に蓄積し、この輝尽性蛍光体層に被写体Ｍの放射線画像情報を記録する。また、カセッテ１２表面には、カセッテ１２を個々に識別するためのカセッテＩＤ（カセッテ識別情報）が貼付されている。

読取装置２０は、輝尽性蛍光体シートを収容したカセッテ１２が装填されると、装填されたカセッテ１２に記録された放射線画像を読み取る装置であり、プログラムに従って装置全体を制御するとともに各種処理を実行するＣＰＵ（Central Processing Unit）、各種プログラムやデータを記憶するＲＯＭ(Read Only Memory)及びＲＡＭ（Random Access Memory）等からなる制御部、操作部、装填されたカセッテ１２に収容された或いは内蔵された輝尽性蛍光体シートに記録された被写体Ｍの放射線画像を読み取って画像データを得る読取部、読取部による読取後に輝尽性蛍光体上の残像を消去する消去部、装填されたカセッテ１２に貼付されたバーコードを読み取るバーコードリーダ、通信ネットワークＮ１上の各装置とのデータ送受信を制御する通信部等を備えて構成される。輝尽性蛍光体シートは、放射線画像を記録する記録媒体である。

読取部は、撮影により被写体Ｍを透過した放射線エネルギーが蓄積された輝尽性蛍光体シートを励起光で２次元的に走査して輝尽発光を生じさせ、この輝尽発光光を光電的に読み取って画像信号を得、これをＡ／Ｄ変換してデジタルの画像データを取得するものである。例えば、レーザービームを励起光として用い、レーザービームを主走査方向に偏向して輝尽性蛍光体シートに入射するとともに、輝尽性蛍光体シートを副走査方向に走査搬送することにより励起光で輝尽性蛍光体シートを二次元的に走査し、生じた輝尽性発光光を光伝播器で伝播して、ＣＣＤ（Charge Coupled Device）やフォトマルチプライヤ等の光検出器で検出することにより、輝尽性蛍光体に記録された放射線画像を読み取る。

コンソール３０は、撮影メニュー（品質評価処理用の撮影を行うＱＣメニューと、患者を撮影する一般撮影メニューとがある）の選択、撮影に使用するカセッテ１２等の登録、一般撮影における撮影対象となる患者の患者情報や撮影部位等の撮影オーダ情報の登録、読取装置２０で読み取られた画像データと撮影オーダ情報との対応付け及び保存、画像データに対する画像処理、画像データのイメージャ４０への送信及び出力制御等を行い、更に、本実施の形態においては、品質評価処理を行なった際に生成される各種レポートをサービスセンタ５０に送信するための装置である。

図２に、コンソール３０の機能構成例を示す。図２に示すように、コンソール３０は、ＣＰＵ３１、操作部３２、表示部３３、第１通信部３４、第２通信部３５、バーコードリーダ３６、ＲＡＭ３７、記憶部３８、ＨＤＤ（Hard Disk Drive）３９等を備えて構成され、各部はバスＮ３により接続されている。

ＣＰＵ３１は、記憶部３８に記憶されている各種プログラムを読み出してＲＡＭ３７内に形成されたワークエリアに展開し、該プログラムに従った処理を実行することにより各部の動作を集中制御する。

操作部３２は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号をＣＰＵ３１に出力する。また、操作部３２は、表示部３３の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号をＣＰＵ３１に出力する。

表示部３３は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ等のモニタにより構成され、ＣＰＵ３１から入力される表示信号の指示に従って、操作部３２からの入力指示やデータ等を表示する。また表示部３３は、ＣＰＵ３１から入力される表示信号の指示に従って、後述する品質評価処理において生成された各種レポートを画面上に出力する。

第１通信部３４は、ＬＡＮアダプタやルータやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮ１に接続された各装置との間の通信を制御する。即ち、第１通信部３４は、通信ネットワークＮ１を介して読取装置２０で得られた画像データが入力され、さらに読取装置２０の装置ＩＤが通信ネットワークＮ１を介して入力される。

第２通信部３５は、ＬＡＮアダプタやルータやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮ２に接続されたサービスセンタ５０との間の通信を制御する。即ち、第２通信部３５は、生成した各種レポートがサービスセンタ５０に出力され、また、サービスセンタ５０に蓄積され各種レポートが通信ネットワークＮ２を介して入力される。

バーコードリーダ３６は、カセッテ１２上に貼付されたバーコードの表すカセッテＩＤを読み取ってＣＰＵ３１に出力する。

ＲＡＭ３７は、ＣＰＵ３１により実行制御される各種処理において、記憶部３８から読み出されたＣＰＵ３１で実行可能な各種プログラム、入力若しくは出力データ、及びパラメータ等の一時的に記憶するワークエリアを形成する。

記憶部３８は、不揮発性の半導体メモリ等により構成され、ＣＰＵ３１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、ＣＰＵ３１は、当該プログラムコードに従った動作を逐次実行する。

ＨＤＤ（Hard Disc Drive）３９は、読取装置２０から送信された一般撮影の画像データを撮影オーダ情報と対応付けて記憶し保存する。また、ＨＤＤ３９は、読取装置２０から送信されたＱＣファントムを撮影して得られた放射線画像の画像データに、当該画像データから取得したＱＣ識別情報、読み取りを行った装置の装置ＩＤ、撮影に使用したカセッテＩＤ、読取時刻（日付を含む）情報等を対応付けて記憶し保存する。
ＱＣ識別情報とは、ＱＣファントムを個別的に識別する識別情報である。

また、ＨＤＤ３９は、ＱＣ撮影管理情報ファイル３９１、一般撮影管理情報ファイル３９２、画像ファイル記憶領域３９３を有している。
ＱＣ撮影管理情報ファイル３９１は、ＱＣファントムの画像データ毎の、管理ＮＯ、装置ＩＤ、カセッテＩＤ、ＱＣ識別情報、読取時刻情報、処理済みフラグ（デフォルトはＯＦＦ）、画像ＩＤを対応付けて１レコードとして記憶するファイルである。
一般撮影管理情報ファイル３９２は、登録された撮影オーダ情報と当該撮影オーダ情報に従って撮影された画像データの画像ＩＤ等を対応付けて記憶するファイルである。
画像ファイル記憶領域３９３は、読取装置２０から受信した画像データを記憶、保存するために設けられた領域である。

なお、読取装置２０から送信される画像データは、ＤＩＣＯＭ形式の画像ファイルであり、画像ＩＤを含む付帯情報が付帯されている。本実施の形態において、画像データは、この画像ＩＤにより、ＱＣ撮影管理情報ファイル３９１に保存されている当該画像データに関する情報や、一般撮影管理情報ファイル３９２に記憶されている当該画像データに関する撮影オーダ情報と対応付けられている。

ＨＤＤ３９の記憶容量は、有る程度余裕を持って準備する必要がある。
通常、撮影が行われ、読取装置２０等から画像データが受信された後、遅くとも１週間以内には医師による読影が開始され、又は、後述の品質評価処理が実行される。従って、画像データが受信されてから少なくとも１週間は保存が可能な記憶容量としなければならない。

例えば、コンソール３０が取り扱う画像データが一画像につき２０００画素×２５００画素で、１日１００画像程度が生成されるとすると、容量４０ＧＢのＨＤＤ（画像保存用の記憶容量は２０ＧＢとする）であれば２０００画像を約２０日保存することができる。

放射線画像読取システムＮＳにおいては、少なくとも一台のコンソール３０に、ＱＣツール３０１が搭載されている。ＱＣツール３０１は、品質評価用の演算プログラムであり、ＱＣツール３０１が搭載されたコンソール３０のＣＰＵ３１は、ＱＣツール３０１との協働によるソフトウエア処理により、品質評価処理を実現する。また、ＱＣツール３０１を搭載したコンソール３０は、品質評価処理により生成される各種レポートを記憶するＱＣ評価結果ファイル３９４を有している。以下、ＱＣツール３０１が搭載されたコンソール３０をコンソール３０ａとして説明する。

ＱＣツール３０１により品質評価処理を実行するコンソール３０ａは、評価部材と設置判別部材と識別部材とが内在されたＱＣファントム（試験用模擬被写体模型）が撮影された放射線画像を読取装置２０から取得し、当該放射線画像に対応するＱＣ識別情報を取得し、取得した放射線画像に基づいて設置判別部材の位置を特定し、特定した設置判別部材の位置に基づいて、ＱＣファントムの設置情報を算出する。また、取得した放射線画像に基づいて評価部材毎に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する。そして、評価結果情報、ＱＣ識別情報、設置情報を対応付けた画像品質評価情報（評価レポート）を作成し、サービスセンタ５０に送信して記憶させる。

図３に、ＱＣファントムＰの一例を示す。図３に示すＱＣファントムＰは、放射線源１１に対向する規定面（以下、表面）から内部を透過して見た図である。
ＱＣファントムＰは、２枚のアクリル板などの放射線透過性の基板の間に、互いに異なる大きさ又は形状を有し放射線を遮蔽する材料や吸収する材料によって読取装置の線形性、表示寸法性精度、鮮鋭性、読取レーザ走査性、低コントラスト分解能等の各画質評価項目用の評価部材が複数封入された品質評価用の被写体である。
また、本実施の形態のＱＣファントムＰは、設置判別部材６と識別部材７とが基板の間に封入されている。

評価部材としては、表示寸法精度評価用の金属円板１ａ〜１ｃ、鮮鋭性評価用のエッジパターン２、線形性評価用のステップウエッジ３ａ〜３ｄ、読取レーザ走査性評価用のジッタ４ａ、４ｂ、低コントラスト分解能評価用のバーガーファントム５ａ〜５ｊとする。
本実施の形態においては、金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂを第１の評価部材とし、第１の評価部材よりも低コントラストであるバーガーファントム５ａ〜５ｊを第２の評価部材とする。

金属円板１ａ〜１ｃは、円状の銅板で構成されており、図３に示すように、ＱＣファントムＰの三隅に配置されている。金属円板１ａ〜１ｃの画像データは、読取装置の表示寸法精度を定量評価するために用いられる。

エッジパターン２は、辺縁が鮮鋭なエッジ部が形成されたタングステン板で形成され、エッジ部を読取装置２０による読み取りの主走査方向及び副走査方向に対して、若干、傾けた状態で配置されている。エッジパターン２の画像データは、読み取り装置の鮮鋭度（ＭＴＦ；Modulation Transfer Function）を定量評価するために用いられる。

ステップウエッジ３ａ〜３ｄは、段階的に互いに厚さの異なる４つの矩形の銅板で校正されている。図３に示すステップウエッジ３ａ〜３ｄは、ステップウエッジ３ａからステップウエッジ３ｄに向かって、順次、厚さが薄くなるものであり、放射線透過量のステップウエッジとなっている。ステップウエッジ３ａ〜３ｄの画像データは、読取装置の線形性やダイナミックレンジを定量評価するために用いられる。

ジッタ４ａ、４ｂは、直線状のエッジ部が形成された長方形状の銅板で形成されており、ジッタ４ａ、４ｂのいずれか一方の長辺は、読取装置２０による読み取りの主走査方向Ｘに平行になるように配置され、他方の長辺は、副走査方向Ｙに平行になるように配置されて構成されている。ジッタ４ａ、４ｂの画像データは、各走査方向についての読取レーザの走査性能（スキャニング精度）を評価するために用いられる。

バーガーファントム５ａ〜５ｊは、所定の直径を有する互いに高さ（厚さ）が異なるアクリル製の円柱で構成されており、図３に示すように、例えば、時計周り方向に、順次、高さ（厚さ）が段階的に低く（薄く）なり、直径が小さくなるように並べて配置されている。バーガーファントム５ａ〜５ｊの画像は、読取装置の低コントラスト分解能を目視評価又は定量評価するために用いられる。

設置判別部材６は、左右非対称及び上下非対称の形状を有する銅板で形成されている。本実施の形態においては、設置判別部材６として、ＱＣファントムの表面において、主走査方向Ｘが副走査方向Ｙよりも短いＬ字状の銅板とする。

識別部材７は、識別部材７が封入されたＱＣファントムＰのＱＣ識別情報を示す銅板で形成されている。本実施の形態においては、識別部材７としてＱＣ識別情報を示す番号が放射線画像として表示される銅板である。

なお、本実施の形態において利用可能なＱＣファントムは、図３に示すＱＣファントムＰに限定されない。

イメージャ４０は、コンソール３０から受信された画像データに基づきフィルム等の記録媒体上に画像を記録し、可視像として再生したハードコピーを出力する装置である。イメージャ４０としては、例えば、画像データに基づいて、透過型記録媒体（フィルム）にレーザ露光することによって潜像を記録し、熱現象処理により潜像を可視化する光熱銀塩方式のイメージャや、インクジェット方式或いはレーザ方式で反射型記録媒体（紙媒体、シール等）に記録するイメージャ等を用いることができる。

サービスセンタ５０は、コンソール３０から通信ネットワークＮ２を介して入力される各種レポートを蓄積し、蓄積した複数のレポートに基づいてＱＣファントムの取扱判定処理を行う装置である。

図４に、サービスセンタ５０の機能構成例を示す。図４に示すように、サービスセンタ５０は、ＣＰＵ５１、操作部５２、表示部５３、通信部５４、ＲＡＭ５５、記憶部５６、ＨＤＤ（Hard Disk Drive）５７等を備えて構成され、各部はバスＮ４により接続されている。

ＣＰＵ５１は、記憶部５６に記憶されている各種プログラムを読み出してＲＡＭ５５内に形成されたワークエリアに展開し、該プログラムに従った処理を実行することにより各部の動作を集中制御する。

操作部５２は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号をＣＰＵ５１に出力する。また、操作部５２は、表示部５３の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号をＣＰＵ５１に出力する。

表示部５３は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ等のモニタにより構成され、ＣＰＵ５１から入力される表示信号の指示に従って、操作部５２からの入力指示やデータ等を表示する。また、本実施の形態において、表示部５３は、ＣＰＵ５１から入力される表示信号の指示に従って、後述する取扱判定処理において生成されたメッセージを画面上に出力する。

通信部５４は、ＬＡＮアダプタやルータやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮ２に接続された各コンソール３０との間の通信を制御する。即ち、通信部５４は、コンソール３０から送信される各種レポートが通信ネットワークＮ２を介して入力される。

ＲＡＭ５５は、ＣＰＵ５１により実行制御される各種処理において、記憶部５６から読み出されたＣＰＵ５１で実行可能な各種プログラム、入力若しくは出力データ、及びパラメータ等の一時的に記憶するワークエリアを形成する。

記憶部５６は、不揮発性の半導体メモリ等により構成され、ＣＰＵ５１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、ＣＰＵ５１は、当該プログラムコードに従った動作を逐次実行する。また、本実施の形態においては、取扱判定処理の実行に必要なプログラムやデータ等を記憶する。

ＨＤＤ（Hard Disc Drive）５７は、レポートファイル５７１を有している。
レポートファイル５７１は、コンソール３０から送信された各種レポートを施設毎に記憶し保存するファイルである。

次に、上記放射線画像読取システムＮＳにおいて、カセッテ１２及び読取装置２０の組み合わせを品質評価対象とした場合の品質評価処理に用いる画像データを取得する画像取得処理を説明する。

ここで、放射線画像読取システムＮＳにおいて得られる画像データの品質は、放射線画像を記録する輝尽性蛍光体シートやＦＰＤ等の記録媒体及び記録媒体から画像データを読み取る読取装置に依存する。即ち、撮影済みのカセッテ１２を読取装置２０に装填して読み取りを行う場合、得られる画像データの品質は使用するカセッテ１２及び読取装置２０の組み合わせによって決まるものである。

画像取得処理は、読取装置２０とコンソール３０とにより実行される。読取装置２０側処理は、読取装置２０のＣＰＵとＲＯＭに記憶されているプログラムとの協働によるソフトウエア処理により実現され、コンソール３０側処理は、ＣＰＵ３１と、記憶部３８に記憶されているプログラムとの協働によるソフトウエア処理により実現される処理である。

コンソール３０において、オペレータ（撮影技師等）による操作部３２からの操作に応じて表示部３３に撮影メニュー画面が表示され、当該撮影メニュー画面からＱＣメニューが選択される。そして、品質評価対象となるカセッテ１２のバーコードがバーコードリーダ３６から読み取られてカセッテＩＤが入力される。入力されたカセッテＩＤと品質評価処理用の撮影であることを示す情報が登録され、自己のコンソールＩＤ等が対応付けられ、通信ネットワークＮ１に接続された全ての読取装置２０に通知（送信）される。

読取装置２０において、通信部を介してコンソール３０からカセッテＩＤ等が通知されると、通知されたコンソールＩＤ、カセッテＩＤ等が対応付けられてＲＡＭに記憶される。オペレータにより、コンソール３０で登録された品質評価対象のカセッテ１２に冶具を付けたり、予め定められた位置に配置したりして、ＱＣファントムを予め定められた位置に固定し、放射線撮影が行われる。

品質評価対象とする何れかの読取装置２０にカセッテ１２が装填されると、読取装置２０のバーコードリーダにより装填されたカセッテ１２のカセッテＩＤが読み取られ、ＲＡＭにおいて読み取られたカセッテＩＤに対応付けて記憶されているコンソールＩＤ等が検索される。
そして、検索により得られた情報に従って、装填されたカセッテ１２に記録されたＱＣ用の画像の読み取りが行われ、得られた画像データに、自己の装置ＩＤ、読み取ったカセッテ１２のカセッテＩＤ、読取時刻情報等が付帯情報として付帯され、通信部を介してカセッテ１２を登録したコンソール３０（検索により得られたコンソールＩＤのコンソール３０）に送信される。

コンソール３０において、第１通信部３４を介して読取装置２０から画像データ及びその付帯情報が受信されると、受信された画像データにその付帯情報が対応付けて保存され、本処理は終了する。

なお、上述した画像取得処理において、カセッテ１２として校正用として定められたカセッテを用いることとすれば、読取装置２０を品質評価対象とし、読取装置２０の品質変動等を評価するための画像データを取得することができる。また、カセッテ１２を読み取る装置として校正用として定められた読取装置２０を用いることとすれば、カセッテ１２を品質評価対象とし、カセッテ１２の品質変動等を評価するための画像データを取得することができる。

図５、図６を参照して、ＱＣツール３０１が搭載されたコンソール３０ａにおいて実行される品質評価処理について説明する。図５、図６に、品質評価処理のフローチャートを示す。当該品質評価処理は、操作部３２の所定の操作により品質評価処理の実行が指示された際に、ＣＰＵ３１とＱＣツール３０１等との協働によるソフトウエア処理により実現される処理である。

まず、第１通信部３４を介して読取装置２０から撮影されたＱＣファントムの画像データが取得される（ステップＳ１）。
即ち、ステップＳ１は、設置判別部材６、識別部材７、各画質評価項目に対応した複数の評価部材が内在されたＱＣファントムが撮影された放射線画像としての画像データを取得する画像取得工程に相当する。

ステップＳ１において取得された画像データにＯＣＲ処理（光学式文字読取処理）が施され、識別部材７に対応する画像から文字情報が取得される。そして、取得された文字情報は、画像データに対応するＱＣファントムの識別情報として取得される（ステップＳ２）。
即ち、ステップＳ２は、取得された画像データに基づいて識別部材が示す識別情報（ＱＣ識別情報）を取得する識別情報取得工程に相当する。

ステップＳ１において取得された画像データから、品質評価対象となるＱＣファントムが写された領域の画像データ（以下、ＱＣ画像データ）が抽出される（ステップＳ３）。なお、ステップＳ３において抽出される領域は、ＱＣファントムが撮影時に所定位置に固定して撮影されることから、予め設定することができるものである。

抽出されたＱＣ画像データの信号値の変換が行われ、画像データからＱＣ画像データが抽出できたか否かが判別され（ステップＳ４）、ＱＣ画像データが抽出できない場合（ステップＳ４；Ｎｏ）、ＱＣ画像データが抽出できなかった旨を示すエラーレポート（異常発生情報）が生成され、当該エラーレポートがＱＣ評価結果ファイル３９４に記憶されると共に、第２通信部３５を介してサービスセンタ５０に送信され（ステップＳ５）、本処理は終了される。

ＱＣ画像データが抽出できた場合（ステップＳ４；Ｙｅｓ）、信号値の変換が行なわれたＱＣ画像データからヒストグラムが作成される（ステップＳ６）。作成されたヒストグラムに基づいて閾値が決定され（ステップＳ７）、決定された閾値を基準としてＱＣ画像データが二値化され、第１の二値化画像が生成される（ステップＳ８）。

即ち、各画素のピクセル値が閾値のピクセル値を超えているか否かに基づいて、各画素に２つのバイナリ・レベルのうち一方が割り当てられる。
例えば、閾値のピクセル値が１４５である場合、ピクセル値が２００（閾値よりも暗いピクセル値）の画素に対しては１が設定され、一方、ピクセル値が１００（閾値よりも明るいピクセル値）の画素に対しては０が設定される。これにより、ＱＣ画像データ内の画素を、設置判別部材、識別部材、評価部材が存在しない領域（１に設定された画素の領域）と、設置判別部材、識別部材、評価部材が存在する領域（０に設定された画素の領域）と、の２つの領域に分類することができる。

第１の二値化画像の０に設定された画素がカウントされ、０に設定された画素の総数、即ち、評価部材が存在する領域に分類された画素の総面積が算出される。この算出された総面積と、予め定められた金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７の平面状の総面積（基準面積）とが比較され、算出された総面積が基準面積の予め定められた範囲内であるか否かが判別される（ステップＳ９）。

ステップＳ８において、基準面積にバーガーファントム５ａ〜５ｊの平面状の総面積が加えられていないのは、バーガーファントム５ａ〜５ｊは、他の部材よりもコントラストが非常に低いため、バーガーファントム５ａ〜５ｊに対応する画素は、ステップＳ８おいてＱＣ画像データを二値化した際に１に設定、即ち、存在しないと判定されるからである。

算出された総面積は基準面積の所定範囲内でない場合（ステップＳ９；Ｎｏ）、ＱＣ画像データが適切でない等を示すエラーレポートが生成され、当該エラーレポートがＱＣ評価結果ファイル３９４に記憶されると共に、第２通信部３５を介してサービスセンタ５０に送信され（ステップＳ５）、本処理は終了される。

ステップＳ９において算出された総面積と基準面積とを比較することにより、二値化処理が成功したか否か、即ち、検出されるべき部材（金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７）の領域が検出されたか否かが判別できる。例えば、算出された総面積が基準面積の所定範囲内よりも小さい場合、検出されるべき部材のいずれかが検出されていないと判別され、ＱＣファントムの撮影時の設置不良やＱＣファントムの破損が発生していると判断される。また、算出された総面積が基準面積の所定範囲内よりも大きい場合、ＱＣ画像データに異物が写りこんでいたりＱＣファントムの破損が発生していたりすることが判断される。このように、算出された総面積が基準面積の所定範囲外である場合には、ＱＣ画像データは正常な品質評価を行うことができるデータではないため、以後の処理は実行されず、当該判断結果がエラーレポートとして生成される。

算出された総面積は基準面積の所定範囲内である場合（ステップＳ９；Ｙｅｓ）、第１の二値化画像に点在する０に設定された画素が隣接する画素と連結した複数の領域を第１の二値化領域として抽出してラベリングを行なうラベリング処理が実行される（ステップＳ１０）。

図７に、ラベリング処理のフローチャートを示す。当該ラベリング処理は、ＣＰＵ３１とＱＣツール３０１等との協働によるソフトウエア処理により実現される処理である。
ラベリング処理は、まず、第１の二値化画像が走査され、ラベル（例えば、各第１の二値化領域を識別する番号）が付されていない０に設定された画素を検出し（ステップＳ３１；Ｙｅｓ）、当該画素に対してラベルを付す（ステップＳ３２）。そして、当該ラベルを付した画素に隣接する画素のうち０に設定されている画素があるか否かが判別され（ステップＳ３３）、０に設定されている隣接する画素がある場合（ステップＳ３３；Ｙｅｓ）、当該隣接する画素に対して同一のラベルを付す（ステップＳ３４）。更に、このラベルを付した画素に着目して、同様に隣接する画素のうち０に設定された画素の有無が判別されラベル付けが実行される。このように同一のラベルを付すべき画素が無くなるまで（ステップＳ３５；Ｙｅｓ）繰り返し実行されることにより、同一のラベルが付された画素の集まり、即ち、一つの第１の二値化領域が形成される。
以後、同様の処理を繰り返すことにより、二値化処理により抽出された０に設定された画素の集まり毎に異なるラベルが付された複数の第１の二値化領域が抽出される。

ステップＳ１０後、ラベル毎に各ラベルに対応する第１の二値化領域の特徴量が算出され、当該特徴量から重心座標が算出される（ステップＳ１１）。特徴量としては、各ラベルに対応する第１の二値化領域の大きさ及び形状であり、形状としては、周囲長、円形度又はアスペクト比等である。

まず、第１の二値化領域の大きさとして面積が算出される。面積は、同一のラベルが付された画素の総数が計数されることにより求められる。次に、周囲長が算出される。周囲長は、第１の二値化領域の最も外側に位置する画素の数を計数することにより求められる。そして、算出された面積及び周囲長に基づいて円形の度合いを示す円形度（０〜１に正規化された値）が算出される。円形度ｅは、面積をＳ、周囲長をＬとして、
ｅ＝４πＳ／Ｌ^２
の算出式により求められる。円形度は、例えば、円では１．０の値、正方形では０．７９、長方形では０．７９未満の値を示し、円に近いほど大きくなり１に近づく。
また、アスペクト比が算出される。更に、算出された特徴量に基づいて、第１の二値化領域の重心座標が算出される。

ラベル毎の特徴量と重心座標が算出された後（ステップＳ１１後）、ラベル毎の特徴量と、部材（金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７）毎に予め設定されている特徴量とが相互に比較される（ステップＳ１２）。

例えば、特徴量と評価部材との関係は、表１のようになる。ステップＳ１２において、全ての特徴量が比較されるものに限らず、評価部材によって特定の特徴量が比較されるようにすることが好ましい。

部材（金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７）毎に予め設定されている特徴量の所定範囲内の特徴量を有するラベルの有無が判別され、全ての各部材（金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７）に対応するラベルが有るか否かが判別される（ステップＳ１３）。なお、ステップＳ１３における判別にあたっては、ラベルの特徴量は二値化した際の閾値の影響を受けるため、適切な所定範囲を設定することが好ましい。

いずれかの部材に対応するラベルが無い場合、即ち、全ての各部材に対応するラベルが判別されなかった場合（ステップＳ１３；Ｎｏ）、ＱＣ画像データが適切でない等を示すエラーレポートが生成され、当該エラーレポートがＱＣ評価結果ファイル３９４に記憶されると共に、第２通信部３５を介してサービスセンタ５０に送信され（ステップＳ５）、本処理は終了される。

全ての各部材に対応するラベルが判別された場合、即ち、ラベル毎に対応する部材が特定された場合（ステップＳ１３；Ｙｅｓ）、ステップＳ１１において算出された重心座標が各評価部材の中心座標として設定される。このことにより、各部材（金属円板１ａ〜１ｃ、エッジパターン２、ステップウエッジ３ａ〜３ｄ、ジッタ４ａ、４ｂ、設置判別部材６、識別部材７）の位置が特定されたこととなる。そして、各部材の中心座標に基づいて各評価部材の評価を行なう際の基準となる解析中心座標が算出される（ステップＳ１４）。

従って、ステップＳ３〜１４は、取得された画像データに基づいて設置判別部材６の位置を特定する位置特定工程に相当する。

図８（ａ）にエッジパターン２の中心座標と解析中心座標と解析領域の位置関係の例、図８（ｂ）にステップウエッジ３ａ〜３ｄの中心座標と解析中心座標と解析領域の位置関係の例、図８（ｃ）にジッタ４ａ、４ｂの中心座標と解析中心座標と解析領域の位置関係の例を示す。なお、中心座標を黒丸、解析中心座標を黒星、解析領域を破線で囲む領域として示す。

図８（ａ）〜（ｃ）に示すように、各評価部材の中心座標と解析中心座標とは異なる座標となっており、解析中心座標を基準として予め定められた範囲が解析領域となっている。中心座標と解析中心座標との位置関係は、評価部材の大きさや配置方法（傾き）が予め定められているため、中心座標が決定されることにより、当該中心座標に基づいて解析中心座標を算出することができる。

ステップＳ１４において設置判別部材６の位置が特定された後、ＱＣファントムの撮影時におけるＱＣファントムの基準位置からの傾き量を算出する傾き量算出処理が実行される（ステップＳ１５）。

ステップＳ１５において実行される傾き量算出処理は、まず、ステップＳ１４において特定された重心位置と、二値化されたＱＣ画像データとに基づいて、設置判別部材６の予め定められた３点の座標が検出される。

図９に、設置判別部材６に予め定められた３点（Ａ、Ｂ、Ｃ）の例を示す。
図９に示すように、設置判別部材６は、読取装置２０の主走査方向Ｘに延在する辺（短辺）と主走査方向Ｘに延在する辺よりも長い副走査方向Ｙに延在する辺（長辺）とにより直角が形成されるように、両辺の端部が接続されたＬ字状の部材である。本実施の形態における予め定められた３点は、両辺の端部が接続された箇所に設定されたＢ点、Ｂ点から副走査方向Ｙに延在する長辺の端部に設定されたＡ点、Ｂ点から主走査方向Ｘに延在する短辺の端部に設定されたＣ点である。また、Ｂ点からＡ点までの長さをｌとする。

次に、検出された３点は、Ｂ点の座標（ｘ１、ｙ１）が基準座標（０，０）の位置になるように平行移動される。そして、平行移動された後のＡ点の座標（ｘ２、ｙ２）が検出される。ＱＣファントムの傾き量をθとした場合、
（ｘ２、ｙ２）＝（ｌｃｏｓθ、ｌｓｉｎθ）
の関係が成立する。この関係式に基づいて、θが算出される。

ステップＳ１５後、ＱＣファントムの撮影時においてＱＣファントムが放射線源１１と対抗する規定面（表面）と対向する面（裏面）に設置されて撮影されたことを示す非規定面情報を生成する裏向き判定処理が実行される（ステップＳ１６）。

ステップＳ１６において実行される裏向き判定処理は、まず、ステップＳ１４において特定された重心位置と、ＱＣ画像データとに基づいて、設置判別部材６の予め定められた３点の座標が検出される。

次に、検出された３点は、Ｂ点の座標（ｘ１、ｙ１）が基準座標（０，０）の位置になるように平行移動される。Ｂ点が基準座標（０，０）に固定された状態で、Ａ点の座標（ｘ２、ｙ２）が（０，ｌ）となるようにＡ点、Ｃ点が回転移動される。回転移動された後のＣ点の座標が検出され、当該Ｃ点のＸ座標が０よりも小さい場合、ＱＣファントムは、撮影時に非規定面（裏面）に設定されたと判別する非規定面情報が生成される。一方、Ｃ点のＸ座標が０以上である場合、ＱＣファントムは、撮影時に表面に設定されたと判別される。

従って、ステップＳ１５及びＳ１６は、設置判別部材６の位置に基づいて、ＱＣファントムの設置情報（傾き量、非規定面情報）を算出する設置情報算出工程に相当する。

次に、低コントラスト分解能の評価部材であるバーガーファントム５ａ〜５ｊの位置認識についての処理が実行される。
バーガーファントム５ａ〜５ｊは、その材質や画質評価項目の性質上、周辺領域とのコントラストが低く、金属性の第１の評価部材と同時に位置認識を行なうことは困難である。そこで、第１の評価部材とは分けてバーガーファントム５ａ〜５ｊの位置認識を行なうものとする。

まず、ステップＳ１４において中心座標が設定された部材うち、いずれかの部材の中心座標から予め定められた位置且つ予め定められた領域（小領域）の画像データが、ステップＳ３において抽出したＱＣ画像データから抽出される（ステップＳ１７）。

ステップＳ１７において、例えば、金属円板１ａとバーガーファントム５ａとの位置関係は、ＱＣファントム上において大きく位置がずれることが少ない。そこで、金属円板１ａからバーガーファントム５ａが存在すべき位置を中心とし、バーガーファントム５ａの平面的な面積を含む予め設定された領域の画像データがステップＳ２において抽出したＱＣ画像データから抽出される。

このように、小領域の画像データが抽出されることにより、この小領域にはバーガーファントムに対応する画像と当該バーガーファントム周辺の画像データのみが抽出される、即ち、他の部材に対応する画像データが含まれない画像データが抽出される。従って、二値化する際、バーガーファントムを抽出可能な閾値の設定を行なうことができる。

ステップＳ１７において抽出された小領域の画像データからヒストグラムが作成される（ステップＳ１８）。作成されたヒストグラムに基づいて閾値が決定され（ステップＳ１９）、決定された閾値を基準として小領域の画像データが二値化され、第２の二値化画像が生成される（ステップＳ２０）。

そして、第２の二値化画像に存在する０に設定された画素が連結した領域を第２の二値化領域として抽出してラベリングを行なうラベリング処理が実行される（ステップＳ２１）。

ステップＳ２１後、ラベルに対応する第２の二値化領域の特徴量が算出され、当該特徴量から重心座標が算出される（ステップＳ２２）。

ラベルの特徴量と重心座標が算出された後（ステップＳ２２後）、当該ラベルの特徴量と、予め設定されているバーガーファントムの特徴量とが相互に比較される（ステップＳ２３）。

ラベルの特徴量は、予め設定されている第２の評価部材としてのバーガーファントムの特徴量の所定範囲内か否かが判別される（ステップＳ２４）。なお、ステップＳ２４における判別にあたっては、ラベルの特徴量は二値化した際の閾値の影響を受けるため、適切な所定範囲を設定することが好ましい。

ラベルの特徴量は、予め設定されているバーガーファントムの特徴量の所定範囲内でない場合（ステップＳ２４；Ｎｏ）、ＱＣ画像データが適切でない等を示すエラーレポートが生成され、当該エラーレポートがＱＣ評価結果ファイル３９４に記憶されると共に、第２通信部３５を介してサービスセンタ５０に送信され（ステップＳ５）、本処理は終了される。

ラベルの特徴量は、バーガーファントムの特徴量の所定範囲内である場合、即ち、ラベルが付された第２の二値化領域がバーガーファントムが存在する領域として特定された場合（ステップＳ２４；Ｙｅｓ）、ステップＳ２２において算出された重心座標がバーガーファントムの中心座標として設定される。このことにより、バーガーファントムの位置が特定されたこととなる。そして、バーガーファントムの中心座標に基づいてバーガーファントムの評価を行なう際の基準となる解析中心座標が算出される（ステップＳ２５）。

そして、ステップＳ１４及びステップＳ２５において算出された各評価部材の解析中心座標に基づいて、各評価に用いる解析領域が抽出される（ステップＳ２６）。抽出された各解析領域に対して各画質評価項目の演算（ＱＣ演算）が実行される（ステップＳ２７）。
ステップＳ２７において実行される各画質評価項目としては、表示寸法精度評価、鮮鋭性評価、線形性評価、読取レーザ走査性評価、低コントラスト分解能評価である。

表示寸法精度評価は、金属円板１ａ〜１ｃを用いて、読取装置２０の主走査方向及び副走査方向の表示寸法精度、即ち、予め定められている金属円板１ａ〜１ｃの相互間の距離との誤差（評価結果値）が算出される。
例えば、算出された金属円板１ａ〜１ｃの夫々の中心座標に基づいて、まず、ＱＣ画像データ上の金属円板１ａ〜１ｂ間の距離が演算され、得られた値と実際にＱＣファントムを測定して得られた金属円板１ａ〜１ｂ間の距離（予め記憶されている）の実測値とに基づき、予め定められた演算式により評価結果値（％）が算出される。また、ＱＣ画像データ上の金属円板１ｂ〜１ｃ間の距離が演算され、得られた値と実際にＱＣファントムＰを測定して得られた金属円板１ｂ〜１ｃ間の距離の実測値とに基づき、予め定められた演算式により評価結果値が算出される。なお、ＱＣ画像データ上の２点間の距離は、画像データの２点間の画素数×読取画素サイズで表される。

鮮鋭性評価は、エッジパターン２を用いて、読取装置の主走査方向及び副走査方向の鮮鋭度の指標となるＭＴＦが算出される。
このＭＦＴは、エッジパターン２の解析中心位置を基に抽出された主走査方向のエッジ部の解析領域の画像と副走査方向のエッジ部の解析領域の画像が微分されて得られるラインスプレッドファンクションがフーリエ変換されることにより算出される。

線形性評価は、ステップウエッジ３ａ〜３ｄを用いて、読取装置の線形性が算出される。例えば、ＱＣ画像データ上のステップウエッジ３ａ〜３ｄのそれぞれに相当する領域の画像データの各画素の信号値の平均値が算出され、当該平均値から換算される透過線量と、ステップウエッジ３ａ〜３ｄの理論透過線量との関係が評価されることにより、線形性が算出される。

読取レーザ走査性評価は、ジッタ４ａ、４ｂを用いて、読取装置の主走査方向及び副走査方向の画素ズレ量が評価される。
この画素ズレ量は、例えば、ＱＣ画像データ上のジッタ４ａのエッジ部分から最小二乗誤差直線が算出され、ジッタ４ｂの最小二乗誤差直線の理論値との画素ズレ量が算出され、また、ＱＣ画像データ上のジッタ４ｂのエッジ部分から最小二乗誤差直線が算出され、ジッタ４ｂの最小二乗誤差直線の理論値との画素ズレ量が算出される。

低コントラスト分解能評価は、バーガーファントム５ａ〜５ｊを用いて、読取装置の低コントラスト分解能が評価される。
例えば、バーガーファントム内の平均信号値と周辺領域の平均信号値とのコントラスト及びノイズ量に基づいて、コントラスト対のノイズ比が算出され、算出されたコントラスト対ノイズ比に基づいて低コントラスト分解能が評価される。

ステップＳ２７の各画質評価項目のＱＣ演算終了後、画質評価項目毎に算出された評価結果値と予め画質評価項目毎に定められた評価基準値との比較が行なわれ、予め画質評価項目毎に設定されている許容範囲と照合し、許容範囲内であれば合格、許容範囲外であれば不合格の各画質評価項目に対する合否の判定が行なわれる（ステップＳ２８）。

本処理が実行されたＱＣ画像データに対するステップＳ２８において判定された結果（評価結果情報）に、ステップＳ２において取得された識別情報と、ステップＳ１５、１６において生成されたＱＣ設置情報（傾き量、非規定面情報）とが対応付けられた評価レポート（画像品質評価情報）が生成され、当該評価レポートがＱＣ評価結果ファイル３９４に記憶されると共に、第２通信部３５を介してサービスセンタ５０に送信され（ステップＳ２９）、本処理は終了される。サービスセンタ５０に送信された評価レポートは、レポートファイル５７１に保存され記憶される。

従って、ステップＳ２７〜２９は、取得された画像データに基づいて評価部材毎に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する評価演算工程に相当する。また、ステップＳ２９は、評価結果情報、識別情報、設置情報が対応付けられた画像品質評価情報（評価レポート）が記憶される記憶工程に相当する。

図１０を参照して、サービスセンタ５０において実行されるＱＣファントムの取扱判定処理について説明する。図１０に、取扱判定処理のフローチャートを示す。当該取扱判定処理は、操作部５２の所定の操作により取扱判定処理の実行が指示された際に、ＣＰＵ５１と記憶部５６等との協働によるソフトウエア処理により実現される処理である。

まず、操作部５２を介して指定された施設内におけるコンソール３０から送信された全ての評価レポートがレポートファイル５７１から読み出される（ステップＳ４１）。
従って、ステップＳ４１は、レポートファイル５７１に記憶された複数の評価レポートを読み出す読出工程に相当する。

読み出された評価レポートは、評価レポートに対応付けられたＱＣ識別情報毎に分別される。そして、任意のＱＣ識別情報に対応するＱＣファントムが撮影された画像データに基づいて生成された複数の評価レポートの誤設置頻度解析処理が実行される（ステップＳ４２）。

ステップＳ４２において実行される誤設置頻度解析処理は、所定のＱＣファントムに対する複数の評価レポートに含まれる設置情報に基づいて、ＱＣファントムの撮影時における誤設置の頻度（誤設置頻度）が算出される処理である。

このステップＳ４２における誤設置頻度は、所定のＱＣファントムに対する複数の評価レポートにおいて、当該ＱＣファントムの評価レポートの総数に対する傾き量が所定値以上である評価レポートの数の割合（以下、ＱＣ傾き頻度）、又は、当該ＱＣファントムの評価レポートの総数に対する非規定面情報を含む評価レポートの数の割合（以下、ＱＣ非規定面頻度）である。

ステップＳ４２の誤設置頻度解析処理により算出された誤設置頻度が所定値以上であるか否かが判別される（ステップＳ４３）。ステップＳ４３における判別は、ＱＣ傾き頻度が予め設定された頻度以上であるか否か、又は、ＱＣ非規定面頻度が予め設定された頻度以上であるか否かが判別される。

なお、ステップＳ４２、４３においては、ＱＣ傾き頻度又はＱＣ非規定面頻度を算出して判別しているが、所定のＱＣファントムの評価レポートにおいて、傾き量が所定値以上である評価レポートの数、又は、非規定面情報を含む評価レポートの数を算出し、傾き量が所定値以上である評価レポートの数が予め設定された数以上、又は、前記非規定面情報を含む評価レポートの数が予め設定された数以上か否かを判別してもよい。

ステップＳ４２の誤設置頻度解析処理により算出された誤設置頻度が所定値以上でない場合、即ち、ＱＣ傾き量頻度が予め設定された頻度未満、かつ、ＱＣ非規定面頻度が予め設定された頻度未満である場合（ステップＳ４３；Ｎｏ）、ＱＣ識別情報毎に分別された全ての各ＱＣファントムの複数の評価レポートに対する誤設置頻度解析処理が完了しているか否か判別される（ステップＳ４４）。

誤設置頻度解析処理が完了していないＱＣファントムの評価レポートがある場合（ステップＳ４４；Ｎｏ）、ステップＳ４２に戻り、誤設置頻度解析処理が完了していないＱＣファントムの評価レポートがない場合（ステップＳ４４；Ｙｅｓ）、本処理は終了される。

ステップＳ４２の誤設置頻度解析処理により算出された誤設置頻度が所定値以上である場合、即ち、ＱＣ傾き量頻度が予め設定された頻度以上、又は、ＱＣ非規定面頻度が予め設定された頻度以上である場合（ステップＳ４３；Ｙｅｓ）、ステップＳ４１において読み出された指定された施設内の全ての評価レポートの誤差設置頻度解析が実行される（ステップＳ４５）。

このステップＳ４５における誤設置頻度は、施設内の全ての評価レポートにおいて、当該施設内の全ての評価レポートの総数に対する傾き量が所定値以上である評価レポートの数の割合（施設傾き頻度）、又は、当該施設内の全ての評価レポートの総数に対する非規定面情報を含む評価レポートの数の割合（施設非規定面頻度）である。

ステップＳ４５の誤設置頻度解析処理により算出された誤設置頻度が所定値以上であるか否かが判別される（ステップＳ４６）。ステップＳ４６における判別は、施設傾き頻度が予め設定された頻度以上であるか否か、又は、施設非規定面頻度が予め設定された頻度以上であるか否かが判別される。

なお、ステップＳ４５、４６においては、施設傾き頻度又は施設非規定面頻度を算出して判別しているが、施設内の全評価レポートにおいて、傾き量が所定値以上である評価レポートの数、又は、非規定面情報を有する評価レポートの数を算出し、傾き量が所定値以上である評価レポートの数が予め設定された数以上、又は、前記非規定面情報を含む評価レポートの数が予め設定された数以上か否かを判別してもよい。

ステップＳ４５の誤設置頻度解析処理により算出された誤設置頻度が所定値以上でない場合、即ち、施設傾き量頻度が予め設定された頻度未満、かつ、施設非規定面頻度が予め設定された頻度未満である場合（ステップＳ４６；Ｎｏ）、ＱＣファントム位置決めタブが欠損したり位置決めマークが消失したりする等の破損が生じているものと判断して、ＱＣファントムの破損確認を指示する注意喚起情報を示すメッセージが表示部５３に表示されることにより報知され（ステップＳ４７）、本処理は終了される。

ステップＳ４７において表示されるメッセージに従い、サービスセンタ５０から当該施設のサービスマンにＱＣファントム交換の指示が出され、サービスマンは、ＱＣファントムを交換する。

ステップＳ４５の誤設置頻度解析処理により算出された誤設置頻度が所定値以上である場合、即ち、施設傾き量頻度が予め設定された頻度以上、又は、施設非規定面頻度が予め設定された頻度以上である場合（ステップＳ４６；Ｙｅｓ）、当該設置施設におけるＱＣファントムの設置・取り扱い手順に不備があると判断して、ＱＣファントムの撮影時の取り扱い手順の見直しを指示する注意喚起情報を示すメッセージが表示部５３に表示されることにより報知され（ステップＳ４８）、本処理は終了される。

ステップＳ４８においてメッセージの報知を受けたサービスマンは、当該施設におもむいて、ＱＣファントムの設置状況を検分し、必要に応じて取り扱いの改善指導、あるいは位置決めタブや位置決めマークの付与されたＱＣファントムへ交換を行なう。

以上のように、本実施の形態によれば、放射線画像の被写体であるＱＣファントムを識別するＱＣ識別情報と当該ＱＣファントムの設置情報とを当該放射線画像の評価結果情報に対応付けて評価レポートとしてサービスセンタ５０のレポートファイル５７１に記憶することができるため、各ＱＣファントムの撮影時における設置傾向を把握することが可能となり、ＱＣファントムの取り扱いに起因する画質評価の不良において、ＱＣファントムの取り扱いに対する対策を講じさせ得ることができる。

更に、評価レポートに含まれるＱＣ設置情報に基づいて判定された注意喚起情報としてのメッセージを表示して報知することができるため、当該メッセージが報知されたユーザは、ＱＣファントムの取り扱いに対する適切な対策を講じ得ることができる。

メッセージとしては、ＱＣ傾き量頻度が所定頻度以上、又は、ＱＣ非規定面頻度が所定頻度以上である場合、且つ、施設傾き量頻度が所定頻度以上、又は、施設非規定面頻度が所定頻度以上である場合、ＱＣファントムの撮影時の取り扱い手順の見直しを指示する注意喚起情報を報知することができ、ＱＣファントムの撮影時の取り扱いによる画像品質の低下の是正を促すことができる。また、ＱＣ傾き量頻度が所定頻度以上、又は、ＱＣ非規定面頻度が所定頻度以上である場合、且つ、施設傾き量頻度が所定頻度未満、及び、施設非規定面頻度が所定頻度未満の場合、ＱＣファントムの破損確認を指示する注意喚起情報を報知することができ、ＱＣファントムの破損による画像品質の低下の是正を促すことができる。

また、ＱＣファントムに自己の識別情報を示す識別部材７が内在されているため、放射線画像に基づいてＱＣ識別情報を取得することにより、ユーザがＱＣファントムのＱＣ識別情報の入力する手間を省き作業性を向上できると共に、誤って他のＱＣファントムのＱＣ識別情報の入力が行なわれることを防止できるため、評価結果情報とＱＣ識別情報との対応付けの信頼性を向上できる。従って、ＱＣファントムの取り扱いに対する正確な対策を講じさせ得ることができる。

なお、本実施の形態においては、ＱＣファントム毎、施設毎に誤設置頻度を算出してメッセージを報知するものであるが、ＱＣファントムの撮影時において、当該ＱＣファントムを撮影した技師の識別情報を当該放射線画像に対応付けて記憶しておくことにより、技師毎に誤設置頻度を算出して各技師に対してメッセージを報知することにより、ＱＣファントムの取り扱いの改善を図ることができる。

また、本発明は、上記実施の形態の内容に限定されるものではなく、本発明の主旨を逸脱しない範囲で適宜変更可能である。

放射線画像読取システムの全体構成図である。 コンソールの機能構成例を示す図である。 ＱＣファントムの一例を示す図である。 サービスセンタの機能構成例を示す図である。 品質評価処理のフローチャートである。 品質評価処理のフローチャートである（図５の続き）。 ラベリング処理のフローチャートである。 図８（ａ）は、エッジパターンの中心座標と解析中心座標と解析領域の位置関係の例、図８（ｂ）は、ステップウエッジの中心座標と解析中心座標と解析領域の位置関係の例、図８（ｃ）は、ジッタの中心座標と解析中心座標と解析領域の位置関係の例を示す図である。 設置判別部材６に予め定められた３点（Ａ、Ｂ、Ｃ）の例を示す図である。 取扱判定処理のフローチャートである。

符号の説明

１ａ、１ｂ 金属円板
２ エッジパターン
３ａ〜３ｄ ステップウエッジ
４ａ、４ｂ ジッタ
５ａ〜５ｊ バーガーファントム
６ 設置判別部材
７ 識別部材
１０ 撮影装置
１１ 放射線源
１２ カセッテ
２０ 読取装置
３０、３０ａ コンソール
３１ ＣＰＵ
３２ 操作部
３３ 表示部
３４ 第１通信部
３５ 第２通信部
３６ バーコードリーダ
３７ ＲＡＭ
３８ 記憶部
３９ ＨＤＤ
３９１ ＱＣ撮影管理情報ファイル
３９２ 一般撮影管理情報ファイル
３９３ 画像ファイル記憶領域
３０１ ＱＣツール
３９４ ＱＣ評価結果ファイル
４０ イメージャ
５０ サービスセンタ
５１ ＣＰＵ
５２ 操作部
５３ 表示部
５４ 通信部
５５ ＲＡＭ
５６ 記憶部
５７ ＨＤＤ
５７１ レポートファイル
Ｍ 被写体
Ｎ１、Ｎ２ 通信ネットワーク
Ｎ３、Ｎ４ バス
ＮＳ 放射線画像読取システム
Ｐ ＱＣファントム
Ｘ 主走査方向
Ｙ 副走査方向

Claims (4)

1. 左右非対称及び上下非対称の形状を有する設置判別部材と、画質評価項目に対応した大きさ又は形状を有する評価用部材とが内在されたＱＣファントムが撮影された放射線画像を取得する画像取得工程と、
   前記ＱＣファントムを個別的に識別する識別情報を取得する識別情報取得工程と、
   前記画像取得工程により取得された放射線画像に基づいて前記設置判別部材の位置を特定する位置特定工程と、
   前記位置特定工程により特定された前記設置判別部材の位置に基づいて、前記ＱＣファントムの設置情報を算出する設置情報算出工程と、
   前記画像取得工程により取得された放射線画像に基づいて前記評価用部材に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する評価演算工程と、
   前記評価演算工程により生成された評価結果情報と、前記識別情報取得工程により取得された前記識別情報と、前記設置情報算出工程により算出された前記設置情報とが対応付けられた画像品質評価情報を記憶する記憶工程と、
   前記記憶工程により記憶された複数の前記画像品質評価情報を読み出す読出工程と、
   読み出された前記複数の画像品質評価情報に含まれる前記設置情報に基づいて、前記ＱＣファントムに対する取り扱いの注意喚起情報を報知するか否かを判定する報知判定工程と、
   前記報知判定工程により前記注意喚起情報を報知すると判定された場合、当該注意喚起情報を報知する報知工程と、
   を含み、
   前記設置情報は、前記ＱＣファントムの撮影時における前記ＱＣファントムの基準位置からの傾き量、又は、前記ＱＣファントムの撮影時において前記ＱＣファントムが照射線源に対向する規定面と対向する面に設置されて撮影されたことを示す非規定面情報を含み、
   前記報知判定工程は、
   所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、
   且つ、
   前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、
   前記ＱＣファントムの撮影時の取り扱い手順の見直しを指示する前記注意喚起情報を報知すると判定すること、
   を特徴とする画像品質管理支援方法。
2. 前記報知判定工程は、
   所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、
   且つ、
   前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数未満、及び、前記非規定面情報を含む画像品質評価情報の数が予め設定された数未満の場合、
   前記ＱＣファントムの破損確認を指示する前記注意喚起情報を報知すると判定すること、
   を特徴とする請求項１に記載の画像品質管理支援方法。
3. 左右非対称及び上下非対称の形状を有する設置判別部材と、画質評価項目に対応した大きさ又は形状を有する評価用部材とが内在されたＱＣファントムが撮影された放射線画像を取得する画像取得工程と、
   前記ＱＣファントムを個別的に識別する識別情報を取得する識別情報取得工程と、
   前記画像取得工程により取得された放射線画像に基づいて前記設置判別部材の位置を特定する位置特定工程と、
   前記位置特定工程により特定された前記設置判別部材の位置に基づいて、前記ＱＣファントムの設置情報を算出する設置情報算出工程と、
   前記画像取得工程により取得された放射線画像に基づいて前記評価用部材に対応する評価演算を行い、当該評価演算に基づく評価結果情報を生成する評価演算工程と、
   前記評価演算工程により生成された評価結果情報と、前記識別情報取得工程により取得された前記識別情報と、前記設置情報算出工程により算出された前記設置情報とが対応付けられた画像品質評価情報を記憶する記憶工程と、
   前記記憶工程により記憶された複数の前記画像品質評価情報を読み出す読出工程と、
   読み出された前記複数の画像品質評価情報に含まれる前記設置情報に基づいて、前記ＱＣファントムに対する取り扱いの注意喚起情報を報知するか否かを判定する報知判定工程と、
   前記報知判定工程により前記注意喚起情報を報知すると判定された場合、当該注意喚起情報を報知する報知工程と、
   を含み、
   前記設置情報は、前記ＱＣファントムの撮影時における前記ＱＣファントムの基準位置からの傾き量、又は、前記ＱＣファントムの撮影時において前記ＱＣファントムが照射線源に対向する規定面と対向する面に設置されて撮影されたことを示す非規定面情報を含み、
   前記報知判定工程は、
   所定の前記ＱＣファントムに対する複数の画像品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数以上、又は、前記非規定面情報を含む画像品質評価情報の数が予め設定された数以上の場合、
   且つ、
   前記記憶工程により記憶された複数の画総品質評価情報において、前記傾き量が所定値以上である画像品質評価情報の数が予め設定された数未満、及び、前記非規定面情報を含む画像品質評価情報の数が予め設定された数未満の場合、
   前記ＱＣファントムの破損確認を指示する前記注意喚起情報を報知すると判定すること、
   を特徴とする画像品質管理支援方法。
4. 前記ＱＣファントムは、自己の識別情報を示す識別部材を内在しており、
   前記識別情報取得工程は、
   画像取得工程により取得された放射線画像に基づいて前記識別部材が示す識別情報を取得すること、
   を特徴とする請求項１から３のいずれか一項に記載の画像品質管理支援方法。

Images