JP6483826B2 - 画像解析装置、画像解析方法、画像解析システム、画像解析プログラム、および記録媒体 - Google Patents

Description

本発明は、組織を撮像した画像を解析して、組織の構成を判定するための情報を出力する画像解析装置および画像解析方法に関する。

病理診断は治療方針を決定するための重要な情報を提供する。病理診断において、診断の誤りの発生を回避するためには、二重チェックシステムを適用することが好ましいとされている。この二重チェックシステムとは、病理医Ａによる診断に加え、病理医Ｂによる診断も行い、両者の所見に基づいて病理診断の結果を決定するシステムである。しかし、現状では、一部の先進国を除いて病理診断を行う病理医の絶対数は不足している。また、病理医は大都市部などの病院に偏在しており、常勤の病理医を持たない病院も多い。それゆえ、複数の病理医を常に確保して迅速に病理診断を実施することは、日本のみならず世界的にも容易ではない。その上、日本において病理医１人が１年間に行う病理診断において観察する組織スライドの枚数は平均約１万枚ともいわれており、各病理医の負担も大きくなっている。このような現状に鑑み、病理診断において病理医が担う役割を計算機技術が代行することによって、病理診断を支援する技術の開発が望まれている。

病理医は、病理診断において顕微鏡観察による組織・細胞の形態学的所見に基づいた客観的な判断が求められる。しかし、形態学的所見のみからでは正しく分類し、峻別することが困難な病変もあり得るため、病変の看過など病理診断の誤りを生じる虞がある。また、病理医の経験に基づいた主観的な判断に基づいた病理診断の場合、病理医ごとに診断結果が異なる可能性があり、病理診断の信頼性を損ねることに繋がる。診断ミスを生じさせることが無く、かつ客観的な診断基準を提供できることが、病理診断を支援する技術には要求される。

病理診断に用いられる組織スライドなどの画像を、計算機技術を用いて解析するために、位相幾何学的な考え方の導入が提案されている。病理診断を支援する技術の例として、特許文献１には、病理画像に基づいて癌病変部を抽出する画像解析装置が開示されている。特許文献１に記載の画像解析装置は、病理画像における単位面積当たりのホモロジーを領域ごとに算出することにより、各領域が注目領域（例えば、癌病変部を含む領域）であるか否かを判断する。

また、非特許文献１には、肝臓を撮像した画像を、パーシステントホモロジー（persistent homology）の概念を適用して解析し、肝臓の病変を分類する方法について記載されている。

ＷＯ２０１０／０８７１１２号（２０１０年 ８月 ５日公開）

Aaron Adcock et al.，「Classification of Hepatic Lesions using the Matching Metric」，Computer Vision and Image Understanding, Vol. 121, p.36-42，2014.

組織を構成する構成要素間の関係を適切に評価することは、当該組織の物性および構造を理解し、その組織の構成を判定する上で重要である。

例えば、生体組織の場合、生体組織に生じた癌の診断に有効な画像解析の場合を例に挙げれば、ある組織に生じた癌細胞が異常増殖するに伴い、細胞の形状は歪になり、組織が有している本来の構成は失われ、癌細胞が無秩序に散在する状態へと推移する。病理診断では、このような組織・細胞の形態に関する病理医の所見に基づいて癌の「分化度」が規定されている。

国際疾病分類、腫瘍学、第３版（ＩＣＤ-Ｏ-３）では、大腸癌はその分化度に応じて「（高・中・低）分化癌」と表記されている場合、「腺癌」として扱われ、形態コード「８１４０／３」が採用される。そして癌（腺腫）の「分化度」は、例えば腺癌の場合、８１４０／３の後ろに分化度を示す数字１〜３を付して規定される。なお、従来同様、正常→腺癌（高分化）→腺癌（中分化）→腺癌（低分化）の順に、組織の細胞が有している特徴が失われた状態を示している。
腺癌（高分化）：８１４０／３１
腺癌（中分化）：８１４０／３２
腺癌（低分化）：８１４０／３３。

病理診断には、癌の分化度の判別が含まれる。組織の物性や構造を理解し、その構成を判定することが重要であることは、生体組織の癌を判定する場合に限らず、さまざまな組織においても同様である。

本発明の目的は、上記の問題点を解決するためになされたもので、その目的は、組織を構成する構成要素間の関係を適切に評価するための情報を生成することができる画像解析装置等を提供することである。

本発明に係る画像解析装置は、上記の課題を解決するために、組織を撮像した１つの撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化部と、上記二値化部が生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出部と、上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定部と、上記特定部が特定した上記特徴数に対応した出力情報を生成する出力情報生成部とを備えている。

また、本発明に係る画像解析方法は、上記の課題を解決するために、組織を撮像した撮像画像を解析する画像解析装置における画像解析方法であって、１つの上記撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化工程と、上記二値化工程にて生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出工程と、上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定工程と、上記特定工程にて特定した上記特徴数に対応した出力情報を生成する出力情報生成工程とを含んでいる。

本発明は、組織の画像を解析することにより、当該組織の構成要素間の関係が異なることを客観的かつ適切に評価できるという効果を奏する。

本発明の実施例１に係る画像解析装置の構成の一例を示すブロック図である。 本発明に係る画像解析装置に適用可能な１次元ベッチ数について説明する図である。 本発明の実施例１に係る画像解析装置における処理の流れの一例を示すフローチャートである。 本発明の実施例１に係る画像解析装置における二値化画像を生成する処理の流れの一例を示すフローチャートである。 （ａ）は元の組織画像（Ｂの画素値の最大値は１７１）であり、（ｂ）、（ｃ）、（ｄ）はＢの画素値の最大値にそれぞれ１、０．７、０．５を乗算した値を用いて二値化された組織画像の一例を示す図である。 画像解析装置によって解析した結果の一例を示す図である。 画像解析装置によって解析した結果の一例を示す図である。 画像解析装置によって解析した結果の一例を示す図である。 近似式の二次の係数ａの正/負の割合を示すグラフである。 （ａ）は近似式として二次関数の式を用いる場合の近似式の一例を示し、（ｂ）は近似式の二次の係数の対数と、近似式の軸の位置を示すｂとを座標として有する点をプロットした分布図である。 判定用画像の一例を示す図である。 本発明の実施例２に係る画像解析装置の構成の一例を示すブロック図である。 本発明に係る画像解析装置を含む画像解析システムの構成例を示す概略図である。 （ａ）は酸性条件でゲル化したシリコーンゲルの分子構造を示す図であり、（ｂ）は酸性条件でゲル化したシリコーンゲルを、透過型電子顕微鏡を用いて撮像した画像であり、（ｃ）は塩基性条件でゲル化したシリコーンゲルの分子構造を示す図であり、（ｄ）は塩基性条件でゲル化したシリコーンゲルを、透過型電子顕微鏡を用いて撮像した画像である。 （ａ）〜（ｇ）は、図１４の（ｂ）の画像を、二値化の基準値を異ならせながら二値化した二値化画像である。 （ａ）〜（ｇ）は、図１４の（ｄ）の画像を、二値化の基準値を異ならせながら二値化した二値化画像である。 シリコーンゲルの構成を解析した結果の一例を示す図である。

以下、本発明の実施の一形態について、図１〜１０を用いて説明すれば、以下の通りである。

（本発明の技術思想）
まず、本発明の技術思想について以下に説明する。本発明の目的は、組織を撮像した撮像画像を解析することにより、当該撮像画像に写っている癌の分化度を判別するために有用な情報を出力することである。

