JP2016221111A - 画像処理装置及び画像処理方法 - Google Patents

Abstract

【課題】網膜等の深さ方向の情報の取得条件によらず、良好な二次元画像を得ることを可能する画像処理装置を提供する。
【解決手段】走査された測定光の被測定物体からの戻り光と測定光に対応する参照光とを干渉させた干渉光に基づいて被測定物体の深さ方向に並ぶ画素値列を取得する取得手段により得られた画像を処理する画像処理装置において、複数の前記画素値列のそれぞれから所定の選択基準に従って選択された画素値に基づいて二次元画像を生成する生成手段と、複数の画素値列の各々において、深さ方向において選択される画素値の選択範囲を設定する設定手段と、設定された選択範囲に応じて所定の選択基準を変更する基準変更手段と、を配する。
【選択図】図３

Description

本発明は、光干渉で得られた断層画像の情報から二次元の画像を得る画像処理装置及び画像処理方法に関するものである。

多波長光波干渉を利用した光コヒーレンストモグラフィ（ＯＣＴ：Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）による撮影装置（以降ＯＣＴ装置と称する。）が、知られている。ＯＣＴ装置は、例えば内視鏡での内蔵情報や、眼科装置での網膜の情報を得るために用いられ、人体に対する適用分野を広げつつある。特に眼に適用したＯＣＴ装置は、眼科用機器として網膜の専門外来では必要不可欠な装置になりつつある。

ＯＣＴ装置では、低コヒーレント光である測定光をサンプルに照射し、干渉系を用いてそのサンプルからの後方散乱光を測定することで、画像の取得を可能にしている。そして、眼に適用した場合には、測定光を被検眼上で走査することにより被検眼の断層画像を高解像度に撮像することが可能である。

また、複数の断層画像から疑似的に眼底を正面から見た二次元画像（以下、「二次元画像」と呼ぶ。）を生成する技術が知られている。二次元画像を生成する方法として、１回のＡスキャンによって得られた深さ方向の画素値列の中から画素値の大きさの順に基づき画素値を選択する。そして、それら選択された画素値を全てのＡスキャンに対して得ることで、断層画像のみで網膜の平面画像と類似した二次元画像（Ｐｒｏｊｅｃｔｉｏｎ画像）を生成している（特許文献１）。

また、特許文献１には、網膜の所定層を選択した場合に、当該所定層に関する二次元画像を生成する方法も開示されている。当該方法では、その層内における画素値列から画素値の大きさの順に基づき画素値を選択し、二次元画像を生成する。

特開２０１４−４５８６９号公報

ところで、網膜の深さ方向の全ての画素値列と、深さ方向の所定の層を包含する範囲内の画素値列とでは、画素値列に占める被測定対象物の割合が大きく異なる。そのため、深さ方向の全ての画素値列から一つの画素値を選択する場合と、所定範囲内の画素値列から一つの画素値を選択する場合とで、同じ基準に基づいて画素値を選択すると、良好な二次元画像が得られないことがある。

また、例えば病変を有する層に関しての二次元画像を得ようとした場合、該病変の影響を受けた組織に対応する画素値がその周辺部と大きく異なることが生じ得る。この場合も、画一的な基準を用いた場合に、該病変に対応する画素値を削除してしまう恐れも存在する。

本発明は上記課題に鑑み、網膜等の深さ方向の情報の取得条件によらず、良好な二次元画像を得ることを可能とする画像処理装置及び画像処理方法の提供を目的とする。

上記目的を達成するために、本発明に係る画像処理装置は、
走査された測定光の被測定物体からの戻り光と前記測定光に対応する参照光とを干渉させた干渉光に基づいて前記被測定物体の深さ方向に並ぶ画素値列を取得する取得手段と、
複数の前記画素値列のそれぞれから所定の選択基準に従って選択された画素値に基づいて二次元画像を生成する生成手段と、
前記複数の画素値列の各々において、前記深さ方向において前記選択される画素値の選択範囲を設定する設定手段と、
前記設定された選択範囲に応じて前記所定の選択基準を変更する基準変更手段と、を有することを特徴とする。

本発明によれば、網膜等の深さ方向の情報の取得条件によらず、良好な二次元画像を得ることが可能となる。

本発明の第１の実施形態に係る画像処理装置を有する撮影システムの構成図である。 図１に例示した撮影システムを構成する眼科装置の側面図である。 本発明の第１の実施形態における画像処理装置の構成を示すブロック図である。 図１に示す撮影システムに供せられる撮影装置の光学系の構成図である。 第１の実施形態の第一のモードにおける二次元画像生成の流れを示すフローチャートである。 Ａスキャン画像の画素値列を示す説明図である。 画素値の並べ替えの一例を示す説明図である。 第１の実施形態の第二のモードにおける二次元画像生成の流れを示すフローチャートである。 深度範囲の設定例を示す説明図である。 図１に示す撮影システムにおける断層画像及び二次元画像の表示例を示す図である。 図１に示す撮影システムにおけるＰｒｏｊｅｃｔｉｏｎ画像及びＥｎＦａｃｅ画像の表示例を示す図である。

以下、本発明の実施形態について、図面を参照して説明する。なお、以下の実施の形態は特許請求の範囲に関わる本発明を限定するものではなく、また、本実施の形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。

［第１の実施形態］
図１は、本発明の第１の実施形態に係る画像処理装置１００及びそれと接続される撮影装置１を有する撮影システム１０００の構成例を示す図である。本形態に係る画像処理装置１００は、中央演算処理装置（ＣＰＵ）１０、主メモリ１１、磁気ディスク１２、及び表示メモリ１３で構成される。また、撮影システム１０００は、モニタ９２８、マウス９２９−１、及びキーボード９２９−２を有する。なお、モニタ９２８がタッチパネルを有することとしてもよい。

