JP6301621B2 - 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法 - Google Patents

2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法 Download PDF

Info

Publication number
JP6301621B2
JP6301621B2 JP2013202023A JP2013202023A JP6301621B2 JP 6301621 B2 JP6301621 B2 JP 6301621B2 JP 2013202023 A JP2013202023 A JP 2013202023A JP 2013202023 A JP2013202023 A JP 2013202023A JP 6301621 B2 JP6301621 B2 JP 6301621B2
Authority
JP
Japan
Prior art keywords
edge line
dimensional tomographic
tomographic image
image processing
anterior segment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2013202023A
Other languages
English (en)
Other versions
JP2015066084A (ja
Inventor
金姫 陳
金姫 陳
圭一郎 岡本
圭一郎 岡本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tomey Corp
Original Assignee
Tomey Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tomey Corp filed Critical Tomey Corp
Priority to JP2013202023A priority Critical patent/JP6301621B2/ja
Priority to US14/499,062 priority patent/US9265411B2/en
Priority to EP14306491.3A priority patent/EP2865324B1/en
Publication of JP2015066084A publication Critical patent/JP2015066084A/ja
Application granted granted Critical
Publication of JP6301621B2 publication Critical patent/JP6301621B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Description

本発明は、被検眼の前眼部の隅角を含む2次元断層画像を処理する2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法に関する。
近年、眼科検査のために用いられる検査装置として、光干渉断層法(Optical Coherence Tomography:OCT)により被検者の眼球(被検眼)における前眼部の断層画像を撮影する光干渉断層撮影装置(以下「前眼部OCT」という)が供されてきている。
具体的には、前眼部OCTは、例えば、緑内障診療のために用いられるようになっており、主に、疑いを含む原発閉塞隅角症や原発閉塞隅角緑内障を含む狭隅角眼における隅角解析が中心になっている(例えば、非特許文献1参照)。
一般に前眼部OCTでは、測定光を被検眼に対して1次元走査することで1つのスライス面の2次元断層画像を取得し(B−スキャン)、さらに、被検眼に対して測定光の走査位置をずらしながら(換言すれば、スライス面を変えながら)2次元断層画像を繰り返し取得する(C−スキャン)ことで、前眼部3次元画像を得る。
スキャンの方法として、例えば、図4(a)に示すようなラスタースキャンと称される方法がある。このラスタースキャンは、水平方向に延びる走査線に沿って1次元走査(B−スキャン)することを、垂直方向にずらせながら繰り返すこと(C−スキャン)で、図4(b)に示すように、各走査線に沿う2次元断層画像を得ることができる。
また例えば、図5(a)に示すようなラジアルスキャンと称される方法もある。このラジアルスキャンは、放射方向に延びる走査線に沿って1次元走査(B−スキャン)することを、円周方向にずらせながら繰り返すこと(C−スキャン)で、図5(b)に示すように、各走査線に沿う2次元断層画像を得ることができる。
そして、従来の画像処理装置では、こうして得られる各スライス面の2次元断層画像において、強膜岬の位置(SS位置)を検者にポイント入力させることにより、SS位置を超えて閉塞している隅角部分(角膜後面と虹彩前面とが接触している部分)を虹彩線維柱帯接触部(iridotrabecular contact:ITC)としてチャート式に表示したりすることを可能としていた。
三嶋弘一著、「前眼部OCTの緑内障への応用:現在」、あたらしい眼科Vol.28 No.6 P.763〜768(2011年6月号)
しかしながら、従来の画像処理装置では、各2次元断層画像において比較的視認しやすい態様で表示される角膜後面や虹彩前面の位置を手掛かりに、経験的に推定され得るSS位置を検者にポイント入力させる構成であったため、入力に手間がかかったり、入力誤差が生じたりする問題があった。
本発明は、こうした問題に鑑みてなされたものであり、前眼部OCTを用いた隅角解析において、検者の負担を軽減可能とする2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法の提供を目的とする。
(1)本発明の2次元断層画像処理装置は、画像取得手段と、候補抽出手段と、特定手段と、を備える。
画像取得手段は、光干渉断層撮影装置(前眼部OCT)を用いて被検眼の前眼部の隅角を含む2次元断層画像を取得する。
候補抽出手段は、画像取得手段により取得された2次元断層画像の所定領域から、前眼部における強膜とぶどう膜との境界を示す強膜−ぶどう膜エッジ線の候補となる候補エッジ線を抽出する。
特定手段は、候補抽出手段により抽出された候補エッジ線の輝度勾配情報に基づいて、強膜−ぶどう膜エッジ線を特定する。
このように構成された2次元断層画像処理装置によれば、例えば、特定された強膜−ぶどう膜エッジ線を表示することにより、各2次元断層画像において比較的視認しやすい態様で表示される隅角周辺の角膜後面や虹彩前面と、表示された強膜−ぶどう膜エッジ線と、の交点を、SS位置として検者にポイント入力させることが可能となる。
したがって、本発明によれば、検者による入力の手間を軽減させたり、入力誤差を低減させたりすることができる。よって、本発明は、前眼部OCTを用いた隅角解析において、検者の負担を軽減可能とすることができる。
(2)また、本発明の2次元断層画像処理装置は、さらに、角膜後面エッジ線検出手段と、虹彩前面エッジ検出手段と、領域抽出手段と、を備えてもよい。
角膜後面エッジ線検出手段は、2次元断層画像において前眼部の角膜後面を示す角膜後面エッジ線を検出する。
虹彩前面エッジ検出手段は、2次元断層画像において前眼部の虹彩前面を示す虹彩前面エッジ線を検出する。
領域抽出手段は、角膜後面エッジ線検出手段により検出された角膜後面エッジ線と、虹彩前面エッジ線検出手段により検出された虹彩前面エッジ線と、に基づいて、前記の所定領域を抽出する。
このような構成によれば、2次元断層画像から角膜後面エッジ線と虹彩前面エッジ線とを容易に検出することが可能であるため、例えば角膜後面エッジ線と虹彩前面エッジ線との交点と、強膜−ぶどう膜エッジ線と、を含む領域を所定領域として、好適に抽出することが可能となる。そして、このように候補エッジ線を抽出するための領域(所定領域)を2次元断層画像上において限定することにより、候補エッジ線の数を不要に増やさずに済ませることができる。これにより、強膜−ぶどう膜エッジ線を特定する際に、2次元断層画像上の全てのエッジ線の輝度勾配を確認する場合と比べて、処理負担を減らすことができ、さらには強膜−ぶどう膜エッジ線の特定を容易にすることができる。
(3)また、本発明の2次元断層画像処理装置では、角膜後面エッジ線検出手段を備える構成において、さらに、SS位置特定手段を備えてもよい。このSS位置特定手段は、特定手段により特定された強膜−ぶどう膜エッジ線と、角膜後面エッジ線検出手段により検出された角膜後面エッジ線と、の交点を、2次元断層画像において前眼部の強膜岬を示すSS位置として特定する。
このように構成された2次元断層画像処理装置によれば、検者にポイント入力させることなく、2次元断層画像におけるSS位置が自動特定される。すなわち、従来、各2次元断層画像において、SS位置を検者にポイント入力させる構成であったため、前眼部OCTによって複数の2次元断層画像を得ることができたとしても、例えばITCを示すチャートの作成を開始するまでに多大な時間がかかり、臨床上で用いるのが困難であるという問題があった。これに対し、少なくとも検者が全ての2次元断層画像においてSS位置を特定する必要がなくなるため、例えばITCを示すチャートの作成を開始するまでの時間を短縮することができる。よって、このような構成によれば、前眼部OCTを用いた隅角解析において、自動化できる工程を増やすことにより、臨床上で有効に活用可能とすることができる。
また、SS位置の特定方法としては、幾つかの方法が採用され得る。
(4)具体的には、SS位置特定手段は、例えば、特定手段により特定された強膜−ぶどう膜エッジ線と、角膜後面エッジ線検出手段により検出された角膜後面エッジ線と、を繋いで形成されるエッジ線の湾曲点を、SS位置(上記交点)として特定してもよい。このような構成によれば、SS位置を好適に特定することができる。すなわち、2次元断層画像において、強膜−ぶどう膜エッジ線の傾きと、角膜後面エッジ線の傾きと、が異なることを利用し、例えば、上記エッジ線において傾きが大きく変化する点をSS位置として特定することができる。その結果、SS位置を好適に自動特定することができる。
(5)また例えば、SS位置特定手段は、特定手段により特定された強膜−ぶどう膜エッジ線と、角膜後面エッジ線検出手段により検出された角膜後面エッジ線と、を繋いで形成されるエッジ線の輝度勾配情報に基づいて、SS位置(上記交点)を特定してもよい。このような構成によっても、SS位置を好適に特定することができる。すなわち、2次元断層画像の上記所定範囲において、角膜後面エッジ線上の輝度勾配の大きさ(エッジの強さ)が、強膜−ぶどう膜エッジ線上の輝度勾配の大きさを上回ることを利用し、例えば、上記エッジ線において輝度勾配が大きく変化する点をSS位置として特定することができる。その結果、SS位置を好適に自動特定することができる。
(6)また、本発明の2次元断層画像処理装置では、虹彩前面エッジ線検出手段を備える構成において、SS位置特定手段が、特定手段により特定された強膜−ぶどう膜エッジ線と、角膜後面エッジ線検出手段により検出された角膜後面エッジ線と、虹彩前面エッジ線検出手段により検出された虹彩前面エッジ線と、の交点を、SS位置として特定してもよい。
このように構成された2次元断層画像処理装置によれば、例えば、虹彩線維柱帯接触部(ITC)が存在しない2次元断層画像において、3つのエッジ線を用いることにより、上記交点を特定しやすくなるため、SS位置の自動特定における精度を向上させることができる。
(7)また、本発明は、プログラムとして市場に流通させることができる。具体的には、コンピュータを、少なくとも、上記の画像取得手段、候補抽出手段および特定手段として機能させるためのプログラムである。
このプログラムは、1ないし複数のコンピュータに組み込まれることにより、本発明の2次元断層画像処理装置によって奏する効果と同等の効果を得ることができる。なお、本発明のプログラムは、コンピュータに組み込まれるROMやフラッシュメモリ等に記憶され、これらROMやフラッシュメモリ等からコンピュータにロードされて用いられてもよいし、ネットワークを介してコンピュータにロードされて用いられてもよい。
また、上記のプログラムは、コンピュータにて読み取り可能なあらゆる形態の記録媒体に記録されて用いられてもよい。この記録媒体としては、例えば、持ち運び可能な半導体メモリ(例えばUSBメモリ)等が含まれる。
(8)また、本発明は、画像取得手段に相当する工程(画像取得工程)と、候補抽出手段に相当する工程(候補抽出工程)と、特定手段に相当する工程(特定工程)と、を備える2次元断層画像処理方法として表すことができる。この方法を適用することにより、本発明の2次元断層画像処理装置によって奏する効果と同等の効果を得ることができる。
前眼部OCT1の光学系の構成を例示する構成図である。 前眼部OCT1の電気的構成を概略的に例示するブロック図である。 制御装置3が実行するアライメント処理の説明を補足する図である。 ラスタースキャン方式を説明するための図である。 ラジアルスキャン方式を説明するための図である。 第1実施形態において画像処理部100が実行する前眼部3次元画像処理(メイン処理)を説明するためのフローチャートである。 第1実施形態における位置ずれ調整処理の説明を補足する図である。 メイン処理において実施される第1のSS位置特定処理を説明するためのフローチャートである。 第1のSS位置特定処理の説明を補足する図である。 メイン処理において実施される基準真円の算出および第2のSS位置特定処理の説明を補足する図である。 隅角解析の一態様(ITCを示す解析画像)を例示するチャート図である。 第2実施形態において画像処理部100が実行する前眼部3次元画像処理(メイン処理)を説明するためのフローチャートである。 第2実施形態における位置ずれ調整処理の説明を補足する図である。
以下に、本発明の実施形態を図面と共に説明する。
なお、本発明は、下記の実施形態によって何ら限定して解釈されない。また、下記の実施形態の一部を、課題を解決できる限りにおいて省略した態様も本発明の実施形態である。また、特許請求の範囲に記載した文言のみによって特定される発明の本質を逸脱しない限度において考え得るあらゆる態様も本発明の実施形態である。
<第1実施形態>
第1実施形態の前眼部光干渉断層撮影装置は、隅角解析、角膜曲率、角膜厚分布、前房深度の測定等の、被検者の眼球(被検眼E)の前眼部Ec(図1参照)の眼科検査のために用いられる装置であり、光干渉断層法(Optical Coherence Tomography:OCT)により被検眼Eの前眼部Ecの2次元断層画像を撮影することで、3次元画像を得るものである。以下、この前眼部光干渉断層撮影装置を「前眼部OCT1」と称する。
ここで、図示は省略するが、前眼部OCT1の装置本体は、保持台に対して、X方向(左右方向)及びY方向(上下方向)並びにZ方向(前後方向)に移動可能に支持されている。装置本体の前面側(被検者側)には、顎受け部及び額当て部が、上記保持台に対して固定的に設けられている。被検者が、上記顎受け部に顎を載せると共に額当て部に額を当てることにより、被検者の眼(被検眼E)が、装置本体の前面に設けられた撮影用の(光の出入りが行われる)検査窓の正面に配置されるようになっている。
このとき、図2に示すように、この前眼部OCT1には、上記装置本体を保持台に対して、X方向、Y方向、Z方向に夫々自在に移動させるための本体駆動部2が設けられている。この本体駆動部2は、X方向移動モータ、Y方向移動モータ、Z方向移動モータなどを備えた周知構成を備えており、制御装置3により制御されるようになっている。
上記装置本体には、図2に示すように、CPUやメモリ等からなるマイクロコンピュータを含んで構成され全体の制御を行う制御装置3、複数の2次元断層画像によって構成される前眼部Ecの3次元画像(以下「前眼部3次元画像」という)を取得するOCTシステム5、被検眼Eの正面画像を撮影する前眼部撮像系6、アライメント光学系4などが設けられている。
さらに、装置本体には、後面側(検者側)に位置して、被検眼Eの正面画像P(図1参照)等を表示する表示装置としてのモニタ7や、検者(オペレータ)が各種操作を行うための操作部8が設けられている。操作部8には、測定開始スイッチ、測定領域指定スイッチ、モニタ7の画面上に配設されたタッチパネル9、キーボードやマウスなどが設けられている。なお、制御装置3には、撮影された前眼部3次元画像の画像データ等を記憶する記憶部10と、記憶されたデータの画像処理等を行う画像処理部100(前眼部3次元画像処理装置の主要部)が接続されている。記憶部10は、例えばハードディスクドライブを含んで構成されている。
OCTシステム5は、光干渉断層法により前眼部3次元画像を得るものであり、本実施形態では時間的に波長を変化させて操作する波長走査光源11(図1参照)を用いたフーリエドメイン(光周波数掃引)方式が採用されている。
すなわち、図1に示すように、波長走査光源11から出力された光は、光ファイバ12aを通して第1のファイバーカプラ13に入力され、この第1のファイバーカプラ13において、例えば1:99の比率で、参照光と測定光とに分岐され出力される。そのうち参照光は、光ファイバ12bを通って第1のサーキュレータ14の入力部に入力され、更にこの第1のサーキュレータ14の入出力部から光ファイバ12cを通ってその端部から出力され、複数個のコリメータレンズ15を通って参照ミラー16に入射される。
そして、参照ミラー16にて反射された参照光が、再び、複数個のコリメータレンズ15を通って光ファイバ12cの端部から入力され、光ファイバ12cを通って第1のサーキュレータ14の入出力部から入力される。そして、第1のサーキュレータ14の出力部から出力された参照光は、光ファイバ12dを通って第2のファイバーカプラ17の第1の入力部に入力される。
一方、第1のファイバーカプラ13から出力された測定光は、光ファイバ12eを通って第2のサーキュレータ18の入力部に入力され、さらにこの第2のサーキュレータ18の入出力部から光ファイバ12fを通ってその端部から出力される。光ファイバ12fの端部から出力された測定光は、コリメータレンズ19を通ってガルバノスキャナ20に入力される。ガルバノスキャナ20は、測定光を走査させるためのもので、ガルバノドライバ21により駆動されるようになっている。
ガルバノスキャナ20から出力された測定光は、波長が長い側の光を反射させ短い側の光を透過させるホットミラー22により90度の角度で反射され、対物レンズ23を通して上記検査窓から出射され、被検眼Eに入射される。被検眼Eに入射された測定光は、前眼部Ecの各組織部分(角膜、前房、虹彩、水晶体、ぶどう膜、強膜等)にて反射し、その反射光が、検査窓から入射され、上記とは逆に、対物レンズ23、ホットミラー22、ガルバノスキャナ20、コリメータレンズ19を順に通って、光ファイバ12fの端部から入力される。そして、その反射光は、光ファイバ12fを通って第2のサーキュレータ18の入出力部から入力され、第2のサーキュレータ18の出力部から出力され、光ファイバ12gを通って第2のファイバーカプラ17の第2の入力部に入力される。
この第2のファイバーカプラ17において、前眼部Ecからの反射光と、光ファイバ12dを通って入力された参照光とが、例えば50:50の比率で合波され、その信号が光ファイバ12h、12iを介して検出器24に入力される。検出器24においては、波長毎の干渉が計測され、計測された干渉信号が、制御装置3に設けられたADボード25に入力される。さらに、制御装置3に設けられた演算部26において、干渉信号に対するフーリエ変換などの処理が行われ、もって走査線に沿う前眼部Ecの断層画像(2次元断層画像)が得られるのである。
このとき、ガルバノスキャナ20による測定光のスキャンパターン、換言すれば、走査線(B−スキャン)の方向は、制御装置3において設定されるようになっている。そして、制御装置3(演算部26)からの指令信号に基づいてガルバノドライバ21がガルバノスキャナ20を制御するようになっている。なお、このように得られた2次元断層画像の画像データは、記憶部10に記憶される。この2次元断層画像の画像データには、少なくとも各画素の輝度を示す情報が含まれている。また、図1に模式的に示しているように、その断層画像Tをモニタ7に表示させることができる。
次に、前眼部撮像系6は、照明光源27,27、対物レンズ23、ホットミラー22、コールドミラー28、結像レンズ29、CCDカメラ30、光学制御部31を備えて構成される。照明光源27,27は、被検眼Eの正面に可視光領域の照明光を照射するようになっており、被検眼Eからの反射光が、上記検査窓から対物レンズ23、ホットミラー22、コールドミラー28、結像レンズ29を通って、CCDカメラ30に入力される。これにて、被検眼Eの正面画像Pが撮影され、撮影された画像データは、光学制御部31によって画像処理が行われて、モニタ7に表示されるようになる。
そして、アライメント光学系4は、被検者が固視灯を見つめることにより眼球(被検眼E)を動かさないようにさせるための固視灯光学系、被検眼Eの角膜頂点のXY方向の位置(装置本体に対する上下左右の位置ずれ)を検出するためのXY方向位置検出系、被検眼Eの角膜頂点の前後方向(Z方向)の位置を検出するためのZ方向位置検出系を含んで構成されている。
そのうち固視灯光学系は、固視灯32、コールドミラー33、リレーレンズ34、ハーフミラー35、コールドミラー28、ホットミラー22、対物レンズ23などから構成されている。これにて、固視灯32から出力された光(例えば緑色の光)は、コールドミラー33、リレーレンズ34、ハーフミラー35、コールドミラー28、ホットミラー22、レンズ23を順に介して、検査窓から被検眼Eに向けて出力されるようになっている。
上記XY方向位置検出系は、XY位置検出光源36、コールドミラー33、リレーレンズ34、ハーフミラー35、コールドミラー28、ホットミラー22、対物レンズ23、結像レンズ37、位置センサ38などを備えて構成されている。これにて、XY位置検出光源36からは、位置検出用のアライメント光が出力され、コールドミラー33、リレーレンズ34、ハーフミラー35、コールドミラー28、ホットミラー22、対物レンズ23を介して、検査窓から被検眼Eの前眼部Ec(角膜)に向けて出射される。
このとき、被検眼Eの角膜表面が球面状をなすことにより、アライメント光は、被検眼Eの角膜頂点の内側で輝点像を形成するようにして角膜表面で反射され、その反射光が、検査窓から入射されるようになっている。角膜頂点からの反射光(輝点)は、対物レンズ23、ホットミラー22、コールドミラー28、ハーフミラー35、結像レンズ37を介して位置センサ38に入力される。位置センサ38によってその輝点の位置が検出されることにより、角膜頂点の位置(X方向及びY方向の位置)が検出されるようになっている(図3(a)参照)。なお、上記輝点は、CCDカメラ30の撮影画像(モニタ7の表示画像)にも写り込むものとなる。
位置センサ38の検出信号は、光学制御部31を介して制御装置3(演算部26)に入力される。この制御装置3の演算部26には、本実施形態では前眼部3次元画像処理装置としての一部機能を実現するためのプログラムがメモリまたは記憶部10に実装されており、演算部26において、CPUがこのプログラムに従ってアライメント処理を実施する。すなわち、このアライメント処理では、位置センサ38の検出結果に基づいて、角膜頂点の所定(正規)の画像取得位置に対する、検出された角膜頂点(輝点)のX方向及びY方向の位置ずれ量ΔX,ΔYを求める。
Z方向位置検出系は、Z方向位置検出光源39、結像レンズ40、ラインセンサ41を備えて構成されている。Z方向位置検出光源39は、被検眼Eに対して斜め方向から検出用の光(スリット光又はスポット光)を照射し、角膜からの斜め方向の反射光が、結像レンズ40を介してラインセンサ41に入射されるようになっている。このとき、装置本体に対する被検眼Eの前後方向(Z方向)の位置によって、ラインセンサ41に入射される反射光の入射位置が異なるようになるので、被検眼Eの装置本体に対するZ方向の位置(距離)が検出されるのである(図3(b)参照)。
ラインセンサ41の検出信号は、光学制御部31を介して制御装置3(演算部26)に入力される。このとき、被検眼Eの角膜頂点の装置本体に対する適切なZ方向位置(距離)が予め設定されており、制御装置3の演算部26は、上記アライメント処理において、ラインセンサ41の検出結果に基づいて、被検眼Eの適切な位置としての角膜頂点の位置に対するZ方向のずれ量ΔZを求める。
そして、アライメント処理では、制御装置3の演算部26は、XY方向位置検出系により検出された角膜頂点のX方向及びY方向の位置ずれ量ΔX,ΔY、並びに、Z方向位置検出系により検出された被検眼EのZ方向の位置ずれ量ΔZを、アライメント情報として、対応する2次元断層画像の画像データが識別可能となる保存形式により記憶部10に記憶する。
つまり、制御装置3の演算部26は、ガルバノスキャナ20を制御して、測定光を被検眼Eに対して1次元走査することで1つのスライス面の2次元断層画像を取得し(B−スキャン)、さらに、被検眼Eに対して測定光の走査位置をずらしながら(換言すれば、スライス面を変えながら)2次元断層画像を繰り返し取得する(C−スキャン)ことで得られた前眼部3次元画像を記憶部10に記憶する。さらに、この前眼部3次元画像を構成する夫々の2次元断層画像に関する上記のアライメント情報を記憶部10に記憶するようになっている。
なお、スキャンの方法として、既述のとおり、図4に示すラスタースキャンと称される方法や、図5に示すラジアルスキャンと称される方法があり、操作部8を介して検者により測定対象が選択された結果に応じて適切な方法が選択される。本実施形態では、測定対象として隅角解析が選択されると、制御装置3の演算部26は、スキャンパターンとしてラジアルスキャンを採用し、具体的には、被検眼Eの角膜頂点を中心とする放射方向をB−スキャン方向、被検眼Eの前眼部Ecの表面円周方向をC−スキャン方向として、各スライス面の2次元断層画像の取り込みを行うようになっている。以下、こうして取り込まれて記憶部10に記憶される各スライス面の2次元断層画像には、前眼部Ecの隅角が2箇所含まれているものとして説明する。
<前眼部3次元画像処理(第1実施形態)>
画像処理部100は、CPUやメモリ等からなるマイクロコンピュータを含んで構成され、前眼部3次元画像処理装置としての主要機能を実現するためのプログラムがメモリまたは記憶部10に実装されており、CPUがこのプログラムに従って図6に示す前眼部3次元画像処理におけるメイン処理を実施する。すなわち、このメイン処理では、S10において、画像処理部100が、記憶部10から前眼部3次元画像を構成する各スライス面の2次元断層画像を取得する。なお、各スライス面は、測定光の光軸を基準として隣接するスライス面とラジアルスキャンのC−スキャン方向に所定の角度をなすように予め設定されている。本実施形態の設定角度は、11.25度に設定されている。つまり、B−スキャン方向の数は32方向となり、16枚の2次元断層画像が取得されることになる。
続いて、S20では、画像処理部100が、S10において取得した16枚の各スライス面の2次元断層画像について、2次元断層画像上における空間座標位置のずれの有無を判定する処理(位置ずれ判定処理)を実施する。この位置ずれ判定処理では、記憶部10に記憶されているアライメント情報を用いて判定する他、各2次元断層画像間で後述する角膜前面曲線のずれが存在するか否か(または、ずれが大きいか否か)に応じて判定する。すなわち、各スライス面の2次元断層画像について、記憶部10に記憶されているアライメント情報に基づき、位置ずれ量ΔX,ΔY,ΔZの少なくとも一つが所定の許容閾値を上回るか否かに応じて、許容閾値を上回る場合には、空間座標位置のずれが存在すると判定し、許容閾値以下である場合には、角膜前面曲線のずれに対する閾値判断を行う。そして、角膜前面曲線のずれが有ると判断した場合には、空間座標位置のずれが存在すると判定し、角膜前面曲線のずれが無いと判断した場合には、空間座標位置のずれが存在しないと判定する。なお、本実施形態では、このように、アライメント情報を用いて判定する手法と、角膜前面曲線のずれに応じた判定を行う手法と、の両方の手法を用いているが、いずれか一方の手法を用いるだけでもよい。
次に、S30では、画像処理部100が、S20における判定結果に応じて処理を分岐する。すなわち、S20において、空間座標位置のずれが存在すると判定した場合、S40に移行し、空間座標位置のずれが存在しないと判定した場合、S50に移行する。
S40では、S20において空間座標位置のずれが存在すると判定した2次元断層画像について、その空間座標位置のずれを調整する処理(位置ずれ調整処理)を実施する。本実施形態の位置ずれ調整処理では、上記のアライメント情報に基づく位置ずれ量ΔX,ΔY,ΔZに対し、2次元断層画像上における空間座標位置のオフセット量をΔX´,ΔZとすると、例えば以下の関係式(1)を満たすように、オフセット量ΔX´を求める。なお、オフセット量ΔX´は、図7(a)に示すように、Z方向に対し、2次元断層画像上において垂直な方向をX´方向とした場合の、X´方向における空間座標位置の補正量である。また、θscanとは、図7(b)に示すように、X方向に対し、ラジアルスキャンのB−スキャン方向がなす角度のことである。
なお、上記関係式(1)は、オフセット量ΔX´が微小(例えば、300μm以下)である場合に近似式として用いられるものである。また、この位置ずれ調整処理では、上記のように記憶部10に記憶されているアライメント情報を用いて調整(補正)する他、各2次元断層画像間で後述する角膜前面曲線のずれを補正する(第2実施形態参照)。S40では、このように、記憶部10に記憶されている全ての2次元断層画像について、位置ずれ調整処理を行うことにより、各画像の空間座標位置を合わせ、前眼部3次元画像を再構築する。なお、本実施形態では、上記のように、アライメント情報を用いて補正する手法と、角膜前面曲線のずれを補正する手法と、の両方の手法を用いているが、いずれか一方の手法を用いるだけでもよい。但し、両方の手法を用いることにより、相互に性質の異なる誤差を補完することができる。例えば、アライメント情報を用いて補正する手法では、被検眼E(眼球)の回旋運動を考慮できないことによる誤差が生じる可能性があり、角膜前面曲線のずれを補正する手法では、眼球が大きく動いてしまった場合に誤差が生じる可能性があり、このように性質の異なる誤差を補完することができる。
次に、S50では、画像処理部100が、S10において取得した16枚の各スライス面の2次元断層画像のうち、4枚の2次元断層画像を代表画像とし、各代表画像において前眼部Ecの強膜岬の空間座標位置を示すSS位置を2箇所ずつ特定する処理(以下「第1のSS位置特定処理」という)を行う。つまり、本実施形態では、4枚の代表画像から8箇所のSS位置を特定する。なお、本実施形態では、この4枚の代表画像として、互いにスライス面のなす角度がある所定角度(例えば30度)以上となる4枚の2次元断層画像を選択する。
具体的に、このS50の処理内容を図8に示す。図8に示す第1のSS位置特定処理では、S110において、画像処理部100が、上記代表画像から前眼部Ecの隅角近傍を局所的に含む画像(以下「局所画像」という:図9(a)参照)を抽出する。
続いて、S120では、画像処理部100が、S110において抽出された局所画像の画像データの各画素について、例えばZ方向に隣接する画素との輝度の差分等を求めることにより輝度勾配を算出する。
次に、画像処理部100は、上記の局所画像における輝度勾配に基づき、S130において、前眼部Ecにおける角膜後面を示すエッジ線(以下「角膜後面エッジ線」という)を検出し、S140において、前眼部Ecにおける虹彩前面を示すエッジ線(以下「虹彩前面エッジ線」という)を検出する。すなわち、上記局所画像の画像データにおいては、角膜後面エッジ線および虹彩前面エッジ線上の画素の輝度勾配が最も高い。このため、例えば、輝度勾配の閾値を適宜設定することにより、局所画像から角膜後面エッジ線や虹彩前面エッジ線を抽出(検出)することができる。
続いて、画像処理部100は、S150において、局所画像の画像データにおいて上記閾値を可変設定することにより、角膜後面エッジ線や虹彩前面エッジ線の他、前眼部Ecにおける各部位を規定する可能性のあるエッジ線を含む画像データ(以下「エッジ画像」という:図9(b)参照)を生成する。そして、S155において、このエッジ画像から、S130において検出された角膜後面エッジ線上にて、SS位置を含むと推定される所定領域を抽出する。例えば、このS155では、S140において虹彩前面エッジ線が検出できた場合、角膜後面エッジ線と虹彩前面エッジ線との2つのエッジ線を繋いで形成されるエッジ線の湾曲点(前眼部Ecの隅角底に相当)を基にして、上記エッジ画像から上記所定領域を限定する。
そして、S160では、画像処理部100が、上記エッジ画像のうち、S150において抽出された所定領域外のエッジ線を不要エッジ線として除去する。例えば、所定領域外の角膜後面エッジ線から分岐している不要エッジ線や、所定領域外の虹彩前面エッジ線から分岐している不要エッジ線を除去する(図9(c)参照)。
S170では、画像処理部100が、S160における不要エッジ線の除去により、前眼部Ecにおける強膜とぶどう膜との境界を示すエッジ線(以下「強膜−ぶどう膜エッジ線」という)の候補となる候補エッジ線を抽出する。
続いて、S180では、画像処理部100が、S170において抽出した各候補エッジ線について、交差方向の輝度勾配の大きさ(エッジの強さ)を算出し、各候補エッジ線のうち、この輝度勾配が最大となるエッジ線を、強膜−ぶどう膜エッジ線として特定する。
そして、S190では、画像処理部100が、S140において虹彩前面エッジ線を検出できたか否かを判定し、その判定結果に応じて処理を分岐する。つまり、被検眼Eによっては、隅角が閉塞されている場合があり、このような場合に、虹彩前面エッジ線が角膜後面エッジ線と一体化されたように映し出され、虹彩前面エッジ線が検出されない可能性があるためである。ここで、画像処理部100が、虹彩前面エッジ線を検出できたと判定した場合、S200に移行し、虹彩前面エッジ線を検出できなかったと判定した場合、S210に移行する。
S200では、画像処理部100が、S180において特定した強膜−ぶどう膜エッジ線と、S130において検出した角膜後面エッジ線と、S140において検出した虹彩前面エッジ線と、の交点を示す空間座標位置を、SS位置として特定し、S220に移行する。
一方、S210では、画像処理部100が、S180において特定した強膜−ぶどう膜エッジ線と、S130において検出した角膜後面エッジ線と、の交点を示す空間座標位置を、SS位置として特定し、S200に移行する。具体的には、強膜−ぶどう膜エッジ線と角膜後面エッジ線との交点(つまりSS位置)の特定方法としては、幾つかの方法が採用され得る。例えば、強膜−ぶどう膜エッジ線と角膜後面エッジ線との両方のエッジ線を繋いで形成されるエッジ線(以下「対象エッジ線」という)の形状に基づいて、SS位置を特定することが可能である。すなわち、エッジ画像において、強膜−ぶどう膜エッジ線の傾きと、角膜後面エッジ線の傾きと、が異なることを利用し、例えば、上記対象エッジ線において湾曲するように傾きが大きく変化する点(湾曲点)をSS位置として特定することができる。また例えば、上記対象エッジ線上の輝度勾配の情報に基づいて、SS位置を特定することも可能である。すなわち、エッジ画像において、角膜後面エッジ線上の輝度勾配が、強膜−ぶどう膜エッジ線上の輝度勾配よりも高いことを利用し、例えば、上記対象エッジ線において輝度勾配が大きく変化する点をSS位置として特定することができる。
S220では、画像処理部100が、全ての代表画像(本実施形態では4枚の代表画像)からそれぞれ所定数のSS位置(本実施形態では2箇所のSS位置)を特定できたか否かを判定し、全てのSS位置(本実施形態では計8箇所のSS位置)を特定できた場合には、メイン処理(S60)に戻り、代表画像において未特定のSS位置が存在する場合には、S110に戻って第1のSS位置特定処理を継続する。
ここでメイン処理に戻ると、S60では、画像処理部100が、S50において特定された複数(本実施形態では8点)のSS位置のうち、少なくとも3点のSS位置を通る基準真円(図10参照)を空間座標上に示す関数を算出する。具体的には、本実施形態では、8点のSS位置のうち、少なくとも3点のSS位置を通り、残りのSS位置との距離が最小となる(換言すれば、残りのSS位置が上記関数上に最も近似するように配置される)空間平面上の基準真円が求められる。このように8点のSS位置のうち、残りのSS位置が近似的に配置されるように、基準真円を求めることにより、画像間での誤差を適度に分散することが可能となり、SS位置の自動特定における精度を向上させることができる。なお、基準真円としては、一般に真円が採用されるが、完全な真円の他、真円に近いものが採用される可能性もある。
S70では、画像処理部100が、前眼部3次元画像を構成する複数(本実施形態では16枚)の2次元断層画像のうち、複数(本実施形態では4枚)の代表画像以外の複数(本実施形態では12枚)の画像(以下「非代表画像」という)について、S60において算出した基準真円の関数に基づいて、残りのSS位置を特定する処理(以下「第2のSS位置特定処理」ともいう)を行う。具体的には、S60において求めた基準真円上に、各非代表画像においてB−スキャン方向に対応する各点を、それぞれの非代表画像におけるSS位置として特定する。これにより、画像処理部100は、メイン処理を終了する。
なお、画像処理部100は、このように全てのスライス面において求められたSS位置を用いることにより、例えば、SS位置を超えて閉塞している隅角部分EP(角膜後面と虹彩前面とが接触している部分)を虹彩線維柱帯接触部(ITC)としてチャート式に示す解析画像(図11参照)を生成したりすることができる。そして、これらの画像は、操作部8を介した検者による操作指示に応じてモニタ7に出力される。
<第2実施形態>
次に、本発明の第2実施形態について説明する。なお、第2実施形態は、第1実施形態に対して、画像処理部100が実施するメイン処理(前眼部3次元画像処理)の内容が異なるだけであるため、その他の内容については説明を省略する。具体的には、第1実施形態の前眼部3次元画像処理では、位置ずれ調整処理(S40)を行うことにより、各2次元断層画像の空間座標位置を合わせ、前眼部3次元画像を再構築し、再構築された前眼部3次元画像を構成する各2次元断層画像に対してSS位置を特定していた(S50〜S70)。これに対して、第2実施形態の前眼部3次元画像処理では、前眼部3次元画像の再構築を行うことなく、各2次元断層画像の空間座標位置のずれを調整するために算出したパラメータを用いてSS位置の決定を行う点で異なる。第2実施形態の前眼部3次元画像処理によれば、前眼部3次元画像を再構築しないで済む分、隅角解析を含む全体の処理速度を向上することができる。
<前眼部3次元画像処理(第2実施形態)>
すなわち、図12に示す第2実施形態のメイン処理では、S1000において、第1実施形態と同様、画像処理部100が、記憶部10から前眼部3次元画像を構成する各スライス面の2次元断層画像を取得する。なお、第2実施形態において各スライス面は、測定光の光軸を基準として隣接するスライス面とラジアルスキャンのC−スキャン方向に所定の角度をなすように予め設定されている。本実施形態の設定角度は、5.625度に設定されている。つまり、B−スキャン方向の数は64方向となり、32枚の2次元断層画像が取得されることになる。
続いて、S2000では、画像処理部100が、S1000において取得した32枚の各スライス面の2次元断層画像について、1枚毎に、隣接する他の1枚のスライス面の2次元断層画像との間で、相互に空間座標位置を角膜前面の位置を基準にして合わせるための動き行列Vを算出する。具体的には、まず、1枚の2次元断層画像とこれに隣接する1枚の2次元断層画像との夫々から例えばパターンマッチング等の周知技術を用いて角膜前面形状を示す曲線(以下「角膜前面曲線」という)を抽出する。そして、例えば、図13に示すように、抽出された夫々の角膜前面曲線について、一方を他方側へ互いに並進および回転させ、その並進距離Tおよび回転角Rが最小となる基準の角膜前面曲線(以下「基準曲線」という)を求め、この基準曲線を出力値とする各2次元断層画像の角膜前面曲線の入力式を夫々の動き行列Vとして算出する。例えば、この動き行列Vの算出を全ての2次元断層画像について行うとともに、全ての2次元断層画像に関する平均の基準曲線を求め、この平均の基準曲線を基にして各2次元断層画像に関する動き行列Vを補正する。こうして算出された動き行列Vは、夫々対応するスライス面の2次元断層画像に対応付けられてメモリに一時記憶される。なお、このS2000は、全ての2次元断層画像について、各2次元断層画像間で角膜前面曲線のずれが存在するか否か(または、ずれが大きいか否か)を判定し、角膜前面曲線のずれが存在する場合に行われるものとする。
次に、S3000では、画像処理部100が、S1000において取得し、S2000において動き行列Vが算出された32枚の各スライス面の2次元断層画像のうち、3枚の2次元断層画像を代表画像とし、各代表画像において前眼部Ecの強膜岬の空間座標位置を示すSS位置を2箇所ずつ特定する処理(第1のSS位置特定処理)を行う。つまり、第2実施形態では、3枚の代表画像から6箇所のSS位置を特定する。なお、第2実施形態では、この3枚の代表画像として、互いにスライス面のなす角度がある所定角度(例えば45度)以上となる3枚の2次元断層画像を選択する。この第1のSS位置特定処理の内容については、第1実施形態と同様であるため、説明を省略する。
S4000では、画像処理部100が、S3000(第1のSS位置特定処理)において特定された複数(本実施形態では6点)のSS位置について、S2000において算出した動き行列Vを用いてそれぞれ空間座標位置を調整(補正)する。
S5000では、画像処理部100が、S4000において補正された複数(本実施形態では6点)のSS位置(以下「SS´位置」という)のうち、少なくとも3点のSS´位置を通る基準真円(図10参照)を空間座標上に示す関数を算出する。具体的には、第2実施形態では、1枚の代表画像により特定された2点のSS´位置の距離を直径とし、少なくとも残り1点のSS´位置を通る空間平面上の基準真円が求められる。すなわち、上記直径を構成する2点のSS´位置の他に、SS´位置を少なくとも1点特定するだけで、空間平面上の基準真円が定まるため、基準真円を求めるために用いる2次元断層画像(代表画像)の数を少なく済ませることが可能となり、自動化できる工程をより増やすことができる。
S6000では、画像処理部100が、前眼部3次元画像を構成する複数(本実施形態では32枚)の2次元断層画像のうち、複数(本実施形態では3枚)の代表画像以外の複数(本実施形態では29枚)の画像(以下「非代表画像」という)について、S5000において算出した基準真円の関数に基づいて、残りのSS´位置を特定する処理(第2のSS位置特定処理)を行う。具体的には、S5000において求めた基準真円上に、各非代表画像においてB−スキャン方向に対応する各点を、それぞれの非代表画像におけるSS´位置として特定する。
そして、S7000では、画像処理部100が、こうして特定された全ての2次元断層画像におけるSS´位置について、S2000において算出した動き行列Vを用いて、補正前のSS位置に戻すことにより、全ての2次元断層画像におけるSS位置を算出(特定)する。これにより、画像処理部100は、メイン処理を終了する。
<主要な効果>
以上説明したように、前眼部OCT1では、第1のSS位置特定処理において、2次元断層画像から前眼部Ecの隅角を含む局所画像(図9(a)参照)を取得し、この局所画像の所定領域から、前眼部における強膜とぶどう膜との境界を示す強膜−ぶどう膜エッジ線の候補となる候補エッジ線を抽出する(図9(b)〜(c)参照)。そして、抽出された候補エッジ線の輝度勾配情報に基づいて、強膜−ぶどう膜エッジ線を特定する(図9(d)参照)。
また、前眼部OCT1では、第1のSS位置特定処理において、前眼部Ecの角膜後面を示す角膜後面エッジ線を検出(特定)し、特定された強膜−ぶどう膜エッジ線と角膜後面エッジ線との交点を、上記2次元断層画像において前眼部Ecの強膜岬を示すSS位置として特定する。
このため、前眼部OCT1によれば、検者にポイント入力させることなく、2次元断層画像におけるSS位置が自動特定されることから、例えばITCを示すチャートの作成を開始するまでの時間を短縮することができる。よって、このような構成によれば、前眼部OCTを用いた隅角解析において、自動化できる工程を増やすことにより、臨床上で有効に活用可能とすることができる。
また、前眼部OCT1では、第1のSS位置特定処理において、前眼部Ecの虹彩前面を示す虹彩前面エッジ線を検出(特定)できた場合、特定された強膜−ぶどう膜エッジ線と角膜後面エッジ線と虹彩前面エッジ線との交点を、上記SS位置として特定する。このため、例えば、虹彩線維柱帯接触部(ITC)が存在しない2次元断層画像において、3つのエッジ線を用いることにより、上記交点を抽出しやすくなるため、SS位置の自動特定における精度を向上させることができる。
また、前眼部OCT1では、メイン処理において、空間座標位置のずれの有無の判定(例えばS10〜S20)を経た2次元断層画像のうち、2つ以上の代表画像を用いて、3点以上のSS位置の特定を自動的に受け付け(例えばS50)、少なくとも3点のSS位置を通る基準真円を空間座標上に示す関数を算出する(例えばS60)。そして、上記代表画像以外の2次元断層画像(非代表画像)におけるSS位置等(残りのSS位置)を、上記基準真円の関数に基づいて特定する(例えばS70)。
このため、前眼部OCT1によれば、検者がSS位置のポイント入力を全く行わなくてもよいことから、例えばITCを示すチャートの作成を開始するまでの時間を大幅に短縮することができる。従って、前眼部OCT1によれば、前眼部OCTを用いた隅角解析において、全工程を自動化することにより、臨床上で有効に活用可能とすることができる。
<他の実施形態>
以上、本発明の第1および第2実施形態について説明したが、本発明はこれらの実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において、様々な態様にて実施することが可能である。
例えば、上記実施形態では、メイン処理において、第1のSS位置特定処理を行うことにより、SS位置を自動特定しているが、これに代えて、ここでは、操作部8を介した検者によるSS位置のポイント入力を受け付けるようにしてもよい。具体的には、第1のSS位置特定処理において、特定された強膜−ぶどう膜エッジ線をモニタ7に表示させることにより、各2次元断層画像において比較的視認しやすい態様で表示される隅角周辺の角膜後面や虹彩前面と、表示された強膜−ぶどう膜エッジ線と、の交点を、SS位置とし、操作部8を介して検者にポイント入力させてもよい。このような構成によっても、検者による入力の手間を軽減させたり、入力誤差を低減させたりすることができるため、前眼部OCT1を用いた隅角解析において、検者の負担を軽減可能とすることができる。
また、上記第2実施形態では、前眼部3次元画像処理(メイン処理)において、角膜前面曲線を用いたSS位置の補正(S2000,S4000)を行って全てのSS´位置を特定した後に、補正前のSS位置に戻す処理(S7000)を行っているが、必ずしもこのようにしなければならないわけではない。例えば、記憶部10に記憶されている全ての2次元断層画像について、このような補正を行うことにより、各画像の空間座標位置を合わせ、前眼部3次元画像を再構築することもできる。
なお、上記第1実施形態では、制御装置3が前眼部3次元画像を構成する各スライス面の2次元断層画像を記憶部10に記憶させるように構成されているが、記憶部10として、例えばインターネット上のサーバに記憶させるようにしてもよい。また、上記第1実施形態では、前眼部OCT1において、制御装置3と画像処理部100とが別体に構成されているが、画像処理部100が行う処理を制御装置3が行うようにしてもよいし、制御装置3が行う処理を画像処理部100が行うようにしてもよい。さらには、前眼部OCT1とは別体に、画像処理部100を備える装置が構成され、この装置が前眼部OCT1や上記サーバとの間で通信可能に接続されることで、各種処理を行うようにしてもよい。また、この装置に各種処理を実行させるためのプログラムは、記憶部10や上記サーバに記憶されており、画像処理部100がこのプログラムをロードして各種処理を実行するようにしてもよい。
1…前眼部OCT、2…本体駆動部、3…制御装置、4…アライメント光学系、5…OCTシステム、6…前眼部撮像系、7…モニタ、8…操作部、9…タッチパネル、10…記憶部、11…波長走査光源、12a〜12h…光ファイバ、13…第1のファイバーカプラ、14…第1のサーキュレータ、15…コリメータレンズ、16…参照ミラー、17…第2のファイバーカプラ、18…第2のサーキュレータ、19…コリメータレンズ、20…ガルバノスキャナ、21…ガルバノドライバ、22…ホットミラー、23…対物レンズ、24…検出器、25…ADボード、26…演算部、27…照明光源、28…コールドミラー、29…結像レンズ、30…CCDカメラ、31…光学制御部、32…固視灯、33…コールドミラー、34…リレーレンズ、35…ハーフミラー、36…位置検出光源、37…結像レンズ、38…位置センサ、39…方向位置検出光源、40…結像レンズ、41…ラインセンサ、100…画像処理部、E…被検眼、Ec…前眼部、P…正面画像、T…断層画像。

Claims (6)

  1. 光干渉断層撮影装置を用いて被検眼の前眼部の隅角を含む2次元断層画像を取得する画像取得手段と、
    前記画像取得手段により取得された前記2次元断層画像の所定領域から、前記前眼部における強膜とぶどう膜との境界を示す強膜−ぶどう膜エッジ線の候補となる候補エッジ線を抽出する候補抽出手段と、
    前記候補抽出手段により抽出された前記候補エッジ線の輝度勾配情報に基づいて、前記強膜−ぶどう膜エッジ線を特定する特定手段と、
    前記2次元断層画像において前記前眼部の角膜後面を示す角膜後面エッジ線を検出する角膜後面エッジ線検出手段と、
    前記特定手段により特定された前記強膜−ぶどう膜エッジ線と、前記角膜後面エッジ線検出手段により検出された前記角膜後面エッジ線と、の交点を、前記2次元断層画像において前記前眼部の強膜岬を示すSS位置として特定するSS位置特定手段と、
    を備えることを特徴とする2次元断層画像処理装置。
  2. 記2次元断層画像において前記前眼部の虹彩前面を示す虹彩前面エッジ線を検出する虹彩前面エッジ線検出手段と、
    前記角膜後面エッジ線検出手段により検出された前記角膜後面エッジ線と、前記虹彩前面エッジ線検出手段により検出された前記虹彩前面エッジ線と、に基づいて、前記所定領域を抽出する領域抽出手段と、
    を備えることを特徴とする請求項1に記載の2次元断層画像処理装置。
  3. 前記SS位置特定手段は、前記特定手段により特定された前記強膜−ぶどう膜エッジ線と、前記角膜後面エッジ線検出手段により検出された前記角膜後面エッジ線と、を繋いで形成されるエッジ線において輝度勾配が大きく変化する点を前記SS位置として特定することを特徴とする請求項または請求項に記載の2次元断層画像処理装置。
  4. 前記2次元断層画像において前記前眼部の虹彩前面を示す虹彩前面エッジ線を検出する虹彩前面エッジ線検出手段を備え、
    前記SS位置特定手段は、前記特定手段により特定された前記強膜−ぶどう膜エッジ線と、前記角膜後面エッジ線検出手段により検出された前記角膜後面エッジ線と、前記虹彩前面エッジ線検出手段により検出された前記虹彩前面エッジ線と、の交点を、前記SS位置として特定することを特徴とする請求項ないし請求項のいずれか1項に記載の2次元断層画像処理装置。
  5. コンピュータを、少なくとも請求項1に記載の前記画像取得手段、前記候補抽出手段前記特定手段、前記角膜後面エッジ線検出手段および前記SS位置特定手段として機能させるためのプログラム。
  6. 2次元断層画像処理装置が実行する2次元断層画像処理方法であって、
    光干渉断層撮影装置に記憶された被検眼の前眼部の隅角を含む2次元断層画像を取得する画像取得工程と、
    前記画像取得工程により取得された前記2次元断層画像の所定領域から、前記前眼部における強膜とぶどう膜との境界を示す強膜−ぶどう膜エッジ線の候補となる候補エッジ線を抽出する候補抽出工程と、
    前記候補抽出工程により抽出された前記候補エッジ線の輝度勾配情報に基づいて、前記強膜−ぶどう膜エッジ線を特定する特定工程と、
    前記2次元断層画像において前記前眼部の角膜後面を示す角膜後面エッジ線を検出する角膜後面エッジ線検出工程と、
    前記特定工程により特定された前記強膜−ぶどう膜エッジ線と、前記角膜後面エッジ線検出工程により検出された前記角膜後面エッジ線と、の交点を、前記2次元断層画像において前記前眼部の強膜岬を示すSS位置として特定するSS位置特定工程と、
    を備えることを特徴とする2次元断層画像処理方法。
JP2013202023A 2013-09-27 2013-09-27 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法 Active JP6301621B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2013202023A JP6301621B2 (ja) 2013-09-27 2013-09-27 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法
US14/499,062 US9265411B2 (en) 2013-09-27 2014-09-26 Anterior segment three-dimensional image processing apparatus, and anterior segment three-dimensional image processing method
EP14306491.3A EP2865324B1 (en) 2013-09-27 2014-09-26 Anterior segment three-dimensional image processing apparatus, and anterior segment three-dimensional image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013202023A JP6301621B2 (ja) 2013-09-27 2013-09-27 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法

Publications (2)

Publication Number Publication Date
JP2015066084A JP2015066084A (ja) 2015-04-13
JP6301621B2 true JP6301621B2 (ja) 2018-03-28

Family

ID=52833520

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013202023A Active JP6301621B2 (ja) 2013-09-27 2013-09-27 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法

Country Status (1)

Country Link
JP (1) JP6301621B2 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6436888B2 (ja) * 2015-10-19 2018-12-12 株式会社トーメーコーポレーション 眼内レンズ度数決定装置
JP6739725B2 (ja) 2015-11-17 2020-08-12 株式会社トーメーコーポレーション 前眼部3次元画像処理装置および前眼部3次元画像処理方法
US10779729B2 (en) 2016-12-05 2020-09-22 Nidek Co., Ltd. Ophthalmic apparatus
JP7091627B2 (ja) * 2016-12-05 2022-06-28 株式会社ニデック 眼科装置
JP2023017402A (ja) * 2021-07-26 2023-02-07 株式会社トプコン 前眼部解析装置、前眼部解析方法、およびプログラム
JP2023111652A (ja) 2022-01-31 2023-08-10 株式会社トーメーコーポレーション 断層画像処理装置及びプログラム

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004006751A2 (en) * 2002-07-12 2004-01-22 Volker Westphal Method and device for quantitative image correction for optical coherence tomography
JP2009066325A (ja) * 2007-09-14 2009-04-02 Tomey Corporation 前眼部断面画像解析システム及び前眼部断面画像撮影装置
JP4933413B2 (ja) * 2007-12-11 2012-05-16 株式会社トーメーコーポレーション 前眼部光干渉断層撮影装置及び前眼部光干渉断層撮影方法
JP5175781B2 (ja) * 2009-03-23 2013-04-03 株式会社トーメーコーポレーション 前眼部3次元画像処理装置
JP5410954B2 (ja) * 2009-12-29 2014-02-05 株式会社ニデック 眼軸長測定装置
JP5421136B2 (ja) * 2010-01-21 2014-02-19 株式会社ニデック 前眼部断面画像解析方法及び前眼部断面画像解析装置
JP2014500096A (ja) * 2010-12-03 2014-01-09 オプトビュー,インコーポレーテッド 光干渉断層法を用いたスキャン及び処理
US8687866B2 (en) * 2011-03-04 2014-04-01 Nanyang Technological University Methods and systems for processing images of the anterior chamber angle of an eye
US9095281B2 (en) * 2012-02-10 2015-08-04 Carl Zeiss Meditec, Inc. Segmentation and enhanced visualization techniques for full-range fourier domain optical coherence tomography

Also Published As

Publication number Publication date
JP2015066084A (ja) 2015-04-13

Similar Documents

Publication Publication Date Title
JP6367534B2 (ja) 前眼部3次元画像処理装置、プログラムおよび前眼部3次元画像処理方法
EP2865324B1 (en) Anterior segment three-dimensional image processing apparatus, and anterior segment three-dimensional image processing method
JP6217185B2 (ja) 眼科撮影装置及び眼科画像処理プログラム
JP6301621B2 (ja) 2次元断層画像処理装置、プログラムおよび2次元断層画像処理方法
JP5836634B2 (ja) 画像処理装置及び方法
EP2822448B1 (en) Apparatus for optical coherence tomography of an eye and method for optical coherence tomography of an eye
US10165939B2 (en) Ophthalmologic apparatus and ophthalmologic apparatus control method
US20160066779A1 (en) Image processing apparatus and control method of image processing apparatus
JP7104516B2 (ja) 断層画像撮影装置
JP2022027879A (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JP7164679B2 (ja) 眼科装置、及びその制御方法
WO2020050308A1 (ja) 画像処理装置、画像処理方法及びプログラム
JP2019042376A (ja) 画像処理装置、画像処理方法およびプログラム
JP6739725B2 (ja) 前眼部3次元画像処理装置および前眼部3次元画像処理方法
JP7005382B2 (ja) 情報処理装置、情報処理方法およびプログラム
JP7297952B2 (ja) 情報処理装置、情報処理方法およびプログラム
JP7130989B2 (ja) 眼科画像処理装置、および眼科画像処理プログラム
JP2018122137A (ja) 前眼部3次元画像処理装置、プログラムおよび前眼部3次元画像処理方法
JP6936676B2 (ja) 眼科装置
JP6942627B2 (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JP7158860B2 (ja) 画像処理装置、画像処理方法及びプログラム
JP6954831B2 (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JP6942626B2 (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JPWO2015129718A1 (ja) 画像処理装置、画像処理方法、及び画像処理プログラム
JP2023029593A (ja) 眼科撮影装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160810

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170412

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170425

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170616

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170912

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171010

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180206

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180301

R150 Certificate of patent or registration of utility model

Ref document number: 6301621

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250