WO2014021175A1

WO2014021175A1 - 壊死細胞領域検出装置及びその方法、コンピュータにより処理可能な壊死細胞領域検出プログラムを記憶する記憶媒体

Info

Publication number: WO2014021175A1
Application number: PCT/JP2013/070105
Authority: WO
Inventors: 新垣　英哉
Original assignee: オリンパス株式会社
Priority date: 2012-07-31
Filing date: 2013-07-24
Publication date: 2014-02-06
Also published as: US9336429B2; US20150131889A1; JP5984559B2; JP2014029287A

Abstract

　壊死細胞領域検出装置は、生細胞の細胞画像を取得する画像取得手段と、前記細胞画像を局所的画質特性が均一になるように複数領域に分割する領域分割手段と、前記細胞画像を低周波成分の低域画像と高周波成分の高域画像とを含む画像に分解する帯域分割手段と、前記各領域毎に前記高域画像からテクスチャ特徴量を算出する特徴量算出手段と、前記各領域毎に前記低域画像から輝度平均値を算出する輝度算出手段と、前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞による領域か否かを判別する判別手段とを含む。

Description

壊死細胞領域検出装置及びその方法、コンピュータにより処理可能な壊死細胞領域検出プログラムを記憶する記憶媒体

　本発明は、例えば明視野顕微鏡等を用いた撮像により取得された細胞画像から壊死細胞により形成される領域を検出する壊死細胞領域検出装置及びその方法、コンピュータにより処理可能な壊死細胞領域検出プログラムを記憶する記憶媒体に関する。

　従来から医療・ライフサイエンス分野では、顕微鏡を通して撮影された細胞画像を用いた様々な細胞解析が行われている。例えば、ＥＳ細胞、ｉＰＳ細胞等、幹細胞の研究においては、細胞分化メカニズムの解明や、創薬開発等を目的として、時系列で撮影された複数の細胞画像から細胞の分化過程や形態的特徴変化を観察し、細胞毎の性質の違いを調べるといった作業が一般的に行われている。

　細胞画像の解析に関しては、従来目視で行われていた個々の細胞のスクリーニング等の煩雑な作業が画像認識等の画像処理技術を応用することで自動化することが可能になりつつある。かかる画像処理技術を応用すれば、細胞画像に含まれる個々の細胞の輪郭線を算出し、細胞の形態情報や個体数等、或いは細胞の動きやその移動距離等を検出することで、細胞活性度合いを計測することが出来る。

　細胞の活動は、細胞分裂による細胞の増殖のみならず、アポトーシス（自然死）及びネクローシス（壊死）と呼ばれる細胞死滅のメカニズムの働きによりバランスを保ちながら成り立っている。　
　アポトーシスでは、初めに細胞核で変化が起こり、細胞が縮小してアポトーシス小体を形成した後、免疫細胞等に貧食され痕跡を残さず消化される。一方、ネクローシスは、細胞全体が次第に膨張し、細胞質が変化した後、最終的に細胞膜が破裂する。その際、周辺に細胞内容物が残留し、炎症などの原因となる（細胞溶解）。

　ネクローシスにより細胞内容物が残留した領域（以降、壊死領域と称する）は、細胞認識における個々の細胞の抽出において生細胞領域との区別を困難とし、正確な細胞数の計測、又は当該細胞数の計測による細胞活性度合いの評価等に悪影響を与える可能性が高いため、予め壊死領域を正確に把握する必要がある。

　特許文献１は、薬剤スクリーニングを目的とし、顕微鏡により取得した細胞画像中の各細胞の形態的な特徴を表す形態特徴量に基づき生細胞又は死細胞のいずれかに判定する方法を開示する。具体的に特許文献１は、細胞の形態を表現する形態特徴量として細胞の輪郭形状の円又は球からのずれ量を用いた判定を行う。この判定の結果、ずれ量の大きな細胞は、円形又は球形を標準としてかなりの歪みを生じているものとし、活発に活動を行っている生細胞と判定する。これに対してずれ量が小さく円形又は球形に近い細胞は、活性を失い死滅に近づいたと判定する。

特開２００７－２０４４９号公報

　特許文献１は、死細胞を円形又は球形に近い形態を成すことを想定しているが、壊死により細胞膜が破壊され、細胞内容物が残留するケースでは、必ずしも円形状を成すとは限らない。

　特に位相差顕微鏡による細胞画像に対する処理において顕著であるが、細胞画像中から個々の細胞を抽出する際に、細胞輪郭線と細胞内部の構造とに基づくエッジ成分（例えば細胞内部に位置する細胞核周囲のエッジ等）とが混在して細胞の区別がつきにくい場合と、或いは複数の細胞が密集し、細胞輪郭線の不明瞭な部分が多く見られる場合とがある。これらの場合には、本来複数の細胞領域であるものを一つの領域として抽出する（未分割）、或いは１つの細胞領域を複数の領域として抽出してしまう（過分割）などの理由で、正確な細胞形態を検出することは難しい。　
　従って、ずれ量を算出するための前提となる円の当てはめも、必ずしも精度良く当てはめが行われるとは限らず、また一つ一つの細胞個別に当てはめられるとも限らない。

　本発明は、上記問題を鑑みてなされたものであり、壊死細胞領域を高精度に検出することが可能な壊死細胞領域検出装置及びその方法、コンピュータにより処理可能な壊死細胞領域検出プログラムを記憶する記憶媒体を提供することにある。

　本発明の主要な局面に係る壊死細胞領域検出装置は、経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得する画像取得部と、前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割する領域分割部と、前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解する帯域分割部と、前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出する特徴量算出部と、前記複数の領域毎に、前記低域画像から輝度平均値を算出する輝度算出部と、前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別する判別部と、を具備する。

　本発明の主要な局面に係る壊死細胞領域検出方法は、コンピュータの処理によって、経時的に変化する生細胞を複数の時点でそれぞれ撮像を行って取得した複数の細胞画像から成る細胞画像群を取得し、前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割し、前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解し、前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出し、前記複数の領域毎に、前記低域画像から輝度平均値を算出し、前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別する。