ＣＰＵ１０は、主として、画像処理装置１００の各構成要素の動作を制御する。主メモリ１１は、ＣＰＵ１０が実行する制御プログラムを格納したり、ＣＰＵ１０によるプログラム実行時の作業領域を提供する。磁気ディスク１２は、オペレーティングシステム（ＯＳ）、周辺機器のデバイスドライバ、後述する再構成処理等を行うためのプログラムを含む各種アプリケーションソフト等を格納する。表示メモリ１３は、モニタ９２８のための表示用データを一時記憶する。モニタ９２８は、例えば、ＣＲＴモニタや液晶モニタ等であり、表示メモリ１３からのデータに基づいて画像を表示する。マウス９２９−１及びキーボード９２９−２はユーザによるポインティング入力及び文字等の入力をそれぞれ行う。上記の各構成要素は共通バス１７により互いに通信可能に接続されている。

画像処理装置１００は、ローカルエリアネットワーク（ＬＡＮ）を介して撮影装置１と接続されており、撮影装置１から画像データを取得できる。なお、本発明の形態はこれに限定されず、これらの機器との接続は、例えば、ＵＳＢやＩＥＥＥ１３９４等の他のインターフェイスを介して行ってもよい。また、これらのデータを管理するデータサーバ等の外部装置３から、ＬＡＮ等を介して必要なデータを読み込む構成であってもよい。また、画像処理装置１００に記憶装置、例えば、ＦＤＤ、ＣＤ−ＲＷドライブ、ＭＯドライブ、ＺＩＰドライブ等を接続し、それらのドライブから必要なデータを読み込むようにしてもよい。

図２は、以上の構成により眼科装置の態様にて構築された該撮影装置１の側面図である。同図に示す撮影装置１は、外部固視標３２４、顎台３２３、ベース部９５１、ステージ部９５０、及び光学ヘッド９００を有する。光学ヘッド９００は、被検眼の前眼画像、及び眼底の表面画像及び断層画像を撮像するための測定光学系である。ステージ部９５０は、光学ヘッド９００を図中ｘｙｚ方向に不図示のモータを用いて移動可能とした移動部を構成する。ベース部９５１は、後述の分光器を内蔵する。

パーソナルコンピュータ（以下「パソコン」と称する。）９２５はステージ部の制御部を兼ね、画像処理装置１００を有する。顎台３２３は、被検者の顎と額とを固定することで、被検者の眼（被検眼）の固定を促す。外部固視標３２４は、被検者の眼を固視させるのに使用する。また、画像処理装置１００を光学ヘッド９００、あるいはステージ部９５０に組み込むこともできる。この場合には撮影装置１と画像処理装置１００が撮影装置として一体的に構成される。

図３は、画像処理装置１００の機能構成を示すブロック図であり、該画像処理装置１００は再構成部１１００、生成部１２００、及び層認識部１３００を有する。再構成部１１００は再構成手段として機能し、測定光の被測定物体からの戻り光と参照光とを干渉させて得た干渉光に基づいて被測定物体の所定範囲の断層画像を得る。より詳細には、センサからの出力値に対して波数変換と高速フーリエ変換（ＦＦＴ）処理とを行い、被検眼眼底上の一点における深さ方向の断層画像（Ａスキャン画像）として再構成する。生成部１２００は、再構成部１１００で得られた断層画像の深さ方向のそれぞれの画素値列から所定の画素をそれぞれの画素値列毎に選択して二次元画像を生成する生成手段として機能する。層認識部１３００は、再構成部１１００で得られた断層画像において被測定対象物（網膜）の層構造を抽出し、それぞれの層境界の形状を特定する。

図４は撮影装置１の測定光学系および分光器の構成を説明する図である。光学ヘッド９００部の内部、及びベース部９５１の内部に配置される光学要素について、以下に説明する。

なお、本実施形態では、被測定物体としての被検眼１０７を例示して以下の説明を行う。光学ヘッド９００では、該被検眼１０７に対向して対物レンズ１３５−１が設置される。該対物レンズ１３５−１の光軸上には、第１ダイクロイックミラー１３２−１及び第２ダイクロイックミラー１３２−２が配置される。被検眼１０７に至る光路は、第１ダイクロイックミラー１３２−１によって前眼観察用の光路３５３へ、波長帯域に応じて分岐される。また分岐後の光路は更に、第２ダイクロイックミラー１３２−２によって固視灯用の光路３５２と、眼底観察のためのＯＣＴ光学系の光路３５１とに、波長帯域ごとに分岐される。

固視灯用光路３５２上には、被検眼１０７側より順に、合焦レンズ１３５−３、レンズ１３５−４、穴あきミラー３０３、第３ダイクロイックミラー１３２−３が配置される。第３ダイクロイックミラー１３２−３の反射方向にはＣＣＤ１７２が、透過方向には固視灯１９１が配置される。ここで、合焦レンズ１３５−３は、固視標１９１および眼底観察用のＣＣＤ１７２の合焦調整のため不図示のモータによって光路３５２に沿った方向に駆動される。

穴あきミラー３０３は、レンズ１３５−４と第３ダイクロイックミラー１３２−３の間に配置され、光路３５２から光路３５４への光路の分岐を為す。

光路３５４は、被検眼１０７の眼底を照明する照明光学系を形成している。被検眼１０７の位置合わせに用いられる眼底観察用照明光源であるＬＥＤ光源３１６、被検眼１０７の眼底の撮像に用いるストロボ管３１４が設置されている。該照明光学系は、最奥のＬＥＤ光源３１６より順に配置される、コンデンサレンズ３１５、ストロボ管３１４、コンデンサレンズ３１３、ミラー３１７、リングスリット３１２、レンズ３１１及びレンズ３０９を有する。

