JP2004007458A

JP2004007458A - Method and device for correcting sharpness in image reader

Info

Publication number: JP2004007458A
Application number: JP2003060120A
Authority: JP
Inventors: Hiroaki Yasuda; 安田　裕昭
Original assignee: Fuji Photo Film Co Ltd
Current assignee: Fujifilm Holdings Corp
Priority date: 2002-03-22
Filing date: 2003-03-06
Publication date: 2004-01-08
Anticipated expiration: 2023-03-06
Also published as: JP4031720B2

Abstract

<P>PROBLEM TO BE SOLVED: To obtain fine image information having approximately the same sharpness in a main scanning direction and a sub-scanning direction. <P>SOLUTION: A shading correction part 82 corrects the shading of radiation image information obtained by reading out a stimulable phosphor sheet IP by line sensors 48a-48c and then a sharpness correction part 86 prepares radiation image information to which sharpness correction processing is applied independently in the main scanning direction and the sub-scanning direction by using a spatial filter selected from a spacial filter storage part 96 by a spatial filter selection part 98 in accordance with a recording condition or a reading condition. <P>COPYRIGHT: (C)2004,JPO

Description

【０００１】
【発明の属する技術分野】
本発明は、主走査方向に複数の光電変換部が配列されてなるセンサを、記録媒体に対して前記主走査方向と略直交する副走査方向に相対的に移動させ、前記記録媒体に記録された画像情報を前記センサにより２次元的に読み取る画像読取装置における鮮鋭度補正方法および装置に関する。
【０００２】
【従来の技術】
例えば、放射線（Ｘ線、α線、β線、γ線、電子線、紫外線等）を照射すると、この放射線エネルギの一部が蓄積され、その後、可視光等の励起光を照射すると蓄積されたエネルギ強度に応じて輝尽発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、人体等の被写体の放射線画像情報をシート状の蓄積性蛍光体層を備えてなる蓄積性蛍光体シートに一旦記録した後、この蓄積性蛍光体シートにレーザ光等の励起光を照射して得られる輝尽発光光を光電変換器を用いて電気信号に変換し、得られた画像信号を写真感光材料等の記録媒体、あるいは、ＣＲＴ（Ｃａｔｈｏｄｅ　Ｒａｙ　Ｔｕｂｅ）等の表示装置に被写体の放射線可視画像として出力させる放射線画像情報記録読取システムが知られている（特許文献１参照）。
【０００３】
この場合、副走査搬送される蓄積性蛍光体シートに対して励起光を主走査方向に移動させながら照射し、得られる輝尽発光光を集光ガイドを介してフォトマルチプライア等からなる単一の光電変換器に集光させ、励起光の走査位置に同期させて読み取る、いわゆる、ポイントスキャン方式のシステムにおいては、蓄積性蛍光体シートから光電変換器に至る光学系の読取条件が主走査方向および副走査方向で略同じである。従って、輝尽発光光が蓄積性蛍光体シートから主走査方向および副走査方向に等方的に出力されるものとすると、光電変換して得られた画像信号の主走査方向および副走査方向に対する空間周波数特性が略等しく、鮮鋭度も主走査方向および副走査方向で略同じになる。
【０００４】
これに対して、主走査方向に複数の光電変換部を配列してなるセンサを副走査方向に移動させながら蓄積性蛍光体シートの放射線画像情報を読み取る、いわゆる、ラインスキャン方式のシステムでは、主走査方向と副走査方向とで読取条件が異なるため、光電変換して得られた画像信号の主走査方向および副走査方向に対する空間周波数特性は等しくならない。この場合、主走査方向に対する画像信号の空間周波数特性は、センサに輝尽発光光を導く集光光学系の主走査方向の結像特性と、センサを構成する各光電変換部の主走査方向の間隔とに依存する。また、副走査方向に対する画像信号の空間周波数特性は、前記集光光学系の副走査方向の結像特性と、前記光電変換部の副走査方向の読取間隔とに依存する。
【０００５】
例えば、集光光学系が主走査方向に多数のレンズを配列してなるセルフォックレンズ等のレンズアレイからなる場合、各光電変換部により得られた画像信号のレンズアレイに対する主走査方向の位置関係と、各光電変換部により得られた画像信号のレンズアレイに対する副走査方向の位置関係とが異なるため、画像信号の空間周波数が主走査方向と副走査方向とで一致しなくなってしまう。従って、このような画像信号から得られる画像は、主走査方向と副走査方向とで鮮鋭度に差異のある不適切なものとなってしまう。
【０００６】
【特許文献１】
特開昭５５−１２４２９号公報
【０００７】
【発明が解決しようとする課題】
本発明は、前記の不具合を解消するためになされたもので、主走査方向および副走査方向に対して略同一の鮮鋭度を有する良好な画像情報を得ることのできる画像読取装置における鮮鋭度補正方法および装置を提供することを目的とする。
【０００８】
【課題を解決するための手段】
前記の課題を解決するために、本発明では、主走査方向に複数の光電変換部が配列されてなるセンサを、記録媒体に対して前記主走査方向と略直交する副走査方向に相対的に移動させ、前記記録媒体に記録された画像情報を前記光電変換部により２次元的に読み取る画像読取装置における鮮鋭度補正方法および装置において、読み取って得られた前記画像情報の前記主走査方向に対する鮮鋭度と、前記画像情報の前記副走査方向に対する鮮鋭度とが略同一となるよう、前記主走査方向および前記副走査方向に対して夫々所望の鮮鋭度補正処理を行うことを特徴とする。
【０００９】
この場合、鮮鋭度補正処理は、例えば、空間周波数特性に影響する記録媒体に対する画像情報の記録条件または読取条件に従い、主走査方向および副走査方向に対して夫々設定した１次元の空間フィルタを選択し、夫々に対して１次元フィルタリング処理を行うことにより、等方的な鮮鋭度からなる画像情報を得ることができる。また、主走査方向に対する空間フィルタと副走査方向に対する空間フィルタとを合成してなる２次元の空間フィルタを選択し、あるいは、予め設定されている２次元のフィルタを選択し、画像情報に対する２次元フィルタリング処理を行うようにしてもよい。
【００１０】
【発明の実施の形態】
図１は、本発明の画像読取装置における鮮鋭度補正方法および装置が適用される放射線画像情報記録読取装置１０の概略構成を示す。
【００１１】
放射線画像情報記録読取装置１０は、人体等の被写体１８の放射線画像を蓄積性蛍光体シートＩＰ（記録媒体）に一旦記録する放射線記録部１２と、前記蓄積性蛍光体シートＩＰに励起光Ｌを照射し、放射線の強度に応じて発光する輝尽発光光Ｒを光電的に読み取る読取部１４と、蓄積性蛍光体シートＩＰに残存している放射線エネルギを放出させる消去光源部１６とを備える。
【００１２】
放射線記録部１２は、放射線Ｓを発生させる放射線源２０と、被写体１８を所定位置に保持する放射線透過性の撮影台２２とを備える。撮影台２２の後方には、フォトタイマ２３と、散乱線除去用のグリッド２４とが配置され、前記グリッド２４のさらに後方には、蓄積性蛍光体シートＩＰが所定の撮影位置に固定される。
【００１３】
読取部１４は、鉛直方向（矢印Ｙ方向）に延在するガイドレール６６に沿って移動可能に構成される。読取部１４を構成する筐体４２は、ガイドレール６６に係合しており、上下両端部に配設される駆動プーリ５８および従動プーリ６０に架け渡されたベルト６４によって副走査方向（矢印Ｙ方向）に移動する。
【００１４】
読取部１４を構成する筐体４２の上部には、励起光Ｌを出力するレーザダイオードアレイ３６が配設される。レーザダイオードアレイ３６から出力された励起光Ｌは、筐体４２に固定されたミラー３７によって反射され、蓄積性蛍光体シートＩＰの主走査方向（矢印Ｙ方向と直交する方向）に照射される。
【００１５】
また、読取部１４を構成する筐体４２の下部には、励起光Ｌによって蓄積性蛍光体シートＩＰから取り出された輝尽発光光Ｒを集光レンズアレイ４４およびプリズム４６を介して受光し、電気信号に変換するＣＣＤ等のラインセンサ４８ａ〜４８ｃが配設される。集光レンズアレイ４４は、例えば、主走査方向に多数のレンズを配列してなるセルフォックレンズによって構成することができる。プリズム４６は、筐体４２の側部に配設されるラインセンサ４８ａ、４８ｃと、筐体４２の下部に配設されるラインセンサ４８ｂとに輝尽発光光Ｒを振り分ける。千鳥状に配設された各ラインセンサ４８ａ〜４８ｃは、夫々多数の光電変換部を主走査方向に配列して構成される。
【００１６】
なお、ラインセンサ４８ａ〜４８ｃは、主走査方向（矢印Ｘ方向）に複数の光電変換部を配列したものであれば、副走査方向（矢印Ｙ方向）に複数ラインを有するセンサや、主走査方向に長尺なエリアセンサを用いることもできる。また、ラインセンサ４８ａ〜４８ｃは、必ずしも３本で構成する必要はない。ラインセンサ４８ａ〜４８ｃには、蓄積性蛍光体シートＩＰから発せられる輝尽発光光Ｒを透過する一方、励起光Ｌを含む外光を遮断する図示しないフィルタが装着される。
【００１７】
消去光源部１６は、読取部１４を挟んで蓄積性蛍光体シートＩＰに対向して配設される。消去光源部１６は、蓄積性蛍光体シートＩＰに消去光Ｑを照射し、蓄積性蛍光体シートＩＰに残存する放射線エネルギを放出させる複数の蛍光灯６８を備える。
【００１８】
図２は、放射線画像情報記録読取装置１０の読取部１４における処理回路ブロックを示す。
【００１９】
読取部１４の処理回路は、ラインセンサ４８ａ〜４８ｃによって光電変換して得られた電気信号を増幅する増幅器７０と、増幅された電気信号をデジタル信号としての画像信号に変換するＡ／Ｄ（Ａｎａｌｏｇ／Ｄｉｇｉｔａｌ）変換器７２と、ラインセンサ４８ａ〜４８ｃのオフセット補正を行うオフセット補正部７４と、ラインセンサ４８ａ〜４８ｃによって光電変換された電気信号の非線形性を補正して線形信号に変換する線形変換部７６と、ラインセンサ４８ａ〜４８ｃの接続部分における信号調整を行うつなぎ処理部７８と、画像信号を対数変換する対数変換部８０と、対数変換された画像信号のシェーディング補正を行うシェーディング補正部８２と、シェーディング補正された画像信号の主走査方向および副走査方向の鮮鋭度を補正する鮮鋭度補正部８６とを備える。
【００２０】
Ａ／Ｄ変換器７２には、オフセット補正部７４に供給するオフセット補正データを演算するオフセット補正データ演算部８８が接続される。オフセット補正データ演算部８８によって演算されたオフセット補正データは、オフセット補正データ記憶部９０に記憶される。この場合、オフセット補正データとは、光を受光しない状態での各ラインセンサ４８ａ〜４８ｃの光電変換部からの出力（暗時出力）を補正する補正データである。
【００２１】
対数変換部８０には、シェーディング補正部８２に供給するシェーディング補正データを演算するシェーディング補正データ演算部９２が接続される。シェーディング補正データ演算部９２によって演算されたシェーディング補正データは、シェーディング補正データ記憶部９４に記憶される。この場合、シェーディング補正データとは、例えば、主走査方向に配列されたレーザダイオードアレイ３６、集光レンズアレイ４４、ラインセンサ４８ａ〜４８ｃを構成する光電変換部の１次元方向（主走査方向）の特性ばらつきを補正する高周波成分を含む補正データや、放射線源２０から出力される放射線Ｓの照射むら、蓄積性蛍光体シートＩＰの蓄積性蛍光体の塗布むら、グリッド２４の散乱線除去特性のむら等を補正する低周波成分を主体とした補正データである。
【００２２】
鮮鋭度補正部８６には、鮮鋭度補正処理を行うための複数種類の空間フィルタを構成するフィルタリング係数を記憶する空間フィルタ記憶部９６が空間フィルタ選択部９８を介して接続される。この空間フィルタ記憶部９６に記憶された空間フィルタは、蓄積性蛍光体シートＩＰに対する放射線画像情報の記録条件、あるいは、蓄積性蛍光体シートＩＰに記録された放射線画像情報の読取条件に従い、空間フィルタ選択部９８によって選択される。この場合、記録条件とは、例えば、放射線画像情報記録読取装置１０による被写体１８に対する撮影メニューによって決定される撮影部位、放射線量、放射線画像情報の記録される蓄積性蛍光体シートＩＰの種類等であり、読取条件とは、例えば、蓄積性蛍光体シートＩＰに記録された放射線画像情報の読取モードによって決定される読取速度等である。
【００２３】
本実施形態の放射線画像情報記録読取装置１０は、基本的には以上のように構成されるものであり、次に、その動作並びに作用効果について説明する。
【００２４】
先ず、図３に示すフローチャートに基づき、放射線画像情報記録読取装置１０におけるオフセット補正データおよびシェーディング補正データの作成方法について説明する。なお、以下の説明では、図４に示すように、ラインセンサ４８ａ〜４８ｃの各光電変換部により蓄積性蛍光体シートＩＰを主走査方向（矢印Ｘ方向）にｎ画素に分割するとともに、副走査方向（矢印Ｙ方向）にｍ画素に分割し、矢印Ｘ方向および矢印Ｙ方向を座標軸とした各画素の座標を（ｘ，ｙ）で表すものとする。
【００２５】
読取部１４のレーザダイオードアレイ３６および消去光源部１６の蛍光灯６８を消灯状態とし、この状態でのラインセンサ４８ａ〜４８ｃからの電気信号を増幅器７０およびＡ／Ｄ変換器７２を介して検出する。オフセット補正データ演算部８８は、検出された電気信号を各ラインセンサ４８ａ〜４８ｃにおける各光電変換部のオフセット補正データＯ（ｘ）とし、オフセット補正データ記憶部９０に記憶させる（ステップＳ１）。
【００２６】
次に、被写体１８を介在させることなく、放射線源２０から出力される放射線Ｓをフォトタイマ２３およびグリッド２４を介して蓄積性蛍光体シートＩＰに一様に照射する。読取部１４は、このようにして一様露光された蓄積性蛍光体シートＩＰの読み取りを行う（ステップＳ２）。この場合、一様露光された蓄積性蛍光体シートＩＰに記録された画像情報には、放射線Ｓの照射むらや蓄積性蛍光体シートＩＰにおける蓄積性蛍光体の塗布むらに基づく変化の小さい低周波成分からなるむらが含まれている。
【００２７】
レーザダイオードアレイ３６から出力された励起光Ｌは、ミラー３７によって反射された後、蓄積性蛍光体シートＩＰの主走査方向に照射される。蓄積性蛍光体シートＩＰに蓄積された放射線エネルギは、輝尽発光光Ｒとして出力され、集光レンズアレイ４４およびプリズム４６を介してラインセンサ４８ａ〜４８ｃに導入される。ラインセンサ４８ａ〜４８ｃは、輝尽発光光Ｒを電気信号に変換する。この電気信号は、増幅器７０によって増幅された後、Ａ／Ｄ変換器７２によってデジタル信号としての補正データ取得画像データＤｏｒｇ（ｘ，ｙ）に変換される。なお、補正データ取得画像データＤｏｒｇ（ｘ，ｙ）には、レーザダイオードアレイ３６、集光レンズアレイ４４、ラインセンサ４８ａ〜４８ｃを構成する光電変換部の個々の特性ばらつきに起因するランダムな高周波成分からなるむらが含まれている。
【００２８】
そこで、オフセット補正部７４は、ステップＳ１で作成されたオフセット補正データＯ（ｘ）をオフセット補正データ記憶部９０から読み込み、
Ｄｏｆｆ（ｘ，ｙ）＝Ｄｏｒｇ（ｘ，ｙ）−Ｏ（ｘ）
として補正データ取得画像データＤｏｒｇ（ｘ，ｙ）のオフセット補正を行い、オフセット補正画像データＤｏｆｆ（ｘ，ｙ）を作成する（ステップＳ３）。
【００２９】
オフセット補正画像データＤｏｆｆ（ｘ，ｙ）は、次に、線形変換部７６において線形変換された後（ステップＳ４）、つなぎ処理部７８において各ラインセンサ４８ａ〜４８ｃ間のつなぎ処理が行われる（ステップＳ５）。つなぎ処理されたオフセット補正画像データＤｏｆｆ（ｘ，ｙ）は、対数変換部８０において、加算演算処理が可能な対数変換データに変換される（ステップＳ６）。
【００３０】
シェーディング補正データ演算部９２は、対数変換されたオフセット補正画像データＤｌｏｇ（ｘ，ｙ）＝Ｌｏｇ（Ｄｏｆｆ（ｘ，ｙ））を各光電変換部毎に矢印Ｙ方向に加算平均し、平均オフセット補正画像データＤｓ（ｘ）を、
Ｄｓ（ｘ）＝１／ｍ・ΣＤｌｏｇ（ｘ，ｙ）
として求める（ステップＳ７）。上式のΣは、オフセット補正画像データＤｌｏｇ（ｘ，ｙ）を矢印Ｙ方向にｍ画素分加算することを意味するものとする。
【００３１】
次に、平均オフセット補正画像データＤｓ（ｘ）を矢印Ｘ方向に加算平均し、平均値Ｄｓ＿ａｖｅを、
Ｄｓ＿ａｖｅ＝１／ｎ・ΣＤｓ（ｘ）
として求める（ステップＳ８）。上式のΣは、平均オフセット補正画像データＤｓ（ｘ）を矢印Ｘ方向にｎ画素分加算することを意味するものとする。
【００３２】
最後に、平均値Ｄｓ＿ａｖｅと平均オフセット補正画像データＤｓ（ｘ）とから、１次元補正データＳ１（ｘ）を、
Ｓ１（ｘ）＝Ｄｓ＿ａｖｅ−Ｄｓ（ｘ）
として求める（ステップＳ９）。このようにして求められた１次元補正データＳ１（ｘ）を用いることにより、レーザダイオードアレイ３６、集光レンズアレイ４４、ラインセンサ４８ａ〜４８ｃを構成する各光電変換部の個々の特性ばらつきに起因するランダムな高周波成分からなるシェーディングを補正することができる。算出された１次元補正データＳ１（ｘ）は、シェーディング補正データ記憶部９４に記憶される。
【００３３】
次に、１次元補正データＳ１（ｘ）を用いて対数変換されたオフセット補正画像データＤｌｏｇ（ｘ，ｙ）の１次元補正を行い、１次元補正画像データＤ１（ｘ，ｙ）を、
Ｄ１（ｘ，ｙ）＝Ｄｌｏｇ（ｘ，ｙ）＋Ｓ１（ｘ）
として求める（ステップＳ１０）。
【００３４】
さらに、１次元補正画像データＤ１（ｘ，ｙ）に対してアンシャープマスキング処理を施し（ステップＳ１１）、低周波補正画像データＤｌｏｗ（ｘ，ｙ）を求める。この低周波補正画像データＤｌｏｗ（ｘ，ｙ）は、例えば、注目する１次元補正画像データＤ１（ｘ，ｙ）の周囲の複数の１次元補正画像データＤ１（ｘ，ｙ）を加算平均することで得ることができる。すなわち、矢印Ｘ方向の処理範囲の画素数をｊ、矢印Ｙ方向の処理範囲の画素数をｋとして、低周波補正画像データＤｌｏｗ（ｘ，ｙ）は、
Ｄｌｏｗ（ｘ，ｙ）＝１／（ｊ・ｋ）・ΣΣＤ１（ｘ，ｙ）
として求めることができる。上式のΣは、夫々１次元補正画像データＤ１（ｘ，ｙ）を矢印Ｘ方向にｊ画素分、矢印Ｙ方向にｋ画素分加算することを意味するものとする。このようにして１次元補正画像データＤ１（ｘ，ｙ）を平滑化することにより、高周波成分からなるノイズの影響を小さくすることができる。
【００３５】
次いで、低周波補正画像データＤｌｏｗ（ｘ，ｙ）を矢印Ｘ方向および矢印Ｙ方向に加算平均し、平均値Ｄａｖｅを、
Ｄａｖｅ＝１／（ｎ・ｍ）・ΣΣＤｌｏｗ（ｘ，ｙ）
として求める（ステップＳ１２）。上式のΣΣは、夫々低周波補正画像データＤｌｏｗ（ｘ，ｙ）を矢印Ｘ方向にｎ画素分、矢印Ｙ方向にｍ画素分加算することを意味するものとする。
【００３６】
平均値Ｄａｖｅと低周波補正画像データＤｌｏｗ（ｘ，ｙ）とから、平滑化された２次元補正データＳ２（ｘ，ｙ）を、
Ｓ２（ｘ，ｙ）＝Ｄａｖｅ−Ｄｌｏｗ（ｘ，ｙ）
として求める（ステップＳ１３）。このようにして求められた２次元補正データＳ２（ｘ，ｙ）を用いることにより、放射線Ｓの照射むらや蓄積性蛍光体シートＩＰにおける蓄積性蛍光体の塗布むらに基づく変化の小さい低周波成分からなるシェーディングを補正することができる。
【００３７】
平滑化された２次元補正データＳ２（ｘ，ｙ）は、低周波成分からなるむらを補正するものであって、高周波成分を含んでいないため、隣接する２次元補正データＳ２（ｘ，ｙ）間で大きく変動することはない。そこで、２次元補正データＳ２（ｘ，ｙ）を間引き処理することにより、データ数の少ない２次元補正データＳ２（ｘ，ｙ）とし（ステップＳ１４）、１次元補正データＳ１（ｘ）とともにシェーディング補正データ記憶部９４に記憶する。
【００３８】
次に、図５に示すフローチャートに基づき、放射線画像情報記録読取装置１０における放射線画像情報の読取処理について説明する。
【００３９】
撮影部位、撮影管電圧、撮影線量等を所定の撮影メニューに従って設定し、被写体１８を放射線画像情報記録読取装置１０における放射線記録部１２の所定位置に位置決めした後、放射線源２０を駆動して放射線Ｓを被写体１８に照射する。被写体１８を透過した放射線Ｓは、フォトタイマ２３およびグリッド２４を介して蓄積性蛍光体シートＩＰに到達し、被写体１８の放射線画像情報が記録される（ステップＳ２１）。
【００４０】
次に、蓄積性蛍光体シートＩＰに対する放射線画像情報の読み取りに先立ち、上述したステップＳ１の場合と同様にして、オフセット補正データＯ（ｘ）の作成を行う（ステップＳ２２）。すなわち、オフセット補正データＯ（ｘ）は、放射線画像情報記録読取装置１０の温度等に依存して変動するものであるため、一定時間周期あるいは各読取動作の直前に行うことが望ましい。作成されたオフセット補正データＯ（ｘ）は、オフセット補正データ記憶部９０に記憶される。
【００４１】
次に、低速走査による高画質モード、高速走査による高速読取モード等の読取モードを設定し、その読取モードに従って読取部１４が駆動され、蓄積性蛍光体シートＩＰに記録された放射線画像情報の読み取りが行われる（ステップＳ２３）。蓄積性蛍光体シートＩＰから読み取られた放射線画像情報に係る画像データＤｏｒｇ（ｘ，ｙ）は、オフセット補正部７４において、オフセット補正データ記憶部９０から取得したオフセット補正データＯ（ｘ）を用いてオフセット補正処理が行われ、オフセット補正画像データＤｏｆｆ（ｘ，ｙ）が、
Ｄｏｆｆ（ｘ，ｙ）＝Ｄｏｒｇ（ｘ，ｙ）−Ｏ（ｘ）
として求められる（ステップＳ２４）。
【００４２】
オフセット補正画像データＤｏｆｆ（ｘ，ｙ）は、線形変換部７６で線形変換された後（ステップＳ２５）、つなぎ処理部７８でラインセンサ４８ａ〜４８ｃ間のつなぎ処理が施され（ステップＳ２６）、さらに、対数変換部８０で対数変換され（ステップＳ２７）、
Ｄｌｏｇ（ｘ，ｙ）＝Ｌｏｇ（Ｄｏｆｆ（ｘ，ｙ））
としてオフセット補正画像データＤｌｏｇ（ｘ，ｙ）が得られる。
【００４３】
対数変換されたオフセット補正画像データＤｌｏｇ（ｘ，ｙ）は、シェーディング補正部８２において、シェーディング補正データ記憶部９４に記憶されている１次元補正データＳ１（ｘ）を用いて、
Ｄ１（ｘ，ｙ）＝Ｄｌｏｇ（ｘ，ｙ）＋Ｓ１（ｘ）
として１次元補正され、１次元補正画像データＤ１（ｘ，ｙ）が得られる（ステップＳ２８）。
【００４４】
次いで、シェーディング補正データ記憶部９４に記憶されている２次元補正データＳ２（ｘ，ｙ）を用いて、１次元補正画像データＤ１（ｘ，ｙ）が、
Ｄ２（ｘ，ｙ）＝Ｄ１（ｘ，ｙ）＋Ｓ２（ｘ，ｙ）
として補正され、２次元補正された補正画像データＤ２（ｘ，ｙ）が得られる（ステップＳ２９）。なお、シェーディング補正データ記憶部９４には、データ量を少なくするため、ステップＳ１４の処理で間引かれた２次元補正データＳ２（ｘ，ｙ）が記憶されているため、補正画像データＤ２（ｘ，ｙ）の算出に当たっては、補間した２次元補正データＳ２（ｘ，ｙ）を用いて補正画像データＤ２（ｘ，ｙ）の算出を行う。
【００４５】
１次元補正および２次元補正された補正画像データＤ２（ｘ，ｙ）は、鮮鋭度補正部８６において、主走査方向に対する鮮鋭度補正処理が施され（ステップＳ３０）、次いで、副走査方向に対する鮮鋭度補正処理が施される（ステップＳ３１）。
【００４６】
ここで、例えば、被写体１８に照射される放射線量は、撮影メニューに従い放射線源２０の撮影管電圧を撮影部位によって制御することで調整される。また、被写体１８に照射される放射線量が同じであっても、蓄積性蛍光体シートＩＰに到達する放射線量は、撮影部位によって異なる。一方、蓄積性蛍光体シートＩＰの種類によって、放出される輝尽発光光Ｒの広がり特性が異なる。さらに、高速走査するか、低速走査するかといった読取モードの違いにより蓄積性蛍光体シートＩＰに照射される励起光Ｌのエネルギが異なるため、その励起光Ｌによって放出される輝尽発光光Ｒの広がり特性も異なってくる。この場合、蓄積性蛍光体シートＩＰから得られる放射線画像情報の空間周波数特性は、蓄積性蛍光体シートＩＰに照射される放射線量に応じて蓄積された放射線エネルギやその放出特性に依存する。
【００４７】
一方、放射線画像情報の空間周波数特性は、ラインセンサ４８ａ〜４８ｃに輝尽発光光Ｒを導く集光レンズアレイ４４およびプリズム４６の主走査方向および副走査方向の特性に依存するとともに、ラインセンサ４８ａ〜４８ｃを構成する各光電変換部の主走査方向の間隔および副走査方向の読取間隔に依存する。
【００４８】
従って、鮮鋭度補正前の放射線画像情報である補正画像データＤ２（ｘ，ｙ）の主走査方向および副走査方向に対する鮮鋭度およびそのバランスは、撮影メニュー、蓄積性蛍光体シートＩＰの種類、読取モード等によって異なることになる。
【００４９】
そこで、空間フィルタ選択部９８は、オペレータによって設定された撮影メニュー、蓄積性蛍光体シートＩＰの種類、読取モード等の記録条件または読取条件に従い、鮮鋭度およびそのバランスに対応した空間フィルタを空間フィルタ記憶部９６から選択し、鮮鋭度補正部８６に供給する。鮮鋭度補正部８６は、供給された空間フィルタを用いて、補正画像データＤ２（ｘ，ｙ）に対する鮮鋭度補正処理を行う。
【００５０】
すなわち、主走査方向の鮮鋭度を補正する空間フィルタを構成するフィルタリング係数をｆ（ａ）とすると、補正画像データＤ２（ｘ，ｙ）は、
Ｄｘ（ｘ，ｙ）＝１／Ｎ・Σｆ（ａ）・Ｄ２（ｘ＋ａ，ｙ）
として、主走査方向に補正された画像データＤｘ（ｘ，ｙ）が求められる（ステップＳ３０）。ａは、フィルタリング処理を行うために利用される画素範囲を表す整数であり、Σは、画素範囲ａでの加算処理を表す。また、Ｎ＝Σｆ（ａ）である。
【００５１】
図６は、フィルタリング係数ｆ（ａ）を、
ｆ（ａ）＝（−２　３　１０　３　−２）
としたときの主走査方向に対する空間周波数（サイクル／ｍｍ）とレスポンスとの関係を表したものである。このフィルタリング係数ｆ（ａ）を用いた場合、鮮鋭度補正前の補正画像データＤ２（ｘ，ｙ）の鮮鋭度を主走査方向に強める効果が得られる。
【００５２】
また、鮮鋭度補正部８６は、空間フィルタ記憶部９６から選択した副走査方向の鮮鋭度を補正する空間フィルタを構成するフィルタリング係数をｇ（ｂ）として、主走査方向の鮮鋭度が補正された画像データＤｘ（ｘ，ｙ）を、
Ｄｙ（ｘ，ｙ）＝１／Ｍ・Σｇ（ｂ）・Ｄｘ（ｘ，ｙ＋ｂ）
として補正することにより、副走査方向に補正された画像データＤｙ（ｘ，ｙ）を求める（ステップＳ３１）。ｂは、フィルタリング処理を行うために利用される画素範囲を表す整数であり、Σは、画素範囲ｂでの加算処理を表す。また、Ｍ＝Σｇ（ｂ）である。
【００５３】
図７は、フィルタリング係数ｇ（ｂ）を、
【００５４】
【数１】