　本発明の主要な局面に係るコンピュータにより処理可能な記憶媒体は、経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得させる画像取得機能と、前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割させる領域分割機能と、前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解させる帯域分割機能と、前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出させる特徴量算出機能と、前記複数の領域毎に、前記低域画像から輝度平均値を算出させる輝度算出機能と、前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別させる判別機能と、を実現させる壊死細胞領域検出プログラムを記憶する。

　本発明によれば、壊死細胞領域を高精度に検出することが可能な壊死細胞領域検出装置及びその方法、コンピュータにより処理可能な壊死細胞領域検出プログラムを記憶する記憶媒体を提供できる。

図１は本発明に係る壊死細胞領域検出装置の第１の実施の形態を示す構成図である。同装置における帯域分割部を示す具体的な構成図である。同装置における撮像部により取得された細胞画像と帯域分割部により生成される低域画像及び高域画像の各輝度値分布を一次元的に表した一例とを示す図である。同装置における撮像部により取得された細胞画像を示す模式図である。同装置における領域分割部により領域分割するときの細胞の境界線部分を示す模式図である。同装置に適用する注目領域内の画素値の一例を示す図である。同装置により記録される隣接した画素対に関して出現頻度をカウントして算出される同時生起行列の一例を示す図である。同装置の判別部によるエントロピーと輝度平均値との空間で線形判別関数を用いての壊死領域・非壊死領域の判別を説明するための図である。同装置の壊死細胞領域検出装置フローチャートである。本発明に係る壊死細胞領域検出装置の第２の実施の形態を示す構成図である。同装置の壊死細胞領域検出装置フローチャートである。

[第１の実施の形態]
　以下、本発明の第１の実施の形態について図面を参照して説明する。　
　図１は壊死細胞領域検出装置の構成図を示す。本装置は、撮像部１００と、帯域分割部１１０と、領域分割部１２０と、テクスチャ特徴量算出部１３０と、輝度平均算出部１４０と、判別部１５０と、記録部１６０とを含む。　
　撮像部１００は、帯域分割部１１０と領域分割部１２０とに接続されている。帯域分割部１１０は、テクスチャ特徴量算出部１３０と、輝度平均算出部１４０とに接続されている。領域分割部１２０は、テクスチャ特徴量算出部１３０と、輝度平均算出部１４０と、判別部１５０とに接続されている。テクスチャ特徴量算出部１３０と輝度平均算出部１４０とは、それぞれ判別部１５０に接続されている。判別部１５０は、記録部１６０に接続されている。各部１００、１１０、…、１６０は、それぞれ例えばシステムコントローラに接続されて動作制御される。

　各部１００、１１０、…、１６０は、例えばＣＰＵ（中央演算処理装置）と演算プログラムを格納するＲＡＭ、ＲＯＭ等の記憶装置などから構成してもよい。ＲＯＭには、演算プログラムとしての壊死細胞検出プログラムが格納されている。　
　壊死領域判別装置プログラムは、コンピュータとしてのＣＰＵに、経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得させる画像取得機能と、細胞画像群中における所定時点に取得した細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割させる領域分割機能と、所定時点で取得した細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解させる帯域分割機能と、複数の領域毎に、高域画像からテクスチャ特徴量を算出させる特徴量算出機能と、複数の領域毎に、低域画像から輝度平均値を算出させる輝度算出機能と、複数の領域毎に対して少なくともテクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別させる判別機能と、を実現させる。

　壊死領域判別装置プログラムは、領域分割機能により分割された各領域毎に、所定時点に取得した細胞画像と、所定時点よりも以前又は以降の時点のうちいずれか一方又は両方に取得した細胞画像との間における局所的な輝度分布の類似度を算出させる類似度算出機能を含む。　
　判別機能は、各領域毎に対してテクスチャ特徴量と輝度平均値と類似度とに基づいて壊死細胞により形成されている領域か否かを判別させる。

　撮像部１００は、観察対象の細胞群を撮像して細胞画像を取得する。この撮像部１００は、経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得する画像取得部としての機能を有する。撮像部１００は、例えばＣＣＤ等の撮像素子と、Ａ／Ｄ変換器とを有する。撮像部１００は、例えば位相差顕微鏡（Phase contrast microscope)にカメラを取り付けたもので、このカメラにより位相差顕微鏡で拡大された細胞群の位相差像を撮像する。撮像部１００は、位相差顕微鏡にカメラを取り付けたものに限定されることはなく、例えば微分干渉顕微鏡（Differential interference contrast microscope：ＤＩＣ）等、他の明視野顕微鏡に対して適用することも可能である。

　撮像部１００は、位相差顕微鏡により撮影された細胞群の位相差像を撮像素子とＡ／Ｄ変換器等とを介してデジタル信号に変換し、例えば８ビット（２５６階調）のモノクロ原画像信号Ｆとして出力する。モノクロ原画像信号Ｆは、帯域分割部１１０と領域分割部１２０に転送される。

　位相差顕微鏡は、光の回折現象を利用し、異なる屈折率を持つ物質間を透過する光の位相差（光路差）をコントラストとして得ることができるので、透明な細胞や微生物等の対象物を観察するのに適している。位相差顕微鏡により撮影した画像は、背景領域と試料との境界線、細胞内部構造の細部ディテール周辺、或いは死滅した細胞の残骸等のごみ領域周辺においてハロ（アーティファクト）と呼ばれる強いコントラストを発生するという特徴を含む。ハロは、特に背景と個々の細胞との境界部分にオーラ状の光として強く表れる。　
　本実施の形態において位相差顕微鏡による位相差像は、見かけ上、背景領域が明るく、細胞領域が相対的に暗く撮影されるポジティブコントラスト像である。位相差顕微鏡による位相差像は、ポジティブコントラスト像に限定されるものではなく、ネガティブコントラスト像の場合においても、階調を反転させポジティブコントラスト像のように扱うことで処理することが可能である。

　帯域分割部１１０は、所定時点で取得した細胞画像から低周波成分により成る低域画像と、高周波成分により成る高域画像とを含む複数の帯域画像に分解する帯域分割部としての機能を有する。具体的に帯域分割部１１０は、所定の多重解像度分解処理を行うことによって、モノクロ原画像信号Ｆを異なる周波数帯域成分を含む複数の帯域画像に分解する。帯域分割部１１０は、モノクロ原画像信号Ｆの中の低周波成分を含む低域画像Ｌと、モノクロ原画像信号Ｆの中の高周波成分を多く含む高域画像Ｈとの２つの成分画像に分解する。

