(57)【要約】【課題】正常組織からがん組織への遷移部分を検出する。【解決手段】観察装置１００は、対象物からの光の偏光状態を検出する偏光検出部４０と、偏光状態を示す値を主成分分析することにより、対象物の状態を判定する判定部２３ａとを備える。【選択図】図２
2020110291.tif 000002
【特許請求の範囲】
【請求項１】
対象物からの光の偏光状態を検出する偏光検出部と、
前記偏光状態を示す値を主成分分析することにより、前記対象物の状態を判定する判定部と、
を備えることを特徴とする観察装置。
【請求項２】
前記偏光状態を示す値は特定次元の値であり、
前記判定部は、主成分分析により前記特定次元の値をこの特定次元よりも低い所定次元の値に変換する
ことを特徴とする請求項１に記載の観察装置。
【請求項３】
前記偏光状態を示す値は前記特定次元の行列値であることを特徴とする請求項２に記載の観察装置。
【請求項４】
前記偏光状態を示す値はミュラー行列の行列値であることを特徴とする請求項３に記載の観察装置。
【請求項５】
前記判定部は、前記所定次元の値とこの所定次元よりも高い次元の値との差分を算出し、前記差分に基づいて前記対象物の状態を判定することを特徴とする請求項２〜４のいずれかに記載の観察装置。
【請求項６】
前記判定部は、前記差分と閾値との大小関係に基づいて前記対象物の状態を判定することを特徴とする請求項５に記載の観察装置。
【請求項７】
前記偏光検出部は、
それぞれ異なる偏光状態を有する光を生成する複数の偏光生成部と、
前記複数の偏光生成部が生成した前記光が照射された前記対象物からの光を検出する検出部と
を有することを特徴とする請求項１〜６のいずれかに記載の観察装置。
【請求項８】
対象物からの光の偏光状態を検出し、
前記偏光状態を示す値を主成分分析することにより、前記対象物の状態を判定する
ことを特徴とする観察方法。
【請求項９】
被検体の体腔内に挿入される管部材を有する挿入部と、
前記挿入部を操作する操作部とを備え、
前記請求項１〜７のいずれかに記載の前記観察装置の前記偏光検出部が前記挿入部に配置されていることを特徴とする内視鏡装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、観察装置、観察方法及び内視鏡装置に関するものである。
【背景技術】
【０００２】
偏光した光を使用し、消化器、特に、胃壁の偏光異方性を持つ粘膜層から戻る戻り光の非偏光の光の割合、すなわち戻り光の偏光度に基づいて粘膜層の厚みを算出することにより、胃壁の粘膜層の厚さの変化を検出し、これにより、がんの浸潤度を診断できる可能性があることが知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第５５０１１５５号公報
【０００４】
しかしながら、上述の特許文献１に開示された技術では、非がん組織、すなわち正常組織からがん組織への遷移過程の状態にある遷移部分を検出することが困難であった。
【発明の概要】
【０００５】
第１の態様によれば、観察装置は、対象物からの光の偏光状態を検出する偏光検出部と、偏光状態を示す値を主成分分析することにより、対象物の状態を判定する判定部とを備える。
第２の態様によれば、観察方法は、対象物からの光の偏光状態を検出し、偏光状態を示す値を主成分分析することにより、対象物の状態を判定する。
【図面の簡単な説明】
【０００６】
【図１】実施の形態である観察装置が適用される内視鏡装置の要部を示す斜視図である。
【図２】実施の形態である観察装置が適用される内視鏡システムの概略構成を示すブロック図である。
【図３】実施の形態である観察装置の要部を示す断面図である。
【図４】実施の形態である観察装置の要部を示す斜視図である。
【図５】実施の形態である観察装置の要部を示す斜視図である。
【図６】実施の形態に係る内視鏡装置の挿入部の先端部の概略構成を示す図である。
【図７】実施の形態である観察装置における、偏光生成部から出射する光の偏光状態の関係を示す図である。
【図８】正常組織、がん組織及び遷移部分の関係を示す図である。
【図９】実施の形態である観察装置による主成分分析の手法を説明するための図である。
【図１０】実施の形態である観察装置による主成分分析の手法を説明するための図である。
【図１１】実施の形態である観察装置による主成分分析の手法を説明するための図である。
【図１２】実施の形態である観察装置の動作を説明するためのフローチャートである。
【図１３】主成分分析した結果である第２主成分、第３主成分及び第４主成分で表現される３次元空間をメルカトル図法により展開した図である。
【発明を実施するための形態】
【０００７】
図１及び図２を参照して、実施の形態である観察装置が適用される内視鏡装置及び内視鏡システムについて説明する。
【０００８】
図１は実施の形態である観察装置が適用される内視鏡装置の要部を示す斜視図であり、図２は実施の形態である観察装置が適用される内視鏡システムの概略構成を示すブロック図である。
【０００９】
これらの図において、本実施の形態に係る内視鏡システム１は、内視鏡装置１０および本体部２０を有する。
【００１０】
内視鏡装置１０は、人体等の被検体の体腔内に挿入部１１の先端部１２が挿入され、被検体内の粘膜内面等の対象物に光を照射する。また、内視鏡装置１０は、対象物からの戻り光を撮像し、対象物像の画像信号を取得する。
【００１１】
本体部２０は、内視鏡装置１０の先端部１２から照射される光を発生する光源２６を有し、内視鏡装置１０が取得した画像信号に対して各種画像処理を行う。
【００１２】
図１に示すように、内視鏡装置１０は、挿入部１１、操作部１３、および、コネクタ部１４を有する。挿入部１１は、可撓性（つまり軟性）を有する管部材１１ａを有し、細長形状を有する。操作部１３は、挿入部１１の先端部１２が行う湾曲操作や観察操作を操作する。コネクタ部１４は、内視鏡装置１０を本体部２０に着脱自在に接続する。
【００１３】
図２に示すように、本体部２０は、制御部２１と、記憶部２２と、画像処理部２３と、検出部ドライバ２４と、光源ドライバ２５と、光源２６と、集光レンズ２７と、ディスプレイコントローラ２８と、入力インタフェース（Ｉ／Ｆ）２９とを有する。本体部２０は、光学コネクタ３０及び電気コネクタ３１により内視鏡装置１０に接続されている。
【００１４】
制御部２１はＣＰＵ等の演算素子を有し、記憶部２２に格納されている図略の制御用プログラムが起動時に読み出されてこの制御部２１において実行される。これにより、制御部２１は、画像処理部２３、検出部ドライバ２４、光源ドライバ２５、ディスプレイコントローラ２８及び入力インタフェース２９を含む内視鏡システム１全体の制御を行う。
【００１５】
記憶部２２はハードディスクドライブ等の大容量記憶媒体及びＲＯＭ、ＲＡＭ等の半導体記憶媒体を備える。この記憶部２２には上述の制御用プログラムが格納されているとともに、制御部２１の制御動作時に必要とされる各種データが一時的に格納される。また、記憶部２２には、画像処理部２３によって処理がされた画像データが格納される。
【００１６】
画像処理部２３は、内視鏡装置１０の検出部５０により撮像された対象物像の画像信号を取得し、この画像信号に対して各種画像処理を行う。画像処理部２３は、画像処理結果をディスプレイコントローラ２８に出力する。
【００１７】
画像処理部２３は判定部２３ａを有する。判定部２３ａは、検出部５０により撮像された対象物像の画像信号に基づいて、対象物からの光の偏光状態を示す値を主成分分析する。そして、判定部２３ａは、この主成分分析の結果に基づいて対象物の状態を判定する。
【００１８】
検出部５０及び画像処理部２３により行われる画像処理の詳細については後述する。
【００１９】
検出部ドライバ２４は、検出部５０に対して、この検出部５０を駆動する信号を発生して検出部５０に供給する。
【００２０】
光源２６は内視鏡装置１０の先端部１２から照射される光を発生するものであり、光源ドライバ２５により駆動される。本実施の形態に係る光源２６は狭帯域でかつ空間的コヒーレンスが低く、スペックルが除去された光を発生する。光源２６が発生する光の波長帯域は目的に応じて３５０ｎｍから８５０ｎｍの間で１波長または複数波長が選択されうる。
【００２１】
但し、生体組織に含まれるヘモグロビンの光吸収スペクトルが４００ｎｍ付近で特有の山を有することが知られている。このため、光源２６からの光は４５０ｎｍ以上の波長帯域の光であることが好ましい。
【００２２】
光源２６が発生する光の帯域幅は５ｎｍ以下、好ましくは３ｎｍ程度とされる。帯域幅の関係で、このような光源２６はレーザー光源から構成されることが好ましい。一方、ＬＥＤ光源にバンドパスフィルタを付加したものであっても同様に狭帯域な光を発生することができる。
【００２３】
光源２６が発生する光の波長帯域及び帯域幅は、検出部５０のイメージセンサ５７上に発生させる後述する干渉縞の本数及び幅によって定められうる。
【００２４】
光源２６により発生された光は、集光レンズ２７により一定幅の光束に収束され、光学コネクタ３０を介して内視鏡装置１０の光ファイバ１５内に導光される。
【００２５】
ここで、本実施の形態に係る内視鏡システム１の本体部２０には光源２６が１つのみ設けられる。一方、後述するように、光ファイバ１５の本数は、後述する偏光検出部４０の偏光生成部６０の個数に等しい本数である。本実施の形態では、偏光生成部６０の個数は６つであるため、光ファイバ１５の本数は６本である。このため、光源２６と集光レンズ２７との間には切替部２６ａが設けられている。切替部２６ａは、光源２６からの光をいずれか一本の光ファイバ１５内に導光する。また、切替部２６ａは、光を導光する光ファイバ１５以外の光ファイバ１５への光を遮断する。
【００２６】
ディスプレイコントローラ２８は、画像処理部２３からの出力である画像信号を、本体部２０の外に設けられたディスプレイ３２の画面に表示させるためのディスプレイ駆動信号を発生し、これをディスプレイ３２に供給する。
【００２７】
入力インタフェース２９は、本体部２０の外に設けられたキーボード等の入力装置３３からの操作入力信号を受け入れ、これを制御部２１に供給する。
【００２８】
そして、本実施の形態である観察装置の偏光検出部４０は、図３及び図６に示すように、内視鏡装置１０の挿入部１１の先端部１２に設けられている。
【００２９】
以下、図３〜図７を参照して偏光検出部４０について説明する。なお、以下の説明では主に光学素子及び光路について説明を行う関係で、図３〜図７においては検出部５０と本体部２０とを結ぶ電気回路や内視鏡装置１０の被覆等他の構成についての詳細な説明は省略する。
【００３０】
偏光検出部４０は、対象物からの光を検出する検出部５０と、この検出部５０の周囲に配置され、それぞれ異なる偏光状態を有する光を生成して対象物に照射する６つの偏光生成部６０とを有する。これら検出部５０及び偏光生成部６０は、図３及び図６に示すように、内視鏡装置１０の挿入部１１の管部材１１ａ内に収納されている。
【００３１】
なお、図３以降において図示するように、検出部５０の光軸方向にｘ軸を取っている。そして、ｘ軸の正方向を対象物に向かう方向とする。
【００３２】
検出部５０は、図３に示すように、前群レンズ５１、絞り５２、偏光分離部５３、結像レンズ５４、偏光板５５、保護ガラス５６及びイメージセンサ５７を有する。
【００３３】
前群レンズ５１は、検出部５０において対象物側（図３において右端）に配置されている。絞り５２は円環状に形成され、前群レンズ５１の出射側に配置されている。偏光分離部５３には絞り５２から出射される光が入射される。結像レンズ５４は、偏光分離部５３から出射される光の干渉縞を後述するイメージセンサ５７の撮像面５７ａ上に結像させる後群レンズである。検光子である偏光板５５及び保護ガラス５６は、イメージセンサ５７の撮像面５７ａの前方（図６において右方）に順に配置されている。イメージセンサ５７は、複数の画素を有する撮像素子であり、その撮像面５７ａに結像された干渉縞を撮像する。
【００３４】
検出部５０は、後に詳述する６つの偏光生成部６０による対象物の照明範囲の共通部分の少なくとも一部をその検出範囲とするように設計されており、より好ましくは、６つの偏光生成部６０による全照明範囲よりも撮像範囲が広くなるように設定されている。あるいは、検出部５０の検出範囲に応じて、この検出範囲を偏光生成部６０による対象物の照明範囲の共通部分の少なくとも一部とするように偏光生成部６０を設計してもよい。
【００３５】
検出部５０の検出範囲に生成した偏光が照射されるように、複数の偏光生成部６０のそれぞれの出射角度を設定してもよい。この場合、複数の偏光生成部６０のそれぞれの配置状態を調整してもよい。このように偏光生成部６０による対象物の照明範囲と検出部５０の検出範囲とを設定することで、対象物の照射範囲からの戻り光を検出部５０により確実に検出することができる。
【００３６】
偏光分離部５３は、２枚のサバール板５３ｃと、これらサバール板５３ｃの間に配置された１／２波長板５３ｄとを有する。本実施の形態におけるサバール板５３ｃは、複屈折性を有する一軸性結晶からなる２枚の平行平面板５３ａ（例えばＹＶＯ４）と、これら平行平面板５３ａの間に配置された１／２波長板５３ｂとを有する。
【００３７】
サバール板５３ｃは、複屈折性を有する一対の平行平面板５３ａを、その光軸が９０°異なるように貼り合わせたものである。サバール板５３ｃに入射する光が異なる偏光状態が重なり合った光であった場合、サバール板５３ｃは異なる偏光状態の光を分離して出射させる。
【００３８】
偏光分離部５３から出射された光は結像レンズ５４により収束され、偏光板５５及び保護ガラス５６を通してイメージセンサ５７の撮像面５７ａ上で結像される。そして、異なる偏光状態の光がこの撮像面５７ａ上で干渉して干渉縞を形成する。
【００３９】
イメージセンサ５７は、その撮像面５７ａ上に結像した光を撮像し、その結果を画像信号として出力する。イメージセンサ５７は、上述のように撮像面５７ａ上に形成される干渉縞を撮像するため、微細な干渉縞が撮像可能な解像度を有している。撮像面５７ａ上に形成される干渉縞の本数及び幅は、既に説明したように、光源２６が発生する光の波長帯域及び帯域幅に依存する。
【００４０】
複数の偏光生成部６０はそれぞれ、図４及び図５に示すように、１枚の偏光板（偏光変換素子）６１ａ〜６１ｆと、この偏光板６１ａ〜６１ｆの出射側に配置された平凹レンズ（発散光学素子）６２とを有する。また、一部の偏光板６１ａ、６１ｂと平凹レンズ６２との間には１／４波長板（偏光光学素子）６３が配置されている。１／４波長板６３が配置されている場合、１枚の偏光板６１ａ、６１ｂ、１／４波長板６３及び平凹レンズ６２により偏光生成部６０が構成される。なお、以下の説明及び図示において、偏光板６１ａ〜６１ｆの位置を特定せずに一般的に示す場合は、符号６１により代表して説明する。
【００４１】
図４及び図５に示すように、６つの偏光生成部６０は二群に分けられて検出部５０の周囲に配置されている。より詳細には、３つの偏光生成部６０が偏光分離部５３の図中上方に配置され、３つの偏光生成部６０が偏光分離部５３の図中下方に配置されている。また、６つの偏光生成部６０は、検出部５０の光軸方向であるｘ軸に直交する同一面内に配置されている。
【００４２】
また、図４及び図５に示すように、ｘ軸に沿って見た場合、複数の偏光生成部６０のそれぞれの一部がイメージセンサ５７の一部と重複している。これにより、検出部５０の光軸方向から見た本実施の形態の偏光検出部４０の外径をよりコンパクトにすることができる。この結果、偏光検出部４０全体の小型化に寄与することができる。複数の偏光生成部６０のうちの少なくとも一つの少なくとも一部が、検出部５０の光軸方向にイメージセンサ５７に重複するように配置されていればよい。
【００４３】
６つの偏光生成部６０は、それぞれ異なる偏光状態を有する光を対象物に照射する。この点について以下詳細に説明する。
【００４４】
光の偏光状態はストークスベクトルで記述することができ、このストークスベクトルは次式のように、４行１列のストークス行列で書き表すことができる。
【数１】
2020110291.tif 000003
ここに、ストークス行列の各成分は、ｓ０が光強度、ｓ１がｘ−ｙ直線偏光、ｓ２が４５°直線偏光、ｓ３が円偏光である。ｓ０は強度を表すため正の値を取る。ｓ１＝ｓ０とは０°直線偏光であることを意味し、ｓ１＝−ｓ０は９０°直線偏光であることを意味する。また、ｓ２＝ｓ０は４５°直線偏光であることを意味し、ｓ２＝−ｓ０は−４５°直線偏光であることを意味する。さらに、ｓ３＝ｓ０は右回り円偏光であることを意味し、ｓ３＝−ｓ０は左回り円偏光であることを意味する。
【００４５】
それぞれの偏光生成部６０、より詳細には偏光板６１ａ〜６１ｆ及び１／４波長板６３のそれぞれの組み合わせは、ｓ１＝±１、ｓ２＝±１、ｓ３＝±１の６種類の偏光状態のいずれか一つの偏光状態を有する光を生成するように設定されている。
【００４６】
具体的には、６つの偏光生成部６０では、図７に示すように、各々の偏光生成部６０から出射される光のストークス成分がそれぞれ図示する成分となるように、偏光板６１ａ〜６１ｆが選択され、また、１／４波長板６３の配置位置が定められている。
【００４７】
各々の偏光生成部６０から出射される光の偏光状態をどのように異ならせるか、言い換えればどのようなストークス成分を有する光を出射させるかは任意に決定可能である。
【００４８】
ここで、ストークス成分ｓ１は角度による偏光状態の変化が他の成分ｓ２、ｓ３よりも小さいので、ストークス成分ｓ１＝＋１を有する光を照射する偏光板６１ｃ及びストークス成分ｓ１＝−１を有する光を照射する偏光板６１ｆが離れて配置されても問題ない。
【００４９】
加えて、ストークス成分ｓ３＝±１を有する光を照射する偏光板６１ａ、６１ｂ及び１／４波長板６３については、これら偏光板６１ａ、６１ｂの出射側に配置される１／４波長板６３を共通化でき、また、１／４波長板６３を保持する図略のホルダを共通化できる観点から、これら偏光板６１ａ、６１ｂを隣り合わせに配置するメリットがある。
【００５０】
レンズ６２は、偏光板６１及び１／４波長板６３から出射される光を発散し、対象物の所定の領域を照明する平凹レンズである。ここで、平凹レンズ６２から出射される光の半画角が１５度以内となるように、平凹レンズ６２の光学設計がなされている。これは、半画角が１５度以内であれば、偏光生成部６０から照射された光の偏光状態は大きく変化しないからである。また、後に詳述する光ファイバ１５の開口数（ＮＡ：numerical aperture）を考慮しても、この光ファイバ１５から出射する光の半画角も１５度以内とされる。
【００５１】
偏光生成部６０に用いられる発散光学素子としての平凹レンズ６２は、同様に光の発散効果を有する他の光学素子、例えばボールレンズと比較して、偏光性能が良いことから、偏光生成部６０として好適である。
【００５２】
偏光生成部６０には、図３及び図６で詳細に示すように、内視鏡装置１０の光ファイバ１５により光が導光される。光ファイバ１５は、図３に示すように、内視鏡装置１０の挿入部１１の管部材１１ａ内を挿入部１１の長手方向に沿って延び、この挿入部１１の先端部１２にまで至ってその出射端１５ａが偏光板６１の入射端（図３において左端）近傍に配置されている。
【００５３】
光ファイバ１５は、光源２６からの光の偏光状態を解消して（偏光スクランブルして）偏光生成部６０まで導光する。本実施の形態では、図６に示すように、光ファイバ１５はコア１５ｂの断面が略正方形でクラッド１５ｃの断面外形が略円形のマルチモード光ファイバである。このような構成の光ファイバ１５は、偏光板６１の入射面に偏光解消された安定な光を供給でき、したがって偏光板６１を通過した偏光状態は空間的及び時間的に安定な光強度を保つことができる。
【００５４】
本実施の形態では、内視鏡装置１０に偏光生成部６０と同数、つまり６本の光ファイバ１５が設けられており（図３では２本のみ図示している）、これら光ファイバ１５には単一の光源２６から共通に光が導光される。従って、偏光生成部６０からは、時分割的に偏光状態の異なる光が対象物に照射されることになる。
【００５５】
なお、本実施の形態である偏光検出部４０を構成する各要素は、高温環境下における滅菌処理を行っても高い偏光計測精度が維持できるような材料で形成されている。具体的には、１／２波長板５３ｂ、５３ｄ及び１／４波長板６３は水晶から構成されている。また、平凹レンズ６２は光学ガラスから構成されている。さらに、偏光板６１は無機材料またはワイヤーグリッドから構成されている。そして、サバール板５３ｃを構成する平行平面板５３ａはＹＶＯ４から構成されている。
【００５６】
ここに、滅菌処理のオートクレーブは、一例として１１５℃で３０分間、１２１℃で２０分間、１２６℃で１５分間、１３４℃で１０分間のいずれかの条件で行われる。
【００５７】
そして、本実施の形態である観察装置１００は、偏光検出部４０と画像処理部２３、特に判定部２３ａとから構成される。
【００５８】
次に、偏光検出部４０を用いた内視鏡システム１によるストークス成分ｓ０〜ｓ３の測定方法の原理について説明する（K. Oka and N. Saito, "Snapshot complete imaging polarimeter using Savart plates", SPIE 6295-9, 1 (2006)参照）。
【００５９】
イメージセンサ５７により撮像された光の強度分布をＩ（ｘ，ｙ）とする。測定光に含まれるストークス成分の２次元分布をそれぞれｓ０（ｘ，ｙ）、ｓ１（ｘ，ｙ）、ｓ２（ｘ，ｙ）、ｓ３（ｘ，ｙ）とすると、これらストークス成分の２次元分布を用いて光強度分布を表すと次式のようになる。
【数２】
2020110291.tif 000004
ここに、ａｒｇは複素数の偏角を示す関数であり、Ｕ１及びＵ２は、それぞれ２枚のサバール板５３ｃにより導入される空間キャリア周波数である。
【００６０】
上式におけるｓ０（ｘ，ｙ）、ｓ２（ｘ，ｙ）及びｓ１３（ｘ，ｙ）は、これら（特にｓ１３（ｘ，ｙ）については実数成分及び虚数成分）がそれぞれ異なる４つの空間キャリア周波数ｆｙ＝０、Ｕ２、Ｕ２−Ｕ１、Ｕ２＋Ｕ１を有することから、光強度分布Ｉ（ｘ，ｙ）を空間周波数フィルタリングすることにより得られる。そして、これらストークス成分の２次元分布は、抽出された成分の振幅及び位相から得ることができる。この際、空間周波数フィルタリング及び振幅、位相の変調は、ストークス成分の２次元分布の変調に適した形にされたフーリエ変換技術により一度に行うことができる。
【００６１】
そこで、内視鏡システム１の画像処理部２３（含む判定部２３ａ）は、検出部５０から出力される画像信号の強度の２次元分布を取得し、この強度の２次元分布をフーリエ変換することで、対象物からの戻り光の偏光状態、具体的にはストークス成分ｓ０〜ｓ３の２次元分布を得る。
【００６２】
そして、画像処理部２３の判定部２３ａは、このストークス成分ｓ０〜ｓ３の２次元分布を用いて、対象物の偏光度を求める。そして、判定部２３ａは、この偏光度に基づいて対象物の粘膜層の厚みを算出し、これにより、がんの浸潤度を診断する。
【００６３】
対象物への入射光のストークス行列をＳ＝（ｓ０，ｓ１，ｓ２，ｓ３）とし、この対象物からの戻り光のストークス行列をＳ′＝（ｓ′０，ｓ′１，ｓ′２，ｓ′３）とすると、これらストークス行列の関係は４行４列のミュラー行列Ｍにより表される。すなわち、
【数３】
2020110291.tif 000005
【００６４】
ここで、ミュラー行列Ｍの全１６の要素ｍ００〜ｍ３３の各要素と偏光の物理的特性との厳密な対応は難しいが、おおまかな関係としては、要素ｍ００は輝度を表し、全１６の要素ｍ００〜ｍ３３は偏光度を表し、要素ｍ０１、ｍ０２、ｍ１０及びｍ２０は二色性（直線複吸収）を表し、要素ｍ０３及びｍ３０は円二色性（円複吸収）を表し、要素ｍ１１、ｍ１２、ｍ２１及びｍ２２は旋光性（円複屈折）を表し、要素ｍ１１〜ｍ１３、ｍ２１〜ｍ２３及びｍ３１〜ｍ３３は複屈折性（直線複屈折）を表すものである。
【００６５】
既に説明したように、本実施の形態である偏光検出部４０によれば、６種類の互いに偏光状態の異なる光、言い換えれば、互いに異なるストークス成分（このストークス成分は全て既知である）を有するストークス行列（これをＳ１〜Ｓ６と置く）により表される光を対象物に入射し、それぞれの光が反射して得られる対象物からの戻り光のストークス成分を検出することができる（このストークス成分からなるストークス行列をＳ′１〜Ｓ′６と置く）。そして、対象物のミュラー行列Ｍは同一であるので、これらストークス行列Ｓ１〜Ｓ６、Ｓ′１〜Ｓ′６からミュラー行列Ｍの各成分を求めることができる。
【００６６】
すなわち、本実施の形態である偏光検出部４０において、６つの偏光生成部６０のうち４つを用いた場合、４つの行列式を解くことでミュラー行列Ｍの各成分は一意に定まる。また、６つの偏光生成部６０を全て用いれば、最小二乗法を用いてミュラー行列Ｍの各成分を求めることができる。
【００６７】
そして、戻り光のストークス行列Ｓ′＝（ｓ′０，ｓ′１，ｓ′２，ｓ′３）のストークス成分の２次元分布がわかっているので、対象物のミュラー行列Ｍの各成分の２次元分布も求めることができる。これにより、対象物の偏光度の２次元分布を求めることができる。そして、この偏光度の２次元分布に基づいて対象物の粘膜層の厚みの２次元分布を算出し、これにより、対象物の状態の１つであるがんの浸潤度を診断することができる。
【００６８】
加えて、本実施の形態である観察装置１００を構成する判定部２３ａによれば、上述した原理に基づいて正常組織及びがん組織を判別するのみならず、正常組織からがん組織に遷移している状態にある遷移部分を検出することができる。以下、この検出原理について説明する。
【００６９】
図８は、対象物として消化管を例に取った、がん組織、正常組織及び遷移部分の関係を示す図である。図８（ａ）、（ｂ）はいずれも消化管を長さ方向に切開して展開した状態を示す。
【００７０】
図８（ａ）に示すように、消化管２００の内面にはがん組織が存在する領域２０１があるものとする。医師は、がん組織を目視で確認し、あるいは、消化管２００の内面の生体組織をサンプルとして切り取り、病理検査を行った上で確定診断を行う。図中、がん組織のサンプル採取箇所を符号２０２で示す。
【００７１】
がん組織は、がん細胞の成長に伴って周辺に拡大する。医師は、がん組織が存在する領域２０１を確認したら、領域２０１の拡大を阻止し、さらに消化管以外への転移を未然に防ぐ必要がある。そこで、医師は、領域２０１の周囲（図中左右方向）に存在する正常組織を含めた切除領域２０３を設定し、この切除領域２０３の消化管を切除する。
【００７２】
図８（ａ）において、正常組織であると診断されたサンプル採取箇所を符号２０４で、遷移部分であると診断されたサンプル採取箇所を符号２０５で示す。
【００７３】
正常組織及びがん組織であるか否かの判定は、上述したように、ミュラー行列Ｍの１６要素の値、言い換えれば各成分に基づいて行うことができる。一方、遷移部分であるか否かの判定は、ミュラー行列Ｍの１６要素の値から直接判定することができなかった。
【００７４】
そこで、発明者らは、正常組織、がん組織及び遷移部分から実際に算出したミュラー行列Ｍの１６要素の値に基づいて種々解析を行った結果、所定の手順に基づけば遷移部分か否かの判定が可能であるとの知見に至った。
【００７５】
より詳細には、対象物のうち所定の領域に、６つの偏光生成部６０からそれぞれ偏光を照射し、その照射領域からの光をイメージセンサ５７で受光する。判定部２３ａは、イメージセンサ５７の複数の画像（本実施例では６つの画像）を用いて画素ごとに、ミュラー行列Ｍを算出する。これにより、対象物の照射領域からの光を受けた画素の数だけミュラー行列Ｍが取得される。判定部２３ａは、取得したミュラー行列Ｍのそれぞれの１６個の要素を１６次のデータとみなす。これにより、光を受けた画素の数と同数の１６次のデータが得られる。判定部２３ａは、この１６次のデータ群を用いて主成分分析することにより、１６個の主成分ベクトルの線形結合で各１６次元データ（ミュラー行列）を表す。
【００７６】
以下、判定部２３ａが行う主成分分析の概要について説明する。なお、主成分分析自体は公知の手法であるので、概要のみ説明する。
【００７７】
主成分分析とは、多次元（ｎ次元）のデータ
【数４】
2020110291.tif 000006
があったとき各成分の分散共分散行列
【数５】
2020110291.tif 000007
【数６】
2020110291.tif 000008
が対角行列（無相関）になるように座標変換を行う多変量解析の方法である。
【００７８】
適当な直交行列
【数７】
2020110291.tif 000009
でｘを
【数８】
2020110291.tif 000010
によって新しい変数に変換する。このとき、各データは
【数９】
2020110291.tif 000011
とあらわされている。あたらしい変数ｙの分散共分散行列
【数１０】
2020110291.tif 000012
【数１１】
2020110291.tif 000013
が対角行列になるような直交行列Ｐを見つけることが主成分分析である。このとき、分散共分散行列の性質（対称行列）により対角成分（ε１１, …，εｎｎ）（分散共分散行列の固有値である）は０または正の数となる。
【数１２】
2020110291.tif 000014
を分散共分散行列の固有値の大きな方から第１主成分(ベクトル)、…、第ｎ主成分(ベクトル)と呼ぶ。
【００７９】
新しい変数
【数１３】
2020110291.tif 000015
を主成分（ベクトル）のスコアと呼ぶ。
【００８０】
分散共分散行列の固有値は新しい変数の分散なので、値が小さいということは全データの対する寄与が小さいことを意味する。従って、他に比べて小さなスコアは無視してもデータの持つ情報はほとんど失われない。
【００８１】
計測データの場合、固有の値が測定誤差の程度であれば、その変数はノイズとみなすことができる。
【００８２】
そこで、情報の欠落が少ない適当な次元ｍ＜ｎを決めて、
【数１４】
2020110291.tif 000016
によって近似すれば、ｎ次元のデータをｍ次元に次元圧縮しても、もとのデータの情報はほとんど失われない。
【００８３】
なお、主成分分析で圧縮したミュラー行列に、光学的意味があることを確認することが必要である。確認の手法には、後述する手法が用いられる。また、次元圧縮の対象として、ミュラー行列から変換できる物理量、例えば分散共分散（Covariant）行列もしくはコヒーレント（Coherent）行列、あるいはその行列要素の絶対値からなる行列、なども有効である。
【００８４】
そこで、本実施の形態である観察装置１００の判定部２３ａは、データ（ミュラー行列）を主成分分析により得られた１６個の主成分ベクトルの一部を使って表現することにより、１６次元よりも低次元の空間に射影する。すなわち、主成分分析により１６次元よりも低次元の特定の次元の値に変換する。これにより、元々ミュラー行列の次元である１６次元の値をより低次元の値に圧縮する。本実施の形態では、上述の説明においてｎ＝１６であり、ｍ＝４であり、つまり、１６次元の値を４次元の値に圧縮する。次いで、判定部２３ａは、上記の圧縮に用いなかった主成分ベクトル（本実施例では第５〜１６主成分ベクトル）を使ってデータを表現する。すなわち、判定部２３ａは、１〜１６次元の値から１〜４次元の値を除いた値である５次元以上の値（以下、これを残差と称する）を算出する。ここで、１〜４次元の値とは、数式９の右辺の第１主成分を含む第１項から第４主成分を含む第４項までを足し合わせた値である。そして、判定部２３ａは、この残差を予め定めた閾値と比較し、残差が閾値以上になった場合、その残差の算出に用いられたミュラー行列Ｍの１６要素が得られた領域が遷移部分であると判定する。
【００８５】
尚、判定部２３ａによる遷移部分か否かの判定は、対象物の照射領域全体が遷移部分か否かを判定するのに代えて、照射領域からの光を受けた複数の画素ごとに残差を求め、画素ごとに残差と閾値との比較を行い、照射領域のうち、残差が閾値以上になった画素に対応する部分が遷移部分であると判定してもよい。
【００８６】
この４次元の値には、がん組織であるか正常組織であるかを判定するための偏光情報が平均で９５％以上が含まれることが、発明者らの実験結果により明らかになった。一方、１〜ｎ次元（ｎ＞４）の値から１〜４次元の値を除いた残差を正常組織、がん組織及び遷移部分とで比較すると、遷移部分についてのみこの残差が大きいことも、発明者らの実験結果により明らかになった。従って、主成分分析の結果得られた１〜４次元の値に基づき残差を算出し、この残差が閾値以上である組織を遷移部分であると判定することが可能になった。
【００８７】
尚、遷移部分であるか否かの判定に用いられる閾値は、例えば、以下のように設定してもよい。イメージセンサ５７の複数の画素ごとに算出したミュラー行列Ｍの各要素のヒストグラムを取る。そして、このヒストグラムの最大幅の、正常（非ガン）のときのヒストグラムの最大幅に対する倍数を、閾値とする。倍数は、例えば２倍もしくは３倍である。ヒストグラムは、ミュラー行列Ｍの複数の要素のそれぞれについて作成され、要素の値を横軸とし、頻度を縦軸としたヒストグラムであり、各要素において、主成分分析に用いた画素の数と同数もしくは特定の画像に含まれる画素数と同数の値が用いられる。
【００８８】
上述した手法によれば、判定部２３ａを含む画像処理部２３はミュラー行列Ｍの１６要素の２次元分布を得ることができるので、判定部２３ａは、遷移部分の２次元分布も得ることができる。これにより、図８（ｂ）に示すように、遷移部分が存在する遷移領域２０６をイメージとして判定し、検出することができる。
【００８９】
本実施の形態である観察装置１００の判定部２３ａの作用について、図９〜図１１を用いて説明する。但し、４次元の値を３次元で表現することは難しいので、３次元の値を主成分分析により２次元の値に圧縮した場合について説明する。
【００９０】
図９は、ミュラー行列Ｍの１６要素のうち、ｍ１１、ｍ２２、ｍ３３の３要素の成分の分布を示す図である。なお、図９〜図１１に示す例は仮想的なものであり、実際の値とは異なる。図９において成分３００は図示するような分布を有するものとする。また、図中に第１主成分、第２主成分及び第３主成分の軸も図示している。
【００９１】
図９に示す分布を、それぞれ直交する第１主成分、第２主成分及び第３主成分の軸で書き直すと図１０（ａ）のようになる。図１０（ａ）は立体的（３次元的）な図示であるので、第１主成分の軸と第２主成分の軸とを含む平面が図の平面に一致するような平面的（２次元的）な図示に書き換えると図１０（ｂ）のようになる。さらに、第１主成分の軸と第３主成分の軸とを含む平面が図の平面に一致するような平面的な図示に書き換えると図１０（ｃ）のようになる。
【００９２】
図１０（ｃ）は、いわば主成分分析により３次元から２次元にその次元を圧縮した表現であると考えられる。従って、図中矢印で挟まれた成分が２次元以上の値の残差に相当する。
【００９３】
発明者らの解析によれば、図１１（ａ）に示すように、正常組織及びがん組織についての２次元以上の残差（図中では２次元残差と表示している）は小さい。一方、図１１（ｂ）に示すように、遷移部分についての２次元残差は大きい。そこで、判定部２３ａは、図１１（ａ）に示す２次元残差と図１１（ｂ）に示す２次元残差とを峻別しうる閾値を設定し、２次元残差が閾値以上であるか否かに基づいて遷移部分を判定する。
【００９４】
上述した手法により遷移部分を判定できる理由について、発明者らは以下のように推測する。
【００９５】
正常組織を構成する細胞の核の大きさはほぼ一定である。また、がん組織を構成する細胞の核の大きさは、正常組織の核より大きいものの、これもほぼ一定である。従って、正常組織及びがん組織はそれぞれ特定の偏光状態を有する偏光子として作用すると考えられる。そして、正常組織とがん組織とではその偏光状態が異なるので、これらをミュラー行列Ｍの１６要素の値に基づいて判定することができる。
【００９６】
一方、遷移部分においては正常組織とみなせる細胞とがん組織とみなせる細胞とが混在していると考えられる。従って、遷移部分に光を当てると複雑な散乱光が発生すると推測される。この散乱は、細胞の核が光の波長の１０倍程度であるので、いわゆるミー（Ｍｉｅ）散乱と考えられる。ミー散乱は粒子径（ここでは核の大きさ）にその散乱強度が依存する。従って、異なる径を有する核が不均一に分布する系からの散乱強度を理論的に求めるのは容易ではない。
【００９７】
しかも、正常組織とみなせる細胞とがん組織とみなせる細胞とが不均一に分布するため、散乱も一様ではない。従って、ミュラー行列Ｍの１６要素の値のバラツキが大きくなり、１６要素の値そのままでは遷移部分であるとの判定は難しい。
【００９８】
そこで、この値のバラつきを残差という形で取り出し、バラつきが大きいことにより遷移部分であるか否かの判定を行えば、比較的簡易に遷移部分の判定が行えると考えられる。
【００９９】
次に、図１２のフローチャートを参照して、本実施の形態である観察装置１００を含む内視鏡システム１の動作について説明する。
【０１００】
まず、ステップＳ１では、光源ドライバ２５により光源２６を駆動し、この光源２６から光を発生させ、偏光生成部６０から所定の偏光状態の光を発生させて、対象物を照明する。ステップＳ２では、対象物からの戻り光を検出部５０のイメージセンサ５７により撮像する。
【０１０１】
ステップＳ３では、ステップＳ２でイメージセンサ５７が撮像した撮像信号に基づいて、画像処理部２３がイメージセンサ５７の画素単位でのミュラー行列Ｍの１６要素の２次元分布を算出する。
【０１０２】
ステップＳ４では、判定部２３ａが、ステップＳ３で得られたミュラー行列Ｍの１６要素の値を１６次元空間内にプロットする。ステップＳ５では、判定部２３ａが、ステップＳ４でプロットされたミュラー行列Ｍの１６要素の値を主成分分析により４次元の値に変換する。ステップＳ６では、判定部２３ａが、ステップＳ５で得られた４次元の値に基づいて残差を算出する。
【０１０３】
ステップＳ７では、判定部２３ａが、ステップＳ７における判定用の閾値を設定する。なお、この閾値は事前に設定されてもよい。ステップＳ８では、ステップＳ６で得られた残差と閾値とを比較し、閾値以下の情報を除去し、閾値以下を除去して残った値すなわち閾値以上の画素の情報を得る。この閾値以上の画素は遷移領域に対応する。ステップＳ９では、画像処理部２３が、ステップＳ８で得られた値を画像化する。ステップＳ９で得られた画像は、遷移部分の分布を示す画像である。遷移部分の分布を示す画像は、例えば、元の画像において、残差が閾値以上である画素を赤く表示することで遷移領域を赤で示すことにより形成される。この場合、遷移領域の境界を判断して輪郭を描き、赤い部分を半透明にしてもよいし、閾値以下の情報を除去して残った値の大きさに応じて色の濃さを変化させてもよい。
【０１０４】
以上のように構成された本実施の形態である観察装置１００は、対象物からの光の偏光状態を検出する偏光検出部４０と、偏光状態を示す値を主成分分析することにより、対象物の状態を判定する判定部２３ａとを有する。
【０１０５】
従って、判定部２３ａの判定結果を用いれば、遷移部分を正確に判定、検出することができる。
【０１０６】
また、本実施の形態である観察装置１００によれば、複雑な機構を用いることなく対象物の偏光状態を観察することができる。さらに、本実施の形態である観察装置１００では、装置全体の小型化を実現できる。
【０１０７】
なお、発明者らは、ミュラー行列Ｍの１６要素の値に基づいて光学的なエントロピーを算出し、このエントロピーの値により正常組織とがん組織とを判別することができるとの知見に至った。
【０１０８】
ミュラー行列をコヒーレント行列に変換した後、固有値分解して得られた固有値を使ってエントロピーを計算する手順は以下の通りである。
【０１０９】
まず、ミュラー行列Ｍ（４×４実数行列）をコヒーレント行列Ｃ（４×４複素数行列）に変換する。以下の説明では、ミュラー行列Ｍの１６要素は次式のように表現するものとする。
【数１５】
2020110291.tif 000017
【０１１０】
コヒーレント行列Ｃは次式で表現される。
【数１６】
2020110291.tif 000018
【０１１１】
コヒーレント行列Ｃの１６要素はミュラー行列Ｍの１６要素を用いて次式により与えられる。
【数１７】
2020110291.tif 000019
【０１１２】
次に、コヒーレント行列Ｃを対角化して固有値分解する。
【数１８】
2020110291.tif 000020
【０１１３】
そして、この正方行列Λの成分λ１〜λ４を用いて、コヒーレント行列ＣのエントロピーＨＴを次式により計算する。
【数１９】
2020110291.tif 000021
【０１１４】
なお、コヒーレント行列Ｃはエルミート行列なので、その固有値は必ず実数になる。また、λｉ＜０の場合はλｉ＝０とする。負の固有値を持つコヒーレント行列が得られた場合、元のミュラー行列Ｍには光学的な意味がない。この性質は、ミュラー行列Ｍが光学的な意味を持つものかどうかの判断に使われる
【０１１５】
上述の説明に基づいて、ミュラー行列ＭからエントロピーＨＴを算出したところ、がん組織から得られるミュラー行列Ｍの１６要素に基づくエントロピーＨＴは、それ以外の組織から得られるミュラー行列Ｍの１６要素に基づくエントロピーＨＴよりも低いという実験結果が得られた。
【０１１６】
また、がん組織のエントロピーの平均値と不偏標準偏差との和が、正常組織のエントロピーの平均値と不偏標準偏差との差より小さくなることが分かった。ここで、エントロピーの平均値は、各画素で計算されたエントロピーの画像毎の平均値であり、（[各画素のエントロピー]の総和）／（総画素数）で画像毎に計算される。また、不偏標準偏差は、データの広がり具合の一つの指標であり、｛（([各画素のエントロピー] - [エントロピーの平均値])2の総和）／（総画素数-1）｝の平方根として画像毎に計算される。
【０１１７】
従って、本実施の形態である観察装置１００において、判定部２３ａにより残差を算出して遷移部分を判定、検出し、エントロピーＨＴも算出してがん組織を判定、検出することができる。
【０１１８】
さらに、発明者らは、正常組織、がん組織及び遷移部分についてそれぞれ得られたミュラー行列Ｍの１６要素を主成分分析して得られた結果のうち、第２主成分、第３主成分及び第４主成分を３次元空間上にプロットしたところ、がん組織と遷移部分とを判別することができるとの知見に至った。
【０１１９】
まず、既に説明したように、正常組織、がん組織及び遷移部分についてそれぞれ得られたミュラー行列Ｍの１６要素を主成分分析して４次元に次元圧縮した。次いで、この４次元の値からなる４次元空間に対してクラスタリング等の手法を用いて探索をした。その結果、分散共分散行列の固有値である第１主成分が、ほぼ単位ジョーンズ（Jones）行列、つまり偏光の強度情報であることが判明した。
【０１２０】
そこで、発明者らは、４次元の値のうち第１主成分を除いた３次元の値、つまり、偏光変換作用成分だけが豊富な値についての探索を行った。
【０１２１】
図１３は、主成分分析の結果である第２主成分（主軸）、第３主成分（主軸）及び第４主成分（主軸）で表現される３次元空間をメルカトル図法により２次元に展開したものである。より詳細には、ある患者における遷移部分及びがん組織から得られた主成分分析の値を、同一の患者における正常組織の値により正規化し、正規化した値を３次元球面上にプロットしたものである。図において縦軸、つまり球面上の緯度は第２主成分の値であり、図中上半分がブロッホ（Bloch）球において｜０〉方向、下半分がブロッホ球において｜１〉方向である。また、図において横軸、つまり球面上の経度は第３主成分及び第４主成分の値である。
【０１２２】
図中、白丸はがん組織の値、黒丸は遷移部分の値である。図に示すように、がん組織の値はその８割近くが北半球の高緯度（４０°以上）に存在することがわかった。従って、図１３に示すような図を描き、緯度を閾値とすることで、がん組織と遷移部分とを区別して判定、検出することができる。
【０１２３】
以上、図面を参照して、実施の形態を詳述してきたが、具体的な構成は、この実施の形態及び実施例に限らず、その要旨を逸脱しない程度の設計的変更は、本開示に含まれる。
【０１２４】
一例として、上述の実施の形態である偏光検出部４０には６つの偏光生成部６０を設けたが、これに代えて、入射光のストークス成分のうちｓ１＝±１、ｓ２＝１、ｓ３＝１をそれぞれ生成する４つの偏光生成部６０を設けてもよい。そして、これら４つの偏光生成部６０のうち、１つの偏光生成部６０は１／４波長板６３を有する。
【０１２５】
さらに、複数の偏光生成部６０には光源２６から順に光が切り替えられて入射されていたが、各々の偏光生成部６０に周波数時分割した光を入射させてもよい。但し、あまり時分割の周期が長いと対象物の偏光度が変化してしまう可能性があるので、時分割の周期は人間が視認できない程度のものであることが望ましい。
【０１２６】
また、偏光生成部６０の配置位置は上述の各実施の形態に限定されず、検出部５０の周囲に配置されていればよい。但し、イメージセンサ５７の撮像面５７ａは通常矩形に形成されているので、観察装置１００全体の小型化を考慮すると、上述の実施の形態のように検出部５０の上下に、あるいは環状に偏光生成部６０を配置することが好ましい。
【０１２７】
さらに、上述の実施の形態における偏光生成部６０において、偏光板６１及び発散レンズ６２、さらには１／４波長板６３を互いに接着させて一体化してもよい。この場合、一部の偏光板６１（上述の実施の形態では６１ｃ〜６１ｆ）の対象物側には１／４波長板６３が設けられていない。そこで、１／４波長板６３の分だけ偏光板６１ｃ〜６１ｆの厚みを加えれば、それぞれの偏光生成部６０の光路長を等しくすることができる。
【０１２８】
同様に、光ファイバ１５の出射端１５ａと偏光板６１とを接着してもよい。
【０１２９】
また、偏光生成部６０の偏光変換素子は、観察装置１００の小型化という観点からすると、この偏光変換素子により得られる偏光状態が変化しない、つまりパッシブな偏光変換素子であることが好ましい。この観点から、上述の実施の形態では偏光板６１を用いている。しかし、対象物に光を照射する時点だけ偏光状態が変化しなければ足りるので、例えば液晶素子を電気駆動して偏光状態を所定の状態に維持したものも偏光変換素子として用いることができる。
【０１３０】
また、内視鏡装置１０が判定部２３ａを備えている例を示したが、これに代えて、内視鏡装置１０の一部ではなく外部装置として判定部２３ａを設けてもよい。この場合、内視鏡装置１０の画像処理部２３で処理された画像信号を判定部２３ａに送信する送信部を内視鏡装置１０に設けてもよい。
【０１３１】
また、観察の対象物の偏光状態を表す値としてミュラー行列を用いた例を示したが、これに代えて、ジョーンズ行列を用いてもよい。この場合、偏光生成部６０で生成する偏光の種類の数や、判定部２３ａでの主成分分析により圧縮される次元の値、残差の値、および、閾値を適宜設定することができる。
【０１３２】
また、判定部２３ａにおいて、高次のデータ群を用いて主成分分析することにより、複数の主成分ベクトルの線形結合で高次データをより低次元の空間で表す例を示したが、これに代えて、カーネル主成分分析を用いてもよい。この場合、例えば、判定部２３ａは、まず主成分分析により次元削減を行い、圧縮する前の次元における情報に対する圧縮された後の次元における情報の割合が所定の割合以下であったときに、カーネル主成分分析を行うと判断してもよい。所定の割合は、必要とする情報量に応じて設定され、例えば９０％である。カーネル主成分分析に用いられる正定値カーネルとして、例えば、線形カーネル、ガウスカーネル、多項式カーネル等が用いられる。
【０１３３】
複数の結果変数からなる多変量データを統計的に扱う手法であれば、主成分分析以外の多変量解析を用いてがん組織、非がん組織、および、遷移部分を区別して判定してもよい。主成分分析以外の多変量解析手法として、主成分分析およびカーネル主成分分析と数学的にほぼ等価な回帰分析（非線形を含む）および因子分析等が挙げられる。
【０１３４】
判定部２３ａは、以下のように、がん組織、非がん組織、および、遷移部分を判定してもよい。判定部２３ａは、例えば異なる照射方向ごとに、観察対象物の所定の領域からの光を撮像したイメージセンサ５７から出力された画像信号を用いて、対象物のミュラー行列を算出する。判定部２３ａは、算出した対象物のミュラー行列と、異常を有するサンプルを予め実測して得られたミュラー行列とを比較して、比較結果に基づいて、対象物の異常の有無を判定する。この場合、サンプルのミュラー行列（偏光特性）とサンプルの状態（異常の種類や異常の程度）とが予め関連付けされている。判定部２３ａは、比較の結果、対象物を実測して得られた実測結果のミュラー行列（偏光特性）との類似度が高いサンプルミュラー行列に関連付けされた状態を、実測した対象物の状態として判定する。
【０１３５】
対象物の状態の判定の基準となるミュラー行列を算出する際、判定部２３ａは、がん組織、非がん組織、または、遷移部分であると予め分かっている上記のサンプルの複数のミュラー行列を、それぞれニューラルネットワークに入力して学習させ、それぞれの状態の特徴量を抽出してもよい。この場合、がん組織の特徴を有するミュラー行列または行列要素、非がん組織の特徴を有するミュラー行列または行列要素、および、遷移部分の特徴を有するミュラー行列または行列要素をそれぞれ算出し、状態と対応付けて記憶してもよい。また、この場合、がん組織、非がん組織、および、遷移部分を示す複数の画像をニューラルネットワークに入力して学習させ、それぞれの状態に対応する画像の特徴量を抽出し、ミュラー行列または行列要素と画像とを対応付けて記憶してもよい。
【０１３６】
判定部２３ａは、状態判定の対象となる対象物のミュラー行列またはその行列要素と、ニューラルネットワークにより算出されたミュラー行列またはその行列要素とを比較することにより、対象物の状態を判定することができる。
【０１３７】
なお、上記のサンプルのミュラー行列は、異なる照射方位、異なる照射角度ごとに取得されたミュラー行列であり、予め記憶部に記憶されてもよい。この場合、判定部２３ａは、複数のサンプルミュラー行列の中から、算出したミュラー行列の照射方位と照射角度とが同一のサンプルミュラー行列を読み出し、比較を行う。制御部２１は、判定部２３ａによる判定の結果（すなわち対象物の異常の結果）をディスプレイ３２に表示させたり、記憶媒体に記憶させたりしてもよい。 サンプルのミュラー行列を一つの照射方位および照射角度で取得し、実測した複数のミュラー行列のそれぞれとサンプルミュラー行列とを比較することにより、対象物の異常の有無を判定してもよい。
【符号の説明】
【０１３８】
１    内視鏡システム
１０    内視鏡装置
１１    挿入部
１１ａ  管部材
１３    操作部
１５    光ファイバ
２０    本体部
２６    光源
４０    偏光検出部
５０    検出部
５７    イメージセンサ
６０    偏光生成部
１００  観察装置
【図１】
2020110291.tif 000022
【図２】
2020110291.tif 000023
【図３】
2020110291.tif 000024
【図４】
2020110291.tif 000025
【図５】
2020110291.tif 000026
【図６】
2020110291.tif 000027
【図７】
2020110291.tif 000028
【図８】
2020110291.tif 000029
【図９】
2020110291.tif 000030
【図１０】
2020110291.tif 000031
【図１１】
2020110291.tif 000032
【図１２】
2020110291.tif 000033
【図１３】
2020110291.tif 000034