【００５５】
としたときの主走査方向に対する空間周波数（サイクル／ｍｍ）とレスポンスとの関係を表したものである。このフィルタリング係数ｇ（ｂ）を用いた場合、画像データＤｘ（ｘ，ｙ）の鮮鋭度を副走査方向に弱める効果が得られる。
【００５６】
以上のようにして、適切なフィルタリング係数ｆ（ａ）、ｇ（ｂ）を有する空間フィルタを用いることにより、主走査方向および副走査方向に対して略同一の鮮鋭度補正処理が施された画像データＤｙ（ｘ，ｙ）を得ることができる。
【００５７】
次いで、画像データＤｙ（ｘ，ｙ）の規格化処理を行った後（ステップＳ３２）、階調補正、カラーコレクション等の画像処理が施される（ステップＳ３３）。この場合、主走査方向および副走査方向の鮮鋭度の差を考慮することなく、所望の画像処理を行うことができる。このようにして画像処理された放射線画像情報は、外部の診断装置等に提供される。なお、ステップＳ３３での画像処理は、放射線画像情報記録読取装置１０で行うのではなく、外部の画像処理装置で行うようにしてもよい。
【００５８】
一方、蓄積性蛍光体シートＩＰからの放射線画像情報の読取処理が終了すると、消去光源部１６が駆動され、蛍光灯６８から出力される消去光Ｑが蓄積性蛍光体シートＩＰの全面に照射されることにより、残存する放射線エネルギが除去され、次の撮影に供される。
【００５９】
ここで、上述した実施形態では、主走査方向に対して鮮鋭度を強める１次元フィルタリング処理を行うとともに、副走査方向に対して鮮鋭度を弱める１次元フィルタリング処理を行っているが、主走査方向または副走査方向の一方の鮮鋭度を基準として固定しておき、他方の鮮鋭度のみを一方の鮮鋭度に合わせるように調整するようにしてもよい。
【００６０】
また、主走査方向のフィルタリング係数ｆ（ａ）と副走査方向のフィルタリング係数ｇ（ｂ）との内積であるフィルタリング係数ｆｇ（ａ，ｂ）を、例えば、
【００６１】
【数２】

【００６２】
として求め、このフィルタリング係数ｆｇ（ａ，ｂ）を用いて、画像データＤ（ｘ，ｙ）を、
Ｄ（ｘ，ｙ）＝１／Ｐ・ΣΣｆｇ（ａ，ｂ）・Ｄ２（ｘ＋ａ，ｙ＋ｂ）
として２次元フィルタリング処理することにより、主走査方向および副走査方向の鮮鋭度補正処理を同時に行うこともできる。ΣΣは、主走査方向および副走査方向のフィルタリング処理の画素範囲ａ、ｂでの加算処理を表す。また、Ｐ＝ΣΣｆｇ（ａ，ｂ）である。
【００６３】
なお、上述した実施形態における放射線画像情報記録読取装置１０では、ラインセンサ４８ａ〜４８ｃを蓄積性蛍光体シートＩＰに対して副走査方向に移動させて放射線画像情報を読み取るようにしているが、固定されたラインセンサ４８ａ〜４８ｃに対して蓄積性蛍光体シートＩＰを副走査方向に移動させるように構成した放射線画像情報記録読取装置に適用することもできる。
【００６４】
【発明の効果】
以上のように、本発明によれば、主走査方向に複数の光電変換部が配列されてなるセンサを用いて読み取られた画像情報から、主走査方向および副走査方向の夫々に対して所望の鮮鋭度補正処理を行うことができ、これによって、主走査方向および副走査方向に略同一の鮮鋭度を有する良好な画像情報を得ることができる。
【図面の簡単な説明】
【図１】本実施形態の放射線画像情報記録読取装置の概略構成図である。
【図２】本実施形態の放射線画像情報記録読取装置の読取部における処理回路ブロック図である。
【図３】本実施形態の放射線画像情報記録読取装置の補正データの作成手順を示すフローチャートである。
【図４】蓄積性蛍光体シートに設定した座標系の説明図である。
【図５】放射線画像情報のシェーディング補正処理および鮮鋭度補正処理のフローチャートである。
【図６】鮮鋭度を強めることのできるフィルタにおける空間周波数とレスポンスとの関係説明図である。
【図７】鮮鋭度を弱めることのできるフィルタにおける空間周波数とレスポンスとの関係説明図である。
【符号の説明】
１０…放射線画像情報記録読取装置　１２…放射線記録部
１４…読取部　　　　　　　　　　　１６…消去光源部
１８…被写体　　　　　　　　　　　２０…放射線源
４８ａ〜４８ｃ…ラインセンサ　　　７４…オフセット補正部
８２…シェーディング補正部　　　　８６…鮮鋭度補正部
８８…オフセット補正データ演算部　９０…オフセット補正データ記憶部
９２…シェーディング補正データ演算部
９４…シェーディング補正データ記憶部
９６…空間フィルタ記憶部　　　　　９８…空間フィルタ選択部
ＩＰ…蓄積性蛍光体シート[0001]
BACKGROUND OF THE INVENTION
In the present invention, a sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction is moved relative to the recording medium in the sub-scanning direction substantially orthogonal to the main scanning direction, and is recorded on the recording medium. The present invention relates to a sharpness correction method and apparatus in an image reading apparatus that reads image information two-dimensionally by the sensor.
[0002]
[Prior art]
For example, when radiation (X-rays, α-rays, β-rays, γ-rays, electron beams, ultraviolet rays, etc.) is irradiated, a part of this radiation energy is accumulated, and thereafter, it is accumulated when irradiated with excitation light such as visible light. Using a stimulable phosphor (stimulable phosphor) that exhibits stimulated emission according to energy intensity, radiation image information of a subject such as a human body is provided with a sheet-like stimulable phosphor layer. Once recorded on the body sheet, the stimulable phosphor light obtained by irradiating the stimulable phosphor sheet with excitation light such as laser light is converted into an electrical signal using a photoelectric converter, and the obtained image signal is converted into an electrical signal. 2. Description of the Related Art A radiation image information recording / reading system that outputs a recording medium such as a photographic light-sensitive material or a display device such as a CRT (Cathode Ray Tube) as a radiation visible image of a subject is known (see Patent Document 1).
[0003]
In this case, the stimulable phosphor sheet transported in the sub-scanning direction is irradiated while moving the excitation light in the main scanning direction, and the resulting photostimulated emission light is a single unit made of a photomultiplier or the like through a condensing guide. In the so-called point scan system, which is focused on the photoelectric converter and read in synchronization with the scanning position of the excitation light, the reading condition of the optical system from the stimulable phosphor sheet to the photoelectric converter is the main scanning direction. And substantially the same in the sub-scanning direction. Accordingly, when the photostimulated luminescence light is output isotropically from the stimulable phosphor sheet in the main scanning direction and the sub-scanning direction, the image signal obtained by photoelectric conversion with respect to the main scanning direction and the sub-scanning direction. The spatial frequency characteristics are substantially equal, and the sharpness is substantially the same in the main scanning direction and the sub-scanning direction.
[0004]
On the other hand, in a so-called line scan type system that reads radiation image information of a stimulable phosphor sheet while moving a sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction in the sub scanning direction, Since scanning conditions differ between the scanning direction and the sub-scanning direction, the spatial frequency characteristics of the image signal obtained by photoelectric conversion in the main scanning direction and the sub-scanning direction are not equal. In this case, the spatial frequency characteristics of the image signal with respect to the main scanning direction are the imaging characteristics in the main scanning direction of the condensing optical system that guides the stimulated emission light to the sensor and the main scanning direction of each photoelectric conversion unit constituting the sensor. Depends on the spacing. The spatial frequency characteristic of the image signal with respect to the sub-scanning direction depends on the imaging characteristic of the condensing optical system in the sub-scanning direction and the reading interval of the photoelectric conversion unit in the sub-scanning direction.
[0005]
For example, when the condensing optical system is composed of a lens array such as a SELFOC lens in which a large number of lenses are arranged in the main scanning direction, the positional relationship in the main scanning direction with respect to the lens array of the image signal obtained by each photoelectric conversion unit Since the positional relationship in the sub scanning direction with respect to the lens array of the image signal obtained by each photoelectric conversion unit is different, the spatial frequency of the image signal does not match in the main scanning direction and the sub scanning direction. Accordingly, an image obtained from such an image signal becomes inappropriate with a difference in sharpness between the main scanning direction and the sub-scanning direction.
[0006]
[Patent Document 1]
JP-A-55-12429
[0007]
[Problems to be solved by the invention]
The present invention has been made to solve the above-described problems, and sharpness correction in an image reading apparatus capable of obtaining good image information having substantially the same sharpness in the main scanning direction and the sub-scanning direction. It is an object to provide a method and apparatus.
[0008]
[Means for Solving the Problems]
In order to solve the above-described problem, in the present invention, a sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction is relatively positioned in the sub-scanning direction that is substantially orthogonal to the main scanning direction with respect to the recording medium. In the sharpness correction method and apparatus in an image reading apparatus that moves and reads image information recorded on the recording medium two-dimensionally by the photoelectric conversion unit, the image information obtained by reading is sharpened in the main scanning direction. Desired sharpness correction processing is performed in each of the main scanning direction and the sub-scanning direction so that the degree and the sharpness of the image information in the sub-scanning direction are substantially the same.
[0009]
In this case, the sharpness correction processing selects, for example, a one-dimensional spatial filter set for each of the main scanning direction and the sub-scanning direction according to the recording condition or reading condition of the image information on the recording medium that affects the spatial frequency characteristics. Then, by performing one-dimensional filtering processing on each, image information having isotropic sharpness can be obtained. In addition, a two-dimensional spatial filter obtained by combining a spatial filter for the main scanning direction and a spatial filter for the sub-scanning direction is selected, or a preset two-dimensional filter is selected and two-dimensional for image information is selected. A filtering process may be performed.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
FIG. 1 shows a schematic configuration of a radiation image information recording / reading apparatus 10 to which a sharpness correction method and apparatus in an image reading apparatus of the present invention is applied.
[0011]
The radiation image information recording / reading device 10 temporarily records a radiation image of a subject 18 such as a human body on a stimulable phosphor sheet IP (recording medium), and excitation light L on the stimulable phosphor sheet IP. A reading unit 14 that photoelectrically reads the stimulated emission light R that is irradiated and emits light according to the intensity of the radiation, and an erasing light source unit 16 that releases the radiation energy remaining in the stimulable phosphor sheet IP are provided.
[0012]
The radiation recording unit 12 includes a radiation source 20 that generates radiation S, and a radiation-transmissive imaging table 22 that holds the subject 18 at a predetermined position. A phototimer 23 and a scattered radiation removal grid 24 are arranged behind the imaging table 22, and the stimulable phosphor sheet IP is fixed at a predetermined imaging position further behind the grid 24.
[0013]
The reading unit 14 is configured to be movable along a guide rail 66 extending in the vertical direction (arrow Y direction). The casing 42 constituting the reading unit 14 is engaged with the guide rail 66, and is driven in the sub-scanning direction (arrow Y) by a belt 64 that spans between a driving pulley 58 and a driven pulley 60 that are disposed at both upper and lower ends. Direction).
[0014]
A laser diode array 36 that outputs excitation light L is disposed on the upper portion of the casing 42 that constitutes the reading unit 14. The excitation light L output from the laser diode array 36 is reflected by the mirror 37 fixed to the housing 42 and is irradiated in the main scanning direction (direction orthogonal to the arrow Y direction) of the stimulable phosphor sheet IP.
[0015]
In the lower part of the casing 42 constituting the reading unit 14, the stimulated emission light R extracted from the stimulable phosphor sheet IP by the excitation light L is received through the condenser lens array 44 and the prism 46, Line sensors 48a to 48c such as CCDs for converting into electrical signals are arranged. The condensing lens array 44 can be constituted by, for example, a Selfoc lens formed by arranging a large number of lenses in the main scanning direction. The prism 46 distributes the stimulated emission light R to the