　低域画像Ｌは、細胞画像の中の背景領域上や細胞領域内部に存在する微細構造・ディテール、及びノイズを取り除き、かつ背景領域と細胞領域との輝度変化に差が現れやすい周波数帯域を多く含むことが好ましい。　
　高域画像Ｈは、細胞画像の中で細胞を構成するエッジ・ハロによる高周波成分をできるだけ多く含むことが好ましい。

　図２は帯域分割部１１０の具体的な構成図を示す。帯域分割部１１０は、例えば低域画像生成手段としての低域生成部２００と、高域画像生成手段としての高域生成部２１０とを含む。低域生成部２００と高域生成部２１０との各入力側は、それぞれ撮像部１００に接続されている。低域生成部２００は、高域生成部２１０と輝度平均算出部１４０とに接続されている。高域生成部２１０は、テクスチャ特徴量算出部１３０に接続されている。

　低域生成部２００は、細胞画像に対して平滑化処理を行って低域画像Ｌを生成する。具体的に低域生成部は、撮像部１００から転送されたモノクロ原画像信号Ｆに対して所定の平滑化フィルタ、例えばガウシアンフィルタを適用し、このフィルタ出力を低域画像Ｌとして輝度平均算出部１４０と高域生成部２１０とに転送する。本実施の形態において低域生成部２００は、ガウシアンフィルタを用いて平滑化しているが、これに限らず、低周波成分を抽出する手段であれば如何なるものでも適用可能である。

　図３は撮像部１００により取得された細胞画像の輝度値分布を一次元に簡略化して示した例である。低域生成部２００は、元の細胞画像に対して平滑化処理を行って低域画像Ｌを生成する。

　図４は撮像部１００により取得された細胞画像の模式図を示す。細胞画像には、各生細胞Ｑと、これら生細胞Ｑ間に存在する壊死細胞やその残留物Ｇとを含む。

　高域生成部２１０は、細胞画像から低域画像Ｌを減算して高域画像Ｈを生成する。具体的に高域生成部２１０は、撮像部１００から転送されたモノクロ原画像信号Ｆと、低域生成部２００から転送された低域画像Ｌとの間の対応する画素値間で各差分値を求め、これら差分値を高域画像Ｈとしてテクスチャ特徴量算出部１３０に転送する。図３は細胞画像から低域画像Ｌを減算して生成された高域画像Ｈの輝度値の分布の一例を示す。　
　このように帯域分割部１１０は、モノクロ原画像信号Ｆから低域画像Ｌと高域画像Ｈとを生成する。

　本実施の形態は、モノクロ原画像信号Ｆを低域画像Ｌと高域画像Ｈとの２つの帯域の画像に分けているが、これに限ることがない。本実施の形態は、モノクロ原画像信号Ｆに対して多重解像度分解を行うことにより帯域をより細かく分解し、この分解によりモノクロ原画像信号Ｆを３つ以上の帯域の画像にそれぞれ分割する。本実施の形態は、分割された各帯域画像の中から背景と細胞領域の輝度変化とが顕著に表れやすい帯域画像、及び細胞・死滅細胞等を構成するエッジ・ハロを多く含む帯域画像を所定の条件により選択する。所定の条件は、例えば画素値のコントラスト・分散に基づく闘値処理等である。本実施の形態は、それぞれ選択された低域画像Ｌと高域画像Ｈとして適用することが可能である。

　領域分割部１２０は、細胞画像群中における所定時点に取得した細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割する領域分割部としての機能を含む。領域分割部１２０は、例えば、細胞画像に対して領域分割処理を行い、細胞画像中の細胞毎に個々の領域として切り分ける。細胞画像に対する領域分割処理は、処理対象の細胞画像に対し、当該細胞画像内の特徴が類似し空間的に近接した１つ以上の画素集合により構成される領域に分割するための処理である。

　一般に位相差顕微鏡により撮影した細胞画像では、細胞境界線上の輝度が細胞内部の輝度と比較して相対的に高い。本実施の形態では、かかる細胞画像の特徴を考慮し、領域分割手法の一つであるウォーターシェッド法（分水嶺領域分割法）を用いた領域分割を行い、個々の細胞領域に切り分ける。ウォーターシェッド法は、細胞画像の輝度値勾配に基づく分割を行い、同画像中の輝度値が高く輝度値勾配の高い部分、すなわち図５に示すように細胞の境界線部分を分割線にして分割を行う。分割された個々の領域には、ラベリング処理により領域番号が付与され、当該領域番号を画素値とした領域分割画像Ｋを生成する。図５では３個の細胞領域に分割され、それぞれに対して領域番号１、２、３が付与されている。また、非細胞領域が背景領域となる領域に対しては領域番号０が付与されている。領域分割処理は、必ずしもウォーターシェッド法に限定されるものではなく、細胞領域を的確に分割できればどのような手法でも適用可能である。この生成された領域分割画像Ｋは、テクスチャ特徴量算出部１３０と、輝度平均算出部１４０と、判別部１５０とにそれぞれ転送される。

　テクスチャ特徴量算出部１３０は、領域分割画像Ｋ中の個々の各領域毎に、高域画像Ｈ上の画素値分布に基づくテクスチャ特徴量を算出する特徴量算出部としての機能を含む。テクスチャ特徴量算出部１３０は、高域画像Ｈ上の各画素値に基づき各領域毎に、局所的なテクスチャ特徴量を算出する。テクスチャ特徴量は、高域画像Ｈにおける画素値分布のランダム性に基づく特徴量、又は高域画像Ｈにおける画素値分布の複雑度に基づく特徴量、又は高域画像Ｈにおける画素の同時生起行列に基づく特徴量である。

　本実施の形態では、テクスチャ解析手法として広く知られる同時生起行列に基づくテクスチャ特徴量の一つであるエントロピーを適用する。以下、本実施の形態では、高域画像Ｈ中の所定の注目領域に関して同時生起行列を作成し、テクスチャ特徴量としてエントロピーを算出する方法について説明する。　
　同時生起行列は、画像中のテクスチャの統計的特徴量算出部として広く知られている方法の一つである。同時生起行列は、画像内に含まれる一定の位置関係にある画素対の出現頻度・確率を行列（同時生起行列）の形で表す。同時生起行列からは、様々なテクスチャ特徴量を算出することができる。

