(57)【要約】【課題】精度よくコロニー数を計数する。【解決手段】培養されたコロニーを含む培地領域及びフィルタに記載される格子を含む画像から培地領域の画像を抽出する領域抽出部と、培地領域の画像を所定の大きさに調整する補正部と、培地領域内の格子を検出する格子検出部と、格子検出部で検出された格子を除く培地領域の画像の色を所定の基準の色に調整する色調整部と、色調整部で調整された培地領域の画像の各画素の色について、赤色成分と青色成分との差分値を求め、差分値が周囲の画素の差分値よりも大きい領域を抽出する色検出部と、色調整部で調整された培地領域の画像について、円形の領域を抽出する円形抽出部と、色検出部による抽出結果と円形抽出部による抽出結果とを統合する統合部と、統合部によって統合された抽出結果に基づいて、抽出結果に含まれる領域の数をコロニーとして計数する計数部と、を備えるコロニー数計数装置とする。【選択図】図１
2015154718.tif 000002
【特許請求の範囲】
【請求項１】
培養されたコロニーを含む培地領域を含む画像から前記培地領域の画像を抽出する領域抽出部と、
前記培地領域の画像の色を所定の基準の色に調整する色調整部と、
前記色調整部で調整された前記培地領域の画像の各画素の色について、赤色成分と青色成分との差分値を求め、前記差分値が周囲の画素の前記差分値よりも大きい領域を抽出する色検出部と、
前記色検出部によって検出された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する計数部と、
を備えるコロニー数計数装置。
【請求項２】
前記画像は、前記培地領域内に格子を含み、
前記培地領域内の前記格子を検出する格子検出部と、を備え、
前記色調整部は、前記格子検出部で検出された前記格子を除く前記培地領域の画像の色を所定の基準の色に調整する
請求項１に記載のコロニー数計数装置。
【請求項３】
前記色調整部で調整された前記培地領域の画像について、円形の領域を抽出する円形抽出部と、
前記色検出部による抽出結果と前記円形抽出部による抽出結果とを統合する統合部とを備え、
前記計数部は、前記統合部によって統合された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する
請求項１または２に記載のコロニー数計数装置。
【請求項４】
前記統合部によって統合された抽出結果から、前記格子の格子点である領域を除去した結果を抽出結果とする除去部を備え、
前記計数部は、前記除去部による抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する
請求項３に記載のコロニー数計数装置。
【請求項５】
前記培地領域の画像を所定の大きさに調整する補正部を備える
請求項１から４のいずれか１項に記載のコロニー数計数装置。
【請求項６】
コンピュータが、
培養されたコロニーを含む培地領域を含む画像から前記培地領域の画像を抽出し、
前記培地領域の画像の色を所定の基準の色に調整し、
調整された前記培地領域の画像の各画素の色について、赤色成分と青色成分との差分値を求め、前記差分値が周囲の画素の前記差分値よりも大きい領域を抽出し、
検出された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する、
ことを実行するコロニー数計数方法。
【請求項７】
前記画像は、前記培地領域内に格子を含み、
前記コンピュータが、
前記培地領域内の前記格子を検出し、
検出された前記格子を除く前記培地領域の画像の色を所定の基準の色に調整する、
ことを実行する請求項６に記載のコロニー数計数方法。
【請求項８】
前記コンピュータが、
調整された前記培地領域の画像について、円形の領域を抽出し、
前記差分値が周囲の画素の前記差分値よりも大きい領域を抽出した抽出結果と前記円形の領域を抽出した抽出結果とを統合し、
統合された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する
ことを実行する請求項６または７に記載のコロニー数計数方法。
【請求項９】
前記コンピュータが、
統合された抽出結果から、前記格子の格子点である領域を除去した結果を抽出結果とし、
当該抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する請求項８に記載のコロニー数計数方法。
【請求項１０】
前記コンピュータが、
前記培地領域の画像を所定の大きさに調整する
請求項６から９のいずれか１項に記載のコロニー数計数方法。
【請求項１１】
コンピュータに、
培養されたコロニーを含む培地領域を含む画像から前記培地領域の画像を抽出し、
前記培地領域の画像の色を所定の基準の色に調整し、
調整された前記培地領域の画像の各画素の色について、赤色成分と青色成分との差分値を求め、前記差分値が周囲の画素の前記差分値よりも大きい領域を抽出し、
検出された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する、
ことを実行させるコロニー数計数プログラム。
【請求項１２】
前記画像は、前記培地領域内に格子を含み、
前記コンピュータに、
前記培地領域内の前記格子を検出し、
検出された前記格子を除く前記培地領域の画像の色を所定の基準の色に調整する、
ことを実行させる請求項１１に記載のコロニー数計数プログラム。
【請求項１３】
前記コンピュータに、
調整された前記培地領域の画像について、円形の領域を抽出し、
前記差分値が周囲の画素の前記差分値よりも大きい領域を抽出した抽出結果と前記円形の領域を抽出した抽出結果とを統合し、
統合された抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数する
ことを実行させる請求項１１または１２に記載のコロニー数計数プログラム。
【請求項１４】
前記コンピュータに、
統合された抽出結果から、前記格子の格子点である領域を除去した結果を抽出結果とし、
当該抽出結果に基づいて、当該抽出結果に含まれる領域の数をコロニーとして計数することを実行させる請求項１３に記載のコロニー数計数プログラム。
【請求項１５】
前記コンピュータに、
前記培地領域の画像を所定の大きさに調整する
ことを実行させる請求項１１から１４のいずれか１項に記載のコロニー数計数プログラム
。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コロニー数計数装置、コロニー数計数方法、及びコロニー数計数プログラムに関する。
【背景技術】
【０００２】
水溶液の清浄度の確認が必要な分野としては、医療分野（製薬会社、透析用水など）、精密機械分野（半導体の洗浄用水など）、飲料水（水道局、飲料製造会社など）などがあり、水の清浄度管理を行われている。透析液や水道水などの清浄度を検査する手段として、生菌数試験が通常行われている。例えば透析液においては、２０１０年の診療報酬改正において透析液水質基準確保加算が算定されるなど、その清浄度を定期測定し、定量的に示すことが求められており、生菌数試験は重要な役割を果たしている。
【０００３】
透析液などの水溶液の生菌数試験としては、平板培養法やメンブレンフィルタ法（ＭＦ法：Membrane法）がある。透析医療施設で適用可能な透析液の生菌数試験としては、平板培養法やメンブレンフィルタ法（ＭＦ法）がある。ＭＦ法は各透析装置や透析液供給装置、水処理装置において水溶液をメンブレンフィルタに通過させて生菌を採取する。採取したメンブレンフィルタは培養液につけて生菌を培養する。また平板法では寒天培地に水溶液を付けて生菌を培養する。これらの方法で４日から７日の培養後、形成されたコロニー数のカウントすることで、生菌数はカウントされる。また、水溶液以外の生菌検査においても、平板法と同様に寒天培地に菌を付着させ、培養してコロニー数をカウントする。以上のように、水溶液検査などにおいて、生菌検査は多数の施設においてこれらの事項が行われている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−２２３４６３号公報
【特許文献２】特開２０１２−０７５４０９号公報
【特許文献３】特開２０１１−２１２０１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
透析医療施設等で適用可能な水溶液の生菌数試験としては、平板培養法やメンブレンフィルタ法（ＭＦ法）がある。各種水処理装置から採取したサンプルは、これらの方法で４日から７日の培養後、形成されたコロニー数をカウントする。現在、透析液のコロニーカウントは、肉眼で行われることで清浄度を判定し、さらには形成されたコロニーの形状や色調も記録する。このような生菌数試験が定期的に行われ、こちらに費やす労力は非常に大きい。
【０００６】
生菌数のカウントでは、Ｒ２Ａ寒天培地が広く用いられている。Ｒ２Ａ 寒天培地は１
ｍＬ中の水生菌を培養している。一方、メンブレンフィルタは１００ｍＬ中の水生菌を培養するため１００倍感度が高い検査方法といえる。Ｒ２Ａ寒天培地は、培地として寒天が用いられている。基本的には、寒天は透明であるため、寒天培地の背面に暗い色の物体を置いて撮影した画像を使用することで、適切な閾値による２値化により、基本的な菌を検出可能である。これに対して、メンブレンフィルタでは、培地の色が基本的には白色であり、場合によってさらに肉眼でのカウントのために格子が入っていることもある。
【０００７】
本発明は、精度よくコロニー数を計数するコロニー数計数装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
開示の技術は、上記課題を解決するために、以下の手段を採用する。
【０００９】
即ち、第１の態様は、
培養されたコロニーを含む培地領域及びフィルタに記載される格子を含む画像から前記培地領域の画像を抽出する領域抽出部と、
前記培地領域の画像を所定の大きさに調整する補正部と、
前記培地領域内の前記格子を検出する格子検出部と、
前記格子検出部で検出された前記格子を除く前記培地領域の画像の色を所定の基準の色に調整する色調整部と、
前記色調整部で調整された前記培地領域の画像の各画素の色について、赤色成分と青色成分との差分値をもとめ、差分値が周囲の画素の差分値よりも大きい領域を抽出する色検出部と、
前記色調整部で調整された前記培地領域内の画像について、円形の領域を抽出する円形抽出部と、
前記色検出部による抽出結果と前記円形抽出部による抽出結果とを統合する統合部と、
前記統合部によって統合された抽出結果に基づいて、抽出結果に含まれる領域の数をコロニーとして計数する計数部と、
を備えるコロニー数計数装置である。
【００１０】
開示の態様は、プログラムが情報処理装置によって実行されることによって実現されてもよい。即ち、開示の構成は、上記した態様における各手段が実行する処理を、情報処理装置に対して実行させるためのプログラム、或いは当該プログラムを記録したコンピュータ読み取り可能な記録媒体として特定することができる。また、開示の構成は、上記した各手段が実行する処理を情報処理装置が実行する方法をもって特定されてもよい。開示の構成は、上記した各手段が実行する処理を行う情報処理装置を含むシステムとして特定されてもよい。
【発明の効果】
【００１１】
開示の技術によれば、精度よくコロニー数を計数するコロニー数計数装置を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】図１は、実施形態のコロニー計数装置の機能構成例を示す図である
【図２】図２は、コロニー数計数装置を実現するコンピュータの例を示す図である。
【図３】図３は、コロニー数計数装置の動作フローの例を示す図である。
【図４】図４は、取得したコロニー数のカウント対象の画像の例を示す図である
【図５】図５は、２値化後の画像の例を示す図である。
【図６】図６は、培地検出部による処理後の画像の例を示す図である。
【図７】図７は、カボールフィルタの例を示す図である。
【図８】図８は、画像における線状カボールフィルタの適用結果の例を示す図である。
【図９】図９は、図８において、閾値３以上として２値化した結果を示す図である。
【図１０】図１０は、差分画像及び２値化画像の例を示す図である。
【図１１】図１１は、適応的閾値処理を説明する図である。
【図１２】図１２は、テンプレートマッチングにおけるテンプレートの例を示す図である。
【図１３】図１３は、誤検出除去部が検出した格子点の例を示す図である。
【図１４】図１４は、目視によるカウント数と実施形態の装置による自動カウント数との関係を示す図である。
【発明を実施するための形態】
【００１３】
以下、図面を参照して実施形態について説明する。実施形態の構成は例示であり、開示の構成は、開示の実施形態の具体的構成に限定されない。開示の構成の実施にあたって、実施形態に応じた具体的構成が適宜採用されてもよい。
【００１４】
〔実施形態〕
（構成例）
本実施形態のコロニー数計数装置は、菌が培養されたメンブレンフィルタが撮影された画像を取得し、画像内のメンブレンフィルタに存在する生菌数（コロニー数）をカウントする。メンブレンフィルタが培養液に浸されることによって、メンブレンフィルタ上で菌が培養される。ここでは、菌が培養されたメンブレンフィルタを、培地領域ともいう。
【００１５】
図１は、本実施形態のコロニー計数装置の機能構成例を示す図である。コロニー計数装置１００は、画像取得部１０２、培地検出部１０４、変形補正部１０６、格子検出部１０８、ホワイトバランス調整部１１０、コロニー検出部１１２、格子除去部１１４、コロニー数カウント部１１６を備える。
【００１６】
画像取得部１０２は、コロニー数のカウント対象の画像を取得する。コロニー数のカウント対象の画像は、臨床現場でデジタルスチルカメラやスマートフォンのカメラ等によって、菌が培養されたメンブレンフィルタを含む領域を撮影した画像である。
【００１７】
培地検出部１０４は、画像取得部１０２によって取得された画像から、培地領域を検出する。
【００１８】
変形補正部１０６は、培地検出部１０４で検出された培地領域の大きさを所定の大きさに調整する。
【００１９】
格子検出部１０８は、画像の培地領域に含まれる格子を検出する。
【００２０】
ホワイトバランス調整部１１０は、培地領域の画像のホワイトバランスを調整する。
【００２１】
コロニー検出部１１２は、ホワイトバランス調整部１１０で処理された画像について、暖色系のコロニーや白色系のコロニーを検出し、検出結果を統合して、出力する。
【００２２】
誤検出除去部１１４は、コロニー検出部１１２によって検出されたコロニーから、格子点、格子線、縁等によるコロニーの誤検出を除去する。
【００２３】
コロニー数カウント部１１６は、誤検出除去部１１４によって処理された画像について、コロニー数をカウントする。
【００２４】
画像ＤＢ１２０は、画像取得部１０２で取得されるコロニー数のカウント対象の画像などを格納する。画像ＤＢ１２０には、培地検出部１０４などで処理された画像が格納され得る。画像ＤＢ１２０には、格子検出部１０８などで使用されるフィルタなどが格納され得る。画像ＤＢ１２０には、コロニー数計数装置１００で使用される各種データが格納され得る。
【００２５】
コロニー数計数装置１００は、パーソナルコンピュータ（ＰＣ、Personal Computer）
のような汎用のコンピュータまたはサーバマシンのような専用のコンピュータを使用して実現可能である。また、コロニー数計数装置１００は、ワークステーション（ＷＳ、Work
Station）、ＰＤＡ（Personal Digital Assistant）のような専用または汎用のコンピュータ、あるいは、コンピュータを搭載した電子機器を使用して実現可能である。また、コロニー数計数装置１００は、スマートフォン、携帯電話、カーナビゲーション装置のような専用または汎用のコンピュータ、あるいは、コンピュータを搭載した電子機器を使用して実現可能である。
【００２６】
図２は、コロニー数計数装置を実現するコンピュータの例を示す図である。図２のコンピュータ１０００は、プロセッサ１００２、主記憶装置１００４、及び、二次記憶装置１００６や、通信インタフェース装置のような周辺装置２０００とのインタフェース装置１００８を含む。主記憶装置１００４及び二次記憶装置１００６は、コンピュータ読み取り可能な記録媒体である。
【００２７】
コンピュータ１０００は、プロセッサが記録媒体に記憶されたプログラムを主記憶装置１００４の作業領域にロードして実行し、プログラムの実行を通じて周辺機器が制御されることによって、所定の目的に合致した機能を実現することができる。
【００２８】
プロセッサ１００２は、例えば、ＣＰＵ（Central Processing Unit）やＤＳＰ（Digital Signal Processor）である。主記憶装置１００４は、例えば、ＲＡＭ（Random Access
Memory）やＲＯＭ（Read Only Memory）を含む。
【００２９】
二次記憶装置１００６は、例えば、ＥＰＲＯＭ（Erasable Programmable ROM）、ハー
ドディスクドライブ（ＨＤＤ、Hard Disk Drive）である。また、二次記憶装置１００６
は、リムーバブルメディア、即ち可搬記録媒体を含むことができる。リムーバブルメディアは、例えば、ＵＳＢ（Universal Serial Bus）メモリ、あるいは、ＣＤ（Compact Disc）やＤＶＤ（Digital Versatile Disc）のようなディスク記録媒体である。
【００３０】
インタフェース装置１００８は、例えば、ＬＡＮ（Local Area Network）インタフェースボードや、無線通信のための無線通信回路である。
【００３１】
周辺装置２０００は、上記の二次記憶装置や通信インタフェース装置の他、キーボードやポインティングデバイスのような入力装置や、ディスプレイ装置やプリンタのような出力装置を含む。また、入力装置は、カメラのような映像や画像の入力装置や、マイクロフォンのような音声の入力装置を含むことができる。また、出力装置は、スピーカのような音声の出力装置を含むことができる。
【００３２】
コロニー数計数装置１００を実現するコンピュータは、プロセッサが二次記憶装置に記憶されているプログラムを主記憶装置にロードして実行することによって、画像取得部１０２、培地検出部１０４、変形補正部１０６、格子検出部１０８、ホワイトバランス調整部１１０、コロニー検出部１１２、誤検出除去部１１４、コロニー数カウント部１１６としての機能を実現する。一方、画像ＤＢ１２０は、主記憶装置または二次記憶装置の記憶領域により実現される。
【００３３】
コロニー数計数装置１００の各ユニットは、ハードウェアの構成要素、ソフトウェアの構成要素、又は、これらの組み合わせとして、それぞれ実現され得る。
【００３４】
ハードウェアの構成要素は、ハードウェア回路であり、例えば、ＦＰＧＡ（Field Programmable Gate Array）、特定用途向け集積回路（ＡＳＩＣ）、ゲートアレイ、論理ゲー
トの組み合わせ、アナログ回路等がある。
【００３５】
ソフトウェアの構成要素は、ソフトウェアとして所定の処理を実現する部品である。ソフトウェアの構成要素は、ソフトウェアを実現する言語、開発環境等を限定する概念ではない。
【００３６】
一連の処理は、ハードウェアにより実行させることもできるが、ソフトウェアにより実行させることもできる。
【００３７】
プログラムを記述するステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくても、並列的または個別に実行される処理を含む。
【００３８】
（動作例）
コロニー数計数装置１００の動作例について、説明する。図３は、コロニー数計数装置１００の動作フローの例を示す図である。
【００３９】
ステップＳ１０１では、画像取得部１０２は、コロニー数のカウント対象の画像を取得する。画像取得部１０２は、デジタルスチルカメラやスマートフォンのカメラ等や、記憶装置等から、コロニー数のカウント対象の画像を取得する。
【００４０】
図４は、取得したコロニー数のカウント対象の画像の例を示す図である。図４の画像の中央には、円形のメンブレンフィルタが撮影されている。メンブレンフィルタは、格子を含む。
【００４１】
ステップＳ１０２では、培地検出部１０４は、画像取得部１０２が取得した画像（処理対象の画像）から、培地領域（メンブレンフィルタの領域）を検出する。市販されている多くのメンブレンフィルタは、紙状の円形の白色フィルタである。ここでは、最初に大津の２値化法を用いて培地領域を検出する。大津の２値化法では、培地検出部１０４は、画像をグレースケール画像に変換し、区間［０，２５５］のグレースケール画像のヒストグラムｈ（ｇ）｜ｇ＝０，・・・，２５５を算出する。ｈ（ｇ）は、濃淡値ｇである画素の頻度を表す。ここで、ｈ（ｇ）を正規化したものをｐ（ｇ）＝ｈ（ｇ）／ωとする。ωは、正規化の係数である。ｐ（ｇ）が双峰性を有する場合、２つの峰の中間を選択すれば良好な閾値を得られる。このとき、クラス内分散は、次式で表される。
【００４２】
【数１】
2015154718.tif 000003
【００４３】
また、クラス間分散は、次式で表される。
【００４４】
【数２】
2015154718.tif 000004
【００４５】
ここで、ｍはｐ（ｇ）の平均、σ２はｐ（ｇ）の分散である。ｍ１はある閾値でヒスト
グラムを分割した際の値の小さい方の平均値、ｍ２は値の大きい方の平均値、σ１２は値の小さい方の分散、σ２２は大きいほうの分散を表す。このとき、次式で表される分離度Ｂなる閾値を用いて、グレースケール画像が２値化される。
【００４６】
【数３】
2015154718.tif 000005
【００４７】
図５は、２値化後の画像の例を示す図である。図５において、白色の領域が培地領域の候補とされる領域である。実際の培地領域は、円形であるが、図５のように、撮影時の環境などによって培地領域以外も培地領域の候補として検出される。そこで、培地検出部１０４は、２値化後の画像中で最も領域が大きく、円形度の高い領域を培地領域として検出する。さらに、培地検出部１０４は、元の処理対象の画像において、処理の簡便性を考慮して、培地領域として検出された領域以外を、白色に塗りつぶした画像を生成する。
【００４８】
図６は、培地検出部による処理後の画像の例を示す図である。図６の例では、培地領域として検出された領域は、画像取得部１０２が取得した画像と同一の画像であり、培地領域として検出された領域以外の領域は、白色となっている。
【００４９】
２値化の方法として、大津の２値化法に限らず、K-means法などの他の方法が使用され
てもよい。
【００５０】
白色フィルタの代わりに黒色フィルタが使用されることがある。黒色フィルタが使用される場合は、上記の白色と黒色とを入れ替えることで、同様に処理することができる。
【００５１】
ステップＳ１０３では、変形補正部１０６は、培地検出部１０４で処理された画像に対して、培地領域の大きさを補正する。
【００５２】
画像における培地領域は、カメラを培地に平行にして撮影しないと、歪む可能性がある。通常、臨床現場では培地領域がカメラに平行になるようして撮影される。よって、ここでは、培地領域は、カメラに平行に撮影されていることを前提とする。これに対して、培地領域（メンブレンフィルタ）からカメラまでの距離は、コントロールすることが難しく、同じ大きさの培地領域であっても撮影毎に画像における培地領域の大きさが変化しうる。そこで、変形補正部１０６は、検出された培地領域の画像における直径を、所定の大きさになるように、拡大処理または縮小処理を行う。ここでは、例えば、培地領域の直径が５００ｐｉｘｅｌになるようにする。変形補正部１０６は、例えば、培地検出部１０４で処理された画像における培地領域の直径が３００ｐｉｘｅｌである場合、当該画像を縦横方向に１．６７倍して、培地領域の直径が５００ｐｉｘｅｌになるようにする。培地領域の直径のサイズは、５００ｐｉｘｅｌに限定されるものではなく、装置の負荷等に応じて自由に選択される。
【００５３】
培地領域の直径は、例えば、培地領域の座標のＸ軸方向（横軸方向）の最小値と最大値との差（即ち、Ｘ軸方向の両端の座標の差）によって求められる。培地の直径は、２次のモーメントのような主軸を求める方法によって、算出され得る。
【００５４】
また、培地領域の直径は、モーメント特徴を用いて、主軸方向ｔａｎθを求め、主軸の長さを求めることで求められる。主軸方向は、次式の（ｐ＋ｑ）次のモーメントＭ（ｐ，ｑ）を用いて計算される。
【００５５】
【数４】
2015154718.tif 000006
【００５６】
このとき、ｆｉｊは図形内で１、図形外で０をとる２値画像である。ここでは、培地領域を１、培地領域以外を０とする。主軸方向ｔａｎθは次式で求められる。
【００５７】
【数５】
2015154718.tif 000007
【００５８】
このとき、図形の重心が円の中心（培地領域の中心）である。
【００５９】
ステップＳ１０４では、格子検出部１０８は、変形補正部１０６で処理された画像（格子検出部１０８における処理対象の画像）から、メンブレンフィルタの格子を検出する。メンブレンフィルタでは、格子が入ったフィルタが使用されることがある。格子は、色補正やコロニー検出に影響するため事前に下記のように複数のガボールフィルタ（ここでは、４９画素×４９画素）等を適用し、このうち最大強度となるものを用いて格子検出を行う。
【００６０】
図７は、カボールフィルタの例を示す図である。図７に示すように、格子検出部１０８は、水平方向からほぼ３０度の６方向の線状ガボールフィルタｇｐ（ｉ，ｊ）（ｐ＝１，２，…６）を用意し、処理対象の画像の各画素において、次式のように値を求める。図７のカボールフィルタにおいて、色が暗い部分の値が小さく、色が明るい部分の値が大きい。ここで使用する線状ガボールフィルタは、回転不変性を得るためＭａｘｉｍｕｍ  Ｒｅｓｐｏｎｓｅ８（Varma, M. and Zisserman, A. (2004). “Unifying statistical texture classification frameworks”, Image Vision Computer,22(14):1175-1183）の線状フィルタを用いる。フィルタスケールはσｘ＝４，σｙ＝１２である。
【００６１】
【数６】
2015154718.tif 000008
【００６２】
ここでＩ（ｘ，ｙ）は処理対象の画像の座標（ｘ，ｙ）の画素値（グレースケール）、ｇｐは６方向のガボールフィルタである（０，３０，６０，９０，１２０，１５０度）。この６方向のカボールフィルタを使用すれば、ほぼどの方向の線状構造であっても検出しうる。この式における値が大きい場合、当該画素がいずれかの方向の線状構造に含まれることを意味する。即ち、値が大きい画素は、格子を形成する画素であることを意味する。また、この式における値が小さい場合、当該画素が、線状構造に含まれないことを意味する。即ち、値が小さい画素は、格子を形成する画素でないことを意味する。ここでは、格子が暗い線状であることを前提として、暗い線状構造が抽出される。
【００６３】
図８は、画像における線状カボールフィルタの適用結果の例を示す図である。図８では
、白い画素が値が大きい画素である。図８のように、線状カボールフィルタの適用結果、線状構造の強い画素の値が大きくなる。
【００６４】
図９は、図８において、閾値３以上として２値化した結果を示す図である。２値化の式を次式に示す。
【００６５】
【数７】
2015154718.tif 000009
【００６６】
ＧＲ’（ｘ，y）は、２値化後の座標（ｘ，y）の画素値である。ここで、ＧＲ’（ｘ，y）＝２５５とされる座標（ｘ，y）の画素が、格子検出部１０８で格子として検出される画素である。ＧＲ’（ｘ，y）による画像は、格子であるか否かを示す画像である。格子
であるか否かを示す画像において、格子を示す画素の画素値は２５５であり、格子を示さないがその画素値は０である。画素値２５５を白、画素値０を黒とすると、格子であるか否かを示す画像は、白黒の２値化画像となる。
【００６７】
格子検出部１０８は、カボールフィルタの代わりに、グレースケールブラックハット処理などの他の処理によって、格子を検出してもよい。
【００６８】
ステップＳ１０５では、ホワイトバランス調整部１１０は、変形補正部１０６で処理された画像の培地領域に対し、ホワイトバランスの調整を行う。血液透析の透析液の培養で検出対象となるのは水系桿菌である。水系桿菌の色調は、主に、白色、乳白色、淡黄色、レモン色、淡茶色、赤色、ピンク色の７色である。このため、白色系（白色、乳白色）以外の菌は暖色系の菌であり、暖色系では色情報を用いた検出が有効である。画像において、色情報を扱う際には、画像のホワイトバランスの調整を行うことが好ましい。撮影された環境（光や機材など）が異なると、元の色が同一の色でも画像上では異なる色となることがあるからである。メンブレンフィルタでは、培地の色が白である。ホワイトバランス調整部１１０は、培地の白を基準にして、ホワイトバランスの調整を行う。メンブレンフィルタの色は、一定であるとみなされるからである。具体的には、ホワイトバランス調整部１１０は、変形補正部１０６で処理された画像の培地領域から、格子検出部１０８で検出された格子の画素を除いて、培地領域とする。ここでは、ホワイトバランス調整部１１０は、培地領域の平均輝度値を求め、平均輝度値が基準値（Ｒ＝１８６，Ｇ＝１８３，Ｂ＝１６３）となるよう補正する。基準値はこれに限定されるものではない。基準値として、例えば、所定の光源のもとで撮影されたフィルタの画像における、フィルタ部分（培地領域）のＲＧＢ値（画素値）が採用されてもよい。ホワイトバランス調整部１１０は、培地領域のすべての画素について、次式に基づいて、ホワイトバランスの調整を行う。
【００６９】
【数８】
2015154718.tif 000010
【００７０】
ここで、ｏｍｇ’Ｒは、処理対象の画像の画素の処理後の赤色成分を表す。同様に、ｏｍｇ’Ｇは緑色成分、ｏｍｇ’Ｂは青色成分を表す。また、ＯＲｎは各色成分の重みを表し、次式で求められる。
【００７１】
【数９】
2015154718.tif 000011
【００７２】
ここで、ＷＭｎは、各色成分における格子を除いた培地の平均輝度である。また、右辺の分子の値は、各色成分の基準値を示す。
【００７３】
ホワイトバランスを調整することで、様々な環境で撮影されたメンブレンフィルタにおいて、色情報に基づいて、より正確にコロニーを検出することができる。ホワイトバランスを調整された培地領域の画像は、例えば、画像ＤＢ１２０に格納される。
【００７４】
ステップＳ１０６では、コロニー検出部１１２は、ホワイトバランス調整部１１０で、ホワイトバランスが調整された培地領域の画像（コロニー検出部１１２での処理対象の画像）を用いて、暖色系のコロニー及び白色系のコロニーの検出を行う。
【００７５】
ここで、使用される培地領域の画像では、格子検出部１０８で検出された格子の位置に相当する座標の画素値は、当該座標の周辺の画素に基づいて変更されてもよい。例えば、格子の位置に相当する座標の画素値は、当該座標の周辺の画素（例えば、５画素×５画素）の中央値（例えば、５画素×５画素では画素値が小さい方から１３番目）の画素値に置き換えられてもよい。ここで、周辺の画素から格子の画素は除かれてもよい。
【００７６】
〈暖色系〉
水系桿菌の色調は、白色系、黄色系、赤系、ピンク系に限定される。白色系以外は暖色系に属する。暖色系のコロニーでは、赤色成分が青色成分よりも大きいという特徴がある。そこで、コロニー検出部１１２は、暖色系のコロニーを、処理対象の画像の画素の赤色成分と青色成分の差分（色差）を取ることで検出する。処理対象の画像の各画素の赤色成分と青色成分の差分画像（差分による画像）は、次式で定義される。差分画像の各画素は、１次元の情報を有する。
【００７７】
【数１０】
2015154718.tif 000012
【００７８】
このとき、コロニー検出部１１２は、培地領域のノイズの影響を軽減するため平滑化処理を行ってもよい。平滑化処理は、例えば、処理対象の画素とその画素の周囲の画素の値の平均値を、処理対象の画素の値とする処理である。
【００７９】
この式において、右辺の絶対値をとることで、暖色系および寒色系のコロニーを検出することができる。寒色系のコロニーでは、青色成分が赤色成分よりも大きい。
【００８０】
図１０（１）は、差分画像の例を示す図である。図１０（１）の例では、差分画像の一
部分を示している。また、図１０（１）の例では、明るい部分では差分値が大きく、暗い部分では差分値が小さい。図１０（１）の例では、数箇所の差分値が大きい領域が見られる。
【００８１】
次に、コロニー検出部１１２は、差分画像を２値化することで、コロニーを検出する。ここで、例えば、２値化後の値１の領域がコロニーの領域で、値０の領域がコロニー以外の領域とする。２値化の際、固定の閾値により２値化すると、うまく検出できないことがある。そこで、ここでは、コロニー検出部１１２は、適応的閾値処理を用いて、多値閾値により、差分画像を２値化する。また、適応的閾値処理の代わりに、リング状フィルタ等のピークを検出する処理が使用されてもよい。
【００８２】
図１１は、適応的閾値処理の説明をする図である。図１１のグラフでは、横軸にＸ座標、縦軸に差分画像の色差をとる。ここでは、簡単のために、１次元の領域を用いて説明するが、差分画像のような２次元の領域に対しても同様である。図１１の曲線は、座標毎の色差を示している。ここでは、色差が周囲より大きい領域がコロニーとされる。図１１のグラフでは、３つの山が見られる。図１１のグラフの３つの山は、３つのコロニー（領域）と考えられる。ここで、固定の閾値を使用すると、閾値ａでは、領域の数は１つであり、閾値ｂとすると、領域の数は２つになり、閾値ｃとすると領域の数は３つとなり、閾値ｄとすると、領域の数は１つとなる。従って、固定の閾値では、どの閾値を採用しても領域の数は３つにならない。
【００８３】
そこで、適応的閾値処理を用いる。ある閾値以上の１つの領域において、このある閾値の一段階上の閾値以上の領域が、複数の領域になった場合、当該１つの領域において閾値を一段階上の閾値とする。これを閾値の範囲で繰り返す。
【００８４】
図１１を用いて、より具体的に説明する。ここでは、閾値をａ、ｂ、ｃ、ｄ（ａ＜ｂ＜ｃ＜ｄとする）の順に変化させるとする。最初に領域全体で、閾値を閾値ａとする。このとき、閾値ａ以上の領域は、領域α１の１つである。閾値ａから閾値ｂにしたとき、領域α１において、閾値ｂ以上の領域は、領域β１と領域β２の２つになる。このとき、領域α１において、閾値をｂとする。さらに、領域β１において、閾値ｂから閾値ｃに変化したとき、閾値ｃ以上の領域は、領域γ１と領域γ２の２つになる。ここで、領域β１の閾値をｃにする。このようにして、閾値の下限から上限まで閾値を変化させ、領域ごとに閾値を決めていく。このように閾値を決めると、領域α１を除く全領域の閾値は閾値ａ、領域β１を除く領域α１の閾値は閾値ｂ、領域β１の閾値は閾値ｃとなる。従って、閾値以上の領域は、領域γ１、領域γ２、領域β２の３つの領域となる。
【００８５】
コロニー検出部１１２は、適応的閾値処理を用いて、このように閾値を決定して、差分画像を２値化する。
【００８６】
図１０（２）は、図１０（１）において、閾値を２０としたときの２値化画像である。図１０（３）は、図１０（１）において、閾値を５としたときの２値化画像である。図１０（２）、図１０（３）では見られる領域の数が異なる。即ち、固定の閾値を用いると、複数の領域が１つの領域に見えたり、周囲より値が大きいが絶対的な値が小さい領域が見えなくなったりすることで、領域数を正確にカウントできない恐れがある。図１０（４）は、図１０（１）において、適応的閾値処理を用いて、閾値を決定したときに２値化画像である。図１０（４）では、大きな領域から、小さな領域まで、見えている。ここでは、閾値を、５から４０まで１ずつ変化させた。閾値の範囲は、目視によるカウント結果などに基づいて、暖色系コロニーと培地とを分けるように、決定される。
【００８７】
このようにして、コロニー検出部１１２は、暖色系コロニーを示す画像（２値化画像）
を生成する。暖色系コロニーを示す画像は、例えば、暖色系コロニーを示す画素の画素値が２５５で、暖色系コロニーを示さない画素の画素値が０である画像である。
【００８８】
〈白色系〉
次に白色系のコロニーの検出について説明する。白色系のコロニーは、赤色成分と青色成分との差分が小さいため、暖色系のコロニーと同様に検出することは難しい。そこで、コロニー検出部１１２は、処理対象の画像の画素値のうち緑色成分を用いて、テンプレートマッチング等により、円形で小さなコロニーの強調検出を行うことで、白色系のコロニーを検出する。円形領域の強調についてはマッチトフィルタなどの円形構造を強調する空間フィルタリング処理でも可能である。緑色成分の代わりに、赤色、緑色、青色の各成分から算出されるグレースケールの成分が使用されてもよい。
【００８９】
テンプレートマッチングは、テンプレートと呼ばれる画像領域と同じパターンが画像全体の中に存在するか否かを調べる処理である。テンプレートマッチングでは、テンプレートの画像を画像の中で移動させながら、テンプレートとテンプレートに重なる部分との画像の相関係数を算出する。相関係数は、画像の画素毎に算出され得る。テンプレートとテンプレートに重なる部分との画像の相関係数は、テンプレートと重なる部分の中心の画素の座標における値となる。相関係数が大きいほど、テンプレートとテンプレートに重なる部分の画像が類似している。相関係数が１のとき、テンプレートとテンプレートに重なる部分の画像が一致している。
【００９０】
図１２は、テンプレートマッチングにおけるテンプレートの例を示す図である。図１２のテンプレートは、円形のテンプレートである。
【００９１】
コロニー検出部１１２は、図１２のテンプレートと処理対象の画像の緑成分を用いて、テンプレートマッチングを行う。コロニー検出部１１２は、テンプレートマッチングにおいて、正規化相互相関を用いる。テンプレートマッチングにおける相関係数の算出式は、次式である。
【００９２】
【数１１】
2015154718.tif 000013
【００９３】
ここで、ｆは処理対象の画像の画素値（緑色成分）、ｔはテンプレートの画素値、ｔのバーはテンプレートの平均値、ｆｕ，ｖのバーはｆ（ｕ，ｖ）の平均値である。
【００９４】
これにより、コロニー検出部１１２によって、各画素のテンプレートとの相関係数を示す画像が生成される。相関係数を示す画像の各画素は、１次元の情報（類似度）を有する。ここで、コロニー検出部１１２は、各画素の相関係数を示す画像に対して、適応的閾値処理を適用し、コロニーの検出を行う。閾値は、相関係数０．０５から０．２までで、０．０１間隔とする。閾値の範囲は、目視によるカウント結果などに基づいて、白色系コロニーと培地とを分けるように、決定され得る。適応的閾値処理を利用することで、周囲より相関係数が高い領域を検出することができる。
【００９５】
これにより、コロニー検出部１１２は、白色系コロニーを示す画像（２値化画像）を生成する。白色系コロニーを示す画像は、例えば、白色系コロニーを示す画素の画素値が２５５で、白色系コロニーを示さない画素の画素値が０である画像である。
【００９６】
〈統合〉
コロニー検出部１１２は、差分画像を用いて抽出したコロニー（暖色系コロニー）と、テンプレートマッチングを用いて抽出したコロニー（白色系コロニー）との、抽出結果を統合する。コロニー検出部１１２は、暖色系コロニーを示す画像及び白色系コロニーを示す画像の両方においてコロニーとされない領域をコロニーでない領域（画素）とし、少なくともいずれか一方の画像でコロニーとされた領域（画素）をコロニーである領域（画素）とする。コロニー検出部１１２は、暖色系コロニーの画像と白色系コロニーの画像とから、コロニーの検出結果として、コロニーを示す画像（２値化画像）を生成する。
【００９７】
テンプレートマッチングを使用したコロニーの抽出は、感度が高いため、仮にコロニーのない画像に適用しても、コロニーの誤検出が発生する。そこで、コロニー検出部１１２は、差分画像を用いて抽出したコロニーの数が５個以上である場合に、白色系コロニーの抽出結果と統合してもよい。
【００９８】
コロニー検出部１１２は、暖色系コロニーを検出する部分について色検出部として、白色系コロニーを検出する部分について円形抽出部として、統合する部分について統合部として動作してもよい。
【００９９】
ステップＳ１０７では、誤検出除去部１１４は、コロニー検出部１１２で検出されたコロニーから、格子や縁によるコロニーの誤検出を除去する。
【０１００】
白色系のコロニーを検出する際、円形領域が強調されるため、格子の格子点も強調される傾向がある。コロニーが存在しない格子の格子点は、コロニーではないため、コロニーとして検出されると誤検出となる。
【０１０１】
そこで、誤検出除去部１１４は、格子検出部１０８が検出した格子の画像を細線化する。細線化は、画像を幅１ピクセルの線画像に変換する処理である。誤検出除去部１１４は、細線化した画像からすべての分岐点を、格子点として検出する。検出した格子点は、格子点の画像として画像ＤＢ１２０などに格納されてもよいし、格子点の座標として画像ＤＢ１２０などに格納されてもよい。誤検出除去部１１４は、コロニー検出部１１２が抽出した画像において、コロニーとして抽出した領域の内、格子点と抽出した位置から５画素以内の距離にあり、面積５０画素以下の領域を除去する。面積５０画素以下の領域は、コロニーとしては小さい。このような条件に合う領域は、格子の格子点である可能性が高い。
【０１０２】
図１３は、誤検出除去部１１４が検出した格子点の例を示す図である。図１３の例では、白い点が格子点に相当する。図１３の例では、メンブレンフィルタにおける複数の格子点が検出されている。
【０１０３】
さらに、誤検出除去部１１４は、コロニーとして抽出した領域の内、細線化した格子線から３画素以内で、面積が３０画素以下の領域を除去してもよい。
【０１０４】
臨床では、一度に、複数のメンブレンフィルタが培養される。そのため、メンブレンフィルタの外周にペンなどで識別の記号を記述することがある。また、メンブレンフィルタを撮影する際の光源は、一定ではないためメンブレンフィルタフィルタの外周に影が入り込み、メンブレンフィルタの輝度分布にばらつきが生じることがある。よって、誤検出除去部１１４は、培地領域の外周から３０画素以内の領域を除去する。培地領域の外周付近の領域は、ノイズが多いからである。
【０１０５】
ステップＳ１０８では、コロニー数カウント部１１６は、誤検出除去部１１４が誤検出を除去した画像（処理対象の画像）に対して、コロニー数をカウントする。コロニー数カウント部１１６は、処理対象の画像に対して、８近傍ラベリング処理を施す。
【０１０６】
ラベリング処理は、白黒に２値化された画像において、白の部分が連続した画素に同一番号（同一ラベル）を割り当てる処理である。８近傍ラベリング処理は、ラベリング処理において、白の部分が縦、横、斜めのいずれかに連続している場合に、同一ラベルを割り当てる処理である。同一のラベルを割り当てられた領域が１つの領域となる。従って、処理対象の画像に対するラベリング処理後のラベルの種類の数が、コロニー数を表す。これによって、コロニー数カウント部１１６は、コロニー数をカウントすることができる。
【０１０７】
（その他）
以上の処理で得られた結果で理想的な場合、良好なコロニー数のカウント結果が得られる。しかし、培地が濡れたことによる歪みや影の影響などで細かなノイズが含まれるケースが存在する。このため、誤検出除去部１１４は、面積が１０００画素以上、面積が４画素以下、離心率が０．９以上の小領域を除去してもよい。ここで、面積は各小領域の画素数である。離心率ｅは次式で求められる。
【０１０８】
【数１２】
2015154718.tif 000014
【０１０９】
ここで、２ａは小領域の長径、２ｂは小領域の短径である。ｅは真円に近いほど値が０に近づく。
【０１１０】
（実施形態の作用、効果）
コロニー数計数装置１００は、水性桿菌が培養されたメンブレンフィルタを撮影した画像から培地領域（メンブレンフィルタの領域）を抽出する。コロニー数計数装置１００は、抽出された培地領域の大きさを所定の大きさになるように補正する。コロニー数計数装置１００は、培地領域に含まれる格子を検出する。コロニー数計数装置１００は、培地領域の色が所定の色になるように、ホワイトバランスを調整する。コロニー数計数装置１００は、赤色成分と青色成分との差分値が周囲より大きい領域を抽出することで、暖色系のコロニーを抽出する。コロニー数計数装置１００は、円形の領域を抽出することで、白色系のコロニーを抽出する。コロニー数計数装置１００は、抽出された暖色系コロニーと白色系コロニーとを統合して、コロニーの数を計数することで、メンブレンフィルタ内のコロニー数を計数する。これにより、コロニー数計数装置１００は、水性桿菌が培養され格子が含まれるメンブレンフィルタの画像から、暖色系コロニー及び白色系コロニーの数を計数することができる。
【０１１１】
コロニー数計数装置１００は、暖色系のコロニーを検出する際に色情報を使用し、白色系のコロニーを検出する際に、コロニーの形状の情報を使用する。形状の情報を用いてコロニーを検出することで、画像内の培地領域の明るさの違いに影響されなくなる。また、コロニー数計数装置１００は、格子を検出する際、形状の情報を使用しているので、画像の撮影時の環境に影響されずに格子を検出できる。
【０１１２】
ここでは、メンブレンフィルタを撮影した画像について説明したが、寒天培地を撮影した画像についても同様に適用することができる。例えば、寒天培地を撮影する際には、背景に白色の用紙などを置くことにより同様に処理されうる。
【０１１３】
図１４は、目視によるカウント数と本実施形態の装置による自動カウント数との関係を示す図である。図１４のグラフは、複数の水性桿菌が培養されたメンブレンフィルタにおいて、目視によるカウントと、コロニー数計数装置１００による自動カウントの関係を示すグラフである。ここでは、９４枚のメンブレンフィルタについて、計数を行った。目視によるカウント数と自動カウント数との間の相関係数は、０．９６（Ｐ＜０．００１）である。これは、強い正の相関があることを意味する。また、ここで使用した９４枚のメンブレンフィルタの画像中で、カウント数の平均誤差は０．６５個であった。
【符号の説明】
【０１１４】
１００      コロニー数計数装置
１０２        画像取得部
１０４        培地検出部
１０６        変形補正部
１０８        格子検出部
１１０        ホワイトバランス調整部
１１２        コロニー検出部
１１４        誤検出除去部
１１６        コロニー数カウント部
１２０        画像ＤＢ
【図１】
2015154718.tif 000015
【図２】
2015154718.tif 000016
【図３】
2015154718.tif 000017
【図４】
2015154718.tif 000018
【図５】
2015154718.tif 000019
【図６】
2015154718.tif 000020
【図７】
2015154718.tif 000021
【図８】
2015154718.tif 000022
【図９】
2015154718.tif 000023
【図１０】
2015154718.tif 000024
【図１１】
2015154718.tif 000025
【図１２】
2015154718.tif 000026
【図１３】
2015154718.tif 000027
【図１４】
2015154718.tif 000028