ＬＥＤ光源３１６とストロボ管３１４とからの照明光はリングスリット３１２によってリング状の光束となり、孔あきミラー３０３によって反射され、被検眼１０７の眼底１２７を照明する。ＬＥＤ光源３１６は、７８０ｎｍ付近を中心波長とする光源である。

前述したように、光路３５２上の穴あきミラー３０３以降には、第３ダイクロイックミラー１３２−３が配置される。該第３ダイクロイックミラー１３２−３による眼底観察用のＣＣＤ１７２及び固視灯１９１への光路分岐は、他のダイクロイックミラーと同じく波長帯域ごとにおこなわれる。

ＣＣＤ１７２は眼底観察用照明光であるＬＥＤ光源３１６から発せられる光の中心波長、具体的には７８０ｎｍ付近に感度を持つものであり、ＣＣＤ制御部１０２に接続されている。一方固視標１９１は可視光を発生して被検者の固視を促すものであり、固視標制御部１０３に接続されている。

ＣＣＤ制御部１０２及び固視標制御部１０３は、共に演算部１０４に接続されており、演算部１０４を通じて、データはパソコン９２５へ入出力される。

第１ダイクロイックミラー１３２−１に分岐された光路３５３には、分岐部より順に、レンズ１３５−２及び前眼観察用の赤外線ＣＣＤ１７１が配置される。このＣＣＤ１７１は不図示の前眼観察用照明光の波長、具体的には９７０ｎｍ付近に感度を持つものである。また、光路３５３には、不図示のイメージスプリットプリズムが配置されており、被検眼１０７に対する光学ヘッド９００部のｚ方向の距離（被検眼１０７に対する接近及び離間の方向）を、前眼観察画像中のスプリット像として検出することができる。

光路３５１は前述の通りＯＣＴ光学系を成しており、実際に被検眼１０７の網膜の断層画像を撮像するための構成が配置される。より具体的には、断層画像を形成するための干渉信号を得るための構成が配置される。光路３５１には、第２ダイクロイックミラー１３２−２から順に、ＸＹスキャナ１３４、合焦レンズ１３５−５、及びレンズ１３５−６が配置され、ファイバー１３１−２端部より当該光路３５１に測定光が供給される。

ＸＹスキャナ１３４は、測定光を眼底上で走査するために用いられる。なお、ＸＹスキャナ１３４は一枚のミラーとして図示してあるが、ＸＹ２軸方向の走査を行うものである。合焦レンズ１３５−５は、後述する光ファイバー１３１−２の端部から出射する光源１０１からの光を眼底１２７上に合焦調整をするために不図示のモータによって光路３５１に沿って駆動される。また、この合焦調整によって眼底１２７からの光は、同時に光ファイバー１３１−２先端にスポット状に結像されて入射されることとなる。

次に、光源１０１からの光路と参照光学系、分光器の構成について説明する。

光源１０１は、光ファイバー１３１−１を介して光カプラー１３１に接続される。光カプラー１３１は光ファイバー１３１−１及び１３１−２に加え、光ファイバー１３１−３及び１３１−４に接続されている。これら光ファイバーは光カプラー１３１に接続されて一体化されているシングルモードに光ファイバーである。光源１０１から射出された光は光カプラー１３１において測定光と参照光とに分割され、測定光は光ファイバー１３１−２を介して光路３５１へ、参照光は光ファイバー１３１−３を介して後述する参照光路へ導かれる。参照光路には、該光路奥より順に、ミラー１３２−４、分散補償用ガラス１１５、及びレンズ１３５−７が配置される。分散補償用ガラス１１５は、測定光と参照光の分散を合わせるために参照光路中に挿入される。

該参照光路を経た参照光は光ファイバー１３１−３の端部に入射され、同様に光カプラー１３１に戻った測定光と合波されて光ファイバー１３１−４を経て分光器１８０に導かれる。本実施形態におけるこれらの構成は、マイケルソン干渉系を構成している。

測定光は前述のＯＣＴ光学系光路を通じ、観察対象である被検眼１０７の眼底に照射され、網膜による反射や散乱により同じ光路を通じて前述したように光カプラー１３１に到達する。一方、参照光は光ファイバー１３１−３、レンズ１３５−７、分散補償ガラス１１５を介してミラー１３２−４に到達し反射される。そして同じ光路を戻り光カプラー１３１に到達する。

光カプラー１３１によって、測定光と参照光とは合波され干渉光となる。ここで、測定光の光路長と参照光の光路長がほぼ同一となったときに干渉を生じる。ミラー１３２−４は不図示のモータおよび駆動機構によって光軸方向に調整可能に保持され、被検眼１０７によって変わる測定光の光路長に参照光の光路長を合わせることが可能である。干渉光は光ファイバー１３１−４を介して分光器１８０に導かれる。

また、測定光側の偏光調整部１３９−１が、光ファイバー１３１−２中に設けられる。参照光側の偏光調整部１３９−２も、光ファイバー１３１−３中に設けられる。これらの偏光調整部は光ファイバーをループ状に引き回した部分を幾つか持ち、このループ状の部分をファイバーの長手方向を中心として回動させることでファイバーに捩じりを加える。この操作によって、測定光と参照光の偏光状態を各々調整して合わせることが可能となる。本装置ではあらかじめ測定光と参照光の偏光状態が調整されて固定されている。

次に分光器１８０について述べる。分光器１８０はレンズ１３５−８、レンズ１３５−９、回折格子１８１、及びラインセンサ１８２から構成される。光ファイバー１３１−４から出射された干渉光はレンズ１３５−８を介して略平行光となった後、回折格子１８１で分光され、レンズ１３５−３によってラインセンサ１８２に結像される。ラインセンサ１８２の出力は、パーソナルコンピュータ９２５へと入力される。

