(57)【要約】【課題】  ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別することができる生体組織解析装置、生体組織解析プログラムおよび生体組織解析方法を提供する。【解決手段】  二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析装置であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析装置；（i）ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；（ii）生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出部；（iii）教師なし学習によって生体組織の変質状態を大別する変質状態大別部；（iv）教師あり学習によって生体組織の変質状態を識別する変質状態識別部。【選択図】  図１
2019163981.tif 000002
【特許請求の範囲】
【請求項１】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析装置であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析装置；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【請求項２】
前記ハイパースペクトルデータの取得時に用いた光源のスペクトル情報によって、前記解析対象部位に相当するピクセルの前記スペクトル情報を除算する光源影響除去部を有している、請求項１に記載の生体組織解析装置。
【請求項３】
前記解析対象部位に相当するピクセルの前記スペクトル情報に基づいて、前記解析対象部位の染色濃度によって変化するスペクトル成分を除去する染色影響除去部を有している、請求項１または請求項２に記載の生体組織解析装置。
【請求項４】
前記解析対象部位抽出部は、前記スペクトル情報の類似度に基づいて各ピクセルを解析対象部位とその他の部位とにラベリングし、前記解析対象部位としてラベリングされたピクセルをフィルタリングする、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項５】
前記解析対象部位抽出部は、前記生体組織の二次元画像を画像認識処理によって腺管を特定するとともに、前記腺管上における前記細胞核を解析対象部位として画像認識処理によって特定する、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項６】
前記変質状態は、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）およびガン細胞の４ステージ、またはポリープにおける、軽度異型腺腫（ＬＧＡ：Low Grade Adenoma）、高度異型腺腫（ＨＧＡ：High Grade Adenoma）およびガン細胞の３ステージによって特定される、請求項１から請求項５のいずれかに記載の生体組織解析装置。
【請求項７】
前記ハイパースペクトルデータの波長帯域に可視光が含まれている、請求項１から請求項６のいずれかに記載の生体組織解析装置。
【請求項８】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析プログラムであって、下記（i），（ii）としてコンピュータを機能させるとともに、下記（iii）および／または（iv）としてコンピュータを機能させる、前記生体組織解析プログラム；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【請求項９】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析方法であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析方法；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得ステップ；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出ステップ；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別ステップ；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別ステップ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ハイパースペクトルデータを用いて生体組織を解析する技術に関し、特に、細胞核のガン化状態等のように、生体組織の変質状態を識別するのに好適な生体組織解析装置、生体組織解析プログラムおよび生体組織解析方法に関するものである。
【背景技術】
【０００２】
従来、細胞核のガン化状態等のように、生体組織の変質状態を診断するに際しては、病理専門医による顕微鏡画像を利用した形態学的診断が行われている。しかしながら、そのような形態学的診断は、病理専門医の知識や経験によるところが大きく、一人前になるまでには数年に渡るトレーニングが必要とされるほど困難なものとされている。
【０００３】
そこで、近年、分光スペクトルを用いて生体組織を光学的に分析する方法が検討されている。例えば、特開２０１７−２０３６３７号公報には、試料に含まれる細胞の分光スペクトルに基づいて、統計的手法、機械学習又はパターン認識により、当該細胞が腫瘍細胞であるか否かを判定する分析工程を有する腫瘍細胞検出方法が提案されている（特許文献１）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１７−２０３６３７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に記載された発明においては、分析対象の細胞が腫瘍細胞であるか否かを判定するものに過ぎない。このため、例えば、発症初期の大腸癌等における４段階のガン化状態（非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）、ガン細胞）等のように、生体組織の変質状態を細かく識別することはできないという問題がある。
【０００６】
本発明は、このような問題点を解決するためになされたものであって、ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別することができる生体組織解析装置、生体組織解析プログラムおよび生体組織解析方法を提供することを目的としている。
【課題を解決するための手段】
【０００７】
本発明に係る生体組織解析装置は、ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別するという課題を解決するために、二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析装置であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析装置；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【０００８】
また、本発明の一態様として、生体組織の変質状態を予め大別することで、解析精度を向上するという課題を解決するために、前記変質状態識別部によって識別する前に、前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類する変質状態大別部を有していてもよい。
【０００９】
さらに、本発明の一態様として、解析用データに対する光源の影響を除去し、解析精度を向上するという課題を解決するために、前記ハイパースペクトルデータの取得時に用いた光源のスペクトル情報によって、前記解析対象部位に相当するピクセルの前記スペクトル情報を除算する光源影響除去部を有していてもよい。
【００１０】
また、本発明の一態様として、解析用データに対する染色の影響を除去し、解析精度を向上するという課題を解決するために、前記解析対象部位に相当するピクセルの前記スペクトル情報に基づいて、前記解析対象部位の染色濃度によって変化するスペクトル成分を除去する染色影響除去部を有していてもよい。
【００１１】
さらに、本発明の一態様として、解析対象部位をピクセル単位で抽出するという課題を解決するために、前記解析対象部位抽出部は、前記スペクトル情報の類似度に基づいて各ピクセルを解析対象部位とその他の部位とにラベリングし、前記解析対象部位としてラベリングされたピクセルをフィルタリングしてもよい。
【００１２】
また、本発明の一態様として、解析対象部位を細胞核単位で抽出するという課題を解決するために、前記解析対象部位抽出部は、前記生体組織の二次元画像を画像認識処理によって腺管を特定するとともに、前記腺管上における前記細胞核を解析対象部位として画像認識処理によって特定してもよい。
【００１３】
さらに、本発明の一態様として、潰瘍やポリープにおける細胞核のガン化状態を病理診断に従って識別するという課題を解決するために、前記変質状態は、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）およびガン細胞の４ステージによって特定されてもよく、またはポリープにおける、軽度異型腺腫（ＬＧＡ：Low Grade Adenoma）、高度異型腺腫（ＨＧＡ：High Grade Adenoma）およびガン細胞の３ステージによって特定されてもよい。
【００１４】
また、本発明の一態様として、生体組織の変質状態を解析するのに適したハイパースペクトルデータを取得するという課題を解決するために、前記ハイパースペクトルデータの波長帯域に可視光が含まれていてもよい。
【００１５】
また、本発明に係る生体組織解析プログラムは、ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別するという課題を解決するために、二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析プログラムであって、下記（i），（ii）としてコンピュータを機能させるとともに、下記（iii）および／または（iv）としてコンピュータを機能させる、前記生体組織解析プログラム；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【００１６】
さらに、本発明に係る生体組織解析方法は、ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別するという課題を解決するために、二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析方法であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析方法；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得ステップ；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する解析対象部位抽出ステップ；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別ステップ；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別ステップ。
【発明の効果】
【００１７】
本発明によれば、ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別することができる。
【図面の簡単な説明】
【００１８】
【図１】本発明に係る生体組織解析装置の一実施形態を示すブロック図である。
【図２】ハイパースペクトルデータの一例を示す図である。
【図３】本実施形態の解析対象部位抽出部によって実行される、（ａ）スペクトル情報抽出処理、（ｂ）ラベリング処理および（ｃ）フィルタリング処理を示す図である。
【図４】本実施形態の生体組織解析装置および生体組織解析プログラムによって実行される生体組織解析方法を示すフローチャートである。
【図５】実施例１において、各ステージにおける全特徴量の組み合わせについての識別率を示す表である。
【図６】実施例２において、主成分分析によって特定された各主成分の寄与率を示すグラフである。
【図７】実施例３において、ロジスティック回帰による予測結果の混同行列を示す表である。
【図８】実施例４において、（ａ）サポートベクターマシンのみによる解析結果、および（ｂ）クラスター分析とサポートベクターマシンによる解析結果を示す表である。
【図９】実施例５において、ハイパースペクトルデータの独立成分分析結果を示すグラフである。
【発明を実施するための形態】
【００１９】
以下、本発明に係る生体組織解析装置、生体組織解析プログラムおよび生体組織解析方法の実施形態について図面を用いて説明する。
【００２０】
本実施形態の生体組織解析システムは、図１に示すように、生体組織のハイパースペクトルデータを取得するハイパースペクトルカメラ１０と、このハイパースペクトルカメラ１０によって取得されたハイパースペクトルデータに基づいて、生体組織の変質状態を解析する生体組織解析装置１とから構成されている。以下、各構成について説明する。
【００２１】
なお、本発明において、生体組織とは、人間や動物等を含む生体において、複数の細胞が所定のパターンで集合してなる全ての組織を含む概念である。また、本発明において、変質状態とは、ガンの進行度やガンの悪性度等のように、生体組織がどの程度変質しているかを示す全ての状態を含む概念である。
【００２２】
ハイパースペクトルカメラ１０は、二次元の空間情報と、複数の波長におけるスペクトル情報（ハイパースペクトル情報）とを同時に取得するものである。具体的には、図２に示すように、二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを取得するものである。なお、スペクトル情報は、分光スペクトルとも呼ばれ、各バンド（波長）ごとの光強度の分布を示すものである。
【００２３】
本実施形態において、上記二次元画像を構成する単位ピクセルは、生体組織における解析対象部位の内部に含まれる大きさに設定されており、当該解析対象部位について複数のピクセルを有するハイパースペクトルデータを取得するようになっている。また、本実施形態では、潰瘍やポリープにおけるガン化状態（ガンの進行度）を解析する場合、細胞核が解析対象部位とされる。しかしながら、解析対象部位は、リンパ球、粘液、組織球（白血球）等のように、解析しようとする生体組織の変質状態に応じて適宜選択され得る。
【００２４】
また、本実施形態において、生体組織は、ホルマリン等で固定してからパラフィン等で包埋して切片を作成し、ヘマトキシリン・エオシン染色等を施した上でプレパラートに挟まれている。そして、当該プレパラートを顕微鏡にセットした状態で、可視光源によって照射して得られる透過光を拡大し、当該透過光を接眼レンズに当接させたハイパースペクトルカメラ１０によって計測するようになっている。
【００２５】
なお、一般的な医療機器に使用される赤外線や紫外線は、顕微鏡レンズの透過性（光学コーディング性能）が可視光域に比べて低く、収差の影響も受けやすい。また、顕微鏡レンズに使用されるフッ化カルシウムは、近赤外光の励起で蛍光を発するため、弱いラマン信号の測定を妨害することが知られている。一方、可視光はＸ線に比べてフォトンエネルギーが小さく、人体への影響が少ない。また、可視光の分光情報は、細胞核内のクロマチン濃度の影響で変化することが知られており、ラマン分光と同様にクロマチンだけでなく多くの細胞組織の情報を潜在的に含んでいる可能性が高い。そこで、本実施形態では、ハイパースペクトルデータの波長帯域に可視光が含まれるように可視光源を使用している。
【００２６】
つぎに、生体組織解析装置１は、パーソナルコンピュータ等のコンピュータによって構成されており、図１に示すように、主として、表示手段２と、入力手段３と、記憶手段４と、演算処理手段５とを有している。以下、各構成手段について詳細に説明する。
【００２７】
表示手段２は、液晶ディスプレイ等で構成されており、生体組織解析装置１による解析結果等を表示するものである。入力手段３は、キーボードやマウス等で構成されており、ユーザによる指示や選択を入力するものである。なお、本実施形態では、表示機能のみを備えた表示手段２と入力機能のみを備えた入力手段３とを別々に使用しているが、この構成に限定されるものではなく、表示機能および入力機能を兼ね備えたタッチパネル等の表示入力手段を使用してもよい。
【００２８】
記憶手段４は、各種のデータを記憶するとともに、演算処理手段５が演算処理を行う際のワーキングエリアとして機能するものである。本実施形態において、記憶手段４は、ハードディスク、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、フラッシュメモリ等で構成されており、図１に示すように、プログラム記憶部４１と、教師データ記憶部４２とを有している。以下、各構成部についてより詳細に説明する。
【００２９】
プログラム記憶部４１には、本実施形態の生体組織解析装置１を制御するための生体組織解析プログラム１ａがインストールされている。そして、演算処理手段５が、当該生体組織解析プログラム１ａを実行することにより、生体組織解析装置１としてのコンピュータを後述する各構成部として機能させるようになっている。
【００３０】
なお、生体組織解析プログラム１ａの利用形態は、上記構成に限られるものではない。例えば、ＣＤ−ＲＯＭやＤＶＤ−ＲＯＭ等のように、コンピュータで読み取り可能な記録媒体に生体組織解析プログラム１ａを記憶させておき、当該記録媒体から直接読み出して実行してもよい。また、外部サーバ等からクラウドコンピューティング方式やＡＳＰ（Application Service Provider）方式で利用してもよい。
【００３１】
教師データ記憶部４２は、生体組織に関するスペクトル情報としての教師データを記憶するものである。本実施形態では、潰瘍における細胞核を解析対象部位としているため、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）およびガン細胞という４ステージによって特定される変質状態のそれぞれについて、スペクトル情報が教師データとして登録されている。
【００３２】
なお、教師データは、上記に限定されるものではなく、解析しようとする生体組織の変質状態に応じて適宜選択される。例えば、ポリープにおける細胞核のガン化状態を解析する場合には、ポリープにおける、軽度異型腺腫（ＬＧＡ：Low Grade Adenoma）、高度異型腺腫（ＨＧＡ：High Grade Adenoma）およびガン細胞という３ステージによって特定される変質状態のそれぞれについて、スペクトル情報が登録されることとなる。
【００３３】
演算処理手段５は、ハイパースペクトルカメラ１０からハイパースペクトルデータを取得し、生体組織の解析処理を実行するものである。本実施形態において、演算処理手段５は、ＣＰＵ（Central Processing Unit）等で構成されており、記憶手段４にインストールされた生体組織解析プログラム１ａを実行することにより、図１に示すように、ハイパースペクトルデータ取得部５１と、解析対象部位抽出部５２と、光源影響除去部５３と、染色影響除去部５４と、データ数削減部５５と、変質状態大別部５６と、変質状態識別部５７として機能するようになっている。以下、各構成部についてより詳細に説明する。
【００３４】
ハイパースペクトルデータ取得部５１は、ハイパースペクトルカメラ１０から生体組織のハイパースペクトルデータを取得するものである。本実施形態において、ハイパースペクトルデータ取得部５１は、ハイパースペクトルカメラ１０に制御信号を出力して生体組織をスキャンさせ、撮像素子から出力された画像信号を１ラインごとにハイパースペクトルデータとして取得するようになっている。
【００３５】
解析対象部位抽出部５２は、ハイパースペクトルデータに基づいて、生体組織の二次元画像から解析対象部位に相当するピクセルを抽出するものである。本実施形態において、解析対象部位抽出部５２は、まず、図３（ａ）に示すように、二次元画像を構成する各ピクセルのスペクトル情報を抽出する。つぎに、図３（ｂ）に示すように、抽出したスペクトル情報の類似度に基づいて、各ピクセルを解析対象部位とその他の部位（リンパ球、組織球、粘液等）とにラベリングする。そして、図３（ｃ）に示すように、解析対象部位としてラベリングされたピクセルをフィルタリングする。このフィルタリングされた各ピクセルのスペクトル情報（各バンドごとの光強度）が解析用データとなる。
【００３６】
なお、本実施形態において、解析対象部位抽出部５２は、ベクトル間距離やマハラノビス距離等によってスペクトル情報の類似度を算出するようになっている。また、本実施形態では、解析対象部位およびその他の各部位ごとに類似度の閾値が設定されており、いずれかの閾値以上の類似度を有するピクセルが、当該閾値に対応する部位としてラベリングされるようになっている。
【００３７】
また、本実施形態において、解析対象部位抽出部５２は、スペクトル情報の類似度に基づいて解析対象部位を抽出しているが、この構成に限定されるものではなく、他の抽出方法を採用してもよい。具体的には、潰瘍やポリープは腺管構造を有しており、当該腺管上に解析対象部位となる細胞核（腺上皮細胞）が存在している。このため、本実施形態のように、潰瘍やポリープにおけるガン化状態を解析する場合、解析対象部位抽出部５２は、ハイパースペクトルデータから生体組織の二次元画像（ＲＧＢ画像）を合成し、当該二次元画像を画像認識処理することによって腺管を特定するとともに、当該腺管上における細胞核を解析対象部位として画像認識処理によって特定してもよい。
【００３８】
光源影響除去部５３は、光源が解析用データに与える影響を除去するものである。本実施形態において、光源影響除去部５３は、ハイパースペクトルデータの取得時に用いた光源のスペクトル情報によって、解析対象部位に相当するピクセルのスペクトル情報を除算することにより、光源の影響を除去するようになっている。
【００３９】
具体的には、光源影響除去部５３は、解析用データの各バンドごとに下記式（１）によって分光反射率Ｒｉを算出する。
Ｒｉ＝Ｉｉ／Ｉ０  …式（１）
ただし、各符号は以下を表す。
Ｉｉ：ｉ番目の細胞核のスペクトル情報
Ｉ０：光源のスペクトル情報
【００４０】
なお、本実施形態では、上記のとおり、解析対象部位抽出部５２が、解析対象部位に相当するピクセルを抽出しているため、ピクセルごとに分光反射率Ｒｉを算出している。しかしながら、上記のとおり、解析対象部位抽出部５２が、画像認識処理によって解析対象部位を抽出する場合、当該解析対象部位ごとに分光反射率Ｒｉを算出する。また、ハイパースペクトルデータの取得時における光源の条件が一定であれば、光源影響除去部５３を機能させる必要はなく、スペクトル情報がそのまま解析用データとなる。
【００４１】
染色影響除去部５４は、解析対象部位に対する染色が解析用データに与える影響を除去するものである。本実施形態において、染色影響除去部５４は、解析対象部位に相当するピクセルのスペクトル情報に基づいて、解析対象部位の染色濃度によって変化するスペクトル成分を除去することにより、染色による影響を除去するようになっている。
【００４２】
なお、解析対象部位の染色濃度によって変化するスペクトル成分は、実施例５で後述するとおり、独立成分分析等によって事前に特定される。また、解析対象部位の染色濃度が一定であれば、染色影響除去部５４を機能させる必要はなく、上記スペクトル成分を除去していないスペクトル情報がそのまま解析用データとなる。
【００４３】
データ数削減部５５は、解析用データのデータ数を削減するものである。本実施形態において、データ数削減部５５は、主成分分析によって解析用データの次元を削減することにより、解析速度を向上するようになっている。なお、データ数の削減方法は、主成分分析に限定されるものではなく、解析用データに含まれる任意のバンドのみを選択することにより、データ数を削減してもよい。また、解析速度よりも解析精度を優先する場合は、データ数削減部５５を機能させずに、全ての解析用データについて解析すればよい。
【００４４】
変質状態大別部５６は、解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別するものである。本実施形態において、変質状態大別部５６は、教師なし学習として、異なる性質のものが混ざり合った集団から、互いに似た性質を持つものを集め、クラスターを作るクラスター分析（ウォード法）を採用している。このため、本実施形態のように、潰瘍におけるガン化状態を解析する場合、実施例１で後述するとおり、上述した４ステージのうち、識別しやすい非ガン細胞および軽度異形成（ＬＧＤ）の２ステージが高精度に大別される。
【００４５】
なお、教師なし学習としては、クラスター分析に限定されるものではなく、生体組織の変質状態を大別しうるものであれば適宜採用することができる。また、変質状態大別部５６による分類は、変質状態識別部５７による識別の前に行うことによって識別対象を絞り込み、識別精度および汎用性を向上するものであるが、必ずしも実行する必要はない。さらに、ポリープにおけるガン化状態を解析する場合は、非ガン細胞というステージはないため、変質状態大別部５６による分類は不要である。
【００４６】
変質状態識別部５７は、解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、生体組織の変質状態を識別するものである。本実施形態において、変質状態識別部５７は、教師あり学習として、高い分類機能を有するロジスティック回帰を採用している。このため、実施例３で後述するとおり、変質状態大別部５６では識別が困難であった高度異形成（ＨＧＤ）およびガン細胞の２ステージについても高精度に識別される。
【００４７】
なお、教師あり学習としては、ロジスティック回帰に限定されるものではなく、生体組織の変質状態を識別しうるものであれば、サポートベクターマシン（ＳＶＭ：Support Vector Machine）やランダムフォレスト等を適宜採用することができる。また、変質状態大別部５６によって非ガン細胞および軽度異形成（ＬＧＤ）を予め大別しない場合には、変質状態識別部５７によって直接、上記４ステージのそれぞれを識別してもよい。さらに、変質状態大別部５６による分類のみで十分な場合は、変質状態識別部５７による識別を実行する必要はない。
【００４８】
つぎに、本実施形態の生体組織解析装置１、生体組織解析プログラム１ａおよび生体組織解析方法による作用について説明する。
【００４９】
本実施形態の生体組織解析装置１、生体組織解析プログラム１ａおよび生体組織解析方法を用いて生体組織の変質状態を解析する場合、図４に示すように、まず、ハイパースペクトルデータ取得部５１が、ハイパースペクトルカメラ１０から生体組織のハイパースペクトルデータを取得する（ステップＳ１：ハイパースペクトルデータ取得ステップ）。
【００５０】
このとき、本実施形態では、赤外線や紫外線と比較して顕微鏡レンズの透過性が高く、収差の影響を受けにくい可視光源を使用する。また、可視光は、Ｘ線と比較して人体への影響が少なく、その分光スペクトルには多くの細胞組織に関する情報が潜在的に含まれている。よって、生体組織の変質状態を解析するのに適したハイパースペクトルデータが取得される。
【００５１】
つぎに、解析対象部位抽出部５２が、ハイパースペクトルデータ取得部５１によって取得されたハイパースペクトルデータに基づいて、生体組織の二次元画像から解析対象部位に相当するピクセルを抽出する（ステップＳ２：解析対象部位抽出ステップ）。これにより、当該ピクセルのスペクトル情報が解析用データとして出力されるため、当該解析用データに対して、解析モードに応じた処理が実行される（ステップＳ３）。
【００５２】
具体的には、高精度モードに設定されている場合（ステップＳ３：ＹＥＳ）、光源影響除去部５３が、解析対象部位に相当する各ピクセルのスペクトル情報を光源のスペクトル情報で除算する（ステップＳ４：光源影響除去ステップ）。これにより、解析用データに対する光源の影響が除去されるため、解析精度が向上する。
【００５３】
つづいて、染色影響除去部５４が、解析対象部位の染色濃度によって変化するスペクトル成分を除去する（ステップＳ５：染色影響除去ステップ）。これにより、解析用データに対する染色の影響が除去されるため、解析精度が向上する。
【００５４】
一方、高精度モードに設定されていない場合（ステップＳ３：ＮＯ）、データ数削減部５５が、解析用データのデータ数を削減する（ステップＳ６：データ数削減ステップ）。これにより、実施例２で後述するとおり、一定の解析精度が維持されたまま、解析速度が向上する。
【００５５】
解析モードに応じた処理が実行された後、変質状態大別部５６が、解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類する（ステップＳ７：変質状態大別ステップ）。これにより、生体組織の変質状態が大別され、変質状態識別部５７によって識別すべき対象が絞り込まれるため、解析精度が向上する。
【００５６】
つづいて、変質状態識別部５７が、解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類する（ステップＳ８：変質状態識別ステップ）。これにより、生体組織の変質状態が高精度に識別される。なお、本実施形態において、解析結果は、全ピクセル数に対する、変質状態のそれぞれに分類されたピクセル数の割合によって表示手段２等に表示される（ステップＳ９：解析結果表示ステップ）。
【００５７】
なお、本発明において、上述したステップＳ１およびステップＳ２は、生体組織解析装置１とは別体の解析装置（プログラム）に実行させてもよい。すなわち、別体の解析装置によって、ハイパースペクトルカメラ１０から生体組織のハイパースペクトルデータを取得させ、当該ハイパースペクトルデータに基づいて、解析対象部位に相当するピクセルを抽出させる。そして、当該ピクセルを生体組織解析装置１が取得して解析してもよい。
【００５８】
また、上述したステップＳ３からステップＳ６およびステップＳ９については必須の処理ではなく、必要に応じてまたはユーザの選択に応じて適宜実行される処理である。さらに、上述したステップＳ７およびステップＳ８は、要求される解析精度に応じていずれか一方または双方が実行されるものであり、ステップＳ７のみ、ステップＳ８のみ、ステップＳ７およびステップＳ８の双方の順に解析精度が向上する。
【００５９】
以上のような本実施形態の生体組織解析装置１、生体組織解析プログラム１ａおよび生体組織解析方法によれば、以下のような効果を奏する。
１．ハイパースペクトルデータを用いて生体組織の変質状態を高精度に識別することができる。
２．解析用データに対する光源の影響を除去し、解析精度を向上することができる。
３．解析用データに対する染色の影響を除去し、解析精度を向上することができる。
４．解析対象部位をピクセル単位で抽出することができる。
５．解析対象部位を細胞核単位で抽出することができる。
６．潰瘍やポリープにおける細胞核のガン化状態を病理診断に従って識別することができる。
７．波長帯域に可視光を含ませることにより、生体組織の変質状態を解析するのに適したハイパースペクトルデータを取得することができる。
【００６０】
つぎに、本発明に係る生体組織解析装置１、生体組織解析プログラム１ａおよび生体組織解析方法の具体的な実施例について説明する。
【実施例１】
【００６１】
本実施例１では、上述した本実施形態において、変質状態大別部５６による教師なし学習として採用されたクラスター分析の識別性能を確認する実験を行った。
【００６２】
具体的には、まず、生体組織として潰瘍性大腸炎に関連する検体を用意し、予め病理医によって、非ガン細胞、軽度異形成（ＬＧＤ）、高度異形成（ＨＧＤ）およびガン細胞の４ステージを識別させた。つぎに、識別された４ステージのそれぞれについて、細胞核のハイパースペクトルデータを取得した。
【００６３】
本実施例１において、ハイパースペクトルデータは、ハイパースペクトルカメラ１０（北海道衛星株式会社：ＨＳＣ１７０２）を正立顕微鏡（オリンパス株式会社）の接眼レンズに直接当てて取得した。このとき、顕微鏡に内蔵された白色ＬＥＤを光源とし、顕微鏡の倍率は４００倍（接眼１０倍、対物４０倍）であった。また、計測波長域は３５０ｎｍから１０５０ｎｍの可視光と近赤外域であり、波長サンプリング間隔を５ｎｍとして計１４１バンドのハイパースペクトルデータを取得した。
【００６４】
つぎに、取得したハイパースペクトルデータについてウォード法によるクラスター分析を行い、３〜４つの特徴量をクラスターとして抽出した。そして、各ステージについて抽出した全特徴量の識別率（類似度）を主成分分析と線形判別分析とによって算出し、各ステージの指標となる特徴量がないかどうかを検討した。その結果を図５に示す。
【００６５】
図５に示すように、非ガン細胞の各特徴量である非ガン１〜４は、他の特徴量に対する識別率がほぼ１００％であった。よって、非ガン細胞は、クラスター分析によって容易に識別しうることが示された。また、軽度異形成（ＬＧＤ）の特徴量の一つであるＬＧＤ１は、他の特徴量に対する識別がいずれも８３．３％以上であり、平均すると約９５．５％であった。よって、軽度異形成（ＬＧＤ）は、ガン細胞や高度異形成（ＨＧＤ）が持たない特有の特徴量を有しており、この特徴量を用いることにより軽度異形成（ＬＧＤ）は、クラスター分析によって他のステージから識別できることが示された。
【００６６】
以上の本実施例１によれば、教師なし学習としてのクラスター分析は、潰瘍のガン化状態における４ステージのうち、非ガン細胞および軽度異形成（ＬＧＤ）を高精度に識別できることが示された。
【実施例２】
【００６７】
本実施例２では、上述した本実施形態において、データ数削減部５５によるデータ数の削減方法として採用された主成分分析の性能を確認する実験を行った。
【００６８】
具体的には、まず、クラスター分析では識別が困難であったガン細胞と高度異形成（ＨＧＤ）のサンプルを用意し、実施例１と同様の方法によりハイパースペクトルデータを取得した。つぎに、主成分分析によってハイパースペクトルデータのデータ数を二次元データに削減し、寄与率を算出した。その結果を図６に示す。図６に示すように、第一主成分から第二主成分までの累積寄与率が９０％を超えているため、第三主成分までの次元を削減しても情報量の損失がきわめて少ないことが示された。
【００６９】
以上の本実施例２によれば、主成分分析は、ハイパースペクトルデータの情報量をほとんど損失させることなくデータ数を削減でき、解析速度を向上しうることが示された。
【実施例３】
【００７０】
本実施例３では、上述した本実施形態において、変質状態識別部５７による教師あり学習として採用されたロジスティック回帰の識別性能を確認する実験を行った。
【００７１】
具体的には、まず、実施例２でデータ数を削減したガン細胞と高度異形成（ＨＧＤ）のハイパースペクトルデータを用意した。つぎに、数値解析ソフト（Python）の拡張モジュール（scikit-learn）中にあるLogisticRegression関数を使用し、ロジスティック回帰によって解析を行った。このとき、高度異形成（ＨＧＤ）のサンプル数（４９０個）と、ガン細胞のサンプル数（３２０個）の比率を保ったまま、７：３の割合でトレーニングデータとテストデータに分割し、学習と予測（識別）を行ったところ、その識別率は９５．９％という高い値であった。
【００７２】
また、予測結果について、真陽性、偽陽性、偽陰性および真陰性の個数を表す混合行列を図７に示す。図７の混合行列に基づいて、予測結果の評価尺度の一つであるＦ値（精度と再現率の調和平均）を算出したところ、９６．９％という高いスコアであった。なお、行列内の各数値は、以下の値を示すものである。
左上（真陽性：ＴＰ）：ガンであり、ガンと識別された細胞
右上（偽陽性：ＦＰ）：ＨＧＤであるが、ガンと識別された細胞
左下（偽陰性：ＦＮ）：ガンであるが、ＨＧＤと識別された細胞
右下（真陰性：ＴＮ）：ＨＧＤであり、ＨＧＤと識別された細胞
【００７３】
以上の本実施例３によれば、教師あり学習としてのロジスティック回帰は、潰瘍のガン化状態における４ステージのうち、ガン細胞および高度異形成（ＨＧＤ）を高精度に識別できることが示された。
【実施例４】
【００７４】
本実施例４では、教師あり学習の一つであるサポートベクターマシンによって解析した場合の解析精度と、教師なし学習の一つであるクラスター分析によって予め解析対象を絞り込んでからサポートベクターマシンによって解析した場合の解析精度とを比較する実験を行った。
【００７５】
具体的には、潰瘍性大腸炎に関連する検体を用意し、サポートベクターマシンのみによって、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ）、高度異形成（ＨＧＤ）およびガン細胞の４ステージを識別した。また、同様の検体について、クラスター分析によって予め非ガン細胞を識別した後、サポートベクターマシンによって、残りの軽度異形成（ＬＧＤ）、高度異形成（ＨＧＤ）およびガン細胞の３ステージを識別した。その結果を図８に示す。
【００７６】
図８（ａ）に示すように、サポートベクターマシンのみによって解析した場合、ガン細胞、軽度異形成（ＬＧＤ）、高度異形成（ＨＧＤ）および非ガン細胞の４ステージについての識別率は、それぞれ９２．２％、７８．８％、５５．６％および１００．０％であった。一方、図８（ｂ）に示すように、クラスター分析とサポートベクターマシンとを組み合わせて解析した場合、ガン細胞、軽度異形成（ＬＧＤ）および高度異形成（ＨＧＤ）の４ステージについての識別率は、それぞれ９０．０％、７１．０％および８６．６％であった。
【００７７】
すなわち、クラスター分析とサポートベクターマシンとを組み合わせて解析した場合、軽度異形成（ＬＧＤ）の識別率は、７．８％向上し、高度異形成（ＨＧＤ）の識別率は１５．４％向上した。
【００７８】
以上の本実施例４によれば、潰瘍のガン化状態の識別は、変質状態識別部５７による教師あり学習のみによって行うよりも、変質状態大別部５６による教師なし学習と組み合わせた方が、識別精度が向上することが示された。
【実施例５】
【００７９】
本実施例５では、上述した本実施形態において、染色影響除去部５４によって除去すべき、解析対象部位の染色濃度によって変化するスペクトル成分を特定する実験を行った。
【００８０】
具体的には、まず、ガンの表層上皮性腫瘍でホルマリン固定したものからパラフィン包埋切片を作製し、ヘマトキシリン・エオジン染色を施した病理標本を用意した。なお、ヘマトキシリン・エオジン染色では、核内のクロマチンが濃青紫色に染色される。このため、細胞核内のクロマチン濃度によって光の透過率が変化し、染色液の吸収スペクトルに依存したスペクトル変化がみられることとなる。
【００８１】
つぎに、実施例１と同様のハイパースペクトルカメラ１０を使用し、倍率を２００倍に固定して１つの標本から平均３枚のハイパースペクトルイメージを撮影した。なお、ガン細胞では、細胞核が肥大化し細胞質を確認することが困難になるため、１枚から１０個の細胞核について分析し、１つの細胞核ごとの平均スペクトルを使用した。
【００８２】
つづいて、ハイパースペクトルイメージにおける染色濃度を目視で３段階に分け、染色濃度が低いものから、染色レベル１、染色レベル２、染色レベル３とした。そして、独立成分の数を３つとして、ハイパースペクトルデータの独立成分分析を行った。その結果を図９に示す。なお、分析に際しては、４４５ｎｍ〜７４５ｎｍの波長間を５ｎｍのバンドで区切って計測された反射強度を使用した。
【００８３】
図９に示すように、３つの独立成分（Signal1，Signal2，Signal3）のうち、特にSignal1で示す独立成分は、５４０ｎｍ〜７００ｎｍの波長範囲において段階的に変化していることが確認された。一方、その他の独立成分（Signal2，Signal3）は、全ての染色レベルにおいてほぼ共通していることが確認された。よって、Signal1で示す独立成分は、染色濃度によって変化するスペクトル成分といえる。
【００８４】
以上の本実施例５によれば、ハイパースペクトルデータを独立成分分析することにより、解析対象部位の染色濃度によって変化するスペクトル成分を抽出できることが示された。
【００８５】
なお、本発明に係る生体組織解析装置１、生体組織解析プログラム１ａおよび生体組織解析方法は、前述した実施形態に限定されるものではなく、適宜変更することができる。
【００８６】
例えば、上述した本実施形態では、変質状態大別部５６による教師なし学習として、クラスター分析が採用されているが、これに限定されるものではない。例えば、教師なし学習として主成分分析を採用し、当該主成分分析によって次元を削減したハイパースペクトルデータを線形判別分析によって解析してもよい。
【符号の説明】
【００８７】
１  生体組織解析装置
１ａ  生体組織解析プログラム
２  表示手段
３  入力手段
４  記憶手段
５  演算処理手段
１０  ハイパースペクトルカメラ
４１  プログラム記憶部
４２  教師データ記憶部
５１  ハイパースペクトルデータ取得部
５２  解析対象部位抽出部
５３  光源影響除去部
５４  染色影響除去部
５５  データ数削減部
５６  変質状態大別部
５７  変質状態識別部
【図１】
2019163981.tif 000003
【図２】
2019163981.tif 000004
【図３】
2019163981.tif 000005
【図４】
2019163981.tif 000006
【図５】
2019163981.tif 000007
【図６】
2019163981.tif 000008
【図７】
2019163981.tif 000009
【図８】
2019163981.tif 000010
【図９】
2019163981.tif 000011
【手続補正書】
【提出日】平成30年7月6日(2018.7.6)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析装置であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析装置；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを画像処理によって抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【請求項２】
前記ハイパースペクトルデータの取得時に用いた光源のスペクトル情報によって、前記解析対象部位に相当するピクセルの前記スペクトル情報を除算する光源影響除去部を有している、請求項１に記載の生体組織解析装置。
【請求項３】
前記解析対象部位に相当するピクセルの前記スペクトル情報に基づいて、前記解析対象部位の染色濃度によって変化するスペクトル成分を除去する染色影響除去部を有している、請求項１または請求項２に記載の生体組織解析装置。
【請求項４】
前記解析対象部位抽出部は、前記スペクトル情報の類似度に基づいて各ピクセルを解析対象部位とその他の部位とにラベリングし、前記解析対象部位としてラベリングされたピクセルをフィルタリングする、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項５】
前記解析対象部位抽出部は、前記生体組織の二次元画像を画像認識処理によって腺管を特定するとともに、前記腺管上における細胞核を解析対象部位として画像認識処理によって特定する、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項６】
前記変質状態大別部によって特定される前記変質状態は、潰瘍における、非ガン細胞および軽度異形成（ＬＧＤ：Low Grade Dysplasia）の２ステージであり、
前記変質状態識別部によって識別される前記変質状態は、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）およびガン細胞の４ステージ、またはポリープにおける、軽度異型腺腫（ＬＧＡ：Low Grade Adenoma）、高度異型腺腫（ＨＧＡ：High Grade Adenoma）およびガン細胞の３ステージである、請求項１から請求項５のいずれかに記載の生体組織解析装置。
【請求項７】
前記ハイパースペクトルデータの波長帯域に可視光が含まれている、請求項１から請求項６のいずれかに記載の生体組織解析装置。
【請求項８】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析プログラムであって、下記（i），（ii）としてコンピュータを機能させるとともに、下記（iii）および／または（iv）としてコンピュータを機能させる、前記生体組織解析プログラム；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを画像処理によって抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【手続補正書】
【提出日】平成30年7月6日(2018.7.6)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析装置であって、下記（i），（ii）を有するとともに、下記（iii）および／または（iv）を有する、前記生体組織解析装置；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを画像処理によって抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。
【請求項２】
前記ハイパースペクトルデータの取得時に用いた光源のスペクトル情報によって、前記解析対象部位に相当するピクセルの前記スペクトル情報を除算する光源影響除去部を有している、請求項１に記載の生体組織解析装置。
【請求項３】
前記解析対象部位に相当するピクセルの前記スペクトル情報に基づいて、前記解析対象部位の染色濃度によって変化するスペクトル成分を除去する染色影響除去部を有している、請求項１または請求項２に記載の生体組織解析装置。
【請求項４】
前記解析対象部位抽出部は、前記スペクトル情報の類似度に基づいて各ピクセルを解析対象部位とその他の部位とにラベリングし、前記解析対象部位としてラベリングされたピクセルをフィルタリングする、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項５】
前記解析対象部位抽出部は、前記生体組織の二次元画像を画像認識処理によって腺管を特定するとともに、前記腺管上における細胞核を解析対象部位として画像認識処理によって特定する、請求項１から請求項３のいずれかに記載の生体組織解析装置。
【請求項６】
前記変質状態大別部によって特定される前記変質状態は、潰瘍における、非ガン細胞および軽度異形成（ＬＧＤ：Low Grade Dysplasia）の２ステージであり、
前記変質状態識別部によって識別される前記変質状態は、潰瘍における、非ガン細胞、軽度異形成（ＬＧＤ：Low Grade Dysplasia）、高度異形成（ＨＧＤ：High Grade Dysplasia）およびガン細胞の４ステージ、またはポリープにおける、軽度異型腺腫（ＬＧＡ：Low Grade Adenoma）、高度異型腺腫（ＨＧＡ：High Grade Adenoma）およびガン細胞の３ステージである、請求項１から請求項５のいずれかに記載の生体組織解析装置。
【請求項７】
前記ハイパースペクトルデータの波長帯域に可視光が含まれている、請求項１から請求項６のいずれかに記載の生体組織解析装置。
【請求項８】
二次元画像を構成する各ピクセルごとにスペクトル情報を対応付けてなるハイパースペクトルデータを用いて生体組織を解析する生体組織解析プログラムであって、下記（i），（ii）としてコンピュータを機能させるとともに、下記（iii）および／または（iv）としてコンピュータを機能させる、前記生体組織解析プログラム；
（i）ハイパースペクトルカメラから生体組織の前記ハイパースペクトルデータを取得するハイパースペクトルデータ取得部；
（ii）前記ハイパースペクトルデータに基づいて、前記生体組織の二次元画像から解析対象部位に相当するピクセルを画像処理によって抽出する解析対象部位抽出部；
（iii）前記解析対象部位に相当するピクセルのそれぞれを教師なし学習によって分類することにより、前記生体組織の変質状態を大別する変質状態大別部；
（iv）前記解析対象部位に相当するピクセルのそれぞれを教師あり学習によって分類することにより、前記生体組織の変質状態を識別する変質状態識別部。