(57)【特許請求の範囲】
【請求項１】
受光する光量に対して出力が直線性をもたないセンサを用い、被測定物に光を照射しその被検出部からの光を前記センサで受光して測定する際に、前記センサの受光量を変化させたときの前記センサからの出力が受光量に比例するように前記センサの出力を補正する直線化処理工程を備えているセンサ出力の補正方法において、
前記直線化処理工程は次の工程（Ａ）と（Ｂ）を含んでいることを特徴とするセンサ出力の補正方法。
（Ａ）受光する光量に対して出力が直線性をもつ光検出器を別途設け、前記センサに入射する光をその光検出器にも同時に入射させるとともに、その入射光量を変化させたときの前記センサの出力と前記光検出器の出力との関係を直線化データとして保持しておく工程、及び
（Ｂ）被測定物を測定したときの前記センサの出力を前記直線化データに基づいて前記光検出器の出力に比例するように補正する工程。
【請求項２】
前記センサがエリアセンサであり、前記直線化処理工程をピクセルごとに行なう請求の範囲第１項に記載のセンサ出力の補正方法。
【請求項３】
前記センサがエリアセンサであり、前記直線化処理工程を画像内で最も明るいピクセル付近の幾つかのピクセルを選択して、それらのピクセルの出力の平均値を用いて行なう請求の範囲第１項に記載のセンサ出力の補正方法。
【請求項４】
前記エリアセンサとしてＣＣＤ型又はＣＭＯＳ型センサを用いる請求の範囲第２項又は第３項に記載のセンサ出力の補正方法。
【請求項５】
前記センサの出力は、受光量がゼロのときの出力をダークデータとして差し引いたオフセット処理後の値である請求項第１項から第４項のいずれかに記載のセンサ出力の補正方法。
【請求項６】
被測定物として基準物を測定し、そのときの前記エリアセンサの各ピクセルの出力を前記直線化処理により補正したものが均一になるように各ピクセルの出力を補正する光むら補正処理をさらに含んでいる請求の範囲第２項から第５項のいずれかに記載のセンサ出力の補正方法。
【請求項７】
被測定物として基準物を測定し、前記光むら補正処理はピクセルが飽和光量に達したときの受光量に対する飽和光量に近い一定割合の受光量での画像データに対して行なう請求の範囲第６項に記載のセンサ出力の補正方法。
【請求項８】
前記基準物は面内の濃淡が均一な反射板又はブランクである請求の範囲第６項又は第７項に記載のセンサ出力の補正方法。
【請求項９】
受光する光量に対して出力が直線性をもつ前記光検出器としてフォトダイオードを用いる請求の範囲第１項に記載のセンサ出力の補正方法
【発明の詳細な説明】
技術分野
本発明は被測定物を光照射しその被検出部からの光をセンサで受光して測定する方法において、そこで用いられるセンサが受光量に対して出力の直線性をもたないものである場合におけるそのセンサ出力の補正方法に関するものである。
本発明で測定対象とする光は、反射光、透過光、蛍光、燐光、化学発光など、定量測定や定性測定に使用されるあらゆる光を含んでいる。
背景技術
光検出器にはフォトダイオードのような単素子の光検出素子、フォトダイオードアレイのような光検出素子がライン上に配列されたリニアセンサ、及びＣＣＤセンサやＣＭＯＳセンサのように光検出素子が二次元に配列されたエリアセンサがある。
受光量に対するセンサの出力特性として、フォトダイオードのような直線性をもったものだけでなく、ＣＣＤセンサやＣＭＯＳセンサなどエリアセンサのように受光量の小さい領域と大きい領域で直線から外れた感度特性をもっているものがある。直線性をもたないセンサは定量には適さないと考えられてきた。
そこで、本発明は出力特性が直線性をもたないセンサを検出器に用いた場合に起こる直線性の問題を解決することにより、そのようなセンサを用いても測定を可能にすることを目的とするものである。
発明の開示
本発明は、被測定物に光を照射しその被検出部からの光をセンサで受光して測定する方法において、そのセンサとして受光量に対して出力が直線性をもたないものを使用する。そして、センサの受光量を変化させたときのセンサからの出力が受光量に比例するようにセンサの出力を補正する直線化処理工程を備えていることを特徴とするものである。
第１の局面では、直線化処理工程は次の工程（Ａ）と（Ｂ）を含んでいる。
（Ａ）受光量に対して出力が直線性をもつ光検出器を別途設け、センサに入射する光をその光検出器にも同時に入射させるとともに、その入射光量を変化させたときのセンサの出力とその光検出器の出力との関係を直線化データとして保持しておく工程、及び
（Ｂ）被測定物を測定したときのセンサの出力をその直線化データに基づいて光検出器の出力に比例するように補正する工程。
センサの出力は受光量に対して直線性をもっていないが、直線化データに基づいて、受光量に対して出力が直線性をもつ光検出器の出力に比例するように補正することにより、補正されたセンサ出力は受光量に対して直線性をもつようになる。
第２の局面では、直線化処理工程は次の工程（Ａ）と（Ｂ）を含んでいる。
（Ａ）発生する光が既知で互いに異なる複数の標準板を用意し、それらの標準板を測定したときのセンサの出力と標準板からの光との関係を直線化データとして保持しておく工程、
（Ｂ）被測定物を測定したときのセンサの出力を直線化データに基づいて標準板からの光に比例するように補正する工程。
この場合もセンサの出力は受光量に対して直線性をもっていないが、直線化データに基づいて、標準板からの光に比例するように補正することにより、補正されたセンサ出力は受光量に対して直線性をもつようになる。
第２の局面では、光量モニタ用に受光量に対して出力が直線性をもつ光検出器とその光検出器の計測器を別途設ける必要がない。この方法は、光源を光量調整した後、標準板を交換して測定するだけであるので、作業が簡単である。実際にセンサが受光したデータで補正するので、誤差要因を軽減することができる。
第３の局面では、直線化処理工程は次の工程（Ａ）と（Ｂ）を含んでいる。
（Ａ）センサは露光時間の設定が可変なものであり、１つの基準物を測定する露光時間を複数段階に異ならせて測定したときのセンサの出力と露光時間との関係を、露光時間に比例した基準物からの光に関する直線化データとして保持しておく工程、及び
（Ｂ）被測定物を測定したときのセンサの出力を直線化データに基づいて露光時間から導かれる基準物からの光に比例するように補正する工程。
基準物として、例えば面内の濃淡が均一な反射板又はブランク（被測定物を置かないで測定光を全てイメージセンサに入射させる状態）を用いる。
第３の局面では、センサの受光量は基準物を露光する露光時間に比例している。この場合もセンサの出力は受光量に対して直線性をもっていないが、直線化データに基づいて、露光時間から導かれる基準物からの光に比例するように補正することにより、補正されたセンサ出力は受光量に対して直線性をもつようになる。
第３の局面では、第２の局面で得られる利点の他に、複数の標準板を必要とせず、基準となる１枚の基準物のみで直線化データを構築することができるので、作業がより簡便になる。
この局面では、基準物として例えば基準白板を装置に内蔵しておけば、測定ごと、又は適当な間隔で直線化データを得ることが容易になる。また、自動で直線化データを得ることも容易であるため、測定精度を維持する上で好都合である。
本発明の対象とするセンサの１つの例はエリアセンサである。その場合、上で述べた直線化処理工程はピクセル（画素）ごとに行なうことができる。また、その直線化処理工程を画像内で最も明るいピクセル付近の幾つかのピクセルを選択して、それらのピクセルの出力の平均値を用いて行なうこともできる。
エリアセンサの例は、ＣＣＤ型又はＣＭＯＳ型センサである。
センサの出力は、受光量がゼロのときの出力をダークデータとして差し引いたオフセット処理後の値であることが好ましい。
センサを用いる分析計の１つとして反射率測定装置がある。反射率測定装置のセンサとしては、精度の良さ、コストパフォーマンス、技術的難易度の観点から主にフォトダイオードが用いられている。しかし、フォトダイオードを用いて複数項目の反射率を得ようとした場合、光学系あるいは試験片を移動させる必要がある。また、フォトダイオードはスポット径内の平均化されたデータを得るものであるため、斑発色検知に代表されるような発色具合の詳細検知を行なう場合には不向きである。
これらの問題を解消する一つの手段として、エリアセンサの採用が考えられる。エリアセンサのデータは対象領域の画像情報であるため、１フレームの情報から複数の項目の測定、斑発色の検知、試験片の位置ずれの補正等を行なうことができる。
エリアセンサによって対象物の形状や色彩を検出することは広く知られている。例えば、免疫測定用のテストストリップの画像をＣＣＤカメラで取り込み、画像の面積又は縦と横の長さの比率に基づいて判定を行なう方法が提案されている（特開平９−２５７７０８号公報参照）。そこでは取り込んだ信号を輝度信号として二値化した後に画像の形状を測定しており、その画像内の濃淡を測定するものではない。
エリアセンサを用いて二次元的な測定を行なっている他の例としては、尿測定機がある。そこでは、尿試験紙の発色の濃淡（明度）ではなく、色調（色相）を判定することによって測定するのが一般的であり、カラーＣＣＤが用いられている。
エリアセンサで被測定物の二次元的な濃淡分布を精度よく検出しようとすると、光の照射むら、レンズの収差などのほか、エリアセンサのピクセル間の感度差により面内の光むらが発生する。そのため、精度よく検出を行なうために、機械的な駆動系を用いてセンサ又は被測定物を移動させるのが一般的である。その場合、エリアセンサを用いていても、一次元のリニアセンサとして利用されているにすぎない。
エリアセンサを用いて、被測定物からの光をもとにした定量測定を行おうとすると、非直線性に基づく上のような問題の他に、次のような問題が生じる。すなわち、二次元的な測定を行おうとすると、光の照射むら、レンズの収差、場所によるピクセルの感度のバラツキなどに基づいて面内での光むらが発生するため、二次元的な測定を行なうと定量結果に場所的なバラツキが生じる。
そこで、エリアセンサを用いて、機械的な駆動系を必要としない簡便な二次元の測定方法を実現するためには、被測定物として基準物を測定し、そのときのエリアセンサの各ピクセルの出力を前述の直線化処理により補正したものが均一になるように各ピクセルの出力を補正する光むら補正処理をさらに含んでいることが好ましい。
エリアセンサを用いた測定装置で撮像した生の画像情報は、エリアセンサの各ピクセル感度の個体差、光源（例えばＬＥＤ）の照射むら、レンズのコサイン４乗則（収差）などの影響をうける。「光むら」はこれらの全ての影響を受けて生じたものである。
直線化処理によりエリアセンサからの出力が直線性をもったものになるので、さらに光むら補正処理により面内での光むらがなくなれば、機械的な駆動系を用いなくても二次元の測定を精度よく行なうことができるようになる。
光むら補正処理の１つの方法は、被測定物として基準物を測定し、ピクセルの受光量が飽和光量に達したときの光量に対する一定割合の受光量での画像データに対して行なう方法である。その割合は、例えば０．８というような比較的飽和光量に近い値にすることにより、精度よく光むら補正処理を行なうことができるようになる。
エリアセンサとしては、ＣＣＤ（電荷結合デバイス）型センサ又はＣＭＯＳ型センサを用いることができる。
本発明により直線化処理を施すことにより、受光量に対して出力が直線性をもたないセンサを用いても、定量測定が可能になる。
さらに光むら補正処理まで施せば、エリアセンサを用いて測定を実施すれば、次のような効果を達成することができる。
（１）尿試験などを測定する際に、同じ色相で濃淡が変化する項目についても測定することができる。
（２）機械的な駆動型を必要とする装置に比べて全体画像を取り込む速度が早く、高速測定が可能になる。
（３）機械的な駆動型を必要としないので、安価に実現できる。
【図面の簡単な説明】
第１図は本発明が適用される反射率測定装置の一例を一部ブロック図で示す概略構成図である。第２図は同反射率測定装置における光学系の具体例を示したものであり、（ａ）は光学系の外観図、（ｂ）はその垂直断面図、（ｃ）は（ｂ）における円内の拡大図である。第３図はエリアセンサの出力特性を示す図である。第４図はエリアセンサにて白板を撮像したときの３次元等高面グラフである。第５図は本発明の補正処理の一実施例を示すフローチャート図である。第６図はデータ取得を行なうピクセルの位置を示す平面図である。第７図は補正の一例を示す図であり、（ａ）は白板を被測定物としたときのエリアセンサの各ピクセル出力のグラフ、（ｂ）はＰＤ電圧値変換した後の出力のグラフ、（ｃ）はさらに（ｂ）に対して光むら補正したグラフである。第８図はＬＥＤ光量の暗いものから明るいものまで３段階に分けて白板を撮像し、その画情報を並べて３次元等高面グラフにしたものであり、（ａ）は光むら補正前の状態、（ｂ）は（ａ）のグラフに対して光むら補正を行なったものである。第９図は１枚の画像を１１×１１の領域に分割する例を示す平面図である。第１０図は１枚の画像を１１×１１の領域に分割する例を示す平面図である。第６図は光むら補正を行なうピクセルの位置を示す平面図である。第１１図は第６図で示す各画像の５点を光むら補正してプロットした図である。第１２図は実施例の反射率測定装置の主なユニットの時間経過と温度の関係を示す図である。第１３図は１０分おきに３回ずつ測定した反射率の結果をプロットした図である。第１４図は本発明が適用される反射率測定装置の他の例を一部ブロック図で示す概略構成図である。第１５図は第２の実施例において反射率の異なる標準板を測定した結果をエリアセンサの出力とともに示したものである。第１６図は第１５図の結果をエリアセンサの１つのピクセルについて示したものである。第１７図は同実施例で直線化データを得るための手順を示したフローチャート図である。第１８図は未知試料の反射率測定手順を示したフローチャート図である。第１９図は第３の実施例において露光時間を減少させていったときのエリアセンサの出力と露光時間の関係を示した図である。第２０図は第１９図の結果をエリアセンサの１つのピクセルについて示したものである。第２１図は同実施例で直線化データを得るための手順を示したフローチャート図である。第２２図は未知試料の反射率測定手順を示したフローチャート図である。第２３図は同実施例において補正後のデータの精度を確認するためのデータ取得ピクセル位置を示す平面図である。第２４図はそれぞれのピクセルにおける補正された出力と反射率の関係を示すグラフである。
発明を実施するための最良の形態
以下では反射率を測定する場合を例として説明するが、本発明は反射率に代えて透過率、蛍光、燐光、化学発光などを測定する場合には同様に適用することができる。
［実施例１］
第１の実施例として、センサとしてエリアセンサを使用し、本発明の第１の局面による出力補正方法を適用した二次元の反射率測定装置の一例の概要を第１図に示す。
２は被測定物であり、試料台（図示略）に保持されて所定の位置に載置される。臨床検査などの実際の測定では、被測定物２は尿試験紙や免疫測定用試験紙などの試験紙、化学分析における薄層クロマトなどであるが、エリアセンサを補正する場合には表面での反射率が均一な白板が該当する。被測定物２を照射するために、光源として３つのＬＥＤ（発光ダイオード）４が被測定物２の周囲上方に互いに１２０度間隔で同じ高さに配置され、被測定物２の中心に向って４５度の入射角で被測定物２を照射する。ＬＥＤ４はいずれも発光の中心波長が６３５ｎｍのものである。
被測定物２の上方には結像レンズ６を介してＣＭＯＳエリアセンサ８が配置されている。被測定物２の反射光がそのレンズ６によりエリアセンサ８に結像されることにより、被測定物２の画情報がエリアセンサ８により検出される。
ＬＥＤ４の光量を感知でき、かつエリアセンサ８の画角外の場所に光検出器（ＰＤ）１０が配置されている。光検出器１０はフォトダイオードであり、受光する光量に対して出力が直線性をもっており、被測定物２への照射光量を電圧に変換する。１２はその光検出器１０が受光した光量を電圧に変換する電圧計である。
破線のブロック１４は、ＬＥＤ４、レンズ６、エリアセンサ８及び光検出器１０がこの反射率測定装置の光学系を構成していることを示している。
破線のブロック２０はエリアセンサドライブ回路であり、エリアセンサ８の出力を増幅する増幅器２２、増幅されたアナログ出力をデジタル信号に変換するＡ／Ｄ変換器２４及び取り込んだデジタル信号を一時的に保持するＲＡＭ（ランダム・アクセス・メモリ）２６を備えている。エリアセンサドライブ回路２０は、撮像時間などのレジスタ設定、画像データ取込み等、エリアセンサ８を制御するものである。また、エリアセンサドライブ回路２０は、ＬＥＤ４の光量を調節したり、パソコン（パーソナルコンピュータ）２８とシリアル通信（５６０００ｂｐｓ）を行なったり、パソコン２８からの命令を実行したりする。
パソコン２８は、エリアセンサ８の各種レジスタ設定を行なったり、エリアセンサドライブ回路２０に命令を与えたり、画像情報を取り込み、モニタに表示したりする。また、適宜の形式のフォーマットでデータを保存する。パソコン２８は、上に示したオフセット処理、直線化処理及び光むら補正処理も実現するものである。
第２図は光学系１４の具体例を示したものである。（ａ）は光学系の外観図であり、（ｂ）はその垂直断面図、（ｃ）は（ｂ）における円内の拡大図である。
この光学系は、レンズ６から被測定物２までの距離、及びレンズ６からエリアセンサ８までの距離を自在に微調整できる構造になっており、焦点あわせ、倍率の変更などが容易にできる。また、被測定物２は試料台のベースプレート３ごと交換できるようになっている。
エリアセンサ８として三菱製ＣＭＯＳイメージセンサ（Ｈ６４２８３ＦＰ）を使用して反射率測定を行なった結果について説明する。
まず、エリアセンサ８の補正処理について説明する。
（１）オフセット処理（ダーク処理）：
ＬＥＤ４の電流値を０（ｍＶ）としたときの、エリアセンサ８の出力（Ａ／Ｄカウント値）をダーク（オフセット）とした。以降に記述する全ての演算結果（補正処理、反射率演算等）は、ＬＥＤ４照射時のエリアセンサ８の出力（Ａ／Ｄカウント値）とダーク成分の差をエリアセンサ８の生出力（Ａ／Ｄカウント値）とした。
（２）光量とエリアセンサ出力の関係（直線化補正）：
被測定物２に対してＬＥＤ４が放つ光の量とエリアセンサ８の出力（ＶｏｕｔをＡ／Ｄ変換したカウント値）とは正比例の関係にはない。
第３図のグラフは、白板（ＮＤ（ニュートラル・デンシティ）９．５、反射率実測値８７．００％（反射率実測値は分光測色計（ＭＩＮＯＬＩＡ  ＣＭ−５０３ｃ）により測定した（以下も同じ））を被測定物２とし、ＬＥＤ４の光量を変化させたときのエリアセンサ８の出力を示している。横軸は光学系１４内に配置した光検出器１０の出力（ｍＶ）、縦軸はエリアセンサ８が最も強く光を受ける部分の連続して並ぶ適当な５ピクセルのデータを平均した値である。
光むら補正を行なうためには、エリアセンサ８の出力（Ａ／Ｄカウント値）を光検出器１０の出力（ｍＶ）に換算する処理を前処理として行なう必要がある。光むら補正を行なう前に、２５℃の環境下で第３図の特性を測定し、その結果をもとにエリアセンサ８の各ピクセルの出力を補正演算する。
（３）画像のむら（光むら）補正：
この反射率測定装置を用いて撮像した生の画像情報は、エリアセンサ８の各ピクセル感度の個体差、ＬＥＤ４の照射むら、レンズ６のコサイン４乗則の影響などを受けて光むらを生じている。第４図は、白板（ＮＤ９．５、反射率実測値８７．００％）を撮像し（画角範囲内は全て白板領域）、その画情報を３次元等高面グラフにしたものである。等高面は、画像を１０×１０の領域に分割し、各領域に含まれるピクセルの平均値から生成した。
白板のように濃淡具合の均一な平面を撮像した場合にも、光むらの影響を受けて画角内の濃淡情報がドーム状に変形してしまっていることが第４図のグラフから読み取れる。このドーム状に変形した画像情報を濃淡の一様な水平面に補正する処理は、エリアセンサ８を反射率測定装置の光検出器として用いる場合に必須である。本発明では、測定結果は、すべてこの光むら補正を施したものである。
本発明では、以下の手順によって補正処理を行なった。第５図を参照してその手順を説明する。
＜補正参照データ獲得手順＞
（１）  白板（ＮＤ９．５、反射率実測値８７．００％）を被測定物２とし、撮像した画像内の最も明るいピクセルが飽和光量に達するときの光検出器（ＰＤ）１０の電圧値を求める（ステップＳ１〜Ｓ５）。
（２）  ０（ｍＶ）からピクセルが飽和光量に達したときの光検出器１０の電圧値を２０等分し、２１段階の各電圧値を求め、低い方からＰ０〜Ｐ２１とする（ステップＳ６）。
（３）  ＬＥＤ４の光量を調整し、光検出器１０の電圧値が各段階になるようにセットする。それそれの光量で白板を撮像しデータを記憶する（２１枚の画像データが得られる。０（ｍＶ）のときの画像はダークデータである）（ステップＳ７〜Ｓ１１）。
（４）  全ての画像のデータをオフセット処理する（ピクセルごとに各画像データからダークデータの値を引く）（ステップＳ１２）。
（５）  画像内で最も明るいピクセル付近の連続して並ぶ５ピクセルの値を平均する。この処理を各画像に対して行ない、光検出器１０の電圧値に対するエリアセンサ出力の関係（第３図）を獲得する（ステップＳ１３）。
（６）  ２１枚の画像データのなかで、飽和光量×０．８あたりの画像データを光むら補正参照用白板データとする（ステップＳ１４）。
＜測定画像の光むら補正手順＞
（７）  測定画像の１２８×１２８ピクセル分のＡ／Ｄデータを第３図の関係から光検出器１０の電圧値に変換（ＰＤ電圧値変換：直線化処理）する（ステップＳ１５）。変換は、第３図のグラフのサンプル点間を直線補間することにより行なう。
（８）  （６）で獲得した光むら補正参照用白板データに対しても同様に、ＰＤ電圧値変換を行なう。
（９）  光むら補正参照用白板データ（ＰＤ電圧値変換後）に対する測定画像のデータ（ＰＤ電圧値変換後）の比を１２８×１２８の各ピクセルに対して求める。この比を光むら補正データとする（ステップＳ１６）。
（ピクセル補正例）
前述までの補正方法により、生の各ピクセルの出力から反射光量を導いた例を示す。補正するピクセルは第６図に示す５点であり、それぞれ点１（３２，３２）、点２（９６，３２）、点３（６４，６４）、点４（３２，９６）、点５（９６，９６）とした。
第７図（ａ）のグラフは、白板（ＮＤ９．５、反射率実測値８７．００％）を被測定物２として、ＬＥＤ４の光量を変化させたときの、第６図で示す５点のエリアセンサの出力である。
これを第３図のグラフの関係からエリアセンサ出力（Ａ／Ｄカウント値）を光検出器の電圧値に変換（ＰＤ電圧値変換）すると、第７図（ｂ）のように補正される。第７図（ｂ）では各点の反射光量に、光むらの影響なとを受けて差異が生じているが、各点はＬＥＤ４の光量に対して正比例の関係を持つようになる。
第７図（ｃ）は（ｂ）に対して光むら補正データを用いて光むらを補正したグラフである。各点はほぼ同一線上にプロットされている。（ｃ）でエリアセンサの補正出力が１のときに完全に５点が一致するのは、この明るさの白板データから光むらを補正しているためである。（ｃ）からは光量が下がるほど各点がばらつき、補正精度が悪くなっていることも読み取ることができる。
（領域補正例）
第８図（ａ）のグラフは、ＬＥＤ光量の暗いものから明るいものまで３段階に分けて白板（ＮＤ９．５、反射率実測値８７．００％）を撮像し（画角範囲内は全て白板領域）、その画情報を並べて３次元等高面グラフにしたものである。等高面は、画像を１０×１０の領域に分割し、各領域に含まれるピクセルの出力の平均値から生成している。３つのドーム状の白板データのうち、左端のものが最も光量が少なく、右端のものが最も光量が多い。
右端の白板データは、真中の白板データに比べて光量が多いにもかかわらず最大値と最小値の差が狭くなっている。これは白板の明るい部分のピクセル光量が飽和量に近づいているためである。
第８図（ａ）のグラフに対して光むら補正を行なうと、第８図（ｂ）のグラフのように平坦になる。
（同時再現性１）
同一ピクセルを対象に異なる反射率のＮＤペーパを撮像し、それらの濃淡値の比を反射率として演算した場合の同時再現性を検証した。
手順は次の通りである。
（１）ＬＥＤ４の電流値を０ｍＡにし、ダーク（オフセット）画像を撮像する。
（２）ＮＤ９．５（反射率実測値８７．００％）とＮＤ６．５（反射率実測値３６．２１％）をそれそれベースプレートに貼り付けたものを用意し、これらを交互に１０回撮像する。
（３）各ピクセルを光むら補正した後で、１枚の画像を第９図のように１１×１１の領域（１つの領域は１０×１０＝１００ピクセル）に分割し、各領域の光量の平均値を算出する。この光量の平均値のＮＤ９．５とＮＤ６．５の比を反射率として、各領域について演算する。
表１は、１０回測定したうちの１回分の測定結果である。表の上段はＮＤ９．５を撮像したときの各領域の光量の平均、中段はＮＤ６．５を撮像したときの各領域の光量の平均、下段は各々同じ領域の比を反射率として求めたものである。
0003754989.tif 000002
表中で、ＡＶＥ．は平均値、Ｃ．Ｖ．（％）は変化率で、（標準偏差／平均値）を表わしている。Δは領域内の最大値と最小値との差を表わす。
表２は、１０回測定した結果の各領域の反射率の平均値（上段）と各領域の反射率のばらつき（下段）を示したものである。
0003754989.tif 000003
表１、表２から各領域の反射率を比較すると、レンズ８の光軸付近（あるいはＬＥＤ４の照射光が最も集光している部分）のばらつきが最も少なく、そこから同心円状に離れていくほどばらつきが多くなる傾向がみられる。これは、光軸から遠いほど補正量が多くなるためであると考えられる。
また、エリアセンサを用いて反射率を測定する場合の特徴として、濃淡が一様なはずの被測定物を測定しても、各領域別に得られる反射率にかなりの差異があることがわかる。原因としては、光むら補正精度が場所によって異なることや、もともと被測定物にある濃淡むらの影響が考えられる。
（同時再現性２）
同一画像内に異なる反射率のＮＤペーパを配置し、それらの濃淡値の比を反射率として演算した場合の同時再現性を検証した。手順は次の通りである。
（１）ＬＥＤ４の電流値を０ｍＶにし、ダーク（オフセット）画像を撮像する。
（２）ＮＤ９．５（反射率実測値８７．００％）とＮＤ６．５（反射率実測値３６．２１％）が画角に半分づつ入るようにベースプレートに貼り付けたものを用意し、これを１０回撮像する。
（３）各ピクセルを光むら補正した後で、１枚の画像を第１０図のように１１×１１の領域（１つの領域は１０×１０＝１００ピクセル）に分割し、各領域の光量の平均値を算出する。この光量の平均値のＮＤ９．５とＮＤ６．５の比を反射率として、各領域の反射率を演算する。
表３は、１０回測定したうちの１回分の測定結果である。表３の上段左側ＮＤ９．５、上段右側はＮＤ６．５の各領域の光量の平均である。下段の左側は、画像内でＮＤ９．５とＮＤ６．５の交わる部分を中心線とし、これから線対称に求めた比を反射率としたものである（対称反射率と記述）。また、下段の右側は、その中心線で領域をＮＤ９．５とＮＤ６．５の領域に分割し、それぞれの同じ領域（例：横軸１０の領域と５０の領域、横軸５０の領域と１１０の領域）の比を反射率としたものである（一方向反射率と記述）。
0003754989.tif 000004
表４は、表３下段のような演算を１０測定分行なった場合の平均値（上段）とばらつき（下段）を示している。
0003754989.tif 000005
同時再現性１の結果に比べて、同時再現性２の結果はＣ．Ｖ．（％）（変化率：標準偏差／平均値）で約２倍程度良好である。これは同時再現性１の測定では、測定毎に被写体を手で交換しなければならないのに対して、同時再現性２の測定では全く手を触れる必要がなかったためであると考えられる。つまり、同時再現性２の結果は純粋なＣＭＯＳエリアセンサの撮像再現性に近いものであると考えられる。
（反射率直線性）
精度管理された既存の分光色計測（ＭＩＮＯＬＴＡ  ＣＭ−５０３ｃ：エリアセンサを用いたものではない）にて、異なる反射率のＮＤペーパを複数種類測定し、本発明の実施例の反射率測定装置との相関を検証した。
あらかじめ、分光色測計により対象とする複数種のＮＤペーパの反射率を測定しておく。反射率はＮＤペーパ上の無作為な位置５点について測定し、平均したものを用いる。
（１）ＬＥＤ４の電流値を０ｍＡにし、ダーク（オフセット）画像を撮像する。
（２）あらかじめ分光色測計にて計測しておいたＮＤペーパを撮像する。
（３）第６図に示されるように、各画像の一様に分布した５点（ピクセル）について光むら補正を行なう。
第１１図のグラフは、横軸が分光測計で測定した値、縦軸が本発明により第６図で示す各画像の５点（ピクセル）を光むら補正してプロットしたものである。
第１１図のグラフの結果は領域平均ではなくピクセル単位での結果を示しているが、光軸に近いピクセル点３（９６，３２）などの直線性は良好であるといえる。点１（３２，３２）は５つのピクセルの中では最も暗い値を示す（生データにて）ピクセルであるが、直線性は５つのピクセルの中で最も劣っている。この実験からも、光軸から離れた部分の光むら補正が難しいことが読み取れる。
（温度特性）
実施例の反射率測定装置の温度特性を把握するための測定を行なった。
以下の操作を１０℃、２０℃、３０℃の各環境でシステム（電源ＯＮ状態）を十分になじませた後行なう。ＮＤ９．５（反射率実測値８７．００％）とＮＤ６．５（反射率実測値３６．２１％）が画角に半分づつ入るようにベースプレートに貼り付けたものを被測定物とする。
（１）ＬＥＤ４の電流値を０ｍＡにし、ダーク（オフセット）画像を撮像する。
（２）各環境温度でＬＥＤ４の電流値を（１０℃：１６．５２（ｍＶ）、２０℃：１７．２０（ｍＡ）、３０℃：１７．９５（ｍＡ））にセットし、光検出器１０が検知するＬＥＤ光量が（１０℃：０．７８８（Ｖ）、２０℃：０．７８６（Ｖ）、３０℃：０．７８３（Ｖ））を超えるのを待つ。
（３）（２）の条件を満たした直後、撮像する。以上の作業を１０回繰り返す。
表５は各１０回の測定による各温度での全領域平均反射率の結果を示したものである。Ｓ．Ｄ．は標準偏差を表わす。
0003754989.tif 000006
この結果から、環境温による影響はほとんどなく、およそ０．２８（％／１０℃）の程度の温度傾向であった。
（ドリフト特性）
使用状態（時間、温度を含む）における、実施例の反射率測定装置のドリフト傾向を把握するための測定を行なった。
（１）実施例の反射率測定装置の主なユニット（ドライブ回路２０内、ＬＥＤ４付近、エリアセンサ８付近）に熱電対を装着して温度をモニタできるようにする。
（２）反射率測定装置を電源ＯＦＦの状態で十分に環境になじませる。
（３）ＬＥＤ４の電流値を０ｍＶにし、ダーク（オフセット）画像を撮像する。
（４）ＬＥＤ４の電流値を１７．３（ｍＶ）にセットし、光検出器１０が検知するＬＥＤ光量が０．７８９（Ｖ）を超えるのを待つ。
（５）（２）の条件を満たした直後、３回撮像する。
（６）（３）〜（５）の処理を１０分おきに繰り返し、モニタしているユニット温度が全て平衡になるまで行なう。
第１２図のグラフは、実施例の反射率測定装置の主なユニット（ドライブ回路２０付近、ＬＥＤ４付近、エリアセンサ８付近）の時間経過（１０分おき）と温度の関係を示している。
第１３図のグラフは１０分おきに３回ずつ測定した反射率の結果をプロットしたものである。
第１２図、第１３図の結果から、使用環境（温度、時間を含む）におけるドリフト現象は確認できず、仮りに存在するとしても、同時測定毎に発生するばらつきの中に包含される程度である。
以上の検証の結果、実施例の反射率測定装置は、同時再現性（ｎ＝１０）では、Ｃ．Ｖ．＝０．２３％（反射率４５％付近にて）、温度特性としては反射率４５％付近において、およそ０．２８（％／１０℃）、使用状態（時間、温度を含む）においてドリフトの傾向はほとんど見られないことが確認された。
実施例で使用したＣＭＯＳエリアセンサは、尿試験紙の測定機など半定量レベルの測定には十分適用できることが判明した。
［実施例２］
第２の実施例として、センサとしてエリアセンサを使用し、本発明の第２の局面による出力補正方法を適用した二次元の反射率測定装置の一例の概要を第１４図に示す。
第１図の反射率測定装置と比較すると、光量をモニタする光検出器１０が配置されていない点で異なる。他の構成は基本的に同じである。被測定物２の反射光が反射板５を介してレンズ６によりエリアセンサ８ａに結像される。エリアセンサ８ａは第１図における増幅器２２まで含んだものを示している。エリアセンサ８ａの検出信号はＡ／Ｄ変換器２４を経て演算部２８ａに取り込まれる。演算部２８ａは、第１図におけるＲＡＭ２６とパソコン２８に該当するものである。演算部２８ａには、表示器３０、キーボード３２、プリンタ３４が接続されている。３６は取り込んだ画像データを保存するイメージ保存部であり、例えばハードディスク装置により構成される。演算部２８ａで算出された反射率を濃度に変換するために検量線データ４０がハードディスク装置やフロッピーディスク装置に保存されている。
演算部２８ａのデータ処理結果は外部出力３８として必要な外部機器に取り出すことができる。
この実施例ではエリアセンサ８ａの出力と被測定物２の反射率との関係を直線化データとして取得するために、被測定物２として反射率が既知の標準板を測定する。標準板としてはＮＤペーパーを使用し、反射率が最も大きいものから最も小さいものまで１１段階のものを用意する。
それらの標準板を被測定物２として測定した結果をエリアセンサ８ａの出力とともに示したものが第１５図である。縦軸は出力を表わし、横軸は各標準板を反射率の大きいものから順に並べたものである。各標準板の出力データは光むら補正処理を施していないので湾曲したものとなっている。
エリアセンサ８ａの１つのピクセルについて反射率と出力の関係を示したものが第１６図である。縦軸エリアセンサ８ａの出力を表わし、横軸は標準板の既知の反射率を表している。エリアセンサ８ａは受光量に対して出力が非直線性を持っているので、この曲線はＳ字型を示しており、第３図に示したものと同じことを意味している。
エリアセンサ８ａの各ピクセルについて第１６図に示されるようなデータをピクセルごとの直線化データとして保持しておく。
反射率が未知の試料を測定したとき、各ピクセルについてそれぞれの直線化データを用いてその出力から第１６図中に矢印で示されるように反射率を求める。反射率はこの直線化データの実測点の間を補間することにより得られる。
このようにして得られた未知試料の反射率は、光源による照射むらやレンズ、エリアセンサ８ａの非直線性を含めて補正された反射率データとなり、直線性を持った反射率的データとなる。
この動作を第１７図と第１８図により改めて説明する。
第１７図は直線化データを得るための手順を示したものである。基準板として反射率の異なるＮ種類のものを用意する。反射率は１００％から０まで１０％単位で変化している１１種類である。１つの基準板を被測定物２の位置に置き、エリアセンサ８ａにより撮像する。そのときのその基準板の既知の反射率ｒと撮像データを記憶する。この動作を全ての基準板について繰り返す。
これにより撮像データの各ピクセルの出力と反射率の関係を示す第１６図の直線化データがピクセルごとに得られる。
第１８図の操作では、反射率が未知の試料を被測定物の位置に置き、エリアセンサ８ａで撮像する。その撮像結果から、ピクセルの位置を示す座標（ｘ，ｙ）について、各ピクセルの出力データから第１６図に矢印で示されるように反射率を求める。この動作を全ピクセルについて行う。
［実施例３］
第３の実施例として、本発明の第３の局面による出力補正方法をセンサとしてエリアセンサを使用して説明する。
用いる反射率測定装置は第１４図に示されたものと同じである。
この実施例ではエリアセンサ８ａとして受光する露光時間がプログラム可能なものを使用する。そのようなエリアセンサとしては、例えば、第１図の実施例で使用した三菱製ＣＭＯＳイメージセンサ（Ｈ６４２８３ＦＰ）を使用する。しかし、用いるエリアセンサ８ａは、ＣＭＯＳイメージセンサに限らず、ＣＣＤイメージセンサでも露光時間がプログラム可能なものであれば、使用することができる。
エリアセンサ８ａの出力は受光量に対して直線性を持っていないが、受光量は露光時間に比例する。また、受光量は反射率に比例するので、１つの基準板を使用して露光時間を異ならせることにより、基準板を共通に使用しても露光時間を変えることによって反射率の異なる基準板を用いた測定と等価な結果を得ることができる。
第１４図における被測定物２の位置に基準板となる白板を載置する。まず基準の露光時間で測定する。次に被測定物の白板はそのままにして、露光時間を基準の露光時間の９０％に下げて同じ測定を行う。同様にして８０％、７０％…というように露光時間を変化させていく。
第１９図はそのように露光時間を減少させていったときのエリアセンサ８ａの出力（縦軸）と露光時間（横軸。右側ほど露光時間が短い）を示したものである。この場合もエリアセンサ８ａ内での光むら補正を行っていないので、ピクセル間での出力は変動している。
各ピクセルについて出力と露光時間の関係を図示すると、第２０図のようになり、この結果第１６図及び第３図と同じものである。各ピクセルについて第２０図のデータを直線化データとして記憶しておく。
第２１図はこの直線化データを取得する手順をまとめて示したものである。基準白板を被測定物２としてセットし、基準露光時間ｔをセットする。その露光時間で照射し、エリアセンサ８ａによる撮像を行ない、露光時間ｔと撮像データを記憶する。
次に、露光時間を１０％減少して同じ測定を繰り返す。このように、露光時間が順次減少するように設定して繰り返して測定していき、各ピクセルについてセンサ出力と露光時間の関係を図示したものが第２０図である。横軸の露光時間は反射率に該当している。
第２２図は反射率が未知の試料を測定したときの手順を示したものであり、第１８図に示した手順と同じである。この実施例ではピクセルごとに反射率に該当する露光時間が得られる。
この実施例において補正後のデータの精度を確認した。反射率が既知の複数の被写体を測定した。エリアセンサ８ａの中央部と周辺部の各位置でのピクセルを第２３図のように選び出し、それぞれのピクセルにおける補正された出力と反射率の関係を第２４図に示す。各位置でのピクセルの出力は、１個のピクセルの出力であってもよく、その位置周辺での幾つかのピクセルの出力の平均値であってもよい。直線は全ピクセルの平均値である。第２４図で、横軸は既知の反射率、縦軸は補正された出力である。
第２４図の結果から、エリアセンサ８ａ内のピクセルの位置にかかわらず、照射むらやレンズ、エリアセンサの非直線性を含めて直線的な反射率的値に補正されていることがわかる。
産業上の利用可能性
本発明の補正方法は、試薬部を支持体に設けた試験片を用いるドライケミストリ分析計を初め、種々の分析計におけるセンサ出力の補正方法として、臨床検査、食品分析、化学分析などの多くの分野で利用することができる。
【図１】
0003754989.tif 000007
【図２】
0003754989.tif 000008
【図３】
0003754989.tif 000009
【図４】
0003754989.tif 000010
【図５】
0003754989.tif 000011
【図６】
0003754989.tif 000012
【図７】
0003754989.tif 000013
【図８】
0003754989.tif 000014
【図９】
0003754989.tif 000015
【図１０】
0003754989.tif 000016
【図１１】
0003754989.tif 000017
【図１２】
0003754989.tif 000018
【図１３】
0003754989.tif 000019
【図１４】
0003754989.tif 000020
【図１５】
0003754989.tif 000021
【図１６】
0003754989.tif 000022
【図１７】
0003754989.tif 000023
【図１８】
0003754989.tif 000024
【図１９】
0003754989.tif 000025
【図２０】
0003754989.tif 000026
【図２１】
0003754989.tif 000027
【図２２】
0003754989.tif 000028
【図２３】
0003754989.tif 000029
【図２４】
0003754989.tif 000030