次に、光源１０１の周辺について説明する。光源１０１には、代表的な低コヒーレント光源であるＳＬＤ（Ｓｕｐｅｒ Ｌｕｍｉｎｅｓｃｅｎｔ Ｄｉｏｄｅ）を用いる。該光源１０１から発せられる光の中心波長は８５５ｎｍ、波長バンド幅は約１００ｎｍである。ここで、バンド幅は、得られる断層画像の光軸方向の分解能に影響するため、重要なパラメータである。また、光源の種類としては、ここではＳＬＤを選択したが、低コヒーレント光が出射できればよく、ＡＳＥ（Ａｍｐｌｉｆｉｅｄ Ｓｐｏｎｔａｎｅｏｕｓ Ｅｍｉｓｓｉｏｎ）等も用いることができる。中心波長は眼を測定することを鑑みると、近赤外光が適する。また、中心波長は得られる断層画像の横方向の分解能に影響するため、なるべく短波長であることが望ましい。双方の理由から中心波長を８５５ｎｍとした。

なお、本実施形態では干渉計としてマイケルソン干渉計を用いたが、マッハツェンダー干渉計を用いてもよい。測定光と参照光との光量差に応じて光量差が大きい場合にはマッハツェンダー干渉計を、光量差が比較的小さい場合にはマイケルソン干渉計を用いることが望ましい。

次に本撮影装置を用いた被検眼の撮像方法について説明する。

まず、検者は本実施例に基づく撮影装置の前に患者を着座させ、被検眼眼底の表面画像の撮影を開始する。光源３１６を射出した光は、リングスリット３１２によってリング状の光束となり、孔あきミラー３０３によって被検眼１０７に向けて反射される。反射された光は、被検眼１０７の眼底１２７を照明する。該照明による眼底１２７からの反射光束は孔あきミラー３０３を通過し、ＣＣＤ１７２上において結像される。ＣＣＤ１７２上にて結像された眼底１２７の反射光は、ＣＣＤ制御部１０２にて眼底の表面画像として画像化され、画像処理装置１００へ送信される。

次に、撮影装置１はＸＹスキャナ１３４を制御することで、被検眼１０７の眼底における所望部位の断層画像を撮像する、より詳細には所望部位の各位置での深さ方向の輝度情報を取得する。

被検眼に向かう測定光は光ファイバー１３１−２を通過しファイバー端から射出され、ＸＹスキャナ１３４へ入射する。ＸＹスキャナ１３４により偏向された測定光は光学系１３５−１を経由して被検眼眼底１２７を照射する。そして被検眼で反射した反射光は逆の経路をたどって光カプラー１３１へと戻される。

一方、参照ミラーに向かう参照光は光ファイバー１３１−３を通過しファイバー端から射出され、コリメート光学系１３５−７及び分散補償光学系１１５を通して参照ミラー１３２−４に到達する。参照ミラー１３２−４で反射された参照光は逆の経路をたどって光カプラー１３１へと戻される。

光カプラー１３１に戻ってきた測定光と参照光は相互に干渉し、干渉光となって光ファイバー１３１−４へと入射し、光学系１３５−８により略平行化され回折格子１８１に入射する。回折格子１８１に入力された干渉光は結像レンズ１３５−９によってラインセンサ１８２に結像し、被検眼眼底上の一点における干渉信号を得ることができる。

ラインセンサ１８２の複数の素子で取得された干渉情報を有する画像信号としての出力値を画像処理装置１００に出力する。なお、図４では、ストロボ管３１４の発光で一度に眼底の表面画像を取得する形態を説明したが、ＳＬＤ光源で発光した光を走査するＳＬＯタイプで眼底の表面画像を得る構成としてもよい。

以上述べたＯＣＴ光学系の構成は、本実施形態において、走査された測定光の被測定物体からの戻り光と測定光に対応する参照光とを干渉させた干渉光に基づいて被測定物体の深さ方向に並ぶ画素値列を取得する取得手段に対応する。

次に画像処理装置１００の画像処理方法の流れを説明する。

被検眼眼底上の一点における断層情報を取得した後、撮影装置１は走査手段としてのＸＹスキャナ１３４の駆動により測定光の照射位置をＸ方向に移動させ、被検眼眼底上の別の一点の干渉光を発生させる。該別の一点の干渉光はラインセンサ１８２を経由し、再構成部１１００に入力される。再構成部１１００は、該入力されたデータに基づいて、被検眼眼底上の別の一点における深さ方向の断層画像（Ａスキャン画像）を形成する。該Ａスキャンの干渉信号を撮影したＸＹスキャナ１３４の位置とＡスキャン画像の眼底画像上の座標とは関連付けて記憶されている。

このＸＹスキャナ１３４による測定光照射位置のＸ方向への移動を連続して行うことにより、被検眼眼底の水平方向の一枚の断層画像（Ｂスキャン画像）を再構成部１１００は再構成する。

ＸＹスキャナ１３４により測定光照射位置をＹ方向に一定量移動した後、上述したＸ方向の走査を再び行うことで、被検眼眼底上の別のＹ方向位置における眼底の水平断層画像（Ｂスキャン画像）を再構成部１１００は再構成する。このＸＹスキャナ１３４による測定光照射位置のＹ方向移動を繰り返すことで、眼底１２７の所定範囲をカバーする複数枚の断層画像を形成することができる。本撮影装置１ではＹ方向に１２８回の一定量の微小駆動を行いながらＢスキャン画像の形成を繰り返すことで、１２８枚の断層画像を生成する。また、再構成部１１００は、１２８枚の断層画像から三次元の断層画像を再構成（形成）する。