line sensors

48 a and 48 c disposed on the side of the housing 42 and the line sensor 48 b disposed on the lower portion of the housing 42. Each of the line sensors 48a to 48c arranged in a staggered pattern is configured by arranging a large number of photoelectric conversion units in the main scanning direction.
[0016]
The line sensors 48a to 48c may be sensors having a plurality of lines in the sub-scanning direction (arrow Y direction) or the main scanning direction as long as a plurality of photoelectric conversion units are arranged in the main scanning direction (arrow X direction). Also, a long area sensor can be used. Further, the line sensors 48a to 48c are not necessarily configured by three. The line sensors 48a to 48c are equipped with filters (not shown) that transmit the stimulated emission light R emitted from the stimulable phosphor sheet IP and block external light including the excitation light L.
[0017]
The erasing light source unit 16 is disposed to face the stimulable phosphor sheet IP with the reading unit 14 interposed therebetween. The erasing light source unit 16 includes a plurality of fluorescent lamps 68 that irradiate the storable phosphor sheet IP with erasing light Q and emit radiation energy remaining in the storable phosphor sheet IP.
[0018]
FIG. 2 shows a processing circuit block in the reading unit 14 of the radiation image information recording / reading apparatus 10.
[0019]
The processing circuit of the reading unit 14 includes an amplifier 70 that amplifies an electric signal obtained by photoelectric conversion by the line sensors 48a to 48c, and an A / D (Analog) that converts the amplified electric signal into an image signal as a digital signal. / Digital) converter 72, offset correction unit 74 that performs offset correction of the line sensors 48a to 48c, and linear conversion that corrects the nonlinearity of the electrical signals photoelectrically converted by the line sensors 48a to 48c and converts them into linear signals. Unit 76, a connection processing unit 78 that performs signal adjustment at the connection part of the line sensors 48 a to 48 c, a logarithmic conversion unit 80 that performs logarithmic conversion on the image signal, and a shading correction unit 82 that performs shading correction on the logarithmically converted image signal. And sharpness of shading-corrected image signal in the main and sub-scan directions And a sharpness correcting unit 86 for correcting.
[0020]
The A / D converter 72 is connected to an offset correction data calculation unit 88 that calculates offset correction data to be supplied to the offset correction unit 74. The offset correction data calculated by the offset correction data calculation unit 88 is stored in the offset correction data storage unit 90. In this case, the offset correction data is correction data for correcting the output (dark output) from the photoelectric conversion units of the line sensors 48a to 48c in a state where light is not received.
[0021]
The logarithmic conversion unit 80 is connected to a shading correction data calculation unit 92 that calculates shading correction data to be supplied to the shading correction unit 82. The shading correction data calculated by the shading correction data calculation unit 92 is stored in the shading correction data storage unit 94. In this case, the shading correction data is, for example, the one-dimensional direction (main scanning direction) of the photoelectric conversion units constituting the laser diode array 36, the condensing lens array 44, and the line sensors 48a to 48c arranged in the main scanning direction. Correction data including high-frequency components for correcting characteristic variations, unevenness of irradiation of radiation S output from the radiation source 20, unevenness of application of the stimulable phosphor of the stimulable phosphor sheet IP, unevenness of scattered radiation removal characteristics of the grid 24, and the like Is correction data mainly composed of a low-frequency component for correcting.
[0022]
A spatial filter storage unit 96 that stores filtering coefficients constituting a plurality of types of spatial filters for performing sharpness correction processing is connected to the sharpness correction unit 86 via a spatial filter selection unit 98. The spatial filter stored in the spatial filter storage unit 96 is a spatial filter according to the recording condition of the radiation image information on the stimulable phosphor sheet IP or the reading condition of the radiation image information recorded on the stimulable phosphor sheet IP. It is selected by the selector 98. In this case, the recording conditions are, for example, an imaging region determined by an imaging menu for the subject 18 by the radiation image information recording / reading apparatus 10, a radiation dose, a type of the stimulable phosphor sheet IP on which the radiation image information is recorded, and the like. The reading condition is, for example, a reading speed determined by a reading mode of radiation image information recorded on the stimulable phosphor sheet IP.
[0023]
The radiographic image information recording / reading apparatus 10 of the present embodiment is basically configured as described above. Next, its operation and effects will be described.
[0024]
First, a method of creating offset correction data and shading correction data in the radiation image information recording / reading apparatus 10 will be described based on the flowchart shown in FIG. In the following description, as shown in FIG. 4, each photoelectric conversion unit of the line sensors 48 a to 48 c divides the stimulable phosphor sheet IP into n pixels in the main scanning direction (arrow X direction) and performs sub-scanning. The pixel is divided into m pixels in the direction (arrow Y direction), and the coordinates of each pixel with the arrow X direction and the arrow Y direction as coordinate axes are represented by (x, y).
[0025]
The laser diode array 36 of the reading unit 14 and the fluorescent lamp 68 of the erasing light source unit 16 are turned off, and electric signals from the line sensors 48a to 48c in this state are detected via the amplifier 70 and the A / D converter 72. . The offset correction data calculation unit 88 stores the detected electrical signal in the offset correction data storage unit 90 as offset correction data O (x) of each photoelectric conversion unit in each of the line sensors 48a to 48c (step S1).
[0026]
Next, the stimulable phosphor sheet IP is uniformly irradiated through the phototimer 23 and the grid 24 with the radiation S output from the radiation source 20 without the subject 18 interposed. The reading unit 14 reads the stimulable phosphor sheet IP that has been uniformly exposed in this way (step S2). In this case, the image information recorded on the stimulable phosphor sheet IP that has been uniformly exposed has a low frequency with a small change based on the irradiation unevenness of the radiation S and the application unevenness of the stimulable phosphor on the stimulable phosphor sheet IP. Contains unevenness of ingredients.
[0027]
The excitation light L output from the laser diode array 36 is reflected by the mirror 37 and then irradiated in the main scanning direction of the stimulable phosphor sheet IP. The radiation energy accumulated in the stimulable phosphor sheet IP is output as the stimulated emission light R and introduced into the line sensors 48 a to 48 c via the condenser lens array 44 and the prism 46. The line sensors 48a to 48c convert the stimulated emission light R into an electrical signal. This electric signal is amplified by the amplifier 70 and then converted into correction data acquisition image data Dorg (x, y) as a digital signal by the A / D converter 72. The correction data acquisition image data Dorg (x, y) includes a random high-frequency component caused by individual characteristic variations of the photoelectric conversion units constituting the laser diode array 36, the condenser lens array 44, and the line sensors 48a to 48c. The unevenness which consists of is included.
[0028]
Therefore, the offset correction unit 74 reads the offset correction data O (x) created in step S1 from the offset correction data storage unit 90,
Doff (x, y) = Dorg (x, y) −O (x)
As described above, offset correction of the correction data acquisition image data Dorg (x, y) is performed to create offset correction image data Doff (x, y) (step S3).
[0029]
Next, the offset corrected image data Doff (x, y) is linearly converted by the linear converter 76 (step S4), and then the connecting process between the line sensors 48a to 48c is performed by the connecting processor 78 (step S4). S5). The offset-corrected image data Doff (x, y) subjected to the linkage processing is converted into logarithmic conversion data that can be added by the logarithmic conversion unit 80 (step S6).
[0030]
The shading correction data calculation unit 92 adds and averages the logarithmically converted offset correction image data Dlog (x, y) = Log (Doff (x, y)) in the direction of the arrow Y for each photoelectric conversion unit. The image data Ds (x) is
Ds (x) = 1 / m · ΣDlog (x, y)
(Step S7). In the above equation, Σ means that the offset correction image data Dlog (x, y) is added by m pixels in the arrow Y direction.
[0031]
Next, the average offset corrected image data Ds (x) is averaged in the direction of the arrow X, and the average value Ds_ave is
Ds_ave = 1 / n · ΣDs (x)
(Step S8). Σ in the above expression means that the average offset correction image data Ds (x) is added by n pixels in the arrow X direction.
[0032]
Finally, the one-dimensional correction data S1 (x) is calculated from the average value Ds_ave and the average offset correction image data Ds (x).
S1 (x) = Ds_ave−Ds (x)
(Step S9). By using the one-dimensional correction data S1 (x) obtained in this way, it is caused by individual characteristic variations of the photoelectric conversion units constituting the laser diode array 36, the condenser lens array 44, and the line sensors 48a to 48c. It is possible to correct shading composed of random high-frequency components. The calculated one-dimensional correction data S1 (x) is stored in the shading correction data storage unit 94.
[0033]
Next, the one-dimensional correction of the offset correction image data Dlog (x, y) logarithmically converted using the one-dimensional correction data S1 (x) is performed, and the one-dimensional correction image data D1 (x, y) is
D1 (x, y) = Dlog (x, y) + S1 (x)
(Step S10).
[0034]
Further, unsharp masking processing is performed on the one-dimensional corrected image data D1 (x, y) (step S11) to obtain low frequency corrected image data Dlow (x, y). The low-frequency corrected image data Dlow (x, y) is, for example, an average of a plurality of one-dimensional corrected image data D1 (x, y) around the focused one-dimensional corrected image data D1 (x, y). Can be obtained at That is, assuming that the number of pixels in the processing range in the arrow X direction is j and the number of pixels in the processing range in the arrow Y direction is k, the low frequency corrected image data Dlow (x, y) is
Dlow (x, y) = 1 / (j · k) · ΣΣD1 (x, y)
Can be obtained as Σ in the above equation means that the one-dimensional corrected image data D1 (x, y) is added by j pixels in the arrow X direction and k pixels in the arrow Y direction. By smoothing the one-dimensional corrected image data D1 (x, y) in this way, it is possible to reduce the influence of noise composed of high-frequency components.
[0035]
Next, the low-frequency corrected image data Dlow (x, y) is averaged in the arrow X direction and the arrow Y direction, and the average value Dave is
Dave = 1 / (n · m) · ΣΣDlow (x, y)
(Step S12). ΣΣ in the above equation means that the low-frequency corrected image data Dlow (x, y) is added by n pixels in the arrow X direction and m pixels in the arrow Y direction.
[0036]
From the average value Dave and the low-frequency corrected image data Dlow (x, y), smoothed two-dimensional correction data S2 (x, y)
S2 (x, y) = Dave-Dlow (x, y)
(Step S13). By using the two-dimensional correction data S2 (x, y) obtained in this way, a low-frequency component having a small change based on the irradiation unevenness of the radiation S or the application unevenness of the storage phosphor in the storage phosphor sheet IP. The shading consisting of can be corrected.
[0037]
The smoothed two-dimensional correction data S2 (x, y) corrects unevenness composed of low-frequency components, and does not include high-frequency components, so adjacent two-dimensional correction data S2 (x, y). It does not fluctuate greatly between. Therefore, by thinning out the two-dimensional correction data S2 (x, y), the two-dimensional correction data S2 (x, y) having a small number of data is obtained (step S14), and the shading correction is performed together with the one-dimensional correction data S1 (x). Store in the data storage unit 94.
[0038]
Next, radiation image information reading processing in the radiation image information recording / reading apparatus 10 will be described with reference to the flowchart shown in FIG.
[0039]
An imaging region, an imaging tube voltage, an imaging dose, etc. are set according to a predetermined imaging menu, and the subject 18 is positioned at a predetermined position of the radiation recording unit 12 in the radiation image information recording / reading apparatus 10, and then the radiation source 20 is driven to emit radiation. The object 18 is irradiated with S. The radiation S transmitted through the subject 18 reaches the stimulable phosphor sheet IP via the phototimer 23 and the grid 24, and the radiographic image information of the subject 18 is recorded (step S21).
[0040]
Next, prior to reading of radiation image information from the stimulable phosphor sheet IP, offset correction data O (x) is created in the same manner as in step S1 described above (step S22). That is, the offset correction data O (x) fluctuates depending on the temperature of the radiation image information recording / reading apparatus 10 and is desirably performed at a fixed time period or immediately before each reading operation. The created offset correction data O (x) is stored in the offset correction data storage unit 90.
[0041]
Next, reading modes such as a high-quality mode by low-speed scanning and a high-speed reading mode by high-speed scanning are set, and the reading unit 14 is driven according to the reading mode to read the radiation image information recorded on the stimulable phosphor sheet IP. Is performed (step S23). The image data Dorg (x, y) related to the radiation image information read from the stimulable phosphor sheet IP is used by the offset correction unit 74 using the offset correction data O (x) acquired from the offset correction data storage unit 90. Offset correction processing is performed, and offset correction image data Doff (x, y)
Doff (x, y) = Dorg (x, y) −O (x)
(Step S24).
[0042]
The offset corrected image data Doff (x, y) is linearly converted by the linear conversion unit 76 (step S25), and then the connection processing unit 78 performs connection processing between the line sensors 48a to 48c (step S26). The logarithmic conversion unit 80 performs logarithmic conversion (step S27),
Dlog (x, y) = Log (Doff (x, y))
As a result, offset-corrected image data Dlog (x, y) is obtained.
[0043]
The logarithmically converted offset correction image data Dlog (x, y) is used in the shading correction unit 82 by using the one-dimensional correction data S1 (x) stored in the shading correction data storage unit 94.
D1 (x, y) = Dlog (x, y) + S1 (x)
Is one-dimensionally corrected to obtain one-dimensional corrected image data D1 (x, y) (step S28).
[0044]
Next, using the two-dimensional correction data S2 (x, y) stored in the shading correction data storage unit 94, the one-dimensional correction image data D1 (x, y)
D2 (x, y) = D1 (x, y) + S2 (x, y)
And corrected image data D2 (x, y) corrected two-dimensionally is obtained (step S29). The shading correction data storage unit 94 stores the two-dimensional correction data S2 (x, y) thinned out in the process of step S14 in order to reduce the amount of data, so that the corrected image data D2 (x , Y), the corrected image data D2 (x, y) is calculated using the interpolated two-dimensional correction data S2 (x, y).
[0045]
The corrected image data D2 (x, y) corrected one-dimensionally and two-dimensionally is subjected to sharpness correction processing in the main scanning direction in the sharpness correction unit 86 (step S30), and then sharpened in the sub-scanning direction. A degree correction process is performed (step S31).
[0046]
Here, for example, the amount of radiation applied to the subject 18 is adjusted by controlling the imaging tube voltage of the radiation source 20 by the imaging region in accordance with the imaging menu. Even if the radiation dose to the subject 18 is the same, the radiation dose that reaches the stimulable phosphor sheet IP varies depending on the imaging region. On the other hand, the spreading characteristics of the photostimulated luminescent light R to be emitted differ depending on the type of the stimulable phosphor sheet IP. Further, since the energy of the excitation light L irradiated to the stimulable phosphor sheet IP differs depending on the reading mode such as high speed scanning or low speed scanning, the stimulated emission light R emitted by the excitation light L is different. The spread characteristics are also different. In this case, the spatial frequency characteristic of the radiation image information obtained from the stimulable phosphor sheet IP depends on the radiation energy accumulated according to the radiation dose irradiated to the stimulable phosphor sheet IP and the emission characteristic thereof.
[0047]
On the other hand, the spatial frequency characteristics of the radiation image information depend on the characteristics of the condensing lens array 44 and the prism 46 in the main scanning direction and the sub scanning direction that guide the stimulated emission light R to the line sensors 48a to 48c, and the line sensor 48a. Dependent on the interval in the main scanning direction and the reading interval in the sub-scanning direction of the photoelectric conversion units constituting .about.48c.
[0048]
Accordingly, the sharpness and the balance of the corrected image data D2 (x, y), which is radiation image information before sharpness correction, with respect to the main scanning direction and the sub-scanning direction are determined according to the shooting menu, the type of the stimulable phosphor sheet IP, the reading It depends on the mode.
[0049]
Therefore, the spatial filter selection unit 98 selects the spatial filter corresponding to the sharpness and the balance according to the recording menu or the recording condition such as the reading mode, which is set by the operator, the type of the stimulable phosphor sheet IP, or the reading condition. It selects from the memory | storage part 96 and supplies to the sharpness correction | amendment part 86. FIG. The sharpness correction unit 86 performs a sharpness correction process on the corrected image data D2 (x, y) using the supplied spatial filter.
[0050]
That is, if the filtering coefficient constituting the spatial filter for correcting the sharpness in the main scanning direction is f (a), the corrected image data D2 (x, y) is
Dx (x, y) = 1 / N · Σf (a) · D2 (x + a, y)
As described above, image data Dx (x, y) corrected in the main scanning direction is obtained (step S30). a is an integer representing a pixel range used for performing the filtering process, and Σ represents an addition process in the pixel range a. N = Σf (a).
[0051]
FIG. 6 shows the filtering coefficient f (a) as
f (a) = (-2 3 10 3 -2)
Represents the relationship between the spatial frequency (cycle / mm) and the response to the main scanning direction. When this filtering coefficient f (a) is used, an effect of increasing the sharpness of the corrected image data D2 (x, y) before sharpness correction in the main scanning direction can be obtained.
[0052]
In addition, the sharpness correction unit 86 corrects the sharpness in the main scanning direction with g (b) as the filtering coefficient constituting the spatial filter for correcting the sharpness in the sub-scanning direction selected from the spatial filter storage unit 96. The image data Dx (x, y)
Dy (x, y) = 1 / M · Σg (b) · Dx (x, y + b)
To obtain image data Dy (x, y) corrected in the sub-scanning direction (step S31). b is an integer representing a pixel range used for performing the filtering process, and Σ represents an addition process in the pixel range b. M = Σg (b).
[0053]
FIG. 7 shows the filtering coefficient g (b) as
[0054]
[Expression 1]

