(57)【特許請求の範囲】
【請求項１】
粒子が撮影された静止画像を、前記粒子の存在する領域と背景に分割する領域分割方法において、多次元色濃度座標を持つ空間上に前記静止画像の色濃度分布を表したとき、濃度分布が最も広がる方向と新たな色濃度座標軸の一軸が平行になるように前記多次元色濃度座標系を回転し、前記粒子が撮影された静止画像の色濃度を回転後の座標の座標値に変換し、該座標値に閾値を設け、前記静止画像を前記背景と前記粒子の存在する領域とに分割することを特徴とする粒子画像領域分割方法。
【請求項２】
粒子が撮影された静止画像の、赤成分の濃度（Ｒ濃度）、緑成分の濃度（Ｇ濃度）、青成分の濃度（Ｂ濃度）を用い、前記静止画像を、前記粒子の存在する領域と背景に分割する領域分割方法において、前記Ｒ濃度、Ｇ濃度、Ｂ濃度を座標軸とする色濃度座標空間上で前記静止画像の色濃度分布を表したとき、色濃度分布の中心と黒色を示す点とを結ぶ直線と、新たな色濃度座標の一軸が平行になるように前記色濃度座標系を回転し、前記Ｒ濃度、Ｇ濃度、Ｂ濃度を回転後の座標系の座標値に変換し、前記回転後の座標系の少なくとも一軸上に設けた閾値により前記静止画像を前記背景と前記粒子の存在する領域とに分割することを特徴とする粒子画像領域分割方法。
【請求項３】
粒子が撮影された静止画像の、赤成分の濃度（Ｒ濃度）、緑成分の濃度（Ｇ濃度）、青成分の濃度（Ｂ濃度）を用い、前記静止画像を、前記粒子の存在する領域と背景に分割する領域分割方法において、前記Ｒ濃度、Ｇ濃度、Ｂ濃度を座標軸とする色濃度座標空間上に前記静止画像を表したときの色濃度の最頻値を求め、前記濃度座標中で前記最頻値を示す点と、黒色を示す点とを結ぶ直線と、新たな色濃度座標の一軸が平行になるように前記色濃度座標を回転し、前記Ｒ濃度、Ｇ濃度、Ｂ濃度を回転後の座標系の座標値に変換し、該回転後の座標系の少なくとも一軸上に設けた閾値により前記静止画像を前記背景と前記粒子の存在する領域とに分割することを特徴とする粒子画像領域分割方法。
【請求項４】
請求項１〜３記載の粒子画像領域分割方法のおいて、前記静止画像を前記背景と、前記粒子の存在する領域と、前記背景と前記粒子の存在する領域との中間的な濃度値を有する領域とに分割し、前記中間的な濃度値を有する領域の周囲長と、背景に接する周囲長を求め、該背景に接する周囲長と、前記周囲長との比に閾値を設け、前記背景に接する周囲長と前記周囲長との比が前記閾値より大きい時には前記中間的な濃度値を有する領域を背景とし、前記背景に接する周囲長と前記周囲長との比が前記閾値より小さい時には、前記中間的な濃度値を有する領域を対象領域として分割することを特徴とする粒子画像領域分割方法。
【請求項５】
請求項１〜４記載の粒子画像領域分割方法において、前記濃度座標を回転する角度を示す情報を表示する工程を有することを特徴とする粒子画像領域分割方法。
【請求項６】
請求項１〜４記載の粒子画像領域分割方法において、前記濃度座標を回転させる角度に基準となる範囲を設け、前記角度が前記基準となる範囲から外れた時に、警告を発する工程を有することを特徴とする粒子画像領域分割方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は計算機を用いた画像処理に関し、特に粒子画像を背景と対象領域（粒子の存在する領域）とに分離する領域分割方法に関する。
【０００２】
【従来の技術】
従来、この種の領域分割方法としては、例えば白血球画像の領域分割方法として、「白血球分類の自動化」（Medical Imaging Technology、 Vol.14、 No.1、 January、 1996、 P.14-22)に記載されている技術が知られている。この技術は、白血球が撮影された静止画像の、赤成分の濃度（Ｒ濃度）と、緑成分の濃度（Ｇ濃度）と青成分の濃度（Ｂ濃度）に閾値を設け、背景と対象領域を分離している。この技術を応用すると、例えばＲ濃度、Ｇ濃度をそれぞれ縦軸、横軸とする座標面上において、図２に示すように背景濃度が１０３、対象領域の濃度が１０４のように分布した場合、Ｒ濃度に閾値Ｔ１、Ｔ２を設定し、Ｇ濃度に閾値Ｔ１１、Ｔ１２を設定し、Ｒ濃度がＴ１以上Ｔ２以下であり、かつＧ濃度がＴ１１以上Ｔ１２以下である濃度を有する画素を背景、それ以外の画素を対象領域とすることにより、背景と対象領域とを分離することができる。
【０００３】
【発明が解決しようとする課題】
しかしながら、尿沈渣画像に上述の従来技術を応用した場合、必ずしも全ての尿沈渣画像の領域分割が正確に行えない、という問題があった。背景領域に色むら、濃度むら（明暗のむら）が存在する場合、背景領域の濃度分布は特定の方向に広がる場合がある。例えばＲＧ平面上において、明るさ方向は原点から放射上に広がる方向であり、背景に濃度むらがある場合には図３の１０１に示すように背景濃度分布は原点と背景濃度分布の中心を結ぶ直線方向に広がった形状となる。この時、対象領域の濃度が１０２のように分布したとすると、Ｒ、Ｇ濃度に閾値を設けるだけで背景と対象とを正確に分離することはできない。例えば背景濃度分布１０１を囲むようにＲ濃度に閾値Ｔ２１、Ｔ２２を設定し、Ｇ濃度に閾値Ｔ３１、Ｔ３２を設定し、 Ｒ濃度がＴ２１以上Ｔ２２以下であり、かつＧ濃度がＴ３１以上Ｔ３２以下である濃度を有する画素を背景、それ以外の画素を対象領域とすると、対象領域の濃度分布１０２は背景濃度とみなされ、対象領域が正確に抽出されない。また、対象領域が抽出されるように閾値を設定すると、背景濃度の一部が対象領域の濃度とみなされ、背景の一部が対象として抽出され、正確な領域分割を行うことができない。特に背景と対象粒子との色調差が小さい画像の領域分割において、上記の問題が顕著に現われる。
【０００４】
本発明は上記事情に鑑みてなされたものであり、その目的とするところは、従来の技術における上述の如き問題を解消し、背景に色むらや濃度むらがある粒子画像や、背景と対象粒子との色調差が小さい粒子画像でも正確に背景と対象領域とを分離することができる領域分割方法を提供することにある。
【０００５】
【課題を解決するための手段】
本発明の上述の目的は、色濃度座標中で、粒子画像の背景の濃度分布が最も広がる方向と、色濃度座標の一軸が平行になるように色濃度座標を回転し、粒子画像の濃度値を回転後の座標の座標値に変換し、該座標値に閾値を設け、背景と対象領域とに分割することにより達成される。
【０００６】
また、特に背景の色が一様で、濃度むらのみがある場合には、ＲＧＢ濃度座標における背景の分布の中心（Ｒｍ、Ｇｍ、Ｂｍ）を求め、濃度座標の黒色を示す点（座標原点）と（Ｒｍ、Ｇｍ、Ｂｍ）を通る直線が、色濃度座標の一軸に平行になるように色濃度座標を回転し、Ｒ濃度、Ｇ濃度、Ｂ濃度を回転後の座標系（ＵＶＷ座標系とする）の座標値（Ｕ値、Ｖ値、Ｗ値）に変換し、Ｕ、Ｖ、Ｗ値に閾値を設け、背景と対象領域を分離することができる。図２に示すように背景濃度と対象領域濃度が分布した場合、原点と背景濃度分布の中心を通る軸を有する座標系は例えば図１に示すＵ軸、Ｖ軸により表現することができる。背景濃度、対象領域濃度をＵ値、Ｖ値により表現し、例えばＵ値に閾値Ｕａ、Ｕｃを設け、Ｖ値に閾値Ｖａ、Ｖｃを設け、Ｕ値がＵａ以上Ｕｃ以下であり、かつＶ値がＶａ以上Ｖｃ以下である画素を背景、それ以外の画素を対象領域とすることにより、背景と対象領域を正確に分離することができる。
【０００７】
背景の面積が大きい粒子画像の場合には、粒子画像のＲ濃度、Ｇ濃度、Ｂ濃度の最頻値が背景濃度分布の中心に一致する。そのため、粒子画像のＲ、Ｇ、Ｂ濃度の最頻値を上記（Ｒｍ、Ｇｍ、Ｂｍ）として用いることができる。
【０００８】
また、濃度座標の変換を行った後、更に画像を背景と対象領域と、背景濃度と対象領域濃度の中間的な濃度を持つ領域（中間領域）の３種類の領域に分割し、中間領域の周囲長Ｌと、背景と接する周囲長Ｌｂを求め、Ｌｂ／Ｌに閾値を設定し、Ｌｂ／Ｌが閾値より小さい中間領域を対象領域に統合し、Ｌｂ／Ｌが閾値より大きい中間領域を背景に統合することにより、背景と対象の色調差が小さい画像に対しても正確な領域分割を行うことができる。
【０００９】
また、濃度座標を回転する際の、回転角に関する情報を表示する工程により、操作者が回転角を確認し、背景の色調を知ることができ、背景の色調がおかしい、領域分割が正常に行われていない可能性がある、等の異常を検出することができる。
【００１０】
また、濃度座標の回転角に基準となる範囲を設け、回転角が基準となる範囲から外れた場合には警告を発する工程により、上記のような異常を自動検出することができる。
【００１１】
【発明の実施の形態】
以下、本発明を尿沈渣検査装置に応用した実施例を、図面を参照して詳細に説明する。
【００１２】
図４は本発明を応用した尿沈渣検査装置の構成の概要を示す図である。測定装置本体２００では染色液添加装置２１０により尿試料に染色液を添加し、一定時間後、画像入力装置２２０により尿中の固形成分の拡大静止画像を撮影する。撮影された画像はパターン認識装置２３０に転送され、画像パターン認識を行い、撮影された固形成分の分類が行われ、１つの検体中にどのような固形成分がどのくらいの頻度で出現したかがカウントされる。カウントされた結果は出力装置２４０を通じて操作者に報告される。出力装置２４０としては例えばプリンタ、液晶ディスプレイ等を用いる。画像入力装置２２０により撮影された画像、測定装置本体２００による分類結果、及び測定装置本体２００による画像パターン認識の途中で得られた画像特徴量、等のデータは通信手段２０を通じレビュー装置３０に転送される。レビュー装置３０は、操作者が任意の画像を表示させ、自動分類の修正や、目視による細分類を行うことができる。レビュー装置３０としてはディスプレイ３１、キーボード３２、記憶装置３３を備えた汎用のパーソナルコンピュータを用いる。
【００１３】
図５は測定装置本体２００を構成する画像入力装置２２０の構成例を示す図である。この装置では、フローセル２２１を用い、対物レンズ２２３とパルスランプ２２２との間に幅広で厚さの薄い、扁平な尿試料の流れを作る。フローセル２２１により形成される尿試料の流れには、パルスランプ２２２が瞬間的に照射され、対物レンズ２２３により拡大される尿中の固形成分像がＣＣＤカラーテレビカメラ２２４により静止画像として撮影される。得られた画像はパターン認識装置２３０に転送される。
【００１４】
次に画像処理装置２３０の、より詳細な構成例を図６を用いて説明する。ＣＣＤカラーテレビカメラ２２４より入力されるアナログ画像信号はＡ／Ｄ変換器２３１によりＲ濃度、Ｇ濃度、Ｂ濃度のデジタルデータに変換され、領域分割部２３２に送られる。領域分割部２３２は画像を背景領域と対象領域とに分割する。特徴量計算部２３３では領域分割部２３２で対象領域として分割された領域の画像特徴量を計算する。画像中に複数の対象が存在する時には、それぞれの対象領域について特徴量を計算する。特徴量としては例えば対象領域の面積、周囲長、平均色濃度、等を用いる。
【００１５】
パターン認識部２３４には特徴量計算部２３３で求められた各対象領域の画像特徴量が入力され、その対象がどのような成分であるか分類を行う。パターン認識としては例えばニューラルネットワークを用いることができる。分類項目としては例えば赤血球、白血球、上皮細胞、円柱、結晶等を用いる。
【００１６】
カウント部２３５は分類すべきクラスと同数のカウンタを持ち、パターン認識部２３４により各分類クラスに分類された対象の数をカウントする。一つの検体について測定が終了すると、カウンタの内容が単位体積中の個数に換算され、出力装置２４０に転送され、出力される。各カウンタの値は一検体の測定が終了するごとに０にリセットされる。
【００１７】
次に図６における領域分割部２３２の処理内容をより詳細に説明する。図７は本発明による領域分割部２３２の処理フローの第一の実施例を説明する図である。まず、ステップＳ１００において領域分割部に入力された粒子画像のＲ濃度、Ｇ濃度、Ｂ濃度のヒストグラムを作成する。次にステップＳ１１０においてＲ濃度、Ｇ濃度、Ｂ濃度の各ヒストグラムの最頻値を求め、Ｒｍ、Ｇｍ、Ｂｍとする。本装置で撮影された尿沈渣画像は背景領域が最も大きな面積を占めるため、濃度ヒストグラム中、最も大きなピークが背景濃度分布の中心を表す。次に、ステップＳ１２０で濃度座標の変換式を決定する。背景色が一様で、濃度むらのみが存在する場合には、背景濃度分布は黒色を示す点（原点）から放射状の方向へ広がる。そこで、変換式は、変換後の座標系のうちの一軸が、ＲＧＢ座標における原点（０、０、０）と背景濃度分布の中心を示す点（Ｒｍ、Ｇｍ、Ｂｍ）を結ぶ方向と平行になるように決定する。変換後の座標軸をＵ、Ｖ、Ｗとすると、例えば図８に示すようにＵ軸をＲＧＢ座標の原点と（Ｒｍ、Ｇｍ、Ｂｍ）を結ぶ方向に一致させるとする。また、Ｕ、Ｖ、Ｗ軸は例えばそれぞれ直交する座標系とする。この時の変換式はＶ、Ｗ軸の向きにより無数に存在し、一意には定まらない
【００１８】
が、例えば【数１】に示す式により変換を行うことができる。ただし、ｒ、ｇ、ｂはそれぞれＲ濃度値、Ｇ濃度値、Ｂ濃度値を表し、ｕ、ｖ、ｗは変換後のＵ値、Ｖ値、Ｗ値を表す。
【００１９】
【数１】
0004061760.tif 000002
【００２０】
ここで、θ、ωは図８に示す角度であり、ｓｉｎθ、ｃｏｓθ、ｓｉｎω、ｃ
【００２１】
ｏｓωはそれぞれ【数２】、【数３】、【数４】、【数５】により求めることができる。
【００２２】
【数２】
0004061760.tif 000003
【００２３】
【数３】
0004061760.tif 000004
【００２４】
【数４】
0004061760.tif 000005
【００２５】
【数５】
0004061760.tif 000006
【００２６】
なお、前述のように座標変換は一意には定まらず、ここで述べた変換方法は一例に過ぎない。他にも様々な変換式を用いることが可能である。また、各座標軸は必ずしも直交している必要はなく、極座標等も用いることができる。極座標（ｒ、θ、φ）を用いる場合にはｒ軸を背景濃度分布の広がり方向に設定する。
【００２７】
ステップＳ１２３ではステップＳ１２０で求めた濃度変換式により座標が回転する回転角度を記録する。回転角度は数１による変換式を用いる場合にはθ、ωに相当する。回転角度は通信手段２０（図４）を介してレビュー装置３０に送られ、記憶装置３３に記録される。操作者はディスプレイ３１とキーボード３２を用いて記録装置３３に記録された回転角度を参照することができる。これらの回転角度は画像の背景色の異常の検出等、メンテナンス時の情報として用いることができる。また、ステップＳ１２６では回転角度を予め設定した基準範囲と比較し、回転角度が基準範囲から外れた場合には、異常を示す警告のメッセージを出力する。例えば上記θに閾値θ１、θ２（θ１＜θ２）を設定し、ωに閾値ω１、ω２（ω１＜ω２）を設定し、θ＜θ１またはθ＞θ２またはω＜ω１またはω＞ω２の時に警告のメッセージを出力する。警告のメッセージは図４の出力装置２４０に表示するか、または通信手段２０、レビュー装置３０を介し、ディスプレイ３１に表示する。また、警告は音により出力しても良い。
【００２８】
次にステップＳ１３０において、画像上の各画素のＲＧＢ値よりＵＶＷ値を求め、Ｕ値、Ｖ値、Ｗ値のヒストグラムを作成し、ステップＳ１４０においてＵ値、Ｖ値、Ｗ値に設ける閾値を定める。閾値はヒストグラムを用いて例えば次のように定めることができる。まず、図９に示す最頻値Ｕｍを求め、Ｕｍにおける頻度Ｆｐｕを記憶し、次に頻度がＦｐｕの半分であるＵ値（Ｕｌ、Ｕｈ）を求める
【００２９】
。その後、求まったＵｍ、Ｕｌ、Ｕｈを元に、【数６】、【数７】、【数８】、
【００３０】
【数９】により閾値Ｕａ、Ｕｂ、Ｕｃ、Ｕｄを求める。ただし、Ｋｕ１、Ｋｕ２（Ｋｕ１＜Ｋｕ２）、Ｋｕ３、Ｋｕ４（Ｋｕ３＜Ｋｕ４）は実験的に求められる定数である。Ｖ値、Ｗ値に関しても全く同様にしてＶ値に設ける閾値Ｖａ、Ｖｂ、Ｖｃ、Ｖｄ、及びＷ値に設ける閾値Ｗａ、Ｗｂ、Ｗｃ、Ｗｄを求める。
【００３１】
【数６】
0004061760.tif 000007
【００３２】
【数７】
0004061760.tif 000008
【００３３】
【数８】
0004061760.tif 000009
【００３４】
【数９】
0004061760.tif 000010
【００３５】
ステップＳ１５０では、ステップＳ１４０で定めた閾値を用い、粒子画像の対象領域を抽出する。ステップＳ１５０における対象領域抽出処理フローの詳細を図１１に示す。まずステップＳ３００において、画像を背景、対象領域、中間領域の３種類の領域に分割する。中間領域とは背景と対象領域との中間的な濃度値をもつ領域であり、ステップＳ１４０で定めた閾値を用いて例えば下記のように３種類の領域を定める。但し、以下では注目する画素のＵ値、Ｖ値、Ｗ値をそれぞれｕ、ｖ、ｗとする
（背景）Ｕａ＜ｕ＜ＵｃかつＶａ＜ｖ＜ＶｃかつＷａ＜ｗ＜Ｗｃを満たす画素
（対象領域）ｕ＜Ｕｂまたはｕ＞Ｕｄまたはｖ＜Ｖｂまたはｖ＞Ｖｄまたはｗ＜Ｗｂまたはｗ＞Ｗｄを満たす画素
（中間領域）背景、対象領域以外の画素
これはＵＶＷ座標系において、濃度値がＵ＝Ｕａ、Ｕ＝Ｕｃ、Ｖ＝Ｖａ、Ｖ＝Ｖｃ、Ｗ＝Ｗａ、Ｗ＝Ｗｃという６個の平面に囲まれた直方体（直方体Ａとする）内にある画素を背景とし、Ｕ＝Ｕｂ、Ｕ＝Ｕｄ、Ｖ＝Ｖｂ、Ｖ＝Ｖｄ、Ｗ＝Ｗｂ、Ｗ＝Ｗｄという６個の平面に囲まれた直方体（直方体Ｂとする）の外にある画素を対象領域としている。直方体Ｂは直方体Ａを内部に含み、直方体Ｂに含まれ、直方体Ａの外にある画素が中間領域となる。わかり易くするため、ＵＶ平面でこの関係を図示したのが図１０である。図１０において外側の長方形が上述の直方体Ｂに相当し、内側の長方形が直方体Ａに相当する。内側の長方形の内部に存在する画素は背景、外側の長方形の外側に存在する画素は対象領域、２つの長方形に挟まれた領域に存在する画素が中間領域として分割される。
【００３６】
画像が３種類の領域に分割された後、ステップＳ３１０において、中間領域にラベリングを行う。ラベリングとは、抽出された領域が１個の連結領域でなく、複数であった場合、それぞれの領域に個別の番号を与える処理である。ここでは、中間領域として抽出された領域が複数であった場合に、それぞれの領域に１から始まる個別の番号を与える。
【００３７】
ラベリングが終了した後、ステップＳ３２０でｎに１を代入し、ステップＳ３３０で中間領域ｎ（ステップＳ３１０のラベリング処理においてｎという番号が与えられた中間領域）の周囲長Ｌを求める。さらにステップＳ３４０で中間領域ｎの周囲長のうち、背景領域に接する長さＬｂを求める。次にステップＳ３５０においてＬｂ／Ｌを予め定めた閾値Ｒｌと比較し、Ｌｂ／ＬがＲｌより大きい時にはステップＳ３６０で中間領域ｎを背景に統合し、Ｌｂ／ＬがＲｌより小さい時にはステップＳ３７０で中間領域ｎを対象領域に統合する。さらにステップＳ３８０でｎの値を１つ増加し、ステップＳ３９０で中間領域ｎの有無を調べ、中間領域ｎが存在すればステップＳ３３０以降の処理を繰り返す。中間領域ｎが存在しない場合には対象領域抽出処理を終了する。更に、必要に応じ、対象領域の穴埋め、背景領域の孤立点除去等の修正処理を施す。
【００３８】
ここで述べた、粒子画像を背景、対象領域、中間領域に分割した後、中間領域の周囲長Ｌ、対象領域と接する周囲長Ｌｂを用いて中間領域を背景、対象領域のどちらかに統合する、という対象領域抽出方法は、特に対象粒子と背景が似た色調の場合に有効な方法である。対象粒子が背景と似た色調の場合、背景、対象領域、中間領域の３種類の領域は例えば図１２に示すように分割される。図１２では、中間領域３２０は対象領域３１０の内部にあり、本来対象領域として抽出されるべき領域である。逆に中間領域３３０は対象領域３１０の外側にあり、本来背景として分割されるべき領域である。各中間領域の周囲長Ｌと背景と接する周囲長Ｌｂの比Ｌｂ／Ｌに着目すると、中間領域３２０のように対象領域の内部にある場合、背景に接している部分が無いかまたは非常に少ないので、Ｌｂ／Ｌは小さな値となる。一方、中間領域３３０のように対象領域の外側にある場合は、背景と接している部分が多いため、Ｌｂ／Ｌは大きな値となる。そのため、上述のようにＬｂ／Ｌが大きい中間領域を背景に統合し、Ｌｂ／Ｌの小さい領域を対象領域に統合することにより、対象領域の内部を誤って背景と分割したり、背景の一部を誤って対象領域として抽出することを防ぐことができる。
【００３９】
以上で述べたステップＳ１４０における閾値の設定方法、ステップＳ１５０における対象領域の抽出方法は一例であり、本発明に用いる閾値設定方法、対象領域抽出方法はこれに限定されるものではない。閾値設定方法、対象領域抽出方法は撮影される粒子の性質に応じ、適切な方法を選択すれば良い。
【００４０】
次に図１３を用いて本発明による領域分割部２３２（図６）の処理フローの第２の実施例を説明する。第１の実施例では、粒子が撮影された画像の濃度ヒストグラムから背景の平均濃度（Ｒｍ、Ｇｍ、Ｂｍ）を決定したが、第２の実施例では、粒子画像の撮影に先立ち、まず、ステップＳ２１０において粒子の存在しない背景のみの画像を撮影する。ステップＳ２２０でこの画像の濃度分布の広がりの方向を求める。濃度分布の広がりの方向を求めるためには、例えば多変量解析技術を用いることができる。多変量解析と呼ばれる技術のうち、主成分分析を行うことにより、分布が最も広がる方向を第１主成分として求めることができる。また、ステップＳ２３０において第２主成分、第３主成分を求めることにより、座標変換式を行列式として求めることができる。また、背景濃度の色が一定で、濃度のみが変化していることが予め分かっている場合には、ステップＳ２２０においてＲ濃度、Ｇ濃度、Ｂ濃度の平均値（または最頻値）を求め、それぞれＲｍ、Ｇｍ、Ｂｍとする。（Ｒｍ、Ｇｍ、Ｂｍ）は第１の実施例と同様、背景濃度分布の中心を表す。次にステップＳ２３０において、（Ｒｍ、Ｇｍ、Ｂｍ）を用いてＲ、Ｇ、Ｂ濃度値をＵ、Ｖ、Ｗ値に変換する変換式を決定する。変換式の決定方法は第１の実施例のステップＳ１２０と同じ方法を用いることができる。ステップＳ２３３ではステップＳ２３０で求めた変換式による座標が回転する回転角度を記録し、ステップＳ２３６で回転角度が異常でないかを判定する。ステップＳ２３３及びステップＳ２３６の処理はそれぞれ第１の実施例におけるステップＳ１２３、Ｓ１２６の処理と同じ処理を用いることができる。
【００４１】
変換式が決定した後、粒子画像の入力を開始し（ステップＳ２４０）、粒子画像のＵ、Ｖ、Ｗ値のヒストグラムを作成し（ステップＳ２５０）、閾値を決定し（ステップＳ２６０）、対象領域を抽出する（ステップＳ２７０）。ステップＳ２５０、Ｓ２６０、Ｓ２７０における処理は第１の実施例で述べたステップＳ１３０、Ｓ１４０、Ｓ１５０と同じ処理を用いることができる。ステップＳ２８０では領域分割結果を出力し、ステップＳ２９０で次の画像が入力されるかどうかを判断する。次の画像が入力される場合にはステップＳ２４０以下の処理を繰り返し、画像が入力されない場合には処理を終了する。
【００４２】
第１の実施例では粒子画像１枚毎に背景濃度を求め、最適な変換式を求めているのに対し、第２の実施例では最初に撮影した背景画像を用いて決定した変換式をその後の全ての粒子画像の領域分割に使用するため、第２の実施例の方が高速な処理が可能である。粒子画像１枚毎に背景の濃度が変動するような場合には第１の実施例が有効である。一方、背景の濃度が全ての画像において殆ど変化しない場合には、第２の実施例が有効である。また、第２の実施例において、主成分分析を行い背景の濃度分布の方向を求める場合には、処理時間は増えるが、背景に濃度むらだけでなく色むらがある場合にも望ましい領域分割結果を得る事が可能である。
【００４３】
以上、本発明を尿沈渣検査装置に応用した実施例について述べたが、本発明は尿沈渣検査装置への応用に限定されるものではない。例えば血球画像、微生物画像等、様々な粒子画像の領域分割に応用が可能である。
【００４４】
【発明の効果】
以上述べたように、本発明によれば粒子画像の背景濃度の広がりの方向に軸を有する濃度座標を定め、領域分割を行うことにより、背景に濃度むらや色むらがある画像についても正確な領域分割が可能である。背景濃度分布の広がりの方向は予め粒子の存在しない画像を撮影することにより、求めることができ、この方法は画像毎に背景濃度が変化しない場合や、背景に濃度むらだけでなく色むらがある場合に有効である。また、粒子が撮影された画像の濃度の最頻値から背景濃度の広がりの方向を求めることも可能であり、この場合、画像一枚毎に背景濃度分布の広がりの方向を求めることができるので、画像一枚毎に背景濃度や背景色が変化する場合に有効である。
【００４５】
更に、画像を背景と対象領域と、背景濃度と対象領域濃度との中間的な濃度を持つ領域（中間領域）の３種類に分割し、中間領域の周囲長Ｌと、背景と接する周囲長Ｌｂを求め、Ｌｂ／Ｌに閾値を設定し、Ｌｂ／Ｌが閾値より大きい中間領域を背景、閾値より小さい中間領域を対象領域に統合することにより、背景と対象との色調差が小さい画像に対しても正確な領域分割を行うことができる。
【００４６】
また、濃度座標を回転する際の回転角度を記録し、後から参照することにより、装置の異常を検出することが可能となる。予め回転角度に基準範囲を設定し、回転角度が基準範囲を外れた時に警告を発することにより、装置の異常を自動検出することができる。
【図面の簡単な説明】
【図１】本発明による領域分割の原理を示す図である。
【図２】従来技術による領域分割の原理を示す図である。
【図３】従来技術による領域分割の問題点を示す図である。
【図４】本発明を用いた尿沈渣検査装置の構成例を示す図である。
【図５】本発明を用いた尿沈渣検査装置の画像入力装置の構成例を示す図である。
【図６】本発明を用いた尿沈渣検査装置のパターン認識装置の構成例を示す図である。
【図７】本発明を用いた尿沈渣検査装置の領域分割部における処理フローの第１の実施例を示す図である。
【図８】変換前後の座標軸の関係を示す図である。
【図９】閾値設定基準を説明する図である。
【図１０】背景濃度と、対象領域濃度と、中間領域濃度の関係を説明する図である。
【図１１】本発明を用いた尿沈渣検査装置の領域分割部における領域抽出処理フローを示す図である。
【図１２】粒子画像上における、背景と対象領域と中間領域の関係を説明する図である。
【図１３】本発明を用いた尿沈渣検査装置の領域分割部における処理フローの第２の実施例を示す図である。
【符号の説明】
２０：通信手段
３０：レビュー装置
３１：ディスプレイ
３２：キーボード
３３：記憶装置
１０１：ＲＧ濃度平面上における背景の濃度分布
１０２：ＲＧ濃度平面上における対象の濃度分布
１０３：ＲＧ濃度平面上における背景の濃度分布
１０４：ＲＧ濃度平面上における対象の濃度分布
２００：測定装置本体
２１０：染色液添加装置
２２０：画像入力装置
２２１：フローセル
２２２：パルスランプ
２２３：対物レンズ
２２４：ＣＣＤカラーテレビカメラ
２３０：パターン認識装置
２３１：Ａ／Ｄ変換器
２３２：領域分割部
２３３：特徴量計算部
２３４：パターン認識部
２３５：カウント部
２４０：出力装置
３００：粒子画像
３１０：対象領域
３２０：中間領域
３３０：中間領域。
【図１】
0004061760.tif 000011
【図２】
0004061760.tif 000012
【図３】
0004061760.tif 000013
【図４】
0004061760.tif 000014
【図５】
0004061760.tif 000015
【図６】
0004061760.tif 000016
【図７】
0004061760.tif 000017
【図８】
0004061760.tif 000018
【図９】
0004061760.tif 000019
【図１０】
0004061760.tif 000020
【図１１】
0004061760.tif 000021
【図１２】
0004061760.tif 000022
【図１３】
0004061760.tif 000023