次に再構成部１１００で生成された断層画像から網膜の二次元画像を生成部１２００で生成する。

本実施形態において、生成部１２００は動作モードとして第一のモードと第二のモードとを有する。第一のモードは、Ａスキャン画像のための輝度情報を得た全ての深度範囲に対応する画素から二次元画像を生成する。第二のモードは、輝度情報を得た深度範囲内の所定の深度範囲から二次元画像を生成する。第一のモードで生成される二次元画像（以下、Ｐｒｏｊｅｃｔｉｏｎ画像と呼ぶ）は強度画像とも呼ばれ、眼底の表面画像と類似したものである。一方、第二のモードで生成される二次元画像（以下、ＥｎＦａｃｅ画像と呼ぶ）は、網膜における任意の深さの情報により生成される平面状の画像であって、眼疾患による網膜層構造の変化を可視化するために用いられる。なお、ユーザはマウス９２９−１、キーボード９２９−２およびタッチパネルの少なくとも一つを用いて第一のモードおよび第二モードから一方のモードを選択することが可能である。例えば、ユーザは、モニタに表示された第一のモードを示すＧＵＩおよび第二のモードを示すＧＵＩのうち一方のＧＵＩをマウス９２９−１を用いてクリックすることでモードを選択することができる。モニタ９２８は、第一のモードおよび第二モードをプルダウンにより選択可能に表示することとしてもよい。

ここで、図５を用いて第一のモードによるＰｒｏｊｅｃｔｉｏｎ画像の生成処理について説明する。前述したように、再構成部１１００にて生成された断層画像は生成部１２００へと入力される（Ｓ２０００）。

Ａスキャン画像は被検眼眼底上の一点における深さ方向の断層画像であり、図６に示すように深さ方向の複数の輝度情報から構成されている。

図６の二次元の断層画像はＡスキャン画像の集合である。この二次元の断層画像はＢスキャン画像である場合と、三次元に再構成された断層画像の断面を示す場合がある。なお、図６からわかるように断層画像には網膜領域とそれ以外の部分とが含まれている。網膜領域とは例えばＩＬＭとＲＰＥとにより挟まれた領域である。

例えば、本実施形態で用いた撮影装置１では２０４８の画素を持つラインセンサ１８２を用いており、ＦＦＴ後のＡスキャン画像Ａｉは１１７６個の画素値から画素値列が構成されている。ここでＰ０は深さ方向の一番浅い部分の輝度情報としての画素値の値を色の濃さで示しており、Ｐ１１７５は深さ方向の一番深い部分の輝度情報としての画素値を示している。これらＡスキャン画像から平面画像を得る場合、本撮影装置１は、これら複数の輝度情報の中から一つの輝度情報を選択抽出する。この抽出された輝度情報に対応する画素値は、該Ａスキャン画像を得た被検眼眼底上の一点の代表的な強度信号とする。すなわちＡスキャンで得られる１１７６個の画素値から一つの画素値を選択する。

ここで、生成部１２００は、外部装置３から取得され且つ再構成された断層画像を処理してＰｒｏｊｅｃｔｉｏｎ画像を生成するように構成してもよい。この場合には、輝度情報等は、再構成部１１００を介さずに不図示のデータ取得部から直接入力を受ける態様とすることが好ましい。

生成部１２００は、図７に示すようにＡスキャン毎に対応する断層画像の輝度情報を輝度の大きい順に並べ替えを行う。すなわち１１７６個の画素値列毎に画素値の大小関係に基づいて画素値を順位付けし、画素値の並び替え（ソート処理）を行う（Ｓ２０１０）。なお、同図において、Ｒ０は最も明るい輝度情報を画素値として持つ画素であり、Ｒ１１７５は最も暗い輝度情報を画素値として持つ画素である。輝度は干渉の強さを示すため、画素値も干渉の強さに対応している。生成部１２００は、更に並び替え後の画素値列における所定順位の画素Ｒｘを選択する。ここで所定順位の画素とは、輝度情報の大きい順に並べ替えた後、先頭からｘ番目に位置している画素のことである。

第一のモードでＰｒｏｊｅｃｔｉｏｎ画像を生成する場合、Ａスキャン画像に含まれる全１１７６個の画素値から所定順位の画素値を選択する。網膜の断層画像を撮影した場合、Ａスキャン画像の中で被測定対象物（網膜）の占める割合は約２０％程度に留まる。つまり、Ａスキャン画像を構成する全１１７６個の画素のうち約２０％が明るい網膜部の画素であり、残りの約８０％は被測定対象物以外の暗い画素で構成される。そのため、網膜の平均的な輝度を取得する場合、ｘは総画素数の半分よりも高順位に位置している画素であることが望ましい。具体的には総画素数１１７６の画素値列で構成されるＡスキャン画像を用いる場合、上位１０％の位置に相当する先頭から１１８番目の画素を所定順位の画素Ｒｘとして選択する（Ｓ２０２０）。

生成部１２００は、全てのＡスキャン画像に対して上位１０％の位置に相当する画素Ｒｘを選択し、それらの画素値を用いてＰｒｏｊｅｃｔｉｏｎ画像を生成する（Ｓ２０３０）。

このＰｒｏｊｅｃｔｉｏｎ画像はＣＣＤ１７２にて得られる眼底の表面画像や、他の眼底カメラ、ＳＬＯで得られる眼底画像と類似した画像であり、疑似的に眼底表面を可視化することが可能である。また、複数の輝度情報の中から有効情報のみを選択的に取得するため、Ａスキャン画像に含まれるノイズ成分や干渉の強度が低い暗部の領域に左右されることなく好適な表面画像を得ることが可能である。