[0055]
The relationship between the spatial frequency (cycle / mm) and the response in the main scanning direction is shown. When this filtering coefficient g (b) is used, an effect of weakening the sharpness of the image data Dx (x, y) in the sub-scanning direction can be obtained.
[0056]
As described above, by using a spatial filter having appropriate filtering coefficients f (a) and g (b), an image subjected to substantially the same sharpness correction processing in the main scanning direction and the sub-scanning direction. Data Dy (x, y) can be obtained.
[0057]
Next, after normalizing the image data Dy (x, y) (step S32), image processing such as gradation correction and color correction is performed (step S33). In this case, desired image processing can be performed without considering the difference in sharpness between the main scanning direction and the sub-scanning direction. The radiographic image information subjected to the image processing in this way is provided to an external diagnostic device or the like. Note that the image processing in step S33 may be performed not by the radiation image information recording / reading apparatus 10 but by an external image processing apparatus.
[0058]
On the other hand, when the reading process of the radiation image information from the stimulable phosphor sheet IP is completed, the erasing light source unit 16 is driven, and the erasing light Q output from the fluorescent lamp 68 is irradiated on the entire surface of the stimulable phosphor sheet IP. As a result, the remaining radiation energy is removed, and it is used for the next imaging.
[0059]
Here, in the above-described embodiment, the one-dimensional filtering process for increasing the sharpness in the main scanning direction and the one-dimensional filtering process for decreasing the sharpness in the sub-scanning direction are performed. Alternatively, one sharpness in the sub-scanning direction may be fixed as a reference, and only the other sharpness may be adjusted to match the one sharpness.
[0060]
Further, a filtering coefficient fg (a, b), which is an inner product of the filtering coefficient f (a) in the main scanning direction and the filtering coefficient g (b) in the sub-scanning direction, for example,
[0061]
[Expression 2]