　同時生起行列のサイズを抑えて計算量を小さくするために、予め高域画像Ｈの階調数を所定の数に圧縮した階調圧縮画像が作成される。同時生起行列の大きさは、階調数×階調数の正方行列となる。本実施の形態では、例えば元の２５６階調（０～２５５）を、４階調（画素値＝０～３）に圧縮する。　
　はじめに、前述の階調圧縮画像内の注目画素を中心とした所定の大きさの注目領域が設定される。本実施の形態では、説明を簡略化するために、注目領域のサイズを例えば５×５画素領域とする。領域の形状は、本実施の形態のように矩形である必要は無く、領域分割部１２０における領域分割結果に基づく任意の形状をしていてもかまわない。

　次に注目領域から抽出する画素対の位置関係δが設定される。本実施の形態では、横方向に隣接する（画素間の距離ｄ＝１、角度θ＝０°）画素対を設定し、この画素対を成す左側の画素をｉ、右側の画素をｊとし、それぞれの画素値をＬｉ、Ｌｊとする。なお、ｉ＝０，１，２，３，…，ｎであり、ｊ＝０，１，２，３，…，ｍである。

　注目領域に含まれる全てに隣接した画素対に関し、各画素対に存在する出現頻度をカウントし、このカウント値を同時生起行列Ｐδ（Ｌｉ、Ｌｊ）に記録する。すなわち、注目領域内で、画素対Ｌｉ、Ｌｊが存在する頻度を同時生起行列ＰδのＬｉ行Ｌｊ列要素に記録する。　
　図６Ａは具体的な注目領域内の画素値の一例を示し、図６Ｂはその際に算出される同時生起行列の一例を示す。図６Ａに示す注目領域内ではＬｉ＝３、Ｌｊ＝２の画素対が２対存在しているので、図６Ｂに示す同時生起行列の要素Ｐδ（３，２）＝２となっている。図６Ａに示す注目領域内ではＬｉ＝０、Ｌｊ＝０の画素対が５対存在しているので、図６Ｂに示す同時生起行列の要素Pδ（０，０）＝５となっている。

　画素対Ｌｉ、Ｌｊが存在する全ての頻度を記録した後、出現頻度の総数でＰδ（Ｌｉ、Ｌｊ）の各要素を除算することで正規化する。そして、算出した同時生起行列Ｃに基づくテクスチャ特徴量を算出する。　
　本実施の形態では、テクスチャ特徴量として式（１）により定義されるテクスチャ特徴量であるエントロピー（Ｅntropy）を適用する。ここで、Ｌは、行列の大きさ（階調数）を表す。テクスチャ特徴量であるエントロピーは、画素値分布のでたらめさ、ランダム性を計る指標である。テクスチャ特徴量のエントロピーは、注目領域内に画素値がランダムに含まれるほど値が小さくなる。　

　高域画像Ｈの画素毎に算出したテクスチャ特徴量のエントロピーは、判別部１５０に転送される。　
　同時生起行列Ｃから算出できるテクスチャ特徴量は、様々に定義されている。例えば以下に示す角二次モーメント（Angular second moment）や、分散（Variance）の逆数もテクスチャ特徴量として適用可能である。　
　角二次モーメントは、式（２）に示すように定義される。角二次モーメントは、特定の画素対が多く存在し、一様性が大きいほど値が大きくなる。　

　分散は、式（３）に示すように定義される。分散は、注目領域に含まれる画素値差が大きく、要素のばらつき・複雑さが大きいほど値が大きくなる。分散の逆数は、逆に小さくなる。

　輝度平均算出部１４０は、領域分割画像Ｋ内の領域ごとに、当該各領域内に含まれる低域画像Ｌ上の画素値（輝度値）の輝度平均値を算出する。輝度平均算出部１４０は、低域画像Ｌに基づいて輝度の平均値を算出することで、細胞画像中に含まれるノイズや、輝度の極端な変動に影響されずに、領域分割画像Ｋ内の領域毎の輝度平均値を安定して算出することができる。

　判別部１５０は、領域のテクスチャ特徴量、及び平均輝度値に基づき壊死領域か否かの判別を行う。この判別は、２変数（すなわちテクスチャ特徴量・輝度平均値）に基づき２クラス、すなわち壊死領域と非壊死領域とに分類する判別分析と等価である。本実施形態では、線形判別関数を用いた判別分析に基づき壊死領域を特定する。　
　図７は判別部１５０によるエントロピーと輝度平均値で構成されるベクトル空間における線形判別関数Ｈを用いた壊死領域・非壊死領域判別を説明するための図を示す。判別部１５０は、エントロピーと輝度平均値との空間において線形判別関数Ｈを境に壊死領域と非壊死領域とを判別する。

　この判別分析は、予め与えられた学習用サンプルデータに基づき、複数のクラス（＝グループ）に分類するための判別関数（判別の基準となる関数）を得る。この判別分析では、新しくデータを得た際に、判別関数を用いてどのクラスに属するかを判別する。この判別分析では、領域分割画像Ｋ内の領域毎にテクスチャ特徴量、平均輝度により構成される２変数ベクトルデータを予め求めた線形判別関数に入力することで、領域分割画像Ｋ内の領域が壊死領域か否かを判定する。又、この判別分析では、判別分析手法として一般に広く知られるフイッシャーの線形判別関数を適用する。

　以下処理の詳細について説明する。　
　予め目視確認により領域の属するクラス（以降、壊死領域をクラスＣｌと称し、その他の非壊死領域をクラスＣ２と称する）が確定している学習用サンプルデータを用意する。　
　サンプルデータは、壊死領域、または非壊死領域に対する“テクスチャ特徴量”、“平均輝度”の２変数を要素とした特徴ベクトルＸの集合として表され、各特徴ベクトルが属するクラス（Ｃ１、またはＣ２）が目視により確認・決定されている。これら特徴ベクトルＸの集合を、２変数に基づく特徴空間上で、クラスＣ１、Ｃ２に最も良く分離することが可能な線形判別関数をフイッシャーの判別基準により予め定義しておく。線形判別関数は、式（４）により表される。