なお、ここでは、再構成部１１００で全データを再構成した後に生成部１２００でＰｒｏｊｅｃｔｉｏｎ画像を生成する例で説明した。しかし、Ａスキャン毎に再構成された断層画像を生成部１２００に順次送信する構成や、Ｂスキャン毎に再構成された断層画像を生成部１２００に順次送信する構成としてもよい。

次に図８を用いて第二のモードによるＥｎＦａｃｅ画像の生成処理について説明する。

層認識部１３００は二次元の断層画像から網膜の層構造（層境界）を抽出し、それぞれの層境界の形状を特定する。ここで特定される層境界の形状はＩＬＭ（網膜-硝子体境界）、ＮＦＬ／ＧＣＬ境界、ＧＣＬ／ＩＰＬ境界、ＩＰＬ／ＩＮＬ境界、ＩＳ／ＯＳライン、ＲＰＥ、ＢＭ等が挙げられる。一般的にＥｎＦａｃｅ画像においてはＲＰＥや推定ＲＰＥ（ＢＭ）の形状が用いられる事が多いが、角膜形状など画像解析によって得られた任意の形状を用いても良い。特定された層境界の形状は断層画像と共に生成部１２００へと入力される（Ｓ３０００）。ここで、生成部１２００は、図示しない取得部が外部装置３から取得した層境界の形状及び断層画像を用いてＥｎＦａｃｅ画像を生成するように構成してよい。この場合に、生成部１２００は、再構成部１１００及び層認識部１３００を介さずに取得部２０００から直接入力を受ける。

生成部１２００は、入力された層境界の形状に基づき、断層画像を構成する各Ａスキャン画像においてＥｎＦａｃｅ画像を生成するための深度範囲Ｚｉとして、画素値列の一部分を設定する。網膜全層のＥｎＦａｃｅ画像を得る場合、図９（ａ）に示すようにＩＬＭからＲＰＥまでの深度範囲Ｚを設定する。また、例えば、ＲＰＥ近傍の浮腫等を観察する場合、図９（ｂ）に示すように推定ＲＰＥ形状を有する２本のラインをＲＰＥを挟んで上下に等距離に配置することで深度範囲Ｚｉを設定してもよい。

また、深度範囲Ｚｉを設定する２本のラインはマウス９２９−１等で上下に移動可能に構成することもできる。この場合、推定ＲＰＥ形状を有するラインの１本をＲＰＥ上に設定し、もう１本を脈絡膜下に配置することで、脈絡膜のＥｎＦａｃｅ画像を得るための深度範囲Ｚｉを設定することもできる。なお、本実施形態では、深さ方向の画素値の選択範囲を設定する際に、表示手段において図９（ａ）等に表示される表示形態を利用して深度範囲を設定させている。しかし、深度範囲を設定する際の表示形態はこれに限定されず、カーソルによって個々の層境界をクリックして選択する等種々の表示形態を用いることが可能である。この表示形態の表示は、パソコン９２５において表示制御手段として機能するモジュールにより実行される。例えば、図９（ａ）に表示された断層画像における２つの層境界がクリック等を用いてユーザにより選択された場合、選択された２つの層境界に挟まれた領域が深度範囲Ｚとして設定される。クリック等を用いてユーザにより選択される層境界は、断層画像に重畳表示された層境界を示す表示であってもよいし、断層画像自体に含まれる層境界であってもよい。

さらに、層境界の形状を用いずに深度範囲Ｚを設定することもできる。例えば、硝子体のＥｎＦａｃｅ画像を得るために、断層画像の上端部からＩＬＭ境界までを深度範囲Ｚと設定することができる。また、断層画像の任意の深度位置に２本の直線を設定し、その２本の直線の間を深度範囲Ｚとすることもできる。

そして、生成部１２００は全てのＡスキャン画像Ａｉに対して深度範囲Ｚｉを設定する（Ｓ３０１０）。生成部１２００は、設定された深度範囲に基づきＥｎＦａｃｅ画像を生成する。即ち、生成部１２００における設定手段として機能するモジュールは、選択範囲として深さ方向のそれぞれの画素値列の各々の一部分を深度範囲Ｚｉとして設定する。

前述のように各Ａスキャン画像Ａｉは１１７６個の画素値から画素値列が構成されている。生成部１２００は、各Ａスキャン画像Ａｉに対して設定された深度範囲Ｚｉに基づき、深度範囲内に相当する画素値のみを取得し、深度範囲Ｚｉにおける画素値列を生成する。以下、例えば、深度範囲における画素値列に１９２個の画素値が含まれるものとして説明する。

生成部１２００は、Ａスキャン毎に設定された深度範囲の画素値列内で画素値の大小関係に基づいて画素値を順位付けし、画素値の並び替えを行う（Ｓ３０２０）。そして生成部１２００は、所定順位の画素Ｒｘを選択する。ここで所定順位の画素とは、輝度情報の大きい順に並べ替えた後、先頭からｘ番目に位置している画素のことである。

第二のモードでＥｎＦａｃｅ画像を生成する場合、設定された深度範囲に含まれる全１９２個の画素値から所定順位の画素値を選択する。Ａスキャン画像において深度範囲を設定する場合、深度範囲の中で着目領域（例えば、被測定対象物である網膜）の占める割合はほぼ１００％である。これはＥｎＦａｃｅ画像を作成する場合、着目領域のみが含まれるよう深度範囲を設定するためである。そのため、着目領域の平均的な輝度を取得する場合、ｘはＰｒｏｊｅｃｔｉｏｎ画像と同じ上位１０％の位置では不適切であり、総画素数の半分程度であることが望ましい。具体的には総画素数１９２の画素値列を用いる場合、上位５０％の位置に相当する先頭から８１番目の画素を所定順位の画素Ｒｘとして選択する（Ｓ３０３０）。