[0062]
And using this filtering coefficient fg (a, b), the image data D (x, y) is
D (x, y) = 1 / P · ΣΣfg (a, b) · D2 (x + a, y + b)
As a result, the sharpness correction processing in the main scanning direction and the sub-scanning direction can be performed simultaneously. ΣΣ represents an addition process in the pixel ranges a and b of the filtering process in the main scanning direction and the sub-scanning direction. Further, P = ΣΣfg (a, b).
[0063]
In the radiation image information recording / reading apparatus 10 according to the above-described embodiment, the line sensors 48a to 48c are moved in the sub-scanning direction with respect to the stimulable phosphor sheet IP, but the radiation image information is read. It can also be applied to a radiation image information recording / reading apparatus configured to move the stimulable phosphor sheet IP in the sub-scanning direction with respect to the line sensors 48a to 48c.
[0064]
【The invention's effect】
As described above, according to the present invention, desired information can be obtained for each of the main scanning direction and the sub-scanning direction from image information read using a sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction. Sharpness correction processing can be performed, and thereby good image information having substantially the same sharpness in the main scanning direction and the sub-scanning direction can be obtained.
[Brief description of the drawings]
FIG. 1 is a schematic configuration diagram of a radiation image information recording / reading apparatus according to an embodiment.
FIG. 2 is a block diagram of a processing circuit in a reading unit of the radiation image information recording / reading apparatus according to the embodiment.
FIG. 3 is a flowchart showing a correction data creation procedure of the radiation image information recording / reading apparatus according to the present embodiment;
FIG. 4 is an explanatory diagram of a coordinate system set in the stimulable phosphor sheet.
FIG. 5 is a flowchart of a radiographic image information shading correction process and a sharpness correction process.
FIG. 6 is an explanatory diagram of a relationship between a spatial frequency and a response in a filter capable of enhancing sharpness.
FIG. 7 is an explanatory diagram of a relationship between a spatial frequency and a response in a filter that can weaken the sharpness.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 ... Radiation image information recording / reading apparatus 12 ... Radiation recording part
14 ... Reading unit 16 ... Erasing light source unit
18 ... Subject 20 ... Radiation source
48a to 48c ... line sensor 74 ... offset correction unit
82: Shading correction unit 86 ... Sharpness correction unit
88 ... Offset correction data calculation unit 90 ... Offset correction data storage unit
92: Shading correction data calculation unit
94: Shading correction data storage unit
96 ... Spatial filter storage unit 98 ... Spatial filter selection unit
IP ... Storage phosphor sheet

