(57)【特許請求の範囲】
【請求項１】
観測対象の状態量を検知媒体の呈色により計測する状態量計測方法であって、
少なくとも２つの既知の状態量に対応する呈色部位を含めて観測対象である呈色部位を同時に撮影した画像データを受け付けるステップと、
前記画像データから前記呈色部位のそれぞれに該当する領域を検出するステップと、
前記領域のそれぞれにおいて表色系に基づく指数を算出するステップと、
算出した前記少なくとも２つの既知の状態量に対応する呈色部位の前記指数を用いて、前記状態量と前記観測対象である呈色部位の前記指数との関係を表す関数を求め、該関数に算出した前記観測対象である呈色部位の前記指数を代入し、前記観測対象の状態量を判定するステップと、
を有することを特徴とする状態量計測方法。
【請求項２】
前記指数は、前記領域内の画像データを統計処理することにより算出することを特徴とする請求項１記載の状態量計測方法。
【請求項３】
前記状態量と前記指数との関係は線形であることを特徴とする請求項１又は２記載の状態量計測方法。
【請求項４】
前記表色系は、ＲＧＢ表色系、ＨＳＬ表色系、グレースケールのいずれかであることを特徴とする請求項１乃至３のいずれか１つに記載の状態量計測方法。
【請求項５】
前記観測対象は、検知紙あるいは検知管のいずれかであることを特徴とする請求項１乃至請求項４のいずれか１つに記載の状態量計測方法。
【請求項６】
観測対象の状態量を検知媒体の呈色により計測する状態量計測装置であって、
少なくとも２つの既知の状態量に対応する呈色部位を含めて観測対象である呈色部位を同時に撮影した画像データを受け付ける入力手段と、
前記画像データから前記呈色部位のそれぞれに該当する領域を検出する領域検出手段と、
前記領域のそれぞれにおいて表色系に基づく指数を算出する指数算出手段と、
算出した前記少なくとも２つの既知の状態量に対応する呈色部位の前記指数を用いて、前記状態量と前記観測対象である呈色部位の前記指数との関係を表す関数を求め、該関数に算出した前記観測対象である呈色部位の前記指数を代入し、前記観測対象の状態量を判定する状態量判定手段と、
を有することを特徴とする状態量計測装置。
【請求項７】
前記指数は、前記領域内の画像データを統計処理することにより算出することを特徴とする請求項６記載の状態量計測装置。
【請求項８】
前記状態量と前記指数との関係は線形であることを特徴とする請求項６又は７記載の状態量計測装置。
【請求項９】
前記表色系は、ＲＧＢ表色系、ＨＳＬ表色系、グレースケールのいずれかであることを特徴とする請求項６乃至８のいずれか１つに記載の状態量計測装置。
【請求項１０】
前記観測対象は、検知紙あるいは検知管のいずれかであることを特徴とする請求項６乃至請求項９のいずれか１つに記載の状態量計測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、比色法により観測対象の状態量を計測する技術に関する。
【背景技術】
【０００２】
ガスセンサの分野では、ガス濃度による色の変化を利用して対象のガスを検知する比色法が広く研究されている。例えば、ガス検知管による方法は、色の変化領域と未変化領域の境界における濃度目盛を目視により読み取ることでガス濃度を測定する。しかしながら、目視による読み取り誤差がかなりの濃度範囲で生じてしまうという問題がある。
【０００３】
また、検知紙の色変化をカラーチャートと比較することによりガス濃度を測定する方法も知られている。検知紙を用いる場合も、カラーチャートと比較する際の個人による誤差が大きいという問題がある。
【０００４】
発明者らも、オゾンや二酸化炭素、ホルムアルデヒドなどと反応することで色が変化する高感度の検知素子を製作してきた。これら検知素子は高感度であるため、目視による読み取り誤差をかなり改善することができた。しかし、依然として読み取りの個人差による誤差は生じてしまう。
【０００５】
分光光度計などの光学測定器を用いて光の反射率を測定し、色の変化を数値化することで読み取り誤差を防ぐことができるが、これらの装置は、高額であることや利用方法が限られるなどの制約がある。
【０００６】
カメラなどにより撮影した画像から色を識別する方法が知られているが、撮影時の証明の違いや撮影機材の違いによる撮影画像の色のバラツキが生じるという問題がある（非特許文献１参照）。例えば、特許文献１に記載の技術のように、ビンやペットボトルを色情報から識別する場合は、数種類の色のちがいが判別できればよく、車の色の判別でも青、赤、白など基本色の違いを判定する程度でよい。
【特許文献１】特開平１０−６２２５１号公報
【非特許文献１】"色色雑学"、[online]、コニカミノルタ、［平成１９年１１月９日検索］、インターネット〈ＵＲＬ：http://konicaminolta.jp/instruments/colorknowledge/part2/07.html〉
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、ガスセンサとしての水溶液や検知紙の色の変化は、従来の色識別で想定している色の違いよりも高分解能で比色識別ができなくてはいけない。
【０００８】
さらに、検知紙などの場合は、ある色から別の色に徐々に変化するものであり、その変化の度合いが識別できなければならない。例えば、最初は青色だった検知紙がガスと反応することで少しずつ白くなる場合、その青色が薄くなる程度の違いが識別できなければならない。
【０００９】
本発明は、上記に鑑みてなされたものであり、その課題とするところは、電子画像を用いた比色法において、撮影時の環境のバラツキによる影響を抑制し、より安定して状態量を計測することにある。
【課題を解決するための手段】
【００１０】
第１の本発明に係る状態量計測方法は、観測対象の状態量を検知媒体の呈色により計測する状態量計測方法であって、少なくとも２つの既知の状態量に対応する呈色部位を含めて観測対象である呈色部位を同時に撮影した画像データを受け付けるステップと、画像データから呈色部位のそれぞれに該当する領域を検出するステップと、領域のそれぞれにおいて表色系に基づく指数を算出するステップと、算出した少なくとも２つの既知の状態量に対応する呈色部位の指数を用いて、状態量と観測対象である呈色部位の指数との関係を表す関数を求め、該関数に算出した観測対象である呈色部位の指数を代入し、観測対象の状態量を判定するステップと、を有することを特徴とする。
【００１１】
本発明にあっては、２つの既知の状態量に対応する呈色部位と観測対象である呈色部位とを含む画像データを入力し、呈色部位それぞれにおける表色系指数を算出して、状態量と表色系指数との関係を表す関数に算出した表色系指数を代入することにより、撮影環境などの影響を抑制してより正確に状態量を計測することが可能となる。既知の状態量に対応する呈色部位も観測対象の呈色部位と同様に、撮影環境などの影響により画像データの色がばらつくが、このばらつきは、基準となる呈色部位も観測対象の呈色部位も同じ法則に従うので、これらの呈色部位を同時に撮影した画像データを入力することで、基準となる呈色部位からの変化量に基づいて状態量を計測することができる。
【００１２】
また、安価でさまざまな場面で利用可能な携帯電話のカメラ、デジタルカメラあるいはインターネットカメラなどの撮影機材を使うことができるので、高価で利用方法が限定される分光光度計などの光学測定器が不要となる。
【００１３】
上記状態量計測方法において、指数は、領域内の画像データを統計処理することにより算出することを特徴とする。
【００１４】
上記状態量計測方法において、状態量と指数との関係は線形であることを特徴とする。
【００１５】
上記状態量計測方法において、表色系は、ＲＧＢ表色系、ＨＳＬ表色系、グレースケールのいずれかであることを特徴とする。
【００１６】
上記状態量計測方法において、観測対象は、検知紙あるいは検知管のいずれかであることを特徴とする。
【００１７】
第２の本発明に係る状態量計測装置は、観測対象の状態量を検知媒体の呈色により計測する状態量計測装置であって、少なくとも２つの既知の状態量に対応する呈色部位を含めて観測対象である呈色部位を同時に撮影した画像データを受け付ける入力手段と、画像データから呈色部位のそれぞれに該当する領域を検出する領域検出手段と、領域のそれぞれにおいて表色系に基づく指数を算出する指数算出手段と、算出した少なくとも２つの既知の状態量に対応する呈色部位の指数を用いて、状態量と観測対象である呈色部位の指数との関係を表す関数を求め、該関数に算出した観測対象である呈色部位の指数を代入し、観測対象の状態量を判定する状態量判定手段と、を有することを特徴とする。
【発明の効果】
【００１８】
本発明によれば、電子画像を用いた比色法において、撮影時の環境のバラツキによる影響を抑制し、より安定して状態量を計測することができる。
【発明を実施するための最良の形態】
【００１９】
以下、本発明の実施の形態について図面を用いて説明する。
【００２０】
図１は、本実施の形態における状態量計測システムの構成を示すブロック図である。同図に示す状態量計測システムは、状態量計測装置１０と、撮影装置３０とを有する。状態量計測装置１０と撮影装置３０とはネットワーク２０を介して接続される。撮影装置３０は、状態量が既知の呈色部位４１，４３と状態量が未知の呈色部位４２とを含む台紙４０を撮影し、得られた画像データをネットワーク２０を介して状態量計測装置１０に送信する。画像データを受信した状態量計測装置１０は、呈色部位４１，４２，４３における表色系指数をそれぞれ算出し、呈色部位４１，４３における表色系指数に基づいて呈色部位４２の状態量を判定する。
【００２１】
図２は、状態量計測装置１０の構成を示すブロック図である。同図に示す状態量計測装置１０は、入力部１１、領域検出部１２、指数算出部１３および状態量判定部１４を備える。
【００２２】
入力部１１は、少なくとも２つの状態量が既知の呈色部位４１，４３と観測対象の検知媒体の呈色部位４２を含む画像データの入力を受け付ける。状態量が既知とは、状態量が数値化されていることであり、検知媒体が呈色部位４１，４３と同じ色になるための状態量が分かっている。基準色となる呈色部位４１，４３には、色が変化する前のものと、完全に色が変化したものを用いることが望ましい。
【００２３】
領域検出部１２は、画像データにおける各呈色部位４１，４２，４３の位置、領域を検出する。呈色部位４１，４２，４３の位置や領域を検出するためには、予めそれぞれを示すマーカを設ける方法や、特定の色とそれ以外の色とを分離する処理（色抽出）などを用いる。色抽出は、カラー画像処理に用いられる前処理手法のひとつして一般的に行われており、色抽出処理ソフトウェアモジュールも入手可能である。
【００２４】
指数算出部１３は、各呈色部位４１，４２，４３の領域における所定の表色系に基づく指数を求める。表色系としては、ＲＧＢ表色系、ＨＳＬ表色系、マンセル表色系、オストワルト表色系およびＰＣＣＳ（Practical Color Coordinate System：日本色研配色体型）などが知られている。
【００２５】
ＲＧＢ表色系は、色を赤（Ｒｅｄ）、緑（Ｇｒｅｅｎ）、青（Ｂｌｕｅ）の３色の組み合わせにより表現したものである。コンピュータのディスプレイで用いられる。
【００２６】
ＨＳＬ表色系は、色を色相（Ｈｕｅ）、彩度（Ｓａｔｕｒａｔｉｏｎ）、輝度（Ｌｕｍｉｎａｎｃｅ）で表現したものである。
【００２７】
マンセル表色系は、色を色相（ｈｕｅ）、彩度（ｃｈｒｏｍａ）、明度（ｖａｌｕｅ）の３つの属性で表現する。色相は、赤、黄、緑、青、紫の５色を基本とする。明度は、黒（０）から白（１０）までの１１段階で表現する。各色相、明度ごとに無彩色（灰色）から最も鮮やかな色（単色）までに分ける。ＪＩＳに採用されている。
【００２８】
オストワルト表色系は、全ての色を、白、黒、純色（各色相）の混合として表す。白、黒、純色を３つの頂点とする正三角形すなわち色三角形を無彩色（白、灰、黒）を軸として構成する。
【００２９】
表色系に基づく指数とは、例えば、ＲＧＢ各８ビットのＲＧＢ表色系で表現された画像データを用いる場合、領域内の画素それぞれのＲＧＢの画像信号値（０−２５５）の平均や最頻値など、数値化した色を統計処理することにより算出したものである。
【００３０】
指数算出部１３で利用する表色系としては、表色系に基づく指数と計測対象の状態量とが関数関係で表せるものならば任意の表色系を利用することができる。グレースケールを用いてもよい。
【００３１】
状態量判定部１４は、状態量と指数との関係を表す関数に指数算出部１３が算出した指数を代入し、計測対象の状態量を判定する。基準色となる呈色部位４１，４３における表色系に基づく指数と、計測対象の呈色部位４２における表色系に基づく指数との距離関係を求めることで呈色部位４２の状態量を判定することができる。ここで、距離とは、状態量が既知の呈色部位４１と呈色部位４３の差分量に対する呈色部位４１と状態量を求める呈色部位４２の差分量の割合のことである。呈色部位４２の状態量は、基準となる呈色部位４１からの変化量から計算することができる。表色系に基づく指数は撮影環境など影響されるが、複数の基準となる呈色部位４１，４３と計測対象の呈色部位４２とを同時に撮影することで、求める状態量におけるその影響を抑制することができる。
【００３２】
［実施例１］
次に、上記状態量計測システムにより、状態変化と呈色が線形関係にある青色のシート状の対象物を用いてガス濃度を判定する場合について説明する。観測する対象物は、ガス濃度に応じて青色が脱色するものである。例えば、市販のオゾン検知シートは、最初は濃い青色であるが、空気中に置いておくことによって、空気中に含まれるオゾンと反応して徐々に青色が薄くなる。十分に反応が進むとオゾン検知シートは脱色されて基板色（白色）に近くなる。
【００３３】
図３に、状態量計測システムの動作の流れを示す。まず、撮影装置３０は、複数の呈色部位４１，４２，４３を含む台紙４０を撮影し、呈色部位４１，４２，４３を含む画像データをネットワーク２０を介して状態量計測装置１０へ送信する。入力部１１は、画像データを受信し、メモリなどの記憶装置に記憶させる（Ｓ３０１）。ここで呈色部位４１は、未使用のオゾン検知シートであり、濃い青色である。一方、呈色部位４３は、十分に反応させて脱色したオゾン検知シートである。呈色部位４２は、状態量が未知の観測対象のオゾン検知シートである。
【００３４】
なお、画像データは、必ずしもネットワーク２０を介して状態量計測装置１０へ入力する必要はなく、磁気ディスク、光ディスクあるいは半導体メモリーなどの記憶媒体を介して状態量計測装置１０へ入力するものであってもよいし、状態量計測装置１０が撮影装置３０を内蔵する構成としてもよい。
【００３５】
続いて、領域検出部１２は、記憶装置から画像データを読み出し、各呈色部位４１，４２，４３の領域を検出する（Ｓ３０２）。
【００３６】
続いて、指数算出部１３は、各領域の画像データのＲＧＢ値のうちＲの値の平均を計算して、各呈色部位４１，４２，４３の指数とする（Ｓ３０３）。図４は、ガス濃度（オゾンの暴露量）とオゾン検知シートにおける画像信号ＲＧＢそれぞれの画像信号値との関係を示した図である。同図に示すように、ガス濃度に対応してもっとも変化量が大きい画像信号がＲであり、画像信号Ｒの画像信号値から精度よくガス濃度を数値化できる。したがって、画像信号値Ｒを用いて表色系指数を算出する。本実施例では、Ｒの値を用いたが、これに限定するものではなく、他の値を用いても良いし、別の表色系を用いてもよい。
【００３７】
そして、状態量判定部１４は、測定対象の状態量と表色系指数との関係を表す関数に状態量が既知の呈色部位４１，４３の表色系指数と、観測対象である呈色部位４２の表色系指数とを代入して状態量を求める（Ｓ３０４）。本実施例では、状態変化と呈色が線形関係にあるので、次式（１）で示す線形関数を用いてガス濃度を数値化する。
【数１】
0005039520.tif 000002
【００３８】
ここで、Ｒ（青）は、未使用時のオゾン検知シートから算出した表色系指数であり、Ｒ（白）は、十分に反応させて脱色したオゾン検知シートから算出した表色系指数であり、Ｒ（対）は、観測対象のオゾン検知シートから算出した表色系指数である。また、Ａは、オゾン検知シートの反応飽和状態に対応する濃度である。
【００３９】
図５は、上記の状態量計測システムを使用し、６４０ｐｐｂ×時間の積算量のオゾンで検知紙が白色となるオゾン検知紙を撮影し、オゾン濃度を計測した結果のグラフである。この実験では、外部空気を取り入れるが温湿度を一定に保った実験槽内にオゾン検知紙を２４時間吊し、そのオゾン検知紙を１時間ごとにインターネットカメラで撮影した。実験槽内のオゾン濃度は屋外の濃度変化とほぼ同じである。比較のため、図５のグラフに実験槽内のオゾン濃度を実際に測定した結果も示す。各時間の左側の斜線を付した棒グラフがオゾン濃度を測定した結果を示し、右側の棒グラフが状態量計測システムにより計測した結果である。各時間毎のオゾン濃度は、暴露した積算のオゾン濃度から１時間前の積算のオゾン濃度を引くことで、その１時間におけるオゾン濃度を算出した。図５に示すように、実験槽内で測定した濃度と、検知紙の画像から状態量計測システムが計算した濃度とはよく一致していることがわかる。
【００４０】
［実施例２］
次に、アルコールに反応して色が変化するガス検知管を用いた例を示す。このガス検知管は、呼気などに含まれるアルコールの量に応じて内部の検知剤の色が変化する。ガス検知管の色は、未使用時は黄色であるが、アルコールの量に応じて徐々に緑色に近くなる。
【００４１】
図６は、未使用時の画素信号値ＲＧＢの頻度分布を示すグラフであり、図７は、少量のアルコールに反応させた時の画素信号値ＲＧＢの頻度分布を示すグラフであり、図８は、大量のアルコールに反応させた時の画素信号値ＲＧＢの頻度分布を示すグラフである。横軸は画像信号値、縦軸は画素数を示している。それぞれの図の符号６０１，７０１，８０１は、画素信号値Ｒの分布を示し、符号６０２，７０２，８０２は、画素信号値Ｇの分布を示し、符号６０３，７０３，８０３は、画素信号値Ｂの分布を示している。各図を比較すると、アルコール量に応じて符号６０１，７０１，８０１で示す画素信号値Ｒが最も変化していることが分かる。したがって、ガス検知管の画像信号値Ｒを用いて表色系指数を算出することで、精度よくアルコール濃度を数値化することが可能である。具体的には、次式（２）で示す線形関数を用いてアルコール濃度を検出する。
【数２】
0005039520.tif 000003
【００４２】
ここで、未使用時のガス検知管の表色系指数をＲ（黄）、十分な量のアルコールと反応させた時のガス検知管の表色系指数をＲ（緑）、観測対象のガス検知管の表色系指数をＲ（対）とする。Ｂは、ガス検知管を飽和反応させるのに必要なアルコール濃度である。
【００４３】
以上説明したように、本実施の形態によれば、基準となる呈色部位４１，４３と状態量を計測する呈色部位４２とが同一の画像データに含まれるように撮影し、呈色部位４１，４２，４３それぞれの表色系指数を算出し、算出した表色系指数それぞれに基づいて呈色部位４２の状態量を求める。これにより、基準となる呈色部位４１，４３も測定対象の呈色部位４２と同様に照明や撮影機材の違いで画像データにおける色がばらつくが、そのばらつき方は呈色部位４１，４３と呈色部位４２とで変わらないので、撮影条件の影響を抑制してより正確に状態量を測定することが可能となる。また、安価でさまざまな場面で利用可能な携帯電話のカメラ、デジタルカメラあるいはインターネットカメラなどの撮影機材を使って、状態に応じて色が変化する対象物体の状態を数値化できるので、高価で利用方法が限定される分光光度計などの光学測定器を用いる必要がなくなる。
【００４４】
なお、図９に示すように、入力される画像データ９０は、ｎ個（ｎ＝１，２，３・・・）の計測対象となる呈色部位９２，９３を含むものでもよい。同図に示すように、基準となる呈色部位９１，９４を加えた（ｎ＋２）個の呈色部位９１，９２，９３，９４を含む画像データ９０を入力し、計測対象となる呈色部位９２，９３の状態量をそれぞれ測定することで、複数の呈色部位９２，９３の状態量をまとめて測定することができる。また、入力する画像データは、基準となる呈色部位を３つ以上含むものでもよい。
【００４５】
また、本願発明は、撮影時の照明スペクトル強度、照射角度、撮影機材の性能といった撮影環境によって生じる電子画像上の表色の揺らぎを抑制するために基準となる呈色部位と計測対象の呈色部位とが含まれた１枚の電子画像を用いたが、上記の撮影環境が十分管理され、電子画像上の表色の揺らぎが抑制された環境において、観測対象となる検知媒体をそれぞれ異なる電子画像とした場合であっても、電子画像毎に表色系指数を算出し、本願発明の状態量計測方法により状態量を特定することが可能である。
【図面の簡単な説明】
【００４６】
【図１】一実施の形態における状態量計測システムの構成を示す図である。
【図２】上記状態量計測システムが備えた状態量計測装置の構成を示すブロック図である。
【図３】上記状態量計測システムの動作の流れを示すフローチャートである。
【図４】オゾンの暴露量と観測対象であるオゾン検知シートの色の画像信号値との関係を示した図である。
【図５】実験槽内のオゾン濃度の測定結果と上記状態量計測システムの計測結果を示すグラフである。
【図６】アルコール用ガス検知管の未使用時の画素信号値ＲＧＢの頻度分布を示すグラフである。
【図７】アルコール用ガス検知管に少量のアルコールを反応させたときの画素信号値ＲＧＢの頻度分布を示すグラフである。
【図８】アルコール用ガス検知管に大量のアルコールを反応させたときの画素信号値ＲＧＢの頻度分布を示すグラフである。
【図９】上記状態量計測システムに入力する別の画像データを示す図である。
【符号の説明】
【００４７】
１０…状態量計測装置
１１…入力部
１２…領域検出部
１３…指数算出部
１４…状態量判定部
２０…ネットワーク
３０…撮影装置
４０…台紙
４１，４２，４３…呈色部位
【図１】
0005039520.tif 000004
【図２】
0005039520.tif 000005
【図３】
0005039520.tif 000006
【図４】
0005039520.tif 000007
【図５】
0005039520.tif 000008
【図６】
0005039520.tif 000009
【図７】
0005039520.tif 000010
【図８】
0005039520.tif 000011
【図９】
0005039520.tif 000012