本発明の発明者は、組織標本を撮像した撮像画像を調べ、癌化に伴う癌細胞の大きさ、形状、または各細胞の位置関係を詳細に比較・検討した。その結果、癌を含む組織が写っている組織画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行い、変化させた二値化の基準値と、当該基準値を用いて二値化された画像に含まれる組織の像が示す１次元ベッチ数との関係に基づいて、癌の分化度を判別することが可能となることを見出した。

なお、本発明に係る画像解析装置１は、生体組織、鉱物組織、および樹脂製品の組織など、さまざまな種類の組織の構成の分析に適用することができる。「組織」とは、複数種類の構成要素が、所定のパターンで集合した構造単位を意味している。本発明の解析対象となる組織は、細胞などの集まりから構成される生体組織、結晶粒などの集まりから構成される鉱物組織、およびシリコーンなどの樹脂材料の集まりから構成される樹脂製品の組織などであってもよい。

以下では、生体組織に生じた癌の分化度を判別するために有用な情報を出力する画像解析装置１を例に挙げて説明するが、本発明の解析対象は、生体組織に限定されない。

＜癌の分化度＞
ある組織に生じた癌細胞が異常増殖するに伴い、細胞の形状は歪になり、組織が有している本来の構成は失われ、癌細胞が無秩序に散在する状態へと推移する。病理診断では、このような組織・細胞の形態学的所見に基づいて癌の「分化度」が規定されている。

生体組織は構成要素として、特定の形態および機能を有する細胞などを有している。組織は構成要素が互いに連携することによって、当該組織が有する機能を発揮する。それゆえ、正常部の構成要素は、組織の機能を発揮するために、互いに強い結びつきを有している。

一方、組織が腫瘍化・癌化している部位では、不純物の存在が増加したり、細胞が異常増加したりするために、構成要素間の関係の強さが正常部より弱くなる。さらに、組織の癌化が進めば、組織の構成が失われた、いわゆる脱分化した状態（低分化、中分化）となる。

＜構成要素間の関係を定量化するための数学的表現＞
組織の構成要素間の関係を定量化するために、本発明の発明者らは、ホモロジーの概念、特にパーシステントホモロジーの適用を試みた。ホモロジーとは、図形の形態上の性質を代数的に置き換えて、図形の結合などの解析を容易にする数学の一分野である。本発明の発明者らは、ホモロジーの概念の中でも１次元ベッチ数の利用に着目した。

構成要素の接触を表す数学的な概念であるホモロジーの概念は、病理画像において、適切な二値化の基準値（二値化パラメータともいう）を設定すれば、二値化された画像から単位面積当たりの１次元ベッチ数を算出することにより、その組織の構成要素同士の接触の程度の評価が可能である。

ベッチ数とは、図形（構成要素）の形状には無関係であり、図形の接触と分離とにのみ関係するトポロジカルな示唆数である、ｑ次元特異ホモロジー群が有限生成のとき、このｑ次特異ホモロジー群は、自由アーベル群と有限アーベル群との直和に分けられる。この自由アーベル群の階級をベッチ数という。２次元の場合、０次元ベッチ数は連結成分の数であり、１次元ベッチ数は当該連結成分が外縁となって囲んだ空間（穴形状の領域）の数、すなわち連結成分中に存在する「穴」の数である。

＜１次元ベッチ数＞
１次元ベッチ数は、数学的には以下のように定義される。一般に有限個の線分を繋ぎ合わせて成る図形（１次元複体）Ｋに対して、適当なｒ個の、両端を含まない１次元単体（線分）を図形Ｋから取り去っても図形Ｋの連結成分の個数は増加しないが、任意の（ｒ＋１）個の両端を含まない１次元単体（線分）をＫから取り去った場合にはＫの連結成分の個数が１つ増加するとき、図形Ｋの１次元ベッチ数はｒである。

図２に示す図形を用いて、１次元ベッチ数について説明する。図２は、本発明に係る画像解析装置１に適用可能なベッチ数について説明する図である。図２に示す図は、環Ｃ１と環Ｃ２が互いに交わっている。この図形が有している「穴」の数は３つである。一方、この図形の１次元ベッチ数は３と算出される。すなわち、ある図形について１次元ベッチ数を算出するということは、当該図形が有している「穴」の数を数えていることに等しい。

〔実施形態１〕
（画像解析装置１の構成）
ここでは、画像解析装置１の構成について図１を用いて説明する。図１は、本発明の実施例１に係る画像解析装置１の構成の一例を示すブロック図である。図示のように、画像解析装置１は、画像取得部２（受信部）、記憶部３、画像解析部４（画像解析装置）、表示制御部５、および一次記憶部７を備えている。なお、画像解析装置１が別体として設けられた表示装置６（提示部）と接続されている例を図示しているが、これに限定されない。例えば、画像解析装置１が表示装置６（提示部）を内蔵する構成であってもよいし、複数の表示装置６と無線／有線で接続された構成であってもよい。

画像取得部２は、外部機器８（例えば、顕微鏡に接続された撮像装置、画像データを保存可能な記憶装置など）から組織を撮像した撮像画像（以下、組織画像と称する）を取得し、取得した組織画像を記憶部３に格納する。上記組織画像は、組織を適当な倍率で撮像した画像である。なお、組織画像の撮像倍率は、解析の対象となる組織に応じて、当業者によって適宜設定されればよい。また、組織画像は、組織をそのまま撮像した画像であってもよいし、適当な染色を施した後の組織を撮像した画像であってもよい。

組織を染色する方法は、特に限定されず、例えば、生体組織の場合には、ＨＥ（Hematoxylin-Eosin）染色などを用いることができる。ＨＥ染色された組織では、細胞核および細胞質が染色されるため、組織を構成する細胞間の関係を観察することに適している。

記憶部３は、画像取得部２が取得した組織画像とともに、画像解析部４が実行する（１）各部の制御プログラム、（２）ＯＳプログラム、（３）アプリケーションプログラム、および、（４）これらプログラムを実行するときに読み出す各種データを記録するものである。記憶部３は、ハードディスク、フラッシュメモリなどの不揮発性の記憶装置によって構成される。

一次記憶部７は、上述の各種プログラムを実行する過程でデータを一時的に保持するための作業領域として使用されるものであり、ＲＡＭ（Random Access Memory）などの揮発性の記憶装置によって構成される。

表示制御部５は、表示装置６を制御することにより、画像解析部４から出力される解析結果を示す画像（出力情報）等を表示装置６に表示する。

表示装置６は、画像解析部４から出力される出力情報等を表示する表示装置であり、例えば液晶ディスプレイである。なお、画像解析装置１が専用の表示装置６を備える構成であってもよい。なお、表示装置６の表示画面にタッチセンサが重畳されており、ユーザによる表示面へのタッチ操作を表示装置６が検出できる構成であってもよい。

（画像解析部４の構成）
画像解析部４は、画像取得部２が取得した組織画像を解析することにより、当該組織画像に含まれる癌組織の分化度を示す情報を出力する。この画像解析部４は、二値化部４１（分割部）、ベッチ数算出部４２（空間数算出部）、関数特定部４３（特定部）、および出力情報生成部４４を備えている。

二値化部４１は、組織画像の二値化処理において、二値化の基準値よりも大きい画素値を有する画素を白色の画素に変換し二値化の基準値以下の画素値を有する画素を黒色の画素に変換する。このとき二値化部４１は、組織を撮像した１つの組織画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する。

