JP4092932B2

JP4092932B2 - Image reading method and image reading apparatus

Info

Publication number: JP4092932B2
Application number: JP2002068964A
Authority: JP
Inventors: 智紀儀同
Original assignee: Konica Minolta Inc
Current assignee: Konica Minolta Inc
Priority date: 2002-03-13
Filing date: 2002-03-13
Publication date: 2008-05-28
Anticipated expiration: 2022-03-13
Also published as: JP2003270740A

Description

【０００１】
【発明の属する技術分野】
本発明は、画像情報方法および画像読取装置に係り、特に画像を記録する画像記録媒体の感度ムラを効果的に補正可能な画像読取方法および画像読取装置に関する。
【０００２】
【従来の技術】
画像が記録された記録媒体にレーザなどの光ビームを照射し、その記録媒体からの反射光、透過光あるいは発光光を検出することにより、画像を画像データに変換する画像読取方法や画像読取装置が知られている。このような画像読取装置として、製版用スキャナー、コンピュータやファクシミリの入力装置等のほか、輝尽性蛍光体プレートを使用した放射線画像情報記録再生システムにおいて用いられる放射線画像読取装置がある。
【０００３】
例えば、上記放射線画像読取装置において画像記録媒体として用いられる輝尽性蛍光体プレートで記録される画像には、集光系や光学系による画像ムラのほか、輝尽性蛍光体の２次元的な感度ムラが生じる。そこで、この感度ムラを補正するため、各輝尽性蛍光体プレートの感度ムラをあらかじめ求めておき、この輝尽性蛍光体プレートから読み取った画像データから補正データを加算するなどして、画像データを補正することが行われている。
【０００４】
ここで、輝尽性蛍光体プレートがカセッテ内に収納されたカセッテ方式の画像記録媒体では、カセッテを画像読取装置の読取部にセットする際のセット位置が、０．数ｍｍ〜数ｍｍ程度の範囲内で、カセッテをセットする毎に変動することがある。ここで、画像を読み取る際の精度を高くする場合には、この微小なセット位置のズレが、画像データの補正に影響を与える。
【０００５】
具体的には、被写体がない状態で画像を記録した輝尽性蛍光体プレートに形成された画像を、画像読取装置によりサンプリングピッチをＰで画像データに変換すると、この画像データは、図６（ａ）、図６（ｂ）中に太実線で示すように、上記感度ムラを包含したものとなる。被写体がない状態では画像データは一定値（図中、細実線）になるべきである。したがって、輝尽性蛍光体プレートに記録された被写体の画像を画像読み取り装置で読み込んで得られた画像データを補正するためには、画像データを太実線で除すれば良い。
【０００６】
しかし、輝尽性蛍光体プレートの画像読取装置媒体へのセット位置が、例えば１サンプリングピッチＰ分だけずれた状態（図中、太点線）で読み取られた画像データに対して上記方法により補正を施すと、補正後の画像データは図中に細点線で示すものとなる。このように、読取時のサンプリングピッチが小さい場合には、補正後の画像データの高周波のノイズが最大で補正前の画像データのノイズの２倍まで増幅されてしまう恐れがある。
【０００７】
そこで、画像読取時のサンプリングピッチに応じて、除去するノイズの周波数を考慮しつつ、低周波から高周波にわたるノイズを効果的に除去し、画像データの補正を良好に行うことが考えられる。
【０００８】
【発明が解決しようとする課題】
本発明の課題は、画像記録媒体に記録された画像を読み取る画像読取装置において、画像記録媒体の画像読取装置への装填位置の誤差があっても、画像記録媒体の感度ムラの大部分を補正可能な、画像読取方法、および画像読取装置を提供することである。
【０００９】
【課題を解決するための手段】
請求項１に記載の発明は、例えば図１〜図５に示すように、所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体（輝尽性蛍光体プレート）を、画像読取装置でサンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取方法において、
前記記録手段によって一様の入力がなされた前記画像記録媒体に対して前記画像読取装置で読取走査を行って前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得し、これを補正データ値Ｓc₀（ｘ，ｙ）とし、
前記画像読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc₀（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆ_xまたはｆ_y）に対応する周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とし、
前記記録手段により画像が記録された前記画像記録媒体に対して前記画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ₀（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出すること
を特徴とする。
【００１０】
ここで、上記画像記録媒体、記録手段とは、それぞれ、例えば放射線撮影装置において記録手段として用いられる輝尽性蛍光体プレート、Ｘ線のことである。
また、上記読取走査とは、例えば輝尽性蛍光体プレートをレーザ光などの励起光で走査して輝尽発光光を発生させ、得られた輝尽発光光を光電的に読み出して、画像信号を得ることである。
また、上記画像データ値がリニアな値である場合には、補正前画像データ値Ｓ₀（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を乗算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出する。上記画像データ値がＬｏｇ変換された値である場合には、補正前画像データ値Ｓ₀（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を和算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出する。
【００１１】
請求項１に記載の発明によれば、一様の入力がなされた画像記録媒体に対して画像読取装置でピッチＰで読取走査を行って補正データ値Ｓc₀（ｘ，ｙ）を取得し、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この画像記録媒体に記録された被写体の画像に対して画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ₀（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ₀（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００１２】
請求項２に記載の発明は、請求項１に記載の画像読取方法において、
前記補正データ値Ｓc₀（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆ_x，ｆ_y）の関数である第２の補正データ値Ｓc₂（ｆ_x，ｆ_y）とし、
Ｆ（ｆ_x，ｆ_y）＝１（０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆ_x，ｆ_y）＝０（ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆ_x，ｆ_y）を用いて、第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）
Ｓc₃（ｆ_x，ｆ_y）＝Ｆ（ｆ_x，ｆ_y）・Ｓc₂（ｆ_x，ｆ_y）、
とし、この第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）を逆フーリエ変換して、前記最終補正データＳc（ｘ，ｙ）とすること、
を特徴とする。
【００１３】
請求項２に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）の集合をフーリエ変換して第２の補正データ値Ｓc₂（ｆ_x，ｆ_y）とし、このうちから（１／（４ｎ・Ｐ））を超える周波数成分を除去して第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）とし、これを逆フーリエ変換して、最終補正データＳc（ｘ，ｙ）とするので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、確実に行える。
【００１４】
請求項３に記載の発明は、請求項１に記載の画像読取方法において、
前記補正データ値Ｓc₀（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆ_xまたはｆ_y）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする。
【００１５】
請求項３に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出するので、請求項２に記載の発明のようにフーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを行える。
【００１６】
請求項４に記載の発明は、請求項３に記載の画像読取方法において、
前記補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数を用いること
を特徴とする。
【００１７】
請求項４に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から得られる前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、空間周波数（ｆ_x，ｆ_y）領域で高周波成分を除去するゲート関数を逆フーリエ変換した結果得られるｓｉｎｃ関数を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行える。
【００１８】
請求項５に記載の発明は、請求項３に記載の画像読取方法において、
前記補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いること
を特徴とする。
ここで、移動平均法とは、読取座標（ｉ，ｊ）の近傍の画像データ値の平均を、読取座標（ｉ，ｊ）における画像データ値とする方法のことを言う。
【００１９】
請求項５に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、簡単な計算方法で確実に行える。
【００２０】
請求項６に記載の発明は、請求項１に記載の画像読取方法において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする。
【００２１】
請求項６に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、輝尽性蛍光体プレートを備える画像記録媒体に一様の入力を施す際に用いる放射線の照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc₀（ｘ，ｙ）を、もっぱら画像記録媒体の構造モトルに起因するものとすることができる。
【００２２】
請求項７に記載の発明は、例えば図１〜図５に示すように、所定の記録手段により画像が記録され、また記録された画像を消去して再び画像を記録することが可能な画像記録媒体を、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に読取走査し、前記画像を前記各読取座標（ｘ，ｙ）に対応する画像データ値に変換する画像読取装置において、
前記画像記録媒体のＩＤを検出可能なＩＤ検出部と、
前記画像記録媒体に対して読取走査を行って前記各読取座標に対応する画像データ値を取得し、これを画像データとして記録する画像データ記録手段と、
一様の入力がなされた前記画像記録媒体に対して、前記画像データ記録手段で読取走査を行い、前記各読取座標（ｘ，ｙ）に対応する画像データ値を取得して、これを補正データ値Ｓc₀（ｘ，ｙ）として記録し、
前記読取装置に前記画像記録媒体をセットする時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎについて、前記補正データ値Ｓc₀（ｘ，ｙ）のうち前記空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆ_xまたはｆ_y）に対応する周波数成分を除去し、これを最終補正データ値Ｓc（ｘ，ｙ）として、この画像記録媒体のＩＤと対応して記録する補正データ算出手段と、
前記補正前画像データ値Ｓ₀（ｘ，ｙ）を前記最終補正データ値Ｓc（ｘ，ｙ）により補正して補正後画像データ値Ｓ（ｘ，ｙ）を算出する画像データ補正手段と、を備えること
を特徴とする。
【００２３】
請求項７に記載の発明によれば、請求項１に記載の発明と同様の効果が得られる。
【００２４】
請求項８に記載の発明は、請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc₀（ｘ，ｙ）の集合をフーリエ変換し、空間周波数（ｆ_x，ｆ_y）の関数である第２の補正データ値Ｓc2（ｆ_x，ｆ_y）を算出し、
Ｆ（ｆ_x，ｆ_y）＝１（０＜ｆ≦１／（４ｎ・Ｐ）のとき）
Ｆ（ｆ_x，ｆ_y）＝０（ｆ＞１／（４ｎ・Ｐ）のとき）
である関数Ｆ（ｆ_x，ｆ_y）を用いて、第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）
Ｓc₃（ｆ_x，ｆ_y）＝Ｆ（ｆ_x，ｆ_y）・Ｓc₂（ｆ_x，ｆ_y）、
を算出し、この第３の補正データ値Ｓc3（ｆ_x，ｆ_y）を逆フーリエ変換して、空間（ｘ，ｙ）の関数である最終補正データＳc（ｘ，ｙ）を算出すること
を特徴とする。
【００２５】
請求項８に記載の発明によれば、請求項２に記載の発明と同様の効果が得られる。
【００２６】
請求項９に記載の発明は、請求項７に記載の画像読取装置において、
前記補正データ算出手段は、前記補正データ値Ｓc₀（ｘ，ｙ）のうち空間周波数ｆが、
ｆ＞１／（４ｎ・Ｐ）
である各空間周波数ｆ（ｆ_xまたはｆ_y）に対応する周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出すること
を特徴とする。
【００２７】
請求項９に記載の発明によれば、請求項３に記載の発明と同様の効果が得られる。
【００２８】
請求項１０に記載の発明は、請求項９に記載の画像読取装置において、
前記補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、ｓｉｎｃ関数が用いられていること
を特徴とする。
【００２９】
請求項１０に記載の発明によれば、請求項４に記載の発明と同様の効果が得られる。
【００３０】
請求項１１に記載の発明は、請求項９に記載の画像読取装置において、
前記補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算が用いられていること
を特徴とする。
【００３１】
請求項１１に記載の発明によれば、請求項５に記載の発明と同様の効果が得られる。
【００３２】
請求項１２に記載の発明は、請求項７〜１１のいずれかに記載の画像読取装置において、
前記記録手段は放射線であり、前記画像記録媒体は放射線を蓄積記録する輝尽性蛍光体プレートを備えるものであって、
この画像記録媒体に、照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上の放射線で前記一様の入力を施すこと
を特徴とする。
【００３３】
請求項１２に記載の発明によれば、請求項６に記載の発明と同様の効果が得られる。
【００３４】
【発明の実施の形態】
以下、本発明に係る画像読取方法、および画像読取装置について、図１〜図５を参照して説明する。
本実施の形態の画像読取装置１は、放射線画像記録読取システムにおいて、撮影台３に装填された状態で画像が撮影記録される輝尽性蛍光体プレート（画像記録媒体）２１の画像を読み取り、画像データに変換する放射線画像読取装置である。
輝尽性蛍光体プレート２１は、図２に示すように、カセッテ２に収容されており、撮影台３にセットされた状態でＸ線発生装置４から出射されたＸ線４１のうち人体などの被写体５およびその周囲を通過するＸ線の照射エネルギーを蓄積することで、被写体５のＸ線画像を記録するものである。この輝尽性蛍光体プレート２１に、可視光線などの励起光を照射すると、蓄積されたエネルギーが解放され、輝尽発光が生じる。この現象を利用して、輝尽性蛍光体プレート２１に記録された画像の読み取りや消去を画像読取装置１で行う。
【００３５】
画像読取装置１は、図１、図３に示すように、上記カセッテ２に収容された輝尽性蛍光体プレート２１を引き出し、サンプリングピッチＰの読取座標（ｘ，ｙ）毎にレーザ光などの励起光で輝尽発光させつつ読取走査し、輝尽性蛍光体プレート２１に記録された画像を各読取座標（ｘ，ｙ）に対応する画像データ値に変換するものである。
【００３６】
画像記録装置１には、カセッテ２を装填するためのスタッカ部１１ａ〜１１ｅが備えられ、これら各スタッカ１１ａ〜１１ｅには、装填されたカセッテ２のＩＤを検出するためのＩＤ検出部１２が設けられている。
そして、表示・操作部１７のＣＲＴ表示部１７１前面のタッチパネル１７２を操作すると、その信号がＣＰＵ１６１に伝達される。そして、ＣＰＵ１６１の命令が、読取部・画像入力制御部１６２、メカ制御部１６３を介してカセッテ搬送機構１３に伝達され、このカセッテ搬送機構１３の動作により、スタッカ部１１ａ〜１１ｅに装填されたカセッテ２のうち所望のカセッテ２が、スタッカ部１１ａ〜１１ｅとカセッテ２に記録された画像の読み取りを行う読取位置との間を搬送される。
また、同タッチパネル１７２の操作により、ＣＰＵ１６１の命令が、読取部・画像入力制御部１６２を介してカセッテ搬送機構１３に設けられた消去部１４に伝達され、消去部１４から可視光線がカセッテ２の記録面全面に照射されて、カセッテ２に記録された画像を消去できるようになっている。
【００３７】
タッチパネル１７２で、所望のカセッテ２の読み取りを指定した場合には、そのカセッテ２に収容された輝尽性蛍光体プレート２１が、カセッテ搬送機構１３の動作によって上記読取位置に搬送され、画像読取部１５からレーザ光からなる励起光を照射されて、この励起光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査が行われる。この読取走査によって、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値が取得される。
【００３８】
この画像データ値は、システム制御部１６に備えられたＣＰＵ１６１の命令により、画像データとしてハードディスク（画像データ記録手段）１６４に記録される。
ハードディスク１６４には、後述する補正データ算出処理や、画像データ補正処理を行うためのプログラムや、各輝尽性蛍光体プレート２１についての最終補正データ（後述）も記録されている。これらの画像データやプログラムは、必要に応じてＣＰＵ１６１から読み出される。ＣＰＵ１６１による処理で扱われるデータやプログラムは、一時的にＲＡＭ（Read-Only Memory）１６５に記憶される。
【００３９】
また、ＣＰＵ１６１は、インターフェイス１６６を介して、ＬＡＮ９に接続され、このＬＡＮ９に接続されている、各患者の情報を管理するホストコンピュータ６や、画像データを出力するためのハードコピー装置７、患者登録ターミナル８などとの間で、各種データの送受信を行えるようになっている。
【００４０】
この画像読取装置１を利用した、輝尽性蛍光体プレート２１に記録された画像の読取方法について説明する。ここでは、画像データ値がリニアな値である場合について説明する。
まず、図２に示すように、輝尽性蛍光体プレート２１を収容するカセッテ２を撮影台３にセットし、この状態で撮影台の前に被写体５を置かずに、Ｘ線発生装置４により照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上、好ましくは２．５×１０^-5［Ｃ／ｋｇ］の放射線を撮影台３に照射して、このＸ線４１でカセッテ２内の輝尽性蛍光体プレート２１を励起させる。
次に、このカセッテ２を撮影台３から取り出して、画像読取装置１のスタッカ部１１ａ〜１１ｅのいずれかに装填すると、カセッテ２のＩＤがＩＤ検出部１２によって読み出される。そして、画像読取部１５からレーザ光からなる励起光を照射させて、この励起光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査を行う。この読取走査によって、カセット２サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値を取得し、これが補正データ値Ｓc₀（ｘ，ｙ）として、輝尽性蛍光体プレート２１が収容されていたカセッテ２のＩＤと対応してハードディスク１６４に記録される。
【００４１】
この第１の補正データ値は、表示・操作部１７のタッチパネル１７２を操作することにより、ＣＰＵ１６１によってハードディスク１６４から読み出された二次元ＦＦＴ（Fast Fourier Transformation：高速フーリエ変換）プログラムによりフーリエ変換が施され、各空間周波数に対応する画像データ値が出力される。この画像データ値がＲＡＭ１６５に第２の補正データ値Ｓc₂（ｆ_x，ｆ_y）として一時的に記録される。
【００４２】
画像読取装置１の上記読取位置に輝尽性蛍光体プレート２１がセットされる時の位置ズレ量の絶対値の最大値をｄとし、ｎを自然数として、
（ｎ−１）Ｐ＜ｄ≦ｎ・Ｐ
となるｎを、あらかじめ設定しておく。
そして、ＣＰＵ１６１により、第２の補正データ値Ｓc₂（ｆ_x，ｆ_y）のうち空間周波数ｆがｆc＝１／（４ｎ・Ｐ）を超える各空間周波数ｆ（ｆ_xまたはｆ_y）に対応する前記画像データ値を０に置換し、その結果が第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）としてＲＡＭ１６５に一時的に記録される。
【００４３】
続いて、ＣＰＵ１６１によって、第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）に二次元逆フーリエ変換が施され、再び各読取座標（ｘ，ｙ）に対応する画像データ値が出力される。この画像データ値をその最大値に対する比として規準化したものが、最終補正データ値Ｓc（ｆ_x，ｆ_y）として、ハードディスク１６４に記録される。すなわち、最終補正データ値Ｓc（ｆ_x，ｆ_y）は、０＜Ｓc（ｆ_x，ｆ_y）≦１となる。そして、この輝尽性蛍光体プレート２１に記録された画像を、カセッテ搬送機構１３の消去部１４で消去し、再びカセッテ２に収容し、このカセッテ２を装填したスタッカ１１ａ〜１１ｅのいずれかから取り出す。
以上で、補正データ算出処理が終了する。
【００４４】
次いで、図２に示すように、このカセッテ２を撮影台３に装填し、撮影台３の全面に人体などの被写体５を配置した状態で、その背後よりＸ線発生手段でＸ線を照射し、カセッテ２の輝尽性蛍光体２１を励起させ、被写体５の画像を記録する。
【００４５】
次に、このカセッテ２を撮影台３から取り出して、画像読取装置１のスタッカ部１１ａ〜１１ｅのいずれかに装填すると、カセッテ２のＩＤがＩＤ検出部１２によって読み出される。そして、画像読取部１５からレーザ光からなる励起光を輝尽性蛍光体プレート２１に照射して、この励起光による輝尽発光を画像読取部１５に設けられたセンサ（図示せず）で受光する読取走査を行う。この読取走査によって、サンプリングピッチＰの読取座標（ｘ，ｙ）毎に、各読取座標（ｘ，ｙ）に対応する画像データ値が取得され、これが補正前画像データ値Ｓ₀（ｘ，ｙ）として、輝尽性蛍光体プレート２１が収容されていたカセッテ２のＩＤと対応してＲＡＭ１６５またはハードディスク１６４に記録される。
【００４６】
そして、補正前画像データ値Ｓ₀（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を乗算して、この結果が補正後画像データ値Ｓ（ｘ，ｙ）としてハードディスク１６４に記録される。
以上で、画像データ補正処理が終了する。
【００４７】
本実施の形態に記載の画像読取方法、および画像読取装置１によれば、一様の入力がなされた輝尽性蛍光体プレート２１に対して画像読取装置１でピッチＰで読取走査を行って補正データ値Ｓc₀（ｘ，ｙ）を取得し、画像読取装置１にカセッテ２に収容された輝尽性蛍光体プレート２１をセットする時にズレる可能性のある最大画素数ｎに応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この輝尽性蛍光体プレート２１に記録された被写体の画像に対して画像読取装置１で読取走査を行って取得される補正前画像データ値Ｓ₀（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ₀（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００４８】
また、輝尽性蛍光体プレート２１を備えるカセッテ２に一様の入力を施す際に用いる放射線の照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc₀（ｘ，ｙ）を、もっぱら輝尽性蛍光体プレート２１の構造モトルに起因するものとすることができる。
したがって、カセッテ２の輝尽性蛍光体プレート２１の構造モトルによるノイズを効果的に除去することができる。
【００４９】
なお、本発明の画像読取方法、および画像読取装置は、上記の実施の形態に限定されることなく、本発明の趣旨を逸脱しない範囲において、種々の改良ならびに設計の変更を行ってもよい。
例えば、上記実施の形態においては、補正データ値Ｓc₀（ｘ，ｙ）から最終補正データ値Ｓc（ｘ，ｙ）を算出する方法として、フーリエ変換および逆フーリエ変換を用いる場合について説明したが、ｓｉｎｃ関数や移動平均法を用いて、空間（ｘ，ｙ）領域において補正データ値Ｓc₀（ｘ，ｙ）から最終補正データ値Ｓc（ｘ，ｙ）を得るようにすれば、フーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行うことができる。
また、上記画像データ値がＬｏｇ変換された値である場合には、高周波成分を除去した補正データ値Ｓc₀（ｘ，ｙ）を、その最大値に対する差分で規準化して最終補正データ値Ｓc（ｘ，ｙ）とし、補正前画像データ値Ｓ₀（ｘ，ｙ）に、前記最終補正データ値Ｓc（ｘ，ｙ）を和算して、補正後画像データ値Ｓ（ｘ，ｙ）を算出すれば良い。
また、上記実施の形態においては、放射線画像読取装置について説明したが、本発明は一般のイメージスキャナなどにも適用可能である。
その他、具体的な細部構造などについても適宜に変更可能であることももちろんである。
【００５０】
【発明の効果】
請求項１または請求項７に記載の発明によれば、一様の入力がなされた画像記録媒体に対して画像読取装置でピッチＰで読取走査を行って補正データ値Ｓc₀（ｘ，ｙ）を取得し、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）としておき、この画像記録媒体に記録された被写体の画像に対して画像読取装置で読取走査を行って取得される補正前画像データ値Ｓ₀（ｘ，ｙ）を、前記最終補正データ値Ｓc（ｘ，ｙ）により補正して、補正後画像データ値Ｓ（ｘ，ｙ）を算出するので、補正前画像データ値Ｓ₀（ｘ，ｙ）のうち高周波のノイズが増幅されることなく、低周波から高周波までにわたる二次元的ノイズが効果的に除去された画像データ値を得ることができる。
【００５１】
請求項２または請求項８に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）の集合をフーリエ変換して第２の補正データ値Ｓc₂（ｆ_x，ｆ_y）とし、このうちから（１／（４ｎ・Ｐ））を超える周波数成分を除去して第３の補正データ値Ｓc₃（ｆ_x，ｆ_y）とし、これを逆フーリエ変換して、最終補正データＳc（ｘ，ｙ）とするので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、確実に行える。
【００５２】
請求項３または請求項９に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を、（ｘ，ｙ）空間領域における演算で除去し、前記最終補正データ値Ｓc（ｘ，ｙ）を算出するので、請求項２に記載の発明のようにフーリエ変換および逆フーリエ変換を行わないので、比較的に少ない計算量で、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを行える。
【００５３】
請求項４または請求項１０に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、空間周波数（ｆ_x，ｆ_y）領域で高周波成分を除去するゲート関数を逆フーリエ変換した結果得られるｓｉｎｃ関数を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、少ない計算量で確実に行える。
【００５４】
請求項５または請求項１１に記載の発明によれば、請求項３に記載の発明と同様の効果が得られるとともに、補正データ値Ｓc₀（ｘ，ｙ）から前記最終補正データ値Ｓc（ｘ，ｙ）を算出する際に、移動平均法による演算を用いるので、画像読取装置に画像記録媒体をセットする時にズレる可能性のある最大画素数に応じて、補正データ値Ｓc₀（ｘ，ｙ）から（１／（４ｎ・Ｐ））を超える周波数成分を除去してこれを最終補正データ値Ｓc（ｘ，ｙ）とすることを、簡単な計算方法で確実に行える。
【００５５】
請求項６または請求項１２に記載の発明によれば、請求項１に記載の発明と同様の効果が得られるとともに、輝尽性蛍光体プレートを備える画像記録媒体に一様の入力を施す際に用いる放射線の照射線量が２．５×１０^-6［Ｃ／ｋｇ］以上であるので、補正データ値Ｓc₀（ｘ，ｙ）を、もっぱら画像記録媒体の構造モトルに起因するものとすることができる。
【図面の簡単な説明】
【図１】本発明に係る画像読取装置の一例の概略を示す、正面図である。
【図２】本発明に係る画像読取方法において用いられる画像記録媒体への画像の記録方法を示す、概念図である。
【図３】本発明に係る画像読取装置の一例の概略を示す、正面図である。
【図４】本発明に係る画像読取方法において、第３の補正データ画像記録媒体の算出に用いるフィルタの特性を示す、グラフである。
【図５】本発明に係る画像読取方法、および従来の画像読取方法における、画像記録媒体の感度ムラ補正状況を示す、グラフである。
【図６】従来の画像読取方法における、画像記録媒体の感度ムラ補正方法を示すグラフである。
【符号の説明】
１画像読取装置
１２ＩＤ検出部
１６１補正データ算出手段、画像データ補正手段（ＣＰＵ）
１６４画像データ記録手段（ハードディスク）
２１画像記録媒体（輝尽性蛍光体プレート）[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an image information method and an image reading apparatus, and more particularly to an image reading method and an image reading apparatus capable of effectively correcting sensitivity unevenness of an image recording medium for recording an image.
[0002]
[Prior art]
An image reading method or image reading apparatus that converts an image into image data by irradiating a recording medium on which an image is recorded with a light beam such as a laser and detecting reflected light, transmitted light, or emitted light from the recording medium. It has been known. As such an image reading apparatus, there is a radiographic image reading apparatus used in a radiographic image information recording / reproducing system using a photostimulable phosphor plate, in addition to a plate-making scanner, a computer or facsimile input device, and the like.
[0003]
For example, an image recorded by a stimulable phosphor plate used as an image recording medium in the radiation image reading apparatus includes, in addition to image unevenness due to a condensing system or an optical system, a two-dimensional stimulable phosphor. Sensitivity unevenness occurs. Therefore, in order to correct this sensitivity unevenness, the sensitivity unevenness of each stimulable phosphor plate is obtained in advance, and the correction data is added from the image data read from this stimulable phosphor plate. It has been done to correct.
[0004]
Here, in the cassette type image recording medium in which the photostimulable phosphor plate is accommodated in the cassette, the set position when the cassette is set in the reading unit of the image reading apparatus is 0. Within a range of several millimeters to several millimeters, it may fluctuate every time the cassette is set. Here, in the case where the accuracy in reading an image is increased, this small set position deviation affects the correction of the image data.
[0005]
Specifically, when an image formed on a photostimulable phosphor plate on which an image is recorded in the absence of a subject is converted into image data with a sampling pitch P by an image reader, the image data is shown in FIG. a) As shown by a thick solid line in FIG. 6B, the sensitivity unevenness is included. When there is no subject, the image data should be a constant value (thin solid line in the figure). Therefore, in order to correct the image data obtained by reading the image of the subject recorded on the photostimulable phosphor plate by the image reading device, the image data may be divided by a thick solid line.
[0006]
However, the above method is used to correct image data read in a state where the setting position of the photostimulable phosphor plate on the image reader medium is shifted by, for example, one sampling pitch P (thick dotted line in the figure). When applied, the corrected image data is indicated by a thin dotted line in the figure. Thus, when the sampling pitch at the time of reading is small, there is a possibility that the high-frequency noise of the image data after correction is amplified up to twice the noise of the image data before correction.
[0007]
Therefore, it is conceivable that image data is corrected well by effectively removing noise from low to high frequencies while considering the frequency of noise to be removed in accordance with the sampling pitch at the time of image reading.
[0008]
[Problems to be solved by the invention]
An object of the present invention is to correct most of the sensitivity unevenness of an image recording medium in an image reading apparatus that reads an image recorded on the image recording medium even if there is an error in the loading position of the image recording medium to the image reading apparatus. An object is to provide an image reading method and an image reading apparatus.
[0009]
[Means for Solving the Problems]
According to the first aspect of the present invention, for example, as shown in FIGS. 1 to 5, an image can be recorded by a predetermined recording means, and the recorded image can be erased and recorded again. The medium (stimulable phosphor plate) is read and scanned by the image reading device at each reading coordinate (x, y) of the sampling pitch P, and the image data value corresponding to each reading coordinate (x, y) is read. In the image reading method of converting to
The image reading medium that is uniformly input by the recording means is scanned by the image reading device to obtain image data values corresponding to the respective reading coordinates (x, y), and this is corrected. Data value Sc ₀ (X, y)
The maximum absolute value of the positional deviation amount when setting the image recording medium in the image reading device is d, and n is a natural number.
(N-1) P <d ≦ n · P
The correction data value Sc for n ₀ Of (x, y), the spatial frequency f is
f> 1 / (4n · P)
Each spatial frequency f (f _x Or f _y ) Is removed and this is used as the final correction data value Sc (x, y),
An uncorrected image data value S obtained by scanning the image recording medium on which the image is recorded by the recording means with the image reading device. ₀ (X, y) is corrected by the final correction data value Sc (x, y) to calculate a corrected image data value S (x, y).
It is characterized by.
[0010]
Here, the image recording medium and the recording means are, for example, a stimulable phosphor plate and an X-ray used as recording means in a radiographic apparatus, for example.
Further, the reading scanning is, for example, scanning the photostimulable phosphor plate with excitation light such as laser light to generate stimulated emission light, photoelectrically reading out the obtained stimulated emission light, and image signal. Is to get.
If the image data value is a linear value, the pre-correction image data value S ₀ (X, y) is multiplied by the final corrected data value Sc (x, y) to calculate a corrected image data value S (x, y). When the image data value is a value subjected to Log conversion, the uncorrected image data value S ₀ The final corrected data value Sc (x, y) is added to (x, y) to calculate a corrected image data value S (x, y).
[0011]
According to the first aspect of the present invention, the correction data value Sc is obtained by reading and scanning the image recording medium with uniform input at the pitch P by the image reading apparatus. ₀ (X, y) is acquired, and the correction data value Sc according to the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. ₀ A frequency component exceeding (1 / (4n · P)) is removed from (x, y) and set as a final correction data value Sc (x, y), and the image of the subject recorded on the image recording medium is set. On the other hand, the pre-correction image data value S obtained by performing scanning scanning with the image reading apparatus. ₀ (X, y) is corrected by the final correction data value Sc (x, y), and the corrected image data value S (x, y) is calculated. ₀ It is possible to obtain an image data value from which two-dimensional noise ranging from low frequency to high frequency is effectively removed without amplifying high frequency noise in (x, y).
[0012]
According to a second aspect of the present invention, in the image reading method of the first aspect,
The correction data value Sc ₀ The set of (x, y) is Fourier transformed and the spatial frequency (f _x , F _y ) As a second correction data value Sc ₂ (F _x , F _y )age,
F (f _x , F _y ) = 1 (when 0 <f ≦ 1 / (4n · P))
F (f _x , F _y ) = 0 (when f> 1 / (4n · P))
Is a function F (f _x , F _y ) To obtain the third correction data value Sc _Three (F _x , F _y )
Sc _Three (F _x , F _y ) = F (f _x , F _y ) ・ Sc ₂ (F _x , F _y ),
And this third correction data value Sc _Three (F _x , F _y ) Is subjected to inverse Fourier transform to obtain the final correction data Sc (x, y),
It is characterized by.
[0013]
According to the second aspect of the invention, the same effect as that of the first aspect of the invention can be obtained, and the correction data value Sc ₀ The set of (x, y) is Fourier transformed to obtain the second correction data value Sc ₂ (F _x , F _y ), And a frequency component exceeding (1 / (4n · P)) is removed from the third correction data value Sc _Three (F _x , F _y ), And this is subjected to inverse Fourier transform to obtain final correction data Sc (x, y). Therefore, the correction data value depends on the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. Sc ₀ It is possible to reliably remove the frequency component exceeding (1 / (4n · P)) from (x, y) and set it as the final correction data value Sc (x, y).
[0014]
The invention described in claim 3 is the image reading method according to claim 1,
The correction data value Sc ₀ Of (x, y), the spatial frequency f is
f> 1 / (4n · P)
Each spatial frequency f (f _x Or f _y ) Is removed by computation in the (x, y) space region, and the final correction data value Sc (x, y) is calculated.
It is characterized by.
[0015]
According to the third aspect of the invention, the same effect as that of the first aspect of the invention can be obtained, and the correction data value Sc ₀ Since the frequency component exceeding (1 / (4n · P)) from (x, y) is removed by calculation in the (x, y) space region, the final correction data value Sc (x, y) is calculated. Since the Fourier transform and the inverse Fourier transform are not performed as in the second aspect of the invention, the maximum amount of pixels that may be shifted when the image recording medium is set in the image reading apparatus with a relatively small amount of calculation. Correction data value Sc ₀ It is possible to remove frequency components exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y).
[0016]
The invention according to claim 4 is the image reading method according to claim 3,
The correction data value Sc ₀ Use the sinc function when calculating the final correction data value Sc (x, y) from (x, y).
It is characterized by.
[0017]
According to the invention described in claim 4, the same effect as that of the invention described in claim 3 can be obtained, and the correction data value Sc ₀ When calculating the final correction data value Sc (x, y) obtained from (x, y), the spatial frequency (f _x , F _y ) Since the sinc function obtained as a result of inverse Fourier transform of the gate function for removing high frequency components in the region is used, the correction data value is determined according to the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. Sc ₀ It is possible to reliably remove a frequency component exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y) with a small amount of calculation.
[0018]
The invention according to claim 5 is the image reading method according to claim 3,
The correction data value Sc ₀ When calculating the final correction data value Sc (x, y) from (x, y), use an operation by the moving average method.
It is characterized by.
Here, the moving average method refers to a method in which an average of image data values in the vicinity of the read coordinates (i, j) is used as an image data value at the read coordinates (i, j).
[0019]
According to the fifth aspect of the invention, the same effect as that of the third aspect of the invention can be obtained, and the correction data value Sc ₀ When calculating the final correction data value Sc (x, y) from (x, y), the calculation by the moving average method is used, so that the maximum pixel that may be shifted when the image recording medium is set in the image reading apparatus. Depending on the number, the correction data value Sc ₀ It is possible to reliably remove a frequency component exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y) by a simple calculation method.
[0020]
According to a sixth aspect of the present invention, in the image reading method of the first aspect,
The recording means is radiation, and the image recording medium includes a stimulable phosphor plate for accumulating and recording radiation,
This image recording medium has an irradiation dose of 2.5 × 10 ^-6 Applying the uniform input with radiation of [C / kg] or more
It is characterized by.
[0021]
According to the invention described in claim 6, the same effect as that of the invention described in claim 1 can be obtained, and the radiation used for uniform input to the image recording medium including the stimulable phosphor plate can be obtained. Irradiation dose is 2.5 × 10 ^-6 Since it is equal to or greater than [C / kg], the correction data value Sc ₀ (X, y) can be attributed solely to the structure of the image recording medium.
[0022]
According to the seventh aspect of the present invention, for example, as shown in FIGS. 1 to 5, an image can be recorded by a predetermined recording means, and the recorded image can be erased and recorded again. In an image reading apparatus that scans and scans a medium for each reading coordinate (x, y) at a sampling pitch P and converts the image into an image data value corresponding to each reading coordinate (x, y).
An ID detection unit capable of detecting an ID of the image recording medium;
Image data recording means for scanning the image recording medium to obtain image data values corresponding to the respective reading coordinates, and recording the image data values as image data;
The image recording medium on which uniform input has been made is read and scanned by the image data recording means to obtain image data values corresponding to the respective reading coordinates (x, y), and this is corrected data. Value Sc ₀ Record as (x, y)
Assuming that the maximum absolute value of the positional deviation amount when setting the image recording medium in the reading device is d, n is a natural number,
(N-1) P <d ≦ n · P
The correction data value Sc for n ₀ The spatial frequency f of (x, y) is
f> 1 / (4n · P)
Each spatial frequency f (f _x Or f _y ), And the correction data calculation means for recording this as the final correction data value Sc (x, y) corresponding to the ID of the image recording medium;
Image data value S before correction ₀ Image data correction means for correcting (x, y) with the final correction data value Sc (x, y) and calculating a corrected image data value S (x, y).
It is characterized by.
[0023]
According to the invention described in claim 7, the same effect as that of the invention described in claim 1 can be obtained.
[0024]
The invention according to claim 8 is the image reading apparatus according to claim 7,
The correction data calculation means is the correction data value Sc ₀ The set of (x, y) is Fourier transformed and the spatial frequency (f _x , F _y ) As a second correction data value Sc2 (f _x , F _y )
F (f _x , F _y ) = 1 (when 0 <f ≦ 1 / (4n · P))
F (f _x , F _y ) = 0 (when f> 1 / (4n · P))
Is a function F (f _x , F _y ) To obtain the third correction data value Sc _Three (F _x , F _y )
Sc _Three (F _x , F _y ) = F (f _x , F _y ) ・ Sc ₂ (F _x , F _y ),
And the third correction data value Sc3 (f _x , F _y ) Is subjected to inverse Fourier transform to calculate final correction data Sc (x, y) which is a function of space (x, y).
It is characterized by.
[0025]
According to the eighth aspect of the invention, the same effect as that of the second aspect of the invention can be obtained.
[0026]
The invention according to claim 9 is the image reading apparatus according to claim 7,
The correction data calculation means is the correction data value Sc ₀ Of (x, y), the spatial frequency f is
f> 1 / (4n · P)
Each spatial frequency f (f _x Or f _y ) Is removed by computation in the (x, y) space region, and the final correction data value Sc (x, y) is calculated.
It is characterized by.
[0027]
According to the ninth aspect of the invention, the same effect as that of the third aspect of the invention can be obtained.
[0028]
According to a tenth aspect of the present invention, in the image reading apparatus according to the ninth aspect,
The correction data value Sc ₀ A sinc function is used when calculating the final correction data value Sc (x, y) from (x, y).
It is characterized by.
[0029]
According to the invention described in claim 10, the same effect as that of the invention described in claim 4 can be obtained.
[0030]
The invention according to claim 11 is the image reading apparatus according to claim 9, wherein
The correction data value Sc ₀ When calculating the final correction data value Sc (x, y) from (x, y), an arithmetic operation using a moving average method is used.
It is characterized by.
[0031]
According to the eleventh aspect of the invention, the same effect as that of the fifth aspect of the invention can be obtained.
[0032]
The invention according to claim 12 is the image reading apparatus according to any one of claims 7 to 11,
The recording means is radiation, and the image recording medium includes a stimulable phosphor plate for accumulating and recording radiation,
This image recording medium has an irradiation dose of 2.5 × 10 ^-6 Applying the uniform input with radiation of [C / kg] or more
It is characterized by.
[0033]
According to the twelfth aspect, the same effect as that of the sixth aspect can be obtained.
[0034]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, an image reading method and an image reading apparatus according to the present invention will be described with reference to FIGS.
The image reading apparatus 1 according to the present embodiment reads an image of a photostimulable phosphor plate (image recording medium) 21 on which an image is photographed and recorded in the radiation image recording and reading system in a state of being loaded on the photographing table 3. It is a radiographic image reading apparatus which converts into image data.
As shown in FIG. 2, the photostimulable phosphor plate 21 is accommodated in the cassette 2, and the X-ray 41 emitted from the X-ray generator 4 in a state of being set on the imaging table 3 is a human body or the like. An X-ray image of the subject 5 is recorded by accumulating X-ray irradiation energy passing through the subject 5 and its surroundings. When this photostimulable phosphor plate 21 is irradiated with excitation light such as visible light, the accumulated energy is released and stimulated light emission occurs. Using this phenomenon, the image reading apparatus 1 reads and erases the image recorded on the photostimulable phosphor plate 21.
[0035]
As shown in FIGS. 1 and 3, the image reading apparatus 1 pulls out the stimulable phosphor plate 21 accommodated in the cassette 2, and uses laser light or the like for each reading coordinate (x, y) of the sampling pitch P. Scanning is performed while stimulating light is emitted with excitation light, and an image recorded on the stimulable phosphor plate 21 is converted into image data values corresponding to the respective reading coordinates (x, y).
[0036]
The image recording apparatus 1 includes stackers 11 a to 11 e for loading the cassette 2, and each stacker 11 a to 11 e is provided with an ID detection unit 12 for detecting the ID of the loaded cassette 2. It has been.
When the touch panel 172 on the front surface of the CRT display unit 171 of the display / operation unit 17 is operated, the signal is transmitted to the CPU 161. The instruction of the CPU 161 is transmitted to the cassette transport mechanism 13 via the reading unit / image input control unit 162 and the mechanical control unit 163, and the cassettes loaded in the stacker units 11 a to 11 e are operated by the operation of the cassette transport mechanism 13. 2, a desired cassette 2 is conveyed between the stacker units 11 a to 11 e and a reading position for reading an image recorded on the cassette 2.
Also, by operating the touch panel 172, a command from the CPU 161 is transmitted to the erasing unit 14 provided in the cassette transport mechanism 13 via the reading unit / image input control unit 162, and visible light is transmitted from the erasing unit 14 to the cassette 2. The image recorded on the cassette 2 can be erased by irradiating the entire recording surface.
[0037]
When reading of the desired cassette 2 is designated on the touch panel 172, the photostimulable phosphor plate 21 accommodated in the cassette 2 is conveyed to the reading position by the operation of the cassette conveying mechanism 13, and the image reading unit 15 is irradiated with excitation light composed of laser light, and reading scanning is performed in which the excitation light is received by a sensor (not shown) provided in the image reading unit 15. By this reading scanning, for each reading coordinate (x, y) of the sampling pitch P, an image data value corresponding to each reading coordinate (x, y) is acquired.
[0038]
This image data value is recorded on the hard disk (image data recording means) 164 as image data in accordance with a command from the CPU 161 provided in the system control unit 16.
The hard disk 164 also stores a correction data calculation process (to be described later), a program for performing an image data correction process, and final correction data (described later) for each stimulable phosphor plate 21. These image data and programs are read from the CPU 161 as necessary. Data and programs handled in the processing by the CPU 161 are temporarily stored in a RAM (Read-Only Memory) 165.
[0039]
The CPU 161 is connected to the LAN 9 via the interface 166. The host computer 6 that manages information of each patient connected to the LAN 9, the hard copy device 7 for outputting image data, and patient registration Various data can be transmitted to and received from the terminal 8 or the like.
[0040]
A method of reading an image recorded on the photostimulable phosphor plate 21 using the image reading apparatus 1 will be described. Here, a case where the image data value is a linear value will be described.
First, as shown in FIG. 2, the cassette 2 containing the photostimulable phosphor plate 21 is set on the imaging table 3, and the X-ray generator 4 does not place the subject 5 in front of the imaging table in this state. Irradiation dose is 2.5 × 10 ^-6 [C / kg] or more, preferably 2.5 × 10 ^-Five Radiation of [C / kg] is applied to the imaging table 3, and the stimulable phosphor plate 21 in the cassette 2 is excited by this X-ray 41.
Next, when the cassette 2 is taken out from the imaging table 3 and loaded into any of the stacker units 11 a to 11 e of the image reading apparatus 1, the ID of the cassette 2 is read by the ID detection unit 12. Then, reading scanning is performed in which excitation light composed of laser light is irradiated from the image reading unit 15 and the excitation light is received by a sensor (not shown) provided in the image reading unit 15. By this scanning, an image data value corresponding to each reading coordinate (x, y) is obtained for each reading coordinate (x, y) of the cassette 2 sampling pitch P, and this is the correction data value Sc. ₀ (X, y) is recorded on the hard disk 164 corresponding to the ID of the cassette 2 in which the photostimulable phosphor plate 21 was accommodated.
[0041]
The first correction data value is subjected to Fourier transform by a two-dimensional FFT (Fast Fourier Transformation) program read from the hard disk 164 by the CPU 161 by operating the touch panel 172 of the display / operation unit 17. Then, the image data value corresponding to each spatial frequency is output. This image data value is stored in the RAM 165 as the second correction data value Sc. ₂ (F _x , F _y ) Temporarily recorded.
[0042]
When the maximum value of the absolute value of the positional deviation amount when the photostimulable phosphor plate 21 is set at the reading position of the image reading apparatus 1 is set as d, and n is a natural number,
(N-1) P <d ≦ n · P
N is set in advance.
Then, by the CPU 161, the second correction data value Sc ₂ (F _x , F _y ), Each spatial frequency f (f) where the spatial frequency f exceeds fc = 1 / (4n · P) _x Or f _y ) Is replaced with 0, and the result is the third correction data value Sc. _Three (F _x , F _y ) Is temporarily recorded in the RAM 165.
[0043]
Subsequently, the CPU 161 makes a third correction data value Sc. _Three (F _x , F _y ) Is subjected to two-dimensional inverse Fourier transform, and image data values corresponding to the respective read coordinates (x, y) are output again. This image data value normalized as a ratio to its maximum value is the final corrected data value Sc (f _x , F _y ) Is recorded on the hard disk 164. That is, the final correction data value Sc (f _x , F _y ) 0 <Sc (f _x , F _y ) ≦ 1. Then, the image recorded on the photostimulable phosphor plate 21 is erased by the erasing unit 14 of the cassette transport mechanism 13, stored again in the cassette 2, and from any one of the stackers 11 a to 11 e loaded with the cassette 2. Take out.
This is the end of the correction data calculation process.
[0044]
Next, as shown in FIG. 2, the cassette 2 is loaded on the imaging table 3, and a subject 5 such as a human body is placed on the entire surface of the imaging table 3, and X-rays are emitted from the back by X-ray generation means. The photostimulable phosphor 21 of the cassette 2 is excited and an image of the subject 5 is recorded.
[0045]
Next, when the cassette 2 is taken out from the imaging table 3 and loaded into any of the stacker units 11 a to 11 e of the image reading apparatus 1, the ID of the cassette 2 is read by the ID detection unit 12. Then, the stimulable phosphor plate 21 is irradiated with excitation light composed of laser light from the image reading unit 15, and stimulated light emission by the excitation light is received by a sensor (not shown) provided in the image reading unit 15. Read scanning is performed. By this reading scanning, image data values corresponding to the respective reading coordinates (x, y) are obtained for each reading coordinate (x, y) of the sampling pitch P, and this is the image data value S before correction. ₀ (X, y) is recorded in the RAM 165 or the hard disk 164 corresponding to the ID of the cassette 2 in which the photostimulable phosphor plate 21 was accommodated.
[0046]
The uncorrected image data value S ₀ (X, y) is multiplied by the final correction data value Sc (x, y), and this result is recorded on the hard disk 164 as a corrected image data value S (x, y).
This is the end of the image data correction process.
[0047]
According to the image reading method and the image reading apparatus 1 described in the present embodiment, the image reading apparatus 1 performs reading scanning at the pitch P with respect to the photostimulable phosphor plate 21 to which uniform input has been made. Correction data value Sc ₀ (X, y) is acquired, and the correction data value Sc according to the maximum number of pixels n that may be shifted when the photostimulable phosphor plate 21 accommodated in the cassette 2 is set in the image reading apparatus 1. ₀ A frequency component exceeding (1 / (4n · P)) is removed from (x, y), and this is set as the final correction data value Sc (x, y), which is recorded on the photostimulable phosphor plate 21. The pre-correction image data value S obtained by scanning the image of the subject with the image reading apparatus 1 ₀ (X, y) is corrected by the final correction data value Sc (x, y), and the corrected image data value S (x, y) is calculated. ₀ It is possible to obtain an image data value from which two-dimensional noise ranging from low frequency to high frequency is effectively removed without amplifying high frequency noise in (x, y).
[0048]
Further, the radiation dose used when applying uniform input to the cassette 2 including the photostimulable phosphor plate 21 is 2.5 × 10 6. ^-6 Since it is equal to or greater than [C / kg], the correction data value Sc ₀ (X, y) can be attributed solely to the structural mottle of the photostimulable phosphor plate 21.
Therefore, noise due to the structural mottle of the photostimulable phosphor plate 21 of the cassette 2 can be effectively removed.
[0049]
The image reading method and the image reading apparatus of the present invention are not limited to the above-described embodiments, and various improvements and design changes may be made without departing from the spirit of the present invention.
For example, in the above embodiment, the correction data value Sc ₀ As a method of calculating the final correction data value Sc (x, y) from (x, y), the case of using the Fourier transform and the inverse Fourier transform has been described. However, using the sinc function or the moving average method, the space (x, y y) correction data value Sc in the region ₀ If the final correction data value Sc (x, y) is obtained from (x, y), the Fourier transform and the inverse Fourier transform are not performed. Therefore, the image reading medium can be loaded into the image reading apparatus with a relatively small amount of calculation. The correction data value Sc according to the maximum number of pixels that may be shifted when setting ₀ By removing frequency components exceeding (1 / (4n · P)) from (x, y) and making this as the final correction data value Sc (x, y), it is possible to reliably carry out with a small amount of calculation. .
When the image data value is a value obtained by Log conversion, the corrected data value Sc from which the high frequency component is removed. ₀ (X, y) is normalized by the difference with respect to the maximum value to obtain the final correction data value Sc (x, y), and the pre-correction image data value S ₀ The corrected image data value S (x, y) may be calculated by adding the final correction data value Sc (x, y) to (x, y).
In the above embodiment, the radiation image reading apparatus has been described. However, the present invention can also be applied to a general image scanner or the like.
In addition, it is needless to say that specific detailed structures can be changed as appropriate.
[0050]
【The invention's effect】
According to the first or seventh aspect of the present invention, the correction data value Sc is obtained by reading and scanning the image recording medium with uniform input at the pitch P by the image reading device. ₀ (X, y) is acquired, and the correction data value Sc according to the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. ₀ A frequency component exceeding (1 / (4n · P)) is removed from (x, y) and set as a final correction data value Sc (x, y), and the image of the subject recorded on the image recording medium is set. On the other hand, the pre-correction image data value S obtained by performing scanning scanning with the image reading apparatus. ₀ (X, y) is corrected by the final correction data value Sc (x, y), and the corrected image data value S (x, y) is calculated. ₀ It is possible to obtain an image data value from which two-dimensional noise ranging from low frequency to high frequency is effectively removed without amplifying high frequency noise in (x, y).
[0051]
According to the invention described in claim 2 or 8, the same effect as that of the invention described in claim 1 can be obtained, and the correction data value Sc ₀ The set of (x, y) is Fourier transformed to obtain the second correction data value Sc ₂ (F _x , F _y ), And a frequency component exceeding (1 / (4n · P)) is removed from the third correction data value Sc _Three (F _x , F _y ), And this is subjected to inverse Fourier transform to obtain final correction data Sc (x, y). Therefore, the correction data value depends on the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. Sc ₀ It is possible to reliably remove the frequency component exceeding (1 / (4n · P)) from (x, y) and set it as the final correction data value Sc (x, y).
[0052]
According to the third or ninth aspect of the invention, the same effect as that of the first aspect of the invention can be obtained, and the correction data value Sc ₀ Since the frequency component exceeding (1 / (4n · P)) from (x, y) is removed by calculation in the (x, y) space region, the final correction data value Sc (x, y) is calculated. Since the Fourier transform and the inverse Fourier transform are not performed as in the second aspect of the invention, the maximum amount of pixels that may be shifted when the image recording medium is set in the image reading apparatus with a relatively small amount of calculation. Correction data value Sc ₀ It is possible to remove frequency components exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y).
[0053]
According to the invention described in claim 4 or 10, the same effect as that of the invention described in claim 3 can be obtained, and the correction data value Sc can be obtained. ₀ When calculating the final correction data value Sc (x, y) from (x, y), the spatial frequency (f _x , F _y ) Since the sinc function obtained as a result of inverse Fourier transform of the gate function for removing high frequency components in the region is used, the correction data value is determined according to the maximum number of pixels that may be shifted when the image recording medium is set in the image reading apparatus. Sc ₀ It is possible to reliably remove a frequency component exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y) with a small amount of calculation.
[0054]
According to the invention described in claim 5 or 11, the same effect as that of the invention described in claim 3 can be obtained, and the correction data value Sc can be obtained. ₀ When calculating the final correction data value Sc (x, y) from (x, y), the calculation by the moving average method is used, so that the maximum pixel that may be shifted when the image recording medium is set in the image reading apparatus. Depending on the number, the correction data value Sc ₀ It is possible to reliably remove a frequency component exceeding (1 / (4n · P)) from (x, y) and use this as the final correction data value Sc (x, y) by a simple calculation method.
[0055]
According to the invention described in claim 6 or claim 12, the same effect as that of the invention described in claim 1 can be obtained, and when uniform input is performed on the image recording medium including the photostimulable phosphor plate. The irradiation dose of radiation used for the test is 2.5 × 10 ^-6 Since it is equal to or greater than [C / kg], the correction data value Sc ₀ (X, y) can be attributed solely to the structure of the image recording medium.
[Brief description of the drawings]
FIG. 1 is a front view schematically showing an example of an image reading apparatus according to the present invention.
FIG. 2 is a conceptual diagram showing a method of recording an image on an image recording medium used in the image reading method according to the present invention.
FIG. 3 is a front view schematically showing an example of an image reading apparatus according to the present invention.
FIG. 4 is a graph showing characteristics of a filter used for calculation of a third correction data image recording medium in the image reading method according to the present invention.
FIG. 5 is a graph showing the sensitivity unevenness correction state of the image recording medium in the image reading method according to the present invention and the conventional image reading method.
FIG. 6 is a graph showing a sensitivity unevenness correction method for an image recording medium in a conventional image reading method.
[Explanation of symbols]
1 Image reader
12 ID detector
161 Correction data calculation means, image data correction means (CPU)
164 Image data recording means (hard disk)
21 Image recording medium (stimulable phosphor plate)