　判別対象の領域特徴ベクトルＸｓが入力された際は、式（５）に示す通りその出力値の正負に基づき壊死領域の判別を行う。　
　ｆ（ｘ）＝ｗ・ｘ＋ｂ　　　　　　　　　　　…（４）　
　ｗ・Ｘｓ＋ｂ＞０→Ｃ１　
　ｗ・Ｘｓ＋ｂ＜０→Ｃ２　　　　　　　　　　…（５）　
　ここで、ｗは重みベクトル、ｂはバイアス項である。ｆ（ｘ）＝０が線形判別関数となる。　
　フイッシャーの判別基準は、サンプルデータから、クラス内共分散Ｓｗ（式（６））と、クラス間共分散ＳＢ（式（７））との比に基づく目的関数（式（８））を定め、この目的関数を最大化、すなわちクラス内共分散を小さく、かつクラス間共分散を大きく）する重みベクトルｗを用いた判別関数がサンプルデータの２クラスに属する各データを互いに最も精度良く分離できるものと仮定するものである。

　ここで、式（８）の最大化は、式（９）に示す固有値問題に書き換えることができる。壊死領域の判別は、式（９）の固有値に対応する固有ベクトルに基づき重みベクトルｗを定めた線形判別関数を用いて行われる。判別結果は、記録部１６０に転送される。

　記録部１６０は、領域分割画像Ｋ内の各領域と、各領域に対応した壊死領域か否かの判定結果とを対としたデータを所定のメモリやファイル等の記録媒体に書き込む。

　次に、上記の如く構成された装置の動作について図８に示す壊死細胞領域検出装置フローチャートに従って説明する。　
　撮像部１００は、観察対象の細胞群を撮像し、細胞画像のモノクロ原画像信号Ｆを出力する。具体的に撮像部１００は、経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得し、そのモノクロ原画像信号Ｆを出力する。　
　帯域分割部１１０は、ステップＳ１０において、撮像部１００から出力されるモノクロ原画像信号Ｆ、すなわち細胞画像を入力する。

　帯域分割部１１０は、ステップＳ２０において、細胞画像から低周波成分により成る低域画像と、高周波成分により成る高域画像とを含む複数の帯域画像に分解する。すなわち、低域生成部２００は、撮像部１００から転送されたモノクロ原画像信号Ｆに対して所定の平滑化フィルタ、例えばガウシアンフィルタを適用し、このフィルタ出力を図３に示すような低域画像Ｌを生成して輝度平均算出部１４０と高域生成部２１０とに転送する。　
　これと共に高域生成部２１０は、撮像部１００から転送されたモノクロ原画像信号Ｆと、低域生成部２００から転送された低域画像Ｌとの間の対応する画素値間で各差分値を求め、これら差分値を図３に示すような高域画像Ｈとしてテクスチャ特徴量算出部１３０に転送する。

　領域分割部１２０は、ステップＳ３０において、細胞画像群中における所定時点に取得した細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割する。分割された個々の領域には、ラベリング処理により領域番号が付与され、当該領域番号を画素値とした領域分割画像Ｋを生成する。

　テクスチャ特徴量算出部１３０は、ステップＳ５０において、領域分割画像Ｋ中の個々の領域毎に、高域画像Ｈ上の画素値分布に基づくテクスチャ特徴量を算出する。テクスチャ特徴量の算出では、同時生起行列に基づくテクスチャ特徴量の一つであるエントロピーを適用する。テクスチャ特徴量の算出に際し、テクスチャ特徴量算出部１０２は、ステップＳ４０において、注目領域を設定し、この注目領域に対し、領域内に含まれる高域画像Ｈ上の画素値分布に基づくテクスチャ特徴量の一つであるエントロピーを算出する。

　輝度平均算出部１４０は、ステップＳ６０において、領域分割画像Ｋ内の領域毎に、低域画像Ｌ上の画素値（輝度値）の輝度平均値を算出する。輝度平均算出部１４０は、低域画像Ｌに基づいて平均輝度値を算出することで、細胞画像中に含まれるノイズや、輝度の極端な変動に影響されない安定した領域内平均輝度値を算出する。

　判別部１５０は、ステップＳ７０において、テクスチャ特徴量としてのエントロピーと輝度平均値とから、線形判別に基づいて領域毎に壊死細胞により形成されているか否かを判別する。判別部１５０は、例えば図７に示すようにエントロピー－平均輝度値の特徴空間において線形判別関数Ｈを境に壊死領域と非壊死領域とを判別する。

　記録部１６０は、領域分割画像Ｋ内の各領域と、当該各領域に対応した壊死領域か否かの判定結果とを対としたデータを所定のメモリやファイル等の記録媒体に書き込む。　
　ＣＰＵは、ステップＳ８０において、判別部１５０により細胞画像の全領域に対して判別処理を終了したか否かを判別する。細胞画像内の全領域に対して判別処理が終了すると、ＣＰＵは、処理を終了する。未処理の領域が存在する場合、ＣＰＵは、ステップＳ４０に移行する。

　このように上記第１の実施の形態によれば、撮像部１００により取得された細胞画像群中における所定時点に取得した細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割し、かつ細胞画像に対して高域画像Ｈと低域画像Ｌとに分解し、領域分割画像Ｋ内の複数の領域毎に高域画像Ｈからテクスチャ特徴量を算出し、領域分割画像Ｋ内の複数の領域毎に低域画像Ｌから輝度平均値を算出し、領域分割画像Ｋ内の複数の領域に対してテクスチャ特徴量と輝度平均値に基づく線形判別により壊死細胞により形成されている領域か否かを判別するので、細胞検出の際に妨げとなる、細胞壊死による残留物等で構成された壊死領域を高精度に検出することができる。

　高域画像Ｈに基づきテクスチャ特徴量を求めるので、細胞画像の全体に亘って明るさムラ等の影響を受けず、細胞画像の局所的な特徴をテクスチャ特徴量として求めることができる。低域画像Ｌに基づき領域分割画像Ｋ内の複数の領域毎に輝度平均値を算出するので、細かな輝度変動によらず安定して領域分割画像Ｋ内の領域毎に輝度平均値を算出できる。これらテクスチャ特徴量と輝度平均値とを用いることで壊死細胞により形成されている領域か否かを高精度に判別することができる。