Claims

A sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction is moved relative to the recording medium in a sub-scanning direction substantially orthogonal to the main scanning direction, and image information recorded on the recording medium is transferred. In the sharpness correction method in the image reading device that is read two-dimensionally by the photoelectric conversion unit,
The sharpness in the main scanning direction of the image information obtained by reading and the sharpness in the subscanning direction of the image information are substantially the same with respect to the main scanning direction and the subscanning direction, respectively. A sharpness correction method in an image reading apparatus, wherein a desired sharpness correction process is performed.

The method of claim 1, wherein
The sharpness correction method in the image reading apparatus, wherein the sharpness correction process is a one-dimensional filtering process using a spatial filter set in each of the main scanning direction and the sub-scanning direction.

The method of claim 1, wherein
The sharpness correction method in the image reading apparatus, wherein the sharpness correction process is a two-dimensional filtering process using a spatial filter set in the main scanning direction and the sub-scanning direction.

The method according to claim 2 or 3,
The sharpness correction method in an image reading apparatus, wherein the spatial filter is selected according to a recording condition or a reading condition of the image information with respect to the recording medium.

A sensor in which a plurality of photoelectric conversion units are arranged in the main scanning direction is moved relative to the recording medium in a sub-scanning direction substantially orthogonal to the main scanning direction, and image information recorded on the recording medium is transferred. In the sharpness correction device in the image reading device that is read two-dimensionally by the photoelectric conversion unit,
The sharpness in the main scanning direction of the image information obtained by reading and the sharpness in the subscanning direction of the image information are substantially the same with respect to the main scanning direction and the subscanning direction, respectively. A spatial filter storage unit that stores a spatial filter that performs a desired sharpness correction process;
A spatial filter selection unit that selects the spatial filter from the spatial filter storage unit according to a recording condition or a reading condition of the image information with respect to the recording medium;
A sharpness correction unit that performs a sharpness correction process on the image information by the selected spatial filter;
A sharpness correction apparatus in an image reading apparatus.