(57)【特許請求の範囲】
【請求項１】
所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体を、画像読取装置でサンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取方法において、
前記記録手段によって一様の入力がなされた前記画像記録媒体に対して前記画像読取装置で読取走査を行って前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得し、これを補正データ値Ｓc0（ｘ，ｙ）とし、
前記画像読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とし、
前記記録手段により画像が記録された前記画像記録媒体に対して前記画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ0（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出すること
を特徴とする画像読取方法。
【請求項２】
請求項１に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆx，ｆy）の関数である第２の補正データ値Ｓc2（ｆx，ｆy）とし、
Ｆ（ｆx，ｆy）＝１  （０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆx，ｆy）＝０  （ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆx，ｆy）を用いて、第３の補正データ値Ｓc3（ｆx，ｆy）
Ｓc3（ｆx，ｆy）＝Ｆ（ｆx，ｆy）・Ｓc2（ｆx，ｆy）、
とし、この第３の補正データ値Ｓc3（ｆx，ｆy）を逆フーリエ変換して、前記最終補正データＳc（ｘ，ｙ）とすること、
を特徴とする画像読取方法。
【請求項３】
請求項１に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする画像読取方法。
【請求項４】
請求項３に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数を用いること
を特徴とする画像読取方法。
【請求項５】
請求項３に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いること
を特徴とする画像読取方法。
【請求項６】
請求項１〜５のいずれかに記載の画像読取方法において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする画像読取方法。
【請求項７】
所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体を、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取装置において、
前記画像記録媒体のＩＤを検出可能なＩＤ検出部と、
前記画像記録媒体に対して読取走査を行って前記各読取座標に対応する画像データ値を取得し、これを画像データとして記録する画像データ記録手段と、
一様の入力がなされた前記画像記録媒体に対して、前記画像データ記録手段で読取走査を行い、前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得して、これを補正データ値Ｓc0（ｘ，ｙ）として記録し、
前記読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc0（ｘ，ｙ）のうち前記空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を除去し、これを最終補正データ値Ｓc（ｘ，ｙ）として、この画像記録媒体のＩＤと対応して記録する補正データ算出手段と、
前記補正前画像データ値Ｓ0（ｘ，ｙ）を前記最終補正データ値Ｓc（ｘ，ｙ）により補正して補正後画像データ値Ｓ（ｘ，ｙ）を算出する画像データ補正手段と、を備えること
を特徴とする画像読取装置。
【請求項８】
請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆx，ｆy）の関数である第２の補正データ値Ｓc2（ｆx，ｆy）を算出し、
Ｆ（ｆx，ｆy）＝１  （０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆx，ｆy）＝０  （ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆx，ｆy）を用いて、第３の補正データ値Ｓc3（ｆx，ｆy）
Ｓc3（ｆx，ｆy）＝Ｆ（ｆx，ｆy）・Ｓc2（ｆx，ｆy）、
を算出し、この第３の補正データ値Ｓc3（ｆx，ｆy）を逆フーリエ変換して、空間（ｘ，ｙ）の関数である最終補正データＳc（ｘ，ｙ）を算出すること
を特徴とする画像読取装置。
【請求項９】
請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする画像読取装置。
【請求項１０】
請求項９に記載の画像読取装置において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数が用いられていること
を特徴とする画像読取装置。
【請求項１１】
請求項９に記載の画像読取装置において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算が用いられていること
を特徴とする画像読取装置。
【請求項１２】
請求項７〜１１のいずれかに記載の画像読取装置において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする画像読取装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像情報方法および画像読取装置に係り、特に画像を記録する画像記録媒体の感度ムラを効果的に補正可能な画像読取方法および画像読取装置に関する。
【０００２】
【従来の技術】
画像が記録された記録媒体にレーザなどの光ビームを照射し、その記録媒体からの反射光、透過光あるいは発光光を検出することにより、画像を画像データに変換する画像読取方法や画像読取装置が知られている。このような画像読取装置として、製版用スキャナー、コンピュータやファクシミリの入力装置等のほか、輝尽性蛍光体プレートを使用した放射線画像情報記録再生システムにおいて用いられる放射線画像読取装置がある。
【０００３】
例えば、上記放射線画像読取装置において画像記録媒体として用いられる輝尽性蛍光体プレートで記録される画像には、集光系や光学系による画像ムラのほか、輝尽性蛍光体の２次元的な感度ムラが生じる。そこで、この感度ムラを補正するため、各輝尽性蛍光体プレートの感度ムラをあらかじめ求めておき、この輝尽性蛍光体プレートから読み取った画像データから補正データを加算するなどして、画像データを補正することが行われている。
【０００４】
ここで、輝尽性蛍光体プレートがカセッテ内に収納されたカセッテ方式の画像記録媒体では、カセッテを画像読取装置の読取部にセットする際のセット位置が、０．数ｍｍ〜数ｍｍ程度の範囲内で、カセッテをセットする毎に変動することがある。ここで、画像を読み取る際の精度を高くする場合には、この微小なセット位置のズレが、画像データの補正に影響を与える。
【０００５】
具体的には、被写体がない状態で画像を記録した輝尽性蛍光体プレートに形成された画像を、画像読取装置によりサンプリングピッチをＰで画像データに変換すると、この画像データは、図６（ａ）、図６（ｂ）中に太実線で示すように、上記感度ムラを包含したものとなる。被写体がない状態では画像データは一定値（図中、細実線）になるべきである。したがって、輝尽性蛍光体プレートに記録された被写体の画像を画像読み取り装置で読み込んで得られた画像データを補正するためには、画像データを太実線で除すれば良い。
【０００６】
しかし、輝尽性蛍光体プレートの画像読取装置媒体へのセット位置が、例えば１サンプリングピッチＰ分だけずれた状態（図中、太点線）で読み取られた画像データに対して上記方法により補正を施すと、補正後の画像データは図中に細点線で示すものとなる。このように、読取時のサンプリングピッチが小さい場合には、補正後の画像データの高周波のノイズが最大で補正前の画像データのノイズの２倍まで増幅されてしまう恐れがある。
【０００７】
そこで、画像読取時のサンプリングピッチに応じて、除去するノイズの周波数を考慮しつつ、低周波から高周波にわたるノイズを効果的に除去し、画像データの補正を良好に行うことが考えられる。
【０００８】
【発明が解決しようとする課題】
本発明の課題は、画像記録媒体に記録された画像を読み取る画像読取装置において、画像記録媒体の画像読取装置への装填位置の誤差があっても、画像記録媒体の感度ムラの大部分を補正可能な、画像読取方法、および画像読取装置を提供することである。
【０００９】
【課題を解決するための手段】
請求項１に記載の発明は、例えば図１〜図５に示すように、所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体（輝尽性蛍光体プレート）を、画像読取装置でサンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取方法において、
前記記録手段によって一様の入力がなされた前記画像記録媒体に対して前記画像読取装置で読取走査を行って前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得し、これを補正データ値Ｓc0（ｘ，ｙ）とし、
前記画像読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とし、
前記記録手段により画像が記録された前記画像記録媒体に対して前記画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ0（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出すること
を特徴とする。
【００１０】
ここで、上記画像記録媒体、記録手段とは、それぞれ、例えば放射線撮影装置において記録手段として用いられる輝尽性蛍光体プレート、Ｘ線のことである。
また、上記読取走査とは、例えば輝尽性蛍光体プレートをレーザ光などの励起光で走査して輝尽発光光を発生させ、得られた輝尽発光光を光電的に読み出して、画像信号を得ることである。
また、上記画像データ値がリニアな値である場合には、補正前画像データ値Ｓ0（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を乗算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出する。上記画像データ値がＬｏｇ変換された値である場合には、補正前画像データ値Ｓ0（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を和算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出する。
【００１１】
請求項１に記載の発明によれば、一様の入力がなされた画像記録媒体に対して画像読取装置でピッチＰで読取走査を行って補正データ値Ｓc0（ｘ，ｙ）を取得し、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この画像記録媒体に記録された被写体の画像に対して画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ0（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ0（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００１２】
請求項２に記載の発明は、請求項１に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆx，ｆy）の関数である第２の補正データ値Ｓc2（ｆx，ｆy）とし、
Ｆ（ｆx，ｆy）＝１  （０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆx，ｆy）＝０  （ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆx，ｆy）を用いて、第３の補正データ値Ｓc3（ｆx，ｆy）
Ｓc3（ｆx，ｆy）＝Ｆ（ｆx，ｆy）・Ｓc2（ｆx，ｆy）、
とし、この第３の補正データ値Ｓc3（ｆx，ｆy）を逆フーリエ変換して、前記最終補正データＳc（ｘ，ｙ）とすること、
を特徴とする。
【００１３】
請求項２に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換して第２の補正データ値Ｓc2（ｆx，ｆy）とし、このうちから（１／（４ｎ・Ｐ））を超える周波数成分を除去して第３の補正データ値Ｓc3（ｆx，ｆy）とし、これを逆フーリエ変換して、最終補正データＳc（ｘ，ｙ）とするので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、確実に行える。
【００１４】
請求項３に記載の発明は、請求項１に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする。
【００１５】
請求項３に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出するので、請求項２に記載の発明のようにフーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを行える。
【００１６】
請求項４に記載の発明は、請求項３に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数を用いること
を特徴とする。
【００１７】
請求項４に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から得られる前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、空間周波数（ｆx，ｆy）領域で高周波成分を除去するゲート関数を逆フーリエ変換した結果得られるｓｉｎｃ関数を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行える。
【００１８】
請求項５に記載の発明は、請求項３に記載の画像読取方法において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いること
を特徴とする。
ここで、移動平均法とは、読取座標（ｉ，ｊ）の近傍の画像データ値の平均を、読取座標（ｉ，ｊ）における画像データ値とする方法のことを言う。
【００１９】
請求項５に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、簡単な計算方法で確実に行える。
【００２０】
請求項６に記載の発明は、請求項１に記載の画像読取方法において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする。
【００２１】
請求項６に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、輝尽性蛍光体プレートを備える画像記録媒体に一様の入力を施す際に用いる放射線の照射線量が２．５×１０-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc0（ｘ，ｙ）を、もっぱら画像記録媒体の構造モトルに起因するものとすることができる。
【００２２】
請求項７に記載の発明は、例えば図１〜図５に示すように、所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体を、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取装置において、
前記画像記録媒体のＩＤを検出可能なＩＤ検出部と、
前記画像記録媒体に対して読取走査を行って前記各読取座標に対応する画像データ値を取得し、これを画像データとして記録する画像データ記録手段と、
一様の入力がなされた前記画像記録媒体に対して、前記画像データ記録手段で読取走査を行い、前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得して、これを補正データ値Ｓc0（ｘ，ｙ）として記録し、
前記読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc0（ｘ，ｙ）のうち前記空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を除去し、これを最終補正データ値Ｓc（ｘ，ｙ）として、この画像記録媒体のＩＤと対応して記録する補正データ算出手段と、
前記補正前画像データ値Ｓ0（ｘ，ｙ）を前記最終補正データ値Ｓc（ｘ，ｙ）により補正して補正後画像データ値Ｓ（ｘ，ｙ）を算出する画像データ補正手段と、を備えること
を特徴とする。
【００２３】
請求項７に記載の発明によれば、請求項１に記載の発明と同様の効果が得られる。
【００２４】
請求項８に記載の発明は、請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆx，ｆy）の関数である第２の補正データ値Ｓc2（ｆx，ｆy）を算出し、
Ｆ（ｆx，ｆy）＝１  （０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆx，ｆy）＝０  （ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆx，ｆy）を用いて、第３の補正データ値Ｓc3（ｆx，ｆy）
Ｓc3（ｆx，ｆy）＝Ｆ（ｆx，ｆy）・Ｓc2（ｆx，ｆy）、
を算出し、この第３の補正データ値Ｓc3（ｆx，ｆy）を逆フーリエ変換して、空間（ｘ，ｙ）の関数である最終補正データＳc（ｘ，ｙ）を算出すること
を特徴とする。
【００２５】
請求項８に記載の発明によれば、請求項２に記載の発明と同様の効果が得られる。
【００２６】
請求項９に記載の発明は、請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc0（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆxまたはｆy）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする。
【００２７】
請求項９に記載の発明によれば、請求項３に記載の発明と同様の効果が得られる。
【００２８】
請求項１０に記載の発明は、請求項９に記載の画像読取装置において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数が用いられていること
を特徴とする。
【００２９】
請求項１０に記載の発明によれば、請求項４に記載の発明と同様の効果が得られる。
【００３０】
請求項１１に記載の発明は、請求項９に記載の画像読取装置において、
前記補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算が用いられていること
を特徴とする。
【００３１】
請求項１１に記載の発明によれば、請求項５に記載の発明と同様の効果が得られる。
【００３２】
請求項１２に記載の発明は、請求項７〜１１のいずれかに記載の画像読取装置において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする。
【００３３】
請求項１２に記載の発明によれば、請求項６に記載の発明と同様の効果が得られる。
【００３４】
【発明の実施の形態】
以下、本発明に係る画像読取方法、および画像読取装置について、図１〜図５を参照して説明する。
本実施の形態の画像読取装置１は、放射線画像記録読取システムにおいて、撮影台３に装填された状態で画像が撮影記録される輝尽性蛍光体プレート（画像記録媒体）２１の画像を読み取り、画像データに変換する放射線画像読取装置である。
輝尽性蛍光体プレート２１は、図２に示すように、カセッテ２に収容されており、撮影台３にセットされた状態でＸ線発生装置４から出射されたＸ線４１のうち人体などの被写体５およびその周囲を通過するＸ線の照射エネルギーを蓄積することで、被写体５のＸ線画像を記録するものである。この輝尽性蛍光体プレート２１に、可視光線などの励起光を照射すると、蓄積されたエネルギーが解放され、輝尽発光が生じる。この現象を利用して、輝尽性蛍光体プレート２１に記録された画像の読み取りや消去を画像読取装置１で行う。
【００３５】
画像読取装置１は、図１、図３に示すように、上記カセッテ２に収容された輝尽性蛍光体プレート２１を引き出し、サンプリングピッチＰの読取座標（ｘ，ｙ）毎にレーザ光などの励起光で輝尽発光させつつ読取走査し、輝尽性蛍光体プレート２１に記録された画像を各読取座標（ｘ，ｙ）に対応する画像データ値に変換するものである。
【００３６】
画像記録装置１には、カセッテ２を装填するためのスタッカ部１１ａ〜１１ｅが備えられ、これら各スタッカ１１ａ〜１１ｅには、装填されたカセッテ２のＩＤを検出するためのＩＤ検出部１２が設けられている。
そして、表示・操作部１７のＣＲＴ表示部１７１前面のタッチパネル１７２を操作すると、その信号がＣＰＵ１６１に伝達される。そして、ＣＰＵ１６１の命令が、読取部・画像入力制御部１６２、メカ制御部１６３を介してカセッテ搬送機構１３に伝達され、このカセッテ搬送機構１３の動作により、スタッカ部１１ａ〜１１ｅに装填されたカセッテ２のうち所望のカセッテ２が、スタッカ部１１ａ〜１１ｅとカセッテ２に記録された画像の読み取りを行う読取位置との間を搬送される。
また、同タッチパネル１７２の操作により、ＣＰＵ１６１の命令が、読取部・画像入力制御部１６２を介してカセッテ搬送機構１３に設けられた消去部１４に伝達され、消去部１４から可視光線がカセッテ２の記録面全面に照射されて、カセッテ２に記録された画像を消去できるようになっている。
【００３７】
タッチパネル１７２で、所望のカセッテ２の読み取りを指定した場合には、そのカセッテ２に収容された輝尽性蛍光体プレート２１が、カセッテ搬送機構１３の動作によって上記読取位置に搬送され、画像読取部１５からレーザ光からなる励起光を照射されて、この励起光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査が行われる。この読取走査によって、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値が取得される。
【００３８】
この画像データ値は、システム制御部１６に備えられたＣＰＵ１６１の命令により、画像データとしてハードディスク（画像データ記録手段）１６４に記録される。
ハードディスク１６４には、後述する補正データ算出処理や、画像データ補正処理を行うためのプログラムや、各輝尽性蛍光体プレート２１についての最終補正データ（後述）も記録されている。これらの画像データやプログラムは、必要に応じてＣＰＵ１６１から読み出される。ＣＰＵ１６１による処理で扱われるデータやプログラムは、一時的にＲＡＭ（Read-Only Memory）１６５に記憶される。
【００３９】
また、ＣＰＵ１６１は、インターフェイス１６６を介して、ＬＡＮ９に接続され、このＬＡＮ９に接続されている、各患者の情報を管理するホストコンピュータ６や、画像データを出力するためのハードコピー装置７、患者登録ターミナル８などとの間で、各種データの送受信を行えるようになっている。
【００４０】
この画像読取装置１を利用した、輝尽性蛍光体プレート２１に記録された画像の読取方法について説明する。ここでは、画像データ値がリニアな値である場合について説明する。
まず、図２に示すように、輝尽性蛍光体プレート２１を収容するカセッテ２を撮影台３にセットし、この状態で撮影台の前に被写体５を置かずに、Ｘ線発生装置４により照射線量が２．５×１０-6［Ｃ／ｋｇ］以上、好ましくは２．５×１０-5［Ｃ／ｋｇ］の放射線を撮影台３に照射して、このＸ線４１でカセッテ２内の輝尽性蛍光体プレート２１を励起させる。
次に、このカセッテ２を撮影台３から取り出して、画像読取装置１のスタッカ部１１ａ〜１１ｅのいずれかに装填すると、カセッテ２のＩＤがＩＤ検出部１２によって読み出される。そして、画像読取部１５からレーザ光からなる励起光を照射させて、この励起光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査を行う。この読取走査によって、カセット２サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値を取得し、これが補正データ値Ｓc0（ｘ，ｙ）として、輝尽性蛍光体プレート２１が収容されていたカセッテ２のＩＤと対応してハードディスク１６４に記録される。
【００４１】
この第１の補正データ値は、表示・操作部１７のタッチパネル１７２を操作することにより、ＣＰＵ１６１によってハードディスク１６４から読み出された二次元ＦＦＴ（Fast Fourier Transformation：高速フーリエ変換）プログラムによりフーリエ変換が施され、各空間周波数に対応する画像データ値が出力される。この画像データ値がＲＡＭ１６５に第２の補正データ値Ｓc2（ｆx，ｆy）として一時的に記録される。
【００４２】
画像読取装置１の上記読取位置に輝尽性蛍光体プレート２１がセットされる時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎを、あらかじめ設定しておく。
そして、ＣＰＵ１６１により、第２の補正データ値Ｓc2（ｆx，ｆy）のうち空間周波数ｆがｆc＝１／（４ｎ・Ｐ）を超える各空間周波数ｆ（ｆxまたはｆy）に対応する前記画像データ値を０に置換し、その結果が第３の補正データ値Ｓc3（ｆx，ｆy）としてＲＡＭ１６５に一時的に記録される。
【００４３】
続いて、ＣＰＵ１６１によって、第３の補正データ値Ｓc3（ｆx，ｆy）に二次元逆フーリエ変換が施され、再び各読取座標（ｘ，ｙ）に対応する画像データ値が出力される。この画像データ値をその最大値に対する比として規準化したものが、最終補正データ値Ｓc（ｆx，ｆy）として、ハードディスク１６４に記録される。すなわち、最終補正データ値Ｓc（ｆx，ｆy）は、０＜Ｓc（ｆx，ｆy）≦１となる。そして、この輝尽性蛍光体プレート２１に記録された画像を、カセッテ搬送機構１３の消去部１４で消去し、再びカセッテ２に収容し、このカセッテ２を装填したスタッカ１１ａ〜１１ｅのいずれかから取り出す。
以上で、補正データ算出処理が終了する。
【００４４】
次いで、図２に示すように、このカセッテ２を撮影台３に装填し、撮影台３の全面に人体などの被写体５を配置した状態で、その背後よりＸ線発生手段でＸ線を照射し、カセッテ２の輝尽性蛍光体２１を励起させ、被写体５の画像を記録する。
【００４５】
次に、このカセッテ２を撮影台３から取り出して、画像読取装置１のスタッカ部１１ａ〜１１ｅのいずれかに装填すると、カセッテ２のＩＤがＩＤ検出部１２によって読み出される。そして、画像読取部１５からレーザ光からなる励起光を輝尽性蛍光体プレート２１に照射して、この励起光による輝尽発光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査を行う。この読取走査によって、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値が取得され、これが補正前画像データ値Ｓ0（ｘ，ｙ）として、輝尽性蛍光体プレート２１が収容されていたカセッテ２のＩＤと対応してＲＡＭ１６５またはハードディスク１６４に記録される。
【００４６】
そして、補正前画像データ値Ｓ0（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を乗算して、この結果が補正後画像データ値Ｓ（ｘ，ｙ）としてハードディスク１６４に記録される。
以上で、画像データ補正処理が終了する。
【００４７】
本実施の形態に記載の画像読取方法、および画像読取装置１によれば、一様の入力がなされた輝尽性蛍光体プレート２１に対して画像読取装置１でピッチＰで読取走査を行って補正データ値Ｓc0（ｘ，ｙ）を取得し、画像読取装置１にカセッテ２に収容された輝尽性蛍光体プレート２１をセットする時にズレる可能性のある最大画素数ｎに応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この輝尽性蛍光体プレート２１に記録された被写体の画像に対して画像読取装置１で読取走査を行って取得される補正前画像データ値Ｓ0（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ0（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００４８】
また、輝尽性蛍光体プレート２１を備えるカセッテ２に一様の入力を施す際に用いる放射線の照射線量が２．５×１０-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc0（ｘ，ｙ）を、もっぱら輝尽性蛍光体プレート２１の構造モトルに起因するものとすることができる。
したがって、カセッテ２の輝尽性蛍光体プレート２１の構造モトルによるノイズを効果的に除去することができる。
【００４９】
なお、本発明の画像読取方法、および画像読取装置は、上記の実施の形態に限定されることなく、本発明の趣旨を逸脱しない範囲において、種々の改良ならびに設計の変更を行ってもよい。
例えば、上記実施の形態においては、補正データ値Ｓc0（ｘ，ｙ）から最終補正データ値Ｓc（ｘ，ｙ）を算出する方法として、フーリエ変換および逆フーリエ変換を用いる場合について説明したが、ｓｉｎｃ関数や移動平均法を用いて、空間（ｘ，ｙ）領域において補正データ値Ｓc0（ｘ，ｙ）から最終補正データ値Ｓc（ｘ，ｙ）を得るようにすれば、フーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行うことができる。
また、上記画像データ値がＬｏｇ変換された値である場合には、高周波成分を除去した補正データ値Ｓc0（ｘ，ｙ）を、その最大値に対する差分で規準化して最終補正データ値Ｓc（ｘ，ｙ）とし、補正前画像データ値Ｓ0（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を和算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出すれば良い。
また、上記実施の形態においては、放射線画像読取装置について説明したが、本発明は一般のイメージスキャナなどにも適用可能である。
その他、具体的な細部構造などについても適宜に変更可能であることももちろんである。
【００５０】
【発明の効果】
請求項１または請求項７に記載の発明によれば、一様の入力がなされた画像記録媒体に対して画像読取装置でピッチＰで読取走査を行って補正データ値Ｓc0（ｘ，ｙ）を取得し、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この画像記録媒体に記録された被写体の画像に対して画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ0（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ0（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００５１】
請求項２または請求項８に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）の集合をフーリエ変換して第２の補正データ値Ｓc2（ｆx，ｆy）とし、このうちから（１／（４ｎ・Ｐ））を超える周波数成分を除去して第３の補正データ値Ｓc3（ｆx，ｆy）とし、これを逆フーリエ変換して、最終補正データＳc（ｘ，ｙ）とするので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、確実に行える。
【００５２】
請求項３または請求項９に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出するので、請求項２に記載の発明のようにフーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを行える。
【００５３】
請求項４または請求項１０に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、空間周波数（ｆx，ｆy）領域で高周波成分を除去するゲート関数を逆フーリエ変換した結果得られるｓｉｎｃ関数を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行える。
【００５４】
請求項５または請求項１１に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc0（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc0（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、簡単な計算方法で確実に行える。
【００５５】
請求項６または請求項１２に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、輝尽性蛍光体プレートを備える画像記録媒体に一様の入力を施す際に用いる放射線の照射線量が２．５×１０-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc0（ｘ，ｙ）を、もっぱら画像記録媒体の構造モトルに起因するものとすることができる。
【図面の簡単な説明】
【図１】本発明に係る画像読取装置の一例の概略を示す、正面図である。
【図２】本発明に係る画像読取方法において用いられる画像記録媒体への画像の記録方法を示す、概念図である。
【図３】本発明に係る画像読取装置の一例の概略を示す、正面図である。
【図４】本発明に係る画像読取方法において、第３の補正データ画像記録媒体の算出に用いるフィルタの特性を示す、グラフである。
【図５】本発明に係る画像読取方法、および従来の画像読取方法における、画像記録媒体の感度ムラ補正状況を示す、グラフである。
【図６】従来の画像読取方法における、画像記録媒体の感度ムラ補正方法を示すグラフである。
【符号の説明】
１        画像読取装置
１２      ＩＤ検出部
１６１    補正データ算出手段、画像データ補正手段（ＣＰＵ）
１６４    画像データ記録手段（ハードディスク）
２１      画像記録媒体（輝尽性蛍光体プレート）
【図１】
0004092932.tif 000002
【図２】
0004092932.tif 000003
【図３】
0004092932.tif 000004
【図４】
0004092932.tif 000005
【図５】
0004092932.tif 000006
【図６】
0004092932.tif 000007