　テクスチャ特徴量は、画素値分布のランダム性に基づく特徴量であるので、壊死細胞領域をその領域内の画素値分布のランダム性に基づき精度良く区別することができる。テクスチャ特徴量は、画素値分布の複雑度に基づく特徴量であるので、壊死細胞領域をその領域内の画素値分布の複雑度に基づき精度良く区別することができる。テクスチャ特徴量は、同時生起行列に基づく特徴量であるので、壊死細胞領域の特徴を効率よく的確に表現することができる。

[第２の実施の形態]　
　次に、本発明の第２の実施の形態について図面を参照して説明する。なお、図１と同一部分には同一符号を付してその詳しい説明は省略する。　
　図９は壊死細胞領域検出装置の構成図を示す。本装置は、上記第１の実施の形態と対比して、バッファ３００と、類似度算出部３１０とを付加し、かつテクスチャ特徴量算出部３２０と判別部３３０との各機能を変更している。

　撮像部１００は、バッファ３００を介して帯域分割部１１０と、領域分割部１２０と、類似度算出部３１０とにそれぞれ接続されている。帯域分割部１１０と領域分割部１２０とは、テクスチャ特徴量算出部３２０に接続されている。領域分割部１２０は、類似度算出部３１０に接続されている。テクスチャ特徴量算出部３２０と類似度算出部３１０と輝度平均算出部１４０と領域分割部１２０とは、それぞれ判別部３３０に接続されている。この判別部３３０は、記録部１６０に接続されている。

　これら各部１００、…、３００、３１０、３３０は、例えばＣＰＵ（中央演算処理装置）と演算プログラムを格納するＲＡＭ、ＲＯＭ等の記憶装置などから構成してもよい。ＲＯＭには、演算プログラムとしての上記壊死細胞検出プログラムが格納されている。壊死領域判別装置プログラムは、コンピュータとしてのＣＰＵに、上記画像取得機能と、上記領域分割機能と、上記帯域分割機能と、上記特徴量算出機能と、上記輝度算出機能と、上記判別機能とを実現し、かつ上記領域分割機能により分割された複数の領域毎に、所定時点に取得した細胞画像と、所定時点よりも以前又は以降の時点のうちいずれか一方又は両方に取得した細胞画像との間における局所的な輝度分布の類似度を算出させる類似度算出機能を含む。判別機能は、複数の領域に対してテクスチャ特徴量と輝度平均値と類似度とに基づいて壊死細胞により形成されている領域か否かを判別させる。

　本第２の実施の形態は、細胞の経時的変化を観察するために複数の時点において撮像した細胞画像を対象としている。撮像部１００により過去の時点から現時点に至るまでに撮影された複数の細胞画像は、バッファ３００に記録される。バッファ３００に記録された複数の細胞画像は、帯域分割部１１０と、領域分割部１２０と、類似度算出部３１０とに転送される。

　帯域分割部１１０は、現時点にて撮影された細胞画像に対し、上記第１の実施の形態と同様の手法により、高域画像Ｈと低域画像Ｌとに帯域を分割する。　
　領域分割部１２０は、現時点にて撮影された細胞画像に対し、上記第１の実施の形態と同様の手法により、領域分割画像Ｋを生成する。

　類似度算出部３１０は、領域分割部１２０により分割された領域分割画像Ｋ内の領域毎に、現時点にて撮影された細胞画像（これ以降、現細胞画像と称する）と、過去の時点にて撮影された細胞画像（これ以降、過去細胞画像と称する）との間で類似度を算出し、この類似度を複数の領域分割画像Ｋに対応する類似度として判別部３３０に転送する。　
　生細胞は運動により位置を移動するが、壊死細胞により構成される領域は、自発的に運動することがないために時間が経過してもその位置を移動することがない。従って、領域毎に現細胞画像と過去細胞画像との間の類似度を算出すると、壊死細胞領域に対する類似度が高くなる。

　以下、類似度算出部３１０の詳細について説明する。　
　類似度算出部３１０は、複数の領域分割画像Ｋ中の領域を注目領域として設定する。類似度算出部３１０は、現細胞画像中の注目領域に対応する（同座標上に位置する）領域をテンプレートとして設定する。類似度算出部３１０は、前記テンプレートを用いて過去細胞画像に対してテンプレートマッチングを行い、過去細胞画像中の各画素の位置座標に対してテンプレートとの類似度（輝度値の相関値）を算出する。類似度は、テンプレート内の輝度分布と、過去細胞画像中の画素を中心とした（テンプレートと同サイズの）周辺領域の輝度分布が似ているほど値が高くなるように定義される。過去細胞画像中の全座標に対して算出した類似度の中で最大値を示す（最も類似性の高い）類似度を注目領域の類似度として決定し、判定部へ転送する。

　類似度には、ＳＳＤ（Sum of Squared Difference ）、ＳＡＤ（Sum of Absolute Difference）等の公知の相関値を用いればよい。本実施の形態では、ＳＳＤを用いている。過去細胞画像は、現細胞画像に対して過去の時点において撮影した画像であれば如何なる画像でもよいが、ここでは現時点の１フレーム前に撮影した画像を過去細胞画像として用いている。　
　テンプレートマッチングを行う際に、過去細胞画像中の全ての画素（座標）に対して類似度を算出するが、注目領域が存在する位置を中心とした所定範囲内に含まれる画素（座標）に限定して類似度を算出することで、処理速度を速くすることができる。

　テクスチャ特徴量算出部３２０は、領域分割画像Ｋ中の複数の領域毎に、帯域分割部１１０により生成された高域画像Ｈからテクスチャ特徴量を算出する。テクスチャ特徴量は、領域内の高域画像Ｈにおける画素値（輝度）分布のヒストグラムに基づく特徴量である。すなわち、テクスチャ特徴量算出部３２０は、上記第１の実施の形態におけるテクスチャ特徴量と相違し、高域画像Ｈの輝度ヒストグラムに基づくテクスチャ特徴量を算出する。テクスチャ特徴量算出部３２０は、テクスチャの複雑さを示す尺度として、輝度ヒストグラムの分散を算出する。