二値化の基準値の決定方法の一例を説明すれば以下の通りである。二値化部４１は、組織画像に含まれる所定の色成分の画素値についての分布データ（ヒストグラム）に基づいて当該色成分の最も頻度が高い（最も多く含まれている）画素値を標準画素値として特定する。そして、該標準画素値に所定の範囲の倍率（例えば、０．５５〜０．７２）を乗算した値を、二値化の基準値として決定する。

二値化部４１は、組織画像のグレースケール画像に対して、上述のように決定した複数の二値化の基準値を用いて二値化処理を行い、複数の二値化画像を生成する。

なお、組織画像を二値化するための基準値は、組織画像に含まれる構成要素の外縁を検出して、構成要素間の関係を判断する上で重要である。二値化の基準値より大きい画素値を有する画素は白になり、二値化の基準値より小さい画素値を有する画素は黒になるため、二値化の基準値を大きく設定し過ぎると、二値化された組織画像は境界線が過度に強調され、黒っぽい画像になる。一方、二値化の基準値を小さく設定し過ぎると、二値化された組織画像は、境界線の少ない白っぽい画像になる（図５参照）。組織画像に対して行う二値化処理については、後に詳述する。

また、二値化部４１が、画像取得部２が取得した組織画像を所定の大きさの分割領域（以下において、「分割画像」と表す場合もある）に分割する画像分割部（分割部）としての機能を備えてもよい。あるいは、画像解析部４が、二値化部４１、ベッチ数算出部４２、関数特定部４３、および出力情報生成部４４に加え、画像分割部（図示せず）をさらに備える構成であってもよい。二値化部４１は、分割領域ごとに二値化の基準値を変化させる範囲を決定して二値化することも可能である。しかし、各分割領域を比較して、組織画像に写っている注目領域を判別する場合には、二値化の基準値を変化させる範囲を共通にしてすべての分割領域の二値化処理を行うようにすることが望ましい。

ベッチ数算出部４２は、二値化部４１が生成した複数の二値化画像について、当該二値化画像に含まれる、閉じた外縁をそれぞれ有する複数の構成成分の外縁によって囲まれた穴の数（穴形状の領域の数）を算出する。すなわち、ベッチ数算出部４２は、複数の二値化画像にそれぞれ含まれる組織の像について、１次元ベッチ数を算出する。

上記穴は、１つ以上の構成成分の外縁の少なくとも一部（構成成分が１つの場合には、その構成成分の外縁の全部）をその外縁として有している開口部である。換言すれば、構成成分の外縁によって囲まれた穴は、（１）１つの構成成分がその内部に有する穴と、（２）互いに連結した複数の構成成分がそれぞれ有する外縁によって囲まれた穴とに分類される。

ベッチ数算出部４２として、既存のプログラムを用いてもよい。このようなプログラムの一例としてＣＨｏｍＰを挙げることができる。ＣＨｏｍＰは、ＧＮＵ（General PublicLicense）に準拠したフリーウェアであり、ホームページ（http://chomp.rutgers.edu/）からダウンロードできる。画像に関する１次元ベッチ数を算出できるプログラムであれば、ＣＨｏｍＰ以外のプログラムを用いてもよい。

関数特定部４３は、複数の二値化の基準値と穴の数との対応関係を特徴づける数値（特徴数）を特定する。この数値は少なくとも１つ特定されればよい。例えば、関数特定部４３は、複数の二値化の基準値と穴の数との関係を近似的に示す関数（近似関数）を表す式を導出して、該導出した式に含まれる係数の少なくとも１つを特徴数として特定してもよい。この場合、関数特定部４３は、二値化の基準値の変化に伴う穴の数の変化を近似的に示す関数の特徴を表す１または複数の数値（特徴数）を算出（特定）する。ここで、二値化の基準値の変化に伴う穴の数の変化を近似的に示す関数（近似関数）としては、任意の関係式で表記可能な関数であってもよく、例えば、ｎ次関数（ｎ≧１）（線形近似および対数近似を含む）、指数関数（累乗近似）、および三角関数などであってもよい。以下では、関数特定部４３は、二値化の基準値の変化に伴う穴の数の変化を近似的に示すｎ次関数（ｎ≧１）の式（ｎ次式）に含まれる、ｎ次の係数（特徴数）を含む複数の係数を算出する場合を例に挙げて説明する。

なお、関数特定部４３は、複数の二値化の基準値と穴の数との対応関係を特徴づける複数の特徴数を、近似関数を導出することなく特定するように構成されていてもよい。例えば、複数の二値化の基準値と穴の数との対応関係を特徴づける特徴数を特定するために、関数特定部４３が深層学習（ディープラーニング）などの手法を適用する構成であってもよい。例えば、関数特定部４３は、所定の特徴（癌の分化度など）を有する典型的な画像（組織画像等）と、当該画像に対して複数の二値化の基準値を用いて二値化処理を行ったときの二値化の基準値と上記二値化処理された画像が示す穴の数との対応関係に対応する特徴数との複数の組み合わせを学習データとして用いた深層学習を施されていてもよい。ここで、学習データとして用いることが可能な特徴数は、上記二値化の基準値と穴の数との対応関係に対応する数であればよく、例えば、癌（組織）の分化度を示す任意の数値（例えば、高分化度の場合は１、中分化度の場合は２、低分化度の場合は３など）であってもよい。このような深層学習を採用した場合、関数特定部４３は、二値化の基準値と穴の数との関係を近似的に示す関数を示す式を導出することなく、適切な特徴数を特定することができる。

出力情報生成部４４は、関数特定部４３が算出した複数の特徴数（例えば、近似関数の係数）のうちの少なくとも１つに対応した出力情報を生成する。上記出力情報としては、例えば、関数特定部４３が算出した複数の係数のそれぞれに１つずつが対応している複数の座標軸によって形成される座標空間における、関数特定部４３が算出した複数の係数を座標として有する点（座標点と称する）の位置を反映した情報（出力情報）が挙げられる。また、出力情報生成部４４は、出力情報として、分割領域が示された判定用画像を生成し、撮像画像における分割領域についての座標点の位置に応じて、判定用画像における分割領域の提示態様を異ならせてもよい（図１１の（ｂ）参照）。例えば、低分化の癌を含む分割領域を赤色、腺腫を含む分割領域をオレンジ色、高分化の癌を含む分割領域を黄色で示すなど、癌を含まない正常な分割領域との差を明確にして表示してもよい。

なお、出力情報生成部４４は、出力情報として、（１）上記座標点を示す複数の数値の組み合わせ、（２）上記座標点が示されたグラフ（例えば、図１０の（ｂ）に示す分布図）、（３）上記座標点の位置に基づいて判定された癌の分化度を示す情報（文字、数値、色情報など）を生成してもよい。分割領域が複数存在する場合には、これらの出力情報に、解析対象となった分割領域を特定する情報を含めることが好ましい。

（画像解析装置１の処理の流れ）
次に、画像解析装置１における処理の流れの一例について図３を用いて説明する。図３は、画像解析装置１における処理の流れの一例を示すフローチャートである。ここでは、組織画像は既に記憶部３に格納されているものとして説明する。また、ここでは、組織画像を分割して分割領域を生成する処理を含む例を示す。

まず、二値化部４１は、記憶部３から組織画像を取得し、取得した組織画像を所定の大きさの領域に分割することにより分割領域を生成する（Ｓ１）。分割領域は、元の組織画像を縦３×横３、縦５×横５などに分割した場合の各領域である。

