(57)【要約】      （修正有）【課題】眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得る眼科解析装置において、信頼性のある解析結果を取得する。【解決手段】本発明の眼科解析装置１０は、前記画像データを解析するための解析処理手段を備え、前記解析処理手段は、眼科用ＯＣＴ装置１によって取得され被検眼の偏光特性を示す画像データである第１のＯＣＴ画像データを処理して第１の解析結果を取得し、眼科用ＯＣＴ装置によって取得され前記第１のＯＣＴ画像データとは異なる画像データである第２のＯＣＴデータを処理して第２の解析結果を取得し、前記第１の解析結果と前記第２の解析結果を統合処理することによって統合解析結果を取得する。【選択図】図１
2018068778.tif 000002
【特許請求の範囲】
【請求項１】
眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得る眼科解析装置であって、
前記画像データを解析するための解析処理手段を備え、
前記解析処理手段は、
眼科用ＯＣＴ装置によって取得され被検眼の偏光特性を示す画像データである第１のＯＣＴ画像データを処理して第１の解析結果を取得し、
眼科用ＯＣＴ装置によって取得され前記第１のＯＣＴ画像データとは異なる画像データである第２のＯＣＴデータを処理して第２の解析結果を取得し、
前記第１の解析結果と前記第２の解析結果を統合処理することによって統合解析結果を取得することを特徴とする眼科解析装置。
【請求項２】
前記解析処理手段は、
前記第１のＯＣＴ画像データとして被検眼眼底の篩状板領域の偏光特性を示す画像データを処理して篩状板領域に関する第１の解析結果を取得することを特徴とする請求項１に記載の眼科解析装置。
【請求項３】
前記解析手段は、
前記第１の解析結果として前記篩状板領域の動径方向と円周方向の少なくともいずれかに関する解析パラメータを得ることを特徴とする請求項１〜２のいずれかの眼科解析装置。
【請求項４】
前記解析手段は、
篩状板領域に関する第１の解析結果として、前記篩状板領域の等方性に関する解析結果を得ることを特徴とする請求項１〜３のいずれかの眼科解析装置。
【請求項５】
前記解析手段は、
前記第１の解析結果として被検眼眼底の強膜の偏光解消度を得ることを特徴とする請求項１〜４のいずれかの眼科解析装置。
【請求項６】
前記解析処理手段は、
前記第２のＯＣＴ画像データとして被検眼のＯＣＴアンジオ画像データを処理して第２の解析結果を取得することを特徴とする請求項１〜５のいずれかの眼科解析装置。
【請求項７】
前記解析処理手段は、
前記第２のＯＣＴ画像データとして被検眼の反射強度特性を示すＯＣＴ画像データを処理して第２の解析結果を取得することを特徴とする請求項１〜６のいずれかの眼科ＯＣＴ解析装置。
【請求項８】
前記解析処理手段は、
前記第１のＯＣＴ画像データを処理して前記第１の解析結果として篩状板領域の位置情報を取得し、
前記第２のＯＣＴ画像データを処理して前記第２の解析結果として篩状板領域の位置情報を取得し、
前記第１のＯＣＴ画像データに基づく位置情報と、前記第２のＯＣＴ画像データに基づく位置情報とを統合処理することによって統合解析結果を得ることを特徴とする請求項１〜７のいずれかの眼科ＯＣＴ解析装置。
【請求項９】
前記解析手段は、
さらに、前記第１の解析結果として被検眼眼底の篩状板の複屈折性及び強膜の偏光解消度を取得し、
篩状板の複屈折性に関する第１の解析結果と強膜の偏光解消度に関する第１の解析結果を統合処理することによって第２の統合解析結果を得ることを特徴とする請求項１〜８のいずれかの眼科ＯＣＴ解析装置。
【請求項１０】
前記解析処理手段は、前記第１の解析結果と前記第２の解析結果として互いに共通する解析パラメータを取得する請求項１〜９のいずれかの眼科ＯＣＴ解析装置。
【請求項１１】
前記解析処理手段は、前記第１の解析結果と前記第２の解析結果の一致度に基づいて前記統合解析結果を得ることを特徴とする請求項１〜１０のいずれかの眼科ＯＣＴ解析装置。
【請求項１２】
前記解析処理手段は、前記第１の解析結果及び前記第２の解析結果に関してそれぞれ篩状板領域上に設定された基準位置を中心とする少なくとも４つの分割領域に関する解析パラメータの代表値を取得し、前記第１の解析結果及び前記第２の解析結果に関する前記代表値に基づいて統合解析結果を得ることを特徴とする請求項１〜１１のいずれかの眼科ＯＣＴ解析装置。
【請求項１３】
前記解析処理手段は、前記第１の解析結果と前記第２の解析結果とを表示部に同時に表示することを特徴とする請求項１〜１２のいずれかの眼科ＯＣＴ解析装置。
【請求項１４】
前記解析処理手段は、第１のＯＣＴ画像データとしてRahe gapとRahe gapの周辺の偏光特性を示す画像データを処理してRahe gap又はRahe gap周辺のいずれかに関する第１の解析結果を取得することを特徴とする請求項１〜１３のいずれかの眼科ＯＣＴ解析装置。
【請求項１５】
眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得る眼科解析装置であって、
前記画像データを解析するための解析処理手段を備え、
前記解析処理手段は、
眼科用ＯＣＴ装置によって取得され被検眼の偏光特性を示す画像データである第１のＯＣＴ画像データを処理して、Rahe gap又はRahe gap周辺のいずれかに関する第１の解析結果を取得することを特徴とする眼科ＯＣＴ解析装置。
【請求項１６】
請求項１〜１５のいずれかに記載の眼科ＯＣＴ解析装置の解析処理手段としてコンピュータを機能させるための眼科解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、被検眼を解析するための眼科ＯＣＴ解析装置、眼科解析プログラムに関する。
【背景技術】
【０００２】
ＯＣＴ画像を解析して解析結果を得る装置が知られている。例えば、眼科用ＯＣＴ装置によって得られた被検眼の反射特性を示すＯＣＴ画像から特徴部位（例えば、篩状板）を検出し、解析結果を得る方法が知られている。一方、ＰＳ−ＯＣＴ装置によって得られた被検眼の偏光特性を示すＯＣＴ画像に基づいて解析結果を得る情報が知られている（例えば、非特許文献１）。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Quantitative evaluation of optical properties of optic nerve head tissues by multi-contrast optical coherence tomographyInvest. Ophthalmol. Vis. Sci.. 2016;57(12):828
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、被検眼の反射特性を示すＯＣＴ画像に基づいて解析結果を得る場合、機能的な信号を検出することが難しく、ＯＣＴ画像のみで被検眼の微小な変化を検出することは困難な場合があった。また、ＰＳ−ＯＣＴによって得られた偏光ＯＣＴ画像に基づいて解析を行う場合、被検眼の偏光特性を含む信号を検出できるので、ＯＣＴ画像よりも被検眼の微小な変化を検出できる場合があるが、偏光ＯＣＴ画像は、ＯＣＴ画像よりもノイズによるアーティファクトが生じる可能性が高く、偏光ＯＣＴ画像のみでは、解析結果としての信頼性が必ずしも十分とはいえなかった。
【０００５】
また、ＯＣＴ画像又は偏光ＯＣＴ画像を用いた各解析の手法においても、改善の余地があると考えられる。
【０００６】
本開示は、上記問題点を鑑み、信頼性のある解析結果を取得できる眼科ＯＣＴ解析装置及び眼科解析プログラムを提供することを技術課題とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明は以下のような構成を備えることを特徴とする。
【０００８】
（１）  眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得る眼科解析装置であって、
前記画像データを解析するための解析処理手段を備え、
前記解析処理手段は、
眼科用ＯＣＴ装置によって取得され被検眼の偏光特性を示す画像データである第１のＯＣＴ画像データを処理して第１の解析結果を取得し、
眼科用ＯＣＴ装置によって取得され前記第１のＯＣＴ画像データとは異なる画像データである第２のＯＣＴデータを処理して第２の解析結果を取得し、
前記第１の解析結果と前記第２の解析結果を統合処理することによって統合解析結果を取得することを特徴とする。
（２）  眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得る眼科解析装置であって、
前記画像データを解析するための解析処理手段を備え、
前記解析処理手段は、
眼科用ＯＣＴ装置によって取得され被検眼の偏光特性を示す画像データである第１のＯＣＴ画像データを処理して、Rahe gap又はRahe gap周辺のいずれかに関する第１の解析結果を取得することを特徴とする。
【発明の効果】
【０００９】
本開示によれば、信頼性のある解析結果を取得できる。
【発明を実施するための形態】
【００１０】
以下、本開示における典型的な実施形態について説明する。
【００１１】
＜概要＞
眼科ＯＣＴ解析装置は、眼科用ＯＣＴ装置によって取得された被検眼の画像データを処理して被検眼の解析結果を得るために用いられてもよい（例えば、図１参照）。眼科ＯＣＴ解析装置は、画像データを解析処理するための解析処理部を備えてもよい。さらに、眼科ＯＣＴ解析装置は、眼科用ＯＣＴ装置の装置本体と一体化された装置であってもよいし、眼科用ＯＣＴ装置の装置本体とは別に設けられた装置であってもよい。また、眼科ＯＣＴ解析装置として、汎用のパーソナル・コンピュータが用いられてもよい。また、解析処理部としては、例えば、ＣＰＵ等のプロセッサが用いられてもよい。被検眼の画像データは、例えば、被検眼に照射された測定光と参照光によるＯＣＴ信号に基づいて取得されてもよい。なお、眼科ＯＣＴ解析装置には、眼科用ＯＣＴ装置によって取得された画像データを処理して解析結果を得るための眼科解析プログラムが用いられてもよい。
【００１２】
解析結果としては、被検眼に関する解析パラメータが取得されてもよく、解析パラメータとしては、例えば、画像データを解析して得られる物理量であってもよいし、画像データを解析して得られる物理量を符号化したパラメータであってもよい。
【００１３】
被検眼の画像データとしては、例えば、眼底の画像データであってもよいし、前眼部の画像データであってもよい。なお、眼底の画像データとしては、例えば、眼底の篩状板領域を少なくとも含む乳頭部の画像データであってもよいし、眼底の黄斑部を少なくとも含む画像データ、眼底のRahe gapとRahe gapの周辺を少なくとも含む画像データであってもよい。もちろん、眼底上の部位は、これらに限定されない。
【００１４】
＜第１の解析結果の取得＞
解析処理部は、例えば、眼科用ＯＣＴ装置によって取得された第１のＯＣＴ画像データを処理して第１の解析結果を取得してもよい。第１のＯＣＴ画像データは、被検眼の偏光特性を示す画像データであってもよい。また、第１のＯＣＴ画像データは、眼科用ＯＣＴ装置から出力されるＯＣＴ信号に基づいて取得された画像データであってもよい。
【００１５】
第１のＯＣＴ画像データを取得するための眼科用ＯＣＴ装置としては、例えば、ＰＳ−ＯＣＴ装置であってもよい。なお、ＰＳ−ＯＣＴ装置としては、被検眼に照射された測定光と参照光によるＯＣＴ信号であり互いに直交する偏光成分を持つ複数のＯＣＴ信号を取得するためのＰＳ−ＯＣＴ光学系を備える構成であってもよい。
【００１６】
第１のＯＣＴ画像データとしては、偏光特性を示すＯＣＴデータ（以下、ＰＳＯＣＴデータ）が異なる位置に関して並べられたＢスキャンＰＳＯＣＴ画像データであってもよい。また、第１のＯＣＴ画像データとしては、深さ方向に直交する方向に関する二次元範囲においてＰＳＯＣＴデータが並べられた３次元ＰＳ画像データであってもよい。
【００１７】
第１のＯＣＴ画像データとしては、３次元ＰＳ画像データに基づく正面ＰＳ画像データ（エンフェイスＰＳ画像データ）であってもよく、例えば、３次元ＰＳ画像データの深さ方向全体に基づく正面ＰＳ画像データであってもよいし、３次元ＰＳ画像データの深さ方向の一部に基づく正面ＰＳ画像データであってもよい。演算手法としては、深さ方向における輝度値の積算処理であってもよいし、他の手法（例えば、ヒストグラム演算、最大値抽出等）であってもよい。
【００１８】
解析処理部は、例えば、第１のＯＣＴ画像データとして眼底領域の偏光特性を示す画像データを処理して眼底領域に関する第１の解析結果を取得してもよい。眼底領域としては、例えば、篩状板領域であってもよいし、篩状板を含む乳頭部領域であってもよいし、黄斑領域、Rahe gapとRahe gap周辺の少なくともいずれかの領域であってもよい。もちろん眼底領域としては、これらの領域に限定されてない。
【００１９】
この場合、第１の解析結果としては、偏光特性に関する解析結果であってもよく、例えば、眼底領域の複屈折性に関する情報（複屈折情報）、眼底領域の偏光非一様性（DOPU）情報、眼底領域の偏光軸回転（Axis-Orientation）情報、の少なくともいずれかであってもよい。なお、これらの情報の算出方法については、当業者に周知の技術を採用できるので説明を省略する。また、これに限定されず、第１の解析結果は、眼底領域の偏光特性を示す画像データに基づいて得られる眼底領域の形態に関する解析結果であってもよく、眼底領域の位置情報、眼底領域の少なくとも一つの層の厚み情報の少なくともいずれかであってもよい。また、第１の解析結果は、眼底領域の密度情報であってもよく、例えば、神経線維層の密度情報であってもよい。
【００２０】
解析処理部は、例えば、第１の解析結果として眼底領域の動径方向と円周方向の少なくともいずれかに関する解析パラメータを得てもよい。なお、この場合、眼底領域上に設定された基準位置を中心とする動径方向と円周方向の少なくともいずれかに関する解析パラメータが得られてもよい。なお、異なる動径方向に関して解析パラメータが取得されてもよく、少なくとも上下左右の４方向に関して解析パラメータが取得されてもよい。この場合、異なる円周方向に関して解析パラメータが取得されてもよい。
【００２１】
なお、眼底領域の複屈折情報としては、例えば、動径方向における複屈折性の変化であってもよく、変化量、変化率、又は局所的な変化における曲率であってもよい。また、眼底領域の位置情報としては、眼底領域の形状情報であってもよく、例えば、動径方向における眼底領域の形状変化であってもよく、眼底領域の傾き、形状変化量、形状変化率、又は局所的な変化における曲率であってもよい。形状変化は、例えば、眼底領域の深さ方向の位置変化であってもよい。この場合、位置情報は、絶対な位置情報として算出されてもよいし、相対的な位置情報として算出されてもよい。
【００２２】
解析処理部は、第１の解析結果として眼底領域の等方性に関する解析結果を取得してもよい。この場合、解析処理部は、例えば、眼底領域に関する解析パラメータを取得し、取得された解析パラメータに基づいて眼底領域の等方性（不等方性）を求めてもよい。等方性を求めるベースとなる解析パラメータとしては、眼底領域に関する解析パラメータであってもよい。この場合、解析パラメータは、動径方向と円周方向の少なくともいずれかに関する解析パラメータの変化を示す変化情報として得られてもよい。
【００２３】
なお、眼底領域上の複数の領域に関する解析パラメータが動径方向と円周方向の少なくともいずれかに関して比較されることで、等方性が求められてもよい。なお、複数の領域としては、眼底領域上に設定された基準位置を中心として設定されてもよい。
【００２４】
なお、上記解析は、眼底領域として篩状板領域が設定され、篩状板を解析する場合に特に有利である。篩状板の解析は、被検眼の緑内障診断において重要な部位であり、篩状板領域に対して上記解析が行われることで、緑内障診断のための信頼性の高い結果が得られる。
【００２５】
なお、篩状板領域に関して相対的な位置情報を得る場合、開口部と篩状板との間の動径方向における距離変化であってもよいし、開口部と篩状板との間の距離に関する動径方向での変化率であってもよいし、当該変化率の局所的な曲率であってもよい。
【００２６】
解析処理部は、例えば、第１の解析結果として眼底領域の偏光解消度を取得してもよい。眼底領域としては、例えば、強膜の偏光解消度を得てもよく、乳頭部周辺における強膜の偏光解消度が得られてもよい。また、眼底領域としては、眼底のRahe gap周辺における偏光解消度が得られてもよい。
【００２７】
＜第２の解析結果の取得＞
解析処理部は、例えば、眼科用ＯＣＴ装置によって取得され第２のＯＣＴデータを処理して第２の解析結果を取得してもよい。第２のＯＣＴデータは、第１のＯＣＴ画像データとは異なる画像データであってもよい。第２のＯＣＴ画像データは、眼科用ＯＣＴ装置から出力されるＯＣＴ信号に基づいて取得された画像データであってもよい。
【００２８】
第２のＯＣＴ画像データを取得するための眼科用ＯＣＴ装置としては、例えば、第１のＯＣＴ画像データを取得するための眼科用ＯＣＴ装置と同一の眼科用ＯＣＴ装置（例えば、ＰＳＯＣＴ装置）であってもよい。この場合、第２のＯＣＴデータは、第１のＯＣＴ画像データを得る際に取得された第２のＯＣＴデータ（例えば、第１のＯＣＴ画像データの基礎となるＯＣＴ画像データ）であってもよいし、別タイミングで取得された第２のＯＣＴデータであってもよい。
【００２９】
また、第２のＯＣＴ画像データを取得する眼科用ＯＣＴ装置としては、第１のＯＣＴ画像データを取得するための眼科用ＯＣＴ装置とは別の眼科用ＯＣＴ装置であってもよい。この場合、眼科用ＯＣＴ装置は、被検眼に照射された測定光と参照光によるＯＣＴ信号を取得するためのＯＣＴ光学系を備える眼科用ＯＣＴ装置であってもよい。
【００３０】
第２のＯＣＴ画像データとしては、ＰＳＯＣＴデータとは異なる第２のＯＣＴデータが異なる位置に関して並べられたＢスキャン第２ＯＣＴ画像データであってもよい。また、第２のＯＣＴ画像データとしては、深さ方向に直交する方向に関する二次元範囲において第２のＯＣＴデータが並べられた３次元第２ＯＣＴ画像データであってもよい。
【００３１】
第２のＯＣＴ画像データとしては、３次元ＯＣＴ画像データに基づく正面第２ＯＣＴ画像データ（エンフェイス第２ＯＣＴ画像データ）であってもよく、例えば、３次元第２ＯＣＴ画像データの深さ方向全体に基づく正面第２ＯＣＴ画像データであってもよいし、３次元第２ＯＣＴ画像データの深さ方向の一部に基づく正面第２ＯＣＴ画像データであってもよい。演算手法としては、深さ方向における輝度値の積算処理であってもよいし、他の手法（例えば、ヒストグラム演算、最大値抽出等）であってもよい。
【００３２】
第２のＯＣＴ画像データとしては、例えば、被検眼の反射強度特性を示すＯＣＴ画像データ、被検眼のＯＣＴアンジオ画像データ（例えば、ＯＣＴモーションコントラスト画像データ）、被検眼のドップラー特性を示すドップラーＯＣＴ画像データであってもよい。
【００３３】
ＯＣＴアンジオ画像データは、例えば、同一位置に関して時間的に異なる少なくとも２つのＯＣＴ信号が処理されることによって、取得されてもよい。ＯＣＴアンジオ画像データは、被検眼の血管領域を画像化したＯＣＴ画像データであってもよい。この場合、血管領域が輝度値として表現された画像データであってもよい。また、ドップラーＯＣＴ画像データは、ドップラ―ＯＣＴによって取得されてもよい。
【００３４】
解析処理部は、被検眼のＯＣＴアンジオ画像データを処理して第２の解析結果を取得してもよい。この場合、第２の解析結果として被検眼の血管領域に関する解析パラメータが取得されてもよい。解析パラメータとしては、血管像の強度、血管密度（例えば、３次元密度、特定の層の血管密度）であってもよい。その他、血管面積、血管体積、血管径、血管蛇行度等が、解析パラメータとして取得されてもよい。
【００３５】
なお、ＯＣＴアンジオ画像データに基づく第２の解析結果としては、眼底の強膜に関する解析結果であってもよく、例えば、乳頭部周辺における強膜に関する解析結果が得られてもよい。
【００３６】
＜統合解析結果の取得＞
解析処理部は、例えば、第１の解析結果と第２の解析結果を統合処理することによって統合解析結果を取得してもよい。統合解析結果として、例えば、解析パラメータ、各領域（セクション）での解析パラメータの比較結果、解析マップ、解析チャートが取得されてもよい。なお、統合解析結果は、表示部、プリンタ、外部サーバーの少なくともいずれかに出力されてもよい。
【００３７】
この場合、統合解析結果としては、例えば、第１のＯＣＴ画像データに基づく第１の解析パラメータと第２のＯＣＴ画像データに基づく第２の解析パラメータとを統合した統合解析パラメータ（例えば、第１の解析パラメータと第２の解析パラメータの代表値）であってもよい。
【００３８】
また、統合解析結果としては、例えば、第１のＯＣＴ画像データに基づく第１の解析マップと第２のＯＣＴ画像データに基づく第２の解析マップとを統合した統合解析マップ（例えば、各マップの各位置における解析パラメータの代表値）であってもよい。
【００３９】
また、統合解析結果としては、例えば、第１のＯＣＴ画像データに基づく第１の解析チャートと第２のＯＣＴ画像データに基づく第２の解析チャートとを統合した統合解析チャート（例えば、各解析チャートのセクションでの解析パラメータの代表値）であってもよい。この場合、統合解析チャートにおける各セクションでの解析パラメータを、動径方向と円周方向の少なくともいずれかに関して比較し、大きな変化が検出されたセクション間の境界を強調表示するようにしてもよい。
【００４０】
なお、統合解析結果の取得において、代表値を求める場合、第１の解析結果と第２の解析結果との重み付けの係数が異なってもいてもよいし、同じ係数であってもよい。
【００４１】
上記のようにして、第１のＯＣＴ画像データに基づく第１の解析結果と、第２のＯＣＴ画像データに基づく第２の解析結果を統合した統合解析結果を得ることで、信頼性の高い解析結果を得ることができる。結果として、検者は、統合解析結果を用いて良好な診断を行うことができる。
【００４２】
なお、統合解析は、検出が比較的に難しい篩状板領域の解析において特に有利である。篩状板の解析は、被検眼の緑内障診断において重要な部位であり、篩状板領域に対し統合解析が行われることで、緑内障診断のための信頼性の高い結果が得られる。
【００４３】
なお、上記説明において、解析マップは、例えば、各位置での解析パラメータの二次元的な分布を示す解析マップであってもよい。解析チャートとしては、例えば、複数の領域（セクション）に応じて解析パラメータの代表値を示す解析チャートとして表現されてもよい。
【００４４】
解析処理部は、例えば、第１の解析結果と第２の解析結果として互いに共通する解析パラメータを取得してもよい。共通する解析パラメータとしては、例えば、眼底領域の位置情報であってもよい。
【００４５】
なお、統合解析に際し、解析処理部は、例えば、第１のＯＣＴ画像データを処理して第１の解析結果として眼底領域の位置情報を取得し、第２のＯＣＴ画像データを処理して第２の解析結果として、第１の解析結果を得た眼底領域と同一の眼底領域の位置情報を取得てもよい。この場合、解析処理部は、例えば、第１のＯＣＴ画像データに基づく位置情報と、第２のＯＣＴ画像データに基づく位置情報とを統合処理することによって統合解析結果を得てもよい。同一の眼底領域としては、例えば、篩状板領域であってもよいし、篩状板を含む乳頭部領域であってもよいし、黄斑領域であってもよい。もちろん眼底領域としては、これらの領域に限定されてない。
【００４６】
解析処理部は、例えば、第１の解析結果と第２の解析結果の一致度に基づいて統合解析結果を得てもよい。この場合、例えば、第１の解析結果と第２の解析結果の一致度としては、第１の解析結果を複数の領域に関して取得したときの異なる領域間での相対的な比較結果（例えば、差分、比率）と、第２の解析結果を複数の領域に関して取得したときの異なる領域間での相対的な比較結果（例えば、差分、比率）と、の間の一致度が求められてもよい。この場合、一致度が高いほど、信頼性の高い結果として扱われてもよい。
【００４７】
例えば、動径方向に関して相対的な比較結果が取得されてもよいし、円周方向に関して相対的な比較結果が取得されてもよい。また、相対的な比較結果に限定されず、各解析結果の絶対値の大小に関する一致度が求められてもよい。
【００４８】
なお、統合解析は、第１の解析結果に関して行われてもよい。解析処理部は、例えば、第１の解析結果として眼底の篩状板の複屈折性及び強膜の偏光解消度を取得し、篩状板の複屈折性に関する第１の解析結果と強膜の偏光解消度に関する第１の解析結果を統合処理することによって第２の統合解析結果を得てもよい。これによって、篩状板と強膜を統合した解析が可能となり、信頼性の高い結果が得られる。
【００４９】
なお、解析処理部は、例えば、第１の解析結果と第２の解析結果を同時に表示部に表示してもよい。この場合、第１の解析結果と第２の解析結果の表示形態は、同一の表示形態であってもよいし、異なる表示形態であってもよい。各解析結果の表示形態としては、例えば、解析パラメータ、各領域（セクション）での解析パラメータの比較結果、解析マップ、解析チャートであってもよい。
【００５０】
＜篩状板解析における基準位置の設定＞
解析処理部は、例えば、第１のＯＣＴ画像データと第２のＯＣＴ画像データの少なくともいずれかを処理して眼底の乳頭部の開口部の重心を求めてもよい。解析処理部は、例えば、篩状板領域における解析パラメータの変化を、重心を中心として求めてもよい。これによって、篩状板領域の解析における基準位置を規定できる。
【００５１】
篩状板領域上に設定された基準位置としては、開口部の重心であってもよい。解析パラメータの変化は、例えば、重心を中心とする動径方向と周方向の少なくともいずれかに関して求められてもよい。もちろん、これに限定されず、他の位置が基準位置として設定されてもよい。
【００５２】
＜分割領域に関する代表値の取得＞
解析処理部は、例えば、第１の解析結果及び第２の解析結果に関してそれぞれ眼底領域上に設定された基準位置を中心とする少なくとも４つの分割領域に関する解析パラメータの代表値を取得してもよい。解析処理部は、例えば、分割領域において取得された第１の解析結果及び第２の解析結果に関する代表値に基づいて統合解析結果を得てもよい。この場合、統合解析結果としては、眼底領域上に設定された基準位置を中心とする少なくとも４つの分割領域に関してそれぞれ第１の解析結果（例えば、第１の解析パラメータ）と第２の解析結果（例えば、第２の解析パラメータ）とを統合した統合解析結果（例えば、統合解析パラメータ）が取得されてもよい。
【００５３】
この場合、代表値とは、例えば、基準位置から等距離の円内ないし二つの円に挟まれた領域における平均値、メディアン、分散値、最大値、最小値、PV値のいずれかであってもよい。
【００５４】
この場合、解析処理部は、例えば、少なくとも４つの領域における解析パラメータの代表値を複数算出し、複数の代表値同士で周方向に関する変化パターンが一致しているかに応じて信頼性に関する解析パラメータ（例えば、解析値）を増減してもよい。
【００５５】
この場合、解析処理部は、例えば、得られた篩状板に関する解析パラメータが、多項式にてフィットして解析してもよい。
【００５６】
上記の少なくとも４つの分割領域としては、少なくともNasal、Temporal、Superior、Inferiorの４つの領域であってもよい。この場合、少なくともNasal、Temporalの領域とSuperior、Inferiorの領域では、フィッティングの関数系を異ならせてもよい。また、フィッティングの関数系は、Nasal、Temporalの領域では動径方向に極値を持たず、Superior、Inferiorの領域では動径方向に極値を持っていてもよい。
【００５７】
また、分割領域に関する解析パラメータの代表値は、動径方向に関して比較されてもよいし、円周方向に関して比較されてもよい。比較手法としては、例えば、代表値間の差分、比率であってもよい。
【００５８】
なお、眼底領域としては、例えば、篩状板領域であってもよいし、篩状板を含む乳頭部領域であってもよいし、黄斑領域であってもよい。もちろん眼底領域としては、これらの領域に限定されてない。
【００５９】
＜PS-OCTを用いたRahe gapの解析＞
解析処理部は、例えば、例えば、第１のＯＣＴ画像データとしてRahe gap周辺の偏光特性を示す画像データを処理してRahe gap周辺に関する第１の解析結果を取得してもよい。なお、Rahe gapとは、視神経乳頭から黄斑部に向けて結ばれる直線の延長の近傍に形成される部分（黄斑部よりも耳側）であり、神経線維層の上下方向の境界ともいえる。第１の解析結果は、解析パラメータ、各領域（セクション）での解析パラメータの比較結果、解析マップ、解析チャートが取得されてもよく、これらが出力されてもよい。この場合、Rahe gapを境界とする解釈結果が出力されてもよい。
【００６０】
この場合、例えば、解析処理部は、第１の解析結果として、Rahe gap周辺の偏光特性に関する解析結果を取得してもよい。偏光特性としては、前述したように、複屈折情報、偏光非一様性（DOPU）情報、偏光軸回転情報の少なくともいずれかであってもよい。解析処理部は、第１の解析結果として、Rahe gap周辺での視神経線維層の密度情報であってもよい。
【００６１】
例えば、解析処理部は、Rahe gapを境界として、Rahe gapよりも上側の解析結果と、Rahe gapよりも下側の解析結果とを取得してもよい。さらに、解析処理部は、Rahe gapを境界として、各解析結果を比較処理してもよい。比較処理としては、差分、比率を求めてもよいし、対称性、一致度を求めるようにしてもよい。
【００６２】
例えば、Rahe gapよりも上側の偏光軸回転情報と、Rahe gapよりも下側の偏光軸回転情報の一致度を求めてもよい。この場合、緑内障の早期においては、被検眼眼底の変化が非対称に進行すると言われている。そこで、偏光軸回転情報の一致度を求めることで、眼底の微小な変化を検出可能であり、緑内障の進行を早期に検出できる。例えば、偏光軸の向きが大きく異なった場合、緑内障が早期進行していると判断可能である。
【００６３】
なお、解析処理部は、例えば、例えば、第２のＯＣＴ画像データとしてRahe gap周辺の画像データを処理してRahe gap周辺に関する第２の解析結果を取得してもよい。さらに、解析処理部は、例えば、Rahe gap周辺に関する第１の解析結果と第２の解析結果を統合処理することによって統合解析結果を取得してもよい。
【００６４】
なお、解析処理部は、Rahe gapでの解析結果と、乳頭周辺の形状情報（構造情報）と、隅角情報とを統合して、緑内障の程度を判定してもよい。
【００６５】
なお、解析処理部は、例えば、例えば、第１のＯＣＴ画像データを解析処理してRahe gapを検出し、Rahe gapに関する第１の解析結果を取得してもよい。第１の解析結果としては、例えば、Rahe gapの偏光特性に関する解析結果、Rahe gapの形状に関する解析結果であってもよい。
【００６６】
なお、上記実施形態に示した各手法につき、独立して実施されてもよいし、連動して実施されてもよい。例えば、＜第１の解析結果の取得＞、＜PS-OCTを用いたRahe gapの解析＞に示した少なくとも一つの手法につき、＜第２の解析結果の取得＞、＜統合解析結果の取得＞とは独立して実施可能であり、＜第１の解析結果の取得＞、＜PS-OCTを用いたRahe gapの解析＞に示した他の手法とも独立して実施可能である。
【００６７】
＜実施例＞
以下、本実施例の眼科ＯＣＴ解析装置について図面を用いて説明する。図１に示す眼科ＯＣＴ解析装置（以下、ＯＣＴ解析装置）１０は、ＯＣＴ装置１によって取得された画像データ（例えば、偏光特性を示す第１のＯＣＴ画像データ、偏光特性とは異なる特性を示す第２のＯＣＴ画像データ）を解析処理して被検眼の解析結果を得る。
【００６８】
ＯＣＴ解析装置１は、例えば、制御部７０を備える。制御部７０は、例えば、一般的なＣＰＵ（Central Processing Unit）７１、ＲＯＭ７２、ＲＡＭ７３、等で実現される。ＲＯＭ７２には、例えば、ＯＣＴ装置１によって取得された画像データを処理するための解析処理プログラム、ＯＣＴ装置１の動作を制御して画像データ（例えば、偏光特性を示す第１のＯＣＴ画像データ、偏光特性とは異なる特性を示す第２のＯＣＴ画像データ）を得るためのプログラム、初期値等が記憶される。ＲＡＭ７３は、例えば、各種情報を一時的に記憶する。
【００６９】
制御部７０には、図１に示すように、例えば、記憶部（例えば、不揮発性メモリ）７４、操作部７６、および表示部７５等が電気的に接続されている。記憶部７４は、例えば、電源の供給が遮断されても記憶内容を保持できる非一過性の記憶媒体である。例えば、ハードディスクドライブ、フラッシュＲＯＭ、着脱可能なＵＳＢメモリ等を記憶部７４として使用することができる。
【００７０】
操作部７６には、検者による各種操作指示が入力される。操作部７６は、入力された操作指示に応じた信号をＣＰＵ７１に出力する。操作部７６には、例えば、マウス、ジョイスティック、キーボード、タッチパネル等の少なくともいずれかのユーザーインターフェイスを用いればよい。
【００７１】
表示部７５は、装置１の本体に搭載されたディスプレイであってもよいし、本体に接続されたディスプレイであってもよい。例えば、パーソナルコンピュータ（以下、「ＰＣ」という）のディスプレイを用いてもよい。表示部７５は、例えば、ＯＣＴ装置１によって取得された画像データ、又は解析結果等を表示する。
【００７２】
なお、本実施例のＯＣＴ解析装置１０には、例えば、ＯＣＴ装置１が接続されている。なお、ＯＣＴ解析装置１０は、例えば、ＯＣＴ装置１と同一の筐体に収納された一体的な構成であってもよいし、別々の構成であってもよい。制御部７０は、接続されたＯＣＴ装置１から画像データを取得してもよい。制御部７０は、ＯＣＴ装置１によって取得された画像データを記憶媒体を介して取得してもよい。
【００７３】
図２は、ＯＣＴ装置１の一例を示す図である。被検物は、眼の眼底である。実施例の装置は、干渉信号における互いに直交する偏光成分を異なる検出器にて検出可能な構成を有する。
【００７４】
光コヒーレンストモグラフィー（ＯＣＴ）装置１は、波長掃引式ＯＣＴ（ＳＳ−ＯＣＴＨwepＴＨource-OCT）を基本的構成とし、波長可変光源１０２、干渉光学系（ＯＣＴ光学系）１００、演算制御器（演算制御部）７０と、を含む。その他、ＯＣＴ装置１には、メモリ７２、モニタ７５、図示無き正面像観察系及び固視標投影系が設けられる。演算制御器（以下、制御器（制御部））７０は、波長可変光源１０２、干渉光学系１００、メモリ７２、モニタ７５に接続されている。
【００７５】
ＯＣＴ光学系１００には、ＳＳ−ＯＣＴ方式が用いられ、光源１０２として出射波長を時間的に高速で変化させる波長可変光源（波長走査型光源）が用いられる。光源１０２は、例えば、レーザ媒体、共振器、及び波長選択フィルタによって構成される。そして、波長選択フィルタとして、例えば、回折格子とポリゴンミラーの組み合わせ、ファブリー・ペローエタロンを用いたフィルタが挙げられる。
【００７６】
本実施例では、瞬間輝線幅が短く、共振器長が短い光源としてＡＸＳＵＮ社のTUNABLE LASER が用いられる（例えば、λc=1060nm、Δλ=110nm、δλ=0.055nm、共振器長~14mm）。このような波長可変光源は、例えば、米国公開２００９／００５９９７１号に記載されている。
【００７７】
カップラー（スプリッタ）１０４は、光分割器として用いられ、光源１０２から出射された光を測定光（測定光）と参照光に分割する。サーキュレータ１０３はカップラー１０４からの光を光ファイバー１０５に導光し、光ファイバー１０５からの光を光ファイバー１１９に導光する。なお、サーキュレータ１０３は、カップラーであってもよい。
【００７８】
ＯＣＴ光学系１００は、測定光学系１０６によって測定光を眼Ｅの眼底Ｅｆに導く。ＯＣＴ光学系１００は、参照光学系１１０に参照光を導く。ＯＣＴ光学系１００は、眼底Ｅｆによって反射された測定光と参照光との干渉、によって取得される干渉信号光を検出器（受光素子）１２０に受光させる。
【００７９】
測定光学系１０６には、光遅延路３００、光ファイバー１０５、光スキャナ１０８、及び対物レンズ系が順次設けられている。
【００８０】
光遅延路３００は、基準光路３００ａと迂回光路３００ｂを有し、互いに光路長差を持つ少なくとも２つの光を生成させるために設けられている。例えば、測定光路に光遅延路３００が配置された場合、測定光は、光遅延路３００に形成された基準光路３００ａと迂回光路３００ｂによって、互いに光路長差を持つ少なくとも２つの測定光として形成される。迂回光路３００ｂは、基準光路３００ａより光路長が長いので、迂回光路３００ｂを通過する測定光は、基準光路を通過する測定光に対し光学的遅延（光路長差）が生じる。このようにして、光路長差を持つ複数の測定光が被検物の同一位置に照射される。
【００８１】
光遅延路３００は、例えば、第１偏光ビームスプリッタ３０２、第１光反射部材３０４、第２光反射部材３０６、第２偏光ビームスプリッタ３０８を備え、測定光を２つの光路に分割し、一方の測定光に対して他方の測定光の光路長を遅延させる。第１偏光ビームスプリッタ３０２は、光源１０２からの測定光を基準光路３００ａ（第１測定光路）と迂回光路３００ｂ（第２測定光路）に分割する。第２偏光ビームスプリッタ３０８は、基準光路３００ａと迂回光路３００ｂを合成する（図２参照）。
【００８２】
第１光反射部材３０４、第２反射部材３０６として、例えば、全反射ミラー、プリズムなどの光学部材が用いられる。光遅延路３００を形成する光学部材は、図２のように互いに離れた光学配置であってもよいし、プリズム等により一体化された光学配置であってもよい。
【００８３】
第１偏光ビームスプリッタ３０２は、光源１０２からの光を垂直偏光成分と水平偏光成分に分割し、一方の偏光成分の光を透過させ、他方の偏光成分の光を反射する特性を持つ。第１光反射部材３０４、第２反射部材３０６は、第１偏光ビームスプリッタ３０２によって分割された光の一方を反射し、光結合部材３０８に戻す。第２偏光ビームスプリッタ３０８は、垂直偏光と水平偏光に分割された光を結合させる特性を持つ。第２偏光ビームスプリッタ３０８によって結合された後、結果的に、偏光成分が互いに直交する２つの測定光がそれぞれ被検眼に照射される（光ファイバー１０５等の影響により両者の偏光成分は、変更されるかもしれないが、偏光成分が互いに直交するという関係は変わらない）。
【００８４】
上記のようにして光遅延路３００は、偏光成分が互いに直交すると共に光路長差を持つ２つの測定光を生成する。２つの測定光は、サーキュレータ１０３、光ファイバー１０５を介して光スキャナ１０８に向かう。２つの測定光は、光スキャナ１０８によって反射方向が変えられる。光スキャナ１０８によって偏向された光は、対物レンズ系によって平行ビームとなって眼Ｅに入射し、眼底Ｅｆ上に入射される。２つの測定光は、眼底Ｅｆ上の同一位置に照射される。
【００８５】
光スキャナ１０８は、眼底Ｅｆ上でＸＹ方向（横断方向）に測定光を走査させる。光スキャナ１０８は、瞳孔と略共役な位置に配置される。光スキャナ１０８は、例えば、２つのガルバノミラーであり、その反射角度が駆動機構によって任意に調整される。
【００８６】
光源１０２から出射された光束は、その反射（進行）方向が変化され、眼底上で任意の方向に走査される。光スキャナ１０８としては、反射ミラー（ガルバノミラー、ポリゴンミラー、レゾナントスキャナ）の他、光の進行（偏向）方向を変化させる音響光学素子（ＡＯＭ）等が用いられる。
【００８７】
制御器７０は、光スキャナ１０８の駆動を制御することにより、眼底Ｅｆの深さ方向に対して垂直な方向（横断方向）に測定光を走査させる。各測定光の眼底Ｅｆからの後方散乱光（反射光）は、対物レンズ系、光スキャナ１０８、光ファイバー１０５、サーキュレータ１０３、光ファイバー１１９を経て、ビームスプリッタ３５０に達する。そして、後方散乱光は、ビームスプリッタ３５０にて参照光と合波されて干渉する。
【００８８】
参照光学系１１０は、眼底Ｅｆでの測定光の反射によって取得される反射光と合成される参照光を生成する。参照光学系１１０は、マイケルソンタイプであってもよいし、マッハツェンダタイプであっても良い。参照光学系１１０は、透過光学系（例えば、光ファイバー）によって形成され、カップラー１０４からの光を戻さず透過させることにより検出器１２０へと導く。参照光学系１１０は、例えば、反射光学系（例えば、参照ミラー）によって形成され、カップラー１０４からの光を反射光学系により反射することにより再度カップラー１０４に戻し、検出器１２０に導いてもよい。
【００８９】
本装置は、測定光と参照光との光路長差を調整するためにＯＣＴ光学系１００に配置された光学部材の少なくとも一部を光軸方向に移動させる。例えば、参照光学系１１０は、参照光路中の光学部材を移動させることにより、測定光と参照光との光路長差を調整する構成を有する。光路長差を変更するための構成は、測定光路中に配置されてもよい。測定光路中に配置された光学部材（例えば、光ファイバーの端部）が光軸方向に移動される。
【００９０】
ビームスプリッタ３５０は、干渉信号光を２つに分割する。ビームスプリッタ３５０によって分割された光路の一方には、偏光ビームスプリッタ３６０が配置され、他方には、偏光ビームスプリッタ３６５が配置されている。偏光ビームスプリッタ３６０、３６５は、入射された干渉信号光を、互いに直交する偏光成分（垂直偏光成分、水平偏光成分）に分割する。
【００９１】
検出器１２０は、垂直偏光検出器１２０Ｖと水平偏光検出器１２０Ｈを持ち、スペクトル信号における垂直偏光成分と水平偏光成分を別々に検出可能な構成を持つ。
【００９２】
垂直偏光検出器１２０Ｖと水平偏光検出器１２０Ｈは、それぞれ、第１受光素子（１２０Ｖａ、１２０Ｈａ）と第２受光素子（１２０Ｖｂ、１２０Ｈｂ）からなる平衡検出器（Balanced Detector）にて構成されるのが有利である。検出器１２０（平衡検出器）は、第１受光素子からの干渉信号と第２受光素子からの干渉信号との差分を得て、干渉信号に含まれる不要なノイズを削減する。各受光素子は、受光部が一つのみからなるポイントセンサであって、例えば、アバランシェ・フォト・ダイオードが用いられる。
【００９３】
垂直偏光検出器１２０Ｖは、偏光ビームスプリッタ３６０、３６５によって分割された垂直偏光成分を，第１受光素子１２０Ｖａ、第２受光素子１２０Ｖｂにより平衡検出を行う。水平偏光検出器１２０Ｈは、偏光ビームスプリッタ３６０、３６５によって分割された水平偏光成分を，第１受光素子１２０Ｈａ、第２受光素子１２０Ｈｂにより平衡検出を行う。
【００９４】
垂直偏光検出器１２０Ｖと水平偏光検出器１２０Ｈによって受光される干渉信号光は、それぞれ、偏光成分が互いに直交すると共に光路長差を持つ２つの測定光に対応する干渉信号光を含んでいる。
【００９５】
光源１０２により出射波長が変化されると、これに対応する干渉信号光が検出器１２０に受光され、結果的に、スペクトル信号として検出器１２０によって検出される。制御器７０は、光源１０２からのトリガ信号を得て、取得するスペクトル信号と光スキャナ１０８を制御する。
【００９６】
垂直偏光検出器１２０Ｖと水平偏光検出器１２０Ｈによって検出される各スペクトル信号は、被検物に照射された２つの測定光のうち、垂直偏光成分を持つ測定光に基づいて形成された第１スペクトル信号と、水平偏光成分を持つ測定光に基づいて形成された第２スペクトル信号と、を含む。第１スペクトル信号と第２スペクトル信号は、光路長差を持つため、スペクトルによって形成される干渉縞の粗密が異なる。
【００９７】
制御器７０は、偏光成分が異なる２つのスペクトル信号を処理して、互いに直交する偏光成分に関する深さ情報ＤＶ、ＤＨを得る。
【００９８】
制御器７０は、垂直偏光検出器１２０Ｖによって検出された垂直偏光成分を持つスペクトル信号を処理して垂直深さ情報ＤＶを得る。垂直深さ情報ＤＶは、第１スペクトル信号に対応する第１垂直深さ情報ＤＶ１と、第２スペクトル信号に対応する第２垂直深さ情報ＤＶ２を含む。第１垂直深さ情報ＤＶ１は、互いに偏光成分が直交する測定光における一方の測定光に基づいて形成された深さ情報であり、第２垂直深さ情報ＤＶ２は、互いに偏光成分が直交する測定光における他方の測定光に基づいて形成された深さ情報である。
【００９９】
制御器７０は、水平偏光検出器１２０Ｈによって検出された水平偏光成分を持つスペクトル信号を処理して水平深さ情報ＤＨを得る。水平深さ情報ＤＨは、第１スペクトル信号に対応する第１水平深さ情報ＤＨ１と、第２スペクトル信号に対応する第２水平深さ情報ＤＨ２を含む。第１水平深さ情報ＤＨ１は、互いに偏光成分が直交する測定光における一方の測定光に基づいて形成された深さ情報であり、第２水平深さ情報ＤＨ２は、互いに偏光成分が直交する測定光における他方の測定光に基づいて形成された深さ情報である。
【０１００】
＜断層画像の取得＞
制御器７０は、光スキャナ１０８の駆動を制御し、眼底Ｅｆ上で測定光を横断方向に走査させる。制御器７０は、各走査位置での深さ情報を順次並べることにより眼底断層画像を形成させる。
【０１０１】
図３Ａ、図３Ｂは、多重スペクトル信号に基づいて取得された断層画像データを示す例であり、図３Ａは、垂直偏光成分に関する断層画像データであり、図３Ｂは、水平偏光成分に関する断層画像データである。なお、フーリエ解析によって取得された断層画像データには、実像とミラーイメージ（虚像）が含まれるが、図３Ａ、図３Ｂは、実像のみを抽出した画像である。
【０１０２】
制御器７０は、互いに直交する偏光成分に関する深さ情報ＤＶ、ＤＨに関してそれぞれ、走査方向に関して並べることにより互いに直交する偏光成分に関する断層画像データＴＶ、ＴＨを得る。断層画像データＴＶ、ＴＨは、深さ方向に分離された眼底Ｅｆの複数の断層像を含む。なお、断層画像データは、各深さ情報における実虚成分の絶対値を求めることにより形成される。各断層像は、眼底Ｅｆ上の同一の走査位置に関して取得された断層像である。
【０１０３】
断層画像データＴＶは、第１垂直深さ情報ＤＶ１に基づく第１垂直断層像ＴＶ１、第２垂直深さ情報ＤＶ２に基づく第２垂直断層像ＴＶ２を含む。断層画像データＴＨは、第１水平深さ情報ＤＨ１に基づく第１水平断層像ＴＨ１、第２水平深さ情報ＤＨ２に基づく第２水平断層像ＴＨ２を含む。
【０１０４】
制御器７０は、上記のように取得された断層画像データＴＶ、ＴＨから第１垂直断層像ＴＶ１、第２垂直断層像ＴＶ２、第１水平断層像ＴＨ１、第２水平断層像ＴＨ２のいずれかを抽出し、モニタ７５の画面上に断層像を表示する。制御器７０は、断層画像データＴＶ、ＴＨを連続的に取得し、動画の断層像を表示するようにしてもよい。
【０１０５】
＜偏光検出＞
制御器７０は、垂直深さ情報ＤＶと水平深さ情報ＤＨを用いて眼底Ｅｆの複屈折特性を求める。スペクトル信号をフーリエ解析した後の各深さ情報における実部と虚部の情報が用いられる。
【０１０６】
制御器７０は、垂直深さ情報ＤＶから第１垂直深さ情報ＤＶ１を得ると共に、水平深さ情報ＤＨから第１水平深さ情報ＤＨ１を得る。制御器７０は、偏光成分が互いに直交する第１垂直深さ情報ＤＶ１と第１水平深さ情報ＤＨ１に基づいて第１の偏光状態を得る。
【０１０７】
制御器７０は、垂直深さ情報ＤＶから第２垂直深さ情報ＤＶ２を得ると共に、水平深さ情報ＤＨから第２水平深さ情報ＤＨ２を得る。制御器７０は、偏光成分が互いに直交する第２垂直深さ情報ＤＶ２と第２水平深さ情報ＤＨ２に基づいて第２の偏光状態を得る。
【０１０８】
制御器７０は、第１の偏光状態と第２の偏光状態に基づいて、眼底表面を基準として眼底Ｅｆのある位置における複屈折特性を得る。制御部７０は、複屈折特性を深さ方向に関して求めることにより、深さ方向に関する眼底Ｅｆの複屈折特性分布を示す偏光深さ情報を得る。
【０１０９】
制御器７０は、各位置での偏光深さ情報を走査方向に関して並べることにより、ある切断面での眼底Ｅｆの複屈折分布（例えば、偏光深さ情報画像）を求める。制御部７０は、求められた複屈折分布をモニタ７５上に表示する。
【０１１０】
なお、制御器７０は、光スキャナ１０８の駆動を制御し、眼底Ｅｆ上で測定光を二次元的に走査することにより３次元データを得てもよい。制御部７０は、各位置における偏光深さ情報を得ることにより、眼底Ｅｆ上の二次元的な複屈折分布を示すマップを得る。制御部７０は、得られたマップをモニタ７５上に表示する。なお、偏光特性を求める具体的手法については、例えば、特開２０１３−１４８４８２号公報を参照されたい。
【０１１１】
以下、上記実施形態に示した解析処理の一例について説明する。
【０１１２】
PS-OCTデータを得たとする。これは例えば被検眼に照射される２種類の偏光と、被検眼から散乱され戻る際の２つの偏光に相当した４つの３次元OCT断層像である（例えば、前述の第１垂直断層像ＴＶ１の３次元断層像、第２垂直断層像ＴＶ２の３次元断層像、第１水平断層像ＴＨ１の３次元断層像、第２水平断層像ＴＨ２の３次元断層像）。これらをジョーンズマトリクスOCT断層像と呼ぶ。なお、ジョーンズマトリクスOCT断層像から強度OCT像への変換処理が行われることによって、強度ＯＣＴ断層像が取得されてもよい。
【０１１３】
４つのジョーンズマトリクスOCT断層像からは、各種の演算によって、例えば強度OCT断層像、偏光非一様性（DOPU）断層像、複屈折（RetardationまたはBirefringence）断層像、偏光軸回転（Axis-Orientation）断層像といった量が求められる。これ以外にも、ミュラーマトリクス変換により、偏光解消（Depolarization）を検出することも可能である。
【０１１４】
このようなジョーンズマトリクスOCT断層像の信号解析により、少なくとも篩状板に対し開口を持つ二次元面を同定する。これは次のような方法によればよい。
【０１１５】
まず得られたジョーンズマトリクスOCT断層像を、画像処理によってセグメンテーションを行い各層に分離する。その際、上記のジョーンズマトリクスOCT断層像から得られる各種信号分布を援用してもよい。
【０１１６】
例えば網膜色素上皮（RPE）は組織によって偏光状態がスクランブルされるので、偏光解消度（Depolarization）または偏光の非一様性（DOPU）などによって特徴的に検出されることが知られている。RPEは信号強度が高いため精度が高く検出できる点で有利である。
【０１１７】
例えば、図４のように強度OCT断層像と、DOPU断層像の各点での信号強度を２値化処理によって特定の閾値の上下で分離することで、図５のような断層像を得ることにより、強度OCT断層像上でRPEを層別することができる（図６）。
【０１１８】
或いは、OCT以外の装置例えば眼底カメラなどとの相関を重視して、ブルッフ膜の検出によってもよい。眼底カメラでは視神経乳頭周辺のブルッフ膜は境界が明瞭な白色信号として映り、色の境界部分として開口を認識しやすい。すなわち、OCT像で検出された開口部分と眼底カメラで検出された開口部分とが一致するようになり、しかし通常ＯＣＴ画像では、ブルッフ膜とRPEはいずれも高輝度で表示されるので、強度OCT像からは分離し難い。これを例えば、高輝度領域を強度OCTから閾値処理などによって求め、前記DOPUによってRPEを分離し、残った部分をブルッフ膜として検出するのでも良い。
【０１１９】
或いは、後に述べるように、強膜は複屈折性を持つので強膜を検出し他との境界を求めてもよい。通常、隣り合うRPEが明瞭に映るので、強膜を選択的に特定しなくても層別が行えるが、例えば網膜色素変性症の患者においてはRPEが消失している場合もあり、その際は複屈折断層像を用いるとよい。
【０１２０】
強度OCT断層像から、公知のセグメンテーション処理手法によって層別する方法でも勿論よい。このような方法は、例えば、特開2016-055122に記載された手法が用いられてもよい。或いは複数枚数の断層像を撮影した場合は、それらに対し統計処理をし、信号の強度分布の統計による層別を行っても良い。すなわち、位置を固定で信号強度のヒストグラムをとると、ノイズの影響で信号はある強度の周りに正規分布をするが、異なる組織ではその強度ないしは正規分布の偏差（分散と言い換えても可）が異なっていれば、ある閾値を設けてこれらを分離することが出来る。勿論、これらの方法に限定されず、例えば上記方法の組合せであってもよい。
【０１２１】
このように、少なくとも２つの層が検出できれば、その境界部分として１つの２次元面を定義することができる。そしてそれは、少なくとも網膜神経乳頭周辺には網膜層が存在しないので、篩状板に対し開口を持つことになる。病態によって、特定の層は歪んだり、消失したりすることがあるため、例えば治療前後で変化しないものや、明瞭に検出できるもの、などの基準から、これらは任意に選んでよい。
【０１２２】
深さ方向をZ、それと垂直な方向をそれぞれX、Y方向と定め、同定した面の深さZがX、Yに対して如何に変化するかの関数Z=U(X、Y)として、これを表す。このとき、ノイズの影響を除くため、境界が滑らかに繋がるように適宜スムージング処理を行ってもよい。例えば、X、Yそれぞれで相隣り合う両側の画素同士との間で平均を取り、これを中央の画素での値とすることで実現してもよいし、同じ箇所の信号を多数回撮影し、平均化してもよい。これを３次元的に表示すると、例えば図７のようになる。
【０１２３】
次に、検出された特定の面は少なくとも篩状板に対し開口しているので、開口部分を滑らかに補間するような仮想面を考え、これを元に基準を定める。このとき、再現性を得るために一意的に形状が決まることが望ましく、これは例えば、開口部分の形状を境界条件とする次のポアッソン方程式
【数１】
2018068778.tif 000003
の解として得ることで実現できる。境界条件を決めたときにポアッソン方程式の解が一意であることは、微分方程式の文献（例えば、James Brown and Ruel Churchill, Fourier Series and Boundary Value Problems, 8th Edition (McGraw-Hill Companies, Inc., New York, NY, 2011)）などに記載されており、これはUが極大値と極小値を持たない条件と等価である。
【０１２４】
そしてこれを求める場合には、開口部分をグリッドに区切り、差分方程式の形で数値的に解く方法が考えられる。すなわち、深さ方向にZ軸を取り、面の表面形状の分布をX、Yの値によってZ=U(X、Y)として記述するものとする。開口部分を一辺がΔsの正方形のグリッドに分割することで、開口内の各点での高さはZi,j＝U(Xi、Yj)として指定できる（図８、図９）。差分方程式で記述したポアッソン方程式は、次の形となる。
【数２】
2018068778.tif 000004
或いはZi,jについてこれを解いて
【数３】
2018068778.tif 000005
としてもよい。
【０１２５】
検出した開口部における境界でのZの値を与えると、相隣り合うグリッドでのZの高さに関する条件が順次求まっていくので、これにより点の数ぶんの連立1次方程式が得られる。多変数の連立1次方程式を解く方法は多数知られているが、例えば行列を用いたガウス・ザイデル法や、変分を元に誤差を小さくしていくSOR法などが知られており、これにより各点での形状を決めることができる。勿論、グリッドの切り方もこれに限られることはなく、これ以外の方法によっても良い。例えば、XとYとでグリッドの辺の長さを変えたりしてもよい。
【０１２６】
結果として、図１０に示すような開口内での補間形状分布が得られる。ただし縦方向のスケールは変更し拡大して描いてある。これは例えば２次元断層像に反映して描けば、図１１のような分布となる。
【０１２７】
更に開口での重心を求めるには、次の方法によっても良い。すなわち、適当な選び原点O（X、Y）＝（0，0）とし、前記グリッド内の全ての点について（Xi、Yi）の平均をとったものをG（Xmean、Ymean）として求める。これが重心Ｇであり、これを改めて原点と選び直せば、各点の位置は重心を中心にして再定義したものとなる（図１２参照）。
【０１２８】
例えば固視点を変更できるOCT装置の場合、特定の固視位置に被検眼を誘導して同じ位置で撮影が行えるため、Z軸の方向をOCTのビームの軸方向（Aスキャン方向）と一致しているとして定義してもよいが、例えば上記で求めた補間された面の法線を算出し、X、Y軸とともにZ軸も再定義すれば更によい。例えば上記で検出した原点におけるX、Y方向の傾き
【数４】
2018068778.tif 000006
【数５】
2018068778.tif 000007
を求め、（-fx、-fy）として求めてもよい。前述の断層像が全てこの軸を含む内にあったとして、この軸を反映して描いたものは図１３のようになる。以降は全て、この軸を基準にした座標系で議論を進めることにする。
【０１２９】
このような軸Ｊの取り直しによって画素そのものも影響を受けるため、例えば任意の座標におけるOCT信号値が、軸の取り直し前での二つ以上の画素に部分的に共通部を持ってしまうことが起こりうる。その時も、重なり量に応じた重み付け平均などの処理を求めることで、その点での物理量を近似しうる。
【０１３０】
ところで、NFLの周囲でのOCT信号により求めた物理量は、NFL周囲の組織の機械強度とよく相関することが知られている。例えば、非特許文献[http://www.arvo.org/webs/am2016/sectionpdf/GL/Session_143.pdf]などにその記載があり、例えば篩状板に関する複屈折断層像の値はその例である。
【０１３１】
篩状板を選択的に検出するために、OCT強度断層像によってセグメンテーションを行っても良いし、上記セグメンテーションの部分で述べたのと同様の方法を用い、複屈折断層像を援用しても良い。これは篩状板が繊維状の組織であるため、比較的大きな複屈折性を持つことによる。複屈折断層像をある閾値で2値化すると図１４のような分布が得られる。篩状板以外の領域にも分布してしまうのは、強膜も同様の複屈折性を持つためであるが、強度OCT断層像と合わせ二つの領域の区別は容易に行える（図１５）。すなわち、高輝度に検出されるRPEと接触し連続している領域は強膜であるし、前記補間によって求めた面よりも深い位置に、離れて存在しているのが篩状板であるとして区別できる。或いは、先に定義した軸と、高複屈折値の位置の交点を求めそこから検索をスタートし、中空の領域を複数持つ連続した領域として検出してもよい。結果として、図１６のような分布が３次元的に求められ、これを篩状板として採用してもよい。
【０１３２】
そしてこのように篩状板は元来偏光依存性を持つが、健常眼であればその分布は上記で定めた中心軸に対しほぼ等方的であるのに対し、篩状板が圧迫されるなどにより篩状板に負荷が加わると、この分布には非等方性が現れてくる。負荷によって生じる歪を、前眼部に複屈折断層像を用いることで検出した例には、例えば次の非特許文献[
201610311916231980__________P___________________________APH_5
. 2014 May 1; 5(5): 1391-1402.  Yamanari et al.(2014)]があるが、上記による分布も、篩状板が応力を受けることでの歪を可視化したものと考えることができる。
【０１３３】
先に定めた軸に対し、動径方向、周方向への分布として複屈折断層像の分布の非等方性を求めることで、上記の非等方性がどのような応力の掛かり方によるものなのかを特定することができる。これは例えば、先に定めた軸に対して複屈折断層像の３次元分布がXi、Yj、ZkにおいてRi,j,k＝B(Xi、Yj、Zk)  で求められているとすると、円筒座標への変換
【数６】
2018068778.tif 000008
【数７】
2018068778.tif 000009
によって、各点での複屈折断層像の分布R(r、φ、Zk)を決定することができる。これを、例えば篩状板の軸方向への所定の厚みについて平均化すると、例えば図１７のような分布が得られる。各方向と値は１にスケールしてあるとする。ここで分布が滑らかでないのは、測定信号にはノイズが乗るためで、これは偏光を用いたOCTでは特に顕著である。
【０１３４】
従って、これを多項式フィットして扱いやすいようにしてもよい。多項式フィットの方法としては、円柱座標系での表示に適したものとしてツェルニケ多項式があり、各次数のツェルニケ多項式と上記分布との積を開口面内で積分することで、重ね合わせの係数を求めるようにしてもよい。例えば36次までなど、高次成分を除いてフィットさせることで、スパイクノイズの影響を小さくすることができる。結果として、図１８のような滑らかな分布にすることができる。中心すなわちｒ＝0の近傍は、相対的に少ない画素からの情報となってしまうので、これを除いてフィッティングを行ってもよい。
【０１３５】
鳥瞰図によって示したが、図１９の如く２次元の等高線プロットによって示してもよい。図１９と同じ分布を、白色ほど値が大きく、黒色ほど値が小さくなるようにしている。
【０１３６】
上記の例では、物理量として篩状板の複屈折断層像の値としたが、これに限るものではない。例えば乳頭周辺強膜の血管像の３次元強度であってもよいし、その密度分布であってもよい。これは、NFL周辺が圧迫により強膜形状が歪むことで、例えば血管の凝集度が上がったり、血流の通りやすさが変わったり、或いはその結果として、局所的に新生血管が発生したりするためである。
【０１３７】
強膜領域を選択的に検出するためには、OCT強度断層像によってセグメンテーションを行っても良いし、上記セグメンテーションの部分で述べたのと同様の方法を用い、複屈折断層像を援用しても良い。先の篩状板の検出で述べた方法により、篩状板を先に検出して、残った領域として検出してもよいし、DOPU断層像によりRPEを選択的に検出し、これと隣接する連続領域として検出してもよい。
【０１３８】
強膜領域において、例えば先に定めた軸周りに、特定の半径（1mmなど）内で血管の分布を求めることで、中空円筒様内の領域での血管の分布A（ｒ、φ、Z）が求められる。これを動径方向位置と周方向位置を定めて深さ（Z）方向の総強度を求めることによって、すなわち
【数８】
2018068778.tif 000010
を行って、ｒとφだけの関数としてもよい。ここに、Zmin、Zmaxは強膜信号のある上下限であり、上式は積分で書いているが、実用上これを離散和に直して計算してもよい。また、方向のみが問題である場合には、これを更にｒ方向にも積分してφだけの関数としてもよい。例えば更に周方向を複数に分割し、図２０のようにこの領域に枝分かれした血管１つしかない場合にてこれを計算すれば、この中空円筒様領域を深さ方向に向けて見た時の図２１のような分布が得られ、枝分かれ部で血管の密度が高いという計算結果を得ることができる。
【０１３９】
血管の同定は、例えば同じ箇所を複数回撮影し、求めた振幅や位相あるいはその両方を考慮してOCT信号同士の差分をとっても良いし、複数回撮影しそのスペックルパターンの統計解析によってもよい。血流がある箇所ではスペックルパターンが異なるので、例えばその分散を調べて閾値処理をすることで、血管のみを同定できる。
【０１４０】
強膜領域では、ここに述べた血管の他に、偏光解消度を見ても良い。これは組織が圧迫を受けて繊維構造が乱れることで、偏光がスクランブルされる度合いが大きくなるためと考えることができる。考える領域の定め方、平均の取り方などは上記の血管の場合と同様であるので省略する。
【０１４１】
また、物理量としては、篩状板の形状そのものに注目してもよい。負荷が加わると、篩状板が変形して非等方な反りを生じることが知られている。例えば、[Lamina Cribrosa Morphology in Glaucoma  KI-HO PARK（2016．ISER）]を参照。
【０１４２】
篩状板の表面形状は、先に示したセグメンテーションの方法と同様に検出してもよい。検出後、ノイズの影響を除くためにスムージングを施してもよい。多項式フィットによるものでも勿論よい。特に篩状板は表面に凹凸が多数あるため、反りのようなマクロな変形を問題にするときはフィッティングの効果が大きい。或いは、篩状板の厚みに注目してもよい。篩状板は力が加わることで薄くなる。
【０１４３】
このようにして複数の物理量のｒ、φ依存性を得たが、ｒが0から最大値まで変化するのに伴って物理量が同じ量だけ変化する場合であっても、ｒに比例して大きくなるのと、ｒが小さい時は緩やかな変化で、ｒが大きくなると急激に大きくなるのとでは応力の掛かり方は異なっていると考えられる。
【０１４４】
従って、これら求められた物理量をｒ方向に微分した、「変化率」に注目してもよい。同様の理由から、応力は局所的に曲率が大きいほど大きいと考えられ、物理量の「変化率」を更にｒ方向に微分した、変化の「加速度」に注目してもよい。
【０１４５】
また更に、上記では平均する操作を用いたが、これに限定されない。極めて大きい振幅のスパイクノイズが発生する場合は、これを除去する目的でメディアンを用いた方がよい。或いは、メディアン操作と平均操作の中間、すなわち最小と最大から数点を除いての平均、といった操作でもよい。
【０１４６】
また、例えば上記で多項式フィットした際の、実測とフィッティング形状との間の乖離（分散）に着目するのもよい。特定の方向で分散が大きいということは、そこに特異な性質があるということであって、例えば極めて大きな変化が存在する場合などが考えられる。同様の理由から、最大値と最小値に注目するのもよい。
【０１４７】
或いは、方向によって多項式のフィットの良し悪しが異なる場合は、領域を複数のエリアに区切り、フィッティングの関数系を変更してもよい。例えばフーリエ級数展開によって近似してもよいし、その他の関数系、例えばルジャンドル陪多項式などを用いてもよい。
また、例えば血管検出のときなどに、表面の反射成分が乗ってしまい、方向によっては偽信号がバイアスされてしまうことがありうる。これを除く目的で、PV値＝（最大値―最小値）に注目するのでも良い。このような相対値とすることで、バイアス同士が相殺し影響を取り除くことができる。或いは、各方向の値から最小値やメディアンなどをオフセットすることで、同様の効果を得ても良い。
【０１４８】
例えば早期の緑内障の病変を診断するにあたっては、構造の変化に先んじて機能の変化が生じてしまうため、早期診断が必要とされている。従って、できるだけ多くの情報を、高い確度を持って得ることが重要視される。例えば形状の変化に注目しただけでは変化が微小であった場合にこれを見逃す可能性が高い。例えば強度像の形状に着眼した特許文献[特開2013-153844]はBowing Angleと呼ばれる角度情報に注目したものであるが、一般にはOCTの分解能以上の情報を得ることができない。
【０１４９】
一方で偏光OCTを用いて、位相レベルの情報量を多くしたとしても、一般に偏光OCTではノイズの影響を受けやすく、絶対値による診断の精度すなわち確度という観点からは課題が大きい。
【０１５０】
これを領域ごとの相対値比較つまりは方向による変化のパターンに注目し比較することで、それぞれの欠点をお互いに相殺し、信頼に足る結果を得ることができる。
【０１５１】
以下に、例を用いて説明する。例えば正常眼圧緑内障（NTG）と高眼圧緑内障（HTG）ではNasal方向、Temporal方向、Superior方向、Inferior方向における影響の現れ方が異なることが知られている。例えば篩状板の表面形状を見ると、NTGではSuperior-Inferior方向で陥凹が大きく、Nasal-Temporal方向で陥凹が小さい。HTGではNasal-Temporal方向、Superior-Inferior方向いずれも陥凹が大きい。すなわち、正常眼と、NTG眼、HTG眼を比較すると、図２２〜図２４に示すような、形態の変化が生じている。尚、スケールは誇張して描いている。
【０１５２】
従って、それぞれについて等高線表示すると表面形状は図２５〜図２７のようになり、それぞれの様態に特徴的な変形のパターン・異方性があることが判る。
【０１５３】
ここで、一般に偏光OCTの結果はノイズに弱いが、これは各点ごとに絶対的な値を求めようとした場合であって、場所ごとの相対値を求めるならば比較的正確な結果が得られるという点に着目する。値のバラツキは解の不安定性と言われるもので、値の小さな数値による逆行列演算などを繰り返すことで生じるものだからである。
【０１５４】
先に述べた通り、複数の物理量を測定した結果のこのようなパターンに注目することで、より確度の高い結果を得ることができる。例えば、偏光により複屈折断層像を得て、これも同様に動径方向と周方向とで分布を求めたとする。これは応力による歪を反映したものとなるが、ここで任意の物体に機械的な力が働く時、形状変化が急峻に起こった箇所ほど応力が集中していることは周知の事実である。従って、例えば複屈折断層像より図２８の複屈折分布を得たとすると、先の篩状板の形状変化を動径方向に２回微分した「局所的な曲率」とのパターンの一致を見ることで、結果がよく説明され、信頼性が増す。
本実施例では、図２９に示すNTGの場合の２回微分と比較することで、Superior-Inferior方向で変化が大きく、Nasal-Temporal方向で値が平坦であるというパターンがよく一致することが判る。
【０１５５】
従って、先の篩状板表面形状が同じくNTGのものに近ければ、複屈折断層像による結果は補強され、表面形状だけから議論をするよりも高い確度でNTGということができる。一方で、先の篩状板表面形状が正常眼やHTGのものに近いのに、複屈折断層像から同結果を得たとすれば、篩状板表面形状だけからNTGであると判断するのは危険であると判断することが可能である。
【０１５６】
上記は形状の局所変化と複屈折断層像について論じたが、もちろんこれに限定されない。圧迫されている箇所は強膜の血流の大小が異なると考えられ、同様にSuperior-Inferior方向で変化が大きく、Nasal-Temporal方向で値が平坦であるというパターンが現れるならば確度は更に高いということになる。
【０１５７】
同様に、圧迫により線維走行が乱れることを考え、DOPU断層像によって強膜での偏光のスクランブルの様子を見て、これもSuperior-Inferior方向で変化が大きく、Nasal-Temporal方向で値が平坦であるというパターンが現れるならば確度は更に高いということになる。
【０１５８】
このように結果のパターンの一致を次々と見ていくことで、診断の確度を論じることができ、確度を持って正常／病変の診断をすべきなのか、確度が低いため再度撮影を行った方が良いかの判断尺度とすることができる。
【０１５９】
なお、両者の結果は完全には一致しない。より定量的に一致度を見るならば、次のようにして類似度の大小を論じても良い。
【０１６０】
例えば領域をｎ＝１、２、・・・、Nとし、それぞれの領域での二つの物理量をPn、Qnとすると
【数９】
2018068778.tif 000011
なる計算で相関係数は求められる。これは完全に一致している場合を１に規格化してあるため、この一致度をそれぞれの物理量に対して求め、その結果に応じて結果の信頼度に反映してもよい。例えば複数の結果同士Γの積を取る方法が考えられるが、もちろん、相関の度合いによって重み付けを変えても足し合わせるのでもよい。
【０１６１】
動径方向の中心は特異な点となりやすいため、これを除いて議論しても良い。また、図は円形の領域において全て信号が得られている（ケラレがない）場合を想定したが、このような状況は稀であり、実際には一部がケラれていたりする場合もある。従って、少なくとも信号が得られている領域の内接円と、中心から所定の割合を除いた中抜け領域（図３０）において、一致度を議論するのもよい。
【０１６２】
更に、上記は周方向の分割数は一般の場合を議論したが、注目するのはパターン・異方性であるので、これらの物理量を、ｒ＝０の軸を中心に、例えばNasal、Temporal、Superior、Inferiorの４領域に分けて、代表値にて取り扱うのもよい。勿論、４領域に限定する必要はないが、４領域程度が扱いやすい。例えばNasal方向への水平断方向をφ＝０として、-45°〜45°をNasal領域、45°〜135°をSuperior領域、-45°〜-135°をInferior領域、残りの領域をTemporal領域と定めてもよい。
【０１６３】
そして代表値としては、領域内に含まれる全ての値の平均や、メディアン、分散値、最大・最小値、PV値のいずれであってもよい。
【０１６４】
更に、前述のように方向によりこれらの物理量の変化の仕方は異なっている。物理量は多項式によるフィッティングを行うのが良いが、ある方向には一致が良いが、別の方向では一致が悪いという状況が生じうる。そのため、フィッティングの度合いが良くなるように、各領域でフィッティングの関数を変えるのもよい。緑内障の場合は、図２２〜２６のような変化を生じるので、Superior-Inferior方向は極値を持つ関数で、Nasal-Temporal方向は極値を持たない関数でフィッティングするのもよい。
【０１６５】
このように、パターンが複数の物理量の間で類似するかどうかを見ることで上記結果が理に適ったものかどうかを論じることができ、診断に役立てることができる。これは強度OCT断層像だけからでは不足した情報量を補い、かつPS-OCTの難点であった確度の低さを克服するものである。
【０１６６】
また、篩状板の検出においては再現性が問題であったが、上記によれば３次元像から一意的に定まる面を定義して軸を決定するため、再現性よく上記を行うことができる。
【０１６７】
なお、ＰＳ−ＯＣＴ装置としての光学系としては、特開２０１３−１４８４８２号公報に記載された光学系を用いたが、これに限定されず、被検眼の偏光特性を示すＯＣＴ画像データを取得可能なＰＳ−ＯＣＴ装置であれば、これに限定されない。
【０１６８】
なお、本実施形態の装置を、以下の方法を実施するために用いてもよい。
【０１６９】
（１）偏光を用いた少なくとも２断面からなる３次元眼科用OCTの解析方法であって、OCTの信号から、
篩状板と、
篩状板よりも浅い位置で、少なくとも篩状板に対し開口している２次元面と、
を同定する手段と、
該開口部を補間することで得られる曲面を生成する処理部と、を持ち、
前記処理部は該曲面上の開口部の重心を求める演算処理を行い、眼底における物理量の
前記重心を中心とした動径方向と周方向への変化を解析し表示することを特徴とする。
【０１７０】
（２）  物理量は、強膜の血管像の強さないし血管像の３次元密度、または特定層の血管像の密度であることを特徴とする（１）記載の方法
（３）  物理量は、強膜の偏光解消度であることを特徴とする（１）記載の方法。
【０１７１】
（４）  物理量は、該曲面から断層像の深さ方向に測った篩状板の表面との距離差ないしその動径方向への変化率、ないし変化の局所曲率であることを特徴とする（１）記載の方法。
【０１７２】
（５）  物理量は、篩状板に入射する光線の偏光遅延量差ないしその動径方向への変化率、ないし変化の局所曲率であることを特徴とする（１）記載の方法
（６）  篩状板ないし少なくとも篩状板に対し開口している２次元面の同定が、
強度信号と入射光線の偏光遅延量差の解析によることを特徴とする（１）〜（５）記載の方法
（７）  篩状板ないし少なくとも篩状板に対し開口している２次元面の同定が、強度信号と強膜領域の入射光線の偏光解消度の解析によることを特徴とする（１）〜（５）記載の方法
（８）  前記開口部の補間は、前記開口部を境界として表面形状が極小値ないし極大値を
持たない条件から決定されることを特徴とする（１）〜（７）記載の方法
（９）  前記重心を中心とした少なくとも４つの領域（Nasal、Temporal、Superior、Inferior）に分けて前記４つの領域における物理量の代表値を表示することを特徴とする（１）〜（８）の方法。なお、ここで代表値とは、前記重心から等距離の円内ないし二つの円に挟まれた領域における平均値、メディアン、分散値、最大値、最小値、PV値  のいずれかである。
【０１７３】
（１０）  前記重心を中心とした少なくとも４つの領域（Nasal、Temporal、Superior、Inferior）に分けて前記４つの領域における物理量の代表値を複数算出し、前記複数の代表値同士で周方向に関する変化パターンが一致しているかに応じて信頼性に関する数値を増減し表示に反映することを特徴とする、（１）〜（９）の方法。
【０１７４】
（１１）  得られた篩状板に関する物理量を、多項式にてフィットして解析することを特徴とする、（１）〜（１０）の方法。
【０１７５】
（１２）  前記重心を中心とした少なくとも４つの領域（Nasal、Temporal、Superior、Inferior）に分けて解析する際に、少なくともNasal、Temporalの領域とSuperior、Inferiorの領域とではフィッティングの関数系を異なるものとすることを特徴とする、（１１）の方法。
【０１７６】
（１３）  前記フィッティングの関数系は、Nasal、Temporalの領域では動径方向に極値を持たずSuperior、Inferiorの領域では動径方向に極値を持つものを含むことを特徴とする、（１２）の方法。
【０１７７】
なお、（２）〜（１４）に記載された方法は、必ずしも（１）の方法に限定されるものではなく、独立して実施されてもよい。また、重心に関して篩状板上に設定された基準位置であれば、必ずしも上記重心に限定されない。
【図面の簡単な説明】
【０１７８】
【図１】本実施例に係る眼科ＯＣＴ解析装置の一例を示す図である。
【図２】本実施例に係るＯＣＴ装置１の一例を示す図である。
【図３Ａ】垂直偏光成分に関する断層画像データの一例を示す図である。
【図３Ｂ】水平偏光成分に関する断層画像データの一例を示す図である。
【図４】強度ＯＣＴ断層像の一例を示す図である。
【図５】DOPU断層像にてセグメンテーションを行った場合の一例を示す図である。
【図６】強度OCT断層像上でRPEを層別した場合の一例を示す図である。
【図７】開口部を３次元的に表示した場合の一例を示す図である。
【図８】グリット分割を行った場合の一例を示す図である。
【図９】グリット分割を行った場合の一例を示す図である。
【図１０】開口内での補間形状分布が得られた場合の一例を示す図である。
【図１１】補間形状分布を２次元断層像に反映した場合の一例を示す図である。
【図１２】開口の重心Ｇの一例を示す図である。
【図１３】開口に基づく軸Ｊの一例を示す図である。
【図１４】複屈折断層像を二値化した場合の一例を示す図である。
【図１５】複屈折断層像を強度OCT断層像と合わせた場合の一例を示す図である。
【図１６】篩状板を検出する場合の一例を示す図である。
【図１７】篩状板の軸方向への所定の厚みを示すマップの一例を示す図である。
【図１８】図１７のマップに関してフィッティングを行った場合の一例を示す図である。
【図１９】２次元の等高線プロットにて篩状板の軸方向への所定の厚みをマップにて示す場合の一例を示す図である。
【図２０】強膜領域において血管の分布を求める場合の一例を示す図である。
【図２１】中空円筒様領域を深さ方向に向けて見た時のの一例を示す図である。
【図２２】正常眼での篩状板の表面形状の一例を示す図である。
【図２３】正常眼圧緑内障での篩状板の表面形状の一例を示す図である。
【図２４】高眼圧緑内障での篩状板の表面形状の一例を示す図である。
【図２５】正常眼での篩状板の表面形状の一例を示す図である（等高線表示）。
【図２６】正常眼圧緑内障での篩状板の表面形状の一例を示す図である（等高線表示）。
【図２７】高眼圧緑内障での篩状板の表面形状の一例を示す図である（等高線表示）。
【図２８】複屈折断層像より複屈折分布を得た場合のパターンの一例を示す図である。
【図２９】NTGの場合の一例を示す図である。
【図３０】中抜け領域について説明するための図である。
【符号の説明】
【０１７９】
１  ＯＣＴ装置
１０  ＯＣＴ解析装置
７０  制御部
７５  表示部
【図１】
2018068778.tif 000012
【図２】
2018068778.tif 000013
【図３Ａ】
2018068778.tif 000014
【図３Ｂ】
2018068778.tif 000015
【図４】
2018068778.tif 000016
【図５】
2018068778.tif 000017
【図６】
2018068778.tif 000018
【図７】
2018068778.tif 000019
【図８】
2018068778.tif 000020
【図９】
2018068778.tif 000021
【図１０】
2018068778.tif 000022
【図１１】
2018068778.tif 000023
【図１２】
2018068778.tif 000024
【図１３】
2018068778.tif 000025
【図１４】
2018068778.tif 000026
【図１５】
2018068778.tif 000027
【図１６】
2018068778.tif 000028
【図１７】
2018068778.tif 000029
【図１８】
2018068778.tif 000030
【図１９】
2018068778.tif 000031
【図２０】
2018068778.tif 000032
【図２１】
2018068778.tif 000033
【図２２】
2018068778.tif 000034
【図２３】
2018068778.tif 000035
【図２４】
2018068778.tif 000036
【図２５】
2018068778.tif 000037
【図２６】
2018068778.tif 000038
【図２７】
2018068778.tif 000039
【図２８】
2018068778.tif 000040
【図２９】
2018068778.tif 000041
【図３０】
2018068778.tif 000042