(57)【要約】      （修正有）【課題】培養器内で培養環境にあるカルスの状態がフライアブルであるか、コンパクトであるかについての客観的な指標を与え、簡便で再現性よく客観的にカルスの状態を判別できるようにする。【解決手段】カルスのカラー画像をグレースケール変換した画像と、カルスのカラー画像の色成分画像をグレースケール変換した画像とを生成し、各々のグレースケール変換した画像についてそれぞれ濃度共起行列法によるテクスチャ解析を行って画像のエントロピー値を算出し、カラー画像をグレースケール変換した画像と少なくとも１種類の色成分画像をグレースケール変換した画像とのエントロピー値の相関対比を行う。カルスの画像を良好に取得するには、カルスを収容する培養器外から載置面内の側方の照明光を与えた場合と、斜め前方からの照明光を与えた場合とについて撮影し取得された画像についてそれぞれテクスチャ解析を行う。【選択図】図１０
2008261631.tif 000002
【特許請求の範囲】
【請求項１】
植物培養細胞塊であるカルスを巨視的な大きさで撮影したカルスのカラー画像からカルスの状態がフライアブルであるかコンパクトであるかを判別する植物培養細胞塊の状態を判別する方法であって、前記カルスのカラー画像における背景部分を除去したカルスのカラー画像を生成することと、該背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成することと、前記カラー画像をグレースケール変換した画像と前記カラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換したカルスの画像のエントロピー値を算出することと、算出された該エントロピー値を対比することによりカルスの状態がフライアブルであるかコンパクトであるかを定量的に判別することと、からなることを特徴とする植物培養細胞塊の状態を判別する方法。
【請求項２】
前記植物培養細胞塊の状態を判別する対象となる巨視的な大きさで撮影したカルスの画像を撮影装置により撮影により取得する際に、載置面上に配置された透明な培養器内に培養環境において収容されたカルスに対し載置面内で側方からの照明光を前記培養器外から照射し撮影して側方照明によるカルスのカラー画像を取得するとともに、前記載置面上に配置された透明な培養器内に培養環境において収容されたカルスに対し載置面に斜め前方から照明光を前記培養器外から照射し撮影して斜め方向照明によるカルスのカラー画像を取得し、前記側方照明によるカルスのカラー画像と斜め方向照明によるカルスのカラー画像との両方についてそれぞれ濃度共起行列法によるテクスチャ解析を行うエントロピー値を対比するようにしたことを特徴とする請求項１に記載の植物培養細胞塊の状態を判別する方法。
【請求項３】
培養器載置台と、該培養器載置台上に配置され植物培養細胞塊を収容する培養器に対し側方および斜め方向の照明光を前記培養器の外方から選択的に照射するための照明手段と、前記培養器載置台上に配置された培養器内に収容された植物培養細胞塊を巨視的な大きさで撮影するための撮影部と、撮影により得られたカラー画像を取得しその背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成する画像取得前処理部と、該画像取得前処理部で生成されたカラー画像をグレースケール変換した画像とカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換したカルスの画像のエントロピー値を算出し相関対比を行うテクスチャ解析処理部と、からなることを特徴とする植物培養細胞塊の状態を判別するための装置。
【請求項４】
植物培養細胞塊であるカルスを巨視的な大きさで撮影したカラー画像における背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成し、前記カラー画像をグレースケール変換した画像と前記カラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換した画像のエントロピー値を算出し、前記カラー画像をグレースケール変換した画像と少なくとも１種類の色成分画像をグレースケール変換した画像とのエントロピー値の相関対比を行うことをコンピュータ上で実行するようにしたことを特徴とする植物培養細胞塊の状態を判別するためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、植物培養細胞塊の状態を判別する方法、そのための装置および植物培養細胞塊の状態を判別するためのプログラムに関する。
【背景技術】
【０００２】
植物細胞培養の技術は、医学、農業、食品製造等に広く利用される植物の二次代謝産物の効率的な生産、種苗生産技術への利用、種子のできにくい植物や種子の発芽率が低いため長期保存に適していない植物の保存等に利用されている。植物は二次代謝産物として、アルカロイド、テルペノイド、フラボノイド等さまざまな有用物質を生産するが、生産する量が少ないので、それらの成分を抽出するには大規模な圃場が必要になり環境を破壊することにもなる。そこで、植物細胞培養技術を用いることにより二次代謝産物を効率的に得ることができる。また、植物細胞培養は不定胚生産への利用など種苗大量生産技術としても有用である。
【０００３】
このような植物細胞培養技術を最大限に利用するために、培養細胞を高品質に維持することが必要であり、特に植物の培養細胞は未分化細胞塊であるカルスを誘導し培養することから始められることが多く、このカルスの適正な評価、選抜が培養による増殖効率や代謝物の生産効率に大きな影響を与える。さらに、カルスの状態の評価は培養の再現性を高めるために不可欠である。
【０００４】
カルスの状態を評価するにあたり、カルスを分類するのに使われるフライアブル・カルス（ｆｒｉａｂｌｅ  ｃａｌｌｕｓ）、コンパクト・カルス（ｃｏｍｐａｃｔ  ｃａｌｌｕｓ）という指標が一般的に用いられおり、この指標での分類により評価を行うのが有効であるとされている。カルスの状態を表す「コンパクト」は細胞同士が硬くしっかりとくっついて細胞塊になっている状態であり、一方「フライアブル」は細胞が緩く集まってもろい、ないし柔らかい細胞塊になっている状態を表しており、培養細胞の利用技術においては、細胞塊のこのようなフライアブル、コンパクトの状態を判別することが重要な意義を有している。
【０００５】
しかしながら、実際にフライアブル・カルス、コンパクト・カルスの判別を行うにあたり、現状では目視で判別することになり、そのため評価する者の主観による個人差が避けられず、客観性、一般性をもつ評価を行うのは困難であった。
【０００６】
培養細胞を客観的に評価するのに培養細胞を撮影し得られた画像から培養細胞の質を表す情報を抽出する画像解析の手法が一般的に試みられており、これには顕微鏡撮影によるミクロ画像を用いたものと、巨視的状態で撮影されたマクロ画像を用いたものとがある。ミクロ画像は多くの情報を含んでおり、解析手法も多様な形が利用可能であるが、試料のサンプリングに問題があり、どうしても試料を培養器外に取り出した破壊的な形での評価とならざるを得ない。これに対し、マクロ画像の場合は試料が培養器内にある状態で撮影することができ、サンプリングが容易で、培養細胞を客観的で単純かつ非破壊的に評価できる可能性をもつ。
【０００７】
非特許文献１には、細胞懸濁培養における培養細胞の評価の重要性、細胞懸濁液のマクロ画像のテクスチャ解析がその中に含まれる懸濁培養細胞の評価に有効であることについて記載されており、培養器の底面から蛍光灯からの透過光下で撮影した懸濁培養細胞のマクロ画像における細胞塊に対応する画素の総数、画像中の全画素のＢ成分の値の合計値を用いて細胞懸濁液の圧縮細胞量を推定することができること、マクロ画像から濃度共起行列法により抽出したテクスチャ特徴量と懸濁培養細胞の不定胚形成能力との対応を調べることにより、細胞懸濁液中に不定胚になる能力を持つ比較的大きい細胞塊が多く含まれるか、比較的細かい細胞塊が多く含まれるかという細胞懸濁液中の細胞塊の大きさの分布の違いがテクスチャ特徴量としてのエントロピーの差として現れることが示されている。
【０００８】
しかしながら、非特許文献１において、液体培地中で培養している懸濁培養細胞の評価にテクスチャ解析が有効であり、評価について細胞塊（カルス）の大きさや色の分布の違いによるテクスチャの違いを評価するものではあるが、培養細胞におけるフライアブル、コンパクトの状態の評価を行うものではなく、評価の有効性について十分なものではなかった。
【非特許文献１】Y.Ibaraki, K.Kurata「Application of image analysis to plant cell suspensioncultures」，Computers and Electronics in Agriculture 30（2001）, pp193-203
【発明の開示】
【発明が解決しようとする課題】
【０００９】
前述のように、従来の培養細胞の評価は目視での観察、顕微鏡を用いた観察で行われるか、培養細胞を撮影して得られたミクロ画像あるいはマクロ画像を用いて画像解析を行うという形態で行われており、このうち培養細胞の目視や顕微鏡を用いた観察によるものでは評価する者の主観による個人差が避けられず、客観的で再現性のある評価がなされないものであり、培養細胞の顕微鏡撮影によるミクロ画像を用いた画像解析によるものは撮影のための装置が煩雑であるとともに、試料を培養器外に取り出した破壊的な形での評価とならざるを得ないものであって、培養器内にある状態での細胞の評価ができないものであった。また、マクロ画像を用いた画像解析によるものでは、カルスの大きさや色の分布の違いによるテクスチャの違いの評価は行えるが、培養細胞の評価において重要なフライアブル、コンパクトの状態の評価を行ってはいなかった。このようなことから、カルスの状態に客観的な指標を与え、簡便で再現性よく客観的に評価できるようにすることが望まれていた。
【課題を解決するための手段】
【００１０】
本発明は、前述した課題を解決すべくなしたものであり、本発明による植物培養細胞塊の状態を判別する方法は、植物培養細胞塊であるカルスを巨視的な大きさで撮影したカルスのカラー画像からカルスの状態がフライアブルであるかコンパクトであるかを判別する植物培養細胞塊の状態を判別する方法であって、前記カルスのカラー画像における背景部分を除去したカルスのカラー画像を生成することと、該背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成することと、前記カラー画像をグレースケール変換した画像と前記カラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換したカルスの画像のエントロピー値を算出することと、算出された該エントロピー値を対比することによりカルスの状態がフライアブルであるかコンパクトであるかを定量的に判別することと、からなるものである。
【００１１】
前記植物培養細胞塊の状態を判別する対象となる巨視的な大きさで撮影したカルスの画像を撮影装置により撮影により取得する際に、載置面上に配置された透明な培養器内に培養環境において収容されたカルスに対し載置面内で側方からの照明光を前記培養器外から照射し撮影して側方照明によるカルスのカラー画像を取得するとともに、前記載置面上に配置された透明な培養器内に培養環境において収容されたカルスに対し載置面に斜め前方から照明光を前記培養器外から照射し撮影して斜め方向照明によるカルスのカラー画像を取得し、前記側方照明によるカルスのカラー画像と斜め方向照明によるカルスのカラー画像との両方についてそれぞれ濃度共起行列法によるテクスチャ解析を行うエントロピー値を対比するようにしてもよい。
【００１２】
また、本発明による植物培養細胞塊の状態を判別するための装置は、培養器載置台と、該培養器載置台上に配置され植物培養細胞塊を収容する培養器に対し側方および斜め方向の照明光を前記培養器の外方から選択的に照射するための照明手段と、前記培養器載置台上に配置された培養器内に収容された植物培養細胞塊を巨視的な大きさで撮影するための撮影部と、撮影により得られたカラー画像を取得しその背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成する画像取得前処理部と、該画像取得前処理部で生成されたカラー画像をグレースケール変換した画像とカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換したカルスの画像のエントロピー値を算出し相関対比を行うテクスチャ解析処理部と、からなるものである。
【００１３】
また、本発明による植物培養細胞塊の状態を判別するためのプログラムは、植物培養細胞塊であるカルスを巨視的な大きさで撮影したカラー画像における背景部分を除去したカルスのカラー画像をグレースケール変換した画像と該背景部分を除去したカルスのカラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とを生成し、前記カラー画像をグレースケール変換した画像と前記カラー画像の少なくとも１種類の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法によるテクスチャ解析を行い各々のグレースケール変換した画像のエントロピー値を算出し、前記カラー画像をグレースケール変換した画像と少なくとも１種類の色成分画像をグレースケール変換した画像とのエントロピー値の相関対比を行うことをコンピュータ上で実行するためのものである。
【発明の効果】
【００１４】
本発明では、カルスのカラー画像をグレースケール変換した画像と、カラー画像の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法を用いたテクスチャ解析を行い画像の特徴量としてのエントロピー値を算出し、同じ色成分間または異なる色成分間でのエントロピー値を対比させることによりフライアブル・カルスとコンパクト・カルスとを分離させることができる。また、この手法でカルスの状態の判別を行う上で、十分な明るさの照明光を側方から照明して撮影した画像と斜め方向から照射して撮影した画像との両方についてテクスチャ解析を行うことによりフライアブル・カルスとコンパクト・カルスとの分離が良好になされる。
【発明を実施するための最良の形態】
【００１５】
本発明においては、培養細胞を撮影して得られた画像についてテクスチャ解析を行い、画像として表された培養細胞についてのフライアブル、コンパクトの状態の評価を行う。そこで、テクスチャ解析について概略説明する。
［テクスチャ解析］
画像についてコンピュータで画像処理を行う上で、画像は多数の微小な画素が配列されたものとして扱われ、画像データとしては画像を構成する各画素について濃度値が与えられている。この濃度値を統計的に処理することで、画像の特徴を表すことができる。簡単な特徴量として画像の濃度ヒストグラムの平均、標準偏差があるが、これらは撮影条件の影響を受けやすい。そこで、濃度値の配列に注目した画像の特徴量の計算方法として、テクスチャ解析を用いることを考える。
【００１６】
テクスチャは簡単に言えば模様である。テクスチャには統計的なものと、構造的なものとの二つのレベルがあり、構造性が強いものほど解析しやすいが、自然界に存在するテクスチャは統計的なものがほとんどであり、あらゆるテクスチャの違いを一般的に解析する方法はない。テクスチャ解析手法を大別すると、統計的特徴の抽出、局所的幾何学的特徴による解析、モデル当てはめによる解析、構造的解析に分けられる。本発明ではテクスチャ解析手法として代表的アルゴリズムである濃度共起行列法を用いるが、これは濃度共起行列を用いて種々の特徴量を計算し、テクスチャの識別を行うものである。
【００１７】
濃度共起行列とは、濃度（ｉ）、濃度（ｊ）の画素対の出現確率を距離・方向別に行列Ｐ（ｉ，ｊ）として整理したものであり、画像の濃度階調数をｎ階調とすると、行列はｎ行ｎ列の正方形行列になる。通常、距離ｄは近接から遠方へ数種類、角度θは０°、４５°、９０°、１３５°の４方向をとる。ある画像のテクスチャ情報はその画像中の画像濃度値の空間的相互関係に一様に含まれていると仮定すると、その画像のテクスチャ特徴量は全て濃度共起行列から導き出される。図１（ａ）は濃度階調数ｎ＝３（濃度値０，１，２）とした場合に原画像の３×３画素の濃度値の例を示しており、図１（ｂ）はこの原画像についてｄ＝１，θ＝０°、ｎ＝３とした時の濃度共起行列を示している。
【００１８】
濃度共起行列から画像についてのいくつかの特徴量が求められるが、本発明では画像の乱雑さを反映する特徴量であるエントロピー（以下、Entropy記す）を式（１）により計算する。
【００１９】
【数１】
2008261631.tif 000003
ここでは画像の等方性を仮定し、Entropyは水平方向（θ＝０）のみ計算する。また、距離ｄはｄ＝１，２，５，１０，２０で計算する。距離ｄ＝ｋの時のEntropyをEntｋとし、距離ｄ＝２０におけるEntropyと距離ｄ＝１におけるEntropyとの差をΔEntで表す。また、Entropyは、ＲＧＢカラー画像を輝度によりグレースケール変換した画像（Brig
htness image）、赤成分画像（R image）、緑成分画像（G image）、青成分画像（B image）についてのEntropyをそれぞれEntｋBright、EntｋR、EntｋG、EntｋBと表す。
【００２０】
本発明において、このようなエントロピーに着目したテクスチャ解析によりカルスの状態を定量化する手法を実現する上で、カルスが撮影された画像を用いたテクスチャ解析を行う。テクスチャ解析を行う対象となる画像はカルスが撮影され画像データとして利用可能なものであればよいのであり、本発明では、論文から収集した既存のカルス画像についてテクスチャ解析を行ってカルスの状態を定量化する手法を定式化し、さらに撮影により得られたカルス画像についてテクスチャ解析を行う上で撮影条件がテクスチャ解析に及ぼす影響を考慮してテクスチャ解析によりカルスの状態を判別する上で良好な撮影条件を求めるものである。
〔１〕既存の画像についてのテクスチャ解析
カルスの画像が掲載されておりフライアブル・カルスとコンパクト・カルスとに分類されている文献からカルスの画像を収集し、それらの画像について濃度共起行列法によるテクスチャ解析を行って、フライアブル・カルス、コンパクト・カルスの特徴量を見出し状態判別の指標を得ることを考える。
【００２１】
文献は論文検索サイトのグーグル・スコラーGoogle Scholar）で検索したＰＤＦ形式の論文からカルスの画像が載っている論文を選び出し、それらの画像を画像抽出ソフトウェアで画像データとして取り込み、フライアブル・カルスとコンパクト・カルスとの巨視的外観の違いを定量化するテクスチャ解析のサンプルとして用いる。その際画像解析ソフトウェアによりカルスの画像における背景を自動選択ツールで選択して取り除き、細胞の部分だけをグレースケール画像に変換し、テクスチャ解析を行う。ここでは、
１）ＲＧＢカラー画像を輝度によりグレースケール変換した画像（Brightness image）
２）Ｒ成分画像をグレースケール変換した画像（R image）
３）Ｇ成分画像をグレースケール変換した画像（G image）
４）Ｂ成分画像をグレースケール変換した画像（B image）
の４種類の８ビットのグレースケール画像に変換したものに対してテクスチャ解析を行う。
【００２２】
カルスの画像についてテクスチャ解析を行う形として、（ａ）同一論文上でフライアブル・カルスあるいはコンパクト・カルスとされている画像についてカクスチャ解析を行い、エントロピーを比較すること、（ｂ）取得したフライアブル・カルスおよびコンパクト・カルスをグループ分けし、グループ別にテクスチャ解析を行い対比することの２つに分けて行った。
（ａ）同一論文上でフライアブル・カルスあるいはコンパクト・カルスとされている画像についてのカクスチャ解析
１２の論文からカルスの画像を抽出し、同一論文上でそれぞれフライアブル・カルスあるいはコンパクト・カルスのいずれであるかが規定されているカルスの画像についてテクスチャ解析を行って式（１）によるEntropyを算出しフライアブルとコンパクトとの間でのエントロピーの比較を行い、その結果は表１、表２のようになっている。表１、表２で番号１〜１２はカルスの画像を抽出した論文に付したものである。
【００２３】
【表１】
2008261631.tif 000004
【００２４】
【表２】
2008261631.tif 000005
表１の結果では、１２の論文中１０の論文でコンパクト・カルスの画像の方がエントロピーが高くなっているのがわかる。フライアブル・カルスの画像の方がエントロピーが高くなった２つの論文のうち１つ（フライアブル・カルスの５−１，５−２、コンパクト・カルスの５）はＢ成分画像のみフライアブル・カルスの方がコンパクト・カルスよりエントロピーが高くなっており、この論文のカルスは他の論文のように黄味がかった色ではなく、Ｂ成分を多く含む緑色であったためにエントロピーが増大したと考えられる。もう１つの論文（１０）ではフライアブル・カルスとコンパクト・カルスの撮影条件が一定でなく、この論文のフライアブル・カルスがコンパクト・カルスに近い状態のものであったことが考えられる。このことから、カルスの撮影条件を一定にすること、カルスの色に応じて抽出する色成分を変えることで、他の１０の論文の場合と同様にコンパクト・カルスの画像の方がフライアブル・カルスの画像よりエントロピーが高くなるものと推察される。
（ｂ）グループ分けしたカルスの画像についてのテクスチャ解析
取得されたフライアブル・カルスの画像およびコンパクト・カルスの画像をそれぞれ次の３つのグループに分ける。グループ分けは論文作者によるものと、目視によるものとを合わせた形でおこなう。
フライアブル・カルス
Ｆ１：論文の作者がフライアブルとし、目視により明らかにフライアブルとされるもの
Ｆ２：論文の作者がフライアブルとしたもののうちＦ１でもＦ３でもないもの
Ｆ３：論文の作者がフライアブルとしたもので、目視によりコンパクトとされるもの
コンパクト・カルス
Ｃ１：論文の作者がコンパクトとし、目視により明らかにコンパクトとされるもの
Ｃ２：論文の作者がコンパクトとしたもののうちＣ１でもＣ３でもないもの
Ｃ３：論文の作者がコンパクトとしたもので、目視によりフライアブルとされるもの
このようにグループ分けしたＦ１、Ｆ２、Ｃ１、Ｃ２のカルスの画像についてテクスチャ解析を行って算出されたエントロピーをそれぞれ表３、表４、表５、表６に示す。
【００２５】
【表３】
2008261631.tif 000006
【００２６】
【表４】
2008261631.tif 000007
【００２７】
【表５】
2008261631.tif 000008
【００２８】
【表６】
2008261631.tif 000009
表３〜６においては、B image、G image、Brightness imageでは、ほぼ全てのカルスの画像においてＣ１のグループのEnt20とΔEntとがＦ１のグループのEn20とΔEntとよりも値が大きくなっている。特にB imageとG imageとでその傾向が顕著になっている。ここで、Ｆ１のＮｏ．４およびＮｏ．５ではＦ１のグループにおけるEntが特に高くなっているが、これらのカルスの色は他のカルスと異なり緑色が強く、このことがＢ成分のグレースケール画像におけるエントロピーの増加をもたらしているとも考えられる。また、R imageではＦ１とＣ１のグループの分離が明確になされなくなっているが、これはサンプルとなったＲＧＢカラー画像のカルスの色が白、黄、薄い緑であったために、Ｒ成分でスケール化することによりテクスチャ情報が失われたことが考えられる。このことから、単一成分のグレースケール画像だけでテクスチャ解析を行って算出されたエントロピー値の高低からフライアブルかコンパクトかを判別するのは難しいと言える。
【００２９】
そこで、異なる色成分画像についての相関を見ることを考える。図２（ａ）は表３〜６の結果において同じ色成分画像についての相関の例としてＣ１とＦ１とのグループにおけるΔEntGとEnt20Gとの関係を示しており，また、図２（ｂ）はΔEntBrightとEnt20Brightとの関係を示しているが、Ｆ１とＣ１とが明確に分離されてはいない。一方、図３（ａ）は異なる色成分画像についての相関の例としてＣ１とＦ１とのグループにおけるEnt20RとEnt20Gとの関係を示し、図３（ｂ）はEnt20BrightとEnt20Rとの関係を示しているが、これらの相関図においては、Ｆ１とＣ１とがより明確に分離されている。このことから、異なる色成分画像もとにした複数のグレースケール画像についてのエントロピーを用いることにより、カルスの色による影響を抑えてカルスの状態の判別を行うことができると考えられる。
【００３０】
また、論文の作者によるカルスの状態（フライアブル、コンパクト）について判断が画像を目視したことでは確認しにくいＦ２およびＣ２のグループにおけるカルスの画像を加えて同様にテクスチャ解析を行いエントロピー値を算出した。その結果から、図４（ａ）は同じ色成分画像の場合についての相関の例としてΔEntBとEnt20Bとの関係を示し、図４（ｂ）はEnt20BrightとEnt20Gとの関係を示している。これらの結果から、Ｆ２、Ｃ２を加えるとある程度フライアブル、コンパクトの判別がし難くなるものの、Ｆ１、Ｃ１のみのグループの画像の場合と同様に、異なる色成分画像についてのエントロピーの相関を見ることによりフライアブル、コンパクトがより明確に分離されることは示される。特にBrightness image−G imageのプロット、B image−G imageのプロットではほぼフライアブル、コンパクトの分離がなされることがわかる。
【００３１】
さらに、フライアブル、コンパクトの判別の精度を調べるために、最も明確な分析がなされたＦ１とＣ１とのグループにおけるB image−G imageのプロットについて統計ソフトウエア（ＭＩＮＩＴＡＢ）を用いて線形判別を行ったが、その結果、フライアブル・カルスの１０個の画像のうちフライアブルと判別できなかったものが１つ、コンパクト・カルスの１２個の画像のうちコンパクト・カルスと判別できなかつたものが０個であり、正確率で９５．５％であった。
【００３２】
このように、カルスのカラー画像をグレースケール変換した画像と、カラー画像の色成分画像をグレースケール変換した画像とについてそれぞれ濃度共起行列法を用いたテクスチャ解析を行い画像の特徴量としてのエントロピー値を算出し、同じ色成分間または異なる色成分間でのエントロピー値を対比させることによりフライアブル・カルスとコンパクト・カルスとを分離させることができる。多数の画像におけるカルスはこのようにしてフライアブルかコンパクトかに分離され、このように多数のカルスの例を表す画像についてテクスチャ解析により得られたエントロピー値を集積しておくことにより、コンパクト・カルスかフライアブル・カルスかのエントロピーの域値を設定しておくことにより、特定のカルスについフライアブルであるか、コンパクトであるかを定量的に判別できるものである。
【００３３】
エントロピー値の相関によりコンパクト・カルスとフライアブル・カルスとが明確に分離されるか否かについては、カルスの色が関係することが考えられる。カルスのカラー画像をグレースケール変換した画像と、カルスのカラー画像を色分解した画像とについてテクスチャ解析によりエントロピー値を得て、同じ色成分あるいは異なる色成分のグレースケール画像のエントロピー値の相関をとることにより、種々の色のカルスに対してカルスの状態の判別が良好になされるようになる。
〔２〕カルスの状態を判別するのに良好な画像を取得するための撮影の形態
［１］において、既存のカルスを撮影した画像について濃度共起行列法を用いたテクスチャ解析を行いカルスの画像についてのエントロピー値を算出することにより、カルスの状態がフライアブルであるかコンパクトであるかを定量的に判別できることを示したが、画像によっては、判別しにくい状況は起こり得る。これは、カルスを撮影する条件、特にカルス撮影の際に照明光をどのように当てるかということが関係することが大きいと考えられる。そのため、カルスを撮影した画像をもとにテクスチャ解析によりカルスの状態を判別する上で撮影条件がどのように影響するかについて知るために、実際にカルスを撮影して画像を取得し、その画像についてテクスチャ解析を行ってカルスが判別し易さがどのようになるかについて実験を行い検討した。
【００３４】
実験は、供試植物体を培地で培養し、カルス誘導を行って得られたカルスについて、培養条件において照明条件を変化させ撮影してカルスの画像を取得し、取得された画像について濃度共起行列法を用いたテクスチャ解析を行って評価するという形で行った。
（ａ）供試植物体
供試植物体としてのニンジン（品種名：金港四寸）を用い、ニンジンの種子を蒸留水で洗浄、表面殺菌し、滅菌水で濯いでから発芽促進媒体が入った試験管に１本あたり１０粒程度の滅菌種子を置床し、プラスチックの蓋をした後にアルミホィルで覆って遮光し、２５°Ｃのインキュベータ内で１４日間静置培養した。無菌播種から１４日目のニンジン小植物体の胚軸を無菌的に５ｍｍ程度切り出し、カルス培養誘導媒体が入った試験管１本につき５本ほどの胚軸切片を入れ、アルミホィルで覆って遮光し、２４°Ｃのインキュベータ内で焼く２か月間静置培養した。カルスができ始めた胚軸切片をカルス誘導媒体が入った三角フラスコに移植し、アルミホィルで蓋をして２７°Ｃ、８０ｒｐｍ、連続暗期の振盪培養機内で１０日間振盪培養した。胚軸から分離したカルスのみをカルス誘導媒体が入った三角フラスコに移植し、さらに２週間振盪培養し、細胞懸濁液を得た。
【００３５】
細胞懸濁液が入った試験管からピペットを用いて懸濁液をシャーレに移植し、パラフィルムで封をした後シャーレを２５°Ｃのインキュベータ内で約２か月間培養しコンパクト・カルスを得た。また、３０°Ｃのインキュベータ内で約２か月間培養しフライアブル・カルスを得た。
【００３６】
培養によって得られたニンジンのコンパクト・カルスもしくはフライアブル・カルスの入ったシャーレを水平な台上に載置し、シャーレに対し垂直上方から撮影するようにデジタルカメラを三脚に固定して数種の照明条件で撮影を行った。
【００３７】
図５（ａ），（ｂ）は培養環境にあるカルスを撮影する装置の配置形態を示す図であり、（ａ）が側面図、（ｂ）が平面図であり、小型暗箱１について断面として内部を示している。小型暗箱１内に培養環境でカルスＣを収容するシャーレ２が載置面上に配置されており、カルスＣを照明する光源として載置面上のシャーレ２の側方から照射光を与えるための白色ＬＥＤからなる側方照射光源３Ａと、シャーレ２の上側で斜め方向から照射光を与えるための白色ＬＥＤからなる斜め方向照射光源３Ｂとが配設されており、斜め方向照射光源３Ｂは支柱４に取り付けられたリング状取付部材５に必要な数だけ取り付けられており、また白色ＬＥＤは砲弾型等の指向性の高いものを用い、カルスＣに向かって照明光を照射するように取り付けられる。このように照明光源の配設することにより、側方照明光源３ＡではカルスＣを収容するシャーレ１の配置された面内の方向（０°）の光が照射され、斜め方向照射光源３Ｂではシャーレ１の配置された面に対して斜めの方向（例えば４５°）の光が照射される。また、側方照明光源３Ａと斜め方向照明光源３Ｂとの白色ＬＥＤはスイッチ手段を介して電源部（図示せず）に接続されており、側方照明光源３Ａと斜め方向照明光源３Ｂとのうちの一方を選択的に点灯して照明できるようにしてある。
【００３８】
小型暗箱１の外側に配設されたカメラ支持固定柱１０にカルスＣを収容したシャーレ１を上方から撮影できるようにダジタルカメラ１１が支持固定されている。デジタルカメラ１１のボデー側は操作のため小型暗箱１外とし、レンズ鏡胴は小型暗箱１内に入り込み、小型暗箱１の鏡胴が通り抜ける部分の壁部は抜き取られた形になっており、他の部分は遮光されている。また、図５（ｂ）において、カメラ支持固定柱１０は図示されておらず、デジタルカメラ１１の位置を示してある。
【００３９】
デジタルカメラはクールピクス９９０（ニコン製）で、レンズがマクロレンズ（Ｎｉｃｏｎ８×２０Ｄ）、撮像素子画素数が３００万画素ＣＣＤ、寸法が２４μｍ／画素のものであり、撮影距離は２５ｃｍとした。照明手段は白色発光ダイオード（ＬＥＤ）を用い、照射光強度（ＰＰＦＤ：５μｍｏｌ，５０μｍｏｌ）と照射角度（０°，４５°）を変えて撮影した。また、撮影により取得された画像は２４ビットＲＧＢカラーモード、２０４８×１５３６画素モードで保存した。
【００４０】
デジタルカメラによりカルスを撮影して得られた画像について前処理を行い、濃度共起行列法によるテクスチャ解析を行った。画像の前処理には画像処理ソフトウェア（フォトショップ６．０）を用い、撮影した画像の画素数をバイキュービック補間方式に基づく画像の変換機能により修正し、画像を５１２×３８４画素にした。その画像の背景となるカルス以外の部分を自動変換ツールで選択して取り除きカルスの部分のみの画像としてから、そのＲＧＢカラー画像を〔１〕で説明したように４種類の８ビットのグレースケール画像に変換した。このようにコンパクト・カルスの画像およびフライアブル・カルスの画像のそれぞれについて４種類の８ビットのグレースケール画像を生成し、それらの各々について濃度共起行列法を用いたテクスチャ解析を行ってエントロピー値を算出した。
【００４１】
４つのコンパクト・カルスの画像Ａ，Ｂ，Ｃ，Ｄから生成された各色成分画像（B image，G image，R image，Brightness image）について算出されたエントロピー値の結果を表７〜１０に示しており、照明光の照射条件はそれぞれ表７が４５°，５０μｍｏｌ、表８が０°，５０μｍｏｌ、表９が４５°，５μｍｏｌ、表１０が０°，５μｍｏｌである。また、４つのフライアブル・カルスの画像Ｅ，Ｆ，Ｇ，Ｈから生成された各色成分画像（B image，G image，R image，Brightness image）について算出されたエントロピー値の結果を表１１〜１４に示しており、照明光の照射条件はそれぞれ表１１が４５°，５０μｍｏｌ、表１２が０°，５０μｍｏｌ、表１３が４５°，５μｍｏｌ、表１４が０°，５μｍｏｌである。
【００４２】
【表７】
2008261631.tif 000010
【００４３】
【表８】
2008261631.tif 000011
【００４４】
【表９】
2008261631.tif 000012
【００４５】
【表１０】
2008261631.tif 000013
【００４６】
【表１１】
2008261631.tif 000014
【００４７】
【表１２】
2008261631.tif 000015
【００４８】
【表１３】
2008261631.tif 000016
【００４９】
【表１４】
2008261631.tif 000017
図６〜９は、表７〜１０における色成分ごとのエントロピー値について、同じ色成分の間でのエントロピー値の相関あるいは異なる色成分の間でのエントロピー値の相関をとったものの例を示している。図６（ａ），（ｂ）は照明光の照射条件４５°、５０μｍｏｌでのそれぞれＢ成分およびＲ成分のグレースケール変換した画像のエントロピー値の相関を示し、図７（ａ）（ｂ）は照明条件０°、５０μｍｏｌでのそれぞれＢ成分およびＲ成分のグレースケール変換した画像のエントロピー値の相関を示し、図８（ａ），（ｂ）は照明光の照射条件４５°、５０ｍｏｌでのエントロピー値と照射条件０°５０μｍｏｌでのエントロピー値との相関をそれぞれＢ成分のグレースケール変換した画像とＲ成分のグレースケール変換した画像の場合について示している。図９（ａ），（ｂ）は照明光の照射条件４５°、５μｍｏｌでのそれぞれＧ成分およびカラー画像からのグレースケール変換した画像のエントロピー値の相関を示している。
【００５０】
表７〜１０における結果では、照明光の照射強度を高くすると、エントロピーが高い値になること、照射強度が低いとコンパクト・カルスとフライアブル・カルスのエントロピーの違いがわかりにくくなり、R imageではコンパクト・カルスのエントロピーがフライアブル・カルスより低い値になっていることがわかる。エントロピー値の相関で見た場合に、図６〜８のように５０μｍｏｌの照射光ではフライアブルとコンパクトとが良好に分離されるが、図９のように５ｍｏｌの照射光の場合にはフライアブルとコンパクトとが明確に分離されないことがわかる。このように、照明光の照射強度が十分でないと、コンパクト・カルスとフライアブル・カルスとにおけるエントロピー値の差が明確にならず、エントロピー値の相関を示す図においてコンパクト・カルスとフライアブル・カルスとが明確に分離されないことがわかる。
【００５１】
また、ＰＰＦＤ５０μｍｏｌｍ−２ｓ−１で照射角度０°の時は、コンパクト・カルスの画像とフライアブル・カルスの画像とでエントロピー値に大きな違いが見られないが、ＰＰＦＤ５０μｍｏｌｍ−２ｓ−１で照射角度４５°の時には、コンパクト・カルスの画像とフライアブル・カルスの画像とでエントロピー値に大きな差が生じた。これは、照射角度が０°の時に撮影された画像においてカルス自体の凹凸による影が写ることによりフライアブル・カルスにおいてもウントロピーが増大するためと考えられ、さらに、照射角度０°では反射光だけでなくカルスを透過した光も画像に入り込み、光の透過率が比較的高いと考えられるフライアブル・カルスの画像におけるコントラストが高まり、結果的に画像における乱雑さがましてエントロピー値が高くなることも考えられる。一方、照射角度４５°で撮影すると、カルス自体の影が写り込むこともなく、また、透過光の影響も抑えられるので、コンパクト・カルスの画像とフライアブル・カルスの画像とにおけるエントロピー値の差が大きくなると考えられる。
【００５２】
このように、十分な照射光強度で撮影されたフライアブル・カルスの画像では、照明光の照射角度０°の場合と照射角度４５°の場合とでエントロピー値に大きな違いがあることがわかる。一方、コンパクト・カルスの画像については、照射光の角度における有意な差は得られていない。つまり、判別対象となるカルスがフライアブルであればあるほど光照射角度による影響を受けやすく、逆にコンパクト・カルスであればあるほど光照射角度の影響を受けにくいということがわかる。
【００５３】
なお、実験結果としては特に照射角度が０°と４５°との場合について示しているが、照射角度０°はシャーレの載置されている面内における側方からの照明であり、一方、照射角度４５°はシャーレの載置されている面に対して角度をなす斜めの方向での照明についての代表的な照射角度であって、照射角度については例えば２０°〜５０°の範囲内の他の角度においても、エントロピー値として若干変わってくるが同様な傾向が見られる。
【００５４】
また、カルスの撮影において、台上に載置されたシャーレの上方から撮影しており、カルスは培養環境にあるので、シャーレ内面側、特に上側の部分に水滴が付着していて、照射角度４５°のような斜め方向の照射の場合に、この水滴の影響を受けるため、前述した実験の際にクリーンベンチ内でシャール内の水滴を取り払っている。この状況を避けるためには、シャーレ内のカルスを下側から撮影する形態をとるようにすればよい。この場合、斜め方向の照明光源は図５の場合と対称的に斜め下方からとすることになる。ただし、培地の光透過性についての評価が必要になる。
【００５５】
カルスの撮影のための照明光源として白色ＬＥＤを用いた例を示したが、他の白点状または線状の白色光源を用いてもよい。
【００５６】
以上のように、本発明によるカルスの画像について濃度共起行列法を用いたテクスチャ解析を行いカルスの状態がコンパクトであるか、フライアブルであるかを判別する上で、照明光源を十分な明るさのものとして、台上に載置されたシャーレ内のカルスを側方から照射して撮影したカルスの画像と、斜め方向から照射して撮影したカルスの画像との両方についてテクスチャ解析を行いエントロピー値の相関をとるのがカルスの状態の判別する上で有効であると言える。
【００５７】
濃度共起行列法を用いたテクスチャ解析によりカルスの状態を判別するのに良好な画像を取得する方法について説明したが、このように画像を撮影し取得されたカルスの画像についてテクスチャ解析処理を行う装置の基本的形態について図１０に示す。図１０において、１２は培養環境にあるカルスを収容する透明な培養器（シャーレ）の載置台であり、３Ａ，３ＢはそれぞれＬＥＤ等の側方照明手段、斜め方向照明手段であり、１１は撮像素子を有する撮影部であり、撮影部１１と載置台１２との間の空間は側方照明手段３Ａおよび斜め方向照明手段３Ｂを含んで遮光部材（図示せず）で外部から遮光される状態になっている。側方照明手段３Ａは培養器内のカルスの高さ位置に合わせて位置調整可能に支持されており、また、斜め方向照明手段３Ｂは培養器内のカルスに対する照射角度を調整できるように位置調整可能に支持されている。撮影動作制御部２１は側方照明手段３Ａ，３Ｂを位置調整し、選択的に動作させるとともに、撮影部での撮影動作の制御を行う。
【００５８】
撮影部１１で撮影された画像のデータは画像取得前処理部２２に移送され、ビット数変換、背景除去、グレースケール変換、色分解・グレースケール変換の処理がなされる。これらの処理がなされた画像データは逐次記憶部２３に蓄積保持される。テクスチャ解析処理部２４においては、前処理のなされたカルスのグレースケール変換画像について濃度共起行列法を用いたテクスチャ解析を行い、各グレースケール変換画像のエントロピー値を算出し、相関対比する演算を行い、演算処理の結果を記憶部２３に蓄積保持する。表示部２５は、撮影動作、演算処理、処理結果のデータの表示を行い、中央制御部２６は全体的な制御を行う。
【００５９】
図１０において、画像動作制御部２１、画像取得前処理部２２、テクスチャ解析処理部２３、中央制御部２６は各々の制御ないし演算の処理を行う専用処理回路を備えるものとして構成され、本発明はこのような専用処理回路を備える植物培養細胞塊の状態を判別するための装置の形態とすることができる。
【００６０】
また、画像取得、前処理、テクスチャ解析処理における演算処理を行うプログラムを備えた汎用コンピュータにより濃度共起行列法を用いたカルスの状態の判別を行うことができ、本発明はこのような濃度共起行列法を用いたカルスの状態の判別を行うことをコンピュータ上で実行するためのプログラムとしての形態をも有するものである。
【図面の簡単な説明】
【００６１】
【図１】濃度共起行列の例を示す図である。
【図２】既存のカルスのカラー画像の同じ色成分のグレースケール画像についてのエントロピー値の相関を示す図である。
【図３】既存のカルスのカラー画像の異なる色成分のグレースケール画像についてのエントロピー値の相関を示す図である。
【図４】図２および図３の場合と異なる種類のカルスの画像を含めた場合のエントロピー値の相関を示す図である。
【図５】培養環境においてカルスを撮影するための装置の配置形態を示す図である。
【図６】十分な照明光におけるエントロピー値の相関を示す図であり、（ａ）Ｂ成分のグレースケール変換した画像の場合（ｂ）Ｒ成分のグレースケール変換した画像の場合について示すものである。
【図７】十分な照明光におけるエントロピー値の相関を示す図であり、（ａ）Ｂ成分のグレースケール変換した画像の場合（ｂ）Ｒ成分のグレースケール変換した画像の場合
【図８】十分な照射光におけるエントロピー値の相関を示す図であり、（ａ）照射角４５°でのエントロピー値と照射角０°でのエントロピー値との相関を示す図であり、（ａ）Ｂ成分のグレースケール変換した画像の場合（ｂ）Ｒ成分のグレースケール変換した画像の場合について示すものである。
【図９】輝度の低い照明光におけるエントロピー値の相関を示す図であり、（ａ）Ｇ成分のグレースケール変換した画像の場合（ｂ）カラー画像のグレースケール変換した画像の場合について示すものである。
【図１０】テクスチャ解析処理を行う装置の基本的形態を示すブロック図である。
【符号の説明】
【００６２】
１    小型暗箱
２    シャーレ
３Ａ  白色ＬＥＤ（側方照明手段）
３Ｂ  白色ＬＥＤ（斜め方向照明手段）
４    支柱
５    リング状取付部材
１０  カメラ支持固定柱
１１  デジタルカメラ（撮影部）
１２  培養器載置台
２１  撮影動作
２２  画像取得前処理部
２３  記憶部
２４  テクスチャ解析処理部
２５  表示部
２６  中央制御部
【図１】
2008261631.tif 000018
【図２】
2008261631.tif 000019
【図３】
2008261631.tif 000020
【図４】
2008261631.tif 000021
【図５】
2008261631.tif 000022
【図６】
2008261631.tif 000023
【図７】
2008261631.tif 000024
【図８】
2008261631.tif 000025
【図９】
2008261631.tif 000026
【図１０】
2008261631.tif 000027