続いて、二値化部４１は、二値化の基準値を変化させながら二値化処理を行い、１つの組織画像から複数の二値化画像を生成する（Ｓ２：二値化工程）。なお、二値化部４１が組織画像から二値化画像を生成する処理については、後に説明する。

画像データは画素（ピクセル:pixel）と呼ばれる微小な正方形を多数用いて、その正方形上の光強度を決めることで各種図形を表現している（ビットマップ形式）。二値化するとは、この光強度がある値より大きければ黒（計算機では１と認識）、ある値以下であるならば白（計算機では０と認識）することである。

二値化部４１は、生成した二値化画像をベッチ数算出部４２に送る。

ベッチ数算出部４２は、受け取った二値化画像が示す分割領域ごとに１次元ベッチ数を算出する（Ｓ３：空間数算出工程）。ベッチ数算出部４２は、算出した１次元ベッチ数を関数特定部４３へ出力する。

関数特定部４３は、ベッチ数算出部４２から出力された１次元ベッチ数と二値化の基準値との関係を示す関係式に含まれる複数の係数を算出する（Ｓ４：関数特定工程）。関数特定部４３は、算出した複数の係数を、分割領域を特定する情報と対応付けて出力情報生成部４４へ出力する。

出力情報生成部４４は、分割領域ごとに、関数特定部４３が算出した複数の係数を座標として有する座標点の座標空間における位置を特定し、特定した位置に対応した出力情報（例えば、上述の判定用画像）を生成する（Ｓ５：出力情報生成工程）。出力情報生成部４４は、生成した出力情報を表示制御部５へ出力する。

表示制御部５は、出力情報生成部４４から出力された出力情報を表示装置６に表示させる（Ｓ６）。

（二値化画像を生成する処理の流れ）
次に、二値化部４１が組織画像から二値化画像を生成する処理について、図４を用いて説明する。図４は、本発明の実施例１に係る画像解析装置１における二値化画像を生成する処理の流れの一例を示すフローチャートである。

二値化部４１は、組織画像の画像データに含まれる画素の色成分（例えば、ＲＧＢ）ごとの画素値を取得し、二値化の基準値を設定するための標準画素値を決定する（Ｓ１１）。例えば、組織画像がＨＥ染色された組織の切片を撮像した画像である場合、ＲＧＢの各色成分のうち青（Ｂ）の画素値に基づいて標準画素値を決定する。この標準画素値の決定方法の具体例について以下に説明する。

画素値ＨＥ染色されている組織の場合、細胞核が青く染まる。核に着目して組織画像を画像解析する場合にはＢの画素値に基づいて二値化の基準値を決定することが適している。組織染色を行った場合には、染色ムラが発生することがあるため、標準的に染色された部分を基準として二値化の基準値を決定することが好ましい。そのため、組織画像を構成するＢの画素値のうち、最も高い頻度で含まれているＢの画素値を標準画素値として決定する。このように二値化の基準値を決定することにより、染色の度合いが異なる複数の組織についての複数の撮像画像を用いる場合に、染色度合いを標準化することができる。

具体的には、組織画像における、Ｂの画素値と、当該画素値を有する画素の個数（画素数）との関係を示す分布データにおいて、上記画素数が極大となる画素値である標準画素数を特定する。すなわち、標準画素数とは、組織画像を構成するＢの画素値のうち、最も高い頻度で含まれている画素値である。なお、出現頻度が極大となるＢの画素値が複数存在する場合には、極大値が大きい方の画素値を選択すればよい。

ここでは、Ｂの画素値に基づいて二値化の基準値を決定する場合について説明するが、組織画像の色成分の分布や染色に用いた色素の色調によっては、ＲＧＢの各色成分のうち赤（Ｒ）または緑（Ｇ）の画素値を特定して二値化の基準値を決定してもよい。

続けて、二値化部４１は、Ｓ１１で決定した標準画素値に、所定の範囲内の様々な倍率を乗算して得られる値を二値化の基準値として、当該組織画像のグレースケール画像から二値化画像を生成する（Ｓ１２）。これにより、二値化の基準値をさまざまに異ならせながら、１つの組織画像から複数の二値化画像を生成することができる。なお、１つの組織画像から複数の二値化画像を生成するための二値化の基準値を算出するために用いる倍率は、解析対象である組織の構成に応じて適宜変更し得る。

ただし、二値化の基準値を大きく設定し過ぎると、二値化された組織画像は境界線が過度に強調され、黒っぽい画像になる。一方、二値化の基準値を小さく設定し過ぎると、二値化された組織画像は、境界線の少ない白っぽい画像になる。このことを、図５を用いて説明する。

図５の（ａ）は元の組織画像である。この組織画像では、Ｂの画素値の分布データにおいて出現頻度が最大となる（極大となる）画素値（標準画素値）は１７１である。図５の（ｂ）、（ｃ）、および（ｄ）はＢの標準画素値にそれぞれ１、０．７、０．５を乗算した値を二値化の基準値として用いて二値化された組織画像の一例を示す図である。Ｂに関する標準画素値をそのまま二値化の基準値に用いて（すなわち、倍率＝１）二値化処理をすると、図５の（ａ）に示すように、画像が黒っぽくなり、構成要素間の関係を評価するには不適な画像である。

一方、倍率を０．５にすると、図５の（ｃ）に示すような白っぽい画像となり、この場合も構成要素間の関係を評価するには不適な画像である。図５の（ｂ）は、倍率を０．７にしたときの画像であるが、組織の構成要素間の関係を評価するには適した画像である。黒っぽくなり過ぎたり、白っぽくなり過ぎたりした二値化画像では、１次元ベッチ数が極端に低く算出されてしまい、構成要素間の関係を示す値が適切に算出されない。

このように様々な倍率について二値化画像の適不適を検討した結果、ＨＥ染色された生体組織の組織画像を解析して、当該組織における癌の分化度を判別するためには、Ｂに関する標準画素値に乗算する倍率を０．５５以上０．７２以下の範囲内から選択すればよいことがわかった。

このように、二値化の基準値を変化させる範囲を好ましい範囲に限定することにより、画像解析装置１の処理速度を向上させることができる。

なお、組織の種類や染色方法などに応じて、この二値化の基準値の決め方は適宜変更し得る。標準画素値を特定する色成分はＢ以外の色成分であってもよいし、これに乗算する倍率も適当な範囲の値を用いればよい。

なお、二値化部４１が組織画像を二値化するための二値化の基準値を、すべての画素について共通の（固定の）閾値として決定する方法を適用してもよいし、画素ごとに可変の閾値として決定する方法を適用してもよい。二値化の基準値をすべての画素について共通の（固定の）閾値として決定する方法としては、例えば、以下の（１）〜（４）に示す方法が知られている。
（１）二値化対象の画像を、当該画像の輝度情報を用いてグレースケール画像に変換した後、当該グレースケール画像における画素値の分布データにおいて出現頻度が最大（極大）となる画素値（ピーク値）から、二値化の基準値を決定する方法。
（２）ＲＧＢ画像を、Ｈ（色相）Ｓ（彩度）Ｖ（輝度）変換あるいはＨ（色相）Ｌ（明度）Ｓ（彩度）変換を行い、Ｖ（輝度）、Ｓ（彩度）等のヒストグラムに基づいて、二値化の基準値を決定する方法。
（３）画素値の分布データにおいて出現頻度が最大（極大）となる画素値をそのまま閾値とするのではなく、当該分布データの双峰性を解析し、分離度が最大になる解を求めることで二値化の基準値を決定する方法（例えば、大津の閾値）。
（４）二値化対象の画像の、二値化したい領域の画素数の画像全体に占める割合をパーセント等の比率で指定し、二値化基準値を決定する方法（例：Ｐタイル法）。

