JP4138199B2 - Image information reading method and apparatus - Google Patents
Image information reading method and apparatus Download PDFInfo
- Publication number
- JP4138199B2 JP4138199B2 JP2000082548A JP2000082548A JP4138199B2 JP 4138199 B2 JP4138199 B2 JP 4138199B2 JP 2000082548 A JP2000082548 A JP 2000082548A JP 2000082548 A JP2000082548 A JP 2000082548A JP 4138199 B2 JP4138199 B2 JP 4138199B2
- Authority
- JP
- Japan
- Prior art keywords
- scanning
- sub
- pixel density
- sheet
- image information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Input (AREA)
- Image Processing (AREA)
- Facsimile Scanning Arrangements (AREA)
- Editing Of Facsimile Originals (AREA)
- Facsimile Image Signal Circuits (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
- Radiography Using Non-Light Waves (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は画像情報読取方法および装置に関し、詳細には、読取られた画像信号の画素密度を変更可能とした画像情報読取方法および装置に関するものである。
【0002】
【従来の技術】
放射線(X線、α線、β線、γ線、電子線、紫外線等)を照射すると、この放射線エネルギーの一部が蓄積され、その後可視光等の励起光を照射すると蓄積されたエネルギーに応じて輝尽発光を示す蓄積性蛍光体(輝尽性蛍光体)を利用して、人体等の被写体の放射線画像情報を一旦シート状の蓄積性蛍光体(蓄積性蛍光体シート)に記録し、この蓄積性蛍光体シートにレーザビーム等の光ビームを走査して信号光である画像情報に応じた輝尽発光光を生じせしめ、得られた輝尽発光光を光電読取手段により光電的に読取ってアナログ画像信号を得、このアナログ画像信号を所定の間隔でサンプリングするとともに量子化することにより所定の画素密度のデジタル画像信号を得ることができ、さらにこのデジタル画像信号に基づき写真感光材料等の記録媒体、CRT等の表示装置に被写体の放射線画像を可視像として出力させる放射線画像情報記録再生システムがすでに知られている(特開昭55-12429号、同56-11395号、同56-11397号など)。
【0003】
このシステムは、従来の銀塩写真を用いる放射線写真システムと比較して極めて広い放射線露出域に亘って画像を記録しうるという実用的な利点を有している。
【0004】
また、上述した輝尽発光光を光電的に読取る方法として、蓄積性蛍光体シートの両面に各別に光電読取手段を配設して、蓄積性蛍光体シートの両面または片面にのみ励起光を照射し、この励起光の照射により蓄積性蛍光体シートの両面からそれぞれ発せられる輝尽発光光を各光電読取手段により光電的に検出する放射線画像読取装置が提案されている(例えば、特開昭55-87970号)。このような放射線画像読取装置は、蓄積性蛍光体シートに1つの放射線画像が蓄積記録され、かつ蓄積性蛍光体シートの両面から読取りを行うように光電読取手段を配し、これにより輝尽発光光を検出するため、集光効率が向上し、得られた2つの画像信号を所定の加算比で加算することにより、ノイズ成分の検出位置がシートの表裏で分散するため、片面のみから得られた画像信号よりも相対的にS/Nがより改善された加算画像信号を得ることができる。
【0005】
さらに、シートの一方の面から得られた単一の画像信号またはシートの両面からそれぞれ得られた2つの画像信号に対して、画像信号(加算画像信号を含む)のS/Nを高くするような周波数応答特性を有するフィルタによってフィルタリング処理を施した後に、加算画像信号を得るようにした画像の重ね合わせ方法が提案されている(特開平7-287330号)。この方法によれば、撮影時において被写体に照射された放射線量を求め、この放射線量に基づいて画像信号にフィルタリング処理を施す際のパラメータ(フィルタ係数)を求めているため、被写体に照射される放射線量に応じて最適な画質の画像を表す画像信号または重ね合わせ画像を表す加算画像信号を得ることができる。さらにまた、この方法は、画像信号の全体に周波数特性を変化させる処理を施すようにしたため、フーリエ変換等の周波数変換を行う必要がなくなり、計算量を少なくすることができるものである。
【0006】
【発明が解決しようとする課題】
ところで上記放射線画像情報記録再生システムの画像読取部において、あるいはこの画像読取部が単独で構成された放射線画像情報読取装置において、得られた単一の画像信号または加算画像信号の画素密度を、上述した所定の画素密度とは異なるものとしたいという要望がある。
【0007】
これは例えば、集団検診などのように大量の放射線画像の撮影が行われる場合等においては、撮影された放射線画像の読取り速度を向上することが要求される一方、再生される画像に求められる画質は、再検査を要するか否かの判断を行うことができる程度のものであればよく、過度に高画素密度で画像情報を読取る必要はない等の事情によるものである。勿論これとは反対に、読取り速度は低下しても高画質の画像を再生しうるように高画素密度で画像情報を読み取ることが求められる場合もある。
【0008】
このように予め設定されている所定の画素密度による読取りの他に、他の画素密度での読取りをも行うには、単純には光ビームの主走査速度と副走査速度とを変更することが考えられる。そして、この光ビームの主走査速度を変更するときは、光ビームの主走査を行う走査光学系(例えばポリゴンミラーやガルバノメータミラー等)の駆動速度を変更する必要がある。
【0009】
しかしながら、走査光学系の駆動速度を変更したときは、その慣性力の影響により駆動速度が安定するまで一定の時間を要し、また全ての速度域において安定した駆動が保証されていない場合もある。このため、光ビームの主走査速度を変更させずに画素密度を変更可能とすることが望ましい。また、このように画素密度を変更した場合に、主走査方向および副走査方向の各画素密度が同一の割合で変更される必要もある。
【0010】
また、画素密度を変更する際には、できるだけ高速に画素密度変換処理を行うことが望まれている。
【0011】
さらに、単に画素密度を変更しただけでは、以下のような問題が生じることも考えられる。
【0012】
すなわち、
(I) 画素密度の変更によりシートから発光される1画素当たりの信号光のエネルギーが、元の画素密度のときの1画素当たりの信号光のエネルギーとは異なるものとなるため、変更前後の画素密度比が大きい場合には再生される画像全体の濃度(または輝度)に変化を生じて、診断性能に影響を与える場合がある。
【0013】
(II) アナログ画像信号に対してシェーディング補正を行う場合には、画素密度変更前後で、補正すべきシェーディングの特性に変化を生じて、シェーディングを適正に補正できない場合がある。
【0014】
(III)アナログ画像信号を対数増幅する場合には、画素密度変更前後で、周波数伝達特性に変化を生じる場合がある。
【0015】
(IV) アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合には、画素密度変更前後で、ナイキスト周波数に変化を生じるため、適正にエリアジングノイズをカットできない場合がある。
【0016】
これらの問題は、上述した放射線画像情報記録再生システムによる、放射線画像情報が記録された蓄積性蛍光体シートからの放射線画像情報の読取りにおいてのみ生じる問題ではなく、例えば画像がプリントされた媒体を光ビームで走査して、この媒体にプリントされた画像に応じて発せられる反射光を信号光として読取る方法や装置においても同様に生じうる。
【0017】
本発明は上記事情に鑑みなされたものであって、光ビームの主走査速度を変更させずに画素密度を変更することができる画像情報読取方法および装置を提供することを第1の目的とするものである。
【0018】
また、本発明は、光ビームの主走査速度を変更させずに画素密度を変更することができる両面読取り用の画像情報読取方法および装置を提供することを第2の目的とするものである。
【0019】
【課題を解決するための手段】
本発明の画像情報読取方法および装置は、副走査の速度およびサンプリングの間隔を変更することにより、主走査速度を変更することなく、任意の画素密度のデジタル画像信号を得るものである。
【0020】
すなわち、本発明の第1の画像情報読取方法は、画像情報が記録された被走査シートに光ビームを一定の方向に一定の速度で繰り返し主走査しつつ該主走査方向に対して略直交する方向に副走査することにより、前記光ビームを該被走査シートの略全面に走査せしめ、該光ビームの走査により該被走査シートから発せられる信号光を光電的に検出し、得られたアナログ画像信号をサンプリングおよび量子化して所定の画素密度のデジタル画像信号を得る画像情報読取方法において、
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍することを特徴とするものである。
【0021】
ここで主走査方向の画素数をa/m(a>0)倍、副走査方向の画素数をa/n倍する画素密度変換処理をさらに加えて行なうようにするのが好ましい。画像の縦横比を元の画像と同一にすることができるからである。
【0022】
また本発明の第2の画像情報読取方法は、シートの両面からそれぞれ読み取って得られた2つの画像信号の加算画像信号に関するものであり、画像情報が記録された被走査シートに光ビームを一定の方向に一定の速度で繰り返し主走査しつつ該主走査方向に対して略直交する方向に副走査することにより、前記光ビームを該被走査シートの略全面に走査せしめ、該光ビームの走査により該被走査シートの両面から発せられる信号光を光電的に検出し、得られた2つのアナログ画像信号をサンプリングおよび量子化して2つのデジタル画像信号を得、該2つのデジタル画像信号を加算して所定の画素密度の加算画像信号を得る画像情報読取方法において、
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍し、さらに前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n倍とする画素密度変換処理を行なうことを特徴とするものである。
【0023】
ここで上記第1および第2の発明において、被走査シートには、上述した放射線画像情報記録再生システムにおいて使用される、放射線画像が蓄積記録された蓄積性蛍光体シートのみならず、写真プリント等の反射原稿やフイルムなどの透過原稿等も含む。したがって、被走査シートから発せられる信号光には、光ビームにより励起されて蓄積性蛍光体シートから発せられる輝尽発光光の他、反射原稿から発せられる反射光や透過原稿から発せられる透過光をも含む。
【0024】
また、サンプリングの間隔を変更する方法としては、基準クロックに対して分周回路を用いて分周したクロックを作り、この分周したクロックによりサンプリングし、あるいは基準クロックに対してPLLにより周期の異なるクロックを作成してこの新たなクロックによりサンプリングし、あるいは周期の異なる複数のサンプリング用クロックを切り換えて適用してサンプリングするなどの方法を適用することができる。
【0025】
さらに、加算画像信号を得る方法としては、単純に2つのデジタル画像信号の相対応する画素同士を加算する方法のみならず、上記特開平7-287330号に記載されたように、2つのデジタル画像信号に対して加算画像信号のS/Nを高くするような周波数応答特性を有するフィルタによってフィルタリング処理を施してから加算する方法を採用することができる。
【0026】
一方、画素密度変換処理としては、画像信号の主走査方向に対して1次元マスク演算を施すことによる方法、画素密度に応じて画素を間引く方法、Bスプライン補間演算、Cubic スプライン補間演算等の高次の補間演算による方法(例えば特開平8-16767号、同9-321981号)、線形補間演算による方法(例えば特開平9-50516号)を採用することができる。
【0027】
なお、この画素密度変換を行う際には、m,nに応じて1次元マスク演算のパラメータ、画素を間引く際のパラメータおよび補間演算の演算式のパラメータを変更することが好ましい。ここで、1次元マスク演算のパラメータはマスク係数であり、画素を間引く際のパラメータは画素の間引き間隔である。なお、補間演算の演算式のパラメータは、補間演算を行う際にどのような補間演算(例えばBスプライン補間演算、Cubic スプライン補間演算、線形補間演算)の演算式を用いるかを表すものである。
【0028】
なお、本発明の第2の画像情報読取方法においては、前記2つのデジタル画像信号に対して前記画素密度変換処理を行うことが好ましい。ここで、2つのデジタル画像信号に対して画素密度変換処理を行うとは、2つのデジタル画像信号を加算して加算画像信号を得る前に画素密度変換を行うことをいう。
【0029】
また、本発明の各画像情報読取方法においては、前記m,nに応じて、以下の(1)から(8)の特性のうち少なくとも1つを変更することが好ましい。
【0030】
(1)前記被走査シートを走査する光ビームの径
(2)前記被走査シートを走査する光ビームのパワー
(3)前記信号光を検出する感度
(4)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(5)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(6)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(7)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数
(8)前記画素密度変換処理におけるフィルタリング処理のパラメータ(フィルタリング処理として間引き処理を行なう場合の間引きのパラメータをも含む)
ここで、上記(1)から(8)の変更事項は、少なくとも1つの該当する事項について適用するものであり、必ずしも(1)から(8)までの全てを構成要件として備える必要はない。すなわち、例えばシェーディング補正を行うことがない画像情報読取方法においては、上記(4)および(5)の構成を備える必要はないものである。
【0031】
さらにこれら(1)から(8)の事項について、少なくとも1つ変更するものであればよく、いずれか1つのみを構成要件として備えるもの、これらのうち2つ以上を備えるものであってもよい。
【0032】
なお、上記(1)から(3)の変更事項は、上述した (I) の問題、すなわち画素密度の変更によりシートから発光される1画素当たりの信号光のエネルギーが、元の画素密度のときの1画素当たりの信号光のエネルギーとは異なるものとなるため、変更前後の画素密度比が大きい場合には再生される画像全体の濃度(または輝度)に変化を生じて、診断性能に影響を与える場合において適用する事項であり、具体的には、画素密度を高密度側に変更したときは、(1)光ビーム径を縮小し、(2)光ビームのパワーを弱く、および/または(3)信号光を検出する感度を高くし、一方、画素密度を低密度側に変更したときは、(1)光ビーム径を拡大し、(2)光ビームのパワーを強く、および/または(3)信号光を検出する感度を低くすればよい。
【0033】
同様に、画素密度を高密度側に変更したときは、(6)の事項は、周波数応答特性をより高周波数帯域まで広げるように変更し、(7)の事項も同様にカットオフ周波数を高周波数側に変更すればよく、画素密度を低密度側に変更したときは、(6)の事項は、カットオフ周波数をより低周波数側とするように変更し、(7)の事項も同様にカットオフ周波数を低周波数側に変更すればよい。
【0034】
上記(8)の事項は、例えばフィルタリング処理のパラメータ(マスク演算の係数)を、高密度のときよりも低密度のときの方がマスク演算したときのカットオフ周波数を低周波側にするように変更すればよい。なお、間引き処理のパラメータを適用する場合は、上記フィルタリング処理のパラメータを適用した場合の画素密度の高低に対する変更方向とは反対の方向に変更すればよい。
【0035】
上記(4)の事項は、画素密度に応じて定性的に変動する特性ではないため、予め代表的な画素密度毎にシェーディング補正用データを準備しておき、変更された画素密度に対応するシェーディング補正用データを、これらの準備されたデータのうちから選択し、または一致する画素密度のデータが準備されていない場合は、準備されているデータのうち少なくとも2種類の画素密度に対応するデータを用いた補間処理により、変更された画素密度に対応するシェーディング補正用データを求めて適用すればよい。
【0036】
上記(5)の事項は、デジタル画像信号に対してリアルタイムでシェーディング補正を行う際に適用されるものであり、予め準備されたシェーディング補正用データを画素密度に応じたタイミングで出力して、得られたデジタル画像信号に対してシェーディング補正を施すものである。具体的には、画素密度を高密度側に変更したときには、サンプリング速度が早くなるためこのサンプリング速度に応じてシェーディング補正用データを出力するタイミングを早くすればよい。一方、画素密度を低密度側に変更したときには、サンプリング速度が遅くなるためこのサンプリング速度に応じてシェーディング補正用データを出力するタイミングを遅くすればよい。
【0037】
ここで、シェーディングとは、光ビームを走査するための光偏向器(ポリゴンミラー、ガルバノメータミラー等)の反射面の反射率ムラによる走査ビームの強度ムラ、また光偏向器の偏向速度のバラツキによる光ビームの走査速度ムラ、あるいは主走査方向に沿って配設される光電検出手段の主走査方向における感度ムラによる検出ムラ等により、光電読取手段から得られるアナログ画像信号の変動(光検出効率の部分的な低下)を意味するものである。またシェーディング補正用データとは、基準となる被走査シート(例えば放射線が一様に照射された蓄積性蛍光体シート)を用いる等して予め求められたシェーディングの特性をいうものである(特開昭61-189763 号、同62-47259号,同62-47261号、同64-86759号、特開平2-58973 号等参照)。
【0038】
なお、前述した変更事項(1)、(2)、(3)、(6)、(7)、(8)においてもそれぞれ、予め画素密度毎の特性((1)にあってはビーム径、(2)にあってはビームパワー、(3)にあっては感度、(6)にあっては周波数応答特性、(7)にあってはカットオフ周波数の異なる周波数応答特性、(8)にあってはフィルタリング処理のパラメータ)を準備しておき、変更された画素密度に対応する特性を選択することにより、変更するものであってもよい。
【0039】
また(4)シェーディング補正用データの変更方法としては上述したような、予め準備された複数の補正用データのうち選択する方式の他、予め設定されている画素密度に対応して準備されている1つのシェーディング補正用データを、変更された後の画素密度に応じてサンプリングすることによって求め、この新たに求められたシェーディング補正用データを用いるようにしてもよい。
【0040】
なお、mおよびnの値毎に例えば、「高密度」(例えばm=n=0.5のとき)、「標準密度」(元の所定の画素密度(m=n=1に相当))、「低密度」(m=n=2のとき)という対応関係を設定しておいてもよい。
【0041】
また、選択され得るm,nの値ごとに予め設定されたシェーディング補正データのうち、選択されたm,nの値に応じたシェーディング補正データを選択することにより、シェーディング補正データを変更するようにしてもよく、この場合、下記の2つの方法[I],[II]のうちいずれか一方を適用するのが好ましい。すなわち、
[I]m,nの値ごとに予め設定されたシェーディング補正データを第1の記憶媒体に予め記憶させておき、選択されたm,nの値に対応するシェーディング補正データを、選択の都度、第1の記憶媒体から第2の記憶媒体に転送し、この転送されたシェーディング補正データを第2の記憶媒体から読み出すことにより、シェーディング補正データの選択を行なう。
[II]m,nの値ごとに予め設定されたシェーディング補正データを第1の記憶媒体に予め記憶させておき、起動時など所望とするときに、これらのシェーディング補正データの全てを、第1の記憶媒体から、m,nの値ごとに互いに異なるアドレスに対応した記憶領域を有する第2の記憶媒体にm,nの値に対応させて転送し、選択されたm,nの値に対応するシェーディング補正データを、選択されたm,nの値に対応するアドレスから読み出すことにより、シェーディング補正データの選択を行なう。
【0042】
[I]の方法によれば、第2の記憶媒体(メモリ)の容量を少なくすることができ、ハードウェア構成を簡単にすることができる。一方、[II]の方法によれば、ソフトウェア構成を簡単にすることができ、第2の記憶媒体からの補正データの読出し速度を高速化することができる。
【0043】
本発明の第3の画像情報読取方法は、画像情報が記録された被走査シートに光ビームを一定の方向に一定の速度で繰り返し主走査しつつ該主走査方向に対して略直交する方向に副走査することにより、前記光ビームを該被走査シートの略全面に走査せしめ、該光ビームの走査により該被走査シートから発せられる信号光を光電的に検出し、得られたアナログ画像信号をサンプリングおよび量子化して所定の画素密度のデジタル画像信号を得る画像情報読取方法において、
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍し、これらm,nに応じて、以下の(1)から(5)の特性のうち少なくとも1つを変更して、前記所定の画素密度に対して1/(m×n)倍の画素密度のデジタル画像信号を得ることを特徴とするものである。
【0044】
(1)前記被走査シートを走査する光ビームの径
(2)前記信号光を検出する感度
(3)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(4)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(5)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数
(1)〜(5)の事項は上述した説明と同様である。
【0045】
本発明の第1の画像情報読取装置は、本発明の第1の画像情報読取方法を実施するための装置であって、画像情報が記録された被走査シートに、光源から出射された光ビームを一定の方向に一定の速度で繰り返し主走査しする主走査手段と、この光ビームの主走査の方向に対して略直交する方向に前記シートが相対的に副走査されるように該シートまたは前記光ビームを副走査する副走査手段と、前記光ビームの走査により前記被走査シートから発せられる信号光を光電的に検出する光電検出手段と、検出されたアナログ画像信号をサンプリングおよび量子化して所定の画素密度のデジタル画像信号を得るA/D変換手段とを備えた画像情報読取装置において、
前記副走査手段による前記副走査の速度を、入力されたmに応じてm倍に変更せしめる副走査速度変更手段と、前記A/D変換手段による前記サンプリングの間隔を、入力されたnに応じてn倍に変更せしめるサンプリング間隔変更手段とをさらに備えたことを特徴とするものである。
【0046】
ここで、主走査方向の画素数をa/m倍、副走査方向の画素数をa/n倍とする画素密度変換処理を行なう画素密度変換処理手段をさらに備えた構成とするのが好ましい。
【0047】
また、本発明の第2の画像情報読取装置は、本発明の第2の画像情報読取方法を実施するための装置であって、画像情報が記録された被走査シートに、光源から出射された光ビームを一定の方向に一定の速度で繰り返し主走査する主走査手段と、この光ビームの主走査の方向に対して略直交する方向に前記シートが相対的に副走査されるように該シートまたは前記光ビームを副走査する副走査手段と、前記光ビームの走査により前記被走査シートの両面から発せられる信号光を光電的に検出する光電検出手段と、検出された2つのアナログ画像信号をサンプリングおよび量子化して2つのデジタル画像信号を得るA/D変換手段と、該2つのデジタル画像信号を加算して所定の画素密度の加算画像信号を得る加算手段とを備えた画像情報読取装置において、
前記副走査手段による前記副走査の速度を、入力されたmに応じてm倍に変更せしめる副走査速度変更手段と、
前記A/D変換手段による前記サンプリングの間隔を、入力されたnに応じてn倍に変更せしめるサンプリング間隔変更手段と、
前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n倍とする画素密度変換処理を行なう画素密度変換処理手段とを備えたことを特徴とするものである。
【0048】
なお、本発明の第2の画像情報読取装置においては、前記画素密度変換処理手段は、前記2つのデジタル画像信号に対して画素密度変換処理を行う手段であることが好ましい。
【0049】
また、本発明の各画像情報読取装置においては、前記m,nに応じて、以下の(1)から(8)の特性のうち少なくとも1つを変更せしめる特性変更手段をさらに備えることが好ましい。
【0050】
(1)前記被走査シートを走査する光ビームの径
(2)前記被走査シートを走査する光ビームのパワー
(3)前記信号光を検出する感度
(4)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(5)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(6)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(7)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数
(8)前記画素密度変換処理におけるフィルタリング処理のパラメータ(フィルタリング処理として間引き処理を行なう場合の間引きのパラメータをも含む)
なお、副走査速度変更手段およびサンプリング間隔変更手段は、特性変更手段の一部をなす構成であってもよい。
【0051】
また特性変更手段が、シェーディング補正データを変更するものを含む場合、選択され得るm,nの値ごとに予め設定されたシェーディング補正データのうち、選択されたm,nの値に応じたシェーディング補正データを選択することにより、シェーディング補正データを変更するものとしてもよく、この場合、下記の2つの方法[I],[II]のうちいずれか一方を適用するのが好ましい。すなわち、
[I]m,nの値ごとに予め設定されたシェーディング補正データを第1の記憶媒体に予め記憶させておき、選択されたm,nの値に対応するシェーディング補正データを、選択の都度、第1の記憶媒体から第2の記憶媒体に転送し、この転送されたシェーディング補正データを第2の記憶媒体から読み出すことにより、シェーディング補正データの選択を行なう。
[II]m,nの値ごとに予め設定されたシェーディング補正データを第1の記憶媒体に予め記憶させておき、起動時など所望とするときに、これらのシェーディング補正データの全てを、第1の記憶媒体から、m,nの値ごとに互いに異なるアドレスに対応した記憶領域を有する第2の記憶媒体にm,nの値に対応させて転送し、選択されたm,nの値に対応するシェーディング補正データを、選択されたm,nの値に対応するアドレスから読み出すことにより、シェーディング補正データの選択を行なう。
【0052】
[I]の方法によれば、第2の記憶媒体(メモリ)の容量を少なくすることができ、ハードウェア構成を簡単にすることができる。一方、[II]の方法によれば、ソフトウェア構成を簡単にすることができ、第2の記憶媒体からの補正データの読出し速度を高速化することができる。
【0053】
本発明の第3の画像情報読取装置は、画像情報が記録された被走査シートに、光源から出射された光ビームを一定の方向に一定の速度で繰り返し主走査する主走査手段と、この光ビームの主走査の方向に対して略直交する方向に前記シートが相対的に副走査されるように該シートまたは前記光ビームを副走査する副走査手段と、前記光ビームの走査により前記被走査シートから発せられる信号光を光電的に検出する光電検出手段と、検出されたアナログ画像信号をサンプリングおよび量子化して所定の画素密度のデジタル画像信号を得るA/D変換手段とを備えた画像情報読取装置において、
前記所定の画素密度に対して、主走査方向について1/m倍、副走査方向について1/n倍となる画素密度に対応した値m、nの入力を受ける画素密度入力手段と、
前記副走査手段による前記副走査の速度をm(m>0)倍に変更せしめる副走査速度変更手段と、
前記A/D変換器による前記サンプリングの間隔をn(n>0)倍に変更せしめるせしめるサンプリング間隔変更手段と、
m,nに応じて、以下の(1)から(5)のうち少なくとも1つを変更せしめる特性変更手段とを備えたことを特徴とするものである。
【0054】
(1)前記被走査シートを走査する光ビームの径
(2)前記信号光を検出する感度
(3)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(4)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(5)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数
【0055】
【発明の効果】
本発明の画像情報読取方法および装置によれば、副走査の速度をm(m>0))倍にすることにより被走査シートを走査する光ビームの走査線(主走査線)数が1/m倍に増減するため副走査方向における画素密度を1/m倍にすることができ、またサンプリング間隔をn(n>0)倍にすることにより主走査方向における画素密度を1/n倍にすることができ、したがって最終的に得られる画像信号(両面読みの場合は加算画像信号)の1/(m×n)2 倍の画素密度の画像信号を得ることができる。また、画素密度を変更したときも光ビームの主走査速度を変更する必要がないため、この光ビームを走査する走査光学系(例えばポリゴンミラーやガルバノメータミラー等)の駆動速度も変更する必要がなく、したがって、走査光学系の駆動速度を変更した場合に生じる駆動速度安定までの読取り不能時間を解消することができる。
【0056】
さらに、主走査方向の画素数をa/m(a>0)倍、副走査方向の画素数をa/n倍とする画素密度変換処理を行うようにした発明においては、最終的に得られる画像信号(両面読みの場合は加算画像信号)の主走査方向の画素密度をa/(m×n)倍、副走査方向の画素密度もa/(m×n)倍にすることができ、主走査方向および副走査方向の画素密度の変更割合が同一になり、画像サイズの縦横比(主副比)を維持しつつも、光ビームの主走査速度を変更させることなく、予め設定されている画素密度の画像信号だけでなく、この画素密度の(a/(m×n))2 倍の画素密度の画像信号をも得ることができる。
【0057】
とくに、被走査シートの両面から信号を読み取る両面読取りを行う(本発明の第2の画像情報読取方法および装置の)場合には、光ビームによるエネルギーをシートの裏面側まで十分に付与するために、片面側からのみ読取る方式に比して走査速度を遅くする必要があるが、その場合においても画素密度を粗くすることが許容されるときは副走査速度を速めることができ、これにより1つのシート当たりの走査時間をより短縮することができる。
【0058】
また、両面読みの場合は、2つのデジタル画像信号に対して画素密度変換処理を施すことにより、画素密度変換処理がなされた画像信号を加算して加算画像信号を得ることができる。このため、デジタル画像信号を加算する際の加算演算の演算量を低減することができ、これにより加算演算の演算時間を短縮して処理を高速に行うことができる。
【0059】
さらに本発明の画像情報読取方法および装置によれば、読取りの画素密度を変更したとき、該当する場合に生じる各種の問題を解決することもできる。
【0060】
すなわち、(I)画素密度の変更によりシートから発光される1画素当たりの信号光のエネルギーが、元の画素密度のときの1画素当たりの信号光のエネルギーとは異なるものとなるため、変更前後の画素密度比が大きい場合には再生される画像全体の濃度(または輝度)に変化を生じて、診断性能に影響を与える場合があるが、本発明の画像情報読取方法および装置によれば、画素密度に応じて(1)前記被走査シートを走査する光ビームの径を変更することによりシートの単位面積に付与するエネルギーを変更することができ、また画素密度に応じて(2)光ビームのパワーを変更することによりシートの単位面積に付与するエネルギーを変更することができ、さらに画素密度に応じて(3)前記信号光を検出する感度を変更することにより信号光の有するエネルギーが変化しても得られる画像信号のレベルが変化するのを抑制することができる。
【0061】
また、(II)アナログ画像信号に対してシェーディング補正を行う場合には、画素密度変更前後で、補正すべきシェーディングの特性に変化を生じて、シェーディングを適正に補正できない場合があるが、本発明の画像情報読取方法および装置によれば、画素密度に応じて(4)シェーディング補正用データを変更することにより、シェーディングの特性の変化に対応したシェーディング補正を行うことができ、画素密度に応じて(5)シェーディング補正用データのメモリからの出力のタイミングを変更することにより、サンプリング速度の変化に対応したシェーディング補正を行うことができる。
【0062】
さらにまた、 (III)アナログ画像信号を対数増幅する場合には、画素密度変更前後で、周波数伝達特性に変化を生じる場合があるが、本発明の画像情報読取方法および装置によれば、画素密度に応じて(6)アナログ画像信号を対数増幅するときの周波数伝達特性を変更することにより、得られた画像信号の周波数特性が変化するのを抑制することができる。
【0063】
また、(IV)アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合には、画素密度変更前後で、ナイキスト周波数に変化を生じるため、適正にエリアジングノイズをカットできない場合があるが、本発明の画像情報読取方法および装置によれば、画素密度に応じて(7)エリアジングノイズを除去するフィルタリングのカットオフ周波数を変更することにより、得られた画像信号の周波数特性が変化するのを抑制することができる。
【0064】
さらにまた画素密度変換処理を行なう場合における(8)画素密度変換処理におけるフィルタリング処理のパラメータ(マスク演算の係数)を変更することにより、画素密度変換する際のエリアジングノイズを一定以下に抑制することができる。またスプライン補間演算等の補間処理により生じ得る誤差を低減することもできる。
【0065】
【発明の実施の形態】
以下、本発明の画像情報読取装置の具体的な実施の形態について、図面を用いて説明する。
【0066】
図1は本発明の画像情報読取装置の一実施形態を示す図である。図示の画像情報読取装置は、図示しないモータにより回転せしめられるエンドレスベルト9a上に放射線画像情報が蓄積記録された蓄積性蛍光体シート(以下、シートという)1が配置され、シート1の上方には、シート1を励起するレーザー光11を発するレーザー光源10と、モーター20により回転され、レーザー光11を反射偏向する回転多面鏡12(主走査周波数160Hzに相当)と、回転多面鏡12で反射偏向されたレーザー光Lをシート1上に収束し、かつ等速度で主走査させる走査レンズ21が配されている。
【0067】
さらに、前記レーザー光11が走査されるシート1の直上には、そのレーザー光11による励起で、シート1の上面から発せられる、蓄積記録されている画像情報に応じた輝尽発光光13aを上方より集光する光ガイド14aが近接して配置されている。光ガイド14aには集光した輝尽発光光13aを光電的に検出してアナログ画像信号QAに変換するフォトマルチプライヤ(光電子増倍管、以下、PMTという)15aが接続されている。このフォトマルチプライヤ15aは、電圧印加手段39aにより印加される電圧に応じた感度で輝尽発光光13aを検出するものである。
【0068】
また、フォトマルチプライヤ15aには対数増幅器16aが接続されており、フォトマルチプライヤ15aにより検出されたアナログ画像信号QAを予め設定されている周波数特性にしたがって対数的に増幅し、対数化画像信号QA′を出力する。
【0069】
一方、メモリー41aには、予め設定されているサンプリング間隔に応じたシェーディング補正用データD1が記憶されており、メモリー41aにはこの補正用データD1を予め設定されている基準クロックでアナログ信号D1′に変換するD/A変換器42aが接続されており、D/A変換器42aには、対数増幅器16aから出力された対数化画像信号QA′にこの補正用アナログ信号D1′を加算して、シェーディング補正が施された画像信号Q1を出力する加算器43aが接続されている。
【0070】
さらに加算器43aには、後述するA/D変換によるエリアジングノイズ(折り返しノイズ)を除去するアンチエリアジングフィルター35aが接続されており、その後段には、このフィルタリングされた後の画像信号Q1′を予め設定されている基準クロックでデジタル画像信号S1に変換するA/D変換器36aが設けられている。なお、アンチエリアジングフィルター35aは高密度用フィルターと低密度用フィルターからなり、初期的には高密度用フィルターが選択されているが、後述する特性変更手段60からの入力信号に応じて低密度用フィルターにも切り換えられるものである。
【0071】
また、A/D変換器36aには、デジタル画像信号S1の画素密度を変換して、画像信号S1′を得る画素密度変換手段37aが接続されている。
【0072】
さらにまた、本実施形態の画像情報読取装置は、画像情報の読取り画素密度として、予め初期的に設定されている「高密度」(10 pix/mm)と、この高密度よりも低密度である「低密度」(5 pix/mm)とのうち、オペレーターが選択した密度の入力を受け、「高密度」のときはパラメータ(m,n)をm=n=1とし、「低密度」のときはパラメータ(m,n)をm=n=2として出力する入力手段70と、この入力手段70から出力されたパラメータ(m,n)が入力されて、これらのパラメーター(m,n)に応じて、エンドレスベルト9aを駆動するモータ8の回転速度を変更し、走査レンズ21を光軸方向に移動させてシート1を走査するレーザー光11のビーム径を変更し、電圧印加手段39aを制御してフォトマルチプライヤ15aの感度を変更し、D/A変換器42aのサンプリングクロックを変更してメモリー41aに記憶されているシェーディング補正用データを変更して出力させ、A/D変換器36aのサンプリングクロックを変更し、対数増幅器16aの周波数特性を変更し、アンチエリアジングフィルター35aの高密度用と低密度用とを切り換え、画素密度変換手段37aにおいて画素密度変換を行う際のパラメータを変更する特性変更手段60とを備えている。
【0073】
ここで、画素密度変換手段37aにおいては、パラメータ(m,n)に応じて、デジタル画像信号S1の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換処理を行って、m=n=1の場合よりも、1/(m×n)2 倍の画素密度を有する画像信号S1′を得るものである。具体的には、デジタル画像信号S1の主副両走査方向に対して1次元マスク演算を施すことによる方法、画素密度に応じて画素を間引く方法、Bスプライン補間演算、Cubic スプライン補間演算等の高次の補間演算による方法(例えば特開平8-16767 号、同9-321981号)、線形補間演算による方法(例えば特開平9-50516 号)等により画素密度変換が行われる。また、この際、パラメータ(m,n)に応じて画素密度変換のパラメータが変更される。例えば1次元マスク演算を行う場合にはマスク係数が、間引き処理を行う場合には間引きの間隔が、補間演算を行う場合にはデジタル画像信号S1に施す補間演算の種類が変更される。
【0074】
次に本実施形態の放射線画像情報読取装置の作用について説明する。
【0075】
まず、オペレーターにより選択された「高密度」を表す信号が入力手段70に入力される。「高密度」はこの画像情報読取装置の初期的設定であり、「高密度」に対応するパラメータ(m,n)=(1,1)が入力手段70から出力されて特性変更手段60に入力される。
【0076】
特性変更手段60は、入力されたパラメータ(m,n)=(1,1)により、エンドレスベルト9aを駆動するモータ8の回転速度、走査レンズ21aの位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、A/D変換器36aのサンプリングクロック、対数増幅器16aの周波数特性、アンチエリアジングフィルター35aおよび画素密度変換手段37aを、それぞれ初期設定(高密度用)のものとなるようにセットする。
【0077】
次いで、エンドレスベルト9a上の所定の位置にセットされた蓄積性蛍光体シート1は、エンドレスベルト9aにより、矢印Y方向に上記高密度用の速度で搬送(副走査)される。
【0078】
一方、レーザ光源10から発せられたレーザー光11は、モータ20により駆動され矢印方向に高速回転する回転多面鏡12によって反射偏向され、この偏向されたレーザー光11は走査レンズ21により、エンドレスベルト9a上を搬送されるシート1で収束され、かつ等速度で走査されて、このシート1を矢印X方向に主走査する。
【0079】
シート1の、レーザー光11で走査された部分からは、そこに蓄積記録されている画像情報に応じた光量の輝尽発光光13aが発光し、輝尽発光光13aはレーザー光11の主走査線(X方向に平行)に沿って配された光ガイド14aの入射端面18aからこの光ガイド14a内に入射し、光ガイド14aの内部を全反射を繰り返しつつPMT15aまで導光される。なお、光ガイド14aとPMT15aとの接続部には、輝尽発光光13aとともに光ガイド14aに入射したレーザ光11がPMT15aに入射するのを阻止しつつ、輝尽発光光13aの入射を許容するように帯域設定されたレーザ光カットフィルター17aが設けられており、シート1の表面等で散乱し光ガイド14aに入射したレーザ光11は、このレーザ光カットフィルター17aによりカットされて、PMT15aに入射することはない。
【0080】
PMT15aは、電圧印加手段39aにより、高密度用の感度に対応した高圧が印加されており、入射した輝尽発光光13aを高密度用の感度で検出し、アナログ画像信号QAに光電変換して出力する。このアナログ画像信号QAは対数増幅器16aに入力され、高密度用の周波数特性にセットされた対数増幅器16aにより対数化画像信号QA′に変換されて加算器43aに入力される。
【0081】
一方、メモリー41aに記憶されている高密度用のシェーディング補正用データD1は、D/A変換器42aにより、高密度用のクロックである基準クロックでのサンプリングレートでアナログ信号D1′に変換される。
【0082】
ここで、D/A変換器42aに入力されるクロックは、図2に示すように、特性変更手段60に設けられたセレクター62により切り換えられて出力される基準クロックである。
【0083】
シェーディング補正用のアナログ信号D1′は加算器43aに入力され、対数増幅器16aから入力された対数化画像信号QA′に加算され、これにより対数化画像信号QA′はシェーディング補正がなされた画像信号Q1に変換されてアンチエリアジングフィルター35aに入力される。
【0084】
アンチエリアジングフィルター35aは、特性変更手段60により高密度用フィルターに切り換えられており、入力された画像信号Q1はこの高密度用フィルターによりエリアジングノイズが適切に除去されてA/D変換器36aに入力される。A/D変換器36aは、特性変更手段60のセレクター62により切り換えられた基準クロックによるサンプリングレートで、入力された画像信号Q1′をデジタル画像信号S1に変換して画素密度変換手段37aに出力する。
【0085】
画素密度変換手段37aは、デジタル画像信号S1の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換を行うものであるが、ここではm=n=1であるため、何ら画素密度変換処理を施すことなく、画像信号S1′を得、画像処理装置等に出力する。
【0086】
次に、オペレーターにより「低密度」が選択された場合の作用について説明する。
【0087】
オペレーターにより「低密度」を表す信号が入力手段70に入力されると、図3に示すように、入力手段70は「低密度」に対応するパラメータ(m,n)=(2,2)を特性変更手段60に出力する。特性変更手段60は、入力されたパラメータ(m,n)=(2,2)により、エンドレスベルト9aを駆動するモータの回転速度、走査レンズ21の位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、A/D変換器36aのサンプリングクロック、対数増幅器16aの周波数特性、アンチエリアジングフィルター35aおよび画素密度変換手段37aを、それぞれ低密度用のものとなるようにセットする。具体的には、エンドレスベルト9aを駆動するモータの回転速度を2倍の速度に変更し、走査レンズ21の位置を、レーザー光11のビーム径がシート1上で略2倍となる位置まで移動し、電圧印加手段39aへの制御信号を、フォトマルチプライヤ15aの感度を下げるような信号とし、セレクター62(図2)を分周回路61から出力されるクロック(基準クロックを分周して得られたクロックであって、基準クロックの2倍の周期のクロック)に切り換えて、D/A変換器42aおよびA/D変換器36aに入力されるサンプリングクロックを変更し、対数増幅器16aの周波数特性を変更し、アンチエリアジングフィルター35aを低密度用フィルターに切り換え、画素密度変換手段37aのパラメータを変更する。
【0088】
このように各部の特性が変更されて、前述した高密度読取りの場合と同様に画像情報の読取りがなされる。すなわち、エンドレスベルト9a上の所定の位置にセットされた蓄積性蛍光体シート1は、エンドレスベルト9aにより、矢印Y方向に上記高密度用の速度の2倍の速度で搬送(副走査)され、一方、レーザ光源10から発せられたレーザー光11は、高密度読取りのときと同一の回転速度で高速回転する回転多面鏡12によって反射偏向され、この偏向されたレーザー光11は走査レンズ21により、エンドレスベルト9a上を搬送されるシート1で収束され、かつ等速度で走査されて、このシート1を矢印X方向に主走査する。このときシート1上におけるビーム径は高密度読取りのときの2倍の大きさで収束される。
【0089】
シート1の、レーザー光11で走査された部分から発光した輝尽発光光13aは、光ガイド14aを介してPMT15aまで導光される。
【0090】
PMT15aは、電圧印加手段39aにより、高密度読取り用よりも低い感度となるような電圧が印加されており、入射した輝尽発光光13aをこの低密度用の感度で検出し、アナログ画像信号QAに光電変換して出力する。このアナログ画像信号QAは対数増幅器16aに入力され、低密度用の周波数特性にセットされた対数増幅器16aにより対数化画像信号QA′に変換されて加算器43aに入力される。
【0091】
一方、メモリー41aに記憶されている高密度用のシェーディング補正用データD1は、D/A変換器42aにより、低密度用のクロックである2分周のクロック(サンプリングレート)でアナログ信号D1′に変換される。すなわち、特性変更手段60に設けられた分周回路61が、基準クロックを2分周して基準クロックの2倍の周期のクロックを発生し、セレクター62はこの2分周のクロックを出力するように切り換えられている。この結果、シェーディング補正用のアナログ信号D1′は高密度用のシェーディング補正用データD1′に対して、主走査方向についてその画素数が1/2とされている。
【0092】
シェーディング補正用のアナログ信号D1′は加算器43aに入力され、対数増幅器16aから入力された対数化画像信号QA′に加算され、これにより対数化画像信号QA′はシェーディング補正がなされた画像信号Q1′に変換されてアンチエリアジングフィルター35aに入力される。
【0093】
アンチエリアジングフィルター35aは、特性変更手段60により低密度用フィルターに切り換えられており、入力された画像信号Q1′はこの低密度用フィルターによりエリアジングノイズが適切に除去されてA/D変換器36aに入力される。A/D変換器36aは、特性変更手段60のセレクター62により切り換えられた2分周のクロックによるサンプリングレートで、入力された画像信号Q1′をデジタル画像信号S1に変換する。
【0094】
画素密度変換手段37aは、デジタル画像信号S1の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換を行う。ここで、m=n=2であるため、デジタル画像信号S1の主走査方向および副走査方向についてそれぞれ画素数を1/2とする画素密度変換処理を施す。この画素密度変換処理として、ここでは1次元マスク演算について説明する。この1次元マスク演算に使用する1次元フィルタの例を下記に示す。
【0095】
a(x,1)=(-8/105,-5/105,34/105,63/105,34/105,-5/105,-8/105)
そしてこのフィルタa(x,1)を用いてデジタル画像信号S1の主走査方向について1画素間隔でフィルタリング処理を行った後、副走査方向について1画素間隔でフィルタリング処理を行うことにより画素密度変換が行われる。具体的には、デジタル画像信号S1の画素値をS1(x,y)、主走査方向における画素密度変換後の画素値をS1A(x/2,y)とすると、画素値S1A(x/2,y)は下記の式(1)により算出される。
【0096】
但し、k=1〜N(Nは主走査方向の画素数(画素位置又は画素番号))
l=1〜M(Mは副走査方向の画素数(画素位置又は画素番号))
同様にしき(1)の演算を副走査方向にも行なうことにより、画素密度変換がなされた画像信号S1′が得られる。
【0097】
一方、このような1次元フィルタを用いたマスク演算においては、デジタル画像信号S1の縁部分においては、マスク演算を行うためのデータが不足する。したがってそのような縁部分の画素に対してフィルタリング処理を施す場合は、その画素の外側に、仮想的な画素位置を定め、この仮想的な画素位置の画素にそれぞれ適切な画素値を与えたうえでフィルタリング処理を行えばよい。
【0098】
なお、画素密度変換手段37aにおける画素密度変換は、上記1次元フィルタを用いたマスク演算に限定されるものではなく、画素を間引く処理や補間演算により画素密度変換を行ってもよい。ここで補間演算としては、線形補間の他、滑らかさを重視したBスプライン補間演算、鮮鋭度を重視したCubic スプライン補間演算等の高次の補間演算を適用することもできる。
【0099】
そして、このように画素密度変換手段37aにおいて得られた画像信号S1′は、「低密度」が選択された場合には「高密度」が選択された場合と比較して1/16の画素密度の信号となり、画像処理装置等に出力される。
【0100】
このように本実施形態の画像情報読取装置によれば、レーザー光の主走査速度を変えることなく、高密度読取りの場合に比して、レーザー光Lの主走査方向について1/2、副走査方向について1/2の画素数である1/4の画素密度の低密度画像信号を得、さらにこの画像信号に対して主走査方向について1/2、副走査方向について1/2の画素数となるように画素密度変換処理を施すことで、全体として1/16の画素密度の低密度画像信号を得ることができる。しかも、特性変更手段60が、画素密度を規定するパラメータに応じて、走査レンズ21の位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、対数増幅器16aの周波数特性、アンチエリアジングフィルター35aおよび画素密度変換手段37aのパラメータをそれぞれ変更するため、読取り画素密度に適した、画素当たりのエネルギー付与、光電検出感度の設定、シェーディング補正、エリアジングノイズの抑制等が施された画像信号を得ることができる。
【0101】
なお、本実施形態の画像情報読取装置においては、高密度読取りが初期設定として設定されているものとしたが、本発明の画像情報読取方法および装置は、この態様に限られるものではなく、初期設定として「標準密度読取り」が設定されており、「高密度読取り」、「低密度読取り」にそれぞれ選択的に切り換えられるものとしてもよい。
【0102】
例えば、主走査周波数(主走査速度)を160Hz、画素密度10pix /mmを標準密度読取りとしたときの、主走査サンプリング間隔、副走査画素ピッチ、副走査速度、アンチエリアジングフィルターカットオフ周波数を以下に示す。
【0103】
主走査サンプリング間隔:100μm 副走査画素ピッチ:100μm
主走査サンプリング周期:1.0μsec 副走査速度:16mm/sec
アンチエリアジングフィルターカットオフ周波数:500kHz
(アナログフィルタのため、500kHz以下、例えば400kHzなどを適用するのが望ましい。)
これに対して、画素密度20pix /mmである高密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔:50μm 副走査画素ピッチ:50μm
主走査サンプリング周期:0.5μsec 副走査速度:8mm/sec
アンチエリアジングフィルターカットオフ周波数:1000kHz
となり、一方、画素密度5pix /mmである低密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔:200μm 副走査画素ピッチ:200μm
主走査サンプリング周期:2.0μsec 副走査速度:32mm/sec
アンチエリアジングフィルターカットオフ周波数:250kHz
となる。
【0104】
なお、上述したような「高密度読取り」、「低密度読取り」等、の読取り密度が段階的に予め設定されたものに限らず、読取り密度を規定するパラメータ(m,n)を任意に選択できるものとしてもよい。
【0105】
また特性変更手段60による、メモリ41aに記憶されているシェーディング補正用データの変更方法の詳細について、図4を参照して説明する。
【0106】
図4に示すように、メモリ41aが、選択されうるm,nに応じたシェーディング補正データD1がこれらのm,nに対応付けられてそれぞれ記憶されている大容量のハードディスク41cと、ハードディスク41cから読み出されたシェーディング補正データをシェーディング補正を行なうシェーディング回路(SHD回路)に受け渡すための一時的メモリであるシェーディングメモリ(SHDメモリ)41dとから構成されており、特性変更手段60は、これら2つのメモリ(ハードディスク41cおよびSHDメモリ41d)に対してメモリ制御を行なうものとすればよい。ここでメモリ制御としては例えば、m,nの値ごとに予め設定されたシェーディング補正データをハードディスク41cに予め記憶させておき、選択されたm,nの値に対応するシェーディング補正データを、選択の都度、ハードディスク41cからSHDメモリ41dに転送し、この転送されたシェーディング補正データをSHDメモリ41dから読み出してSHD回路に送るようにする方法や、m,nの値ごとに予め設定されたシェーディング補正データをハードディスク41cに予め記憶させておき、起動時など所望とするときに、これらのシェーディング補正データの全てを、ハードディスク41cから、m,nの値ごとに互いに異なるアドレスに対応した記憶領域を有するSHDメモリ41dにm,nの値に対応させて転送し、選択されたm,nの値に対応するシェーディング補正データを、選択されたm,nの値に対応するSHDメモリのアドレスから読み出してSHD回路に送るようにする方法などを適用することができる。
【0107】
ここでシェーディング補正の都度、特定のシェーディング補正データのみをSHDメモリ41dに転送する方法によれば、SHDメモリ41dの容量を少なくすることができ、ハードウェア構成を簡単にすることができる。一方、全てのシェーディング補正データをハードディスク41cからSHDメモリ41dに読み込む方法によれば、ソフトウェア構成を簡単にすることができ、SHDメモリ41dからの補正データの読出し速度を高速化することができる。
【0108】
図5は本発明の第2の画像情報読取装置の一実施形態を示す図である。図5に示すように、本実施形態の画像情報読取装置は、被写体の放射線画像が蓄積記録された蓄積性蛍光体シートの両面から放射線画像を表す画像信号を得るものである。蓄積性蛍光体シート1が、モータ8により回転せしめられるエンドレスベルト9a,9b上に配置される。このシート1の上方には、励起光としてのレーザ光11を発するレーザ光源10と、レーザ光11を反射偏向してシート1を主走査する、モータ20により回転される回転多面鏡12と、レーザ光11をシート1に結像するための走査レンズ21とが配されている。さらに、レーザ光11が走査される位置の上方には、そのレーザ光11の走査により発せられる輝尽発光光を上方より集光する集光ガイド14aがシート1に近接して配置され、その位置の下方には、輝尽発光光を下方より集光する集光ガイド14bがシート1と略垂直に配置されている。各集光ガイド14a,14bは、それぞれ輝尽発光光を光電的に検出するPMT15a,15bが接続されている。このPMT15a,15bは対数増幅器16a,16bに接続されており、PMT15a,15bにより検出されたアナログ画像信号QA,QBを予め設定されている周波数特性にしたがって対数的に増幅し、対数化画像信号QA′,QB′を出力する。
【0109】
一方、メモリ41a,41bには、予め設定されているサンプリング間隔に応じたシェーディング補正用データD1,D2が記憶されており、メモリ41a,41bにはこの補正用データD1,D2を予め設定されている基準クロックでアナログ信号D1′,D2′に変換するD/A変換器42a,42bが接続されており、D/A変換器42a,42bには、対数増幅器16a,16bから出力された対数化画像信号QA′,QB′にこの補正用アナログ信号D1′、D2′を加算して、シェーディング補正が施された画像信号Q1,Q2を出力する加算器43a,43bが接続されている。
【0110】
さらに加算器43a,43bには、後述するA/D変換によるエリアジングノイズ(折り返しノイズ)を除去するアンチエリアジングフィルタ35a,35bが接続されており、その後段には、このフィルタリングされた後の画像信号Q1′,Q2′を予め設定されている基準クロックでデジタル画像信号S1,S2に変換するA/D変換器36a,36bが設けられている。なお、アンチエリアジングフィルタ35a,35bは高密度用フィルタと低密度用フィルタとからなり、初期的には高密度用フィルタが選択されているが、後述する特性変更手段60からの入力信号に応じて低密度用フィルタにも切り換えられるものである。
【0111】
さらにまた、A/D変換器36a,36bには、デジタル画像信号S1,S2の画素密度を変換して、画像信号S1′,S2′を得る画素密度変換手段37a,37bが接続されている。そして、画素密度変換された画像信号S1′,S2′は加算手段38において加算されて加算画像信号S3が得られる。
【0112】
また、本実施形態の画像情報読取装置は、図6に示すように、画像情報の読取り画素密度として、予め初期的に設定されている「高密度」(10pix/mm)と、この高密度よりも低密度である「低密度」(5pix/mm)とのうち、オペレーターが選択した密度の入力を受け、「高密度」のときはパラメータ(m,n)をm=n=1とし、「低密度」のときはパラメータ(m,n)をm=n=2として出力する入力手段70と、この入力手段70から出力されたパラメータ(m,n)が入力されて、これらのパラメータ(m,n)に応じて、エンドレスベルト9a,9bを駆動するモータ61の回転速度を変更し、レーザ光源10の出力を変更してレーザ光11のパワーを変更し、走査レンズ21を光軸方向に移動させてシート1を走査するレーザ光11のビーム径を変更し、電圧印加手段39a,39bを制御してPMT15a,15bの感度を変更し、D/A変換器42a,42bのサンプリングクロックを変更してメモリ41a,41bに記憶されているシェーディング補正用データD1,D2を変更して出力させ、D/A変換器42a,42bがシェーディング補正用データをメモリ41a,41bから読み出すタイミングを変更し、A/D変換器36a,36bのサンプリングクロックを変更し、対数増幅器16a,16bの周波数特性を変更し、アンチエリアジングフィルタ35a,35bの高密度用と低密度用とを切り換え、画素密度変換手段37a,37bにおいて画素密度変換を行う際のパラメータを変更する特性変更手段60とを備えている。
【0113】
ここで、画素密度変換手段37a,37bにおいては、パラメータ(m,n)に応じて、デジタル画像信号S1,S2の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換処理を行って、m=n=1の場合よりも、1/(m×n)2 倍の画素密度を有する画像信号S1′,S2′を得るものである。具体的には、デジタル画像信号S1,S2の主副両走査方向に対して1次元マスク演算を施すことによる方法、画素密度に応じて画素を間引く方法、Bスプライン補間演算、Cubic スプライン補間演算等の高次の補間演算による方法、線形補間演算による方法等により画素密度変換が行われる。また、この際、パラメータ(m,n)に応じて画素密度変換のパラメータが変更される。例えば1次元マスク演算を行う場合にはマスク係数が、間引き処理を行う場合には間引きの間隔が、補間演算を行う場合にはデジタル画像信号S1,S2に施す補間演算の種類が変更される。
【0114】
なお、加算手段38においては、上記特開平7-287330号に記載されたように、加算画像信号S3のS/Nを高くするような周波数応答特性を有するフィルタによるフィルタリング処理を画像信号S1′,S2′に対して施した後に、加算画像信号S3を得る方法を採用することが好ましいが、単純に画像信号S1′,S2′を加算するものであってもよい。なお、フィルタリング処理を行う場合には、画素密度変換手段37a,37bにおいて行うようにしてもよい。
【0115】
次に本実施形態の放射線画像情報読取装置の作用について説明する。
【0116】
まず、オペレーターにより選択された「高密度」を表す信号が入力手段70に入力される。「高密度」はこの画像情報読取装置の初期的設定であり、「高密度」に対応するパラメータ(m,n)=(1,1)が入力手段70から出力されて特性変更手段60に入力される。
【0117】
特性変更手段60は、入力されたパラメータ(m,n)=(1,1)により、エンドレスベルト9a,9bを駆動するモータ8の回転速度、レーザ光源10の出力、走査レンズ21の位置、電圧印加手段39a,39bへの制御信号、D/A変換器42a,42bのサンプリングクロック、D/A変換器42a,42bがシェーディング補正用データD1,D2をメモリ41a,41bから読み出すタイミング、A/D変換器36a,36bのサンプリングクロック、対数増幅器16a,16bの周波数特性、アンチエリアジングフィルタ35a,35b、および画素密度変換手段37a,37bを、それぞれ初期設定(高密度用)のものとなるようにセットする。
【0118】
被写体の放射線画像が蓄積記録された蓄積性蛍光体シート1がエンドレスベルト9a,9bに配置されると、エンドレスベルト9a,9bにより矢印Y方向に搬送(副走査)される。一方、レーザ光源10から発せられたレーザ光11はモータ20により駆動され矢印方向に高速回転する回転多面鏡12によって反射偏向され、走査レンズ21を通ってシート1に入射し副走査の方向(矢印Y方向)と略垂直な矢印X方向に主走査する。このレーザ光11が照射されたシート1の箇所からは、蓄積記録されている放射線画像情報に応じた光量の輝尽発光光13a,13b(ここで、輝尽発光光13a,13bはそれぞれシート1の上方(表面)、下方(裏面)から発散されたものを示す)が発散される。この輝尽発光光13aは集光ガイド14aによって導かれ、PMT15aによって光電的に検出される。入射端面18aから集光ガイド14a内に入射した輝尽発光光13aは、集光ガイド14aの内部を全反射を繰り返して進み、出射端面から出射してPMT15aに受光され、放射線画像を表す輝尽発光光13aの光量がPMT15aによってアナログ画像信号QAに変換される。同様に、輝尽発光光13bは集光ガイド14bによって導かれ、PMT15bによって光電的に検出されてアナログ画像信号QBに変換される。
【0119】
これらのアナログ画像信号QA,QBは対数増幅器16a,16bに入力され、高密度用の周波数特性にセットされた対数増幅器34により対数化画像信号QA′,QB′に変換されて加算器43a,43bに入力される。
【0120】
一方、メモリ41a,41bに記憶されている高密度用のシェーディング補正用データD1,D2は、D/A変換器42a,42bにより、高密度用のタイミングで読み出されるとともに、高密度用のクロックである基準クロックでのサンプリングレートでアナログ信号D1′,D2′に変換される。
【0121】
ここで、D/A変換器42a,42bに入力されるクロックは、図7に示すように、特性変更手段60に設けられたセレクター62により切り換えられて出力される基準クロックである。
【0122】
シェーディング補正用のアナログ信号D1′,D2′は加算器43a,43bに入力され、対数増幅器16a,16bから入力された対数化画像信号QA′,QB′に加算され、これにより対数化画像信号QA′,QB′はシェーディング補正がなされた画像信号Q1,Q2に変換されてアンチエリアジングフィルタ35a,35bに入力される。
【0123】
アンチエリアジングフィルタ35a,35bは、特性変更手段60により高密度用フィルタに切り換えられており、入力された画像信号Q1,Q2はこの高密度用フィルタによりエリアジングノイズが適切に除去されてA/D変換器36a,36bに入力される。
【0124】
A/D変換器36a,36bは、特性変更手段60のセレクター62により切り換えられた基準クロックによるサンプリングレートで、入力された画像信号Q1′,Q2′をそれぞれデジタル画像信号S1,S2に変換する。
【0125】
画素密度変換手段37a,37bは、デジタル画像信号S1,S2の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換を行うものであるが、ここではm=n=1であるため、何ら画素密度変換処理を施すことなく、画像信号S1′,S2′を得る。
【0126】
加算手段38は、画像信号S1′,S2′に対して、上記特開平7-287330号に記載されたような、加算画像信号S3のS/Nを高くするような周波数応答特性を有するフィルタによるフィルタリング処理を施し、このフィルタリング処理を施した画像信号S1′,S2′を相対応する画素同士で加算して加算画像信号S3を得て、画像処理装置等に出力する。
【0127】
次に、オペレーターにより「低密度」が選択された場合の作用について説明する。
【0128】
オペレーターにより「低密度」を表す信号が入力手段70に入力されると、入力手段70は「低密度」に対応するパラメータ(m,n)=(2,2)を特性変更手段60に出力する。
【0129】
特性変更手段60は、入力されたパラメータ(m,n)=(2,2)により、エンドレスベルト9a,9bを駆動するモータ8の回転速度、レーザ光源10の出力、走査レンズ21の位置、電圧印加手段39a,39bへの制御信号、D/A変換器42a,42bのサンプリングクロック、D/A変換器42a,42bがシェーディング補正用データD1,D2をメモリ41a,41bから読み出すタイミング、A/D変換器36a,36bのサンプリングクロック、対数増幅器16a,16bの周波数特性、アンチエリアジングフィルタ35a,35b、および画素密度変換手段37a,37bを、それぞれ低密度用のものとなるようにセットする。
【0130】
具体的には、エンドレスベルト9a,9bを駆動するモータ8の回転速度を2倍の速度に変更し、レーザ光源10の出力を略2倍とし、走査レンズ21の位置をレーザ光11のビーム径がシート1上で略2倍となる位置まで移動し、電圧印加手段39a,39bへの制御信号をフォトマルチプライヤ15a,15bの感度を下げるような信号とし、D/A変換器42a,42bがシェーディング補正用データD1,D2をメモリ41a,41bから読み出すタイミングを変更し、セレクター62(図7参照)を分周回路61から出力されるクロック(基準クロックを分周して得られたクロックであって、基準クロックの2倍の周期のクロック)に切り換えて、D/A変換器42a,42bおよびA/D変換器36a,36bに入力されるサンプリングクロックを変更し、対数増幅器16a,16bの周波数特性を変更し、アンチエリアジングフィルタ35a,35bを低密度用フィルタに切り換え、画素密度変換手段37a,37bのパラメータを変更する。
【0131】
このように各部の特性が変更されて、前述した高密度読取りの場合と同様に画像情報の読取りがなされる。
【0132】
すなわち、エンドレスベルト9a,9b上の所定の位置にセットされた蓄積性蛍光体シート1は、エンドレスベルト9a,9bにより、矢印Y方向に上記高密度用の速度の2倍の速度で搬送(副走査)され、一方、レーザ光源10から発せられた高密度用の2倍のパワーのレーザ光11は、高密度読取りのときと同一の回転速度で高速回転する回転多面鏡12によって反射偏向され、この偏向されたレーザ光11は走査レンズ21により、エンドレスベルト9a,9b上を搬送されるシート1で収束され、かつ等速度で走査されて、このシート1を矢印X方向に主走査する。このときシート1上におけるビーム径は高密度読取りのときの2倍の大きさで収束される。
【0133】
シート1の、レーザ光11で走査された部分から発光した輝尽発光光13a,13bは、集光ガイド14a,14bを介してPMT15a,15bまで導光される。PMT15a,15bは、電圧印加手段39a,39bにより、高密度読取り用よりも低い感度となるような電圧が印加されており、入射した輝尽発光光13a,13bをこの低密度用の感度で検出し、アナログ画像信号QA,QBに光電変換して出力する。このアナログ画像信号QA,QBは対数増幅器16a,16bに入力され、低密度用の周波数特性にセットされた対数増幅器16a,16bにより対数化画像信号QA′,QB′に変換されて加算器43a,43bに入力される。
【0134】
一方、メモリ41a,41bに記憶されている高密度用のシェーディング補正用データD1,D2は、D/A変換器42a,42bにより、低密度用のタイミングで読み出されるとともに、低密度用のクロックである2分周のクロックでのサンプリングレートでアナログ信号D1′,D2′に変換される。すなわち、特性変更手段60に設けられた分周回路61が、基準クロックを2分周して基準クロックの2倍の周期のクロックを発生し、セレクター62はこの2分周のクロックを出力するように切り換えられている。この結果、シェーディング補正用のアナログ信号D1′,D2′は高密度用のシェーディング補正用データD1,D2に対して、主走査方向についてその画素数が1/2とされている。
【0135】
シェーディング補正用のアナログ信号D1′,D2′は加算器43a,43bに入力され、対数増幅器16a,16bから入力された対数化画像信号QA′,QB′に加算され、これにより対数化画像信号QA′,QB′はシェーディング補正がなされた画像信号Q1,Q2に変換されてアンチエリアジングフィルタ35a,35bに入力される。
【0136】
アンチエリアジングフィルタ35a,35bは、特性変更手段60により低密度用フィルタに切り換えられており、入力された画像信号Q1,Q2はこの低密度用フィルタによりエリアジングノイズが適切に除去されてA/D変換器36a,36bに入力される。
【0137】
A/D変換器36a,36bは、特性変更手段60のセレクター62により切り換えられた2分周のクロックによるサンプリングレートで、入力された画像信号QA′,QB′をデジタル画像信号S1,S2に変換する。
【0138】
画素密度変換手段37a,37bは、デジタル画像信号S1,S2の主走査方向について1/m倍、副走査方向について1/n倍とする画素密度変換を行う。ここで、m=n=2であるため、デジタル画像信号S1,S2の主走査方向および副走査方向についてそれぞれ画素数を1/2とする画素密度変換処理を施す。この画素密度変換処理として、ここでは1次元マスク演算について説明する。この1次元マスク演算に使用する1次元フィルタの例を下記に示す。
【0139】
a(x,1)=(-8/105,-5/105,34/105,63/105,34/105,-5/105,-8/105)
そしてこのフィルタa(x,1)を用いてデジタル画像信号S1,S2の主走査方向について1画素間隔でフィルタリング処理を行った後、副走査方向について1画素間隔でフィルタリング処理を行うことにより画素密度変換が行われる。具体的には、デジタル画像信号S1(S2についても同様であるため、ここではS1についてのみ説明する)の画素値をS1(x,y)、主走査方向における画素密度変換後の画素値をS1A(x/2,y)とすると、画素値S1A(x/2,y)は下記の式(1)により算出される。
【0140】
但し、k=1〜N(Nは主走査方向の画素数(画素位置又は画素番号))
l=1〜M(Mは副走査方向の画素数(画素位置又は画素番号))
同様に式(1)の演算を副走査方向にも行うことにより、画素密度変換がなされた画像信号S1′,S2′が得られる。
【0141】
一方、このような1次元フィルタを用いたマスク演算においては、デジタル画像信号S1,S2の縁部分においては、マスク演算を行うためのデータが不足する。例えば、図8に示す画像信号において、S1(1,1)およびS1(N,1)(斜線部)の画素に対してフィルタリング処理を施す場合は、3画素分のデータが不足する。この場合、仮想的な画素位置S1(−1,1),S1(−2,1),S1(−3,1)およびS1(N+1,1),S1(N+2,1),S1(N+3,1)を定め、この仮想的な画素位置に画素位置S1(1,1)およびS1(N,1)のデータ値をコピーして、仮想的な画素位置においてデータ値が存在するものとしてフィルタリング処理を行えばよい。
【0142】
なお、画素密度変換手段37a,37bにおける画素密度変換は、上記1次元フィルタを用いたマスク演算に限定されるものではなく、画素を間引く処理や補間演算により画素密度変換を行ってもよい。
【0143】
ここで補間演算としては、線形補間の他、滑らかさを重視したBスプライン補間演算、鮮鋭度を重視したCubic スプライン補間演算等の高次の補間演算を適用することができる。
【0144】
ここで、Cubic スプライン補間演算およびBスプライン補間演算について説明する。本実施形態において使用される画像信号S1,S2は、等間隔の周期でサンプリングされた一方向に配列されたサンプリング点(画素)Xk-2 ,Xk-1 ,Xk ,Xk+1 ,Xk+2 ,…にそれぞれ対応した信号値(Sk-2 ,Sk-1 ,Sk ,Sk+1 ,Sk+2 ,…)を有するものする。Cubic スプライン補間演算は、オリジナルのサンプリング点(画素)Xk 〜Xk+1 間に設けられた補間点Xp の補間データY′を表す3次のCubic スプライン補間演算式(2)における補間データYk-1 ,Yk ,Yk+1 ,Yk+2 にそれぞれ対応する補間係数ck-1 ,ck ,ck+1 ,ck+2 を、下記にそれぞれ示す演算により求めるものである。
【0145】
Y′=ck-1 Yk-1 +ck Yk +ck+1 Yk+1 +ck+2 Yk+2 (2)
ck-1 =(−t3 +2t2 −t)/2
ck =(3t3 −5t2 +2)/2
ck+1 =(−3t3 +4t2 +t)/2
ck+2 =(t3 −t2 )/2
(但し、t(0≦t≦1)は格子間隔を1とし、画素Xk を基準としたときの補間点Xp の画素Xk+1 方向への位置を示す。)
Bスプライン補間演算は、オリジナルのサンプリング点Xk 〜Xk+1 間に設けられた補間点Xp の補間データY′を表す3次のBスプライン補間演算式(3)における補間データYk-1 ,Yk ,Yk+1 ,Yk+2 にそれぞれ対応する補間係数bk-1 ,bk ,bk+1 ,bk+2 を、下記にそれぞれ示す演算により求めるものである。
【0146】
Y′=bk-1 Yk-1 +bk Yk +bk+1 Yk+1 +bk+2 Yk+2 (3)
bk-1 =(−t3 +3t2 −3t+1)/6
bk =(3t3 −6t2 +4)/6
bk+1 =(−3t3 +3t2 +3t+1)/6
bk+2 =t3 /6
(但し、t(0≦t≦1)は格子間隔を1とし、画素Xk を基準としたときの補間点Xp の画素Xk+1 方向への位置を示す。)
本実施形態においては、m,nの値に応じて線形補間演算も含めて、補間演算の種類を選択すればよい。
【0147】
そして、このように画素密度変換手段37a,37bにおいて得られた画像信号S1′,S2′は加算手段38に入力される。ここで、「低密度」が選択された場合には「高密度」が選択された場合と比較して、画素密度は1/16となる。
【0148】
加算手段38は、画像信号S1′,S2′に対して、上記特開平7-287330号に記載されたような、加算画像信号S3のS/Nを高くするような周波数応答特性を有するフィルタによるフィルタリング処理を施し、このフィルタリング処理を施した画像信号S1′,S2′を相対応する画素同士で加算して加算画像信号S3を得て、画像処理装置等に出力する。
【0149】
このように本実施形態の画像情報読取装置によれば、レーザ光11の主走査速度を変えることなく、高密度読取りの場合に比して、レーザ光11の主走査方向について1/2、副走査方向について1/2の画素数である、1/4の画素密度の画像信号を得、さらにこの画像信号に対して主走査方向について1/2、副走査方向について1/2の画素数となるように画素密度変換処理を施して、1/16の画素密度の低密度画像信号を得ることができる。しかも、特性変更手段60が、画素密度を規定するパラメータに応じて、走査レンズ21の位置、電圧印加手段39a,39bへの制御信号、D/A変換器42a,42bのサンプリングクロック、対数増幅器16a,16bの周波数特性、アンチエリアジングフィルタ35a,35b、および画素密度変換手段37a,37bのパラメータをそれぞれ変更するため、読取り画素密度に適した、画素当たりのエネルギー付与、光電検出感度の設定、シェーディング補正、エリアジングノイズの抑制等が施された画像信号を得ることができる。
【0150】
とくに、本実施形態のように両面読取りを行う場合には、レーザ光11によるエネルギーをシート1の裏面側まで十分に付与するために、片面側からのみ読取る方式に比して走査速度を遅くする必要があるが、その場合においても画素密度を粗くすることが許容されるときは副走査速度を速めることができ、これにより1つのシート当たりの走査時間を短縮することができる。
【0151】
また、2つのデジタル画像信号S1,S2に対して画素密度変換処理を施すことにより、加算画像信号に対して画素密度変換を施す場合と比較して、加算画像信号を得る際の加算マスク処理を行うための演算量を低減することができ、これにより加算演算の演算時間を短縮して処理を高速に行うことができる。
【0152】
なお、本実施形態の画像情報読取装置においては、高密度読取りが初期設定として設定されているものとしたが、本発明の画像情報読取方法および装置は、この態様に限られるものではなく、初期設定として「標準密度読取り」が、設定されており、「高密度読取り」、「低密度読取り」に切り換えられるものとしてもよい。
【0153】
例えば、主走査周波数(主走査速度)を160Hz、画素密度10pix /mmを標準密度読取りとしたときの、主走査サンプリング間隔、副走査画素ピッチ、副走査速度、アンチエリアジングフィルタカットオフ周波数を以下に示す。
【0154】
主走査サンプリング間隔:100μm 副走査画素ピッチ:100μm
主走査サンプリング周期:1.0μsec 副走査速度:16mm/sec
アンチエリアジングフィルタカットオフ周波数:500kHz
(アナログフィルタのため、カットオフ周波数は500kHz以下であることが好ましく、例えば400kHzなどを適用するのが望ましい。)
これに対して、画素密度20pix /mmである高密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔:50μm 副走査画素ピッチ:50μm
主走査サンプリング周期:0.5μsec 副走査速度:8mm/sec
アンチエリアジングフィルタカットオフ周波数:1000kHz
となり、一方、画素密度5pix /mmである低密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔:200μm 副走査画素ピッチ:200μm
主走査サンプリング周期:2.0μsec 副走査速度:32mm/sec
アンチエリアジングフィルタカットオフ周波数:250kHz
となる。
【0155】
なお、上述したような「高密度読取り」、「低密度読取り」等、の読取り密度が段階的に予め設定されたものに限らず、読取り密度を規定するパラメータ(m,n)を任意に選択できるものとしてもよい。
【0156】
また、上記実施形態においては、デジタル画像信号S1,S2に対して画素密度変換処理を施してるが、デジタル画像信号S1,S2を先に加算して加算画像信号S3を得、この加算画像信号S3に対して画素密度変換処理を施してもよい。
【0157】
さらに、上記実施形態においては、画素密度変換手段37a,37bにおいて、デジタル画像信号S1,S2の主走査方向について1/m倍、副走査方向について1/n倍となるように画素密度変換をしているが、主走査方向についてa/m(a>0)倍、副走査方向についてa/n倍となるように画素密度変換を行うようにしてもよい。この場合、加算画像信号S3の画素密度は(a/(m×n))2 倍となる。
【0158】
図9は本発明の第3の画像情報読取装置の一実施形態を示す図である。図示の画像情報読取装置は、モータ8により回転せしめられるエンドレスベルト9a上に放射線画像情報が蓄積記録された蓄積性蛍光体シート1が配置され、シート1の上方には、シート1を励起するレーザー光11を発するレーザー光源10と、モーター12により回転され、レーザー光11を反射偏向する回転多面鏡20(主走査周波数160Hzに相当)と、回転多面鏡20で反射偏向されたレーザー光11をシート1上に収束し、かつ等速度で主走査させる走査レンズ21が配されている。
【0159】
さらに、前記レーザー光11が走査されるシート1の直上には、そのレーザー光11による励起で、シート1の上面から発せられる、蓄積記録されている画像情報に応じた輝尽発光光13aを上方より集光する光ガイド14aが近接して配置されている。光ガイド14aには集光した輝尽発光光13aを光電的に検出してアナログ画像信号QAに変換するフォトマルチプライヤ(光電子増倍管)15aが接続されている。このフォトマルチプライヤ15aは、電圧印加手段39aにより印加される電圧に応じた感度で輝尽発光光13aを検出するものである。
【0160】
また、フォトマルチプライヤ15aには対数増幅器16aが接続されており、フォトマルチプライヤ15aにより検出されたアナログ画像信号QAを予め設定されている周波数特性にしたがって対数的に増幅し、対数化画像信号QA′を出力する。
【0161】
一方、メモリー41aには、予め設定されているサンプリング間隔に応じたシェーディング補正用データD1が記憶されており、メモリー41aにはこの補正用データD1を予め設定されている基準クロックでアナログ信号D1′に変換するD/A変換器41aが接続されており、D/A変換器42aには、対数増幅器16aから出力された対数化画像信号QA′にこの補正用アナログ信号D1′を加算して、シェーディング補正が施された画像信号Q1を出力する加算器43aが接続されている。
【0162】
さらに加算器43aには、後述するA/D変換によるエリアジングノイズ(折り返しノイズ)を除去するアンチエリアジングフィルター35aが接続されており、その後段には、このフィルタリングされた後の画像信号Q1′を予め設定されている基準クロックでデジタル画像信号S1に変換するA/D変換器36aが設けられている。なお、アンチエリアジングフィルター35aは高密度用フィルターと低密度用フィルターとからなり、初期的には高密度用フィルターが選択されているが、後述する特性変更手段60からの入力信号に応じて低密度用フィルターにも切り換えられるものである。
【0163】
さらにまた、本実施形態の画像情報読取装置は、画像情報の読取り画素密度として、予め初期的に設定されている「高密度」(10 pix/mm)と、この高密度よりも低密度である「低密度」(5 pix/mm)とのうち、オペレーターが選択した密度の入力を受け、「高密度」のときはパラメータ(m,n)をm=n=1とし、「低密度」のときはパラメータ(m,n)をm=n=2として出力する入力手段70と、この入力手段70から出力されたパラメータ(m,n)が入力されて、これらのパラメーター(m,n)に応じて、エンドレスベルト9aを駆動するモータ8の回転速度を変更し、走査レンズ21を光軸方向に移動させてシート1を走査するレーザー光11のビーム径を変更し、電圧印加手段39aを制御してフォトマルチプライヤ15aの感度を変更し、D/A変換器42aのサンプリングクロックを変更してメモリー41aに記憶されているシェーディング補正用データを変更して出力させ、A/D変換器36aのサンプリングクロックを変更し、対数増幅器16aの周波数特性を変更し、アンチエリアジングフィルター35aの高密度用と低密度用とを切り換える特性変更手段60とを備えている。
【0164】
次に本実施形態の放射線画像情報読取装置の作用について説明する。
【0165】
まず、オペレーターにより選択された「高密度」を表す信号が入力手段70に入力される。「高密度」はこの画像情報読取装置の初期的設定であり、「高密度」に対応するパラメータ(m,n)=(1,1)が入力手段70から出力されて特性変更手段60に入力される。
【0166】
特性変更手段70は、入力されたパラメータ(m,n)=(1,1)により、エンドレスベルト9aを駆動するモータ8の回転速度、走査レンズ21の位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、A/D変換器36のサンプリングクロック、対数増幅器16aの周波数特性、およびアンチエリアジングフィルター35aを、それぞれ初期設定(高密度用)のものとなるようにセットする。
【0167】
次いで、エンドレスベルト9a上の所定の位置にセットされた蓄積性蛍光体シート1は、エンドレスベルト9aにより、矢印Y方向に上記高密度用の速度で搬送(副走査)される。
【0168】
一方、レーザ光源10から発せられたレーザー光11は、モータ12により駆動され矢印方向に高速回転する回転多面鏡20によって反射偏向され、この偏向されたレーザー光11は走査レンズ21により、エンドレスベルト9a上を搬送されるシート1で収束され、かつ等速度で走査されて、このシート1を矢印X方向に主走査する。
【0169】
シート1の、レーザー光11で走査された部分からは、そこに蓄積記録されている画像情報に応じた光量の輝尽発光光13aが発光し、輝尽発光光13aはレーザー光11の主走査線に沿って配された光ガイド14aの入射端面18aからこの光ガイド14a内に入射し、光ガイド14aの内部を全反射を繰り返しつつ、レーザ光カットフィルタ17aを介してフォトマルチプライヤ15aまで導光される。
【0170】
フォトマルチプライヤ15aは、電圧印加手段39aにより、高密度用の感度に対応した高圧が印加されており、入射した輝尽発光光13aを高密度用の感度で検出し、アナログ画像信号QAに光電変換して出力する。
【0171】
このアナログ画像信号QAは対数増幅器16aに入力され、高密度用の周波数特性にセットされた対数増幅器16aにより対数化画像信号QA′に変換されて加算器43aに入力される。
【0172】
一方、メモリー41aに記憶されている高密度用のシェーディング補正用データD1は、D/A変換器42aにより、高密度用のクロックである基準クロックでのサンプリングレートでアナログ信号D1′に変換される。
【0173】
ここで、D/A変換器42aに入力されるクロックは、図2に示すように、特性変更手段60に設けられたセレクター62により切り換えられて出力される基準クロックである。
【0174】
シェーディング補正用のアナログ信号D1′は加算器43aに入力され、対数増幅器16aから入力された対数化画像信号QA′に加算され、これにより対数化画像信号QA′はシェーディング補正がなされた画像信号Q1に変換されてアンチエリアジングフィルター35aに入力される。
【0175】
アンチエリアジングフィルター35aは、特性変更手段60により高密度用フィルターに切り換えられており、入力された画像信号Q1はこの高密度用フィルターによりエリアジングノイズが適切に除去された信号Q1′に変換されてA/D変換器36aに入力される。
【0176】
A/D変換器36aは、特性変更手段60のセレクター62により切り換えられた基準クロックによるサンプリングレートで、入力された画像信号Q1′をデジタル画像信号S1に変換して、画像処理装置等に出力する。
【0177】
次に、オペレーターにより「低密度」が選択された場合の作用について説明する。
【0178】
オペレーターにより「低密度」を表す信号が入力手段70に入力されると、入力手段70は「低密度」に対応するパラメータ(m,n)=(2,2)を特性変更手段60に出力する。
【0179】
特性変更手段70は、入力されたパラメータ(m,n)=(2,2)により、エンドレスベルト9aを駆動するモータの回転速度、走査レンズ21の位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、A/D変換器36aのサンプリングクロック、対数増幅器16aの周波数特性、およびアンチエリアジングフィルター35aを、それぞれ低密度用のものとなるようにセットする。具体的には、エンドレスベルト9aを駆動するモータの回転速度を2倍の速度に変更し、走査レンズ21の位置を、レーザー光11のビーム径がシート1上で略2倍となる位置まで移動し、電圧印加手段39aへの制御信号を、フォトマルチプライヤ15aの感度を下げるような信号とし、セレクター62(図2参照)を分周回路61から出力されるクロック(基準クロックを分周し手得られたクロックであって、基準クロックの2倍の周期のクロック)に切り換えて、D/A変換器42aおよびA/D変換器36aに入力されるサンプリングクロックを変更し、対数増幅器16aの周波数特性を変更し、アンチエリアジングフィルター35aを低密度用フィルターに切り換える。
【0180】
このように各部の特性が変更されて、前述した高密度読取りの場合と同様に画像情報の読取りがなされる。
【0181】
すなわち、エンドレスベルト9a上の所定の位置にセットされた蓄積性蛍光体シート1は、エンドレスベルト9aにより、矢印Y方向に上記高密度用の速度の2倍の速度で搬送(副走査)され、一方、レーザ光源10から発せられたレーザー光11は、高密度読取りのときと同一の回転速度で高速回転する回転多面鏡20によって反射偏向され、この偏向されたレーザー光11は走査レンズ21により、エンドレスベルト9a上を搬送されるシート1で収束され、かつ等速度で走査されて、このシート1を矢印X方向に主走査する。このときシート1上におけるビーム径は高密度読取りのときの2倍の大きさで収束される。
【0182】
シート1の、レーザー光11で走査された部分から発光した輝尽発光光13aは、光ガイド14aを通じてレーザ光カットフィルタ17aを介し、フォトマルチプライヤ15aまで導光される。
【0183】
フォトマルチプライヤ15aは、電圧印加手段39aにより、高密度読取り用よりも低い感度となるような電圧が印加されており、入射した輝尽発光光13aをこの低密度用の感度で検出し、アナログ画像信号QAに光電変換して出力する。
【0184】
このアナログ画像信号QAは対数増幅器16aに入力され、低密度用の周波数特性にセットされた対数増幅器16aにより対数化画像信号QA′に変換されて加算器43aに入力される。
【0185】
一方、メモリー41aに記憶されている高密度用のシェーディング補正用データD1は、D/A変換器42aにより、低密度用のクロックである2分周のクロックでのサンプリングレートでアナログ信号D1′に変換される。すなわち、特性変更手段60に設けられた分周回路61が、基準クロックを2分周して基準クロックの2倍の周期のクロックを発生し、セレクター62はこの2分周のクロックを出力するように切り換えられている。この結果、シェーディング補正用のアナログ信号D1′は高密度用のシェーディング補正用データQA′に対して、主走査方向についてその画素数が1/2とされている。
【0186】
シェーディング補正用のアナログ信号D1′は加算器43aに入力され、対数増幅器16aから入力された対数化画像信号QA′に加算され、これにより対数化画像信号QA′はシェーディング補正がなされた画像信号Q1に変換されてアンチエリアジングフィルター35aに入力される。
【0187】
アンチエリアジングフィルター35aは、特性変更手段60により低密度用フィルターに切り換えられており、入力された画像信号Q1はこの低密度用フィルターによりエリアジングノイズが適切に除去された信号Q1′に変換されてA/D変換器36aに入力される。
【0188】
A/D変換器36aは、特性変更手段60のセレクター62により切り換えられた2分周のクロックによるサンプリングレートで、入力された画像信号Q1′をデジタル画像信号S1に変換する。
【0189】
このように本実施形態の画像情報読取装置によれば、レーザー光の主走査速度を変えることなく、高密度読取りの場合に比して、レーザー光Lの主走査方向について1/2、副走査方向について1/2の画素数である、全体として1/4の画素密度の低密度画像信号を得ることができる。しかも、特性変更手段60が、画素密度を規定するパラメータに応じて、走査レンズ21の位置、電圧印加手段39aへの制御信号、D/A変換器42aのサンプリングクロック、対数増幅器16aの周波数特性、およびアンチエリアジングフィルター35aをそれぞれ変更するため、読取り画素密度に適した、画素当たりのエネルギー付与、光電検出感度の設定、シェーディング補正、エリアジングノイズの抑制等が施された画像信号を得ることができる。
【0190】
なお、本実施形態の画像情報読取装置においては、高密度読取りが初期設定として設定されているものとしたが、本発明の画像情報読取方法および装置は、この態様に限られるものではなく、初期設定として「標準密度読取り」が、設定されており、「高密度読取り」、「低密度読取り」に切り換えられるものとしてもよい。
【0191】
例えば、主走査周波数(主走査速度)を160Hz、画素密度10pix /mmを標準密度読取りとしたときの、主走査サンプリング間隔、副走査画素ピッチ、副走査速度、アンチエリアジングフィルターカットオフ周波数を以下に示す。
【0192】
主走査サンプリング間隔: 100μm 副走査画素ピッチ: 100μm
主走査サンプリング周期:1.0μsec 副走査速度:16mm/sec
アンチエリアジングフィルターカットオフ周波数:500kHz
(アナログフィルタのため、正確には500kHz以下であることが好ましく、例えば400kHzなどを適用するのが望ましい。)
これに対して、画素密度20pix /mmである高密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔:50μm 副走査画素ピッチ:50μm
主走査サンプリング周期:0.5μsec 副走査速度:8mm/sec
アンチエリアジングフィルターカットオフ周波数:1000kHz
となり、一方、画素密度5pix /mmである低密度読取りのときは、主走査周波数を160Hz不変として、
主走査サンプリング間隔: 200μm 副走査画素ピッチ: 200μm
主走査サンプリング周期:2.0μsec 副走査速度:32mm/sec
アンチエリアジングフィルターカットオフ周波数:250kHz
となる。
【0193】
なお、上述したような「高密度読取り」、「低密度読取り」等、の読取り密度が段階的に予め設定されたものに限らず、読取り密度を規定するパラメータ(m,n)を任意に選択できるものとしてもよい。
【図面の簡単な説明】
【図1】本発明の第1の画像情報読取装置の一実施形態の構成を示す図
【図2】クロックを変更する構成の一実施形態を示す図
【図3】図1に示した実施形態における特性変更手段の構成を示す図
【図4】図1に示した実施形態におけるメモリ管理の詳細を説明する図
【図5】本発明の第2の画像情報読取装置の一実施形態を示す図
【図6】図5に示した実施形態における特性変更手段の構成を示す図
【図7】クロックを変更する構成の一実施形態を示す図
【図8】フィルタリング処理を説明するための図
【図9】本発明の第3の画像情報読取装置の一実施形態を示す図
【符号の説明】
1 蓄積性蛍光体シート
8 モータ
9a,9b エンドレスベルト
10 レーザ光源
11 レーザ光
12 回転多面鏡
14a,14b 光ガイド
15a,15b フォトマルチプレイヤ
16a,16b 対数増幅器
17a,17b レーザ光カットフィルタ
20 モータ
21 走査レンズ
35a,35b アンチエリアジングフィルタ
36a,36b A/D変換器
37a,37b 画素密度変換手段
38 加算手段
39a,39b 電圧印加手段
41a,41b メモリ
42a,42b A/D変換器
43a,43b 加算器
60 特性変更手段
70 入力手段[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an image information reading method and apparatus, and more particularly to an image information reading method and apparatus that can change the pixel density of a read image signal.
[0002]
[Prior art]
When irradiated with radiation (X-rays, α-rays, β-rays, γ-rays, electron beams, ultraviolet rays, etc.), a part of this radiation energy is accumulated, and when irradiated with excitation light such as visible light after that, depending on the accumulated energy Using a stimulable phosphor (stimulable phosphor) exhibiting stimulating luminescence, radiation image information of a subject such as a human body is temporarily recorded on a sheet-like stimulable phosphor (storable phosphor sheet), The stimulable phosphor sheet is scanned with a light beam such as a laser beam to generate stimulated emission light according to image information as signal light, and the obtained stimulated emission light is photoelectrically read by a photoelectric reading means. An analog image signal is obtained, a digital image signal having a predetermined pixel density can be obtained by sampling and quantizing the analog image signal at a predetermined interval, and a photographic photosensitive material based on the digital image signal. A radiation image information recording / reproducing system for outputting a radiation image of a subject as a visible image on a display device such as a recording medium or a CRT is already known (Japanese Patent Laid-Open Nos. 55-12429, 56-11395, 56). -11397 etc.).
[0003]
This system has a practical advantage that an image can be recorded over a very wide radiation exposure range as compared with a radiographic system using a conventional silver salt photograph.
[0004]
In addition, as a method of photoelectrically reading the above-mentioned stimulated emission light, photoelectric reading means are separately provided on both sides of the stimulable phosphor sheet, and excitation light is irradiated only on both sides or one side of the stimulable phosphor sheet. In addition, a radiation image reading apparatus has been proposed in which the stimulated emission light emitted from both sides of the stimulable phosphor sheet by the irradiation of the excitation light is photoelectrically detected by each photoelectric reading means (for example, JP-A-55). -87970). Such a radiation image reading apparatus has photoelectric reading means arranged so that one radiation image is accumulated and recorded on the stimulable phosphor sheet, and reading is performed from both sides of the stimulable phosphor sheet, and thereby, the light emission is stimulated. Light detection improves light collection efficiency, and by adding the two obtained image signals at a predetermined addition ratio, the detection position of the noise component is dispersed on the front and back of the sheet, so it can be obtained from only one side. Thus, it is possible to obtain an added image signal having a relatively improved S / N compared to the image signal.
[0005]
Furthermore, the S / N of the image signal (including the added image signal) is increased with respect to a single image signal obtained from one side of the sheet or two image signals obtained from both sides of the sheet. An image superposition method has been proposed in which an added image signal is obtained after filtering processing using a filter having a proper frequency response characteristic (Japanese Patent Laid-Open No. 7-287330). According to this method, the amount of radiation applied to the subject at the time of imaging is obtained, and the parameters (filter coefficients) for performing filtering processing on the image signal are obtained based on the amount of radiation, so that the subject is irradiated. An image signal representing an image with an optimum image quality or an added image signal representing a superimposed image can be obtained according to the radiation dose. Furthermore, in this method, since processing for changing the frequency characteristics is performed on the entire image signal, it is not necessary to perform frequency conversion such as Fourier transform, and the amount of calculation can be reduced.
[0006]
[Problems to be solved by the invention]
By the way, in the image reading unit of the radiation image information recording / reproducing system or in the radiation image information reading apparatus in which the image reading unit is configured alone, the pixel density of the obtained single image signal or the added image signal is described above. There is a demand to make the pixel density different from the predetermined pixel density.
[0007]
This is because, for example, when a large amount of radiographic images are taken, such as in a mass screening, it is required to improve the reading speed of the radiographic images taken, while the image quality required for the reproduced image is required. The reason is that it is sufficient to determine whether or not re-examination is necessary, and it is because there is no need to read image information at an excessively high pixel density. Of course, on the contrary, there is a case where it is required to read image information at a high pixel density so that a high-quality image can be reproduced even if the reading speed is lowered.
[0008]
In order to perform reading at other pixel densities in addition to reading at a predetermined pixel density set in advance as described above, the main scanning speed and the sub-scanning speed of the light beam can be simply changed. Conceivable. When changing the main scanning speed of the light beam, it is necessary to change the driving speed of a scanning optical system (for example, a polygon mirror or a galvanometer mirror) that performs the main scanning of the light beam.
[0009]
However, when the drive speed of the scanning optical system is changed, a certain time is required until the drive speed is stabilized due to the influence of the inertial force, and stable drive may not be guaranteed in all speed ranges. . For this reason, it is desirable that the pixel density can be changed without changing the main scanning speed of the light beam. Further, when the pixel density is changed in this way, the pixel densities in the main scanning direction and the sub-scanning direction need to be changed at the same rate.
[0010]
In addition, when changing the pixel density, it is desired to perform the pixel density conversion process as fast as possible.
[0011]
Furthermore, the following problems may occur when the pixel density is simply changed.
[0012]
That is,
(I) Since the energy of the signal light per pixel emitted from the sheet by changing the pixel density is different from the energy of the signal light per pixel at the original pixel density, the pixels before and after the change If the density ratio is large, the density (or luminance) of the entire reproduced image may change, which may affect the diagnostic performance.
[0013]
(II) When shading correction is performed on an analog image signal, there may be a case where the shading characteristics to be corrected are changed before and after the pixel density change, and the shading cannot be corrected appropriately.
[0014]
(III) When the analog image signal is logarithmically amplified, the frequency transfer characteristics may change before and after the pixel density is changed.
[0015]
(IV) When performing filtering that removes aliasing noise prior to sampling the analog image signal, the Nyquist frequency changes before and after the pixel density change, so the aliasing noise may not be cut properly. is there.
[0016]
These problems do not arise only in the reading of radiation image information from the stimulable phosphor sheet on which the radiation image information is recorded by the above-described radiation image information recording / reproducing system. The same may occur in a method or apparatus that scans with a beam and reads reflected light emitted as a signal light according to an image printed on the medium.
[0017]
The present invention has been made in view of the above circumstances, and a first object thereof is to provide an image information reading method and apparatus capable of changing the pixel density without changing the main scanning speed of the light beam. Is.
[0018]
A second object of the present invention is to provide an image information reading method and apparatus for double-sided reading that can change the pixel density without changing the main scanning speed of the light beam.
[0019]
[Means for Solving the Problems]
The image information reading method and apparatus of the present invention obtains a digital image signal having an arbitrary pixel density by changing the sub-scanning speed and sampling interval without changing the main scanning speed.
[0020]
That is, according to the first image information reading method of the present invention, a light beam is repeatedly main-scanned at a constant speed in a constant direction on a scanned sheet on which image information is recorded, and substantially orthogonal to the main-scanning direction. By sub-scanning in the direction, the light beam is scanned over substantially the entire surface of the scanned sheet, and the signal light emitted from the scanned sheet is photoelectrically detected by scanning the light beam, and the obtained analog image In an image information reading method for obtaining a digital image signal having a predetermined pixel density by sampling and quantizing a signal,
The sub-scanning speed is multiplied by m (m> 0), and the sampling interval is multiplied by n (n> 0).
[0021]
Here, it is preferable to further perform a pixel density conversion process for a / m (a> 0) times the number of pixels in the main scanning direction and a / n times the number of pixels in the sub-scanning direction. This is because the aspect ratio of the image can be made the same as that of the original image.
[0022]
The second image information reading method of the present invention relates to an added image signal of two image signals obtained by reading from both sides of the sheet, and a constant light beam is applied to the scanned sheet on which the image information is recorded. The light beam is scanned over substantially the entire surface of the scanned sheet by performing sub-scanning in a direction substantially perpendicular to the main scanning direction while repeatedly performing main scanning at a constant speed in the direction of. The signal light emitted from both sides of the scanned sheet is detected photoelectrically, the two analog image signals obtained are sampled and quantized to obtain two digital image signals, and the two digital image signals are added In an image information reading method for obtaining an added image signal having a predetermined pixel density,
The sub-scanning speed is multiplied by m (m> 0), the sampling interval is multiplied by n (n> 0), and the number of pixels in the main scanning direction is a / m (a> 0) times. The pixel density conversion process is performed to increase the number of pixels in the scanning direction by a / n times.
[0023]
Here, in the first and second inventions, the scanned sheet includes not only the stimulable phosphor sheet on which the radiation image is stored and recorded, but also a photographic print or the like used in the radiation image information recording / reproducing system described above. Also includes transparent originals such as reflective originals and films. Therefore, the signal light emitted from the scanned sheet includes reflected light emitted from the reflective original and transmitted light emitted from the transparent original, in addition to the stimulated light emitted from the storage phosphor sheet by being excited by the light beam. Including.
[0024]
Further, as a method of changing the sampling interval, a clock divided by a frequency dividing circuit is created with respect to the reference clock, and sampling is performed with this divided clock, or the period is different depending on the PLL with respect to the reference clock. It is possible to apply a method such as creating a clock and sampling with the new clock, or switching and applying a plurality of sampling clocks having different periods.
[0025]
Further, as a method of obtaining the added image signal, not only a method of simply adding pixels corresponding to each other of two digital image signals but also two digital images as described in the above-mentioned Japanese Patent Application Laid-Open No. 7-287330. It is possible to employ a method of performing addition after performing a filtering process with a filter having a frequency response characteristic that increases the S / N of the added image signal with respect to the signal.
[0026]
On the other hand, as the pixel density conversion processing, a method by performing a one-dimensional mask operation in the main scanning direction of an image signal, a method of thinning out pixels according to the pixel density, a B-spline interpolation operation, a Cubic spline interpolation operation, and the like The following interpolation calculation methods (for example, JP-A-8-16767 and 9-321981) and linear interpolation calculation methods (for example, JP-A-9-50516) can be employed.
[0027]
When performing this pixel density conversion, it is preferable to change parameters for one-dimensional mask calculation, parameters for thinning out pixels, and parameters for calculation formulas for interpolation calculation according to m and n. Here, the parameter for the one-dimensional mask calculation is a mask coefficient, and the parameter for thinning out pixels is the thinning-out interval of pixels. The parameter of the calculation formula of the interpolation calculation represents what type of interpolation calculation (for example, B-spline interpolation calculation, Cubic spline interpolation calculation, linear interpolation calculation) is used when performing the interpolation calculation.
[0028]
In the second image information reading method of the present invention, it is preferable to perform the pixel density conversion process on the two digital image signals. Here, performing pixel density conversion processing on two digital image signals means performing pixel density conversion before adding two digital image signals to obtain an added image signal.
[0029]
In each image information reading method of the present invention, it is preferable to change at least one of the following characteristics (1) to (8) according to m and n.
[0030]
(1) Diameter of a light beam that scans the scanned sheet
(2) Power of a light beam for scanning the scanned sheet
(3) Sensitivity to detect the signal light
(4) Shading correction data set in advance when shading correction is performed on the analog image signal
(5) Timing of outputting the shading correction data from the memory storing the shading correction data
(6) Frequency transfer characteristics when the analog image signal is logarithmically amplified
(7) Cut-off frequency when performing filtering to remove aliasing noise prior to sampling the analog image signal
(8) Parameters for filtering processing in the pixel density conversion processing (including thinning parameters when performing thinning processing as filtering processing)
Here, the changed items (1) to (8) are applied to at least one applicable item, and it is not always necessary to provide all the items (1) to (8) as configuration requirements. That is, for example, in an image information reading method that does not perform shading correction, it is not necessary to provide the configurations (4) and (5).
[0031]
Furthermore, it is only necessary to change at least one of the items (1) to (8), and only one of them may be included as a configuration requirement, or two or more of these may be included. .
[0032]
The items (1) to (3) are changed when the above-mentioned problem (I), that is, when the energy of the signal light per pixel emitted from the sheet by changing the pixel density is the original pixel density. Therefore, if the pixel density ratio before and after the change is large, the density (or brightness) of the entire image to be reproduced will change, affecting the diagnostic performance. Specifically, when the pixel density is changed to the high density side, (1) the light beam diameter is reduced, (2) the light beam power is reduced, and / or ( 3) When the sensitivity for detecting the signal light is increased while the pixel density is changed to the lower density side, (1) the light beam diameter is enlarged, (2) the power of the light beam is increased, and / or ( 3) Reduce the sensitivity to detect signal light Just do it.
[0033]
Similarly, when the pixel density is changed to the high density side, the item (6) is changed so that the frequency response characteristic is expanded to a higher frequency band, and the item (7) is also set to increase the cutoff frequency. When the pixel density is changed to the lower density side, the item (6) is changed so that the cutoff frequency is set to the lower frequency side, and the item (7) is similarly applied. The cut-off frequency may be changed to the low frequency side.
[0034]
The above item (8) is such that, for example, the filtering frequency (coefficient for mask calculation) is set so that the cutoff frequency when the mask calculation is performed at the low density side is lower than that at the high density side. Change it. Note that, when applying the thinning process parameter, it may be changed in a direction opposite to the change direction with respect to the pixel density level when the filtering process parameter is applied.
[0035]
Since item (4) above is not a characteristic that varies qualitatively according to the pixel density, shading correction data is prepared in advance for each representative pixel density, and shading corresponding to the changed pixel density. Data for correction is selected from these prepared data, or when data having a matching pixel density is not prepared, data corresponding to at least two types of pixel densities is prepared from the prepared data. The shading correction data corresponding to the changed pixel density may be obtained and applied by the interpolation processing used.
[0036]
The above item (5) is applied when shading correction is performed in real time on a digital image signal, and is obtained by outputting shading correction data prepared in advance at a timing according to the pixel density. The digital image signal is subjected to shading correction. Specifically, when the pixel density is changed to the high density side, the sampling speed increases, so the timing for outputting the shading correction data may be advanced according to the sampling speed. On the other hand, when the pixel density is changed to the low density side, the sampling speed is slowed down. Therefore, the timing for outputting the shading correction data may be delayed according to the sampling speed.
[0037]
Here, shading refers to light caused by unevenness in the intensity of the scanning beam due to uneven reflectance on the reflecting surface of an optical deflector (polygon mirror, galvanometer mirror, etc.) for scanning the light beam, and variations in the deflection speed of the optical deflector. Variations in the analog image signal obtained from the photoelectric reading means (part of the light detection efficiency) due to unevenness in the scanning speed of the beam or detection unevenness due to sensitivity unevenness in the main scanning direction of the photoelectric detecting means arranged along the main scanning direction. Mean decrease). The shading correction data refers to shading characteristics obtained in advance using, for example, a reference scanned sheet (for example, a stimulable phosphor sheet that is uniformly irradiated with radiation). (See JP 61-189763, 62-47259, 62-47261, 64-86759, JP-A-2-58973, etc.).
[0038]
In addition, in the above-mentioned changes (1), (2), (3), (6), (7), and (8), the characteristics for each pixel density (the beam diameter, Beam power in (2), sensitivity in (3), frequency response characteristics in (6), frequency response characteristics with different cutoff frequencies in (7), (8) In this case, the parameter may be changed by preparing a filtering parameter) and selecting a characteristic corresponding to the changed pixel density.
[0039]
Further, (4) as a method for changing the shading correction data, a method of selecting from a plurality of correction data prepared in advance as described above, and a method corresponding to a preset pixel density are prepared. One shading correction data may be obtained by sampling according to the changed pixel density, and the newly obtained shading correction data may be used.
[0040]
For each value of m and n, for example, “high density” (for example, when m = n = 0.5), “standard density” (original predetermined pixel density (corresponding to m = n = 1)), A correspondence relationship of “low density” (when m = n = 2) may be set.
[0041]
Further, the shading correction data is changed by selecting the shading correction data corresponding to the selected values of m and n from the shading correction data set in advance for each value of m and n that can be selected. In this case, it is preferable to apply one of the following two methods [I] and [II]. That is,
[I] Shading correction data set in advance for each value of m, n is stored in the first storage medium in advance, and shading correction data corresponding to the selected value of m, n is selected each time. The shading correction data is selected by transferring the shading correction data from the first storage medium to the second storage medium and reading the transferred shading correction data from the second storage medium.
[II] Pre-set shading correction data for each value of m and n is stored in advance in the first storage medium, and when desired, such as at the time of startup, all of these shading correction data are Transfer from the storage medium to the second storage medium having a storage area corresponding to a different address for each value of m, n, corresponding to the value of m, n, and corresponding to the selected value of m, n The shading correction data is selected by reading the shading correction data to be read from the address corresponding to the selected values of m and n.
[0042]
According to the method [I], the capacity of the second storage medium (memory) can be reduced, and the hardware configuration can be simplified. On the other hand, according to the method [II], the software configuration can be simplified and the reading speed of the correction data from the second storage medium can be increased.
[0043]
According to the third image information reading method of the present invention, a light beam is repeatedly main-scanned at a constant speed in a constant direction on a scanned sheet on which image information is recorded, in a direction substantially perpendicular to the main scanning direction. By sub-scanning, the light beam is scanned over substantially the entire surface of the scanned sheet, signal light emitted from the scanned sheet is detected photoelectrically by scanning the light beam, and the resulting analog image signal is In an image information reading method for obtaining a digital image signal having a predetermined pixel density by sampling and quantization,
The sub-scanning speed is multiplied by m (m> 0), the sampling interval is multiplied by n (n> 0), and depending on these m and n, the following characteristics (1) to (5) At least one of them is changed to obtain a digital image signal having a
[0044]
(1) Diameter of a light beam that scans the scanned sheet
(2) Sensitivity to detect the signal light
(3) Shading correction data set in advance when shading correction is performed on the analog image signal
(4) Frequency transfer characteristics when the analog image signal is logarithmically amplified
(5) Cutoff frequency when performing filtering to remove aliasing noise prior to sampling the analog image signal
Items (1) to (5) are the same as described above.
[0045]
A first image information reading apparatus of the present invention is an apparatus for carrying out the first image information reading method of the present invention, and is a light beam emitted from a light source on a scanned sheet on which image information is recorded. Main scanning means for repeatedly performing main scanning repeatedly in a constant direction at a constant speed, and the sheet or the sheet or the sheet so that the sheet is relatively sub-scanned in a direction substantially perpendicular to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the light beam, photoelectric detection means for photoelectrically detecting signal light emitted from the scanned sheet by scanning the light beam, and sampling and quantizing the detected analog image signal In an image information reading apparatus provided with an A / D conversion means for obtaining a digital image signal having a predetermined pixel density,
The sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m times according to the inputted m, and the sampling interval by the A / D converting means according to the inputted n And a sampling interval changing means for changing the sampling interval to n times.
[0046]
Here, it is preferable to further include a pixel density conversion processing unit that performs a pixel density conversion process in which the number of pixels in the main scanning direction is a / m times and the number of pixels in the sub-scanning direction is a / n times.
[0047]
A second image information reading apparatus of the present invention is an apparatus for carrying out the second image information reading method of the present invention, and is emitted from a light source onto a scanned sheet on which image information is recorded. Main scanning means for repeatedly scanning the light beam at a constant speed in a constant direction, and the sheet so that the sheet is relatively sub-scanned in a direction substantially perpendicular to the main scanning direction of the light beam. Alternatively, sub-scanning means for sub-scanning the light beam, photoelectric detection means for photoelectrically detecting signal light emitted from both sides of the scanned sheet by scanning the light beam, and two detected analog image signals Image information reading comprising: A / D conversion means for obtaining two digital image signals by sampling and quantization; and addition means for obtaining an added image signal having a predetermined pixel density by adding the two digital image signals In the location,
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m times according to the input m;
Sampling interval changing means for changing the sampling interval by the A / D conversion means to n times according to the input n;
Pixel density conversion processing means for performing pixel density conversion processing for increasing the number of pixels in the main scanning direction by a / m (a> 0) times and setting the number of pixels in the sub-scanning direction by a / n times. It is what.
[0048]
In the second image information reading apparatus of the present invention, it is preferable that the pixel density conversion processing means is a means for performing pixel density conversion processing on the two digital image signals.
[0049]
Each image information reading apparatus of the present invention preferably further includes a characteristic changing unit that changes at least one of the following characteristics (1) to (8) according to m and n.
[0050]
(1) Diameter of a light beam that scans the scanned sheet
(2) Power of a light beam for scanning the scanned sheet
(3) Sensitivity to detect the signal light
(4) Shading correction data set in advance when shading correction is performed on the analog image signal
(5) Timing of outputting the shading correction data from the memory storing the shading correction data
(6) Frequency transfer characteristics when the analog image signal is logarithmically amplified
(7) Cut-off frequency when performing filtering to remove aliasing noise prior to sampling the analog image signal
(8) Parameters for filtering processing in the pixel density conversion processing (including thinning parameters when performing thinning processing as filtering processing)
The sub-scanning speed changing unit and the sampling interval changing unit may be part of the characteristic changing unit.
[0051]
When the characteristic changing means includes one that changes the shading correction data, among the shading correction data set in advance for each m and n value that can be selected, the shading correction according to the selected m and n values. The shading correction data may be changed by selecting data. In this case, it is preferable to apply one of the following two methods [I] and [II]. That is,
[I] Shading correction data set in advance for each value of m, n is stored in the first storage medium in advance, and shading correction data corresponding to the selected value of m, n is selected each time. The shading correction data is selected by transferring the shading correction data from the first storage medium to the second storage medium and reading the transferred shading correction data from the second storage medium.
[II] Pre-set shading correction data for each value of m and n is stored in advance in the first storage medium, and when desired, such as at the time of startup, all of these shading correction data are Transfer from the storage medium to the second storage medium having a storage area corresponding to a different address for each value of m, n, corresponding to the value of m, n, and corresponding to the selected value of m, n The shading correction data is selected by reading the shading correction data to be read from the address corresponding to the selected values of m and n.
[0052]
According to the method [I], the capacity of the second storage medium (memory) can be reduced, and the hardware configuration can be simplified. On the other hand, according to the method [II], the software configuration can be simplified and the reading speed of the correction data from the second storage medium can be increased.
[0053]
The third image information reading apparatus of the present invention comprises a main scanning unit that repeatedly scans a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction; Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in a direction substantially perpendicular to the main scanning direction of the beam; Image information comprising photoelectric detection means for photoelectrically detecting signal light emitted from the sheet, and A / D conversion means for sampling and quantizing the detected analog image signal to obtain a digital image signal having a predetermined pixel density In the reading device,
Pixel density input means for receiving inputs of values m and n corresponding to a pixel density of 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction with respect to the predetermined pixel density;
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m (m> 0) times;
Sampling interval changing means for changing the sampling interval by the A / D converter to n (n> 0) times;
It is characterized by comprising characteristic changing means for changing at least one of the following (1) to (5) according to m and n.
[0054]
(1) Diameter of a light beam that scans the scanned sheet
(2) Sensitivity to detect the signal light
(3) Shading correction data set in advance when shading correction is performed on the analog image signal
(4) Frequency transfer characteristics when the analog image signal is logarithmically amplified
(5) Cutoff frequency when performing filtering to remove aliasing noise prior to sampling the analog image signal
[0055]
【The invention's effect】
According to the image information reading method and apparatus of the present invention, the number of scanning lines (main scanning lines) of a light beam for scanning a scanned sheet is increased by multiplying the sub-scanning speed by m (m> 0)) times. The pixel density in the sub-scanning direction can be increased to 1 / m times because it increases or decreases to m times, and the pixel density in the main scanning direction is increased to 1 / n times by increasing the sampling interval by n (n> 0) times. Therefore, it is possible to obtain an image signal having a pixel density that is 1 / (m × n) 2 times the finally obtained image signal (added image signal in the case of double-sided reading). Further, since it is not necessary to change the main scanning speed of the light beam even when the pixel density is changed, it is not necessary to change the driving speed of a scanning optical system (for example, a polygon mirror or a galvanometer mirror) that scans the light beam. Therefore, it is possible to eliminate the unreadable time until the drive speed is stabilized, which occurs when the drive speed of the scanning optical system is changed.
[0056]
Further, in the invention in which the pixel density conversion processing is performed in which the number of pixels in the main scanning direction is a / m (a> 0) times and the number of pixels in the sub-scanning direction is a / n times, the pixel density conversion process is finally obtained. The pixel density in the main scanning direction of the image signal (added image signal in the case of double-sided reading) can be a / (m × n) times, and the pixel density in the sub-scanning direction can also be a / (m × n) times. The change rate of the pixel density in the main scanning direction and the sub-scanning direction is the same, and it is set in advance without changing the main scanning speed of the light beam while maintaining the aspect ratio (main / sub ratio) of the image size. (A / (m × n)) of this pixel density 2 An image signal having a double pixel density can also be obtained.
[0057]
In particular, when performing double-sided reading for reading signals from both sides of the scanned sheet (the second image information reading method and apparatus of the present invention), in order to sufficiently apply energy by the light beam to the back side of the sheet. However, it is necessary to slow down the scanning speed as compared with the method of reading only from one side, but even in this case, if it is allowed to reduce the pixel density, the sub-scanning speed can be increased. The scanning time per sheet can be further shortened.
[0058]
In the case of double-sided reading, by performing pixel density conversion processing on two digital image signals, the image signals subjected to pixel density conversion processing can be added to obtain an added image signal. For this reason, it is possible to reduce the amount of addition calculation when adding digital image signals, thereby shortening the calculation time of the addition calculation and performing processing at high speed.
[0059]
Furthermore, according to the image information reading method and apparatus of the present invention, it is possible to solve various problems that occur when the reading pixel density is changed.
[0060]
That is, (I) Since the energy of signal light per pixel emitted from the sheet by changing the pixel density is different from the energy of signal light per pixel at the original pixel density, before and after the change When the pixel density ratio is large, the density (or luminance) of the entire image to be reproduced may change and affect the diagnostic performance. According to the image information reading method and apparatus of the present invention, Depending on the pixel density, (1) the energy applied to the unit area of the sheet can be changed by changing the diameter of the light beam that scans the scanned sheet, and (2) the light beam according to the pixel density. By changing the power, the energy applied to the unit area of the sheet can be changed, and (3) by changing the sensitivity for detecting the signal light according to the pixel density. It can be the level of the image signal is energy of the signal light obtained be varied to suppress a change.
[0061]
(II) When shading correction is performed on an analog image signal, the shading characteristics to be corrected may change before and after the pixel density change, and the shading may not be corrected appropriately. According to this image information reading method and apparatus, (4) by changing the shading correction data according to the pixel density, it is possible to perform the shading correction corresponding to the change in the shading characteristics, and according to the pixel density. (5) By changing the timing of output of the shading correction data from the memory, the shading correction corresponding to the change in the sampling rate can be performed.
[0062]
Furthermore, (III) when the analog image signal is logarithmically amplified, there may be a change in the frequency transfer characteristic before and after the pixel density change, but according to the image information reading method and apparatus of the present invention, the pixel density (6) By changing the frequency transfer characteristic when the analog image signal is logarithmically amplified, it is possible to suppress the change in the frequency characteristic of the obtained image signal.
[0063]
In addition, (IV) When performing filtering to remove aliasing noise prior to sampling the analog image signal, the Nyquist frequency changes before and after the pixel density change, so the aliasing noise cannot be cut properly. In some cases, according to the image information reading method and apparatus of the present invention, the frequency of the image signal obtained by changing the cutoff frequency of (7) filtering to remove aliasing noise according to the pixel density. It can suppress that a characteristic changes.
[0064]
Further, when performing the pixel density conversion process, (8) the aliasing noise at the time of the pixel density conversion is suppressed to a certain level or less by changing the parameter of the filtering process (masking coefficient) in the pixel density conversion process. Can do. In addition, errors that can occur due to interpolation processing such as spline interpolation calculation can be reduced.
[0065]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, specific embodiments of the image information reading apparatus of the present invention will be described with reference to the drawings.
[0066]
FIG. 1 is a diagram showing an embodiment of an image information reading apparatus according to the present invention. In the illustrated image information reading device, a stimulable phosphor sheet (hereinafter referred to as a sheet) 1 on which radiation image information is accumulated and recorded is disposed on an
[0067]
Further, immediately above the
[0068]
Further, a logarithmic amplifier 16a is connected to the photomultiplier 15a, and the analog image signal QA detected by the photomultiplier 15a is logarithmically amplified according to a preset frequency characteristic to obtain a logarithmic image signal QA. ′ Is output.
[0069]
On the other hand, shading correction data D1 corresponding to a preset sampling interval is stored in the
[0070]
Further, an anti-aliasing filter 35a that removes aliasing noise (folding noise) caused by A / D conversion, which will be described later, is connected to the adder 43a, and the filtered image signal Q1 'is connected to the subsequent stage. Is converted to a digital image signal S1 with a preset reference clock. The anti-aliasing filter 35a includes a high-density filter and a low-density filter. The high-density filter is initially selected. However, the anti-aliasing filter 35a has a low density according to an input signal from the characteristic changing
[0071]
The A /
[0072]
Furthermore, the image information reading apparatus of the present embodiment has a “high density” (10 pix / mm) that is initially set as a pixel density for reading image information, and a lower density than this high density. It receives the density selected by the operator from “Low Density” (5 pix / mm). When “High Density”, the parameter (m, n) is set to m = n = 1, and “Low Density” In some cases, the input means 70 outputs the parameter (m, n) as m = n = 2, and the parameter (m, n) output from the input means 70 is input to these parameters (m, n). Accordingly, the rotational speed of the
[0073]
Here, the pixel density conversion unit 37a performs a pixel density conversion process for 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction of the digital image signal S1 in accordance with the parameter (m, n). As a result, an image signal S1 ′ having a
[0074]
Next, the operation of the radiation image information reading apparatus of this embodiment will be described.
[0075]
First, a signal representing “high density” selected by the operator is input to the input means 70. “High density” is an initial setting of the image information reading apparatus, and a parameter (m, n) = (1, 1) corresponding to “high density” is output from the input means 70 and input to the characteristic changing means 60. Is done.
[0076]
The characteristic changing means 60, according to the input parameter (m, n) = (1, 1), the rotational speed of the
[0077]
Next, the
[0078]
On the other hand, the
[0079]
From the portion of the
[0080]
The PMT 15a is applied with a high voltage corresponding to the sensitivity for high density by the
[0081]
On the other hand, the high-density shading correction data D1 stored in the
[0082]
Here, the clock input to the D / A converter 42a is a reference clock that is switched and output by a
[0083]
The shading correction analog signal D1 'is input to the adder 43a and added to the logarithmic image signal QA' input from the logarithmic amplifier 16a, whereby the logarithmic image signal QA 'is subjected to shading correction. And is input to the anti-aliasing filter 35a.
[0084]
The anti-aliasing filter 35a is switched to a high-density filter by the characteristic changing means 60, and the input image signal Q1 is appropriately freed of aliasing noise by the high-density filter, and the A /
[0085]
The pixel density conversion unit 37a performs pixel density conversion by multiplying the digital image signal S1 by 1 / m times in the main scanning direction and 1 / n times in the sub-scanning direction. Here, m = n = 1. Therefore, the image signal S1 ′ is obtained without performing any pixel density conversion processing and is output to the image processing apparatus or the like.
[0086]
Next, an operation when “low density” is selected by the operator will be described.
[0087]
When a signal representing “low density” is input to the input means 70 by the operator, the input means 70 sets a parameter (m, n) = (2, 2) corresponding to “low density” as shown in FIG. Output to the characteristic changing means 60. The characteristic changing means 60, according to the input parameter (m, n) = (2, 2), the rotational speed of the motor that drives the
[0088]
In this way, the characteristics of each part are changed, and image information is read in the same manner as in the case of the high-density reading described above. That is, the
[0089]
The stimulated emission light 13a emitted from the portion of the
[0090]
The voltage is applied to the PMT 15a by the
[0091]
On the other hand, the high-density shading correction data D1 stored in the
[0092]
The shading correction analog signal D1 'is input to the adder 43a and added to the logarithmic image signal QA' input from the logarithmic amplifier 16a, whereby the logarithmic image signal QA 'is subjected to shading correction. ′ And input to the anti-aliasing filter 35a.
[0093]
The anti-aliasing filter 35a is switched to a low-density filter by the characteristic changing means 60, and the input image signal Q1 'is appropriately freed of aliasing noise by the low-density filter so that an A / D converter is used. It is input to 36a. The A /
[0094]
The pixel density conversion unit 37a performs pixel density conversion of 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction of the digital image signal S1. Here, since m = n = 2, pixel density conversion processing is performed in which the number of pixels is halved in the main scanning direction and the sub-scanning direction of the digital image signal S1. As this pixel density conversion process, a one-dimensional mask calculation will be described here. An example of a one-dimensional filter used for this one-dimensional mask operation is shown below.
[0095]
a (x, 1) = (-8 / 105, -5 / 105,34 / 105,63 / 105,34 / 105, -5 / 105, -8 / 105)
The filter a (x, 1) is used to perform filtering processing at intervals of 1 pixel in the main scanning direction of the digital image signal S1, and then perform filtering processing at intervals of 1 pixel in the sub-scanning direction to perform pixel density conversion. Done. Specifically, assuming that the pixel value of the digital image signal S1 is S1 (x, y) and the pixel value after pixel density conversion in the main scanning direction is S1A (x / 2, y), the pixel value S1A (x / 2 , Y) is calculated by the following equation (1).
[0096]
However, k = 1 to N (N is the number of pixels in the main scanning direction (pixel position or pixel number))
l = 1 to M (M is the number of pixels in the sub-scanning direction (pixel position or pixel number))
Similarly, the image signal S1 ′ having undergone pixel density conversion is obtained by performing the calculation of (1) also in the sub-scanning direction.
[0097]
On the other hand, in such a mask calculation using a one-dimensional filter, data for performing the mask calculation is insufficient at the edge portion of the digital image signal S1. Therefore, when performing filtering processing on pixels at such edge portions, a virtual pixel position is determined outside the pixel, and an appropriate pixel value is assigned to each pixel at the virtual pixel position. The filtering process may be performed with
[0098]
The pixel density conversion in the pixel density conversion unit 37a is not limited to the mask calculation using the one-dimensional filter, and the pixel density conversion may be performed by a process of thinning out pixels or an interpolation calculation. Here, as the interpolation calculation, a high-order interpolation calculation such as a B-spline interpolation calculation focusing on smoothness and a Cubic spline interpolation calculation focusing on sharpness can be applied in addition to linear interpolation.
[0099]
The image signal S1 ′ thus obtained in the pixel density conversion means 37a has a pixel density of 1/16 when “low density” is selected, compared to when “high density” is selected. And output to an image processing apparatus or the like.
[0100]
As described above, according to the image information reading apparatus of the present embodiment, the sub-scanning is 1/2 in the main scanning direction of the laser beam L as compared with the case of high-density reading without changing the main scanning speed of the laser beam. A low-density image signal having a pixel density of ¼, which is ½ the number of pixels in the direction, is obtained. By performing the pixel density conversion process as described above, a low-density image signal having a pixel density of 1/16 as a whole can be obtained. In addition, the characteristic changing
[0101]
In the image information reading apparatus according to the present embodiment, high-density reading is set as an initial setting. However, the image information reading method and apparatus according to the present invention are not limited to this mode. “Standard density reading” is set as the setting, and it may be selectively switched between “high density reading” and “low density reading”.
[0102]
For example, when the main scanning frequency (main scanning speed) is 160 Hz and the pixel density is 10 pix / mm, the main scanning sampling interval, sub-scanning pixel pitch, sub-scanning speed, and anti-aliasing filter cutoff frequency are as follows: Shown in
[0103]
Main scanning sampling interval: 100 μm Sub-scanning pixel pitch: 100 μm
Main scanning sampling period: 1.0 μsec Sub scanning speed: 16 mm / sec
Anti-aliasing filter cutoff frequency: 500 kHz
(Because it is an analog filter, it is desirable to apply 500 kHz or less, for example, 400 kHz.)
On the other hand, at the time of high density reading with a pixel density of 20 pix / mm, the main scanning frequency is assumed to be 160 Hz invariant.
Main scanning sampling interval: 50 μm Sub-scanning pixel pitch: 50 μm
Main scanning sampling period: 0.5 μsec Sub scanning speed: 8 mm / sec
Anti-aliasing filter cutoff frequency: 1000 kHz
On the other hand, at the time of low density reading with a pixel density of 5 pix / mm, the main scanning frequency is set to 160 Hz invariant,
Main scanning sampling interval: 200 μm Sub-scanning pixel pitch: 200 μm
Main scanning sampling period: 2.0 μsec Sub scanning speed: 32 mm / sec
Anti-aliasing filter cutoff frequency: 250 kHz
It becomes.
[0104]
Note that the reading density such as “high density reading” and “low density reading” as described above is not limited to those set in advance, and parameters (m, n) that define the reading density are arbitrarily selected. It may be possible.
[0105]
Details of a method for changing the shading correction data stored in the
[0106]
As shown in FIG. 4, the
[0107]
Here, according to the method of transferring only specific shading correction data to the SHD memory 41d for each shading correction, the capacity of the SHD memory 41d can be reduced, and the hardware configuration can be simplified. On the other hand, according to the method of reading all the shading correction data from the
[0108]
FIG. 5 is a diagram showing an embodiment of the second image information reading apparatus of the present invention. As shown in FIG. 5, the image information reading apparatus according to the present embodiment obtains image signals representing a radiation image from both sides of a stimulable phosphor sheet on which a radiation image of a subject is accumulated and recorded. The
[0109]
On the other hand, shading correction data D1 and D2 corresponding to a preset sampling interval are stored in the
[0110]
Furthermore,
[0111]
Furthermore, the A /
[0112]
Further, as shown in FIG. 6, the image information reading apparatus according to the present embodiment has a “high density” (10 pix / mm) that is initially set as a read pixel density of image information, In response to the input of the density selected by the operator among “low density” (5 pix / mm), the parameters (m, n) are set to m = n = 1, In the case of “low density”, the input means 70 that outputs the parameter (m, n) with m = n = 2, and the parameter (m, n) output from the input means 70 are input, and these parameters (m , N), the rotational speed of the
[0113]
Here, in the pixel density conversion means 37a and 37b, the pixels which are 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction of the digital image signals S1 and S2 according to the parameters (m, n). 1 / (m × n) than when m = n = 1 after density conversion processing 2 Image signals S1 'and S2' having double pixel density are obtained. Specifically, a method by performing a one-dimensional mask operation in both the main and sub scanning directions of the digital image signals S1 and S2, a method of thinning out pixels according to the pixel density, a B-spline interpolation operation, a Cubic spline interpolation operation, etc. Pixel density conversion is performed by a high-order interpolation calculation method, a linear interpolation calculation method, or the like. At this time, the pixel density conversion parameter is changed according to the parameter (m, n). For example, the mask coefficient is changed when a one-dimensional mask calculation is performed, the interval between thinnings is performed when a thinning process is performed, and the type of interpolation calculation applied to the digital image signals S1 and S2 is changed when an interpolation calculation is performed.
[0114]
In addition, in the adding means 38, as described in the above-mentioned JP-A-7-287330, filtering processing by a filter having a frequency response characteristic that increases the S / N of the added image signal S3 is performed on the image signal S1 ′, Although it is preferable to adopt a method of obtaining the added image signal S3 after applying to S2 ', the image signals S1' and S2 'may be simply added. In addition, when performing a filtering process, you may make it perform in the pixel density conversion means 37a and 37b.
[0115]
Next, the operation of the radiation image information reading apparatus of this embodiment will be described.
[0116]
First, a signal representing “high density” selected by the operator is input to the input means 70. “High density” is an initial setting of the image information reading apparatus, and a parameter (m, n) = (1, 1) corresponding to “high density” is output from the input means 70 and input to the characteristic changing means 60. Is done.
[0117]
The characteristic changing means 60 determines the rotational speed of the
[0118]
When the
[0119]
These analog image signals QA and QB are input to
[0120]
On the other hand, the high-density shading correction data D1 and D2 stored in the
[0121]
Here, the clocks input to the D /
[0122]
The shading correction analog signals D1 'and D2' are input to the
[0123]
The
[0124]
The A /
[0125]
The pixel
[0126]
The adding means 38 is a filter having a frequency response characteristic for increasing the S / N of the added image signal S3 as described in Japanese Patent Laid-Open No. 7-287330 with respect to the image signals S1 'and S2'. A filtering process is performed, and the image signals S1 ′ and S2 ′ subjected to the filtering process are added together with corresponding pixels to obtain an added image signal S3, which is output to an image processing apparatus or the like.
[0127]
Next, an operation when “low density” is selected by the operator will be described.
[0128]
When a signal representing “low density” is input to the
[0129]
The characteristic changing means 60 determines the rotational speed of the
[0130]
Specifically, the rotational speed of the
[0131]
In this way, the characteristics of each part are changed, and image information is read in the same manner as in the case of the high-density reading described above.
[0132]
That is, the
[0133]
The photostimulated light 13a and 13b emitted from the portion of the
[0134]
On the other hand, the high-density shading correction data D1 and D2 stored in the
[0135]
The shading correction analog signals D1 'and D2' are input to the
[0136]
The
[0137]
The A /
[0138]
The pixel
[0139]
a (x, 1) = (-8 / 105, -5 / 105,34 / 105,63 / 105,34 / 105, -5 / 105, -8 / 105)
Then, after filtering the digital image signals S1 and S2 in the main scanning direction using the filter a (x, 1) at intervals of 1 pixel, the filtering processing is performed at intervals of 1 pixel in the sub-scanning direction. Conversion is performed. More specifically, the pixel value of the digital image signal S1 (S2 is the same, so only S1 will be described here) is set to S1 (x, y), and the pixel value after pixel density conversion in the main scanning direction is set to S1A. If (x / 2, y), the pixel value S1A (x / 2, y) is calculated by the following equation (1).
[0140]
However, k = 1 to N (N is the number of pixels in the main scanning direction (pixel position or pixel number))
l = 1 to M (M is the number of pixels in the sub-scanning direction (pixel position or pixel number))
Similarly, the image signal S1 ′, S2 ′ having undergone pixel density conversion is obtained by performing the calculation of Expression (1) also in the sub-scanning direction.
[0141]
On the other hand, in the mask calculation using such a one-dimensional filter, data for performing the mask calculation is insufficient at the edge portions of the digital image signals S1 and S2. For example, in the image signal shown in FIG. 8, when the filtering process is performed on the pixels of S1 (1, 1) and S1 (N, 1) (shaded area), data for three pixels is insufficient. In this case, the virtual pixel positions S1 (-1,1), S1 (-2,1), S1 (-3,1) and S1 (N + 1,1), S1 (N + 2,1), S1 (N + 3, 1), the data value of the pixel positions S1 (1,1) and S1 (N, 1) is copied to this virtual pixel position, and the filtering process is performed assuming that the data value exists at the virtual pixel position. Can be done.
[0142]
The pixel density conversion in the pixel
[0143]
Here, as the interpolation calculation, a high-order interpolation calculation such as a B-spline interpolation calculation focusing on smoothness and a Cubic spline interpolation calculation focusing on sharpness can be applied in addition to linear interpolation.
[0144]
Here, the Cubic spline interpolation calculation and the B spline interpolation calculation will be described. The image signals S1 and S2 used in this embodiment are sampling points (pixels) X arranged in one direction sampled at equal intervals. k-2 , X k-1 , X k , X k + 1 , X k + 2 ,... Corresponding signal values (S k-2 , S k-1 , S k , S k + 1 , S k + 2 , ...). Cubic spline interpolation operation is the original sampling point (pixel) X k ~ X k + 1 Interpolation point X provided between p Interpolation data Y in the cubic Cubic spline interpolation equation (2) representing the interpolation data Y ′ k-1 , Y k , Y k + 1 , Y k + 2 Interpolation coefficient c corresponding to each k-1 , C k , C k + 1 , C k + 2 Are obtained by the operations shown below.
[0145]
Y '= c k-1 Y k-1 + C k Y k + C k + 1 Y k + 1 + C k + 2 Y k + 2 (2)
c k-1 = (-T Three + 2t 2 -T) / 2
c k = (3t Three -5t 2 +2) / 2
c k + 1 = (-3t Three + 4t 2 + T) / 2
c k + 2 = (T Three -T 2 ) / 2
(However, t (0 ≦ t ≦ 1) is set to 1 for the lattice interval and the pixel X k Interpolation point X with reference to p Pixel X k + 1 Indicates the position in the direction. )
The B-spline interpolation operation uses the original sampling point X k ~ X k + 1 Interpolation point X provided between p Interpolation data Y in the cubic B-spline interpolation equation (3) representing the interpolation data Y ′ k-1 , Y k , Y k + 1 , Y k + 2 Interpolation coefficient b respectively corresponding to k-1 , B k , B k + 1 , B k + 2 Are obtained by the operations shown below.
[0146]
Y '= b k-1 Y k-1 + B k Y k + B k + 1 Y k + 1 + B k + 2 Y k + 2 (3)
b k-1 = (-T Three + 3t 2 -3t + 1) / 6
b k = (3t Three -6t 2 +4) / 6
b k + 1 = (-3t Three + 3t 2 + 3t + 1) / 6
b k + 2 = T Three / 6
(However, t (0 ≦ t ≦ 1) is set to 1 for the lattice interval and the pixel X k Interpolation point X with reference to p Pixel X k + 1 Indicates the position in the direction. )
In the present embodiment, the type of interpolation calculation including linear interpolation calculation may be selected according to the values of m and n.
[0147]
The image signals S1 ′ and S2 ′ thus obtained in the pixel density conversion means 37a and 37b are input to the addition means 38. Here, when “low density” is selected, the pixel density is 1/16 compared to when “high density” is selected.
[0148]
The adding means 38 is a filter having a frequency response characteristic for increasing the S / N of the added image signal S3 as described in Japanese Patent Laid-Open No. 7-287330 with respect to the image signals S1 'and S2'. A filtering process is performed, and the image signals S1 ′ and S2 ′ subjected to the filtering process are added together with corresponding pixels to obtain an added image signal S3, which is output to an image processing apparatus or the like.
[0149]
As described above, according to the image information reading apparatus of this embodiment, the main scanning speed of the
[0150]
In particular, when performing double-sided reading as in this embodiment, in order to sufficiently apply energy from the
[0151]
In addition, by performing pixel density conversion processing on the two digital image signals S1 and S2, compared to the case where pixel density conversion is performed on the addition image signal, addition mask processing when obtaining the addition image signal is performed. The amount of calculation to be performed can be reduced, thereby shortening the calculation time of the addition calculation and performing processing at high speed.
[0152]
In the image information reading apparatus according to the present embodiment, high-density reading is set as an initial setting. However, the image information reading method and apparatus according to the present invention are not limited to this mode. “Standard density reading” is set as the setting, and it may be switched to “high density reading” or “low density reading”.
[0153]
For example, when the main scanning frequency (main scanning speed) is 160 Hz and the pixel density of 10 pix / mm is standard density reading, the main scanning sampling interval, sub scanning pixel pitch, sub scanning speed, and anti-aliasing filter cutoff frequency are as follows: Shown in
[0154]
Main scanning sampling interval: 100 μm Sub-scanning pixel pitch: 100 μm
Main scanning sampling period: 1.0 μsec Sub scanning speed: 16 mm / sec
Anti-aliasing filter cutoff frequency: 500 kHz
(Because it is an analog filter, the cut-off frequency is preferably 500 kHz or less, and preferably 400 kHz, for example).
On the other hand, at the time of high-density reading with a pixel density of 20 pix / mm, the main scanning frequency is assumed to be 160 Hz unchanged,
Main scanning sampling interval: 50 μm Sub-scanning pixel pitch: 50 μm
Main scanning sampling period: 0.5 μsec Sub scanning speed: 8 mm / sec
Anti-aliasing filter cutoff frequency: 1000 kHz
On the other hand, in the case of low density reading with a pixel density of 5 pix / mm, the main scanning frequency is assumed to be 160 Hz unchanged,
Main scanning sampling interval: 200 μm Sub-scanning pixel pitch: 200 μm
Main scanning sampling period: 2.0 μsec Sub scanning speed: 32 mm / sec
Anti-aliasing filter cutoff frequency: 250 kHz
It becomes.
[0155]
Note that the reading density such as “high density reading” and “low density reading” as described above is not limited to those set in advance, and parameters (m, n) that define the reading density are arbitrarily selected. It may be possible.
[0156]
In the above embodiment, pixel density conversion processing is performed on the digital image signals S1 and S2. However, the digital image signals S1 and S2 are added first to obtain an added image signal S3, and this added image signal S3. Pixel density conversion processing may be performed on the.
[0157]
Further, in the above embodiment, the pixel density conversion means 37a, 37b performs pixel density conversion so that the digital image signals S1, S2 are 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction. However, the pixel density conversion may be performed so that a / m (a> 0) times in the main scanning direction and a / n times in the sub-scanning direction. In this case, the pixel density of the added image signal S3 is (a / (m × n)) 2 Doubled.
[0158]
FIG. 9 is a diagram showing an embodiment of the third image information reading apparatus of the present invention. In the illustrated image information reading apparatus, a
[0159]
Further, immediately above the
[0160]
Further, a logarithmic amplifier 16a is connected to the photomultiplier 15a, and the analog image signal QA detected by the photomultiplier 15a is logarithmically amplified according to a preset frequency characteristic to obtain a logarithmic image signal QA. ′ Is output.
[0161]
On the other hand, shading correction data D1 corresponding to a preset sampling interval is stored in the
[0162]
Further, an anti-aliasing filter 35a that removes aliasing noise (folding noise) caused by A / D conversion, which will be described later, is connected to the adder 43a, and the filtered image signal Q1 'is connected to the subsequent stage. Is converted to a digital image signal S1 with a preset reference clock. The anti-aliasing filter 35a includes a high-density filter and a low-density filter, and the high-density filter is initially selected. However, the anti-aliasing filter 35a is low depending on the input signal from the characteristic changing means 60 described later. It can also be switched to a density filter.
[0163]
Furthermore, the image information reading apparatus of the present embodiment has a “high density” (10 pix / mm) that is initially set as a pixel density for reading image information, and a lower density than this high density. It receives the density selected by the operator from “Low Density” (5 pix / mm). When “High Density”, the parameter (m, n) is set to m = n = 1, and “Low Density” In some cases, the input means 70 outputs the parameter (m, n) as m = n = 2, and the parameter (m, n) output from the input means 70 is input to these parameters (m, n). Accordingly, the rotational speed of the
[0164]
Next, the operation of the radiation image information reading apparatus of this embodiment will be described.
[0165]
First, a signal representing “high density” selected by the operator is input to the input means 70. “High density” is an initial setting of the image information reading apparatus, and a parameter (m, n) = (1, 1) corresponding to “high density” is output from the input means 70 and input to the characteristic changing means 60. Is done.
[0166]
The characteristic changing means 70, according to the input parameter (m, n) = (1, 1), the rotational speed of the
[0167]
Next, the
[0168]
On the other hand, the
[0169]
From the portion of the
[0170]
The photomultiplier 15a is applied with a high voltage corresponding to the sensitivity for high density by the
[0171]
The analog image signal QA is input to the logarithmic amplifier 16a, converted into the logarithmic image signal QA 'by the logarithmic amplifier 16a set to the high frequency characteristic, and input to the adder 43a.
[0172]
On the other hand, the high-density shading correction data D1 stored in the
[0173]
Here, the clock input to the D / A converter 42a is a reference clock that is switched and output by a
[0174]
The shading correction analog signal D1 'is input to the adder 43a and added to the logarithmic image signal QA' input from the logarithmic amplifier 16a, whereby the logarithmic image signal QA 'is subjected to shading correction. And is input to the anti-aliasing filter 35a.
[0175]
The anti-aliasing filter 35a is switched to a high-density filter by the characteristic changing means 60, and the input image signal Q1 is converted into a signal Q1 'from which aliasing noise has been appropriately removed by this high-density filter. To the A /
[0176]
The A /
[0177]
Next, an operation when “low density” is selected by the operator will be described.
[0178]
When a signal representing “low density” is input to the
[0179]
The characteristic changing means 70, based on the input parameters (m, n) = (2, 2), the rotational speed of the motor that drives the
[0180]
In this way, the characteristics of each part are changed, and image information is read in the same manner as in the case of the high-density reading described above.
[0181]
That is, the
[0182]
The stimulated emission light 13a emitted from the portion of the
[0183]
The photomultiplier 15a is applied with a voltage that has a lower sensitivity than that for high-density reading by the
[0184]
The analog image signal QA is input to the logarithmic amplifier 16a, converted into the logarithmic image signal QA 'by the logarithmic amplifier 16a set to the frequency characteristic for low density, and input to the adder 43a.
[0185]
On the other hand, the high-density shading correction data D1 stored in the
[0186]
The shading correction analog signal D1 'is input to the adder 43a and added to the logarithmic image signal QA' input from the logarithmic amplifier 16a, whereby the logarithmic image signal QA 'is subjected to shading correction. And is input to the anti-aliasing filter 35a.
[0187]
The anti-aliasing filter 35a is switched to a low density filter by the characteristic changing means 60, and the input image signal Q1 is converted into a signal Q1 'from which aliasing noise has been appropriately removed by the low density filter. To the A /
[0188]
The A /
[0189]
As described above, according to the image information reading apparatus of the present embodiment, the sub-scanning is 1/2 in the main scanning direction of the laser beam L as compared with the case of high-density reading without changing the main scanning speed of the laser beam. A low-density image signal having a pixel density of 1/4 as a whole, which is 1/2 the number of pixels in the direction, can be obtained. In addition, the characteristic changing
[0190]
In the image information reading apparatus according to the present embodiment, high-density reading is set as an initial setting. However, the image information reading method and apparatus according to the present invention are not limited to this mode. “Standard density reading” is set as the setting, and it may be switched to “high density reading” or “low density reading”.
[0191]
For example, when the main scanning frequency (main scanning speed) is 160 Hz and the pixel density is 10 pix / mm, the main scanning sampling interval, the sub scanning pixel pitch, the sub scanning speed, and the anti-aliasing filter cutoff frequency are as follows: Shown in
[0192]
Main scanning sampling interval: 100 μm Sub-scanning pixel pitch: 100 μm
Main scanning sampling period: 1.0 μsec Sub scanning speed: 16 mm / sec
Anti-aliasing filter cutoff frequency: 500 kHz
(Because it is an analog filter, it is preferably 500 kHz or less, and preferably 400 kHz or the like is applied.)
On the other hand, at the time of high-density reading with a pixel density of 20 pix / mm, the main scanning frequency is assumed to be 160 Hz unchanged,
Main scanning sampling interval: 50 μm Sub-scanning pixel pitch: 50 μm
Main scanning sampling period: 0.5 μsec Sub scanning speed: 8 mm / sec
Anti-aliasing filter cutoff frequency: 1000 kHz
On the other hand, in the case of low density reading with a pixel density of 5 pix / mm, the main scanning frequency is assumed to be 160 Hz unchanged,
Main scanning sampling interval: 200 μm Sub-scanning pixel pitch: 200 μm
Main scanning sampling period: 2.0 μsec Sub scanning speed: 32 mm / sec
Anti-aliasing filter cutoff frequency: 250 kHz
It becomes.
[0193]
Note that the reading density such as “high density reading” and “low density reading” as described above is not limited to those set in advance, and parameters (m, n) that define the reading density are arbitrarily selected. It may be possible.
[Brief description of the drawings]
FIG. 1 is a diagram showing a configuration of an embodiment of a first image information reading apparatus of the present invention.
FIG. 2 is a diagram showing an embodiment of a configuration for changing a clock;
3 is a diagram showing a configuration of characteristic changing means in the embodiment shown in FIG. 1;
FIG. 4 is a diagram for explaining details of memory management in the embodiment shown in FIG. 1;
FIG. 5 is a diagram showing an embodiment of a second image information reading apparatus of the present invention.
6 is a diagram showing a configuration of characteristic changing means in the embodiment shown in FIG. 5;
FIG. 7 is a diagram showing an embodiment of a configuration for changing a clock;
FIG. 8 is a diagram for explaining filtering processing;
FIG. 9 is a diagram showing an embodiment of a third image information reading apparatus of the present invention.
[Explanation of symbols]
1 Storage phosphor sheet
8 Motor
9a, 9b Endless belt
10 Laser light source
11 Laser light
12 rotating polygon mirror
14a, 14b Light guide
15a, 15b Photo multiplayer
16a, 16b logarithmic amplifier
17a, 17b Laser light cut filter
20 Motor
21 Scanning lens
35a, 35b Anti-aliasing filter
36a, 36b A / D converter
37a, 37b Pixel density conversion means
38 Adding means
39a, 39b Voltage applying means
41a, 41b memory
42a, 42b A / D converter
43a, 43b Adder
60 Characteristic changing means
70 Input means
Claims (14)
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍するとともに、
前記m,nに応じて、以下の(1)から(4)の特性のうち少なくとも1つを変更することを特徴とする画像情報読取方法。
(1)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(2)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(3)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(4)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数The light beam is sub-scanned in a direction substantially perpendicular to the main scanning direction while repeatedly scanning the light beam at a constant speed in a constant direction on the scanning sheet on which the image information is recorded. A digital image having a predetermined pixel density is obtained by scanning almost the entire surface of the scanning sheet, photoelectrically detecting signal light emitted from the scanned sheet by scanning the light beam, and sampling and quantizing the obtained analog image signal. In an image information reading method for obtaining a signal,
The sub-scan speed is multiplied by m (m> 0), the sampling interval is multiplied by n (n> 0),
An image information reading method, wherein at least one of the following characteristics (1) to ( 4 ) is changed in accordance with m and n.
( 1 ) Preset shading correction data when performing shading correction on the analog image signal ( 2 ) Timing for outputting the shading correction data from the memory storing the shading correction data ( 3 ( 4 ) Frequency transfer characteristics when the analog image signal is logarithmically amplified ( 4 ) Cut-off frequency when filtering to remove aliasing noise is performed prior to sampling the analog image signal
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍し、さらに前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n倍とする画素密度変換処理を行なうとともに、
前記m,nに応じて、前記画素密度変換処理におけるフィルタリング処理のパラメータを変更することを特徴とする画像情報読取方法。The light beam is sub-scanned in a direction substantially perpendicular to the main scanning direction while repeatedly scanning the light beam at a constant speed in a constant direction on the scanning sheet on which the image information is recorded. A digital image having a predetermined pixel density is obtained by scanning almost the entire surface of the scanning sheet, photoelectrically detecting signal light emitted from the scanned sheet by scanning the light beam, and sampling and quantizing the obtained analog image signal. In an image information reading method for obtaining a signal,
The sub-scanning speed is multiplied by m (m> 0), the sampling interval is multiplied by n (n> 0), and the number of pixels in the main scanning direction is a / m (a> 0) times. While performing pixel density conversion processing to increase the number of pixels in the scanning direction by a / n times,
An image information reading method, wherein a parameter of a filtering process in the pixel density conversion process is changed according to m and n.
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍するととともに、
前記m,nに応じて、以下の(1)から(4)の特性のうち少なくとも1つを変更することを特徴とする画像情報読取方法。
(1)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(2)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(3)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(4)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数The light beam is sub-scanned in a direction substantially perpendicular to the main scanning direction while repeatedly scanning the light beam at a constant speed in a constant direction on the scanning sheet on which the image information is recorded. The scanning sheet is scanned over substantially the entire surface, the signal light emitted from both sides of the scanned sheet is photoelectrically detected by scanning the light beam, and the two analog image signals obtained are sampled and quantized to obtain two digital signals. In an image information reading method for obtaining an image signal and adding the two digital image signals to obtain an added image signal having a predetermined pixel density,
Multiplying the sub-scan speed by m (m> 0) and multiplying the sampling interval by n (n>0);
An image information reading method, wherein at least one of the following characteristics (1) to ( 4 ) is changed in accordance with m and n.
( 1 ) Preset shading correction data when performing shading correction on the analog image signal ( 2 ) Timing for outputting the shading correction data from the memory storing the shading correction data ( 3 ( 4 ) Frequency transfer characteristics when the analog image signal is logarithmically amplified ( 4 ) Cut-off frequency when filtering to remove aliasing noise is performed prior to sampling the analog image signal
前記副走査の速度をm(m>0)倍し、前記サンプリングの間隔をn(n>0)倍し、さらに前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n倍とする画素密度変換処理を行なうとともに、
前記m,nに応じて、前記画素密度変換処理におけるフィルタリング処理のパラメータを変更することを特徴とする画像情報読取方法。The light beam is sub-scanned in a direction substantially perpendicular to the main scanning direction while repeatedly scanning the light beam at a constant speed in a constant direction on the scanning sheet on which the image information is recorded. The scanning sheet is scanned over substantially the entire surface, the signal light emitted from both sides of the scanned sheet is photoelectrically detected by scanning the light beam, and the two analog image signals obtained are sampled and quantized to obtain two digital signals. In an image information reading method for obtaining an image signal and adding the two digital image signals to obtain an added image signal having a predetermined pixel density,
The sub-scanning speed is multiplied by m (m> 0), the sampling interval is multiplied by n (n> 0), and the number of pixels in the main scanning direction is a / m (a> 0) times. While performing pixel density conversion processing to increase the number of pixels in the scanning direction by a / n times,
An image information reading method, wherein a parameter of a filtering process in the pixel density conversion process is changed according to m and n.
前記副走査手段による前記副走査の速度を、入力されたm(m>0)に応じてm倍に変更せしめる副走査速度変更手段と、
前記A/D変換手段による前記サンプリングの間隔を、入力されたn(n>0)に応じてn倍に変更せしめるサンプリング間隔変更手段と、
前記m,nに応じて、以下の(1)から(4)の特性のうち少なくとも1つを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。
(1)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(2)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(3)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(4)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in the direction, and signal light emitted from the scanned sheet is detected photoelectrically by scanning the light beam. In an image information reading apparatus comprising photoelectric detection means and A / D conversion means for sampling and quantizing a detected analog image signal to obtain a digital image signal having a predetermined pixel density,
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m times according to the input m (m>0);
Sampling interval changing means for changing the sampling interval by the A / D conversion means to n times according to the input n (n> 0) ;
An image information reading apparatus comprising: a characteristic changing unit that changes at least one of the following characteristics (1) to ( 4 ) according to m and n.
( 1 ) Preset shading correction data when performing shading correction on the analog image signal ( 2 ) Timing for outputting the shading correction data from the memory storing the shading correction data ( 3 ( 4 ) Frequency transfer characteristics when the analog image signal is logarithmically amplified ( 4 ) Cut-off frequency when filtering to remove aliasing noise is performed prior to sampling the analog image signal
前記副走査手段による前記副走査の速度を、入力されたm(m>0)に応じてm倍に変更せしめる副走査速度変更手段と、
前記A/D変換手段による前記サンプリングの間隔を、入力されたn(n>0)に応じてn倍に変更せしめるサンプリング間隔変更手段と、
前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n(a>0)倍とする画素密度変換処理を行う画素密度変換処理手段と、
前記m,nに応じて、前記画素密度変換処理におけるフィルタリング処理のパラメータを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in the direction, and signal light emitted from the scanned sheet is detected photoelectrically by scanning the light beam. In an image information reading apparatus comprising photoelectric detection means and A / D conversion means for sampling and quantizing a detected analog image signal to obtain a digital image signal having a predetermined pixel density,
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m times according to the input m (m>0);
Sampling interval changing means for changing the sampling interval by the A / D conversion means to n times according to the input n (n>0);
Pixel density conversion processing means for performing a pixel density conversion process in which the number of pixels in the main scanning direction is a / m (a> 0) times and the number of pixels in the sub-scanning direction is a / n (a> 0) times;
An image information reading apparatus comprising: characteristic changing means for changing a parameter of filtering processing in the pixel density conversion processing according to m and n.
前記副走査手段による前記副走査の速度を、入力されたmに応じてm(m>0)倍に変更せしめる副走査速度変更手段と、
前記A/D変換手段による前記サンプリングの間隔を、入力されたnに応じてn(n>0)倍に変更せしめるサンプリング間隔変更手段と、
前記m,nに応じて、以下の(1)から(4)の特性のうち少なくとも1つを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。
(1)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(2)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(3)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(4)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in a direction, and signal light emitted from both sides of the scanned sheet by the scanning of the light beam photoelectrically. Photoelectric detection means for detecting; A / D conversion means for obtaining two digital image signals by sampling and quantizing two detected analog image signals; and adding the two digital image signals to obtain a predetermined pixel density In an image information reading apparatus provided with an adding means for obtaining an added image signal,
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m (m> 0) times according to the input m;
Sampling interval changing means for changing the sampling interval by the A / D conversion means to n (n> 0) times according to the input n;
An image information reading apparatus comprising: a characteristic changing unit that changes at least one of the following characteristics (1) to ( 4 ) according to m and n.
( 1 ) Preset shading correction data when performing shading correction on the analog image signal ( 2 ) Timing for outputting the shading correction data from the memory storing the shading correction data ( 3 ( 4 ) Frequency transfer characteristics when the analog image signal is logarithmically amplified ( 4 ) Cut-off frequency when filtering to remove aliasing noise is performed prior to sampling the analog image signal
前記副走査手段による前記副走査の速度を、入力されたmに応じてm(m>0)倍に変更せしめる副走査速度変更手段と、
前記A/D変換手段による前記サンプリングの間隔を、入力されたnに応じてn(n>0)倍に変更せしめるサンプリング間隔変更手段と、
前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n(a>0)倍とする画素密度変換処理を行う画素密度変換処理手段と、
前記m,nに応じて、前記画素密度変換処理におけるフィルタリング処理のパラメータを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in a direction, and signal light emitted from both sides of the scanned sheet by the scanning of the light beam photoelectrically. Photoelectric detection means for detecting; A / D conversion means for obtaining two digital image signals by sampling and quantizing two detected analog image signals; and adding the two digital image signals to obtain a predetermined pixel density In an image information reading apparatus provided with an adding means for obtaining an added image signal,
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m (m> 0) times according to the input m;
Sampling interval changing means for changing the sampling interval by the A / D conversion means to n (n> 0) times according to the input n;
Pixel density conversion processing means for performing a pixel density conversion process in which the number of pixels in the main scanning direction is a / m (a> 0) times and the number of pixels in the sub-scanning direction is a / n (a> 0) times;
An image information reading apparatus comprising: characteristic changing means for changing a parameter of filtering processing in the pixel density conversion processing according to m and n.
前記所定の画素密度に対して、主走査方向について1/m倍、副走査方向について1/n倍となる画素密度に対応した値m、nの入力を受ける画素密度入力手段と、
前記副走査手段による前記副走査の速度をm(m>0)倍に変更せしめる副走査速度変更手段と、
前記A/D変換器による前記サンプリングの間隔をn(n>0)倍に変更せしめるサンプリング間隔変更手段と、
前記m,nに応じて、以下の(1)から(4)の特性のうち少なくとも1つを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。
(1)前記アナログ画像信号に対してシェーディング補正を行う場合における、予め設定されているシェーディング補正用データ
(2)前記シェーディング補正用データを記憶したメモリから該シェーディング補正用データを出力するタイミング
(3)前記アナログ画像信号を対数増幅する場合における周波数伝達特性
(4)前記アナログ画像信号をサンプリングするのに先立ってエリアジングノイズを除去するフィルタリングを行う場合におけるカットオフ周波数A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in the direction, and signal light emitted from the scanned sheet is detected photoelectrically by scanning the light beam. In an image information reading apparatus comprising photoelectric detection means and A / D conversion means for sampling and quantizing a detected analog image signal to obtain a digital image signal having a predetermined pixel density,
Pixel density input means for receiving inputs of values m and n corresponding to a pixel density of 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction with respect to the predetermined pixel density;
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m (m> 0) times;
Sampling interval changing means for changing the sampling interval by the A / D converter to n (n> 0) times;
An image information reading apparatus comprising: a characteristic changing unit that changes at least one of the following characteristics (1) to ( 4 ) according to m and n.
( 1 ) Preset shading correction data when performing shading correction on the analog image signal ( 2 ) Timing for outputting the shading correction data from the memory storing the shading correction data ( 3 ( 4 ) Frequency transfer characteristics when the analog image signal is logarithmically amplified ( 4 ) Cut-off frequency when filtering to remove aliasing noise is performed prior to sampling the analog image signal
前記所定の画素密度に対して、主走査方向について1/m倍、副走査方向について1/n倍となる画素密度に対応した値m、nの入力を受ける画素密度入力手段と、
前記副走査手段による前記副走査の速度をm(m>0)倍に変更せしめる副走査速度変更手段と、
前記A/D変換器による前記サンプリングの間隔をn(n>0)倍に変更せしめるサンプリング間隔変更手段と、
前記主走査方向の画素数をa/m(a>0)倍、前記副走査方向の画素数をa/n(a>0)倍とする画素密度変換処理を行う画素密度変換処理手段と、
前記m,nに応じて、前記画素密度変換処理におけるフィルタリング処理のパラメータを変更せしめる特性変更手段とを備えたことを特徴とする画像情報読取装置。A main scanning means for repeatedly scanning a light beam emitted from a light source on a scanned sheet on which image information is recorded at a constant speed in a constant direction, and substantially orthogonal to the main scanning direction of the light beam. Sub-scanning means for sub-scanning the sheet or the light beam so that the sheet is relatively sub-scanned in the direction, and signal light emitted from the scanned sheet is detected photoelectrically by scanning the light beam. In an image information reading apparatus comprising photoelectric detection means and A / D conversion means for sampling and quantizing a detected analog image signal to obtain a digital image signal having a predetermined pixel density,
Pixel density input means for receiving inputs of values m and n corresponding to a pixel density of 1 / m times in the main scanning direction and 1 / n times in the sub scanning direction with respect to the predetermined pixel density;
Sub-scanning speed changing means for changing the sub-scanning speed by the sub-scanning means to m (m> 0) times;
Sampling interval changing means for changing the sampling interval by the A / D converter to n (n> 0) times;
Pixel density conversion processing means for performing a pixel density conversion process in which the number of pixels in the main scanning direction is a / m (a> 0) times and the number of pixels in the sub-scanning direction is a / n (a> 0) times;
An image information reading apparatus comprising: characteristic changing means for changing a parameter of filtering processing in the pixel density conversion processing according to m and n.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000082548A JP4138199B2 (en) | 1999-03-23 | 2000-03-23 | Image information reading method and apparatus |
Applications Claiming Priority (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7855799 | 1999-03-23 | ||
| JP26831399 | 1999-09-22 | ||
| JP2000-45953 | 2000-02-23 | ||
| JP2000045953 | 2000-02-23 | ||
| JP11-78557 | 2000-02-23 | ||
| JP11-268313 | 2000-02-23 | ||
| JP2000082548A JP4138199B2 (en) | 1999-03-23 | 2000-03-23 | Image information reading method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001313815A JP2001313815A (en) | 2001-11-09 |
| JP4138199B2 true JP4138199B2 (en) | 2008-08-20 |
Family
ID=27466183
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000082548A Expired - Fee Related JP4138199B2 (en) | 1999-03-23 | 2000-03-23 | Image information reading method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4138199B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4715044B2 (en) * | 2001-06-28 | 2011-07-06 | ソニー株式会社 | Additional information decoding apparatus for television signal |
| JP5921511B2 (en) | 2013-09-30 | 2016-05-24 | 富士フイルム株式会社 | Radiation image reading apparatus, radiation image reading program, and radiation image reading method |
| JP6684145B2 (en) * | 2016-04-28 | 2020-04-22 | キヤノン株式会社 | Image forming device |
-
2000
- 2000-03-23 JP JP2000082548A patent/JP4138199B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001313815A (en) | 2001-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6580525B1 (en) | Method of and system reading out image signal | |
| JPH052031B2 (en) | ||
| JP3394375B2 (en) | Radiation image reader | |
| JP2754068B2 (en) | Radiation image signal generation method and radiation image reading device | |
| JP4138199B2 (en) | Image information reading method and apparatus | |
| JP3416052B2 (en) | How to get shading correction data | |
| JP2952418B2 (en) | Radiation image signal generation method | |
| JP2651625B2 (en) | Energy subtraction image generation method | |
| US20020024027A1 (en) | Image signal generating method, apparatus and program | |
| JPH02280099A (en) | Accumulation type fluorescent material sheet and energy subtraction method | |
| JPH0617984B2 (en) | Radiation image information reading method | |
| JP2707369B2 (en) | Image display method | |
| JPH0950516A (en) | Method and device for interpolation operation of picture data | |
| JPH09160149A (en) | Image processing method and device | |
| JP4163401B2 (en) | Image signal generation method, apparatus, and program | |
| US7366358B2 (en) | Method of and system for generating image signal | |
| JPH08116435A (en) | Shading correction method in radiation image reading | |
| JP2003319260A (en) | Image signal generating method and apparatus | |
| JP2707346B2 (en) | Radiation image information reader | |
| JP2867096B2 (en) | Image smoothing method | |
| JPH10232452A (en) | Method and device for reading radiograph information | |
| JPH11341221A (en) | Method and device for image reading | |
| JP3738851B2 (en) | Shading correction method in radiographic image reading | |
| JPH0968767A (en) | Method and device for reading image | |
| JPS62227272A (en) | Method and device for reducing noise of picture signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050912 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061201 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070626 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070710 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070910 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080304 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080425 |
|
| A911 | Transfer of reconsideration by examiner before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20080519 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080603 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080605 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120613 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120613 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130613 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |