JP4250245B2 - Fundus examination device - Google Patents
Fundus examination device Download PDFInfo
- Publication number
- JP4250245B2 JP4250245B2 JP04546699A JP4546699A JP4250245B2 JP 4250245 B2 JP4250245 B2 JP 4250245B2 JP 04546699 A JP04546699 A JP 04546699A JP 4546699 A JP4546699 A JP 4546699A JP 4250245 B2 JP4250245 B2 JP 4250245B2
- Authority
- JP
- Japan
- Prior art keywords
- blood vessel
- filter
- image
- information
- calculating
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Eye Examination Apparatus (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、眼底上の血管を検査する眼底検査装置に関するものである。
【0002】
【従来の技術】
従来から眼底血管の血管径を測定する方法は、眼底写真を顕微鏡観察して径を測定する方法と、光電素子を用いて眼底像を電気信号に変換して血管径を測定する方法が提案されている。前者はリアルタイムで測定できず、測定に要する時間が非常にかかる等の問題がある。また、後者は血管径算出のための信号処理に関しての詳しい開示は、鈴木らの方法:日眼会誌98−1pp92〜97が知られている。
【0003】
更に、He−Neレーザーを測定用光源として、被検眼の眼底血管に投射し、その反射光をCCDセンサで撮像し、この画像信号から眼底血管の径を算出しているものが特開平7−31596号公報で提案されている。
【0004】
【発明が解決しようとする課題】
しかし、光電素子から出力される画像信号をそのまま利用して血管径を算出した場合に、ノイズ、血管周辺の眼底の模様、血管壁の正反射成分等の影響を受けて、血管径の値を正確に算出すことは難しい。また、光電素子から出力される画像信号に対して、移動平均等のフィルタ処理を行い、画像信号に含まれる不要成分を除去した後の画像信号を用いて血管径を算出した場合に、測定対象となる血管の太さにより、血管壁の正反射成分や、血管像信号の特徴点の周波数成分が異なる。
【0005】
例えば、コントラストの強い太い血管を測定対象に選んだ場合に、移動平均点数を小さくすると、その血管や周辺組織の模様、ノイズ等の影響により血管像信号に含まれる特徴点が必要以上に存在することとなり、実際の血管径よりも細く算出されてしまう場合がある。また、細い血管に対して移動平均点数を大きくすると、血管像信号に含まれる必要な特徴点の情報が失われて、実際の血管径よりも太く算出される場合がある。
【0006】
本発明の目的は、上述の問題点を解消し、血管像の特徴点を検出し、正確な血管径を測定し得る眼底検査装置を提供することにある。
【0007】
【課題を解決するための手段】
上記目的を達成するための本発明に係る眼科検査装置は、ターゲットとなる眼底血管を含む領域を照明する照明手段と、前記領域の像を撮像して画像信号を出力する撮像手段と、該撮像手段からの画像信号の不要成分を除去する特性変更が可能なフィルタ手段と、前記画像信号を解析する解析手段と、該解析手段の解析結果を基に前記フィルタ手段のフィルタ特性を変更するフィルタ特性変更手段と、該フィルタ特性変更手段によって特性を変更した前記フィルタ手段の出力信号を基に前記ターゲットとなる眼底血管の情報を算出する血管情報算出手段とを有し、前記解析手段は前記画像信号の複数の極大点、極小点から血管情報に関する特徴点を算出する特徴点算出手段を有し、前記フィルタ特性変更手段は該特徴点算出手段によって検出した特徴点を基に前記フィルタ手段のフィルタ特性を変更することを特徴とする。
【0008】
また、本発明に係る眼底検査装置は、ターゲットとなる眼底血管を合む領域を照明する照明手段と、前記領域の像を撮像して画像信号を出力する撮像手段と、該撮像手段からの画像信号の不要成分を除去する特性変更が可能なフィルタ手段と、前記フィルタ手段のフィルタ特性を変更するフィルタ特性変更手段とを有し、該フィルタ手段の出力信号を基に前記ターゲットとなる眼底血管の情報を算出する第1のステップと、該第1のステップを複数の前記画像信号に対して行う第2のステップと、該第2のステップを実行する毎に前記フィルタ特性変更手段によって前記フィルタ手段の特性を変更し、前記第2のステップを複数回行う第3のステップと、前記第3のステップの出力結果の内の最もばらつきの少ないフィルタ特性を採用することを特徴とする。
【0010】
本発明に係る眼底検査装置は、ターゲットとなる眼底血管を含む領域を照明する照明手段と、前記領域の像を撮像して画像信号を出力する撮像手段と、該撮像手段からの画像信号の不要成分を除去する特性変更が可能なフィルタ手段と、前記照明手段の照明位置を記憶する照明位置記憶手段と、該照明位置記憶手段の記憶された同一部位の前回値の情報を基に前記フィルタ手段のフィルタ特性を変更するフィルタ特性変更手段と、該フィルタ特性変更手段によって特性を変更された前記フィルタ手段の出力信号を基に前記ターゲットとなる眼底血管の情報を算出する血管情報算出手段とを有することを特徴とする。
【0011】
【発明の実施の形態】
本発明を図示の実施例に基づいて詳細に説明する。
図1は本発明を眼底血流計に応用した実施例の構成図であり、白色光を発するタングステンランプ等から成る観察用光源1から被検眼Eと対向する対物レンズ2に至る照明光路上には、ほぼ被検眼Eの眼底と光学的に共役な位置に配され光路に沿って移動自在な固視標表示用素子である透過型液晶板3、リレーレンズ4、孔あきミラー5、黄色域の波長光を透過し他の光束を殆ど反射するバンドパスミラー6が順次に配列されている。孔あきミラー5の背後には眼底観察光学系が構成されており、光路に沿って移動自在な結像レンズ7、接眼レンズ8が順次に配列され、検者眼eに至っている。
【0012】
バンドパスミラー6の反射方向の光路上には、イメージローテータ9、紙面に垂直な回転軸を有し両面研磨されたガルバノメトリックミラー10が配置され、ガルバノメトリックミラー10の下側反射面10aの反射方向には光路に沿って移動自在なリレーレンズ11が配置され、上側反射面10bの反射方向にはレンズ12、光路に沿って移動自在なフォーカスユニット13が配置されている。
【0013】
なお、ガルバノメトリックミラー10は前述の回転軸の下方に切欠部を有しており、またレンズ12の前側焦点面は被検眼Eの瞳孔と共役関係にあり、この焦点面にガルバノメトリックミラー10が配置されている。また、ガルバノメトリックミラー10の後方には、レンズ14及び凹面ミラー15が配され、ガルバノメトリックミラー10の下側反射面10aで反射されず切欠部を通過する光束を、ガルバノメトリックミラー10の上側反射面10bに導くリレー光学系が構成されている。
【0014】
フォーカスユニット13においては、レンズ12と同一光路上にダイクロイックミラー16が配列され、このダイクロイックミラー16の反射方向の光路上には矩形状の絞りを有するマスク板17、ミラー18が配置され、ダイクロイックミラー16の透過方向の光路上にはレンズ19が配置されており、このフォーカスユニット13は一体的に移動ができるようになつている。
【0015】
また、レンズ19の入射方向の光路上には、コリメートされたコヒーレントな例えば赤色光を発するレーザーダイオードなどの測定用光源20が配列されている。更に、ミラー18の入射方向の光路上には、高輝度の他の光源と異なる例えば緑色光を発するヘリウムネオンレーザーなどのトラッキング用光源21が配列されている。
【0016】
ガルバノメトリックミラー10の下側反射面10aの反射方向の光路上には、光路に沿って移動自在なリレーレンズ11、ダイクロイックミラー22、拡大レンズ23、イメージインテンシファイヤ24、一次元CCD25が順次に配列され、血管検出系が構成されている。また、ダイクロイックミラー22の反射方向にはフォトマルチプライヤ26、27が配置され、測定用受光光学系が構成されている。なお、図示の都合上、全ての光路を同一平面上に示したが、ダイクロイックミラー22の反射方向などは紙面に直交している。
【0017】
更に、装置全体を制御するためのシステム制御部30が設けられ、このシステム制御部30には検者が操作する入力手段31、測定結果を表示する表示手段32、フォトマルチプライヤ26、27の出力、イメージインテンシファイヤ24及び一次元CCD25の出力がそれぞれ接続されており、システム制御部30の出力はガルバノメトリックミラー10を制御するガルバノメトリックミラー制御回路33に接続されている。
【0018】
また、前述の透過型液晶板3、結像レンズ7、フォーカスユニット13及びリレーレンズ11は、図示しないフォーカシングノブを操作することにより、被検眼Eの眼底Eaと透過型液晶板3、検者眼eの眼底、マスク板17及びイメージインテンシファイヤ24の受光面とが常に光学的に共役になるように、共に連動して光軸方向に移動するようになっている。
【0019】
図2はシステム制御部30の構成図を示し、一次元CCD25の出力信号はサンプルホールド回路41に入力されサンプルホールドされた後に、増幅器42に入力され適当なゲインで増幅される。増幅器42の出力信号はA/D変換器43に入力され、ディジタル信号に変換されて、血管像を含む画像データとしてメモリ44に格納される。
【0020】
演算部45はメモリ44から血管像を含む画像データを読み出し、この信号に対して不要な周波数成分をカットしノイズ除去処理するディジタルフィルタリング処理、極大・極小点を特徴点として算出、血管径の算出、ディジタルフィルタのフィルタ特性を変更等を行う。
【0021】
ここで、各光源から発せられた光束が被検眼Eに導かれ、その反射散乱光が検者又はフォトマルチプライヤ26、27、イメージインテンシファイヤ24の受光面に導かれる。観察用光源1から発した白色光は、透過型液晶板3を背後から照明し、リレーレンズ4を通って孔あきミラー5で反射され、黄色域の波長光のみがバンドパスミラー6を透過し、対物レンズ2を通り、被検眼Eの瞳孔Ep上で眼底照明光光束像Iとして一旦結像した後に、眼底Eaをほぼ一様に照明する。このとき、透過型液晶板3には図示しない固視標Fが表示されており、照明光により被検眼Eの眼底Eaに投影され、視標像F’として被検眼Eに呈示される。
【0022】
眼底Eaからの反射光は同じ光路を戻り、瞳孔上から眼底観察光束Oとして取り出され、孔あきミラー5の中心の開口部、結像レンズ7を通り、検者眼eにより接眼レンズ8を介して図3に示す眼底像Ea’が観察可能となる。検者はこの眼底像Ea’を観察しながら、接眼レンズ8より装置のアライメントを行う。
【0023】
レーザーダイオードなどの測定用光源20を発したコリメートされた測定光は、レンズ19を通過しダイクロイックミラー16を透過する。一方、ヘリウムネオンレーザーなどのトラッキング用光源21から発したトラッキング光は、ミラー18で反射された後に、マスク板17で所望の形状に整形された後に、ダイクロイックミラー16に反射されて、上述の測定光と重畳される。このとき、測定光はレンズ19によりマスク板17の開口部中心と共役な位置にスポット状に結像されている。更に、測定光とトラッキング照射光はレンズ12を通り、ガルバノメトリックミラー10の上側反射面10bで反射され、レンズ14を一旦通った後に凹面ミラー15で反射され、再びレンズ14を通りガルバノメトリックミラー10の方に戻される。
【0024】
ここで、ガルバノメトリックミラー10は被検眼瞳の共役な位置に配されており、凹面ミラー15とレンズ14は光軸上に同心に配置されかつ共働して、ガルバノメトリックミラー10を−1倍で結像するリレー系の機能が与えられている。
【0025】
従って、ガルバノメトリックミラー10の上側反射面10bで反射された光束は、ガルバノメトリックミラー10の切欠部に戻されることになり、ガルバノメトリックミラー10で反射されることなくイメージローテータ9に向かうことになる。イメージローテータ9を経て、バンドパスミラー6により対物レンズ2へ偏向された両光束は、対物レンズ2を介して被検眼Eの眼底Eaを照射する。
【0026】
このように、測定光とトラッキング照射光はガルバノメトリックミラー10の上側反射面10b内で反射されて、再び戻されるように対物レンズ2の光軸から偏心した状態でガルバノメトリックミラー10に入射され、図4に示すように瞳孔Ep上でスポット像P2又はP2’として結像した後に眼底Eaを点状に照射する。
【0027】
眼底Eaからの測定光とトラッキング照射光の散乱反射光は再び対物レンズ2で集光され、バンドパスミラー6で殆どの光束が反射されてイメージローテータ9を通り、ガルバノメトリックミラー10の下側反射面10aで反射され、リレーレンズ11を通り、ダイクロイックミラー22において測定光とトラッキング光とが分離される。
【0028】
トラッキング光はダイクロイックミラー22を透過し、拡大レンズ12により眼底観察光学系による眼底像Ea’よりも拡大された血管像Ev’として、イメージインテンシファイヤ24の光電面に結像し、増幅された後に一次元CCD25により撮像される。そして、一次元CCD25で撮像された血管像Ev’に基づいて、システム制御部30において血管像Ev’の移動量を表すデータが作成され、ガルバノメトリックミラー制御回路33に血管像Ev’と移動量が出力される。そして、ガルバノメトリックミラー制御回路33がこの移動量を補償するようにガルバノメトリックミラー10を駆動することにより、被測定部の血管のトラッキングが行える。
【0029】
また、測定光はダイクロイックミラー22により反射され、フォトマルチプライヤ26、27に受光される。フォトマルチプライヤ26、27の出力はそれぞれシステム制御部30に出力され、この受光信号は従来例と同様に周波数解析されて眼底Eaの血流速度が求められる。
【0030】
一方、眼底Eaからの測定光とトラッキング照射光の散乱反射光は再び対物レンズ2で集光され、バンドパスミラー6を透過した一部の光束は観察用光源1から発した光束の被検眼Eの眼底Eaからの反射散乱光と同様の光路をたどって、検者眼eに達し、観察眼底像Ea’と共にトラッキング指標像、測定光像として検者が観察できるようになっている。
【0031】
先ず、検者は図示しない操作桿を操作して、被検眼Eの光軸と対物レンズ2の光軸が一致するように位置合わせを行う。次に、眼底像Ea’を観察しながら前述のフォーカスノブを操作して、被検眼Eの眼底Eaにフォーカスを合わせる。すると、前述したように透過型液晶板3の固視標Fと眼底Eaが光学的に共役になり被検眼Eに呈示され、被検者が固視標像F’を固視すると検者は図3に示したような眼底像Ea’を観察できる。そして、検者は第1の被測定部位が観察視野の略中央付近にくるように入力手段31を操作して固視標Fを動かし、被検眼Eを誘導する。
【0032】
次に、入力手段31を操作してトラッキング光を眼底Eaに照射し、更にトラッキング指標像と第1の被測定血管に垂直になるように図示しないローテータ操作ノブを操作し、第1の被測定血管上に測定光が照射されるようにガルバノメトリックミラー10の角度を制御する。トラッキング光で照射された血管Evは、前述のように血管像Ev’としてイメージインテンシファイヤ24の光電面に結像し、増幅された後に一次元CCD25上に撮像され、血管像を含む画像信号として出力される。
【0033】
検者は測定部位を決定した後に、再び入力手段31を操作してトラッキングの開始を入力する。トラッキングが安定したところで、検者は再び入力手段31を操作して測定の開始を入力する。一次元CCD25から出力される画像信号は、サンプルホールド回路41でサンプルホールドされた後に、増幅器42に入力され適当なゲインで増幅される。増幅器42の出力信号はA/D変換器43に入力されてディジタル信号に変換され、血管像を含む画像データとしてメモリ44に格納される。
【0034】
図5のフローチャート図において、システム制御部30の演算部45はステップS2で血管径の仮計算を行う設定とし、ステップS3でメモリ44に記憶された画像データを読み出す。ステップS4で不要な周波数成分を遮断するノイズ成分除去処理を行うフィルタ処理が行われる。このときのフィルタ特性は、ステップS1で設定された仮計算時の初期設定値に設定されている。本実施例では、フィルタ特性は移動平均計算を行うように設定してあり、血管径の仮計算時には測定可能な全ての血管径に対して高周波成分のノイズ成分だけを除去できるような移動平均点数に設定されている。ステップS4でのフィルタ処理後に、ステップS5で入力データの微分処理を行い、複数の極大、極小点を算出する。更にステップS6で、ステップS5で得られた複数の極大、極小点の中から、血管情報に関する特徴点を抽出し、ステップS7ではステップS6で得られた血管惰報に関する特徴点を用いて、仮の血管径を算出する。
【0035】
血管径の仮計算が終了すると、ステップS7で算出された仮の血管径値を基に、フィルタ特性の変更が必要かどうかを判断し、必要な場合はステップS9でフィルタ特性の変更を行って本計算に移行し、再びステップS3〜S7までの処理を行い、最終的な血管径を算出する。必要でない場合には、仮計算の血管径の値を最終的な血管径値とする。なお、本実施例では仮の血管径値を基にフィルタの移動平均点数を決定・変更することで、フィルタの特性を変更している。
【0036】
ステップS6で行う複数の極大、極小点の中からの血管情報に関する特徴点の抽出については、図6の信号波形において、上側反射面の暗い部分つまりMax1〜Max2の間が、ほぼ血管の像を表している。また、この血管像信号はコントラストの強い太い血管を測定対象に選んだ場合に得られたもので、その血管中心には血管壁の正反射が観測されていて明るくなっている部分がある。更に、血管像を含む映像信号波形を微分処理して算出した極大値座標を「△」で、極小値座標を「×」で示している。
【0037】
ステップS61で、極小値の中での極小値Min1を検出する。ステップS62で、極大値の中で極大値Max0を検出する。ステップS63でMin1の次に小さい極小値をMin2を検出する。ただし、極小値Min2の検出には振幅方向に一定の基準値を設ける必要があり、これを設けないと不必要なノイズ成分を誤検出してしまう場合がある。本実施例では、Min1−Max0の幅の1/2(=Ref)よりも小さいことを条件としている。ステップ64で極小値Min1の左側でMin1に最も近い極大値Max1を検出する。ステップS65で極小値Min2の右側でMin2に最も近いMin2に最も近い極大値Max2を検出する。ステップS66で極小値Min1とMin2の間で最も大きい極大値Max2を検出し、ステップS67でMax1、Max2、Max3、Min1、Min2の5点を血管像の特徴点として検出する。
【0038】
なお、コントラストの小さい細い血管を測定対象に選んだ場合には、血管中心には血管壁の正反射成分が観測されない場合があり、この場合には、Max3、Min2は検出されず、Max1、Max2、Min1の3点が血管像の特徴点として検出される。
【0039】
図7(a)は太い血管(血管径=約200μm)を測定対象に選んだ場合のフィルタ処理前の画像信号を示す。本実施例では、血管径の仮計算時には、移動平均点数は10画素に設定されていて、演算部45は図5に示すステップS1〜S6までの処理を行い、仮の血管径=約200μmを算出する。
【0040】
表1は本実施例で設定している仮計算で得られた血管径の値と、ステップS9で設定するフィルタ特性である移動平均の画素数の関係を示す。
【0041】
【0042】
この表1から仮計算で得られた血管径が200μmの場合の移動平均の画素数は10画素であるため、ステップS9のフィルタ特性の再設定は行わずに、血管径は200μmを最終的な血管径の値として算出する。図7(c)に移動平均の画素数は10画素でフィルタ処理を行った後の画像信号を示すが、Max1、Max2、Max3、Min1、Min2の5個の極大・極小点を用いて血管径を算出する。
【0043】
図7(b)は移動平均点数を4画素に設定した場合の例で、血管信号に対して高い周波数成分のノイズが取りきれていないため、図7(c)の場合と比較すると血管像の特徴点Max1、Max2、Max3、Min1、Min2の内、Min2、Max3の検出位置がずれていて、これらの特徴点を用いて血管径を算出すると、血管径は100μmとなり、実際の値よりも細い値を算出してしまう。
【0044】
図7(d)は細い血管(血管径は約100μm)を測定対象に選んだ場合のフィルタ処理前の画像信号を示す。演算部45は図5に示すステップS1〜S6までの処理を行い、仮の血管径を約140μmとして算出する。このときのフィルタ処理後の出力信号を図7(f)に示す。仮計算では、図7(f)に示すMax1、Max2、Min1の3個の極大・極小点が血管像情報の特徴点として検出され、これらを用いて血管径を算出する。
【0045】
表1から、仮計算で得られた血管径が140μmの場合の移動平均の画素数は6画素であるため、ステップS9でフィルタ処理の移動平均点数を6画素に設定して、再度ステップS2からS6までの処理を行い、血管径を100μmとして算出する。このときのフィルタ処理後の画像信号を図7(e)に示すが、Max1、Max2、Max3、Min1、Min2の5個の極大−極小点が血管像情報の特徴点として検出され、これらを用いて血管径を算出する。
【0046】
このように、仮計算で得られた血管径によって本計算に用いる血管像波形に対する移動平均点数を可変とすることにより、血管の情報を損なうことなく、また不要な周波数成分を拾うことなく、血管像の特徴点を正確に検出することができる。
【0047】
本実施例における血管径の仮計算時に、ステップS7ではステップS6で得られた血管像情報の特徴点から、仮の血管径の値を算出し、この値を基にステップS9でフィルタ特性の変更を行っているが、ステップS7では血管径の仮計算時に血管径を算出する代りに、ステップS6で得られた特徴点の内、両端の極大点間距離情報、Max1−Max2間の距離情報を算出する。この距離情報は血管径の値に関連するため、この距離情報を基にステップS9でフィルタ特性を変更し、本計算時に血管径を算出してもよい。
【0048】
図8は第2の実施例の演算部45の処理フローチャート図であり、本実施例では60mSの間に15個の血管像を含む画像信号をメモリ44に取り込み、これらの信号に対して5種類のフィルタ特性のフィルタ処理を行い、血管径を算出している。
【0049】
ステップS10でフィルタ特性を変更する回数カウンタを初期化し、ステップS11でフィルタ特性の初期設定を行う。ステップS12で画像信号を取り込む回数カウンタを初期化する。ステップS13でメモリ44に記憶されたi個目(0<i<14)の画像データは演算部45によって読み出され、ステップS14で不要な周波数成分を連断するノイズ成分除去処理を行うフィルタ処理が行われる。このときのフィルタ特性は、ステップS11で設定された初期設定に設定されている。
【0050】
なお、本実施例においてフィルタは移動平均計算を行うように設定してある。ステップS14でのフィルタ処理後に、ステップS15で入力データの微分処理を行い、複数の極大、極小点を算出する。更に、ステップS16でステップS15で得られた複数の極大、極小点の中から、血管情報に関する特徴点を抽出し、ステップS17ではステップS16で得られた血管情報に関する特徴点を用いて血管径を算出する。
【0051】
ここで、ステップS16での複数の極大・極小点から血管情報の特徴点を抽出する方法の詳細は第1の実施例と同様であり、図5のステップS61〜S67までの処理を行う。15個分の画像信号について血管径の算出が終了するまでステップS13〜S17までの処理を繰り返し、15個分の画像信号について血管径の算出が終了すると、ステップS19で15個分の血管径の平均値・標準偏差を算出する。
【0052】
更に、フィルタ特性を5回変更するまで、ステップS11〜S19までの処理を繰り返す。本実施例において、ステップS23でのフィルタ特性の変更は、移動平均処理を行う画素数を変更する処理を行っていて、2、4、6、8、10画素の5種類の画素数で移動平均を行うように設定されている。ここで、フィルタ特性の変更が5回に達すると、ステップS23ではステップS19で算出した5種類のフィルタ特性での血管径の値の内、標準偏差が最も小さいものを最終的な血管径の値として選択する。
【0053】
このように、複数の血管像を含む画像信号を取り込み、血管径の値を算出し、その平均値を血管径の値とする。更に、血管像を含む画像信号に対して、フィルタ特性を変えた複数のフィルタ処理を行う、本実施例においては、異なる移動平均点数で複数の血管径の値を算出し、この中で最もばらつきの小さいものを最終的な血管径の値とすることにより、血管の情報を損なうことなく、また不要な周波数成分を拾うことなく血管像の特徴点を正確に検出することができ、より正確な血管径の値を算出することが可能となる。
【0054】
本実施例では、血管径の仮測定で得られた値を基に、フィルタ処理の特性を変更していたが、検者が予め測定対象血管の概略値を図示しない入力手段から入力し、この入力情報を基にフィルタの特性を決定してもよい。
【0055】
また、演算部45は血管径の測定部位とフィルタ特性をメモリ44に記憶しておいて、同一部位の血管径を測定する場合に、第1、第2の実施例で求めた前回のフィルタ特性の設定値を用いて血管径を算出してもよい。
【0056】
測定部位は検者が図示しない入力部から入力してもよいし、本装置がイメージローテータ9の角度、ガルバノメトリックミラー10の角度、及び透過型液晶板3に投影される内部固視標の位置等を記憶することもできる。また、眼底画像をテレビカメラで取り込み、画像処理により測定部位を断定してもよい。
【0057】
以上の説明では、眼底血流計、特にその血管径算出に応用した実施例について詳細に説明したが、血管径の計測に関する本発明の構成・効果は血管位置の特定においても同様の効果を有する。その違いは単に血管の走行方向に垂直な情報として血管位置を算出するのか、血管径を算出するかの違いであるため、ここでは詳細には言及しないこととする。
【0058】
【発明の効果】
以上説明したように本発明に係る眼科検査装置は、取り込んだ画像信号を解析した結果に基づいてフィルタの特性を可変とし、このフィルタを通過した画像信号を用いて血管径の値を算出するため、血管の太さ、血管の模様、血管周辺の眼底の模様、血管壁の正反射成分、ノイズの影響を受け難く、画像信号に含まれる血管像の特徴点を正確に検出することができるため、より正確な血管径計測が可能になる。
【図面の簡単な説明】
【図1】第1の実施例の構成図である。
【図2】システム制御部の構成図である。
【図3】観察眼底像の説明図である。
【図4】被検眼Eの瞳面上の各光束の様子の説明図である。
【図5】第1の実施例における演算部の処理を表すフローチャート図である。
【図6】血管像信号波形の説明図である。
【図7】血管像信号波形の説明図である。
【図8】第2の実施例における演算部の処理を表すフローチャート図である。
【符号の説明】
1 観察用光源
2 対物レンズ
5 孔あきミラー
9 イメージローテータ
10 ガルバノメトリックミラー
15 凹面ミラー
20 測定用光源
23 トラッキング用光源
24 イメージインテンシファイア
25 一次元CCD
26、27 フォトマルチプライヤ
30 システム制御部
E 被検眼
Ea’ 被検眼眼底
Ev’ 血管[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a fundus examination apparatus for examining blood vessels on the fundus.
[0002]
[Prior art]
Conventional methods for measuring the diameter of the fundus blood vessel have been proposed: a method of measuring the diameter by observing a fundus photograph with a microscope, and a method of measuring the diameter of a blood vessel by converting a fundus image into an electrical signal using a photoelectric element. ing. The former cannot be measured in real time, and there is a problem that it takes a very long time for measurement. As for the latter, a detailed disclosure regarding signal processing for calculating a blood vessel diameter is known from Suzuki et al .: Nikkankai 98-1pp92-97.
[0003]
Furthermore, the He-Ne laser is used as a measurement light source, projected onto the fundus blood vessel of the eye to be examined, the reflected light is imaged by a CCD sensor, and the diameter of the fundus blood vessel is calculated from this image signal. This is proposed in Japanese Patent No. 31596.
[0004]
[Problems to be solved by the invention]
However, when the blood vessel diameter is calculated using the image signal output from the photoelectric element as it is, the value of the blood vessel diameter is affected by noise, the fundus pattern around the blood vessel, the specular reflection component of the blood vessel wall, etc. It is difficult to calculate accurately. In addition, when the blood vessel diameter is calculated using the image signal after removing unnecessary components included in the image signal by performing filter processing such as moving average on the image signal output from the photoelectric element, the measurement target Depending on the thickness of the blood vessel, the regular reflection component of the blood vessel wall and the frequency component of the characteristic point of the blood vessel image signal are different.
[0005]
For example, when a thick blood vessel with high contrast is selected as the measurement target, if the moving average score is reduced, there are more feature points included in the blood vessel image signal due to the influence of the blood vessel, surrounding tissue pattern, noise, etc. In other words, it may be calculated to be narrower than the actual blood vessel diameter. Further, if the moving average score is increased for a thin blood vessel, information on necessary feature points included in the blood vessel image signal may be lost, and may be calculated to be thicker than the actual blood vessel diameter.
[0006]
An object of the present invention is to provide a fundus examination apparatus that can solve the above-described problems, detect a feature point of a blood vessel image, and measure an accurate blood vessel diameter.
[0007]
[Means for Solving the Problems]
To achieve the above object, an ophthalmic examination apparatus according to the present invention includes an illumination unit that illuminates a region including a target fundus blood vessel, an imaging unit that captures an image of the region and outputs an image signal, and the imaging unit. Filter means capable of changing characteristics for removing unnecessary components of the image signal from the means, analysis means for analyzing the image signal, and filter characteristics for changing the filter characteristics of the filter means based on the analysis result of the analysis means And a blood vessel information calculation unit that calculates information about the target fundus blood vessel based on an output signal of the filter unit whose characteristics have been changed by the filter characteristic changing unit. The analysis means includes feature point calculation means for calculating a feature point related to blood vessel information from a plurality of maximum points and minimum points of the image signal, and the filter characteristic changing means is a feature point detected by the feature point calculation means. The filter characteristics of the filter means are changed based on It is characterized by that.
[0008]
The fundus examination apparatus according to the present invention includes an illuminating unit that illuminates a region that joins a target fundus blood vessel, an imaging unit that captures an image of the region and outputs an image signal, and an image from the imaging unit. Filter means capable of changing characteristics for removing unnecessary components of the signal, and filter characteristic changing means for changing the filter characteristics of the filter means, and based on an output signal of the filter means, the target fundus blood vessel A first step of calculating information; a second step of performing the first step on a plurality of the image signals; and the filter means changing means by the filter characteristic changing unit each time the second step is executed. The third step of performing the second step a plurality of times and the filter characteristic with the least variation among the output results of the third step are adopted. The features.
[0010]
The fundus examination apparatus according to the present invention includes an illumination unit that illuminates a region including a target fundus blood vessel, an imaging unit that captures an image of the region and outputs an image signal, and an image signal from the imaging unit is unnecessary. Filter means capable of changing characteristics for removing components, illumination position storage means for storing the illumination position of the illumination means, and the filter means based on the previous value information of the same part stored in the illumination position storage means Filter characteristic changing means for changing the filter characteristic of the blood vessel, and blood vessel information calculating means for calculating information of the target fundus blood vessel based on the output signal of the filter means whose characteristic has been changed by the filter characteristic changing means It is characterized by that.
[0011]
DETAILED DESCRIPTION OF THE INVENTION
The present invention will be described in detail based on the embodiments shown in the drawings.
FIG. 1 is a block diagram of an embodiment in which the present invention is applied to a fundus blood flow meter. The illumination light path extends from an
[0012]
On the optical path in the reflection direction of the
[0013]
The
[0014]
In the
[0015]
Further, on the optical path in the incident direction of the
[0016]
On the optical path in the reflection direction of the lower
[0017]
Further, a
[0018]
The transmission
[0019]
FIG. 2 is a block diagram of the
[0020]
The
[0021]
Here, the luminous flux emitted from each light source is guided to the eye E, and the reflected and scattered light is guided to the light-receiving surface of the examiner or the
[0022]
Reflected light from the fundus oculi Ea returns along the same optical path, is taken out from the pupil as a fundus oculi observation light beam O, passes through an opening at the center of the
[0023]
The collimated measurement light emitted from the
[0024]
Here, the
[0025]
Accordingly, the light beam reflected by the
[0026]
Thus, the measurement light and the tracking irradiation light are reflected on the
[0027]
The scattered and reflected light of the measurement light and tracking irradiation light from the fundus oculi Ea is condensed again by the
[0028]
The tracking light passes through the
[0029]
Further, the measurement light is reflected by the
[0030]
On the other hand, the scattered and reflected light of the measurement light from the fundus oculi Ea and the tracking irradiation light is collected again by the
[0031]
First, the examiner operates an operating rod (not shown) to perform alignment so that the optical axis of the eye E and the optical axis of the
[0032]
Next, the input means 31 is operated to irradiate the fundus Ea with tracking light, and further, a rotator operation knob (not shown) is operated so as to be perpendicular to the tracking index image and the first blood vessel to be measured. The angle of the
[0033]
After determining the measurement site, the examiner operates the
[0034]
In the flowchart of FIG. 5, the
[0035]
When the temporary calculation of the blood vessel diameter is completed, it is determined whether or not the filter characteristic needs to be changed based on the temporary blood vessel diameter value calculated in step S7, and if necessary, the filter characteristic is changed in step S9. Shifting to this calculation, the processing from step S3 to S7 is performed again, and the final blood vessel diameter is calculated. If not necessary, the provisionally calculated blood vessel diameter value is used as the final blood vessel diameter value. In this embodiment, the filter characteristics are changed by determining and changing the moving average score of the filter based on the provisional blood vessel diameter value.
[0036]
With respect to extraction of feature points relating to blood vessel information from a plurality of maximum and minimum points performed in step S6, in the signal waveform of FIG. 6, a dark portion of the upper reflection surface, that is, between Max1 and Max2, substantially represents a blood vessel image. Represents. Further, this blood vessel image signal is obtained when a thick blood vessel with high contrast is selected as a measurement target, and a regular reflection of the blood vessel wall is observed at the center of the blood vessel, and there is a bright portion. Further, the maximum value coordinate calculated by differentiating the video signal waveform including the blood vessel image is indicated by “Δ”, and the minimum value coordinate is indicated by “x”.
[0037]
In step S61, the minimum value Min1 in the minimum value is detected. In step S62, the maximum value Max0 is detected among the maximum values. In step S63, Min2 is detected as a minimum value next to Min1. However, in order to detect the minimum value Min2, it is necessary to provide a constant reference value in the amplitude direction. If this is not provided, an unnecessary noise component may be erroneously detected. In the present embodiment, the condition is that it is smaller than ½ (= Ref) of the width of Min1−Max0. In step 64, the maximum value Max1 closest to Min1 is detected on the left side of the minimum value Min1. In step S65, the maximum value Max2 closest to Min2 closest to Min2 is detected on the right side of the minimum value Min2. In step S66, the maximum maximum value Max2 between the minimum values Min1 and Min2 is detected, and in step S67, five points of Max1, Max2, Max3, Min1, and Min2 are detected as feature points of the blood vessel image.
[0038]
If a thin blood vessel with low contrast is selected as the measurement target, the regular reflection component of the blood vessel wall may not be observed at the center of the blood vessel. In this case, Max3 and Min2 are not detected, and Max1 and Max2 , Min1 are detected as feature points of the blood vessel image.
[0039]
FIG. 7A shows an image signal before filtering when a thick blood vessel (blood vessel diameter = about 200 μm) is selected as a measurement target. In this embodiment, at the time of provisional calculation of the blood vessel diameter, the moving average score is set to 10 pixels, and the
[0040]
Table 1 shows the relationship between the value of the blood vessel diameter obtained by the provisional calculation set in this embodiment and the number of pixels of the moving average which is the filter characteristic set in step S9.
[0041]
[0042]
When the blood vessel diameter obtained by provisional calculation from Table 1 is 200 μm, the moving average number of pixels is 10 pixels. Therefore, without resetting the filter characteristics in step S9, the blood vessel diameter should be 200 μm at the end. Calculated as the value of the blood vessel diameter. FIG. 7 (c) shows the image signal after filtering with a moving average pixel count of 10 pixels. The maximum and minimum points of Max1, Max2, Max3, Min1, and Min2 are used to determine the vessel diameter. Is calculated.
[0043]
FIG. 7B shows an example in which the moving average score is set to 4 pixels. Since the high frequency component noise is not completely removed from the blood vessel signal, the blood vessel image is compared with the case of FIG. 7C. Among the feature points Max1, Max2, Max3, Min1, and Min2, the detection positions of Min2 and Max3 are misaligned. If the vessel diameter is calculated using these feature points, the vessel diameter is 100 μm, which is smaller than the actual value. The value is calculated.
[0044]
FIG. 7D shows an image signal before filter processing when a thin blood vessel (blood vessel diameter is about 100 μm) is selected as a measurement target. The
[0045]
From Table 1, since the number of moving average pixels when the blood vessel diameter obtained by provisional calculation is 140 μm is 6 pixels, the moving average score of the filtering process is set to 6 pixels in step S9, and again from step S2. The processing up to S6 is performed, and the blood vessel diameter is calculated as 100 μm. FIG. 7E shows the image signal after the filter processing at this time, and five maximum-minimum points of Max1, Max2, Max3, Min1, and Min2 are detected as feature points of the blood vessel image information and are used. To calculate the vessel diameter.
[0046]
In this way, by making the moving average score for the blood vessel image waveform used in this calculation variable according to the blood vessel diameter obtained in the preliminary calculation, the blood vessel can be obtained without damaging blood vessel information and without picking up unnecessary frequency components. It is possible to accurately detect image feature points.
[0047]
In the provisional calculation of the blood vessel diameter in the present embodiment, in step S7, a temporary blood vessel diameter value is calculated from the feature points of the blood vessel image information obtained in step S6, and the filter characteristics are changed in step S9 based on this value. In step S7, instead of calculating the blood vessel diameter at the time of temporary calculation of the blood vessel diameter, the distance information between the maximum points at both ends and the distance information between Max1 and Max2 among the feature points obtained in step S6 are obtained. calculate. Since this distance information is related to the value of the blood vessel diameter, the filter characteristic may be changed in step S9 based on this distance information, and the blood vessel diameter may be calculated during this calculation.
[0048]
FIG. 8 is a processing flowchart of the
[0049]
In step S10, the counter for changing the filter characteristics is initialized. In step S11, the filter characteristics are initialized. In step S12, the number counter for capturing image signals is initialized. The i-th (0 <i <14) image data stored in the
[0050]
In this embodiment, the filter is set to perform moving average calculation. After the filtering process in step S14, the input data is differentiated in step S15 to calculate a plurality of maximum and minimum points. Further, in step S16, feature points related to blood vessel information are extracted from the plurality of maximum and minimum points obtained in step S15. In step S17, the blood vessel diameter is determined using the feature points related to blood vessel information obtained in step S16. calculate.
[0051]
Here, the details of the method of extracting the feature points of the blood vessel information from the plurality of maximum / minimum points in step S16 are the same as in the first embodiment, and the processing from steps S61 to S67 in FIG. 5 is performed. Steps S13 to S17 are repeated until calculation of blood vessel diameters for 15 image signals is completed. When calculation of blood vessel diameters for 15 image signals is completed, the blood vessel diameters for 15 images are calculated in step S19. Calculate the mean and standard deviation.
[0052]
Further, the processing from step S11 to S19 is repeated until the filter characteristic is changed five times. In this embodiment, the filter characteristic is changed in step S23 by changing the number of pixels for which moving average processing is performed, and moving average is performed with five types of pixels of 2, 4, 6, 8, and 10 pixels. Is set to do. Here, when the change of the filter characteristic reaches five times, in step S23, among the blood vessel diameter values in the five types of filter characteristics calculated in step S19, the one having the smallest standard deviation is the final blood vessel diameter value. Choose as.
[0053]
In this manner, an image signal including a plurality of blood vessel images is taken in, a blood vessel diameter value is calculated, and the average value is set as the blood vessel diameter value. Further, a plurality of filter processes with different filter characteristics are performed on an image signal including a blood vessel image. In this embodiment, a plurality of blood vessel diameter values are calculated with different moving average points, and the most varied among them. By making the final value of the blood vessel diameter smaller, the characteristic points of the blood vessel image can be accurately detected without damaging blood vessel information and without picking up unnecessary frequency components. The value of the blood vessel diameter can be calculated.
[0054]
In this embodiment, the characteristics of the filter processing are changed based on the value obtained by the temporary measurement of the blood vessel diameter. However, the examiner inputs in advance an approximate value of the blood vessel to be measured from an input means (not shown). Filter characteristics may be determined based on input information.
[0055]
Further, the
[0056]
The measurement site may be input by the examiner from an input unit (not shown), or the angle of the
[0057]
In the above description, the fundus blood flow meter, in particular, the embodiment applied to the calculation of the blood vessel diameter has been described in detail. However, the configuration and effect of the present invention relating to the measurement of the blood vessel diameter has the same effect in specifying the blood vessel position. . The difference is simply whether the blood vessel position is calculated as information perpendicular to the blood vessel traveling direction or the blood vessel diameter is calculated, and therefore will not be described in detail here.
[0058]
【The invention's effect】
As described above, the ophthalmic examination apparatus according to the present invention makes the characteristics of the filter variable based on the result of analyzing the captured image signal, and calculates the value of the blood vessel diameter using the image signal that has passed through the filter. Because it is difficult to be affected by the thickness of the blood vessel, the blood vessel pattern, the fundus pattern around the blood vessel, the specular reflection component of the blood vessel wall, and noise, the feature points of the blood vessel image included in the image signal can be accurately detected. More accurate blood vessel diameter measurement becomes possible.
[Brief description of the drawings]
FIG. 1 is a configuration diagram of a first embodiment.
FIG. 2 is a configuration diagram of a system control unit.
FIG. 3 is an explanatory diagram of an observation fundus image.
FIG. 4 is an explanatory diagram of the state of each light beam on the pupil plane of the eye E to be examined.
FIG. 5 is a flowchart showing processing of a calculation unit in the first embodiment.
FIG. 6 is an explanatory diagram of a blood vessel image signal waveform.
FIG. 7 is an explanatory diagram of a blood vessel image signal waveform.
FIG. 8 is a flowchart showing processing of a calculation unit in the second embodiment.
[Explanation of symbols]
1 Light source for observation
2 Objective lens
5 Perforated mirror
9 Image rotator
10 Galvanometric mirror
15 Concave mirror
20 Light source for measurement
23 Light source for tracking
24 Image Intensifier
25 One-dimensional CCD
26, 27 Photomultiplier
30 System controller
E Eye to be examined
Ea 'eye fundus
Ev 'blood vessels
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP04546699A JP4250245B2 (en) | 1999-02-23 | 1999-02-23 | Fundus examination device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP04546699A JP4250245B2 (en) | 1999-02-23 | 1999-02-23 | Fundus examination device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000237141A JP2000237141A (en) | 2000-09-05 |
JP4250245B2 true JP4250245B2 (en) | 2009-04-08 |
Family
ID=12720164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP04546699A Expired - Fee Related JP4250245B2 (en) | 1999-02-23 | 1999-02-23 | Fundus examination device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4250245B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4616972B2 (en) * | 2000-08-09 | 2011-01-19 | キヤノン株式会社 | Fundus examination device |
JP2018023478A (en) * | 2016-08-09 | 2018-02-15 | 大日本印刷株式会社 | Fundus image processing device |
DE102018107625A1 (en) * | 2018-03-29 | 2019-10-02 | Imedos Systems GmbH | Apparatus and method for examining retinal vascular endothelial function |
-
1999
- 1999-02-23 JP JP04546699A patent/JP4250245B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2000237141A (en) | 2000-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6569104B2 (en) | Blood vessel detecting apparatus | |
US5894337A (en) | Eye fundus examining apparatus | |
US6535757B2 (en) | Ocular examination system | |
JP3647164B2 (en) | Ophthalmic measuring device | |
JP3679547B2 (en) | Fundus blood vessel observation device | |
JP2003019119A (en) | Eyeground rheometer | |
JPH08215150A (en) | Ophthalmological diagnosis device | |
JP2003019116A (en) | Opthalmologic measuring device | |
US6699198B2 (en) | Ocular-blood-flow meter | |
JP4250245B2 (en) | Fundus examination device | |
JP2000135200A (en) | Optometric apparatus | |
JP3762025B2 (en) | Ophthalmic examination equipment | |
JP3880209B2 (en) | Ophthalmic examination equipment | |
JP4536884B2 (en) | Ophthalmic examination equipment | |
JP3591952B2 (en) | Fundus examination device | |
JP2001112716A (en) | Ophthalmic examination apparatus | |
JP3610139B2 (en) | Fundus examination device | |
JP4838428B2 (en) | Ophthalmic equipment | |
JP2001275976A (en) | Fonduscope | |
JP4616972B2 (en) | Fundus examination device | |
JP3604836B2 (en) | Ophthalmic measurement device | |
JP3624071B2 (en) | Ophthalmic equipment | |
JP4035247B2 (en) | Fundus blood flow meter | |
JP3639658B2 (en) | Beam deflector for ophthalmic examination | |
JP2000037351A (en) | Eyeground examination device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060217 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080118 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081014 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081204 |
|
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: 20081224 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090119 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120123 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4250245 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130123 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140123 Year of fee payment: 5 |
|
LAPS | Cancellation because of no payment of annual fees |