（組織画像の解析例）
続いて、本発明の画像解析装置１を用いて生体組織の組織画像を解析した例を、図６〜図８を用いて説明する。図６〜図８は、画像解析装置１によって解析した結果の一例を示す図であり、図６は正常な組織、図７および図８は癌を含む組織の解析例である。

図６の（ａ）に示す組織画像を二値化して二値化画像を生成するために、Ｂに関する標準画素値に、０．５５〜０．７２の値を０．０１ずつ変更して乗算して得られた値を二値化の基準値とした。こうして得られた二値化の基準値を用いて、組織画像のグレースケール画像から複数の二値化画像を生成し、それぞれの二値化画像について１次元ベッチ数を算出した。

図６の（ａ）に示す組織画像の分割領域Ａおよび分割領域Ｂは、癌を含まない正常な領域である。二値化の基準値を異ならせたときの分割領域Ａおよび分割領域Ｂの１次元ベッチ数の変化をプロットしたグラフを図６の（ｂ）に示す。なお、図６〜図８の（ｂ）において、横軸はＢに関する標準画素値に乗算した倍率であり、縦軸は穴の数（すなわち、１次元ベッチ数）である。図６の（ｂ）に示すように、分割領域Ａおよび分割領域Ｂについて倍率に対する穴の数はいずれも近似式（近似関数の式）として二次関数（図中の点線）を用いて近似することが可能であった。近似の程度を示す相関係数（Ｒ^２）はいずれも０．９９より高く、この近似が非常に良い精度の近似であることを示している。

なお、ここでは近似式として二次関数を用いる例を説明したがこれに限定されず、３次関数、４次関数などの高次関数を用いて近似してもよい。また、組織によっては、１次関数を用いて近似してもよい。回帰曲線の描画および近似式の導出には、マイクロソフト社製のExcel（登録商標）に組み込まれたアルゴリズムなど、任意のソフトウェアを用いればよい。

次に、癌を含む領域の解析例を示す。図７の（ａ）に示す組織画像の分割領域Ｃおよび分割領域Ｄは、癌を含む領域である。二値化の基準値を異ならせたときの分割領域Ｃおよび分割領域Ｄの１次元ベッチ数の変化を、プロットしたグラフを図７の（ｂ）に示す。図７の（ｂ）に示すように、分割領域Ｃおよび分割領域Ｄについて倍率に対する穴の数もまた、二次関数（図中の点線）を用いて非常に良い精度で近似することが可能であった。

また、図８の（ａ）に示す組織画像の分割領域Ｅ〜Ｇは、癌を含む領域である。二値化の基準値を異ならせたときの分割領域Ｅ〜Ｇの１次元ベッチ数の変化を、プロットしたグラフを図８の（ｂ）に示す。図８の（ｂ）に示すように、分割領域Ｅ〜Ｇについて倍率に対する穴の数もまた、二次関数（図中の点線）を用いて非常に良い精度で近似することが可能であった。

（二次の係数が負になる場合について）
ここでは、二次関数の二次の係数が負になる場合について、図９を用いて説明する。多くの組織画像において、二次関数の二次の係数は正であったが、分割領域Ｅのように、二次関数の二次の係数が負となる場合も見られた。

図９は、近似式の二次の係数の正/負の割合を示すグラフである。図示のように、正常な組織、腺腫および高分化の癌において、二次の係数が負となる例が確認された。一方、中分化および低分化の癌では、１例を除き、二次の係数が負となることはなかった。しかも、二次の係数が負になっていた低分化の癌の１例では、組織画像の状態が悪く、細胞核が潰れていた。すなわち、近似式の二次の係数が負になる場合は、腺腫および高分化の癌に多く見られ、中分化・低分化の癌では近似式の二次の係数が負になる場合はほとんど見られないことが分かった。

（癌の分化度の判別への適用例）
画像解析装置１を用いて生体組織の組織画像を解析し、癌の分化度を判別するための情報を抽出する例について説明する。近似式の複数の係数に関する複数の座標軸によって形成される座標空間に、解析された組織画像（あるいは各分割画像）についての近似式の複数の係数を座標として有する点（座標点）をプロットしたグラフの一例を図１０に示す。図１０の（ａ）は近似式として用いた二次関数の式を示し、（ｂ）は近似式の二次の係数ａの対数（Log a）と、近似式の軸の位置を示すｂとを座標として有する座標点をプロットした分布図である。対数の底は１０でもよいし、自然数ｅでもよい。なお、ここでは近似式の二次の係数ａの値が大きいため便宜上対数を用いているが、対数を用いなくてもよい。

図１０の（ｂ）に示すように、中分化および低分化の癌を含む画像についての座標点は、縦軸の値および横軸の値が共に大きい位置に、集まって分布している。一方、高分化の癌、腺腫、正常の組織では、座標点は、それぞれ固有の傾きを有する領域（図中の楕円で示した領域）内に分布しており、高分化の癌、腺腫、正常の組織の順にその傾きが緩くなる傾向があることがわかった。

この座標点の分布の傾向に基づいて、当該座標点を算出するに至った組織画像（または分割画像）に写っている癌がどのような分化度の癌であるか、を表す指標を算出して、この指標に基づいて、癌の判定において医師が注目すべき領域が明示された組織画像（判定用画像）を出力情報として生成することができる。この出力情報は、癌の分化度を判別するために有用な情報である。

（判定用画像について）
ここでは、判定用画像について、図１１を用いて説明する。図１１は判定用画像の一例を示す図である。図１１の（ａ）は解析対象の組織画像、図１１の（ｂ）は当該組織画像に写っている癌の分化度に応じて、分割画像を色分けした判定用画像をそれぞれ示している。

例えば、図１１の（ｂ）において、図１１の（ａ）に含まれる領域Ｈおよび領域Ｉに含まれる各分割領域が、当該分割領域の解析結果（上記座標点の位置）に応じて色分けされている。例えば、低分化の癌を含んでいる可能性のある分割領域を赤色、腺腫を含んでいる可能性のある分割領域をオレンジ色、高分化の癌を含んでいる可能性のある分割領域を黄色で示してもよい。また、分化度を１つに特定できない場合（例えば、図１０の（ｂ）において、破線領域が互いに重なっている領域に座標点が存在する場合）、グラデーションによって分化度を表示してもよい。

（二次の係数が負の中分化・高分化の癌の判別について）
なお、前述したように、中分化・高分化の癌については、近似式の二次の係数が負になる場合がある。近似式の二次の係数が負である場合、図１０の（ｂ）に示すようなグラフにプロットすることができない。

そこで、関数特定部４３が算出した２次の係数が負である場合、出力情報生成部４４は、所定の二値化画像における穴の数の実数に基づいて出力情報を生成してもよい。

具体的には、ベッチ数算出部４２は、Ｂに関する標準画素値に０．５５〜０．７２の値を乗算して得られた値を二値化の基準値として生成された複数の二値化画像のそれぞれについて１次元ベッチ数を算出する。次に、関数特定部４３は、二値化の基準値と穴の数との関係を示す二次の近似式に含まれる、二次の係数を含む複数の係数を算出する。出力情報生成部４４は、二次の係数の符号が正の場合には、穴の数の実数と、予め定められた判別基準値とを比較して、比較結果に基づいて癌の分化度の判別するための上述のような判定用画像を出力情報として生成してもよい。

例えば、近似式の二次の係数が負になる場合には、標準画素値に０．７２を乗算して得られる値を二値化の基準値として生成した二値化画像における、穴の数の実数に着目する。例えば、癌を含まない正常な組織では、穴の数は平均１０３５であるのに対し、腺腫では穴の数は平均２５５０に増え、高分化の癌ではさらに平均６２５５にまで増えた。正常、腺腫、および高分化の癌における穴の数の実数の差異は統計学的に有意であった。

このことから、近似式の二次の係数が負になる場合には、二値化画像における穴の数の実数に基づいて、癌の分化度の判別を行うために有用な情報を提供することができることが分かった。

〔実施形態２〕
（画像解析装置１ａの構成）
次に、図１１を用いて、画像解析装置１ａの構成について説明する。図１１は、本発明の実施例２に係る画像解析装置１ａの構成の一例を示すブロック図である。画像解析装置１ａは、表示制御部５の代わりに解析結果送信部９（送信部）を備えている点で、画像解析装置１と異なっている。解析結果送信部９は、外部機器８から受信した画像データが示す撮像画像を解析した結果を出力情報生成部４４から取得し、提示装置１０に送信する。

なお、外部機器８および提示装置１０の数は複数であってもよい。

（画像解析システム）
ここでは、画像解析装置１ａを含む画像解析システム１００および１００ａの構成例について、図１２を用いて説明する。図１２は、本発明に係る画像解析装置１ａを含む画像解析システム１００および１００ａの構成例を示す概略図である。図１２の（ａ）は、外部機器８と提示装置１０とが離れた場所に設置されている例を示しており、（ｂ）は、提示装置１０が外部機器８ａに接続されている例を示している。

画像解析システム１００は、外部機器８、画像解析装置１ａ、および提示装置１０を有している。外部機器８、画像解析装置１ａ、および提示装置１０は、インターネットなどの情報通信ネットワーク５０に接続されており、相互にデータの送受信が可能である。

外部機器８は、例えば、組織を撮像する機能を備える顕微鏡などの機器であってもよいし、組織を撮像した画像を集約して管理するサーバ（電子カルテサーバ、顕微鏡画像データサーバなど）であってもよい。

提示装置１０は、ユーザに画像解析の結果を提示する機能を有する装置であればよく、例えば、提示装置１０はディスプレイを備える表示装置である。あるいは、医療関係者が携帯するタブレット端末などの通信端末機器であってもよい。

組織を撮像した撮像画像の画像データは、外部機器８から画像解析装置１ａへと送信される。画像データを受信した画像解析装置１ａは、画像解析部４が該画像を解析し、出力情報生成部４４が生成した出力情報を解析結果送信部９から提示装置１０または外部機器８に送信する。

画像解析システム１００ａは、外部機器８ａ、画像解析装置１ａ、および提示装置１０を有している。外部機器８ａおよび画像解析装置１ａは、インターネットなどの情報通信ネットワーク５０に接続されており、相互にデータの送受信が可能である。提示装置１０は、外部機器８ａに接続されている。

すなわち、画像解析装置１ａは、遠隔地で撮像された画像を外部機器８から受信して、画像解析を行い、解析結果である出力情報および判定用画像の少なくともいずれか一方を、提示装置１０および１０ａに送信することができる。提示装置１０および１０ａは、外部機器８に接続された装置であってもよいし、画像解析装置１ａおよび外部機器８と独立している装置であってもよい。

〔実施形態３〕
本発明は、生体組織における癌の判定に限らず、他の組織における組織の構成の解析に適用できる。その一例として、本実施形態では、シリコーンゲル（組織）における構成の解析例について、図１４〜図１７を用いて説明する。

（シリコーンゲルの構成および画像）
シリコーンゲルは、ケイ素原子（Ｓｉ）と酸素原子（Ｏ）とを構成要素として含む非晶質材料である。シリコーンゲルは、酸性条件でゲル化させた場合と、塩基性条件でゲル化させた場合とで、分子間の結合（関係）のパターン（構成）が異なることが知られている。以下では、シリコーンゲルにおける分子間の結合を、本発明の画像解析装置１、１ａを用いて解析した結果の一例を示す。

まず、図１４を用いて、シリコーンゲルのゲル化条件に依存した分子構造（シリコーンゲルを構成する分子間の結合パターン）の違いについて説明する。図１４の（ａ）は酸性条件でゲル化したシリコーンゲルの分子構造を示す図であり、（ｂ）は酸性条件でゲル化したシリコーンゲルを、透過型電子顕微鏡を用いて撮像した画像（組織画像、撮像画像）であり、（ｃ）は塩基性条件でゲル化したシリコーンゲルの分子構造を示す図であり、（ｄ）は塩基性条件でゲル化したシリコーンゲルを、透過型電子顕微鏡を用いて撮像した画像（組織画像、撮像画像）である。なお、図１４の（ｂ）および（ｄ）に示す画像において、白く表示される部分は酸素原子が存在していることを示している。ここでは、図１５の（ｂ）および（ｄ）の画像中に１ｎｍを示すスケールバーを描画している。撮像倍率は解析の対象となる組織に応じて、当業者によって適宜設定されればよい。

画像取得部２は、透過型電子顕微鏡などの外部機器８から、シリコーンゲルを撮像した画像を取得する。取得した画像は、記憶部３に格納される。

（シリコーンゲルを撮像した画像の二値化）
まず、二値化部４１は、記憶部３から図１５の（ｂ）の画像を読み出して、当該画像から、二値化の基準値を異ならせながら二値化した二値化画像を生成する。図１５は、図１４の（ｂ）の画像を、二値化の基準値を異ならせながら二値化した二値化画像である。また、図１６は、図１４の（ｄ）の画像を、二値化の基準値を異ならせながら二値化した二値化画像である。

ここで、上述の実施形態と異なり、シリコーンゲルの画像はＲＧＢなどの色成分を含んでいない。そのため、本実施形態では、二値化部４１は、所定の画素値を二値化の基準値として決定する。

図１５および図１６では、（ａ）２００、（ｂ）１７５、（ｃ）１５０、（ｄ）１２５、（ｅ）１００、（ｆ）７５、（ｇ）５０を二値化の基準値としてそれぞれ設定した場合に生成される二値化画像を示す。

図１５および図１６に示す複数の二値化画像のうち、（ａ）に示す二値化画像は二値化の基準値が大きく設定し過ぎており、一方、（ｆ）および（ｇ）に示す二値化画像は二値化の基準値を小さく設定し過ぎている。

（シリコーンゲルの二値化画像における穴の数を算出）
そこで、本実施形態では、ベッチ数算出部４２は、図１５および図１６の二値化画像（ｂ）〜（ｅ）のそれぞれについて穴の数（穴形状の領域の数）を算出した。図１４の（ｂ）および（ｄ）に示す画像において白く表示される部分は、酸素原子が存在していることを示す部分である。それゆえ、二値化画像について穴の数を算出することは、当該二値化画像に含まれる領域内の酸素原子の数を数えている処理、あるいは当該二値化画像に含まれる領域内に酸素原子が統計的にどの程度分布しているかを調べる処理に相当しているということができる。

算出された穴の数は、図１５の（ｂ）、（ｃ）、（ｄ）、および（ｅ）についてはそれぞれ、５９、４５、３８、２９であり、図１６の（ｂ）、（ｃ）、（ｄ）、および（ｅ）についてはそれぞれ、１６２、１３７、５２、６であった。