Claims

An image recording medium on which an image is recorded by a predetermined recording means, and the recorded image can be erased and recorded again is recorded for each reading coordinate (x, y) of the sampling pitch P by the image reading apparatus. In an image reading method that scans and converts the image into an image data value corresponding to each reading coordinate (x, y),
The image reading medium that is uniformly input by the recording means is scanned by the image reading device to obtain image data values corresponding to the respective reading coordinates (x, y), and this is corrected. Let the data value Sc ₀ (x, y) be
The maximum absolute value of the positional deviation amount when setting the image recording medium in the image reading device is d, and n is a natural number.
(N-1) P <d ≦ n · P
For n, the spatial frequency f of the correction data value Sc ₀ (x, y) is
f> 1 / (4n · P)
This was the final correction data value Sc (x, y) to remove frequency components corresponding to the respective spatial frequency f (f _x or f _y) is,
An uncorrected image data value S ₀ (x, y) obtained by performing scanning scanning with the image reading apparatus on the image recording medium on which an image is recorded by the recording means is used as the final corrected data value Sc ( A method of reading an image, wherein the corrected image data value S (x, y) is calculated by correcting with x, y).

The image reading method according to claim 1,
The correction data value Sc ₀ (x, y) a set of Fourier transform, spatial frequency (f _x, f _y) a second correction data value Sc ₂ (f _x, f _y) is a function of the,
F (f _x , f _y ) = 1 (when 0 <f ≦ 1 / (4n · P))
F (f _x , f _y ) = 0 (when f> 1 / (4n · P))
In a function F (f _x, f _y) by using the third correction data value _{_{_{Sc 3 (f x, f y}}} )
_{_{_{Sc 3 (f x, f y}}} ) = F (f x, f y) · Sc 2 (f x, f y),
And then, the third correction data value _{_{_{Sc 3 (f x, f y}}} ) by inverse Fourier transform, the final correction data Sc (x, y) to be,
An image reading method characterized by the above.

