(57)【要約】  変化に富んだ性質をもつ尿試料において、尿中粒子を撮像した尿中粒子画像の対象粒子の領域を正確に抽出する方法と装置を提供する。  初めに、粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を抽出する。続いて、Ｒ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像の第１対象領域の色濃度分布と大きさを算出し、これらの特徴量に基いて、第１対象領域を所定数の群に分類する。群毎にＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を含む局所領域から第２対象領域を抽出する。この構成により、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割が行えるようになる。
2009125678.tif 000002
【特許請求の範囲】
【請求項１】
尿中粒子画像の領域分割方法であって、
粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像の赤成分の画像（Ｒ画像）、緑成分の画像（Ｇ画像）、青成分の画像（Ｂ画像）の中の１又は複数の画像を使用して、第１対象領域を抽出する工程と、
前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像の前記第１対象領域の色濃度分布と前記第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、前記第１対象領域を所定数の群に分類する工程と、
前記群毎に前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像を使用して、前記画像中の前記第１対象領域を含む局所領域から第２対象領域を抽出する工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項２】
尿中粒子画像の領域分割方法であって、
粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像を使用して、Ｒ濃度、Ｇ濃度、Ｂ濃度の濃度ヒストグラムを作成し、
前記濃度ヒストグラムの形状を表す、１又は複数のパラメータを求める第１の工程と、
前記１又は複数のパラメータと前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を抽出する第２の工程と、
前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像の前記第１対象領域の色濃度分布と前記第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、前記第１対象領域を所定数の群に分類する第３の工程と、
前記群毎に前記１又は複数のパラメータと前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像を使用して、前記画像中の前記第１対象領域を含む局所領域から第２対象領域を抽出する第４の工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項３】
請求項２記載の尿中粒子画像の領域分割方法であって、
前記第２の工程は、
前記第１の工程で求めた前記１又は複数のパラメータから、
前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像の濃度閾値を算出する工程と、
前記閾値を前記画像に適用して、対象領域を抽出する工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項４】
請求項２記載の尿中粒子画像の領域分割方法であって、
前記第２の工程は、
前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像をフィルタ処理して値を得る工程と、
前記値に閾値処理を行い対象領域を抽出する工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項５】
請求項２記載の尿中粒子画像の領域分割方法であって、前記第２の工程は、前記第１の工程で求めた前記１又は複数のパラメータから、前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像の濃度閾値を算出し、前記閾値を前記画像に適用して領域分割を行い抽出された領域分割結果と、前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像をフィルタ処理して値を得て、前記値に閾値処理を行い、領域分割を行い抽出された領域分割結果との論理和を取ることで抽出された領域を前記第１対象領域とすることを特徴とする粒子画像の領域分割方法。
【請求項６】
請求項２記載の尿中粒子画像の領域分割方法であって、
前記第４の工程は、
前記第３の工程で識別した群毎に、前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像を選択する工程と、
前記選択した色画像毎に、前記第１の工程で求めた前記１又は複数のパラメータを用いて閾値を算出する工程と、
前記閾値を利用して前記群毎に領域分割を行う工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項７】
請求項２記載の尿中粒子画像の領域分割方法であって、
前記第４の工程は、
前記第３の工程で識別した群毎に、前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像を選択する工程と、
前記選択した色画像毎に、前記第１の工程で求めた前記１又は複数のパラメータを用いて閾値を算出する工程と、
前記閾値を利用して前記群毎に領域分割を行い、第１領域分割結果を得る工程と、
前記第３の工程で識別した群毎に、前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の１又は複数の画像をフィルタ処理する工程と、
選択した前記フィルタ処理により得られる値に閾値処理を行い、第２領域分割結果を得る工程と、
前記第１領域分割結果と前記第２領域分割結果との論理和を取って、抽出された領域を前記第２対象領域とする工程と
を有することを特徴とする粒子画像の領域分割方法。
【請求項８】
請求項２記載の尿中粒子画像の領域分割方法であって、前記第３の工程で識別した群が所定の一群に分類された群に対しては、前記第４の工程で、前記第２の工程の領域分割結果を出力することを特徴とする粒子画像の領域分割方法。
【請求項９】
尿中粒子画像の領域分割装置であって、
粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を抽出する手段と、
前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像の前記第１対象領域の色濃度分布と前記第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、前記第１対象領域を所定数の群に分類する手段と、
前記群毎に前記Ｒ画像、前記Ｇ画像、前記Ｂ画像の中の１又は複数の画像を使用して、前記画像中の前記第１対象領域を含む局所領域から第２対象領域を抽出する手段と
を有することを特徴とする尿中粒子画像の領域分割装置。
【請求項１０】
尿を被検液とし、被検液を染色する染色手段と、該染色された被検液を撮像する手段と、撮像された画像から前記被検液中の有形成分を識別する識別手段とを備え、該識別手段は、前記画像における前記有形成分を含む第一領域を切り出し、該第一領域中の色濃度分布と該第一領域中の大きさを算出し、該算出された前記色濃度分布と前記大きさを基に前記有形成分を所定数分類することを特徴とする尿沈渣自動分析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、色濃度情報を用いた画像の領域分割方法に関し、特に尿中の粒子像を分割するのに好適な領域分割方法、及び装置に関する。
【背景技術】
【０００２】
尿中の粒子を形態学的に検査するには、従来、採取した尿を遠心分離し、検査技師が尿中粒子を直接顕微鏡観察する方法で検査を行っていた。顕微鏡検査は、（１）技師の習熟度に結果が左右される、（２）検査に時間を要する、等の問題があるため効率化が求められている。【０００３】
近年では、尿沈渣（尿中粒子）検査の自動化が進み、例えば、「粒子分析装置及び粒子分析方法（特許文献１，２）」では、尿試料を特別な形状の流路（フローセル）に通し、試料中の粒子を幅広の撮像領域中に流し、フラッシュランプを点灯して、尿中の粒子の拡大像を静止画像として撮影する方法が示されている。【０００４】
静止画像として撮像された尿中粒子を自動分類するためには、まず、画像上で尿中粒子の領域と背景領域とを分離した後、尿中粒子領域における画像特徴量を求め、これらの特徴量に基き分類を行う。【０００５】
画像上で尿中粒子の領域を背景領域から分離する従来技術としては、例えば「粒子画像の領域分割方法（特許文献３）」に記載がある。この従来技術では、濃度ヒストグラムにより求めた閾値と濃度値の変化の大きさとから、背景と対象領域とを分離する方法が示されている。【０００６】
また、画像特徴量から対象の分類を行う従来技術としては、例えば、「パターン認識装置（特許文献４）」や「血球自動分類装置（特許文献５）」に記載がある。この従来技術には、認識論理として階層型ネットワークを用いることの記載がある。
【特許文献１】特開平５−２９６９１５号公報
【特許文献２】特開昭６３−９４１５６号公報
【特許文献３】特開平８−１４５８７１号公報
【特許文献４】特開平１０−３０２０６７号公報
【特許文献５】特開平３−１３１７５６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
尿沈渣検査装置では、領域分割を行う対象粒子が、以下のように変化に富んだ性質をもっている。【０００８】
（１）尿沈渣検査では、粒子判定を容易に行うために、尿試料を染色することがあるが、充分に染色される（有色）粒子、染色されにくい（充分に染色されない（淡色）、ほとんど染色されない（無色））粒子が混在する。【０００９】
（２）同一種類の粒子でも、充分に染色される粒子と、染色されにくい粒子とがある。【００１０】
（３）色調の異なる複数個の粒子が、同一の画像内に存在する。【００１１】
以上のように、尿沈渣装置での対象粒子は、変化に富んだ性質をもつため、全ての対象に同一の領域分割方法を用いると、必ずしも正確な領域が抽出されない場合があると言う問題があった。【００１２】
例えば、前記の従来の領域分割方法（特許文献３）には、濃度ヒストグラムの最頻値を背景の平均濃度として背景濃度分布を推定し、この分布より暗い部分及び明るい部分の両方を対象領域として抽出する方法が開示されている。染色されにくい粒子は、対象の内部での光の屈折・反射により、対象領域中に背景より明るい濃度値をもつ画素が存在する。そのため、この方法によると、背景より明るい部分も対象として抽出することで、染色されにくい細胞の形状を正確に抽出できるとある。しかし、対象の輪郭付近での光の反射・屈折の影響により、対象に隣接する背景領域に、背景の平均濃度よりもかなり明るい部分が生じることがあり、そのため、輪郭近辺の背景領域を対象として抽出する場合があり、領域分割の精度低下の一因になっている。【００１３】
また、前記の従来の領域分割方法（特許文献３）を用いる場合、染色されにくい粒子の濃度ヒストグラムでは、対象領域と背景領域とが、一部重なっているため、これらの粒子の領域を抽出するためには、背景濃度分布を推定する際、背景領域とする濃度範囲を極力狭くして領域分割を行う必要がある。しかし、同じ閾値を用いて充分に染色される粒子の領域を抽出する場合、本来は背景である領域が対象として抽出される場合があり、領域分割の精度低下の一因になっている。【００１４】
そこで、本発明では、様々な大きさや色調の異なる尿中粒子が混在する尿試料において、粒子毎に安定した領域分割を行うための領域分割方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【００１５】
上記課題を解決するために、本発明の領域分割方法は、粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像の赤成分の画像（以下、Ｒ画像）、緑成分の画像（以下、Ｇ画像）、青成分の画像（以下、Ｂ画像）の中の１又は複数の画像を使用して、第１対象領域を抽出する工程と、Ｒ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像の第１対象領域の色濃度分布と第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、第１対象領域を所定数の群に分類する工程と、前記群毎にＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、画像中の第１対象領域を含む局所領域から第２対象領域を抽出する工程とを有することを特徴としている。この構成により、第１対象領域を粒子像より広めに抽出し、第１対象領域の色調と大きさから所定数の群に分類し、分類結果に基づき、粒子像の特徴に応じて、第２対象領域を抽出できるため、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割が行えるようになる。【００１６】
更に、本発明の領域分割方法は、粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像を使用して、Ｒ濃度、Ｇ濃度、Ｂ濃度の濃度ヒストグラムを作成し、濃度ヒストグラムの形状を表す、１又は複数のパラメータを求める第１の工程と、前記１又は複数のパラメータとＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を抽出する第２の工程と、Ｒ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像の第１対象領域の色濃度分布と第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、第１対象領域を所定数の群に分類する第３の工程と、前記群毎に前記１又は複数のパラメータとＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、画像中の第１対象領域を含む局所領域から第２対象領域を抽出する第４の工程とを有することを特徴としている。濃度ヒストグラムを用いた閾値処理を行うことで、色調が様々な尿中粒子が混在する尿試料においても、安定した閾値処理を行うことが可能となり、粒子像毎に安定した領域分割が行えるようになる。【００１７】
また、本発明の尿中粒子画像の領域分割装置は、粒子画像を入力する画像入力光学系により得た、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、第１対象領域を抽出する手段と、Ｒ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像の第１対象領域の色濃度分布と第１対象領域の大きさを算出し、該色濃度分布と該大きさに基づき、第１対象領域を所定数の群に分類する手段と、前記群毎にＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、画像中の第１対象領域を含む局所領域から第２対象領域を抽出する手段とを有することを特徴としている。
【発明の効果】
【００１８】
本発明の領域分割方法によれば、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割ができ、より正確な２値画像が得られる。その結果、対象領域の特徴量がより正確に求められ、対象粒子の分類の誤り防止が可能となり、各種の尿中粒子の識別精度が向上できる効果がある。
【図面の簡単な説明】
【００１９】
【図１】本発明が適用される尿沈渣検査装置の領域分割処理の一構成例を示す図。
【図２】尿沈渣検査装置の構成例を説明する図。
【図３】尿沈渣検査装置の画像入力装置の構成例を説明する図。
【図４】尿沈渣検査装置の画像処理方法を説明する図。
【図５】領域分割処理の詳細な一構成例を説明する図。
【図６】領域分割処理に画像毎の濃度ヒストグラムを用いた場合の詳細な構成を説明する図。
【図７】尿沈渣検査装置で得られる尿中粒子画像の一例を示す図。
【図８】画像毎の濃度ヒストグラムの形状を表すパラメータを求める詳細な手順を説明する図。
【図９】濃度ヒストグラムの一例を示す図。
【図１０】対象を広めにとる第１領域分割処理についての詳細な手順を示す図。
【図１１】第１対象領域を群分けする処理の詳細な手順を示す図。
【図１２】面積を用いた場合の群分け処理例のフローを説明する図。
【図１３】群分けに用いる色濃度値にＧ濃度とＢ濃度を選択した場合に用いる分布図例を説明する図。
【図１４】群分けに用いる色濃度値にＧ−Ｂ分布図を用いた場合の群分け例のフローを説明する図。
【図１５】群毎の第２領域分割処理についての詳細な手順を示す図。
【図１６】群毎に更なる領域分割の必要の有無を判断する処理例のフローを説明する図。
【図１７】群毎の第２領域分割処理についての詳細な手順を示す図。
【図１８】本発明の第２の実施例の尿中粒子画像の領域分割装置の構成を説明する図。
【符号の説明】
【００２０】
１０１  尿中粒子画像例
１０２  第１対象領域例
１０３  群分け例
１０４  第２対象領域例
２０１  測定装置本体
２０２  染色液添加装置
２０３  画像入力装置
２０４  ディスプレイ
２０５  キーボード
２０６  画像処理装置
３０１  フローセル
３０２  パルスランプ
３０３  対物レンズ
３０４  カメラ
７０１  尿中粒子画像例
７０２  対象粒子例
７０３  第１対象領域例
７０４  第１対象領域を含む局所領域矩形例
１３０１  識別境界線
１３０２  識別境界線
１３０３  Ａ群に属するシンボル
１３０４  Ｂ群に属するシンボル
１３０５  Ｃ群に属するシンボル
１３０６  Ｄ群に属するシンボル
【発明を実施するための最良の形態】
【００２１】
以下、図面を参照して本発明の実施の形態を説明する。
【実施例１】
【００２２】
図２は、本発明が適用される尿沈渣自動分析装置の構成例を説明する図である。測定装置本体２０１では、染色液添加装置２０２により尿試料に染色液を添加し、一定時間後、画像入力装置２０３により尿中の粒子の拡大静止画像を撮像する。撮像された画像は画像処理装置２０６に転送され、画像パターン認識により尿中粒子の分類を行い、１検体中に含まれる尿中粒子の種類とその出現頻度がカウントされる。画像処理装置２０６としては、例えばディスプレイ２０４、キーボード２０５を備えた汎用のパーソナルコンピュータを用いる。カウントされた結果はディスプレイ２０４を通じて操作者に報告される。画像入力装置２０３により撮像された画像、画像処理装置２０６による測定結果、各対象領域の分類結果及び画像パターン認識の途中で得られた画像特徴量等のデータは画像処理装置２０６内の記憶装置に保存される。また、画像処理装置２０６はレビュー機能も兼ね備えている。レビュー機能では、操作者が任意の画像を表示し、自動分類の修正、目視による細分類を実行できる。【００２３】
図３は、測定装置本体２０１を構成する画像入力装置２０３の構成図である。画像入力装置２０３では、フローセル３０１を用い、対物レンズ３０３とパルスランプ３０２との間に、幅広で厚さの薄い偏平な尿試料の流れを作る。フローセル３０１により形成される尿試料の流れに、パルスランプ３０２が瞬間的に照射され、対物レンズ３０３により拡大される尿中の粒子像がカメラ３０４により静止画像として撮像される。カメラ３０４としては、例えばＣＣＤカラーカメラ、ＣＭＯＳカラーカメラ等を用いる。得られた画像は画像処理装置２０６に転送される。【００２４】
図４は、画像処理装置２０６での画像処理方法を説明する図である。カメラ３０４により撮像された尿中粒子画像はデジタルデータとして、画像処理装置２０６に転送される。シェーディング補正処理ステップＳ４０１では、光学系の特性に由来する画像上の濃度むらを除去する。【００２５】
領域分割処理ステップＳ４０２では、尿中粒子が撮像された画像を背景領域と対象領域とに分割し、背景領域を０、対象領域を１とする２値画像を作成する。領域分割補正処理ステップＳ４０３では、対象領域の穴埋め、背景領域のノイズ除去等の２値画像の修正、整形を行う。手段としては、例えば、太め処理、細め処理等のフィルタ処理を含め、周知の従来技術を使用できる。【００２６】
ラベリング処理ステップＳ４０４では、２値画像の連結成分毎にラベリング処理を行い、画像中の複数の対象を一意に特定するための番号を付ける。特徴量計算処理ステップＳ４０５では、番号を付けたそれぞれの対象領域について、面積、周囲長、平均色濃度値等の特徴量を算出する。【００２７】
パターン認識処理ステップＳ４０６では、ステップＳ４０５で求めた各対象領域の画像特徴量を使用し、その対象がどのような成分であるか分類を行う。パターン認識処理手段としては、例えば、ニューラルネットワーク、統計的識別手法等を使用できる。分類項目としては、例えば、赤血球、白血球、扁平上皮細胞、その他の上皮細胞、円柱、結晶、細菌等を用いる。【００２８】
カウント処理ステップＳ４０７では、パターン認識処理ステップＳ４０６による分類結果に従い、各分類クラスに分類された対象数をカウントする。カウントした結果は、尿試料単位体積当たりの個数もしくは濃度に換算され、換算結果がディスプレイ２０４に出力される。【００２９】
なお、図４に示す処理の全て、もしくは一部をハードウェアで処理することも可能である。【００３０】
図５は領域分割処理ステップＳ４０２の詳細な一構成例を説明する図である。【００３１】
ステップＳ５０１では、尿中粒子が撮像された画像から、対象粒子が存在する第１対象領域を抽出する。ここでは、染色されにくい粒子の対象粒子を、連結した１つの領域として抽出することを目的として、対象粒子を広めにとる領域分割を行う。例えば、粒子画像を入力する画像入力光学系により、図１に示す尿中粒子を撮像した画像１０１を得たとする。この画像１０１のＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、予め実験的に定められた固定閾値等を用いて領域分割を行う。その結果、図１の１０２に示すように、対象粒子を広めにとる領域分割結果が得られる。【００３２】
ステップＳ５０２では、各第１対象領域を所定数の群に分ける。群分けには第１対象領域の特徴量を用いる。特徴量には、大きさに関する特徴量と色調に関する特徴量を用いる。大きさに関する特徴量には、例えば面積や周囲長等を用いる。色調に関する特徴量には、例えば、平均色濃度値等を用いる。例えば、図１の１０３では色調や大きさから第１対象領域を３つの群に分けた結果を示す。【００３３】
ステップＳ５０３では、群毎に詳細な領域分割を行うことで、第１対象領域を含む局所領域から第２対象領域を抽出する。領域分割方法には、例えば、尿中粒子を撮像した画像のＲ画像、Ｇ画像、Ｂ画像の中の１又は複数の画像を使用して、予め実験的に定められた固定閾値を用いて行う方法等を用いる。ここで、選択する画像の種類や数、固定閾値が群毎に異なる。例えば代表的なSternheimer(S)染色で染色された尿中粒子の場合、染色された粒子の吸収ピークは５５０ｎｍ付近にあるので、波長が４００〜５００ｎｍのＢ画像よりも、５００〜７００ｎｍのＧ画像、Ｂ画像を用いた方が感度が高い。また、染色されたものの中でも、充分に染色される有色粒子と、充分に染色されない淡色粒子とがあり、その場合は、同じ画像を用いたとしても、固定閾値を色調に応じて変える必要がある。また、領域分割方法には、動的輪郭（スネーク）といった周知の従来技術である輪郭抽出法を使用しても良い。例えば、図１の１０４では、群毎に異なる領域分割を行った結果、正確な対象領域が抽出された結果を示す。【００３４】
この構成により、第１対象領域を粒子像より広めに抽出し、第１対象領域の色調と大きさから所定数の群に分類し、分類結果に基づき、粒子像の特徴に応じて、第２対象領域を抽出できるため、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割が行えるといった効果が得られる。【００３５】
図６は領域分割処理ステップＳ４０２に尿中粒子画像のＲ濃度、Ｇ濃度、Ｂ濃度の濃度ヒストグラムを用いた場合の詳細な構成を説明する図である。【００３６】
ステップＳ６０１では、尿中粒子画像のＲ濃度、Ｇ濃度、Ｂ濃度の濃度ヒストグラムをそれぞれ作成し、各濃度ヒストグラムの形状を表すパラメータを求める。ステップＳ６０２では、ステップＳ６０１で求めたパラメータ等を用いて、画像から対象粒子が存在する第１対象領域を抽出する。ここでは、染色されにくい粒子の対象粒子を、連結した１つの領域として抽出することを目的として、対象粒子を広めにとる領域分割を行う。例えば図７に示すような尿中粒子画像７０１を得た場合、ステップＳ６０２によって得られる第１対象領域は、対象粒子７０２を広めにとった領域７０３になる。【００３７】
図６の群分け処理ステップＳ６０３では第１対象領域の平均色濃度値分布と、大きさから、第１対象領域を所定数の群に分ける。なお、大きさには面積や周囲長等を用いる。第２領域分割処理ステップＳ６０４では、分けられた群ごとに最適な異なる領域分割方法を用いることで、第１対象領域を含む局所領域から、対象粒子の形状がより正確な第２対象領域を抽出する。抽出された第２対象領域はステップＳ４０３で補正処理される。【００３８】
図８は、図６のステップＳ６０１で行われる、画像毎に濃度ヒストグラムの形状を表すパラメータPd(*)、Phl(*)、Phh(*)、dl(*)、dh(*)を求める詳細な手順を示す。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００３９】
ステップＳ８０１では、各画像において、Ｒ画像、Ｇ画像、Ｂ画像の濃度ヒストグラムを作成する。ステップＳ８０２では、各濃度ヒストグラムで、図９に示すように頻度の最大値Pmax(*)をもつ濃度値Pd(*)を求める。*は、Ｒ，Ｇ，Ｂのいずれかを表す。図８のステップＳ８０３では、図９に示すように、各濃度ヒストグラムで、ピークの半値幅Pmax/2(*)を与える濃度、Phl(*)、Phh(*)を求める。*は、Ｒ，Ｇ，Ｂのいずれかを表す。なお、ここではピークの半値を示す濃度値を用いたが、これに限ることはなく、ピークの１／４、１／１０等を示す濃度値でもよい。ヒストグラムの形状から最適な濃度値を用いるようにする。図８のステップＳ８０４では、ステップＳ８０２で求めた図９のPd(*)と、ステップＳ８０３で求めた図９のPhl(*)、Phh(*)を使用して、次の式（１）、式（２）により、図９のdl(*)とdh(*)を算出する。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００４０】
dl(*)＝Pd(*)−Phl(*)          (1)
dh(*)＝Phh(*)−Pd(*)          (2)
図１０は、図６のステップＳ６０２で行われる、対象を広めにとる第１領域分割処理についての詳細な手順を示す。【００４１】
はじめにステップＳ１００１で、濃度ヒストグラムを用いて、対象領域を抽出する手順を示す。ステップＳ１００２では、図８のステップＳ８０２で求めたPd(*)とステップＳ８０４で求めたdl(*)とdh(*)を使用して、次の式（３）、式（４）により、図９に示す閾値T1(*)，T2(*)を算出する。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００４２】
T1(*)＝Pd(*)−dl(*)×k1(*)          (3)
T2(*)＝Pd(*)＋dh(*)×k2(*)          (4)
k1(*)、k2(*)は予め実験的に定める係数で、Ｒ画像、Ｇ画像、Ｂ画像の色毎に異なる最適値を求めておく。なお、閾値の算出には、尿中粒子の色特性とカメラの分光特性から最も感度の高い色濃度値を１つもしくは複数選択して使用する。【００４３】
例えば、Ｒ画像とＧ画像を選択した場合は、以下のようになる。【００４４】
T1(R)＝Pd(R)−dl(R)×k1(R)
T2(R)＝Pd(R)＋dh(R)×k2(R)
T1(G)＝Pd(G)−dl(G)×k1(G)
T2(G)＝Pd(G)＋dh(G)×k2(G)
k1(R)、k2(R）、k1(G）、k2(G)は予め実験的に定める係数で、色毎に異なる最適値を求めておく。【００４５】
ステップＳ１００３では、求めた閾値を使用して対象領域を抽出する。例えば、閾値T1(R)、T2(R)、T1(G)、T2(G)を用いる場合には、次式（５）を満たす画素(x,y)を対象領域として抽出する。【００４６】
{T1(R)≧Ｒ(x,y)}II{T2(R)＜Ｒ(x,y)}II{T1(G)≧Ｇ(x,y)}II{T2(G)＜Ｇ(x,y)}  (5)
ここで、IIは論理和を表す。Ｒ(x,y)は、画素(x,y)におけるＲ濃度値を、Ｇ(x,y)は、画素(x,y)におけるＧ濃度値を表す。【００４７】
図６のステップＳ６０２では、第１対象領域を広めに抽出することが目的であるため、係数kn(*)はやや小さめに設定する。係数kn(*)を小さめに設定することで、図９の対象領域を広げることになる。【００４８】
続けて、ステップＳ１００４で、濃度値の変化の大きさを用いて、対象領域を抽出する手順を示す。ステップＳ１００５では、濃度値の変化の大きさを表す値を算出する。【００４９】
ここでは、濃度値の変化の大きさを表す指標として、局所的な小領域における濃度値の差分を使用する場合について説明する。ステップＳ１００５で局所的な小領域における濃度値の差分を使用する場合は、次の式（６）で定義される濃度差分値を用いる。画像上の画素の位置(x,y)における濃度値を*(x,y)、濃度差分値をr(*)(x,y)とすると、r(*)(x,y)は、
【数１】
2009125678.tif 000003
【００５０】
により、得られる。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００５１】
ステップＳ１００６では、r(*)(x,y)から、次の式（７）を満たす、画素(x,y)を対象領域として抽出する。【００５２】
sn (*)≦|r(*)(x,y)|          (7)
ここで、｜｜は絶対値を表し、sn(*)は予め実験的に定める定数で、色毎に異なる最適値を求めておく。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００５３】
ある画素の差分値を計算する際に用いる近傍画素の個数は２ｎ＋１で表され、この個数をマスクサイズと呼ぶ。例えば、ｎ＝２の時、ある画素の差分値の計算には、その画素を中心とし、ｘ方向もしくはｙ方向の前後２画素が必要であり、マスクサイズ＝５である。マスクサイズ＝１の場合は、差分処理を行わないことを意味する。【００５４】
なお、領域分割は、尿中粒子の色特性とカメラの分光特性から最も感度の高い色画像を１つもしくは複数選択して使用する。また、ｘ方向、ｙ方向と各マスクサイズは、無色粒子の対象領域内で、できるだけ大きな値を取るように設定する。【００５５】
例えば、色画像はＧ画像を選択し、ｘ方向のマスクサイズ＝５、ｙ方向のマスクサイズ＝１を選択した場合は、次の式（８）を満たす、画素(x,y)が対象領域となる。【００５６】
r(G)(x,y)=Ｇ(x+2,y)+Ｇ(x+1,y)-Ｇ(x-2,y)-Ｇ(x-1,y)
s1(G)≦|r(G)(x,y)|                            (8)
｜｜は絶対値を、定数s1(G)は予め実験的に最適値を求めておく。【００５７】
なお、濃度値の変化の大きさを表すのには、濃度値の差分を用いた方法に限定されることはなく、局所的な小領域における濃度値の分散を用いる方法、対象領域に特異的に含まれる周波数成分を強調するフィルタ処理を用いる方法、等他の方法を用いても良い。【００５８】
ステップＳ１００７では、ステップＳ１００１とステップＳ１００４で求めた対象領域を重ね合わせて（論理和を取って）、第１対象領域を抽出する。例えば、ステップＳ１００１ではＲ画像とＧ画像を用いて、式（５）より対象領域を抽出し、ステップＳ１００４では、Ｇ画像でｘ方向のマスクサイズ＝５、ｙ方向のマスクサイズ＝１を組合わせて、式（８）より対象領域を抽出した場合は、次式を満たす画素(x,y)が、第１対象領域として求められる。ここで、IIは論理和を表す。【００５９】
{T1(R)≧Ｒ(x,y)}II{T2(R)＜Ｒ(x,y)}II{T1(G)≧Ｇ(x,y)}II{T2(G)＜Ｇ(x,y)}
IIs1(G)≦|r(G)(x,y)|
ステップＳ１００８では、抽出された第１対象領域の領域分割補正処理を行う。方法は、図４のステップＳ４０３と同様である。ステップＳ１００９では、各第１対象領域にラベリング処理を行う。方法は、図４のステップＳ４０４と同様である。ステップＳ１０１０では、番号を付けた各第１対象領域につき、面積、周囲長、平均色濃度値等の特徴量を算出する。算出した特徴量を元に、図６のステップＳ６０３で群分け処理を行う。【００６０】
図１１は、図６のステップＳ６０３で行われる群分け処理の詳細な手順を示す。ステップＳ１１０１では、各第１対象領域の大きさに関する特徴量から、第１対象領域を大型粒子群と小型粒子群の２群に分ける。大きさに関する特徴量には、面積、周囲長等を用いる。【００６１】
図１２は、例えば、特徴量に面積を用いた場合の群分け処理のフローを説明する図である。ある第１対象領域の面積をＭとすると、ステップＳ１２０１の判定で、Ｍ≧ｍを満たす第１対象領域は大型粒子とし、満たさない場合は小型粒子と分類する。ここで、ｍは予め実験的に最適値を求めておく。【００６２】
図１１のステップＳ１１０２では、各第１対象領域の平均色濃度値から、第１対象領域を色調に応じた所定数の群に分ける。平均色濃度値には、対象領域の色特性とカメラの分光特性から最も感度の高い色画像を１つもしくは複数選択する。選択した平均色濃度値が、色濃度空間のどこに位置するのかで、第１対象領域を所定の群に分ける。ここで、色濃度空間で群に分けるための識別境界線は、予め実験的に定めておく。識別境界線を求める手段としては、例えば、判別分析等、周知の従来技術を使用できる。【００６３】
図１３は、例えば、色濃度値にＧ濃度とＢ濃度を用いた場合の群分け処理に用いる分布図例を説明する図である。シンボル１３０３〜１３０６は、予め分布図を作成するために集められた尿中粒子の平均濃度値座標（ｇ，ｂ）を示す。集められた尿中粒子は、最適な領域分割方法が異なることが実験的に明確にされている。ここでは、尿中粒子を色調により、４種類の群Ａ，Ｂ，Ｃ，Ｄに分けるものとし、シンボル１３０３〜１３０６は、それぞれ、群Ａ，Ｂ，Ｃ，Ｄに属する粒子のＧ−Ｂ空間中での分布を表すものとする。識別境界線１３０１と１３０２は、例えば、線形判別分析を用いて算出する。線形判別分析には、前記シンボルの濃度値データを用いる。算出された識別境界線は次式で示せる。【００６４】
ｂ＝ｊ1ｇ＋ｋ1  （識別境界線１３０１）
ｂ＝ｊ2ｇ＋ｋ2  （識別境界線１３０２）
図１４は、例えば、図１１のステップＳ１１０２の群分け処理に図１３のＧ−Ｂ分布図を用いた場合の群分けのフローを説明する図である。【００６５】
ある第１対象領域のＧ画像、Ｂ画像の色濃度平均値を（ｇ，ｂ）とすると、ステップＳ１４０１の判定で、ｂ≧ｊ1ｇ＋ｋ1を満たし、なおかつステップＳ１４０２でｂ≧ｊ2ｇ＋ｋ2を満たす場合は、第１対象領域は淡色の粒子群Ｂ（１３０４）に属する。ステップＳ１４０１の判定で、ｂ≧ｊ1ｇ＋ｋ1を満たし、なおかつステップ１４０２でｂ≧ｊ2ｇ＋ｋ2を満たさない場合は、第１対象領域は無色の粒子群Ｄ（１３０６）に属する。【００６６】
ステップＳ１４０１の判定で、ｂ≧ｊ1ｇ＋ｋ1を満たさず、なおかつステップ１４０３でｂ≧ｊ2ｇ＋ｋ2を満たす場合は、第１対象領域は有色の粒子群Ａ（１３０３）に属する。ステップＳ１４０１の判定で、ｂ≧ｊ1ｇ＋ｋ1を満たさず、なおかつステップ１４０３でｂ≧ｊ2ｇ＋ｋ2を満たさない場合は、第１対象領域は淡色の粒子群Ｃ（１３０５）に属する。【００６７】
図１１のステップＳ１１０３では、ステップＳ１１０２と同様の方法で群分け処理を行う。但し、分類する群の数や分類方法（識別境界線）は、予め実験的に別途設定しておく。【００６８】
図１５は、図６のステップＳ６０４で行われる、群毎の第２領域分割処理についての詳細な手順を示す。【００６９】
ステップＳ１５０１では、群毎に更なる領域分割が必要か否かの判断を行う。図１６にステップＳ１５０１の詳細なフローを示す。例えば、群Ｎは更なる領域分割が必要でないと予め実験的に定められている群とする。例えば、ある第１対象領域ｘが属する群Ｘが、ステップＳ１６０１で、Ｘ＝Ｎを満たす場合は、第１対象領域ｘは更なる領域分割を必要とせずに、第１対象領域を最終的な領域分割結果とする。ステップＳ１６０１で、Ｘ＝Ｎを満たさない場合は、第１対象領域ｘは更なる領域分割を必要とするので、ステップＳ１５０２に進む。【００７０】
群Ｎは例えば、大型粒子である。大型粒子の場合、図６のステップＳ６０２で対象を広めにとる領域分割を行ったとしても、元々大型であるため、図４のステップＳ４０５で算出される特徴量に大きな差は生じない。そのため、図４のステップＳ４０６でのパターン認識処理も良好に行われ、最終的な尿中粒子の誤分類も生じにくいからである。また、更なる領域分割を必要としないことにより、領域分割処理に要する時間を節約できることになるため、画像処理の処理速度が向上する効果がある。【００７１】
図１５のステップＳ１５０２では、第１対象領域を含む局所領域を抽出する。局所領域としては、第１対象領域を含む局所的な矩形を抽出するようにする。例えば図７に示すような尿中粒子画像７０１を得た場合、図６のステップＳ６０２によって得られた第１対象領域７０３を含む局所領域として矩形７０４を抽出する。【００７２】
ステップＳ１５０３では、群毎に異なる、群特異的な領域分割を行うことで、第１対象領域を含む局所領域から、より正確な第２対象領域を抽出する処理を行う。【００７３】
図１７は、ステップＳ１５０３で行われる、第２領域分割処理についての詳細な手順を示す。ステップＳ１７０１は基本的に図１０のステップＳ１００１と同様だが、ステップＳ１７０２で用いる色画像やパラメータが群毎に異なる。また、閾値の算出に用いる係数kn(*)も群毎に異なるため、閾値の算出方法も群毎に異なる。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００７４】
ステップＳ１７０４は基本的に図１０のステップＳ１００４と同様だが、ステップＳ１７０５で用いる色画像や定数sn(*)、マスクサイズが群毎に異なる。*は、Ｒ，Ｇ，Ｂのいずれかを表す。【００７５】
ステップＳ１７０７では、図６のステップＳ６０４で抽出される第２対象領域を抽出する。第２対象領域は、ステップＳ１７０３とステップＳ１７０６で求めた対象領域の重ね合わせ（論理和を取る）で求める。例えば、予め実験的に分けた群を、有色粒子群（１）、淡色粒子群（２）、無色粒子群（３）の３群とする。【００７６】
ある第１対象領域が、有色粒子群（１）に属する場合は、ステップＳ１７０１ではＲ画像とＧ画像を用いて対象領域を抽出し、ステップＳ１７０４では、Ｇ画像でｘ方向のマスクサイズ＝５、ｙ方向のマスクサイズ＝１を組合わせて、対象領域を抽出したとすると、第２対象領域は、次式により得られる。【００７７】
{T3(R)≧Ｒ(x,y)}II{T4(R)＜Ｒ(x,y)}II{T3(G)≧Ｇ(x,y)}II{T4(G)＜Ｇ(x,y)}
IIs2(G)≦|r(G)(x,y)|
ある第１対象領域が、淡色粒子群（２）に属する場合は、ステップＳ１７０１ではＲ画像とＧ画像を用いて対象領域を抽出し、ステップＳ１７０４では、ｘ，ｙ両方向に対して、マスクサイズ＝１を用いたとすると、第２対象領域は、次式により得られる。【００７８】
{T5(R)≧Ｒ(x,y)}II{T6(R)＜Ｒ(x,y)}II{T5(G)≧Ｇ(x,y)}II{T6(G)＜Ｇ(x,y)}
ある第１対象領域が、無色粒子群（３）に属する場合は、ステップＳ１７０１ではＢ画像を用いて対象領域を抽出し、ステップＳ１７０４では、ｘ，ｙ両方向に対して、マスクサイズ＝１を用いたとすると、第２対象領域は、次式により得られる。【００７９】
{T1(B)≧B(x,yj)}II{T2(B)＜B(x,y)}
本発明の濃度ヒストグラムを用いて閾値処理を行う領域分割方法により、色調が様々な尿中粒子が混在する尿試料においても、安定した閾値処理を行えるといった効果が得られる。【００８０】
また、第１対象領域を粒子像より広めに抽出し、第１対象領域の色調と大きさから所定数の群に分類し、分類結果に基づき、粒子像の特徴に応じて、第２対象領域を抽出できるという本発明の構成により、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割が行えるといった効果が得られる。その結果、対象領域の特徴量がより正確に求められ、対象粒子の分類の誤り防止が可能となり、各種の尿中粒子の識別精度が向上できる効果がある。
【実施例２】
【００８１】
図１８は、本発明の第２の実施例の尿中粒子画像の領域分割装置の構成を説明する図である。【００８２】
カメラ等の画像入力装置に撮像された元画像はメモリ１８０１に転送される。元画像は第１領域分割手段１８０２に転送され、第１領域分割が行われる。第１領域分割方法は実施例１に示した方法を用いれば良い。領域分割された第１対象領域画像はメモリ１８０１に転送される。【００８３】
第１対象領域は次に、群分け手段１８０３に転送され、第１対象領域毎に特徴量が算出され、所定数の群に分けられる。群分け処理方法は実施例１に示した方法を用いれば良い。群分け処理結果はメモリ１８０１に転送される。【００８４】
群毎の詳細な第２領域分割手段１８０４では、メモリ１８０１に保存してある、元画像、第１対象領域、群分け処理結果を用いて、領域分割を行い、対象粒子の形状に正確な第２対象領域を抽出する。群毎の詳細な領域分割方法は実施例１に示した方法を用いれば良い。第２対象領域結果はメモリ１８０１に転送される。ただし、第２領域分割結果を保存する手段としては、メモリ１８０１に限らず、例えば、ＨＤＤもしくはフロッピディスクといった外部記憶媒体を用いても良い。【００８５】
第１対象領域を粒子像より広めに抽出し、第１対象領域の色調と大きさから所定数の群に分類し、分類結果に基づき、粒子像の特徴に応じて、第２対象領域を抽出できるという、本発明の装置構成により、様々な大きさや色調の異なる尿中粒子が混在する尿試料においても、粒子像毎に安定した領域分割が行えるといった効果が得られる。
【図１】
2009125678.tif 000004
【図２】
2009125678.tif 000005
【図３】
2009125678.tif 000006
【図４】
2009125678.tif 000007
【図５】
2009125678.tif 000008
【図６】
2009125678.tif 000009
【図７】
2009125678.tif 000010
【図８】
2009125678.tif 000011
【図９】
2009125678.tif 000012
【図１０】
2009125678.tif 000013
【図１１】
2009125678.tif 000014
【図１２】
2009125678.tif 000015
【図１３】
2009125678.tif 000016
【図１４】
2009125678.tif 000017
【図１５】
2009125678.tif 000018
【図１６】
2009125678.tif 000019
【図１７】
2009125678.tif 000020
【図１８】
2009125678.tif 000021
【国際調査報告】
2009125678.tif 000022
2009125678.tif 000023