（二値化の基準値の変化に応じた穴の数の変化を近似する関数の特定）
次に、関数特定部４３は、二値化の基準値の変化に伴う穴の数の変化を近似的に示す関数の特徴を表す数値を算出する。図１７は、シリコーンゲルの構成を解析した結果の一例を示す図であり、縦軸に穴の数、横軸に二値化の基準値をプロットしたグラフである。

図１７に示すように、図１５の（ｂ）〜（ｅ）、および図１６の（ｂ）〜（ｅ）について、二値化の基準値に対する穴の数の変化はいずれも近似式（近似関数）として一次関数（図中の点線）を用いて近似することが可能であった。図１７において、図１５の（ｂ）〜（ｅ）に対応するデータのプロットはデータＬ、図１６の（ｂ）〜（ｅ）に対応するデータのプロットはデータＫでそれぞれ示す。このことは、酸性条件および塩基性条件でゲル化させたシリコーンゲルにはいずれも一定の割合で酸素原子が含まれていることを示している。なお、近似の程度を示す相関係数（Ｒ^２）はいずれも０．９８より高く、この近似が非常に良い精度の近似であることを示していた。

関数特定部４３は、データＬおよびデータＫのそれぞれを近似的に示す一次関数の一次の係数（特徴数）を算出する。なお、関数特定部４３は、一次関数の一次の項の係数に加えて、定数項の値（特徴数）を算出してもよい。

関数特定部４３が特定した一次関数の特徴を示す一次の係数（特徴数）は、シリコーンゲルが酸性条件でゲル化されたものか、あるいは塩基性条件でゲル化されたものかによって、明確に異なっていた。

（シリコーンゲルの判定）
このことは、関数特定部４３が特定した一次関数の特徴を示す一次の係数に基づいて、シリコーンゲルの分子間の結合のパターン（構成）を判別することが可能であることを示している。

そこで、出力情報生成部４４は、関数特定部４３が特定した一次関数の特徴を示す一次の係数に基づいて判定した判定結果を示す出力情報を生成する。これにより、ユーザは、表示装置６の画面において提示される出力情報に基づいて、シリコーンゲルの分子構造を判別することができる。

なお、シリコーンゲルの分子構造を判別する場合を例に挙げて説明したが、これに限定されない。本発明に係る画像解析装置１、１ａを用いれば、任意の組織について、該組織を構成する複数の構成要素間の結びつき（関係）を、高い精度で評価することが可能である。

〔ソフトウェアによる実現例〕
画像解析装置１、１ａの制御ブロック（特に二値化部４１、ベッチ数算出部４２、関数特定部４３、および出力情報生成部４４）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ＣＰＵ（Central Processing Unit）を用いてソフトウェアによって実現してもよい。

後者の場合、画像解析装置１、１ａは、各機能を実現するソフトウェアであるプログラムの命令を実行するＣＰＵ、上記プログラムおよび各種データがコンピュータ（またはＣＰＵ）で読み取り可能に記録されたＲＯＭ（Read Only Memory）または記憶装置（これらを「記録媒体」と称する）、上記プログラムを展開するＲＡＭ（Random Access Memory）などを備えている。そして、コンピュータ（またはＣＰＵ）が上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。

本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
〔まとめ〕
本発明の態様１に係る画像解析装置１、１ａは、組織を撮像した１つの撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化部４１と、上記二値化部が生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出部（ベッチ数算出部４２）と、上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定部（関数特定部４３）と、上記特定部が特定した上記特徴数に対応した出力情報を生成する出力情報生成部４４とを備えている。

この構成によれば、画像解析装置は、組織を撮像した１つの撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する。ここで、穴形状の領域とは、二次元平面である撮像画像における空間（穴）である。穴形状の領域の数を算出する処理としては、上記複数の二値化画像の各々について１次元ベッチ数を求める処理を行ってもよい。なお、画像からベッチ数を算出するためのプログラムは公知であり、そのプログラムを実行する装置を空間数算出部として用いればよい。

次に、画像解析装置は、複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する。そして、画像解析装置は、特定した上記特徴数に対応した出力情報を生成する。なお、特徴数は少なくとも１つ特定されればよい。

それゆえ、構成要素間の関係を数理的に表現した情報を出力情報として生成することができ、さまざまな組織の構成を判別することが可能となる。なお、ここでいう「組織」とは、１または複数種類の構成要素が、一定のパターンで集合した構造単位を意味しており、細胞などの集まりから構成される生体組織、結晶粒などの集まりから構成される鉱物組織、およびシリコーンなどの樹脂材料の分子の集まりから構成される樹脂製品の組織などであってもよい。

例えば、生体組織に生じた癌の分化度を判定する場合、生体組織を構成する構成要素間の関係（結びつきの強さ）が変化するため、この構成要素間の関係を数値化することにより、癌の分化度を判別するために有用な情報を抽出することができる。

本発明の態様２に係る画像解析装置は、上記態様１において、上記特定部は、上記複数の二値化の基準値と上記穴形状の領域の数との関係を近似的に示す関数の特徴を表す複数の特徴数のうち少なくとも１つを特定してもよい。

すなわち、本態様に係る画像解析装置は、組織を撮像した１つの撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化部４１と、上記二値化部が生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出部（ベッチ数算出部４２）と、上記複数の二値化の基準値と上記穴形状の領域の数との関係を近似的に示す関数の特徴を表す複数の特徴数を特定する関数特定部４３と、上記関数特定部が特定した上記複数の特徴数のうちの少なくとも１つに対応した出力情報を生成する出力情報生成部４４とを備えていてもよい。

これにより、複数の二値化の基準値と穴形状の領域の数との関係を近似的に示す関数の特徴を表す複数の特徴数を、当該複数の二値化の基準値と穴形状の領域の数との対応関係を特徴づける複数の特徴数として利用することができる。

本発明の態様３に係る画像解析装置は、上記態様１または２において、上記出力情報生成部は、上記複数の特徴数のそれぞれに１つずつが対応している複数の座標軸によって規定される座標空間における、上記特定部が特定した複数の特徴数を座標として有する点（座標点）の位置に対応した出力情報を生成してもよい。

これにより、組織の構成要素間の関係を高い精度で評価するために有用な情報を出力することができる。

解析対象が生体組織の場合、上記座標空間における上記座標点の位置によって、当該生体組織における癌の分化度を判定することができる可能性があることを本発明の発明者は見出した。

上記のように構成すれば、上記座標点の位置に対応した出力情報を生成することができ、この出力情報を病理医に提示することにより、病理医による癌の分化度の判定を適切にサポートすることができる。

本発明の態様４に係る画像解析装置は、上記態様１から３のいずれかにおいて、上記撮像画像を所定の大きさの分割領域に分割する分割部（二値化部４１）をさらに備え、上記空間数算出部は、上記分割部が生成した分割領域毎に上記穴形状の領域の数を算出し、上記特定部は、上記分割領域毎に上記複数の特徴数を特定し、上記出力情報生成部は、上記分割領域毎に上記出力情報を生成してもよい。

上記の構成によれば、分割部は、組織の撮像画像を所定の大きさの分割領域に分割し、空間数算出部は、分割部が生成した分割領域ごとに穴形状の領域の数を算出する。そして、特定部は、各分割領域に関して空間数算出部が算出した穴形状の領域の数から分割領域ごとに、関数の特徴を表す複数の特徴数を特定し、出力情報生成部は、分割領域ごとに特定された特徴数を座標として有する点の位置に対応した出力情報を生成する。