生成部１２００は、全てのＡスキャン画像に対して深度範囲Ｚｉ内の画素値列から上位５０％の位置に相当する画素Ｒｘを選択し、それらの画素値を用いてＥｎＦａｃｅ画像を生成する（Ｓ３０４０）。

なお、上記ではＥｎＦａｃｅ画像を生成する際に深度範囲内の画素値列から上位５０％の位置に相当する画素を選択したが、これは上位５０％の位置の画素に限るものではない。所定順位は画素値列の輝度分布（ヒストグラム）に応じて設定することもできるし、深度範囲内に占める着目領域の割合に応じて所定順位を可変とすることもできる。例えば、画素値列の輝度分布（ヒストグラム）に応じて所定順位を設定する際に、深度範囲内に含まれる網膜の占める割合に応じて、上位何％の位置の画素を選択するかを変更することもできる。このような構成では、画像に含まれる網膜の割合に応じて適切な画素が選択されるため、患者ごとの網膜厚の違いによらず良好なＥｎｆａｃｅ画像を生成することができる。例えば、層認識部１３００はＩＬＭとＲＰＥとを断層画像から抽出し、ＩＬＭとＲＰＥとに挟まれる領域を網膜の領域として認識する。すなわち、層認識部１３００は、断層画像における網膜領域の座標を取得する。一方、生成部１２００はユーザにより指定された深度範囲Ｚｉの座標を取得する。そして、生成部１２００は層認識部１３００から断層画像における網膜領域の座標を取得し、深度範囲Ｚｉの座標と比較することで、深度範囲Ｚｉ内において網膜の占める割合を算出することができる。すなわち、生成部１２００は画素値列に含まれる注目領域の占める割合を算出することができる。この割合に基づいて生成部１２００は自動的に上位何％の位置の画素を選択するかを変更することができる。 なお、生成部１２００は、深度範囲Ｚｉ内に網膜が含まれる割合が少ない程、より上位の位置の画素を選択するようにする。言い換えれば、生成部１２００は、深度範囲Ｚｉ内に網膜が含まれる割合が多い程、より下位の位置の画素を選択するようにする。すなわち、基準変更手段の一例である生成部１２００は、第二のモードにおいて、画素値列に含まれる注目領域の占める割合に応じて二次元画像の生成に用いる画素値の順位を変更する。
なお、生成部１２００はＡスキャン画像毎に網膜の占める割合を算出してＡスキャン画像毎に上位何％の位置の画素を選択するかを変更することとしてもよいし、複数Ａスキャン画像における網膜の占める割合の平均値に基づいて上位何％の位置の画素を選択するかを変更することとしてもよい。また、生成部１２００は、第一のモードにおいても、Ａスキャン画像に含まれる網膜の割合に応じて上位何％の位置の画素を選択するかを変更することとしてもよい。

このような本実施形態において所定順位として述べられる設定された選択範囲に応じた所定の選択基準の変更は、ここで述べたように網膜の割合に例示される患者個々の眼の特徴等に応じて行われることが好ましい。また、予め保存されている患者個人のデータより、任意に選択し条件に基づいて行われても良い。この場合、当該処理は同画像処理装置１００中において基準変更手段として機能するモジュールにより実行される。また、この基準変更は、例えば網膜の占める割合の値が閾値を超えるか否かで判断し実行されても良く、検者による入力、指示等に応じて実行されても良い。

このＥｎＦａｃｅ画像はＣＣＤ１７２にて得られる眼底の表面画像とは異なり、深度範囲の設定によって様々な様態の画像を得ることができる。特定の深度範囲内の情報のみで二次元画像を生成するため、深度範囲内における被測定対象物の構造とその構造の変化を顕著に可視化することができる。

次に、本実施形態に係る撮影装置１では、生成した眼底の表面画像、断層画像、Ｐｒｏｊｅｃｔｉｏｎ画像、ＥｎＦａｃｅ画像をモニタ９２８に表示する。図１０に示すように、眼底の表面画像Ｓと、断層画像Ｔｉ、Ｐｒｏｊｅｃｔｉｏｎ画像Ｐはモニタ９２８上に並べて表示される。また、Ｐｒｏｊｅｃｔｉｏｎ画像Ｐ及び表面画像Ｓ上には断層画像Ｔｉの取得位置Ｌｉが重ねて表示される。Ｐｒｏｊｅｃｔｉｏｎ画像Ｐは、不図示の切り替え手段によりＥｎＦａｃｅ画像Ｅに切り替えて表示することもできる。或いは、図１１に示すように、断層画像Ｔｉ、Ｐｒｏｊｅｃｔｉｏｎ画像Ｐ、及びＥｎＦａｃｅ画像Ｅをモニタ９２８上に並べて表示することもできる。

なお、本実施形態における画像処理装置１００では、１２８枚の断層画像を生成する。しかし、モニタ９２８上には選択された１枚の断面としての断層画像Ｔｉ（ｉ＝０〜１２８）、又は三次元で再構成された断層画像の断面画像Ｔｉ（この場合には任意の番号ｉが振られる）が表示される。検者は入力部９２９−１、９２９−２を操作して、表示する断層画像を切り替える切り替え表示可能とすることができる。ここで、断層画像Ｔｉが切り替わると、Ｐｒｏｊｅｃｔｉｏｎ画像Ｐ、ＥｎＦａｃｅ画像Ｅ及び表面画像Ｓ上に表示される断層画像Ｔｉの取得位置Ｌｉの表示位置も更新される。これにより、検者は表示された断層画像Ｔｉが被検眼眼底１２７上のどの位置の断層画像なのかを容易に知ることができる。

