(57)【要約】【課題】ＯＣＴ計測データから、脈絡膜の層内に存在する脈絡膜血管網のみの画像データを選択的に抽出し、該画像データから脈絡膜血管の太さや、脈絡膜血管網の厚さを正確に得て、脈絡膜血管網の定量評価を可能とする。【解決手段】光干渉断層計装置は、光干渉断層計と、光干渉断層計で取得したＯＣＴ計測データに基づいて得た三次元のＯＣＴ断層画像の画像処理を行うコンピュータと、を備え、コンピュータは、三次元のＯＣＴ断層画像から脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段、として機能する。【選択図】図３
2015000131.tif 000002
【特許請求の範囲】
【請求項１】
光干渉断層計と、該光干渉断層計で取得したＯＣＴ計測データに基づき三次元のＯＣＴ断層画像を得て、該三次元のＯＣＴ断層画像の画像処理を行うコンピュータと、を備えた脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置であって、
前記コンピュータは、前記三次元のＯＣＴ断層画像から脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段、として機能し、
前記脈絡膜の血管の画像データを取得する手段は、ＯＣＴ計測データから脈絡膜の断層画像データを抽出し、該脈絡膜の断層画像データから脈絡膜における深さ方向の等間隔ごとの位置でスライスした画像のデータを抽出して、該スライスした画像のデータから脈絡膜の血管の画像データを得る構成であることを特徴とする脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置。                                              
【請求項２】
前記脈絡膜の血管の画像データを取得する手段は、大きさの異なる複数のウィンドウのそれぞれについて、                            
前記スライスした画像の画素の色濃度が、予め定めた所定の閾値以上であるか否かに応じて画素毎に２値化し、血管推定部分が抽出された２値画像を得て、
該血管推定部分が抽出された２値画像において、当該ウィンドウ内における各画素に対する２値データの異なる画素の存在する割合が、予め定めた所定値以上の領域は擬血管として消去するとともに、当該ウィンドウの寸法に対する予め定めた所定の大きさの直径より小さな直径の前記血管推定部分を、ノイズであり非血管として消去して、脈絡膜の血管網の２値画像のデータを得る構成であることを特徴とする請求項１に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置。
【請求項３】
前記脈絡膜の血管の画像データを取得する手段は、前記三次元のＯＣＴ断層画像から眼底におけるブルッフ膜の前側の領域の深さ方向の勾配を検出することで、眼底における視細胞との内節外節結合部の位置のデータを得て、該内節外節結合部及びブルッフ膜の位置のデータを使用して網膜色素上皮周辺の高反射構造体を抽出し、該高反射構造体における光強度を平均することによって、高反射構造体における網膜血管による陰を特定する画像データを得て、該画像データを、線状を強調するフィルターであるフランジフィルターにかけて得られた画像データを２値化し、網膜血管の陰の２値画像のデータを得る構成であることを特徴とする請求項２に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置。
【請求項４】
前記脈絡膜の血管の画像データを取得する手段は、前記脈絡膜の断層画像データに基づき、血管の部分に相当する全ての画素について、光強度の大小により、中小の血管と大きな血管に分類した分類データを作成し、該分類データに基づき、大きさの異なる複数のウィンドウのそれぞれについて得た前記脈絡膜の血管網の２値画像のデータを選択的に組み合わせるとともに、前記網膜血管の陰の２値画像のデータに基づき、脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータを形成する構成であることを特徴とする請求項３に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置。
【請求項５】
前記脈絡膜の血管の形状の定量評価を行うためのデータを得る手段は、前記脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータに基づき、脈絡膜の血管の直径を見積もり、該見積もりで得られたデータを基に脈絡膜の血管の直径分布データを作成するとともに、脈絡膜血管網の厚さを計測し、該計測で得られたデータを基に脈絡膜血管網の厚さマップを作成する構成であることを特徴とする請求項４に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置。
【請求項６】
光干渉断層計と、該光干渉断層計で取得したＯＣＴ計測データに基づき三次元のＯＣＴ断層画像を得て、該三次元のＯＣＴ断層画像の画像処理を行うコンピュータと、を備えた脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラムであって、
前記コンピュータを、前記三次元のＯＣＴ断層画像から脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段、として機能させるものであり、
前記脈絡膜の血管の画像データを取得する手段は、ＯＣＴ計測データから脈絡膜の断層画像データを抽出し、該脈絡膜の断層画像データから脈絡膜における深さ方向の等間隔ごとの位置でスライスした画像のデータを抽出して、該スライスした画像のデータから脈絡膜の血管の画像データを得ることを特徴とする脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラム。
【請求項７】
前記脈絡膜の血管の画像データを取得する手段は、大きさの異なる複数のウィンドウのそれぞれについて、                            
前記スライスした画像の画素の色濃度が、予め定めた所定の閾値以上であるか否かに応じて画素毎に２値化し血管推定部分が抽出された２値画像を得て、
該血管推定部分が抽出された２値画像において、当該ウィンドウ内における各画素に対する２値データの異なる画素の存在する割合が、予め定めた所定値以上の領域は擬血管として消去するとともに、当該ウィンドウの寸法に対する予め定めた所定の大きさの直径より小さな直径の前記血管推定部分を、ノイズであり非血管として消去して、脈絡膜の血管網の２値画像のデータを得ることを特徴とする請求項６に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラム。
【請求項８】
前記脈絡膜の血管の画像データを取得する手段は、前記三次元のＯＣＴ断層画像から眼底におけるブルッフ膜の前側の領域の深さ方向の勾配を検出することで、眼底における視細胞との内節外節結合部の位置のデータを得て、該内節外節結合部及びブルッフ膜の位置のデータを使用して網膜色素上皮周辺の高反射構造体を抽出し、該高反射構造体における光強度を平均することによって、高反射構造体における網膜血管による陰を特定する画像データを得て、該画像データを、線状を強調するフィルターであるフランジフィルターにかけて得られた画像データを２値化し、網膜血管の陰の２値画像のデータを得ることを特徴とする請求項７に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラム。
【請求項９】
前記脈絡膜の血管の画像データを取得する手段は、前記脈絡膜の断層画像データに基づき、血管の部分に相当する全ての画素について、光強度の大小により、中小の血管と大きな血管に分類した分類データを作成し、該分類データに基づき、大きさの異なる複数のウィンドウのそれぞれについて得た前記脈絡膜の血管網の２値画像のデータを選択的に組み合わせるとともに、前記網膜血管の陰の２値画像のデータに基づき、脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータを形成することを特徴とする請求項８に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラム。
【請求項１０】
前記脈絡膜の血管の形状の定量評価を行うためのデータを得る手段は、前記脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータに基づき、脈絡膜の血管の直径を見積もり、該見積もりで得られたデータを基に脈絡膜の血管の直径分布データを作成するとともに、脈絡膜血管網の厚さを計測し、該計測で得られたデータを基に脈絡膜血管網の厚さマップを作成することを特徴とする請求項９に記載の脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光干渉計測技術に分野であって、光干渉断層計（光コヒーレンストモグラフィー：Optical Coherence Tomography、略して「ＯＣＴ」という）で得られた画像に対する画像処理の技術に関し、特に、特に、目の眼底における脈絡膜の血管の可視化し解析する光干渉断層計装置及びその画像処理プログラムに関する。
【０００２】
上記本発明の画像処理プログラムは、概要としては、３次元のＯＣＴ断層画像を解析し擬似的な眼底造影画像を構築する擬似造影画像作成のためのアルゴリズムと、擬似的な眼底造影画像をさらに処理することで、眼底疾患のため診断脈絡血管構造から定量評価するために有用な値（具体的には、脈絡膜の直径、脈絡膜網の厚さ）を取得するアルゴリズムと、を有するものである。
【背景技術】
【０００３】
従来、医療分野等で用いられる非侵襲の断層計測技術の１つとして、時間的に低コヒーレンス光をプローブ（探針）として用いる光断層画像化法「光コヒーレンストモグラフィー」（ＯＣＴ）がある（特許文献１参照）。ＯＣＴは、光を計測プローブとして用いるため、被計測物体の屈折率分布、分光情報、偏光情報（複屈折率分布）等が計測できるという利点がある。
【０００４】
基本的なＯＣＴ４３は、マイケルソン干渉計を基本としており、その原理を図８で説明する。光源４４から射出された光は、コリメートレンズ４５で平行化された後に、ビームスプリッター４６により参照光と物体光に分割される。物体光は、物体アーム内の対物レンズ４７によって被計測物体４８に集光され、そこで散乱・反射された後に再び対物レンズ４７、ビームスプリッター４６に戻る。
【０００５】
一方、参照光は参照アーム内の対物レンズ４９を通過した後に参照鏡５０によって反射され、再び対物レンズ４９を通してビームスプリッター４６に戻る。このようにビームスプリッター４６に戻った物体光と参照光は、物体光とともに集光レンズ５１に入射し光検出器５２（フォトダイオード等）に集光される。
【０００６】
ＯＣＴの光源４４は、時間的に低コヒーレンス光（異なった時刻に光源から出た光同士は極めて干渉しにくい光）の光源を利用する。時間的低コヒーレンス光を光源としたマイケルソン型の干渉計では、参照アームと物体アームの距離がほぼ等しいときにのみ干渉信号が現れる。
【０００７】
この結果、参照アームと物体アームの光路長差（τ）を変化させながら、光検出器５２で干渉信号の強度を計測すると、光路長差に対する干渉信号（インターフェログラム）が得られる。
【０００８】
そのインターフェログラムの形状が、被計測物体４８の奥行き方向の反射率分布を示しており、１次元の軸方向走査により被計測物体４８の奥行き方向の構造を得ることができる。このように、ＯＣＴ４３では、光路長走査により、被計測物体４８の奥行き方向の構造を計測できる。
【０００９】
このような軸方向の走査のほかに、横方向の機械的走査を加え、２次元の走査を行うことで被計測物体の２次元断面画像が得られる。この横方向の走査を行う走査装置としては、被計測物体を直接移動させる構成、物体は固定したままで対物レンズをシフトさせる構成、被計測物体も対物レンズも固定したままで、対物レンズの瞳面付近においたガルバノミラーの角度を回転させる構成等が用いられている。
【００１０】
以上の基本的なＯＣＴが発展したものとして、分光器を用いてスペクトル信号を得るスペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）と、光源の波長を走査してスペクトル干渉信号を得る波長走査型ＯＣＴ（Swept Source OCT、略して「ＳＳ−ＯＣＴ」という。）がある。ＳＤ−ＯＣＴには、フーリエドメインＯＣＴ（Fourier Domain OCT、略して「ＦＤ−ＯＣＴ」という。特許文献２参照）、及び偏光感受型ＯＣＴ（Polarization-Sensitive OCT、略して「ＰＳ−ＯＣＴ」という。特許文献３参照）がある。
【００１１】
ＦＤ−ＯＣＴは、被計測物体からの反射光の波長スペクトルを、スペクトロメーター（スペクトル分光器）で取得し、このスペクトル強度分布に対してフーリエ変換することで、実空間（ＯＣＴ信号空間）上での信号を取り出すことを特徴とするものであり、このＦＤ−ＯＣＴは、奥行き方向の走査を行う必要がなく、ｘ軸方向の走査を行うことで被計測物体の断面構造を計測可能である。
【００１２】
ＳＳ−ＯＣＴは、高速波長スキャニングレーザーにより光源の波長を変え、スペクトル信号と同期取得された光源走査信号を用いて干渉信号を最配列し、信号処理を加えることで３次元光断層画像を得るものである。なお、光源の波長を変える手段として、モノクロメーターを利用したものでも、ＳＳ−ＯＣＴとして利用可能である。
【００１３】
ドップラー光コヒーレンストモグラフィー（ドップラーＯＣＴ）による、網膜の血流の分布の計測が知られている。これは、上記ＦＤ−ＯＣＴ等を用いて、網膜の血流の分布が計測を行うことのできる手段であり、同様に、スペクトルドメインＯＣＴを使うことによって、横断面網膜血流画像形成が得られ、また次元の網膜の脈管構造も観察することができる。
【００１４】
本発明者らは、従来、ドップラーＯＣＴに着目し、生体内、特に眼底の血流を非侵襲で計測する手法の研究・開発を行って来た。従来、本発明者らはＳＤ−ＯＣＴを技術基盤上にドップラーＯＣＴを実装し、眼底血流検査に成功してきたが、近年ＯＣＴの技術基板はＳＤ−ＯＣＴから、次世代であるＳＳ−ＯＣＴに移行しつつある。
【００１５】
ところで、近年、眼科分野における病気として、緑内障や加齢黄斑変性症（ＡＭＤ）が多いことが社会的な問題となっており、その効果的な診断、治療のために、眼底の断層画像の計測データに基づき、定量的に評価する技術が重要となっている。
【００１６】
眼底において、脈絡膜は、網膜色素上皮（ＲＰＥ）と強膜との間にあり、多くの血管（脈絡膜血管）を有する層である。そして、脈絡膜は、網膜の外層に酸素及び栄養を供給する機能を有する。脈絡膜血管の形態の変化は、緑内障や加齢黄斑変性症（ＡＭＤ）等、循環異常に係る眼病の診断に役立つものである。そこで、脈絡膜血管の画像は、眼科学ではきわめて有用である。                                 
【００１７】
従来、脈絡膜血管の画像化手段として、通常、インドシアニングリーン蛍光眼底血管造影（ＩＣＧＡ）が使用されており、近赤外線及び近赤外線検出器を使用すると、蛍光造影より透過性を向上することができる。
【００１８】
ところで、上記緑内障や加齢黄斑変性症（ＡＭＤ）等の病気の診断に、血管を分離して識別できるようにすることは必要であり、そのために、２次元画像又は３次元画像の血管を分離するために多くのアルゴリズムが提案されている（非特許文献１、２参照）。
【００１９】
脈絡膜の血管の分離画像は、脈絡膜の血管の直感的な視覚像だけでなく、形態面における定量的な情報を得ることができ、これは、血流に関連する眼の病気の病理を理解するために有用なことである。
【００２０】
血管の識別は、通常は、追尾に基づく手段（非特許文献３参照）と、ウィンドウに基づく手段（非特許文献４、５参照）の２つの異なる手段で行われている。追尾に基づく手段は、血管を抽出する必要の部位をなぞる手段であり時間がかかり、他方、ウィンドウに基づく手段は、フィルターを利用して血管の境界等を検出する手段でありさらに画像の精度を上げる必要がある。
【００２１】
最近、病因検出及び領域成長手法を利用した脈絡膜の血管の自動抽出方法が開発され、３Ｄ画像を得て、脈絡膜の血管を抽出された（非特許文献６参照）
【００２２】
また、マルチスケールの３Ｄエッジフィルターを利用した自動血管分離手段が提案されている（非特許文献７参照）
【先行技術文献】
【特許文献】
【００２３】
【特許文献１】特開２００２−３１０８９７号公報
【特許文献２】特開平１１−３２５８４９号公報
【特許文献３】特開２００４−０２８９７０号公報
【非特許文献】
【００２４】
【非特許文献１】M. Martinez-Perez, A. Hughes, A.Stanton, S. Thom, A. Bharath, and K. Parker,“ Retinal blood vessel segmentation by means of scale-space analysis and region growing, " in“Medical Image Computing and Computer-Assisted Intervention lMICCAll9," ,vol. 1679 of Lecture Notes in Computer Science, C. Tay1or and A. Colchester, eds.(Springer Berlin Heidelberg, 1999), pp. 90-97.
【非特許文献２】J.Staa1, M. Abramoff, M. Niemeijer, M. Viergever, and B. van Ginneken,“Ridge-based vessel segmentation in color images of the retina,"IEEE Trans.Med.Imag. 501-509(2004).
【非特許文献３】W.Cai and A.C.S. Chung,“Multi-resolution vessel segmentation using normalized cuts in retinal images," Lecture Notes in Computer Science 4191, 928-936(2006).
【非特許文献４】Q. Li, J. You, L. Zhang, and P. Bhattacharya, “A multiscale approach to retinal vessel segmentation using gabor filters and scale multiplication, "Systems, Man and Cybernetics, 2006. SMC '06, IEEE International Conference on 4, 3521-35279(2006).
【非特許文献５】M.D. Saleh, C. Eswaran, and A. Mueen,“An automated blood vessel segmentation algorithm using histogram equalization and automatic threshold selection, " J. Digit. Imag. 24, 564-572(20011).
【非特許文献６】L.Zhang, K. Lee, M.NiemeiJer, R.F. Mu11ins, M.Sonka, and M.D. Abrmoff,“Automated segmentation of the choroid from clinical sd-oct, "Invest. Ophthalmo1. Vis.Sci. 53. 7510-7519(2012).
【非特許文献７】V. Kajic, M.Esmaeelpour, C. Glittenberg, M.F. Kraus, J. Honegger, R. Othara, S. Binder, J.G. Fujimoto,
【非特許文献８】N．Otsu，"A threshold selection method from gray-1evel histograms, "Systems，Man and Cybemetics，IEEE Transaction on 9，62?66(1979）．
【発明の概要】
【発明が解決しようとする課題】
【００２５】
現在、眼科診療において、眼底血管可視化検査ＩＣＧＡが広く用いられている。この検査は眼底疾患の確定診断に不可欠なものであるが、染色濃度の正面画像しか提供することができず、深さ方向については良い分解度は得られないという問題がある。
【００２６】
要するに、深さ方向の異なる層で得られる画像信号がオーバーラップしてしまい、特に、毛細血管の密度が多い部位では、血管の輪郭が不鮮明となってしまう。また、ＩＣＧＡは、染料を投入する侵襲検査であるために、患者の身体に負担をかけ、同時に、副作用のリスクの伴うという問題がある。
【００２７】
−方、本発明のベースとなるＯＣＴは非侵襲の眼底構造可視化技術であるが、ＯＣＴで得られる脈絡膜血管の画像の画像を得るための画像処理において、次のような問題が生じる。
（１）脈絡膜より表面側にある網膜の血管が陰となって脈絡膜血管の画像中に投影されて脈絡膜血管のように誤認される偽血管の画像が生じる。
【００２８】
（２）脈絡膜血管の画像を得るためには、脈絡膜の層のみ取り込む必要があるが、脈絡膜より奥側に位置する強膜との境界がＯＣＴの断層画像でははっきりせず、境界に近い強膜の一部も脈絡膜血管のように誤認される偽血管の画像が生じる。
【００２９】
（３）ＯＣＴで得られた計測データの画像化処理において、画像を構成する画素（ピクセル）ごとの白黒の２値化処理の際に生じる、例えば、２値データ（１又は０）が同じ複数の画素が大きな密度で存在する領域が、脈絡膜血管ではないのに脈絡膜血管のように見えてしまうという偽血管の画像が生じる。
【００３０】
なお、上記非特許文献１〜７に示される方法は、確かに文献の記載からすると良い結果が得られるようであるが、数学的モデル及び分離プログラムは比較的複雑である。
【００３１】
本発明は、従来、眼科診断に必須であるにもかかわらず侵襲的にしか得られなかった眼底造影と近い画像情報を非侵襲に得ることが可能とすることを課題とするとともに、ＯＣＴで得られた計測データから、脈絡膜の層内に存在する脈絡膜血管のみの画像を得て、それにより、脈絡膜血管の太さや、血管網の厚さ等の形態を正確に得て、定量評価を可能とし、眼科診断の的確性を向上させる光干渉断層計装置及びその画像処理プログラムを実現する。
【課題を解決するための手段】
【００３２】
本発明は上記課題を解決するために、光干渉断層計と、該光干渉断層計で取得したＯＣＴ計測データに基づき三次元のＯＣＴ断層画像を得て、該三次元のＯＣＴ断層画像の画像処理を行うコンピュータと、を備えた脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置であって、前記コンピュータは、前記三次元のＯＣＴ断層画像から脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段、として機能し、前記脈絡膜の血管の画像データを取得する手段は、ＯＣＴ計測データから脈絡膜の断層画像データを抽出し、該脈絡膜の断層画像データから脈絡膜における深さ方向の等間隔ごとの位置でスライスした画像のデータを抽出して、該スライスした画像のデータから脈絡膜の血管の画像データを得る構成であることを特徴とする脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置を提供する。                        
【００３３】
前記脈絡膜の血管の画像データを取得する手段は、大きさの異なる複数のウィンドウのそれぞれについて、前記スライスした画像の画素の色濃度が、予め定めた所定の閾値以上であるか否かに応じて画素毎に２値化し、血管推定部分が抽出された２値画像を得て、該血管推定部分が抽出された２値画像において、当該ウィンドウ内における各画素に対する２値データの異なる画素の存在する割合が、予め定めた所定値以上の領域は擬血管として消去するとともに、当該ウィンドウの寸法に対する予め定めた所定の大きさの直径より小さな直径の前記血管推定部分を、ノイズであり非血管として消去して、脈絡膜の血管網の２値画像のデータを得る構成であることが好ましい。
【００３４】
前記脈絡膜の血管の画像データを取得する手段は、前記三次元のＯＣＴ断層画像から眼底におけるブルッフ膜の前側の領域の深さ方向の勾配を検出することで、眼底における視細胞との内節外節結合部の位置のデータを得て、該内節外節結合部及びブルッフ膜の位置のデータを使用して網膜色素上皮周辺の高反射構造体を抽出し、該高反射構造体における光強度を平均することによって、高反射構造体における網膜血管による陰を特定する画像データを得て、該画像データを、線状を強調するフィルターであるフランジフィルターにかけて得られた画像データを２値化し、網膜血管の陰の２値画像のデータを得る構成であることが好ましい。
【００３５】
前記脈絡膜の血管の画像データを取得する手段は、前記脈絡膜の断層画像データに基づき、血管の部分に相当する全ての画素について、光強度の大小により、中小の血管と大きな血管に分類した分類データを作成し、該分類データに基づき、大きさの異なる複数のウィンドウのそれぞれについて得た前記脈絡膜の血管網の２値画像のデータを選択的に組み合わせるとともに、前記網膜血管の陰の２値画像のデータに基づき、脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータを形成する構成であることが好ましい。
【００３６】
前記脈絡膜の血管の形状の定量評価を行うためのデータを得る手段は、前記脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータに基づき、脈絡膜の血管の直径を見積もり、該見積もりで得られたデータを基に脈絡膜の血管の直径分布データを作成するとともに、脈絡膜血管網の厚さを計測し、該計測で得られたデータを基に脈絡膜血管網の厚さマップを作成する構成であることが好ましい。
【００３７】
本発明は上記課題を解決するために、光干渉断層計と、該光干渉断層計で取得したＯＣＴ計測データに基づき三次元のＯＣＴ断層画像を得て、該三次元のＯＣＴ断層画像の画像処理を行うコンピュータと、を備えた脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラムであって、前記コンピュータを、前記三次元のＯＣＴ断層画像から脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段、として機能させるものであり、前記脈絡膜の血管の画像データを取得する手段は、ＯＣＴ計測データから脈絡膜の断層画像データを抽出し、該脈絡膜の断層画像データから脈絡膜における深さ方向の等間隔ごとの位置でスライスした画像のデータを抽出して、該スライスした画像のデータから脈絡膜の血管の画像データを得ることを特徴とする脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置におけるコンピュータに搭載される画像処理プログラムを提供する。
【００３８】
前記画像処理プログラムによって、前記脈絡膜の血管の画像データを取得する手段は、大きさの異なる複数のウィンドウのそれぞれについて、前記スライスした画像の画素の色濃度が、予め定めた所定の閾値以上であるか否かに応じて画素毎に２値化し血管推定部分が抽出された２値画像を得て、該血管推定部分が抽出された２値画像において、当該ウィンドウ内における各画素に対する２値データの異なる画素の存在する割合が、予め定めた所定値以上の領域は擬血管として消去するとともに、当該ウィンドウの寸法に対する予め定めた所定の大きさの直径より小さな直径の前記血管推定部分を、ノイズであり非血管として消去して、脈絡膜の血管網の２値画像のデータを得ることが好ましい。
【００３９】
前記画像処理プログラムによって、前記脈絡膜の血管の画像データを取得する手段は、前記三次元のＯＣＴ断層画像から眼底におけるブルッフ膜の前側の領域の深さ方向の勾配を検出することで、眼底における視細胞との内節外節結合部の位置のデータを得て、該内節外節結合部及びブルッフ膜の位置のデータを使用して網膜色素上皮周辺の高反射構造体を抽出し、該高反射構造体における光強度を平均することによって、高反射構造体における網膜血管による陰を特定する画像データを得て、該画像データを、線状を強調するフィルターであるフランジフィルターにかけて得られた画像データを２値化し、網膜血管の陰の２値画像のデータを得ることが好ましい。
【００４０】
前記画像処理プログラムによって、前記脈絡膜の血管の画像データを取得する手段は、前記脈絡膜の断層画像データに基づき、血管の部分に相当する全ての画素について、光強度の大小により、中小の血管と大きな血管に分類した分類データを作成し、該分類データに基づき、大きさの異なる複数のウィンドウのそれぞれについて得た前記脈絡膜の血管網の２値画像のデータを選択的に組み合わせるとともに、前記網膜血管の陰の２値画像のデータに基づき、脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータを形成することが好ましい。
【００４１】
前記画像処理プログラムによって、前記脈絡膜の血管の形状の定量評価を行うためのデータを得る手段は、前記脈絡膜の血管のみを分離した脈絡膜網の血管の２値画像のデータに基づき、脈絡膜の血管の直径を見積もり、該見積もりで得られたデータを基に脈絡膜の血管の直径分布データを作成するとともに、脈絡膜血管網の厚さを計測し、該計測で得られたデータを基に脈絡膜血管網の厚さマップを作成することが好ましい。
【発明の効果】
【００４２】
本発明に係る眼底血管の可視化機能を有する光干渉断層計及びそのプログラムによれば、次のような効果が生じる。
（１）従来、眼科診断に必須であるにもかかわらず侵襲的にしか得られなかった眼底造影と近しい画像情報を非侵襲に得ることが可能とする。
【００４３】
（２）ドップラー機能を有しないＯＣＴ装置で得られた画像を解析することで、ドップラー機能を有するＯＣＴ（ドップラーＯＣＴ）を用いて得られる擬似眼底造影像と同様、もしくは、さらに高コントラストな画像をドップラ一拡張されていない通常のＯＣＴ画像から得ることができる。
【００４４】
（３）ドップラー機能を有しない通常のＯＣＴ装置で得られた画像の画像処理で得られたデータを利用することにより、眼底血管の血管系などの定量評価を可能とする。なお、本発明を、ドップラーＯＣＴ画像に適用すれば、従来、さらに高コントラストな擬似眼底造影像を得ることができる。
【図面の簡単な説明】
【００４５】
【図１】本発明に係る光干渉断層計装置の実施例の全体構成を示す図である。
【図２】上記実施例の光干渉断層計装置の画像処理を行うコンピュータを示す図である。
【図３】上記実施例の光干渉断層計装置に搭載された画像処理プログラムによるフローを示す図である。
【図４】上記実施例の光干渉断層計で取得される断層像を示す写真である。
【図５】（ａ）は上記実施例の光干渉断層計装置で取得される断層像をスライスした画像を示す写真であり、（ｂ）は（ａ）を２値化した白黒画像である。
【図６】（ａ）は上記実施例の光干渉断層計装置で取得される網膜のスライスした画像であり、（ｂ）は（ａ）を２値化した白黒画像であり、（ｃ）は（ｂ）の網膜血管を強調できるように画像処理した図である。
【図７】（ａ）〜（ｄ）は、それぞれ網膜色素上皮から２５μｍ、１００μｍ、１７５μｍ、２５０μｍだけ、深さ方向の内部にある正面像（en-face像）の血管径の分布を示す図である。
【図８】従来の基本的なＯＣＴを説明する図である。
【発明を実施するための最良の形態】
【００４６】
本発明に係る脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置及びその画像処理プログラムを実施するための実施するための形態を実施例に基づいて図面を参照して、以下に説明する。
【実施例】
【００４７】
本発明に係る光干渉断層計装置１は、図１に示すように、光干渉断層計２と、該光干渉断層計で取得したＯＣＴ計測データを画像処理するコンピュータ１７とを備えている。
【００４８】
光干渉断層計（ＯＣＴ）２は、三次元のＯＣＴ断層画像を得るためのＯＣＴ計測データが取得できれば、どのような種類のものでもよいが、この実施例では、光干渉断層計として、波長走査型光干渉断層計（ＳＳ−ＯＣＴ）を使用する。このＳＳ−ＯＣＴ２について、その概要を、まず説明する。
【００４９】
図１において、波長走査型光源３から出射された出力光を、ファイバー４を通してファイバーカップラー５に送る。この出力光を、ファイバーカップラー５において、ファイバー６を通して被計測物体７への照射する物体光と、ファイバー８を通して固定参照鏡９に照射する参照光に分割する。
【００５０】
物体光は、ファイバー６、レンズ１０、角度が可変な走査鏡１１及びレンズ１２を介して、被計測物体７に照射、反射され、同じルートでファイバーカップラー５に戻る。参照光は、ファイバー８、レンズ１３及びレンズ１４を介して固定参照鏡９に照射、反射されて同じルートでファイバーカップラー５に戻る。
【００５１】
そして、これらの物体光と参照光はファイバーカップラー５で重ねられ、ファイバー１５を通して光検知器１６（ＰＤ（フォトダイオード）等のポイントセンサが使用される。）に送られ、スペクトル干渉信号として検出され、スペクトル干渉信号がＯＣＴ計測データとして、コンピュータ１７に取り込まれる。
【００５２】
コンピュータ１７において、光検知器１６における検知出力であるＯＣＴ計測データに基づいて、被計測物体７の奥行き方向（Ａ方向）と走査鏡の走査方向であるＢ方向及びＣ方向の３次元のＯＣＴ断層画像が形成される。コンピュータ１７には、ディスプレー１８が接続されている。
【００５３】
波長走査型光源３は、時間的に波長を変化させて走査する光源であり、即ち波長が時間依存性を有する光源である。これにより、Ａ−スキャン（被計測物体の奥行き方向への走査。このスキャンの結果Ａラインが得られる）をするために、固定参照鏡９を走査（光軸方向に移動して行う光軸方向の走査）することなく、被計測物体７の奥行き方向の反射率分布を得て奥行き方向の構造を取得することができる。
【００５４】
従って、Ｂ−スキャン（光軸方向と直交する１次方向の走査）及びＣ−スキャン（光軸方向と直交する２次方向の走査）をするだけで、３次元のＯＣＴ断層画像を得るためのＯＣＴ計測データを取得することができる。
【００５５】
そして、光検知器１６における検知出力であるＯＣＴ計測データに基づいて、コンピュータ１７によって３次元のＯＣＴ断層画像が形成され、コンピュータ１７に接続されたディスプレー１８で表示することができる。
【００５６】
本発明に係る画像処理プログラムは、光干渉断層計装置を構成するコンピュータに搭載され、そのコンピュータを、光干渉断層計で取得されるＯＣＴ計測データを画像処理して、脈絡膜血管の太さや、血管網の厚さ等の形態を正確に取得する手段として機能させる画像処理プログラムに関する。      
【００５７】
（画像処理手段、画像処理プログラム）                                            
光干渉断層計で得られたＯＣＴ計測データの処理を、コンピュータ１７を用いて画像処理する画像処理手段の構成、及びコンピュータ１７に搭載される本発明に係る画像処理プログラムについて、以下説明する。              
【００５８】
上記構成の光干渉断層計で得られたＯＣＴ計測データは、画像処理装置として使用するコンピュータ１７に入力される。このコンピュータ１７は通常のコンピュータであり、図２に示すように、入力部２１、出力部２２、ＣＰＵ２３、記憶装置２４及びデータバス２５を備えている。コンピュータ１７の記憶装置２４は、本発明に係る画像処理プログラムが記憶されて搭載されている。
【００５９】
コンピュータ１７に入力された光干渉断層計で計測されたＯＣＴ計測データに基づき  本発明に係る画像処理プログラムは、図３のフローに示すように、ＯＣＴ断層画像を形成し、このＯＣＴ断層画像を処理して、脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段、及び脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段（脈絡膜血管の太さや、血管網の厚さ等の形態を正確に取得する手段として機能させる手段）として、コンピュータ１７を機能させる画像処理プログラムである。
【００６０】
脈絡膜の血管の画像のみを選択的に分離して脈絡膜の血管の画像データを取得する手段は、図３のフローにおける（１）〜（９）の脈絡膜の血管分離プロセスに示す各機能を行うものであり、具体的には、以下の（１）〜（９）で説明する手段を含む。
【００６１】
また、脈絡膜の血管の画像データに基づき脈絡膜の血管の形状の定量評価を行うためのデータを得る手段は、図３のフローにおける（１０）〜（１３）の脈絡膜の血管定量評価用データ作成プロセスに示す各機能を行うものであり、具体的には、以下の（１０）〜（１３）で説明する手段を含む。
【００６２】
以下、ＯＣＴ断層画像を画像処理するために、画像処理プログラムによって、コンピュータが機能する手段を、図３のフローに従って、順次説明する。
【００６３】
（１）脈絡膜の断層画像についての画像データを抽出する手段（図３の（１））  
ＯＣＴ計測データにより、被計測物体である眼底の３次元の断層画像が得られるが、画像処理プログラムは、コンピュータを、この３次元の断層画像のうち、大まかに脈絡膜の領域に相当する断層画像についての画像データを抽出するように機能させる。
【００６４】
このように脈絡膜領域を抽出するために、ブルッフ膜（脈絡膜と網膜色素上皮を隔てる構造）を検出して自動的に切り分ける。この切り出しは、脈絡膜及びブルッフ膜のＡ−ラインについて、網膜色素上皮から深さ方向に手前から奥に向かって信号の二乗強度を計算し、その勾配が最小値となる画素をもってブルッフ膜と認定しブルッフ膜は自動的に除去する。
【００６５】
（２）脈絡膜の３次元の断層画像を等間隔ごとの深さ位置でスライスした複数の画像のデータを抽出する手段（図３の（２））
画像処理プログラムは、コンピュータを、網膜色素上皮から、脈絡膜の深さ方向に等間隔の複数の深さ位置において、脈絡膜の３次元の断層画像をスライス（より詳しくは、深さ方向に直交して厚みを薄く切り取ること）し、各深さ位置における眼の正面から見た２次元の画像（正面像（en-face像））を抽出するように機能させる。
【００６６】
なお、断層画像をスライスした脈絡膜の領域は、例えば、網膜色素上皮と脈絡膜の間に位置するブルッフ膜から例えば３７０μｍの深さまでの領域について選択し、この領域には、脈絡膜の奥にある強膜も一部含まれる。
【００６７】
即ち、ブルッフ膜のすぐ下の領域について脈絡膜を十分に含む領域（例えば３７０μｍ）を選択し、解析領域記として平坦化（薄い平坦面としてスライス）する。このとき、黄斑の周りの典型的な最大脈絡膜の厚さを考慮に入れて、ボリュームを選択するため、ほとんどの場合、脈絡膜のすべてと強膜のいくらかを含む領域となる。
【００６８】
また、スペックルを減らすため適当な平均フィルター（２０×２０画素：水平方向に５８．６μｍ、深さ方向に５８．６μｍなど）をかけ、画質の向上を図る。
【００６９】
図４（ａ）〜（ｄ）は、それぞれ脈絡膜の深さ方向に異なる位置でスライス（水平面に２０×２０画素を含み、５８．６μｍ四方の面積であり、厚さは５８．６μｍでスライス）して抽出した画像を示し、網膜色素上皮からの深さが、（ａ）は２５μｍ、（ｂ）は１００μｍ、（ｃ）は１７５μｍ、（ｄ）は２５０μｍである。
【００７０】
次に、画像処理プログラムは、上記（２）でスライスした複数の画像のデータについて、正面視で４７μｍウィンドウ（４７μｍ四方の窓）、９４μｍウィンドウ（９４μｍ四方の窓）、１８８μｍウィンドウ（１８８μｍ四方の窓）、３７５μｍウィンドウ（３７５μｍ四方の窓）のそれぞれについて、以下の（３）、（４）、（５）の手段を順次行うように、コンピュータを機能させる。
【００７１】
（３）血管推定部分を抽出し２値画像を形成する手段（図３の（３））
画像処理プログラムは、コンピュータを、上記（２）でスライスした複数の画像のデータについて、血管推定部分を抽出するように機能させる。ここで、血管推定部分とは、「脈絡膜の血管であろうと推定できる程度の部分」であって、実際は、血管以外の部分、陰等、非血管部分を含んでいる。
【００７２】
血管推定部分を抽出し２値画像を形成する手段は、血管と非血管の画像データの色濃度の違いに着目し、上記ウィンドウ内に含まれる複数の画素毎に、色濃度が予め定めた所定の閾値以上の場合は「１」とし、該閾値未満の場合は「０」として、画素毎に２値化する。
【００７３】
そして、「１」とされた画素は血管推定部分とし黒を割り当て、「０」とされた画素は非血管部分として白を割り当てて、白黒の２値画像を形成する。このような動作を、１枚の画像データの全域について、上記ウィンドウをそれぞれずらすようにして行い、１枚の画像データの全域についての白黒の２値画像を形成する。
【００７４】
図５（ａ）はスライスした画像の一例を示す図であるが、血管推定部分を抽出し２値画像を形成する手段で、画素毎に２値化することで、図５（ｂ）に示すような白黒の２値画像が得られる。
【００７５】
以上の手段をより詳細に説明すると次のとおりである。正面像（en-face像）を部分に分割しそのウィンドウ（領域）について次の式（数１）で決まる閾値（アダプティブ閾値という）で２値化する。
【００７６】
【数１】
2015000131.tif 000003
【００７７】
ここで、ｘ０とｘ１は横方向の座標、ｙ０とｙ１は縦方向の座標である。和Σは（ｘ０、ｙ０） を中心とするウィンドウ（領域）内（例えば、ｗｘμｍ×ｗｙμｍ等）の画素についておこなう。Ｉ（ｘｉ、ｙｉ）は、OCT像強度の対数をとったものであり、ｇ（ｘｉ、ｙｉ）はエッジ強度とよばれ、次の式(数２)で定義される。
【００７８】
【数２】
2015000131.tif 000004
【００７９】
このようなフィルター（アダプティブフィルターという）により血管領域の画素を「１」、非血管領域の画素を「０」とした２値画像が得られる。このような本発明の手段は、高速にかつ簡単な手順でグレースケール画像を２値化する手段である。
【００８０】
（４）ビジネスフィルターにより疑血管を消去する手段（図３の（４））
上記（３）の血管推定部分を抽出し２値画像を形成する手段により、図５（ｂ）に示すような白黒の２値画像が得られるが、この２値画像に、血管でなく単なる背景部分（バックグラウンド）であっても血管と識別される非血管部分（これを「疑血管」と言う）が含まれてしまう。
【００８１】
このような場合、ウィンドウ（領域）内がバックグラウンドのみで血管を含んでいない場合、適切な２値化ができずバックグラウンドも血管と分類してしまう。その場合、連続した血管ではなく、細かなエッジの集合のようになる。本発明に係る画像処理プログラムは、このような疑血管をビジネスフィルターにより消去する手段として、コンピュータを機能させる。
【００８２】
ビジネスフィルターにより血管と擬血管を識別する原理は、次のとおりである。ウィンドウ内の「或る１つの画素Ａ」（その２値データは例えば「０」）に対して、「２値データの異なる画素Ｂ」（その２値データは例えば「１」）が該ウィンドウ内に存在する場合、上記「或る画素Ａ」はエッジピクセルとして定義される。
【００８３】
このようなエッジピクセルが２０％以上存在する領域にある血管推定部分は、擬血管として消去する。エッジピクセルの比率をビジネス値という。具体的には、図５（ｂ）に示すような白黒の２値画像の全領域について、ウィンドウをずらしてエッジピクセルの比率であるビジネス値を取得する。
【００８４】
真の血管の部分の画素の２値データは「１」とし、２値画像において黒が割り当てられ、非血管の部分の画素の２値データは「０」とし、２値画像において白が割り当てられている場合に、非血管の部分の領域は、ビジネス値が高くなる。
【００８５】
このビジネス値が２０％以上存在する領域は、擬血管として、消去する。擬血管が消去された後の２値画像（ここでは図示しない）では、例えば、図５（ｂ）の中程の高さの右側の領域の偽血管の部分は消去される。                  
【００８６】
さらに分かり易く説明すると、ビジネスフィルターは、２値画像のある特定の画素に着目し、その周囲８画素を用いてビジーであるかどうか定義する。周囲がすべて自分と同じ値であればビジーではなく、逆に、１画素でも自分と異なる値であればビジーとする。
【００８７】
上記アダプティブフィルターで用いたウィンドウ（領域）と同じウィンドウ（領域）において、そのビジネス値（ビジー度）は、ビジー画素の数（エッジ画素という）の数と、本ウィンドウ（領域）の全画素数の比の値で定義される。一般にランダムな画像であればビジネス値は上がり、ウィンドウ（領域）がはっきり分かれているような画像ではビジネス値は小さくなる。
【００８８】
例えば、ウィンドウ（領域）の大きさがｗｘ×ｗｙ＝３７５μｍのとき、のとき、ビジネス値が０．２ より大きなウィンドウ（領域）は血管を含まないウィンドウ（領域）であるとして非血管領域として「０」をセットする。
【００８９】
（５）非血管画素の消去する手段（図３の（５））
ウィンドウの大きさが血管の直径に対してかなり大きい場合は、血管におけるコントラストの変化部分が、血管と非血管を識別の支障となる。
【００９０】
即ち、ウィンドウ（領域）ｗの大きさが血管径よりずっと大きい場合十分なコントラストが得られず、血管構造の検出がうまくいかない（後記するが経験的にはウィンドウの寸法の５分の１より小さい場合など）。逆にいえば、検出された小さいサイズ（ウィンドウの寸法の５分の１より小さい）の血管はアーチファクト（偽像、ノイズ）といえ、これを消す必要がある。
【００９１】
そのために、本発明の画像処理プログラムは、上記（４）で得られた２値画像から、さらに非血管部分を構成する画素を消去する手段として、コンピュータを機能させる。以下、その機能の詳細を説明する。
【００９２】
まず、上記（４）で得られた２値画像を、逆転する。即ち、血管を表示する画素は「１」で黒であったが、それを「０」で白とし、血管以外の部分を表示する画素は「０」で白であったが、それを「１」で黒とする。
【００９３】
そして、２つの画素より小さな直径の粒部分は消去する。このような動作の目的は、大きな血管の小さな孔の部分を埋めることにある。このようなことをしないと、次に説明する形態的な分析（morphlogical analysys）において、上記粒部分は、２つの隣接する薄い血管として認識されてしまうから、その前処理として行うのである。
【００９４】
上記前処理を行ってから、形態的な分析（morphlogical analysys）を行う。本発明者らは、本発明の研究開発の過程で、上記（４）に示すビジネスフィルターにより疑血管を消去する手段では、ウィンドウの寸法より小さな直径（ウィンドウの寸法に対する予め定めた所定の大きさの直径より小さな直径、具体的には、ウィンドウの寸法の１／５より小さな直径）の血管は識別できないという知見を得た。        
【００９５】
換言すると、ウィンドウの寸法の１／５より小さな直径の血管推定部分は、ノイズであって非血管であるから消去する。このような知見に基づいて、形態的な分析（morphlogical analysys）の手始めとして、ウィンドウの寸法より１／５より小さな直径の血管推定部分は、ノイズであって非血管であるから消去する。
【００９６】
次に、形態的な分析（morphlogical analysys）として、ヘイウッドの円形ファクタ（Heywood Circularity Factor：「有る形状の外周の長さ」／「円周の長さ」）を用いて、血管とは異なる形状である、円形又は円形に近い形状の部分を構成する画素を消去する。このヘイウッドの円形ファクタ（Heywood Circularity Factor）自体は、周知であって、「有る形状の外周の長さ」／「円周の長さ」）で定義される。
【００９７】
具体的には、ヘイウッドの円形ファクタ＜１．５の部分を構成する画素を消去する。このファクタが１．５未満であると、円形又は円形に近い形状で長細い血管ではないと認識され、１．５以上あると、長細いもの、即ち、血管と認識される。このようにして、非血管部分を構成する画素を消去する。
【００９８】
以上整理すると、まず、２値画像を反転する（血管領域が「０」、それ以外が「１」）。２画素以下の「１」の領域を除去する（「０」とする）。これは大血管中の小さな穴（ノイズ）を除去するためである。
【００９９】
次いで除去したい構造の円（たとえば直径ｗ／５）とのコンボリューション（たたみ込み演算）をおこない、「１」の値を持つ領域を広げる（ダイレーション）。次いで領域の縮小（エロージョン）をおこなうことにより、小さい構造は除去されるが、大きな構造は維持される。
【０１００】
（６）網膜色素上皮周辺の高反射構造体を分離する手段（図３の（６））
脈絡膜のＯＣＴ断層画像には、その手前にある網膜の血管の陰が投影されており、これが脈絡膜の血管として誤認されてしまう。この誤認は、脈絡膜の厚さや脈絡膜の血管の太さを正確に定量評価する際に、阻害要因となる。この陰を取り除く必要がある。
【０１０１】
その解決の手始めとして、本発明に係る画像処理プログラムは、コンピュータを、網膜色素上皮周辺の高反射構造体を分離する手段として機能させる。この手段は、図３のフロー図に示すように、上記（１）の脈絡膜を抽出する動作と並行して、下記のような動作を行う。
【０１０２】
即ち、網膜色素上皮と視細胞（光受容体）との内節外節結合部（ＩＳ／ＯＳライン）は、高い反射率を持つ高反射構造体であり、この内節外節結合部に網膜血管が位置している。そのため、網膜血管の投影像はこの高反射構造体である内節外節結合部に最もはっきり投影されて見える。
【０１０３】
そのため、ＯＣＴ断層画像の光強度の正の最大勾配位置から高反射構造体に投影された網膜血管の投影像が抽出できる。具体的には、ブルッフ膜の前側の領域の深さ方向の勾配（変化度）を計測し光強度の正の最大勾配位置検出することで、網膜色素上皮と視細胞（光受容体）との内節外節結合部の位置のデータを得る。
【０１０４】
そして、この内節外節結合部及びブルッフ膜の位置のデータを使用することで、網膜色素上皮周辺の高反射構造体を分離して光強度を抽出することができる。分離された高反射構造体における光強度を平均することによって、図６（ａ）に示すように、高反射構造体における画像を得る。このようなデータは、次に説明する網膜血管の陰の信号を特定する
【０１０５】
（７）網膜血管による陰を分離する手段（図３の（７））
上記（６）で説明したとおり、網膜血管は、脈絡膜の画像データに陰として入り込み、スライスした画像のデータ内で脈絡膜の血管として誤認されてしまう。この誤認は、脈絡膜の厚さや脈絡膜の血管の太さを正確に定量評価する際に、阻害要因となる。
【０１０６】
そのため、上記（６）で得られた高反射構造体における画像を利用して、網膜血管による陰を脈絡膜の血管とは分離する必要がある。本発明に係る画像処理プログラムは、コンピュータを、網膜血管による陰を分離する手段として機能させ、次のような動作を行わせる。               
【０１０７】
高反射構造体における画像データを、線状を強調するフィルターであるフランジフィルターにかけて、適当な閾値を持つ非線形２値化処理（フランジフィルターなど）をおこない、血管部分を「１」、それ以外を「０」として、図６（ｂ）に示すような画像データを得る。
【０１０８】
その後、血管の連続性を確保するためにモルフォロジカルクロージング処理（「１」の値を持つ領域の縮小（エロージョン）の後で領域を広げる（ダイレーション）処理）をおこない、これにより、実際より若干太めに抽出された網膜血管画像の影は脈絡膜血管画像から差し引かれる。これにより、網膜血管像は完全に除去され、図６（ｃ）に示すような２値画像のデータを得ることができる。
【０１０９】
（８）光強度に基づき画素を分類する手段（図３の（８））
前記したとおり、大きさの異なる４つのウィンドウについて、上記（３）〜（５）のプロセスを行うが、後記する理由によって、４つのウィンドウについて得られたデータについて、適切に組み合わせて、最終的に、脈絡膜の血管の厚さ、直径等の定量的評価を行う。
【０１１０】
この４つのウィンドウについて得られたデータの適切な組み合わせのために、事前に大きな血管に相当する画素と、中小の血管に相当する画素を、次のように分類しておく。本発明に係る画像処理プログラムは、コンピュータを、ＯＣＴ断層画像における各画素の光強度に基づき画素を分類する手段として機能させる。
【０１１１】
結合プロセスでは、上記（１）で得られた脈絡膜の断面像部分における全ての画素を、まず、大きな血管を構成する画素と、中小の血管を構成する画素に分ける。計測光は、脈絡膜を通過し、血液で散乱するために、厚い脈絡膜の血管のＯＣＴで得られるスペクトル干渉信号の強度、即ちＯＣＴ断層画像の光強度は、網膜色素上皮近くに存在する薄い中小の血管よりも小さい。
【０１１２】
このようなＯＣＴ断層画像の光強度の違いによって、主に薄い中小の血管に対応する明るい画素と、主に厚い血管に相当する暗い画素に分類する。実際は、この分類は、周知のオオツ閾値（非特許文献８参照）によって行われ、全ての画素について、同時に、ＯＣＴ断層画像のその光強度がオオツ閾値より、大きい場合は中小の血管に分類し、小さい場合は大きな血管に分類し、大まかな分類であるが分類データを生成する。
【０１１３】
（９）異なるウィンドウに基づく画像処理データを結合する手段（図３の（９））
観測している全視野（たとえば６ｍｍ四方）を複数の領域に分割する場合、いくつかのスケールに分割し、それぞれから血管構造を抽出し、合成する手段が有効である。
【０１１４】
例えば、６ｍｍ四方を複数の、たとえば、１／１２８、１／６４、１／３２、１／１６の４つのスケールに分割すると、その領域の大きさは、それぞれ、ｗ（＝ｗｘ＝ｗｙ）＝４７μｍ、９４μｍ、１８８μｍ、３７５μｍ となる。
【０１１５】
血管径は数ミクロンから数百ミクロンに及んでいるため、アダプティブ閾値によって、それぞれの分割はそれに最適な直径の血管を抽出することができる。よって、上記複数の（４つの）スケールの結果を組み合わせて血管を抽出する手段が有効となる。
【０１１６】
この点についてさらに分かり易く説明する。上記（５）に示す手段によって、ウィンドウの寸法より１／５より小さな直径の血管推定部分は、ノイズであって非血管であるから消去した。しかしながら、それは、真の小さな血管が分離されることではない。小さなウィンドウは、適切な分離に適切ではない。
【０１１７】
その理由は次のとおりである。血管推定部分を抽出し２値画像を形成する手段は、ウィンドウ内に含まれる複数の画素毎に、色濃度が所定の閾値以上の場合は「１」とし、該閾値未満の場合は「０」として、画素毎に２値化した。
【０１１８】
もし、血管と背景の間の境界における真の縁部の画素の数が、当該ウィンドウ内できわめて少ない場合は、上記（５）に示す手段は適切な閾値を提供することができない。それ故、小さなウィンドウは、大きな直径を有する血管を識別できない。
【０１１９】
脈絡膜の血管の直径は、数μｍから数百μｍの範囲まであるので、直径において広いレインジで血管を識別することが必要である。そのために、本発明では、前記したとおり、複数の大きさのウィンドウを使用して画像処理を行い、マルチスケールで広い直径範囲について、血管の識別を可能とするのである。
【０１２０】
本実施例では、前記のとおり、ウィンドウの縦横の大きさｗを、４７μｍ、９４μｍ、１８８μｍ、３７５μｍとしたが、これは、６mmの画像の１／１２８、１／６４、１
／３２、１／１６に相当する大きさである。
【０１２１】
本発明に係る画像処理プログラムは、コンピュータを、４つの異なる大きさウィンドウに基づいて得られた画像処理データを結合する手段（どの大きさのウィンドウに基づく画像処理データを使用するか選択する手段）として機能させる。
【０１２２】
この画像処理データを結合する手段は、次のようにしてどの大きさのウィンドウに基づく画像処理データＶ（ｘｉ,ｙｉ,ｚｉ,ｗ）を使用するか選択する。ここで、ｗの大きさのウィンドウに基づく２値画像のデータ（血管画素は「１」、非血管画素は「０」）である。
【０１２３】
（ア）Ｉ（ｘｉ,ｙｉ,ｚｉ）≧ｋ＊の場合は、Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝４７）∪Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝９４）の画像処理データが使用できる。
【０１２４】
ここで、Ｉ（ｘｉ,ｙｉ,ｚｉ）は、ＯＣＴ断層画像のｘｉ,ｙｉ,ｚｉの画素の光強度（対数値）であり、ｋ＊はオオツ閾値である。Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝４７）は、ｗ＝４７μｍのウィンドウに基づく２値画像のデータである。Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝９４）は、ｗ＝９４μｍのウィンドウに基づく２値画像のデータである。∪は論理和演算を意味する。
【０１２５】
要するに、上記（８）においてオオツ閾値ｋ＊以上で分類された画素に対応して、４７μｍ又は９４μｍの大きさのウィンドウに基づく２値画像処理データを採用する。これは、小さな血管については、小さなウィンドウに基づく２値画像処理データを採用するという意味である。
【０１２６】
（イ）Ｉ（ｘｉ,ｙｉ,ｚｉ）≧ｋ＊ではない場合は、Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝９４）∪Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝１８８）∪Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝３７５）の２値画像のデータが使用できる。
【０１２７】
ここで、Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝１８８）はｗ＝１８８μｍのウィンドウに基づく２値画像のデータであり、Ｖ（ｘｉ,ｙｉ,ｚｉ,ｗ＝３７５）はｗ＝３７５μｍのウィンドウに基づく２値画像のデータである。∪は論理和演算を意味する。
【０１２８】
要するに、上記（８）においてオオツ閾値ｋ＊以上ではないと分類された画素に対応して、９４μｍ、１８８μｍ又は３７５μｍの大きさのウィンドウに基づく２値画像のデータを採用する。これは、大きな血管については、大きなウィンドウに基づく２値画像のデータを採用するという意味である。
【０１２９】
以上のようにして、本発明に係る画像処理プログラムは、コンピュータを、断層画像の画素の光強度に対応して、大きさの異なるウィンドウに基づく２値画像のデータを組み合わせる作業を行う機能をさせることとなる。
【０１３０】
なお、上記（６）における網膜色素上皮周辺の高反射構造体を分離する手段と、上記（７）における網膜血管による陰を分離する手段による２値画像のデータは、この（９）のプロセスにおいて、異なる大きさのウィンドウに基づく２値画像のデータに反映される。
【０１３１】
（１０）脈絡膜の血管の直径を見積もる手段（図３の（１０））
本発明に係る画像処理プログラムは、上記（９）において最終的に得られた２値画像のデータに基づき、白黒の２値画像の各部の血管の直径を見積もる手段として機能させる。
【０１３２】
具体的には、血管の長手方向に直交する方向の長さをもって直径とする。円形パターンによるモルフォロジカルオープニング（「１」の値を持つ領域の領域を広げる（ダイレーション）の後で縮小（エロージョン）処理）によってパターンより大きな直径を持つ血管が抽出される。
【０１３３】
従って、抽出したい最小の直径をｄ１とし、最大の直径をｄｉｍａｘとすると、ｄｉ＝ｉ×ｄ１（ｉ＝１……ｉｍａｘ） の径の円でモルフォロジカルオープニングすることで血管径を分類することができる。
【０１３４】
（１１）脈絡膜の血管の直径の分布データを作成する手段（図３の（１１））
本発明に係る画像処理プログラムは、上記（１０）において得られた２値画像の各部の血管の直径の見積もりデータに基づき、脈絡膜の血管の直径の分布データを作成する手段として機能させる。
【０１３５】
上記（１０）における血管径を分類する処理を、異なる深さの正面像（en-face像）に施すことにより、図７に示すような血管径の分布を求めることができる。
【０１３６】
図７（ａ）〜（ｄ）は、それぞれ網膜色素上皮（ＲＰＥ）から２５μｍ、１００μｍ、１７５μｍ、２５０μｍだけ、深さ方向の内部にある正面像（en-face像）の血管径の分布を示す図であり、明るさは、血管径を示しており、明るさと血管径の大きさの指標は図７の下に示す。但し、脈絡膜血管のない領域と網膜血管の領域は、明るさとは関係なく、黒く見える部分である。
【０１３７】
ところで、上記のとおり血管径の分布を求める手段は、円形パターンが小さくなると多くの時間を要する。円形パターンの直径を最大の検出したい血管径に固定し（例えば、ｄ１＝ｄｆ（＝ｄｉｍａｘ））、その代わり、血管像の正面像（en-face像）をＮ／ｉ× Ｎ／ｉにリサイズ（画素数を減らし画像サイズをｉ 分の１にすること）しながらモルフォロジカルオープニングをおこなうことで、高速に血管径の分類をおこなうことができる。
【０１３８】
（１２）脈絡膜の血管網の厚さを計測する手段（図３の（１２））
本発明に係る画像処理プログラムは、上記（１０）において得られた２値画像の各部の血管の直径の見積もりデータに基づき、脈絡膜の血管網の厚さを計測する手段として機能させる。
【０１３９】
この脈絡膜の血管網の厚さを計測する手段は、具体的には、次のとおりである。脈絡膜血管網の厚さは、網膜色素上皮（ＲＰＥ）と脈絡膜血管の後ろ側（奥側）の包絡面をなめらかに結んだエンベロプ（いわば、血管網の領域に風呂敷をかけたような面）の領域に風呂敷をの間の厚さである。
【０１４０】
このエンベロープを能動的変形表面モデル（active deformable surface model）を用いて抽出する。まず、網膜色素上皮（ＲＰＥ）を平坦化し裏返し、上部（眼底の奥側）から変形可能なメッシュをかける。個々での説明では、ｚ座標（眼底の深さ方向の座標）について、手前側（眼底の奥側）を正（ｚ）とする。
【０１４１】
メッシュにおける格子点を制御点とし、この制御点に相互作用（力）を与えると表面の変形がおこり、なめらかで、最適なエンベロープとなる。この作業をすべての制御点に繰り返しおこなう。相互作用（力）Ｆｊは下記の式で定義される。
【０１４２】
Ｆｊ（τ） ＝αＦｊ（τ）＋βＲｊ（τ）＋G
ここで、τは繰り返し数、ｊは制御点の番号であり、Ｇは負の定数（重力のようなもの）でメッシュ（制御点）全体に下向きの力をかける。Ｒｊ（τ）は変形表面の局所的な硬直性で、次の式（数３）で表される。
【０１４３】
【数３】
2015000131.tif 000005
【０１４４】
ここでＳ（ｘｉ、ｙｉ）は制御点間を、例えば、２次元のバイキュービック補間（求める位置の周辺の４×４（１６点）の高さ情報を用いて面を３次関数で表し補間する方法）等適当な補間方法によって得られたなめらかな表面であり、（ｕｊ、νｊ）はｊ番目の制御点の座標である。Ｐｊ（τ）は局所的な場所における平均血管径から定義される斥力である。
【０１４５】
局所的な場所は、例えば、制御点を中心としてメッシュの大きさにとることができる。αとβは正の定数で上記３つの力のバランスをとる。制御点の位置（ｚ座標）は繰り返し操作により更新される。つまり、制御点の位置は、次の式（数４）に示される。
【０１４６】
【数４】
2015000131.tif 000006
【０１４７】
ここで、ｚｊ（τ）はτ番目の繰り返しにおけるｊ番目の制御点の位置であり、Ｆθは予め決められている正の閾値である。繰り返しはすべての制御点が安定（収束）するまで繰り返される。収束したら制御点について、例えば、２次元のバイキュービック補間等、適当な補間をおこない血管網のエンベロープを求める。
【０１４８】
（１３）脈絡膜の血管網の厚さマップを作成する手段（図３の（１３））
本発明に係る画像処理プログラムは、上記（１２）において得られた脈絡膜の血管網の厚さを計測データに基づき、脈絡膜の血管網の厚さマップを作成する手段として機能させる。
【０１４９】
以上、脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置及びその画像処理プログラムを実施するための形態を実施例に基づいて説明したが、本発明はこのような実施例に限定されることなく、特許請求の範囲記載の技術的事項の範囲内で、いろいろな実施例があることは言うまでもない。
【産業上の利用可能性】
【０１５０】
本発明に係る脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置及びその画像処理プログラムを使用すれば、ドップラ一拡張されていない通常の光干渉断層計装置で取得されたＯＣＴ画像から、脈絡膜の血管網を他の層及び非血管から分離して抽出し、定量評価を可能とするので、眼科診断装置への利用が最適である。
【０１５１】
例えば、擬似血管造影検査、脈絡膜血管の定量評価による緑内障、糖尿病網膜症、加齢黄斑変性症等の超早期診断に、きわめて有用であり、さらに、光干渉断層計装置ドップラ一拡張されたＯＣＴ画像へも適用可能である。
【符号の説明】
【０１５２】
１  光干渉断層計装置
２  波長走査型ＯＣＴ         
３  波長走査型光源           
４、６、８、１５  ファイバー 
５  ファイバーカップラー     
７  被計測物体               
９  固定参照鏡               
１１  走査鏡                 
１０、１２、１３、１４  レンズ     
１６  光検知器               
１７  コンピュータ           
１８  ディスプレー 
２１  入力部
２２  出力部
２３  ＣＰＵ
２４  記憶装置
２５  データバス
４３  ＯＣＴ
４４  光源 
４５  コリメートレンズ 
４６  ビームスプリッター 
４７  物体アーム内の対物レンズ 
４８  被計測物体
４９  参照アーム内の対物レンズ 
５０  参照鏡
５１  集光レンズ 
５２  光検出器（フォトダイオード等） 
【図１】
2015000131.tif 000007
【図２】
2015000131.tif 000008
【図３】
2015000131.tif 000009
【図８】
2015000131.tif 000010
【図４】
2015000131.tif 000011
【図５】
2015000131.tif 000012
【図６】
2015000131.tif 000013
【図７】
2015000131.tif 000014
【手続補正書】
【提出日】平成26年2月6日(2014.2.6)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００２４
【補正方法】変更
【補正の内容】
【００２４】
【非特許文献１】M. Martinez-Perez, A. Hughes, A.Stanton, S. Thom, A. Bharath, and K. Parker,“ Retinal blood vessel segmentation by means of scale-space analysis and region growing, " in“Medical Image Computing and Computer-Assisted Intervention lMICCAll9," ,vol. 1679 of Lecture Notes in Computer Science, C. Tay1or and A. Colchester, eds.(Springer Berlin Heidelberg, 1999), pp. 90-97.
【非特許文献２】J.Staa1, M. Abramoff, M. Niemeijer, M. Viergever, and B. van Ginneken,“Ridge-based vessel segmentation in color images of the retina,"IEEE Trans.Med.Imag. 23, 501-509(2004).
【非特許文献３】W.Cai and A.C.S. Chung,“Multi-resolution vessel segmentation using normalized cuts in retinal images," Lecture Notes in Computer Science 4191, 928-936(2006).
【非特許文献４】Q. Li, J. You, L. Zhang, and P. Bhattacharya, “A multiscale approach to retinal vessel segmentation using gabor filters and scale multiplication, "Systems, Man and Cybernetics, 2006. SMC '06, IEEE International Conference on 4, 3521-3527(2006).
【非特許文献５】M.D. Saleh, C. Eswaran, and A. Mueen,“An automated blood vessel segmentation algorithm using histogram equalization and automatic threshold selection, " J. Digit. Imag. 24, 564-572(20011).
【非特許文献６】L.Zhang, K. Lee, M.Niemeijer, R.F. Mu11ins, M.Sonka, and M.D. Abrmoff,“Automated segmentation of the choroid from clinical sd-oct, "Invest. Ophthalmo1. Vis.Sci. 53, 7510-7519(2012).
【非特許文献７】V. Kajic, M.Esmaeelpour, C. Glittenberg, M.F. Kraus, J. Honegger, R. Othara, S. Binder, J.G. Fujimoto, and W. Drexler, "Automated three-dimensional choroidal vessel segmentation of 3d 1060 nm oct retinal data," Biomed. Opt. Express 4, 134-150(2013).
【非特許文献８】N．Otsu，"A threshold selection method from gray-1evel histograms, "Systems，Man and Cybernetics，IEEE Transactions on 9，62-66(1979）．