それゆえ、組織の撮像画像全体を解析する場合よりも組織の構成要素間の関係を数値化する精度を高めることができるとともに、撮像画像のどの分割領域が注目すべき領域であるかを特定することができる。

本発明の態様５に係る画像解析装置は、上記態様４において、上記出力情報生成部は、上記出力情報として、上記分割領域を含む判定用画像を生成し、上記撮像画像における複数の分割領域のそれぞれについて、上記複数の特徴数に応じた態様で、上記判定用画像における当該分割領域の画像を生成してもよい。

これにより、ユーザに対して、組織の構成を判定するために必要な情報を含み、かつ直観的で利便性の高い判定用画像を提示することができる。

本発明の態様６に係る画像解析装置は、上記態様１から５のいずれかにおいて、上記撮像画像に含まれる画素の画素値に基づいて上記二値化の基準値を決定してもよい。

これにより、組織を撮像した撮像画像に写っている構成要素間の関係をより明確に特定できる画素値に基づいて、二値化画像を生成する際の二値化の基準値を決定することができる。

なお、生体組織の場合、多くは染色された組織を切片化して撮像される。組織の染色に用いた色素に応じて、二値化画像を生成するための二値化の基準値を適切に決定することができる。この場合、撮像画像を構成する所定の色（例えば青色（Ｂ））の画素値のうち、最も高い頻度で含まれているＢの画素値を標準画素値として決定し、この標準画素値に基づいて二値化の基準値を決定すればよい。

本発明の態様７に係る画像解析装置は、上記態様１から６のいずれかにおいて、上記特定部は、上記二値化の基準値と上記穴形状の領域の数との関係を示すｎ次式（ｎ≧１）に含まれる、ｎ次の係数を含む複数の係数を特定してもよい。

二値化の基準値と穴形状の領域の数とを座標として有する点をグラフにプロットした場合、二値化の基準値を変えたときの穴形状の領域の数の変化は、ｎ次関数の式で近似することができる。すなわち、二値化の基準値と上記穴形状の領域の数との関係性はｎ次関数の式で近似できることが多い。

上記の構成によれば、ｎ次式の係数を含む複数の係数を特定する。これにより、二値化の基準値と上記穴形状の領域の数との関係を、客観的で簡便な数学的処理に基づいて抽出することができる。この場合、算出される係数は関数の特徴を表す複数の特徴数と捉えることができる。

本発明の態様８に係る画像解析システムは、上記態様１から７のいずれかに記載の画像解析装置と、上記撮像画像の画像データを、上記画像解析装置へ送信する外部機器と、上記画像解析装置において生成された上記出力情報を取得して、該出力情報を提示する提示装置とを含んでいる。

これにより、例えば遠隔地の外部機器を用いて撮像された撮像画像を受信して、当該画像を解析することができ、かつ撮像画像および出力情報を遠隔地のユーザに提示することができる。

また、本発明の態様９に係る画像解析方法は、上記の課題を解決するために、組織を撮像した撮像画像を解析する画像解析装置１、１ａにおける画像解析方法であって、１つの上記撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化工程（Ｓ２）と、上記二値化工程にて生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出工程（Ｓ３）と、上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定工程（Ｓ４）と、上記特定工程にて特定した上記特徴数に対応した出力情報を生成する出力情報生成工程（Ｓ５）とを含んでいる。この構成によれば、上記態様１と同様の効果を奏する。

また、画像解析装置を動作させるための画像解析プログラムであって、コンピュータを上記各部として機能させるための画像解析プログラムおよび当該画像解析プログラムを記録したコンピュータ読み取り可能な記録媒体も本発明の技術的範囲に含まれる。

１、１ａ 画像解析装置
２ 画像取得部（受信部）
４ 画像解析部（画像解析装置）
５ 表示制御部
６ 表示装置（提示部）
８ 外部機器
９ 解析結果送信部（送信部）
１０ 提示装置
４１ 二値化部(分割部)
４２ ベッチ数算出部（空間数算出部）
４３ 関数特定部（特定部）
４４ 出力情報生成部
１００、１００ａ 画像解析システム
Ｓ２ 二値化工程
Ｓ３ 空間数算出工程
Ｓ４ 特定工程
Ｓ５ 出力情報生成工程

Claims (11)

1. 組織を撮像した１つの撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化部と、
   上記二値化部が生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出部と、
   上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定部と、
   上記特定部が特定した上記特徴数に対応した出力情報を生成する出力情報生成部とを備えることを特徴とする画像解析装置。
2. 上記特定部は、上記複数の二値化の基準値と上記穴形状の領域の数との関係を近似的に示す関数の特徴を表す複数の特徴数のうち少なくとも１つを特定することを特徴とする請求項１に記載の画像解析装置。
3. 上記出力情報生成部は、上記複数の特徴数のそれぞれに１つずつが対応している複数の座標軸によって規定される座標空間における、上記特定部が算出した複数の特徴数を座標として有する点の位置に対応した出力情報を生成することを特徴とする請求項１または２に記載の画像解析装置。
4. 上記撮像画像を所定の大きさの分割領域に分割する分割部をさらに備え、
   上記空間数算出部は、上記分割部が生成した分割領域毎に上記穴形状の領域の数を算出し、
   上記特定部は、上記分割領域毎に上記複数の特徴数を算出し、
   上記出力情報生成部は、上記分割領域毎に上記出力情報を生成することを特徴とする請求項１から３のいずれか１項に記載の画像解析装置。
5. 上記出力情報生成部は、上記出力情報として、上記分割領域を含む判定用画像を生成し、
   上記撮像画像における複数の分割領域のそれぞれについて、上記複数の特徴数に応じた態様で、上記判定用画像における当該分割領域の画像を生成することを特徴とする請求項４に記載の画像解析装置。
6. 上記二値化部は、上記撮像画像に含まれる画素の画素値に基づいて上記二値化の基準値を決定することを特徴とする請求項１から５のいずれか１項に記載の画像解析装置。
7. 上記特定部は、上記二値化の基準値と上記穴形状の領域の数との関係を示すｎ次式（ｎ≧１）に含まれる、ｎ次の係数を含む複数の係数を算出することを特徴とする請求項１から６のいずれか１項に記載の画像解析装置。
8. 請求項１から７のいずれか１項に記載の画像解析装置と、
   上記撮像画像の画像データを、上記画像解析装置へ送信する外部機器と、
   上記画像解析装置において生成された上記出力情報を取得して、該出力情報を提示する提示装置とを含むことを特徴とする画像解析システム。
9. 組織を撮像した撮像画像を解析する画像解析装置における画像解析方法であって、
   １つの上記撮像画像に対して、二値化の基準値を異ならせながら複数回の二値化処理を行うことにより、複数の二値化画像を生成する二値化工程と、
   上記二値化工程にて生成した複数の二値化画像のそれぞれについて、穴形状の領域の数を算出する空間数算出工程と、
   上記複数の二値化の基準値と上記穴形状の領域の数との対応関係を特徴づける特徴数を特定する特定工程と、
   上記特定工程にて特定した上記特徴数に対応した出力情報を生成する出力情報生成工程とを含むことを特徴とする画像解析方法。
10. 請求項１に記載の画像解析装置としてコンピュータを機能させるための画像解析プログラムであって、上記二値化部、上記空間数算出部、上記特定部、および上記出力情報生成部としてコンピュータを機能させるための画像解析プログラム。
11. 請求項１０に記載の画像解析プログラムを記録したコンピュータ読み取り可能な記録媒体。