そして、取得位置ＬｉがＥｎＦａｃｅ画像Ｅ上におかれた場合、特定の深度範囲における被測定対象物の構造変化と断層画像Ｔｉを対比しながら観察することも可能である。

なお、本実施例では被検眼眼底の断層画像に基づいて網膜のＰｒｏｊｅｃｔｉｏｎ画像及びＥｎＦａｃｅ画像を生成しているが、被検眼前眼部の断層画像に基づいて前眼部のＰｒｏｊｅｃｔｉｏｎ画像及びＥｎＦａｃｅ画像を生成しても良い。

以上述べたＰｒｏｊｅｃｉｏｎ画像及びＥｎＦａｃｅ画像を生成する生成部１２００は、本実施形態において、複数の前記画素値列のそれぞれから所定の選択基準に従って選択された画素値に基づいて二次元画像を生成する。この場合、所定の選択基準は前述した上位何パーセントという基準が例示される。

また、この所定の選択基準は、画素値選択時に用いるＡスキャン画像中の画素値列として用いる画素の数に応じて変更される。この場合、複数の画素値列の各々において、深さ方向において幾つの或いはどのような条件の範囲に含まれる画素値を用いるかは、選択される画素値の選択範囲を設定する設定手段として機能する画像処理装置１００中のモジュールにより実行される。また、設定された選択範囲に応じた所定の選択基準の変更は、同画像処理装置１００中において基準変更手段として機能するモジュールにより実行される。

前述したように、本発明によれば、深さ方向の全ての情報に基づき二次元画像を生成する場合と、深さ方向の所定範囲内の情報に基づき二次元画像を生成する場合の双方で、良好な二次元画像を得ることが可能となる。

［その他の実施形態］
なお、本件は上記の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲内において、種々の変形、変更して実施することができる。例えば、上記の実施形態では、被測定物体が眼、特に眼底の場合について述べているが、眼以外の皮膚や臓器等の被測定物体に本発明を適用することも可能である。この場合、本発明は眼科装置以外の、例えば内視鏡等の医療機器としての態様を有する。従って、本発明は眼科装置に例示される検査装置のための画像処理装置として把握され、被検眼は被検査物体の一態様として把握されることが望ましい。

また、本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。

Claims (12)

1. 走査された測定光の被測定物体からの戻り光と前記測定光に対応する参照光とを干渉させた干渉光に基づいて前記被測定物体の深さ方向に並ぶ画素値列を取得する取得手段と、
   複数の前記画素値列のそれぞれから所定の選択基準に従って選択された画素値に基づいて二次元画像を生成する生成手段と、
   前記複数の画素値列の各々において、前記深さ方向において前記選択される画素値の選択範囲を設定する設定手段と、
   前記設定された選択範囲に応じて前記所定の選択基準を変更する基準変更手段と、を有することを特徴とする画像処理装置。
2. 前記生成手段は、前記選択範囲が前記深さ方向の画素値列の全てであって、前記所定の選択基準として前記画素値の大きさの順に基づき選択された画素値を用いて前記二次元画像を生成する第一のモードと、
   前記選択範囲として前記深さ方向のそれぞれの画素値列の各々の一部分が前記設定手段に設定され、前記所定の選択基準として前記選択範囲の内で画素値の大きさの順に基づき選択された画素値を用いて二次元画像を生成する第二のモードとを備えることを特徴とする請求項１に記載の画像処理装置。
3. 前記生成手段は画素値列毎に画素値の大小関係に基づいて画素値を順位付けし、
   前記順位づけされた画素値列における所定の順位の画素値に基づいて前記二次元画像を生成し、
   第一のモードと第二のモードで所定の順位が異なることを特徴とする請求項２に記載の画像処理装置。
4. 前記基準変更手段は、前記第二のモードにおいて、前記画素値列に含まれる注目領域の占める割合に応じて前記所定の順位を変更することを特徴とする、請求項３に記載の画像処理装置。
5. 前記第一のモードにおける所定の順位は前記第二のモードにおける所定の順位よりも高順位であることを特徴とする請求項３又は４に記載の画像処理装置。
6. 前記生成手段で生成された前記二次元画像を表示するための表示手段を更に備えることを特徴とする請求項１乃至５のいずれか一項に記載の画像処理装置。
7. 前記設定手段が前記深さ方向の画素値の選択範囲を設定する際に、前記表示手段において前記設定手段に対して選択範囲を設定させる表示形態を表示させる表示制御手段を有することを特徴とする請求項６に記載の画像処理装置。
8. 前記表示手段は、前記生成手段で生成された前記第一のモードによる二次元画像と前記第二のモードによる二次元画像とを切り替えて表示可能であることを特徴とする請求項６又は７に記載の画像処理装置。
9. 前記表示手段は、前記生成手段で生成された前記第一のモードによる二次元画像と前記第二のモードによる二次元画像とを並べて表示可能であることを特徴とする請求項６又は７に記載の画像処理装置。
10. 前記被測定物体は、眼であることを特徴とする請求項１乃至９の何れか１項に記載の画像処理装置。
11. 走査された測定光の前記被測定物体からの戻り光と前記測定光に対応する参照光とを干渉させた干渉光に基づいて被測定物体の深さ方向に並ぶ画素値列を取得する取得工程と、
    複数の前記画素値列のそれぞれから所定の選択基準に従って選択された画素値に基づいて二次元画像を生成する生成工程と、
    前記複数の画素値列の各々において、前記深さ方向において前記選択される画素値の選択範囲を設定する設定工程と、
    前記設定された選択範囲に応じて前記所定の選択基準を変更する基準変更工程と、を有することを特徴とする画像処理方法。
12. 請求項１１に記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。