　具体的にテクスチャ特徴量算出部３２０は、領域分割画像Ｋ中の領域を注目領域として設定する、テクスチャ特徴量算出部３２０は、注目領域に対して輝度ヒストグラムＨist[Ｌｖ]を算出する。Ｌｖは、輝度値（高域画像の画素値）を表し、例えば「０－２５５」の範囲で値をとる。　
　テクスチャ特徴量算出部３２０は、注目領域内の画素値平均Ａｖｅを算出し、式（１０）に従って輝度ヒストグラムの分散（Complexity）を算出する。

　ＰＣallは、注目領域内の画素数を表す。輝度ヒストグラムの分散（Complexity）が大きい程、注目領域内に含まれるテクスチャの複雑度が高いので、テクスチャ特徴量算出部３２０は、輝度ヒストグラムの分散（Complexity）の逆数（－１／輝度ヒストグラムの分散（Complexity））を算出して判別部３３０へ転送する。

　判別部３３０は、領域分割画像Ｋ内の各領域に対してテクスチャ特徴量と輝度平均値と類似度との３次元特徴空間内で壊死細胞により形成されている領域か否かを判別する。

　上記第１の実施の形態の判別部１５０は、線形判別分析による判別を行ったが、本第２の実施の形態の判別部３３０は、線形サポートベクタマシン（以下ＳＶＭと称する）及びカーネル関数を用いた非線形判別分析による判別を行う。線形ＳＶＭは、線形判別分析手法であるが、カーネル関数を用いた高次元空間への非線形写像と組み合わせることで、非線形判別分析を行うことが出来る。

　線形ＳＶＭについて説明する。線形ＳＶＭは、学習サンプルデータに基づき、特徴べクトル（本実施の形態では３変数）を２クラスに分類するための超平面（式（１１））を定め、判別対象の領域特徴ベクトルＸｓが入力された際、この領域特徴ベクトルＸｓを式（１２）に示す通りその出力値の正負に基づき壊死領域の判別を行う。

　ｆ（ｘ）＝ｗ・ｘ＋ｂ　　　　　　　　　　　…（１１）　
　ｗ・Ｘｓ＋ｂ＞０→Ｃ１　
　ｗ・Ｘｓ＋ｂ＜０→Ｃ２　　　　　　　　　　…（１２）　
　ここで、ｗは重みベクトル、ｂはバイアス項である。ｆ（ｘ）＝０が識別境界となる超平面となる。重みベクトルｗとバイアス項bとを求めるためには、式（１３）で示される目的関数を最小化すればよい。

　但し、ｙｉ（ｗ・ｘｉ＋ｂ）≧１－εｉ，εｉ≧０　　　…（１３）　
　ここで、ｙｉはラベルであり、１又は－１の値をとる。εはスラック変数であり、超平面で２群（群Ｃ_１又はＣ_２）を分離することができないときに、ある程度の誤判別を認めるためのパラメータである。Ｃはどの程度誤判別を認めるかを示すパラメータであり、ＳＶＭの使用時に実験的に設定する。

　ラグランジェの乗数αに対するラグランジェの未定乗数法を用いて、式（１４）で示されるように目的関数を変形することもできる。　
　maxΣα－（１／２）Σα_ｉα_ｊｙ_ｉｙ_ｊｘ_ｉ ^Ｔｘ_ｊ　
　但し、０≦α_ｉ≦Ｃ,Σα_ｉｙ_ｉ＝０　　　　　　　　　…（１４）
　カーネル関数Φにより高次元特徴空間ηへ写像を行うことにより、線形ＳＶＭを非線形ＳＶＭへと拡張することができる。カーネル関数Φにより写像した高次元特徴空問η内での重みｗは、係数αを用いて、式（１５）のように表される。　
　ｗ＝Σα_ｉｙ_ｉΦ（ｘ_ｉ）　　　　　　　　　　　　　…（１５）　
　識別関数は、式（１６）のように表される。　
　ｆ（Φ（ｘ））＝Σα_ｉｙ_ｉΦ（ｘ）^ＴΦ（ｘ_ｉ）＋ｂ　
　　　　　　　＝Σα_ｉｙ_ｉｋ（ｘ，ｘ_ｉ）＋ｂ　　　　…（１６）　
　このとき、目的関数は、式（１７）のように表される。　
　maxΣα－（１／２）Σα_ｉα_ｊｙ_ｉｙ_ｊｋ（ｘ，ｘ_ｉ）　
　但し、０≦α_ｉ≦Ｃ,Σα_ｉｙ_ｉ＝０　　　　　　　　　…（１７）　
　ここでは、カーネル関数として、一般によく利用されるGaussian カーネルを用いている。

　次に、上記の如く構成された装置の動作について図１０に示す壊死細胞領域検出装置フローチャートに従って説明する。なお、図８に示すステップと同一ステップには同一ステップの符号を付してその詳しい説明は省略する。　
　ＣＰＵは、ステップＳ９０において、撮像部１００により過去の時点から現時点に至るまでに撮影された複数の細胞画像をバッファ３００に記録する。　
　帯域分割部１１０は、ステップＳ２０において、細胞画像から低周波成分により成る低域画像と、高周波成分により成る高域画像とを含む複数の帯域画像に分解する。　
　領域分割部１２０は、ステップＳ３０において、細胞画像群中における所定時点に取得した細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割する。

　テクスチャ特徴量算出部３２０は、ステップＳ１００において、複数の領域分割画像Ｋ中の領域毎に、帯域分割部１１０により生成された高域画像Ｈからテクスチャ特徴量を算出する。テクスチャ特徴量算出部３２０は、上記第１の実施の形態におけるテクスチャ特徴量と相違し、高域画像Ｈの輝度ヒストグラムに基づくテクスチャ特徴量を算出する。ここでは、テクスチャ特徴量算出部３２０は、テクスチャの複雑さを示す尺度として、輝度ヒストグラムの分散を算出する。

　輝度平均算出部１４０は、ステップＳ６０において、領域分割画像Ｋ内の領域ごとに、当該領域に含まれる低域画像Ｌ上の画素値（輝度値）の輝度平均値（平均輝度値）を算出する。　
　ＣＰＵは、ステップＳ１１０において、撮像部１００から入力してバッファ３００に記録された細胞画像が１枚目であるか否かを判定する。この判定の結果、１枚目であれば、ステップＳ１３０に移り、判別部１５０は、領域の類似度は用いずテクスチャ特徴量と輝度平均値だけに基づいて壊死細胞により形成されている領域か否かを判別する。