The image reading method according to claim 1,
Of the corrected data value Sc ₀ (x, y), the spatial frequency f is
f> 1 / (4n · P)
The corresponding frequency components in the spatial frequency f (f _x or f _y) is, characterized by calculating the (x, y) is removed in operation in the spatial domain, the final correction data value Sc (x, y) An image reading method.

The image reading method according to claim 3.
An image reading method using a sinc function when calculating the final correction data value Sc (x, y) from the correction data value Sc ₀ (x, y).

The image reading method according to claim 3.
An image reading method, wherein a calculation using a moving average method is used when calculating the final correction data value Sc (x, y) from the correction data value Sc ₀ (x, y).

In the image reading method according to any one of claims 1 to 5,
The recording means is radiation, and the image recording medium includes a stimulable phosphor plate for accumulating and recording radiation,
An image reading method characterized by applying the uniform input to the image recording medium with radiation having an irradiation dose of 2.5 × 10 ⁻⁶ [C / kg] or more.

An image recording medium on which an image is recorded by a predetermined recording means and the recorded image can be erased and recorded again is read and scanned at every reading coordinate (x, y) of the sampling pitch P, In the image reading apparatus for converting the image into an image data value corresponding to each reading coordinate (x, y),
An ID detection unit capable of detecting an ID of the image recording medium;
Image data recording means for scanning the image recording medium to obtain image data values corresponding to the respective reading coordinates, and recording the image data values as image data;
The image recording medium on which uniform input has been made is read and scanned by the image data recording means to obtain image data values corresponding to the respective reading coordinates (x, y), and this is corrected data. Record as the value Sc ₀ (x, y),
Assuming that the maximum absolute value of the positional deviation amount when setting the image recording medium in the reading device is d, n is a natural number,
(N-1) P <d ≦ n · P
For n, the spatial frequency f of the correction data value Sc ₀ (x, y) is
f> 1 / (4n · P)
Removing the frequency components corresponding to the respective spatial frequency f (f _x or f _y) is, which as a final correction data value Sc (x, y), the correction data to be recorded corresponding to the ID of the image recording medium A calculation means;
Image data correction means for correcting the uncorrected image data value S ₀ (x, y) with the final corrected data value Sc (x, y) and calculating the corrected image data value S (x, y); An image reading apparatus comprising:

The image reading apparatus according to claim 7.
The correction data calculating means, the correction data value Sc ₀ (x, y) a set of Fourier transform, the spatial frequency (f _x, f _y) a second correction data value which is a function of Sc2 (f _x, f _y )
F (f _x , f _y ) = 1 (when 0 <f ≦ 1 / (4n · P))
F (f _x , f _y ) = 0 (when f> 1 / (4n · P))
In a function F (f _x, f _y) by using the third correction data value _{_{_{Sc 3 (f x, f y}}} )
_{_{_{Sc 3 (f x, f y}}} ) = F (f x, f y) · Sc 2 (f x, f y),
Is calculated, and the third correction data value Sc3 (f _x, f _y) by inverse Fourier transform, to calculate the spatial (x, y) is a function of the final correction data Sc (x, y) A characteristic image reading apparatus.

The image reading apparatus according to claim 7.
The correction data calculating means determines that the spatial frequency f of the correction data value Sc ₀ (x, y) is
f> 1 / (4n · P)
The corresponding frequency components in the spatial frequency f (f _x or f _y) is, characterized by calculating the (x, y) is removed in operation in the spatial domain, the final correction data value Sc (x, y) An image reading apparatus.

The image reading apparatus according to claim 9.
An image reading apparatus, wherein a sinc function is used in calculating the final correction data value Sc (x, y) from the correction data value Sc ₀ (x, y).

The image reading apparatus according to claim 9.
An image reading apparatus, wherein a calculation based on a moving average method is used when calculating the final correction data value Sc (x, y) from the correction data value Sc ₀ (x, y).

The image reading apparatus according to any one of claims 7 to 11,
The recording means is radiation, and the image recording medium includes a stimulable phosphor plate for accumulating and recording radiation,
An image reading apparatus characterized by applying the uniform input to the image recording medium with radiation having an irradiation dose of 2.5 × 10 ⁻⁶ [C / kg] or more.