(57)【要約】  時系列画像から細胞塊の複層化の状態を判別して出力する手段を提供する。  画像処理プログラム（ＧＰ）は、時系列画像を取得するステップ（Ｓ１０）と、取得した各画像から細胞塊を抽出するステップ（Ｓ２０）と、抽出された細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を細胞塊全域について算出するステップ（Ｓ５０）と、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別するステップ（Ｓ７０）と、判別結果を出力するステップ（Ｓ８０）とを備え、細胞塊の時系列画像から細胞塊の複層化の状態を判別して出力するように構成される。
2010143420.tif 000002
【特許請求の範囲】
【請求項１】
撮像装置により撮影された細胞塊の時系列画像を取得し、
取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出し、
算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別することを特徴とする細胞塊の状態判別手法。
【請求項２】
複層化の状態判別は、前記小領域内で空間的な輝度値の変化が閾値以上であり、かつ時間的な輝度値の変化が閾値以上である場合に、複層化されていると判別することを特徴とする請求項１に記載の細胞塊の状態判別手法。
【請求項３】
前記空間的な輝度値の変化は、前記小領域内の輝度値の分散または輝度値微分の総和であることを特徴とする請求項２に記載の細胞塊の状態判別手法。
【請求項４】
前記時間的な輝度値の変化は、時系列画像において対応する前記小領域の輝度値の分散または微分値であることを特徴とする請求項２または３に記載の細胞塊の状態判別手法。
【請求項５】
コンピュータにより読み取り可能であり、撮像装置により撮影された画像を取得して画像処理する画像処理装置としてコンピュータを機能させるための画像処理プログラムであって、
撮像装置により撮影された細胞塊の時系列画像を取得するステップと、
取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出するステップと、
算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別するステップと、
判別結果を出力するステップとを備え、
細胞塊の時系列画像から細胞塊の複層化の状態を判別して出力するようにコンピュータを機能させるための画像処理プログラム。
【請求項６】
前記細胞塊の複層化の状態を判別するステップは、
前記小領域内における空間的な輝度値の変化が閾値以上であり、かつ時間的な輝度値の変化が閾値以上である場合に、複層化されていると判別するように構成されることを特徴とする請求項５に記載の画像処理プログラム。
【請求項７】
前記空間的な輝度値の変化は、前記小領域内の輝度値の分散または輝度値微分の総和であることを特徴とする請求項６に記載の画像処理プログラム。
【請求項８】
前記時間的な輝度値の変化は、時系列画像において対応する前記小領域の輝度値の分散または微分値であることを特徴とする請求項６または７に記載の画像処理プログラム。
【請求項９】
前記判別結果を出力するステップは、前記細胞塊の複層化の状態を判別するステップにおいて複層化されていると判断された部位の前記細胞塊における位置情報を出力するように構成されることを特徴とする請求項５〜８のいずれかに記載の画像処理プログラム。
【請求項１０】
前記判別結果を出力するステップは、前記細胞塊の複層化の状態を判別するステップにおいて複層化されていると判断された部位の前記細胞塊に占める大きさの情報を出力するように構成されることを特徴とする請求項５〜９のいずれかに記載の画像処理プログラム。
【請求項１１】
前記時系列の各画像中に複数の細胞塊が含まれる場合において、
前記細胞塊の複層化の状態を判別するステップは、細胞塊ごとに複層化の状態を判別するとともに、複層化部位を持つ細胞塊と複層化部位をもたない細胞塊とに分別し、
前記判別結果を出力するステップは、前記分別された判別結果を出力するように構成されることを特徴とする請求項５〜１１のいずれかに記載の画像処理プログラム。
【請求項１２】
撮像装置により撮影された細胞塊の時系列画像を取得し画像を解析する画像解析部と、前記画像解析部による解析結果を出力する出力部とを備えた画像処理装置であって、
前記画像解析部は、取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出し、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別し、
前記出力部が、前記画像解析部により判別された細胞塊の複層化の状態を出力するように構成したことを特徴とする画像処理装置。
【請求項１３】
前記画像解析部は、前記小領域内における空間的な輝度値の変化が閾値以上であり、かつ時間的な輝度値の変化が閾値以上である場合に、複層化されていると判別するように構成されることを特徴とする請求項１２に記載の画像処理装置。
【請求項１４】
前記空間的な輝度値の変化は、前記小領域内の輝度値の分散または輝度値微分の総和であることを特徴とする請求項１３に記載の画像処理装置。
【請求項１５】
前記時間的な輝度値の変化は、時系列画像において対応する前記小領域の輝度値の分散または微分値であることを特徴とする請求項１３または１４に記載の画像処理装置。
【請求項１６】
前記画像解析部は、複層化されていると判断された部位の前記細胞塊における位置を算出し、前記出力部は、前記画像解析部により算出された複層化の位置情報を出力するように構成したことを特徴とする請求項１２〜１５のいずれかに記載の画像処理装置。
【請求項１７】
前記画像解析部は、複層化されていると判断された部位の前記細胞塊に占める大きさの情報を算出し、前記出力部は、画像解析部により算出された複層化の大きさの情報を出力するように構成したことを特徴とする請求項１２〜１６のいずれかに記載の画像処理装置。
【請求項１８】
前記時系列の各画像中に複数の細胞塊が含まれる場合において、前記画像解析部は、細胞塊ごとに複層化の状態を判別するとともに、複層化部位を持つ細胞塊と複層化部位をもたない細胞塊とに分別し、
前記出力部は、前記分別された判別結果を出力するように構成したことを特徴とする請求項１２〜１７のいずれかに記載の画像処理装置。
【請求項１９】
細胞を培養する培養ステップと、
前記培養ステップにおいて培養される細胞を、請求項１２〜１８のいずれかに記載の画像処理装置を用いて観察し、培養により変化する前記細胞における細胞塊の複層化の状態を判別する判別ステップと、
を備えることを特徴とする細胞塊の製造方法。
【請求項２０】
細胞を培養する培養ステップと、
前記培養ステップにおいて培養される細胞を、撮像装置により撮影し、培養により変化する前記細胞における細胞塊の時系列画像を取得する取得ステップと、
前記取得ステップにおいて取得した前記時系列画像の中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出する算出ステップと、
前記算出ステップにおいて算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別する判別ステップと、
を備えることを特徴とする細胞塊の製造方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、細胞観察において取得された時系列画像から細胞塊の複層化の状態を判別する状態判別手法に関するものである。
【背景技術】
【０００２】
細胞を培養しながら観察する装置の一例として培養顕微鏡が挙げられる。培養顕微鏡は、細胞の培養に好適な環境を形成する培養装置と、培養容器内の細胞を顕微観察する顕微観察系とを備え、生きた細胞を培養しながら、細胞の変化や分裂などを観察できるように構成される（例えば特許文献１を参照）。生細胞の培養過程では、細胞分裂の進行により細胞塊が形成される。細胞分裂の初期過程では、分裂した細胞が単層状態で培地内を水平方向に広がるが、細胞分裂が活発化して細胞塊が成熟してくると、細胞が泡立つように上下方向にも広がり、いわゆる複層化が進行する。【０００３】
培養顕微鏡を用いた従来の細胞観察手法では、細胞塊の複層化の状態判定を、顕微観察画像を目視観察して判別する目視判定や、試薬を投与して着色状態等から判別する試薬判定により行っていた。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−２２９６１９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来行われてきた目視判定の手法では、多数の時系列画像から細胞塊を抽出し複層化された状態を判別するのに、一定の経験を有する知見者が時間をかけて判別する必要があった。特に、各時刻の観察画像に多数の細胞塊が含まれる場合に、個々の細胞塊を識別しながら複層化の状態判別を行うことは大変煩雑な作業であった。また、目視判定では細胞塊において複層化された部位の位置や大きさ（面積や細胞塊に占める比率等）を定量的に捉えることが困難であるという課題があった。さらに、試薬判定の手法では、試薬の投与により細胞に与える化学的・物理的影響や、培養された細胞を利用する際の制約が大きいという課題があった。【０００６】
本発明は、このような課題に鑑みてなされたものであり、細胞に試薬投与によるダメージを与えることなく、撮像装置により撮影された時系列画像から細胞塊の複層化の状態を判別可能な手段を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明を例示する第１の態様に従えば、撮像装置により撮影された細胞塊の時系列画像を取得し、取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出し、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別することを特徴とする細胞塊の状態判別手法が提供される。【０００８】
本発明を例示する第２の態様に従えば、コンピュータにより読み取り可能であり、撮像装置により撮影された画像を取得して画像処理する画像処理装置としてコンピュータを機能させるための画像処理プログラムであって、撮像装置により撮影された細胞塊の時系列画像を取得するステップと、取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出するステップと、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別するステップと、判別結果を出力するステップとを備え、細胞塊の時系列画像から細胞塊の複層化の状態を判別して出力するようにコンピュータを機能させるための画像処理プログラムが提供される。【０００９】
本発明を例示する第３の態様に従えば、撮像装置により撮影された細胞塊の時系列画像を取得し画像を解析する画像解析部と、画像解析部による解析結果を出力する出力部とを備えた画像処理装置が構成される。この画像処理装置において、前記画像解析部は、取得された画像中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出し、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別し、前記出力部が、画像解析部により判別された細胞塊の複層化の状態を出力するように構成される。【００１０】
なお、以上の本発明において、複層化の状態判別は、前記小領域内で空間的な輝度値の変化が閾値以上であり、かつ時間的な輝度値の変化が閾値以上である場合に、複層化されていると判別することが好ましい。この場合において、上記空間的な輝度値の変化は、小領域内の輝度値の分散または輝度値微分の総和であることが好ましく、上記時間的な輝度値の変化は、時系列画像において対応する小領域の輝度値の分散または微分値であることが好ましい。【００１１】
本発明の画像処理プログラムまたは画像処理装置において、複層化されていると判断された部位の細胞塊における位置情報を出力するように構成することが好ましく、複層化されていると判断された部位の細胞塊に占める大きさの情報（面積、体積、これらの比率等）を出力するように構成することが好ましい。また、時系列の各画像中に複数の細胞塊が含まれる場合に、細胞塊ごとに複層化の状態を判別するとともに、複層化部位を持つ細胞塊と複層化部位をもたない細胞塊とに分別し、分別された判別結果を出力するように構成することが望ましい。【００１２】
また、本発明を例示する第４の態様に従えば、細胞を培養する培養ステップと、前記培養ステップにおいて培養される細胞を、上述の画像処理装置を用いて観察し、培養により変化する前記細胞における細胞塊の複層化の状態を判別する判別ステップと、を備えて細胞塊の製造方法が提供される。【００１３】
また、本発明を例示する第５の態様に従えば、細胞を培養する培養ステップと、前記培養ステップにおいて培養される細胞を、撮像装置により撮影し、培養により変化する前記細胞における細胞塊の時系列画像を取得する取得ステップと、前記取得ステップにおいて取得した前記時系列画像の中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出する算出ステップと、前記算出ステップにおいて算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別する判別ステップと、を備えて細胞塊の製造方法が提供される。
【発明の効果】
【００１４】
本発明の細胞塊の状態判別手法、画像処理プログラム及び画像処理装置おいては、細胞塊の小領域の空間的な輝度分布及び時間的な輝度変化から細胞塊の複層化の状態が判別される。従って、本発明によれば、試薬の投与によって細胞にダメージを与えるようなことがなく、撮像装置により撮影された時系列画像から細胞塊の複層化の状態を判別可能な手段を提供することができる。
【図面の簡単な説明】
【００１５】
【図１】画像処理プログラムの概要構成を示すフローチャートである。
【図２】本発明の適用例として示す培養観察システムの概要構成図である。
【図３】上記培養観察システムのブロック図である。
【図４】画像処理装置の概要構成を例示するブロック図である。
【図５】縦軸に小領域の輝度値の空間的分散、横軸に時間をとった場合の空間的分散の時系列変化の様子を示すグラフである。
【図６】セグメンテーション、ラベリングされた細胞塊の状況を例示する模式図である。
【図７】ラベルｎの細胞塊の時系列取得を例示する模式図である。
【図８】画像解析により判別された複層化部位（位置、大きさ、範囲）の出力形態を例示する模式図である。
【図９】細胞塊の製造方法を示すフローチャートである。
【発明を実施するための形態】
【００１６】
以下、本発明を実施するための形態について、図面を参照しながら説明する。本発明の画像処理装置を適用したシステムの一例として、培養観察システムの概要構成図及びブロック図を図２及び図３に示しており、まず培養観察システムＢＳの全体構成について概要説明する。
【００１７】
培養観察システムＢＳは、大別的には、筐体１の上部に設けられた培養室２、複数の培養容器１０を収容保持するストッカー３、培養容器１０内の試料を観察する観察ユニット５、培養容器１０を搬送する搬送ユニット４、システムの作動を制御する制御ユニット６、及び画像表示装置を備えた操作盤７などから構成される。
【００１８】
培養室２は、培養環境を形成する部屋であり、この培養室２に付随して温度調整装置２１、加湿器２２、ＣＯ2ガスやＮ2ガス等のガスを供給するガス供給装置２３、循環ファン２４、培養室２の温度や湿度等を検出する環境センサ２５などが設けられている。ストッカー３は、前後及び上下に仕切られた棚状に形成され各棚に固有の番地が設定される。培養容器１０は、培養する細胞の種別や目的に応じて適宜選択され、例えばディッシュタイプの培養容器に細胞試料が液体培地とともに注入保持される。各培養容器１０にはコード番号が付与され、ストッカー３の指定番地に対応づけて収容される。搬送ユニット４は、培養室２の内部に設けられて上下移動可能なＺステージ４１、前後移動可能なＹステージ４２、左右移動可能なＸステージ４３などからなり、Ｘステージ４３の先端側に培養容器１０を持ちあげ支持する支持アーム４５が設けられている。
【００１９】
観察ユニット５は、試料台１５の下側から試料を照明する第１照明部５１、顕微観察系５５の光軸に沿って試料台１５の上方から試料を照明する第２照明部５２及び下方から試料を照明する第３照明部５３、試料のマクロ観察を行うマクロ観察系５４、試料のミクロ観察を行う顕微観察系５５、及び画像処理装置１００などから構成される。試料台１５には顕微観察系５５の観察領域に透明な窓部１６が設けられている。
【００２０】
マクロ観察系５４は、観察光学系５４ａと、観察光学系により結像された試料の像を撮影するＣＣＤカメラ等の撮像装置５４ｃとを有して構成され、第１照明部５１によりバックライト照明された培養容器１０の上方からの全体観察画像（マクロ像）が取得される。顕微観察系５５は、対物レンズや中間変倍レンズ、蛍光フィルタ等からなる観察光学系５５ａと、観察光学系５５ａにより結像された試料の像を撮影する冷却ＣＣＤカメラ等の撮像装置５５ｃとを有して構成される。対物レンズ及び中間変倍レンズは各々複数設けられ、レンズの組み合わせを変化させることにより任意の観察倍率に設定可能に構成される。顕微観察系５５では、第２照明部５２により照明された細胞の透過像、第３照明部５３により照明された細胞の反射像、第３照明部５３により照明された細胞の蛍光像など、培養容器１０内の細胞を顕微鏡観察した顕微観察像（ミクロ像）が取得される。
【００２１】
画像処理装置１００は、マクロ観察系５４の撮像装置５４ｃ、顕微観察系５５の撮像装置５５ｃにより撮影され、これらの撮像装置から入力された信号を処理して全体観察画像や顕微観察画像などの画像を生成する。また、画像処理装置１００は、これらの観察画像（画像データ）に画像解析を施し、タイムラプス画像の生成、細胞の移動方向予測、細胞の運動状態の解析、細胞塊の複層化の状態の解析などを行う。なお、画像処理装置１００については、後に詳述する。
【００２２】
制御ユニット６は、処理を実行するＣＰＵ６１、培養観察システムＢＳの制御プログラムや制御データ等が設定記憶されたＲＯＭ６２、ハードディスクやＤＶＤ等の補助記憶装置を含み観察条件や画像データ等を一時記憶するＲＡＭ６３などを有し、培養観察システムＢＳの作動を制御する。そのため、図３に示すように、培養室２、搬送ユニット４、観察ユニット５、操作盤７の各構成機器が制御ユニット６に接続されている。ＲＡＭ６３には、観察プログラムに応じた培養室２の環境条件や、観察スケジュール、観察ユニット５における観察種別や観察位置、観察倍率等が設定され記憶される。また、ＲＡＭ６３には、観察ユニット５により撮影された画像データを記録する画像データ記憶領域が設けられ、培養容器１０のコード番号や撮影日時等を含むインデックス・データと画像データとが対応付けて記録される。
【００２３】
操作盤７には、キーボードやスイッチ等の入出力機器が設けられた操作パネル７１、操作画面や観察画像、解析結果等を表示する表示パネル７２が設けられ、操作パネル７１において観察プログラムの設定や条件選択、動作指令等の入力が行われる。通信部６５は有線または無線の通信規格に準拠して構成されており、この通信部６５に外部接続されるコンピュータ等との間でデータの送受信が可能になっている。
【００２４】
このように概要構成される培養観察システムＢＳは、操作盤７において設定された観察プログラムに従ってＣＰＵ６１が各部の作動を制御し、培養容器１０内の試料の撮影を自動的に実行する。観察プログラムがスタートされると、ＣＰＵ６１はＲＡＭ６３に記憶された環境条件に基づいて温度調整装置２１、加湿器２２等の作動を制御する。また、ＲＡＭ６３に記憶された観察条件を読み込み、観察スケジュールに基づいてＸ，Ｙ，Ｚステージ４３，４２，４１を作動させてストッカー３から観察対象の培養容器１０を試料台１５に搬送し、観察ユニット５による観察を開始させる。例えば、観察プログラムにおいて設定された観察が細胞のミクロ観察である場合には、該当する培養容器１０を顕微観察系５５の光軸上に位置決めし、第２照明部５２または第３照明部５３の光源を点灯させて、顕微観察像を撮像装置５５ｃに撮影させる。
【００２５】
以上のように構成される培養観察システムＢＳにおいて、画像処理装置１００は、撮像装置（５４ｃ、５５ｃ）により撮影された時系列画像を取得し、画像に含まれる細胞塊の複層化の状態を判別する機能を有しており、例えば、ｉＰＳ細胞やＥＳ細胞等の解析に好適に利用される。画像処理装置１００は、細胞塊を撮影した時系列画像について、輝度分布の空間的、時間的特長から細胞塊が複層化しているか否かを判断する。この手法は、複層化した部位の画像が、以下の二つの特徴をもつことを利用する。
【００２６】
（Ａ）細胞塊中、複層化が進んだ部分の輝度分布は、複層化していない部分の輝度分布と比較して、空間的な輝度変化（輝度値のばらつき）が大きい。これは、細胞が複層化していない単層領域では、細胞塊が水平方向に広がるため小領域内の空間的な輝度変化がさほど大きくないのに対し、複層化領域では細胞が泡立つように上方に盛り上がるため小領域内の空間的な輝度変化が大きくなるからと考えられる。
【００２７】
（Ｂ）複層化が進んだ部分で生じる輝度変化の時間スケールは、複層化していない部分の輝度変化や、細胞塊境界の変化に伴う輝度変化とくらべて短い。上記のように、複層化領域では、細胞が泡立つように広がるため、時系列画像において、小領域の輝度が比較的短時間で大きく変動するからと考えられる。
【００２８】
具体的に、培養中の細胞塊を撮影した時系列画像について、細胞塊の大きさに比較して十分に小さい小領域を指定し、小領域の空間的な輝度変化とその時間的変化をプロットすると、図５に示すようになる。この図は、縦軸に小領域の空間的な輝度値の分散（例えば縦横２０×２０画素（細胞２〜３個分の大きさ）の小領域の輝度値の分散）、横軸に時間をとり、小領域の空間的な輝度値分散の変化状態を、複層化された状態の小領域（実線）と、他の小領域（点線および二点鎖線）とについて模式的に示したものである。
【００２９】
図５中の実線は、複層化された部位の小領域の特性を示す。複層化された部位の小領域では、泡立つような細胞の増殖により、空間的な輝度値の変化（輝度分散値）が大きく、かつ、比較的短い周期で大きく波打つように変化する。
【００３０】
図５中の点線は、細胞塊の境界（輪郭部）を含む小領域の特性を示す。細胞塊の境界では輪郭線を挟んで輝度が大きく変化するため、空間的な輝度値の変化（輝度分散値）は比較的大きくなる。しかし、細胞塊が成長または縮退して細胞塊の境界が移動したとしても、その移動速度は遅く、輝度値の時間的な変化は緩慢かつ小さなものとなる。
【００３１】
図５中の一点鎖線は、複層化していない部位の小領域の特性を示す。単層状態の部位の小領域では、細胞塊が広がっても空間的な輝度値の変化は小さく、時間的にも輝度値の変化は小さいものとなる。
【００３２】
本発明は、上記のような複層化領域の画像上の特徴に着目し、細胞塊の時系列画像を画像処理することにより複層化の状態を判別する。画像処理装置１００は、細胞塊内部の小領域の空間的な輝度分布及び時間的な輝度変化を細胞塊全域（解析範囲がマウス等により指定されている場合には指定された解析範囲内）について算出し、算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別する。
【００３３】
図４に画像処理装置１００のブロック図を示し、図１に複層化の判別処理を行う画像処理プログラムＧＰのフローチャートを示す。
【００３４】
画像処理装置１００は、撮像装置（５５ｃ，５４ｃ）により撮影された細胞塊の時系列画像を取得して画像を解析する画像解析部１２０と、画像解析部１２０により解析された解析結果を出力する出力部１３０とを備えて構成され、画像解析部１２０による解析結果、例えば複層化が進んでいると判断された部位の位置や大きさの情報（面積、体積、比率等）、複層化された部位を含む細胞塊と含まない細胞塊との判別などが出力部１３０から出力され表示パネル７２等に表示されるように構成される。
【００３５】
画像処理装置１００は、ＲＯＭ６２に設定記憶された画像処理プログラムＧＰがＣＰＵ６１に読み込まれ、ＣＰＵ６１によって画像処理プログラムＧＰに基づく処理が順次実行されることによって構成される。換言すれば、画像処理プログラムＧＰはハードウェア資源であるＣＰＵ６１（コンピュータ）を画像処理装置１００として機能させるためのソフトウェアである。
【００３６】
画像解析部１２０は、画像処理プログラムＧＰに基づき、撮像装置（説明ではミクロ系の撮像装置５５ｃとする）により撮影され、ＲＡＭ６３に記録された細胞塊の時系列画像を以下のように画像処理する。なお、撮像装置５５ｃにより撮影された画像を順次リアルタイムで処理するようにしてもよい。
【００３７】
画像解析部１２０は、ＲＡＭ６３に記憶された細胞塊の時系列画像を順次取得し（ステップＳ１０）、時刻ｔ1・・・ｔnの各画像に対し、レベルセット（Level Set）、及び分散フィルタによって細胞塊のセグメンテーションを行う（ステップＳ２０）。次いで、図６に示すように、セグメント化された細胞塊ＭＣに対してラベリングを行い（ステップＳ３０）、各画像間での細胞塊の対応をとる。例えば、各画像においてラベリング１，２，３…を施した細胞塊ＭＣについて、像が重なるラベル同士を同じ細胞塊として対応づける。
【００３８】
次いで、細胞塊ＭＣが移動する場合の効果を小さくするため、同一ラベルの細胞塊の位置合わせを行う。位置合わせの基準には、細胞塊の重心位置や、外接矩形の頂点位置などを用い、画像間での移動量を算出して位置合わせすることができ、また図形モーメントの相関が最大（差分が最小）となるように細胞塊の回転角度を揃えることができる。時刻ｔ1の第１像と時刻ｔ2の第２像との差分が最小（相関値が最大）となる位置及び角度で位置合わせを行うようにしてもよい。
【００３９】
このように対応付けされ、位置合わせされた同一ラベルの細胞塊ＭＣの時系列画像（図７はラベルｎの細胞塊の時系列画像を示す）に対して、ステップＳ４０において、画像を形成する画素を中心とする多数の「小領域」Ａに分割する。小領域Ａは、細胞塊ＭＣの大きさに比べて十分に小さい局所領域として設定され、例えば、１０×１０〜３０×３０画素程度（細胞２〜５個程度の大きさ）に設定される。
【００４０】
そして、ステップＳ５０において画像の解析が行われる。画像解析は、小領域の空間的な解析５０Ａと、小領域の時系列解析５０Ｂとにより構成される。
【００４１】
（５０Ａ：小領域の空間的な解析）
小領域の空間的な解析では、各時刻画像における各小領域について、小領域Ａの輝度分布から空間的な輝度値の変化（ばらつき）を算出する。空間的な輝度値変化の評価指標として、小領域Ａの輝度値の分散や、空間方向に対する輝度値の微分の総和などが例示され、この値（スコア）ａによって空間的な輝度値変化の特徴が定量的に表現される。スコアａが大きな小領域は、輝度値が空間的に大きく変化している領域であり、細胞塊中で複層化している部位や、細胞塊内外の境界の部分が該当する。
【００４２】
（５０Ｂ：小領域の時系列解析）
小領域の時系列解析では、時系列画像における小領域Ａの時間的な輝度値の変化（ばらつき）を算出する。時間的な輝度値変化の評価指標として、時系列画像間での小領域Ａの輝度値の分散や、微分値などが例示され、この値（スコア）ｂによって時間的な輝度値変化の特徴を定量的に表現することができる。スコアｂが大きな小領域は輝度値が時間的に大きく変化している領域であり、細胞塊中で複層化している部分が該当する。
【００４３】
スコアｂを算出する場合の「小領域の輝度値」は、５０Ａで算出したスコアａ（輝度値の分散、微分の総和）を用いることができ、このような構成により、ＣＰＵ６１の処理負担を軽減して演算処理を高速化することができる。なお、小領域の輝度値は、小領域Ａを形成する複数画素の少なくともいずれかの画素の輝度値、あるいは複数画素の適宜な平均輝度値を用いることもできる。
【００４４】
ステップＳ６０では、上記のようにして得られた空間的な輝度値の変化のスコアａ及びスコアｂを、細胞塊ＭＣの各小領域に割り当て、細胞塊全域について、スコアａ及びｂを付したマップを形成する。
【００４５】
そしてステップＳ７０において、細胞塊ＭＣの各小領域について、スコアａが所定の閾値以上であり、かつスコアｂが所定の閾値以上であるか否かを判断し、ステップＳ７５において、スコアａ及びｂがともに閾値以上である領域を複層化された領域、そうでない領域をその他の領域に分別する（図５を併せて参照）。なお、ステップＳ７０において判断する際の小領域のスコアａは、時系列画像における各時刻の小領域のスコアａ1・・・ａnの総和Σ（ａ1・・・ａn）や平均値等を用いることができる。上記各閾値は、観察対象となる細胞の種別や観察条件、解析条件などに応じて適宜な値が設定される。また複層化領域の判断は、スコアａ×ｂが所定の閾値以上である領域を複層化領域としてもよい。
【００４６】
ステップＳ８０では、ステップＳ７５の分別結果に基づき、細胞塊の複層化の状態を出力部１３０から出力して表示パネル７２等に表示させる。細胞塊の複層化の状態は、例えば図８に模式図を示すように、ラベルｎの細胞塊ＭＣの画像中に、複層化されたと判断された領域をハッチングや色分け等により識別可能に表示する。また、当該ラベルｎの細胞塊に占める複層化領域の大きさの情報（面積、体積、これらの比率等）を数値データで表示する。
【００４７】
図８では、マウス等により選択された特定の（ラベルｎの）細胞塊を拡大表示した構成例を示すが、図６のように画像中に複数の細胞塊ＭＣが含まれるような場合には、表示画面を全体像に切り換えることにより、複層化された部位を含む細胞塊と含まない細胞塊とを識別可能に表示する。
【００４８】
このような全体画像における識別表示の具体的な形態として、例えば、画像中の各細胞塊について図８と同様の表示を行い、あるいは、複層化された部位を有する細胞塊を黄色、複層化部位の無い内細胞塊を青色のように分別表示し、または、各細胞塊に占める複層化部位の面積比に応じて、面積比が高い細胞塊ほど赤く、面積比が低くなるにつれて黄色〜緑〜青のように識別表示する表示形態が例示される。なお、解析結果をプリンタやＲＡＭ６３、磁気記録媒体等に出力して記録させ、あるいは通信部６５を介してシステム外部に出力するように構成してもよい。
【００４９】
これにより、観察者は画像に含まれる細胞塊の複層化の状態を定量的に目視判断することができる。このように、画像処理装置１００においては、細胞塊の小領域の空間的な輝度分布及び時間的な変化から細胞塊の複層化の状態が判別される。従って、画像処理装置１００による細胞塊の状態判別手法によれば、試薬の投与によって細胞にダメージを与えるようなことがなく、撮像装置（５５ｃ，５４ｃ）により撮影された時系列画像から細胞塊の複層化の状態を判別可能な手段を提供することができる。
【００５０】
なお、以上説明した実施形態では、培養観察システムＢＳにおいて撮像装置に撮影され、ＲＡＭ６３に記憶された時系列画像（画像データ）を読み出して複層化状態を解析する構成を例示したが、撮像装置により撮影された画像を逐次リアルタイムで解析するように構成してもよく、また、他の観察システムにおいて撮影され磁気記憶媒体等に記録された時系列画像を読み出して複層化状態を解析するように構成してもよい。また、オペレータが時系列画像に含まれる細胞塊の所定範囲（特定の細胞塊や、細胞塊における特定部位）をマウス等により解析範囲として設定し、設定された解析範囲について画像処理装置が複層化状態の解析を実行するように構成してもよい。
【００５１】
次に、本発明の実施形態に係る細胞塊の製造方法について、図９を参照して説明する。この製造方法は、基本的には、細胞を培養する培養ステップ（Ｓ１１０）と、この培養ステップにおいて培養される細胞を、前述した画像処理装置を用いて観察し、培養により変化する細胞における細胞塊の複層化の状態を判別する判別ステップ（Ｓ１２０−Ｓ１４０）とを備える。
【００５２】
より具体的には、細胞を培養する培養ステップ（Ｓ１１０）と、この培養ステップにおいて培養される細胞を、撮像装置により撮影し、培養により変化する細胞における細胞塊の時系列画像を取得する取得ステップ（Ｓ１２０）と、この取得ステップにおいて取得した時系列画像の中の細胞塊における小領域の空間的な輝度値の分布及び時間的な輝度値の変化を算出する算出ステップ（Ｓ１３０）と、この算出ステップにおいて算出された輝度値の空間分布及び時系列変化に基づいて細胞塊の複層化の状態を判別する判別ステップ（Ｓ１４０）と、所定の基準に基づいて、細胞塊を選別する選別ステップ（Ｓ１５０）と、選別した細胞塊を採取、保存する採取保存ステップ（Ｓ１６０）とを備えて、製造方法が構成される。なお、培養される細胞は、ヒト、ウシ、ウマ、ブタ、マウス等の動物由来の細胞であってもよいし、植物由来の細胞であってもよい。また、細胞塊の保存は凍結保存であってもよい。
【符号の説明】
【００５３】
Ａ  小領域                          ＢＳ  培養観察システム
ＧＰ  画像処理プログラム            ＭＣ  細胞塊
５  観察ユニット                    ６  制御ユニット
５４  マクロ観察系                  ５４ｃ  撮像装置
５５  顕微観察系                    ５５ｃ  撮像装置
６１  ＣＰＵ（コンピュータ）        ６２  ＲＯＭ
６３  ＲＡＭ                        １００  画像処理装置
１２０  画像解析部                  １３０  出力部
【図１】
2010143420.tif 000003
【図２】
2010143420.tif 000004
【図３】
2010143420.tif 000005
【図４】
2010143420.tif 000006
【図５】
2010143420.tif 000007
【図６】
2010143420.tif 000008
【図７】
2010143420.tif 000009
【図８】
2010143420.tif 000010
【図９】
2010143420.tif 000011
【国際調査報告】
2010143420.tif 000012
2010143420.tif 000013