　細胞画像が２枚目以降であれば、類似度算出部３１０は、ステップＳ１２０において、領域分割部１２０により分割された複数の領域分割画像Ｋ内の領域毎に、現時点にて撮影された現細胞画像と、過去の時点にて撮影された過去細胞画像との間で類似度を算出し、この類似度を領域分割画像Ｋ内の領域に対応する類似度として判別部３３０に転送する。

　判別部３３０は、ステップＳ１３０において、領域分割画像Ｋ内の領域に対して図１０に示すようにテクスチャ特徴量と輝度平均値と類似度との３次元空間内で壊死細胞により形成されている領域か否かを判別する。

　このように上記第２の実施の形態によれば、複数の領域分割画像Ｋ内の領域毎に、現時点にて撮影された現細胞画像と、過去の時点にて撮影された過去細胞画像との間で類似度を算出し、領域分割画像Ｋ内の領域に対してテクスチャ特徴量と輝度平均値と類似度に基づき壊死細胞により形成されている領域か否かを判別するので、上記第１の実施の形態と同様の効果を奏することができる。エントロピーと平均輝度値と類似度に基づく非線形判別を行うことで、壊死領域・非壊死領域を高精度に判別することができる。

Claims

　経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得する画像取得部と、
　前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割する領域分割部と、
　前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解する帯域分割部と、
　前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出する特徴量算出部と、
　前記複数の領域毎に、前記低域画像から輝度平均値を算出する輝度算出部と、
　前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別する判別部と、
を具備することを特徴とする壊死細胞領域検出装置。
　前記領域分割部により分割された前記複数の領域毎に、前記所定時点に取得した前記細胞画像と、前記所定時点よりも以前又は以降の時点のうちいずれか一方又は両方に取得した前記細胞画像との間における局所的な輝度分布の類似度を算出する類似度算出部を具備し、
　前記判別部は、前記複数の領域毎に対してそれぞれ前記テクスチャ特徴量と前記輝度平均値と前記類似度とに基づいて前記壊死細胞により形成されているか否かを判別する、
ことを特徴とする請求項１に記載の壊死細胞領域検出装置。
　前記類似度算出部は、前記複数の領域毎に、当該各領域に対してそれぞれ周辺を対象とした前記細胞画像間における平均二乗誤差を類似度として算出することを特徴とする請求項２に記載の壊死細胞領域検出装置。
　前記帯域分割部は、前記細胞画像に対して平滑化処理を行って前記低域画像を生成する低域画像生成部と、
　前記細胞画像から前記低域画像を減算して前記高域画像を生成する高域画像生成部と、
を含むことを特徴とする請求項１乃至３のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記テクスチャ特徴量は、画素値分布のランダム性に基づく特徴量であることを特徴とする請求項１乃至４のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記テクスチャ特徴量は、画素値分布の複雑度に基づく特徴量であることを特徴とする請求項１乃至４のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記テクスチャ特徴量は、同時生起行列に基づく特徴量であることを特徴とする請求項１乃至６のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記テクスチャ特徴量は、画素値のヒストグラムに基づく特徴量であることを特徴とする請求項１乃至６のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記判別部は、線形判別関数に基づく線形判別処理により前記複数の領域がそれぞれ前記壊死細胞により形成されている領域か否かを判別することを特徴とする請求項１乃至８のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記判別部は、非線形判別処理により前記複数の領域毎がそれぞれ前記壊死細胞により形成されている領域か否かを判別することを特徴とする請求項１乃至８のうちいずれか１項に記載の壊死細胞領域検出装置。
　前記細胞画像群は、明視野顕微鏡により取得されることを特徴とする請求項１乃至１０のうちいずれか１項に記載の壊死細胞領域検出装置。
　コンピュータの処理によって、
　経時的に変化する生細胞を複数の時点でそれぞれ撮像を行って取得した複数の細胞画像から成る細胞画像群を取得し、
　前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割し、
　前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解し、
　前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出し、
　前記複数の領域毎に、前記低域画像から輝度平均値を算出し、
　前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別する、
ことを特徴とする壊死細胞領域検出方法。
　前記複数の領域の分割では、前記各領域毎に、前記所定時点に取得した前記細胞画像と、前記所定時点よりも以前又は以降の時点のうちいずれか一方又は両方に取得した前記細胞画像との間における局所的な輝度分布の類似度を算出し、
　前記領域か否かを判別では、前記各領域に対してそれぞれ前記テクスチャ特徴量と前記輝度平均値と前記類似度とに基づいて前記壊死細胞により形成されている領域か否かを判別する、
ことを特徴とする請求項１２に記載の壊死細胞領域検出方法。
　経時的に変化する生細胞を複数の時点でそれぞれ撮像して取得した複数の細胞画像から成る細胞画像群を取得させる画像取得機能と、
　前記細胞画像群中における所定時点に取得した前記細胞画像に対してそれぞれ局所的画質特性が均一になるように複数の領域に分割させる領域分割機能と、
　前記所定時点で取得した前記細胞画像に対して低周波成分により成る低域画像と高周波成分により成る高域画像とを含む複数の帯域画像に分解させる帯域分割機能と、
　前記複数の領域毎に、前記高域画像からテクスチャ特徴量を算出させる特徴量算出機能と、
　前記複数の領域毎に、前記低域画像から輝度平均値を算出させる輝度算出機能と、
　前記複数の領域毎に対して少なくとも前記テクスチャ特徴量と前記輝度平均値とに基づいて壊死細胞により形成されている領域か否かを判別させる判別機能と、
を実現させる壊死細胞領域検出プログラムを記憶するコンピュータにより処理可能な記憶媒体。
　前記領域分割機能により分割された前記各領域毎に、前記所定時点に取得した前記細胞画像と、前記所定時点よりも以前又は以降の時点のうちいずれか一方又は両方に取得した前記細胞画像との間における局所的な輝度分布の類似度を算出させる類似度算出機能を含み、
　前記判別機能は、前記各領域毎に対して前記テクスチャ特徴量と前記輝度平均値と前記類似度とに基づいて前記壊死細胞により形成されている領域か否かを判別させる、
ことを実現させるための請求項１４に記載の壊死細胞領域検出プログラムを記憶するコンピュータにより処理可能な記憶媒体。