JP4478286B2 - 眼底血流計及び眼底血流速度測定方法 - Google Patents
眼底血流計及び眼底血流速度測定方法 Download PDFInfo
- Publication number
- JP4478286B2 JP4478286B2 JP2000102240A JP2000102240A JP4478286B2 JP 4478286 B2 JP4478286 B2 JP 4478286B2 JP 2000102240 A JP2000102240 A JP 2000102240A JP 2000102240 A JP2000102240 A JP 2000102240A JP 4478286 B2 JP4478286 B2 JP 4478286B2
- Authority
- JP
- Japan
- Prior art keywords
- blood flow
- frequency
- light receiving
- flow velocity
- light
- 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
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Eye Examination Apparatus (AREA)
Description
【発明の属する技術分野】
本発明は、ドップラ現象を利用して、被検眼の血管内の血流速度を計測する眼底血流計及び眼底血流速度測定方法に関するものである。
【0002】
【従来の技術】
従来からドップラ方式の眼血流計では、被測定血管に測定光を照射し、血管内を流れる赤血球等の粒子によってドップラシフトされた散乱信号光と、血管壁や周辺組織からのドップラシフトされていない散乱参照光とがミキシングされた光ビート信号を、2方向から2個の受光器により受光し、それらのFFT(高速フーリエ変換)波形の解析を行っている。このとき、図8に示すように血管内の血流の流れをポアゼイユの流れと仮定して、図9、図10に示すように、血管中心の最大流速に対応する最大ドップラシフト量であるカットオフ周波数Δfmax1、Δfmax2を求め、この値から最大血流速度を求めている。
【0003】
このカットオフ周波数の決定は、古くはオペレータが目視で判断して決定しているが、APPLIED OPTICS、Vol.27、No.6、pp.l126・1134(1988)「Retinal laser Doppler velocimetry: toward its computer-assisted clinical use」(B.L.Petrig、C.E,Riva)には、FFT波形のパワースペクトルがカットオフ周波数の個所で不連続に垂直に落ちる理想的なモデルを考慮し、カットオフ周波数を自動的に求める方法が開示されている。
【0004】
【発明が解決しようとする課題】
しかしながら上述の従来例では、実際のFFT波形のカットオフ部が垂直ではないために、精度良くカットオフ周波数を求めることができず、このため血流速を高精度で求めることができない。
【0005】
このため、本発明者は精度良くカットオフ周波数を求めることで、高精度に血流速を求めることができる眼底血流計を得るために、既に受光手段の受光面積を制限するための受光瞳の形状に関する情報、又は受光瞳の形状を基に算出した受光信号の理論的パワースペクトル形状の情報を記憶し、その情報を使用して受光信号を解析することにより、実際の周波数解析結果に近い形状を使用した解析を行う眼底血流計を提案している。
【0006】
本発明の目的は、更に精度良くカットオフ周波数を求め、高精度に血流速を測定できる眼底血流計及び眼底血流速度測定方法を提供することにある。
【0007】
【課題を解決するための手段】
上記目的を達成するための請求項1に係る本発明は、測定光を被検眼の血管に照射する照射手段と、前記測定光により血管内粒子から生ずる散乱光を異なる方向から受光する2個の受光手段と、該受光手段からの出力信号を解析して血流速度を算出するデータ処理手段とを有する眼底血流計において、前記データ処理手段は前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする眼底血流計である。
【0008】
請求項2に係る本発明は、前記受光手段の受光瞳を形成する形成部材と、受光瞳の形状に関する情報又は受光瞳の形状を基に算出した理論的パワースペクトル形状の情報を記憶する記憶手段とを有し、前記データ処理手段は前記記憶手段に記憶された情報を解析して、前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする請求項1に記載の眼底血流計である。
【0009】
請求項3に係る本発明は、前記データ処理手段は、前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、基準となる軸と前記血管とが成す角度に対してもフィッティングを行うことにより、血流速度を算出することを特徴とする請求項2に記載の眼底血流計である。
【0010】
請求項4に係る本発明は、測定光を被検眼の血管に照射し、前記測定光により血管内粒子から生ずる散乱光を異なる2方向から受光し、前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティッングして血流速度を算出することを特徴とする眼底血流速度測定方法である。
【0011】
請求項5に係る本発明は、記憶手段が記憶する受光瞳の形状に関する情報又は受光瞳の形状を基に算出した理論的パワースペクトル形状の情報を解析して、前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする請求項4に記載の眼底血流速度測定方法である。
【0012】
請求項6に係る本発明は、前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、基準となる軸と前記血管とが成す角度に対してもフィッティッングを行うことにより、血流速度を算出することを特徴とする請求項4に記載の眼底血流速度測定方法である。
【0013】
【発明の実施の形態】
本発明を図1〜図7に図示の実施例に基づいて詳細に説明する。
図1は実施例の眼底血流計の構成図を示し、白色光を発するタングステンランプ等から成る観察用光源1から被検眼Eと対向する対物レンズ2に至る照明光路上には、コンデンサレンズ3、例えば黄色域の波長光のみを透過するバンドパスフィルタ4、被検眼Eの瞳孔とほぼ共役な位置に設けられたリングスリット5、光路に沿って移動自在な固視標表示用素子である透過型液晶板6、リレーレンズ7、孔あきミラー8、黄色域の波長光を透過し他の光束を殆ど反射するバンドパスミラー9が順次に配列されている。なお、リングスリット5は被検眼Eの前眼部において、眼底照明光と眼底観察光を分離するためのものであり、必要な遮光領域を形成するものであれば、その形状や数は問題とならない。
【0014】
孔あきミラー8の背後には眼底観察光学系が構成されており、光路に沿って移動自在なフォーカスレンズ10、リレーレンズ11、スケール板12、接眼レンズ13が順次に配列され、検者眼eに至っている。バンドパスミラー9の反射方向の光路上には、イメージローテータ14、紙面に垂直な回転軸を有する両面研磨されたガルバノメトリックミラー15が配置され、ガルバノメトリックミラー15の下側反射面15aの反射方向には、光路に沿って移動自在なフォーカシングレンズ16が配置され、上側反射面15bの反射方向にはレンズ17、フォーカスユニット18が配置されている。
【0015】
なお、レンズ17の前側焦点面は被検眼Eの瞳孔と共役関係にあり、この焦点面に瞳孔上において非対称な形状とされたガルバノメトリックミラー15の後方には、凹面ミラー19が光軸上に同心的に配置され、ガルバノメトリックミラー15の上側反射面15bで反射されたレーザービームが、ガルバノメトリックミラー15の切欠部を通過するようにするために、ガルバノメトリックミラー15の上側反射面15bと下側反射面15aとを−1倍で結像するリレー光学系が構成されている。
【0016】
フォーカスユニット18においては、レンズ17と同一光路上にダイクロイックミラー20、集光レンズ21、レーザーダイオード等の測定用光源22が順次に配列され、ダイクロイックミラー20の反射方向の光路上にはマスク23、ミラー24が配置されており、このフォーカスユニット18は一体的に矢印で示す方向に移動可能とされている。更に、ミラー24の入射方向の光路上には、他の光源と異なる高輝度の例えば緑色光を発するトラッキング用光源25が配列されている。
【0017】
ガルバノメトリックミラー15の下側反射面15aの反射方向の光路上には、フォーカスレンズ16の後方に、ダイクロイックミラー26、拡大レンズ27、イメージインテンシファイヤ付の二次元撮像素子28が順次に配列され、血管検出系が構成されている。また、ダイクロイックミラー26の反射方向の光路上には、受光瞳を形成するミラー29a、29bを介してフォトマルチプライヤ30a、30bが配置され、測定用受光光学系が構成されている。なお、図示の都合上、全ての光路を同一平面上に示したが、ミラー29a、29b、フォトマルチプライヤ30a、30bはそれぞれ紙面に直交した方向に配置されている。
【0018】
二次元撮像素子28の出力はトラッキング制御部31に接続されており、トラッキング制御部31の出力はガルバノメトリックミラー15に接続されており、更にトラッキング制御部31は装置全体を制御するシステム制御部32に接続されている。また、システム制御部32にはフォトマルチプライヤ30a、30b、操作部33の出力が接続され、システム制御部32の出力はデータ処理部34に接続されている。
【0019】
観察用光源1から発した白色光はコンデンサレンズ3を通り、バンドパスフイルタ4により黄色の波長光のみが透過され、リングスリット5を通過した光束が透過型液晶板6を背後から照明し、リレーレンズ7を通って孔あきミラー8で反射される。その後に、黄色域の光のみがバンドパスミラー9を透過し、対物レンズ2を通り、被検眼Eの瞳孔上でリングスリット像として一旦結像した後に、眼底Eaをほぼ一様に照明する。このとき、透過型液晶板6には固視標が表示されており、照明光により眼底Eaに投影され、視標像として被検眼Eに呈示される。
【0020】
眼底Eaからの反射光は同じ光路を戻り、瞳孔上から眼底観察光光束として取り出され、孔あきミラー8の中心の開口部、フォーカスレンズ10、リレーレンズ11を通り、スケール板12で眼底像Ea’として結像した後に、接眼レンズ13を介して検者眼eにより観察される。検者はこの眼底像Ea’を観察しながら、装置のアライメントを行う。
【0021】
測定用光源22を発した測定光は集光レンズ21の上方を偏心して通過し、ダイクロイックミラー20を透過する。一方、トラッキング用光源25から発したトラッキング光はミラー24で反射された後に、マスク23で所望の形状に整形され、更にダイクロイックミラー20で反射されて、集光レンズ21によりマスク23の開口部中心と共役な位置にスポット状に結像している測定用光源22からの測定光と重畳される。
【0022】
測定光とトラッキング光はレンズ17を通り、ガルバノメトリックミラー15の上側反射面15bで一旦反射され、更に凹面ミラー19で反射され、再びガルバノメトリックミラー15の方に戻される。ここで、リレー光学系の機能により、ガルバノメトリックミラー15の上側反射面15bで反射された両光束は、ガルバノメトリックミラー15の切欠部の位置に戻されることになり、ガルバノメトリックミラー15に反射されることなく、イメージローテータ14に向かう。イメージローテータ14を経て、バンドパスミラー9により対物レンズ2の方向に偏向された両光束は、対物レンズ2を介して眼底Eaに照射される。
【0023】
このとき、トラッキング光はマスク23により、測定点を含みその血管をカバーする長方形の領域を照明するように、その大きさが血管走行方向300〜500μm程度、血管直角方向に500〜1200μm程度に整形されており、また測定光は測定する血管の太さ程度の50〜120μmの円形スポット、又は血管走行方向に長手方向を有する楕円形状とされている。
【0024】
眼底Eaでの散乱反射光は再び対物レンズ2で集光され、バンドパスミラー9で反射されてイメージローテータ14を通り、ガルバノメトリックミラー15の下側反射面15aで反射され、フォーカスレンズ16を通りダイクロイックミラー26において、測定光とトラッキング光とが分離される。
【0025】
そして、トラッキング光はダイクロイックミラー26を透過し、拡大レンズ27により二次元撮像素子28上で眼底観察光学系による眼底像Ea’よりも拡大された血管像として結像する。このときの撮像範囲はトラッキング光の照射範囲とほぼ同一の大きさである。この血管像信号はトラッキング制御部31に入力され、血管の位置信号に変換される。トラッキング制御部31はこの信号を使用して、ガルバノメトリックミラー15の回転角を制御し血管のトラッキングを行う。
【0026】
また、測定光とトラッキング光による眼底Eaでの散乱反射光の一部はバンドパスミラー9を透過し、孔あきミラー8の背後の眼底観察光学系に導かれ、トラッキング光はスケール板12上に棒状のインジケータとして結像し、測定光はこのインジケータの中心部にスポット像として結像する。これらの像は接眼レンズ13を介して眼底像及び視標像と共に観察される。このとき、インジケータの中心には測定ビームのスポット像が重畳して観察される。インジケータは操作部33によってガルバノメトリックミラー15を回転することにより、眼底Ea上を一次元に移動することができる。
【0027】
測定に際して、検者は先ず眼底像のピント合わせを行う。操作部33のフォーカスノブを調整すると、図示しない駆動手段により透過型液晶板6、フォーカスレンズ10、16、フォーカスユニット18が連動して光路に沿って移動する。眼底像のピントが合うと、透過型液晶板6、スケール板12、二次元撮像素子28は同時に眼底Eaと共役になる。
【0028】
検者は眼底像のピントを合わせた後に、被検眼Eの視線を誘導して観察領域を変更し、測定対象とする血管を適当な位置に移動するために操作部33を操作する。システム制御部32は透過型液晶板6を制御し視標像を移動し、イメージローテータ14を回転して測定対象とする血管の走行方向に対して、フォトマルチプライヤ30a、30bの中心を結んだ線が平行になるように操作する。このとき、ガルバノメトリックミラー15を回転することにより、二次元撮像素子28の画素配列の垂直方向と測定ビームの移動方向は、同時にこれと直角の血管に対して垂直な方向に調整される。
【0029】
検者はトラッキングを開始してその良否を確認した後に、操作部33の測定スイッチを押して測定を開始する。この測定の間は、測定ビームはトラッキング制御部31の働きにより血管上に保持されるが、その散乱反射光はダイクロイックミラー26、ミラー29a、29bにより反射され、フォトマルチプライヤ30a、30bに受光される。フォトマルチプライヤ30a、30bの出力はそれぞれシステム制御部32に入力され、FFT処理などの周波数解析が行われる。これにより得られたFFT波形はデータ処理部34において解析され、眼底Eaの血流速度が求められる。
【0030】
図2はデータ処理部34の動作のフローチャート図を示し、先ずステップS1で、フォトマルチプライヤ30a、30bからの信号をFFT処理したそれぞれのFFT波形に対して、スムージングなどの処理を行う。
【0031】
図3は測定対象とする血管Evの測定部位V、測定用光源22、フォトマルチプライヤ30a、30b、そしてミラー29a、29bが形成する受光瞳A1、A2の関係を示したものである。受光瞳A1、A2は反射ミラーに限らず、絞り或いはフォトマルチプライヤのセンサ部等で形成してもよい。角度の基準とする眼底血流計の光軸Oと測定対象とする血管Evの走行方向との成す角、即ち血管角度をθ、光軸Oと測定用光源22の測定部位Vへの入射方向との成す角つまり入射角度をγ0、光軸Oと測定部位Vから受光瞳A1、A2への受光方向との成す角、即ち受光角度をβとする。
【0032】
また、図4は受光瞳A1、A2の位置関係と平面図を示し、本実施例では半円形のミラー29a、29bを用い、受光瞳A1、A2は半径Rの半円形状であるとし、これら受光瞳A1、A2の形状に関する情報は、システム制御部32に記憶されている。β0は受光瞳A1、A2の重心G1、G2の受光角度である。ここで、図3の血管角度θは血管Evの走行方向によって異なる変数であり、図3の受光角度βは図4の受光瞳A1、A2についてβ1〜β2の角度の中で受光方向によって異なる変数である。なお、これらの角度θ、γ0、βは人眼中での角度に換算したものを使用する。
【0033】
血管Ev内を流れる赤血球等の散乱粒子が速度ベクトルVeで進んでいるものとし、入射光及び散乱光の進行方向の波数ベクトルをKl、Ksとすると、ドップラ効果の原理からドップラシフト量Δfは次式で与えられる。
Δf=(Ks−Kl)・Ve/(2π) …(1)
【0034】
図2に戻り、血管角度θとして予め定められた例えば75度〜105度のような取り得る角度の範囲θs〜θeの内の初めの角度θsを、仮に定めた血管角度とする。血管角度θが与えられると、受光瞳A2、A1の重心G2、G1の受光方向の光ビート信号によるFFTのそれぞれのカットオフ周波数の理論的な比率Ratioは、(1)式から、角度θ、γ0、β0により次式のように表され、ステップS2でこのRatioの値を計算する。
Ratio=|{cos(θ+β0)+cos(θ−γ0)}/{cos(θ−β0)+cos(θ−γ0)}|
…(2)
【0035】
次に、予め定められた周波数N1の範囲Fs〜Feの内の初めの周波数Fsを、重心G1での仮に定めたカットオフ周波数とし、ステップS3で次式により、重心G2で仮に定めたカットオフ周波数N2を求める。
N2=N1×Ratio …(3)
【0036】
ステップS4で、ステップS1で処理したフォトマルチプライヤ30a、30bからのFFT波形のそれぞれについて、モデル形状であるフィッティング曲線を求める。ここで受光瞳が点であれば、FFT波形の理想的なモデル形状は、図5の点線FLに示すような周波数Fsにおいてパワースペクトルが不連続に落ちる形状になるが、実際には受光瞳が面積を持っているために、周波数Fsで不連続にはならず、実線SLのFFT波形のように、或る曲線で急激に落ちる形状となる。
【0037】
フォトマルチプライヤ30a、30bで受光される光ビート信号の内、受光角度βの受光方向の光ビート信号によるFFTのカットオフ周波数Δf1、Δf2は、(1)式からそれぞれの受光角度β0におけるカットオフ周波数Δfo1、Δfo2と角度θ、γ0、β、β0により、次式のように表される。
Δf1=Δfo1・|{cos(θ−β)+cos(θ−γ0)}
/{cos(θ−β0)+cos(θ−γ0)}| …(4)
Δf2=Δfo2・{cos(θ+β)+cos(θ−γ0)}
/{cos(θ+β0)+cos(θ−γ0)}| …(5)
【0038】
また、受光瞳A1、A2の内の受光角度βの部分から形成される光ビート信号のFFTのパワースペクトルPは、その部分の面積に比例し、次の関係式が成立する。なお、Lは受光瞳A1、A2が作る平面と測定部位Vとの距離である。
P∝[R2−{L(tanβ−tanβ0)}2]1/2 …(6)
【0039】
受光瞳A1、A2の全体で形成される光ビート信号のFFTのパワースペクトルPは、図4の受光角度β1 〜β2 までのパワースペクトルPの重ね合わせとなり、FFT波形のパワースペクトルPがカットオフ周波数のところで連続的に落ちる曲線の形状、つまり理論カットオフ形状を求めることができる。
【0040】
なお本実施例では、計算により理論カットオフ形状を求めているが、各血管角度に対する理論カットオフ形状の数値を記憶しておく方法でもよい。また精度は落ちるが、簡易的には図5の破線FLに示すようなパワースペクトルが不連続に落ちる形状であるFFT波形の理想的なモデル形状を用いることもできる。
【0041】
図6はフォトマルチプライヤ30aで受光される光ビート信号によるFFT波形SLとフィッティング曲線の例を示す。破線FL1は受光出力信号の周波数解析結果のフィッティングのために、θsを仮に定めた血管角度、Fsを仮に定めたカットオフ周波数としたときのフィッティング曲線を計算により求めたものである。
【0042】
理論カットオフ形状の曲線部は上述の重ね合わせで求め、理論カットオフ形状の落ち始めの周波数をFaとすると、周波数Faよりも低い個所のパワースペクトルPS1は、実線SLに示す実際のFFT波形の周波数Faよりも低い個所のパワースペクトルPの平均値を使用している。また、理論カットオフ形状の終端の周波数をFbとすると、周波数Fbよりも高い個所のパワースペクトルPS2は、実際のFFT波形の内の周波数の十分高いノイズによる部分の平均値を使用している。これらの平均値を求める方法は計算が単純なので、短時間で計算が可能である。
【0043】
このように、理論カットオフ形状の終端がPS2になるように、理論カットオフ形状にノイズ成分としてPS2を上乗せし、更に理論カットオフ形状の高さがPS1−PS2になるように、拡大又は縮小を行ってフィッティング曲線を作成する。このフィッティング曲線のカットオフ周波数部は、実際のFFT波形に近い形状であり、より精度良くカットオフ周波数を求めることができ、高精度の血流速度を得ることができる。
【0044】
図7はフィッティング曲線を求める他の方法を示したものである。理論カットオフ形状の曲線部と周波数Fbよりも高い個所のパワースペクトルPS2の求め方は図6と同様であるが、周波数Faよりも低い個所のパワースペクトルPS1の求め方が異なる。ここでは、実線SLに示す実際のFFT波形のパワースペクトルPからノイズ成分と考えられるパワースペクトルPS2を差し引いた部分を全周波数に渡って積分した値と、フィッティング曲線のパワースペクトルPからパワースペクトルPS2を差し引いた部分を、周波数Fbよりも低い個所について積分した値が等しくなるように、即ちパワースペクトルPの合計が等しくなるようにパワースペクトルPS1を決めている。
【0045】
この場合には、仮のカットオフ周波数Fsが実際のFFT波形のカットオフ周波数から離れていると、フィッティング曲線の形状は実際のFFT波形と大きく異なるが、実際のカットオフ周波数に近付くにつれてフィッティング曲線の形状も実際のFFT波形に近付く。即ち、この方法では図2のステップS6で数値が最小となる周波数を求める際に収束性が良くなり、より容易にカットオフ周波数を求めることができる。
【0046】
図2のフローチャート図に戻り、ステップS4でフォトマルチプライヤ30bで受光される光ビート信号によるFFT波形のフィッティング曲線も同様に求める。そしてステップS5において、上述のようにして求めたフォトマルチプライヤ30a、30bからのFFT波形それぞれについてのフィッティング曲線と、実際のFFT波形とのパワースペクトルPの値の差の二乗和を計算し、両者の合計Toを求める。ステップS3〜S5を、予め定められた間隔で周波数Fsの次の周波数についても実行し、最終の周波数Feまで繰り返す。即ち、周波数に対するフィッティングを行う。
【0047】
次に、予め定められた間隔で血管角度θsの次の血管角度についてもステップS2と、Fs〜Fe間で予め定められた間隔の各周波数について、ステップS3〜S5を実行し、最終の血管角度θeまで繰り返し、血管角度に対するフィッティングを行う。最終的に、血管角度θと周波数N1を変数とする二乗和Toの値が得られる。
【0048】
その後に、ステップS6で二乗和Toの最小値を求め、そのときの周波数N1、N2をそれぞれフォトマルチプライヤ30a、30bで受光される受光角度β0におけるカットオフ周波数Δfmax1、Δfmax2とする。このときのフィッティング曲線はフォトマルチプライヤ30aについての図6の実線FL0のようにFFT波形SLにより近い形状となり、カットオフ周波数が精度良く求められていることが分かる。フォトマルチプライヤ30bについても同様である。
【0049】
次に、ステップS7で次式により最大血流速度Vmaxを求める。ここで、λは測定光束の波長、nは測定部位の屈折率である。
Vmax={λ/(n・2β0)}・|Δfmax1−Δfmax2| …(7)
【0050】
このようにして、測定光により血管内粒子から生ずる散乱光を異なる方向から受光する2個の受光手段からの出力信号の周波数解析結果を、カットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして、カットオフ周波数を求めることにより、2倍のデータの出力信号に対してフィッティングしてカットオフ周波数を求めることになるので、容易に精度良くカットオフ周波数が求められ、高精度に最大血流速度Vmaxを得ることができる。
【0051】
なお、本実施例では血管角度θをもパラメータとしてフィッティングを行い、フィッティングが最適となるようにしているので、より精度良くカットオフ周波数が求められているが、簡単のためには特開平10−80398号公報に記載されているような方法で角度θを求め、その値を使用して周波数軸に対してフィッティングして血流速度を算出してもよく、この場合は精度は落ちるが計算量が少なくなる。
【0052】
また、両フォトマルチプライヤ30a、30bの受光出力信号のS/Nが十分に良好な場合には、ステップS6で二乗和Toが最小値となるときの血管角度θを測定部位Vの血管角度θvとし、上述のようにして求めたΔfmax1、Δfmax2を用いて、次式のように等しくなるべき2個の最大血流速度を求め、それらを平均するなどすれば、より精度良く最大血流速度を求めることができる。
Vmax=λ・Δfmax1/{n|cos(θv−β0)+cos(θv−γ0)|} …(8)
Vmax=λ・Δfmax2/{n|cos(θv+β0)+cos(θv−γ0)|} …(9)
【0053】
以上の実施例では、眼底血管の血流を測定する眼底血流計を例にしているが、光学系を多少変更した強膜上の血管の血流を測定する装置にも応用可能である。
【0054】
【発明の効果】
以上説明したように、請求項1及び請求項4に係る眼底血流計及び眼底血流速度測定方法は、測定光により血管内粒子から生ずる散乱光を異なる方向から受光する2個の受光手段からの出力信号の周波数解析結果を、カットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングしてカットオフ周波数を求めることにより、2倍のデータの出力信号に対してフィッティングしてカットオフ周波数を求めることになるので、容易に精度良くカットオフ周波数、血流速を求めることが可能となる。
【0055】
また、請求項2及び請求項5に係る本発明は、データ処理手段が情報記憶手段に記憶した受光瞳の形状に関する情報又は受光瞳の形状を基に算出した受光信号の理論的パワースペクトル形状の情報を解析して、2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、周波数軸に対してフィッティングして血流速度を算出することにより、更に精度良くカットオフ周波数、血流速を求めることが可能となる。
【0056】
請求項3及び請求項6に係る本発明は、データ処理手段は2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、更に眼底血流計の基準軸と前記血管とが成す角度に対してもフィッティングを行うことにより、容易により精度の良いカットオフ周波数、血流速を求めることができる。
【図面の簡単な説明】
【図1】本実施例の眼底血流計の構成図である。
【図2】データ処理動作のフローチャート図である。
【図3】測定部位、測定用光源、フォトマルチプライヤ、受光瞳の関係の説明図である。
【図4】受光瞳の位置関係の説明図である。
【図5】FFT波形のモデル形状のグラフ図である。
【図6】フィッティング曲線のグラフ図である。
【図7】他のフィッティング曲線のグラフ図である。
【図8】血管内の流速分布の説明図である。
【図9】光ビート信号のFFT波形のグラフ図である。
【図10】光ビート信号のFFT波形のグラフ図である。
【符号の説明】
1 観察用光源
2 対物レンズ
6 透過型液晶板
9 バンドパスミラー
14 イメージローテータ
15 ガルバノメトリックミラー
18 フォーカスユニット
20 ダイクロイックミラー
22 測定用光源
25 トラッキング用光源
28 二次元撮像素子
30a、30b フォトマルチプライヤ
31 トラッキング制御部
32 システム制御部
33 操作部
34 データ処理部
Claims (6)
- 測定光を被検眼の血管に照射する照射手段と、前記測定光により血管内粒子から生ずる散乱光を異なる方向から受光する2個の受光手段と、該受光手段からの出力信号を解析して血流速度を算出するデータ処理手段とを有する眼底血流計において、前記データ処理手段は前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする眼底血流計。
- 前記受光手段の受光瞳を形成する形成部材と、受光瞳の形状に関する情報又は受光瞳の形状を基に算出した理論的パワースペクトル形状の情報を記憶する記憶手段とを有し、前記データ処理手段は前記記憶手段に記憶された情報を解析して、前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする請求項1に記載の眼底血流計。
- 前記データ処理手段は、前記2個の受光手段からの出力信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、基準となる軸と前記血管とが成す角度に対してもフィッティングを行うことにより、血流速度を算出することを特徴とする請求項2に記載の眼底血流計。
- 測定光を被検眼の血管に照射し、前記測定光により血管内粒子から生ずる散乱光を異なる2方向から受光し、前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティッングして血流速度を算出することを特徴とする眼底血流速度測定方法。
- 記憶手段が記憶する受光瞳の形状に関する情報又は受光瞳の形状を基に算出した理論的パワースペクトル形状の情報を解析して、前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように周波数軸に対してフィッティングして血流速度を算出することを特徴とする請求項4に記載の眼底血流速度測定方法。
- 前記2方向から受光した信号の周波数解析結果をカットオフ周波数の理論的な比率を保つように、基準となる軸と前記血管とが成す角度に対してもフィッティッングを行うことにより、血流速度を算出することを特徴とする請求項4に記載の眼底血流速度測定方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000102240A JP4478286B2 (ja) | 2000-04-04 | 2000-04-04 | 眼底血流計及び眼底血流速度測定方法 |
US09/650,348 US6454722B1 (en) | 1999-08-31 | 2000-08-29 | Doppler velocimeter for blood flow |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000102240A JP4478286B2 (ja) | 2000-04-04 | 2000-04-04 | 眼底血流計及び眼底血流速度測定方法 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2001276002A JP2001276002A (ja) | 2001-10-09 |
JP2001276002A5 JP2001276002A5 (ja) | 2007-06-07 |
JP4478286B2 true JP4478286B2 (ja) | 2010-06-09 |
Family
ID=18616161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000102240A Expired - Fee Related JP4478286B2 (ja) | 1999-08-31 | 2000-04-04 | 眼底血流計及び眼底血流速度測定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4478286B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5438744B2 (ja) * | 2011-11-25 | 2014-03-12 | 国立大学法人 東京大学 | 血流可視化診断装置、及びプログラム |
-
2000
- 2000-04-04 JP JP2000102240A patent/JP4478286B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2001276002A (ja) | 2001-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4402601A (en) | Fundus camera-based retinal laser doppler velocimeter | |
US5615683A (en) | Ophthalmologic measuring apparatus | |
JP3647164B2 (ja) | 眼科測定装置 | |
JPH01262835A (ja) | 眼科診断方法及び装置 | |
JPH04242628A (ja) | 眼科測定装置 | |
JP2002034921A (ja) | 眼底検査装置 | |
US6454722B1 (en) | Doppler velocimeter for blood flow | |
US6332683B1 (en) | Fundus examination apparatus | |
JPH08215150A (ja) | 眼科診断装置 | |
US6607272B1 (en) | Retinal blood flow measuring apparatus using a laser beam | |
JP4478286B2 (ja) | 眼底血流計及び眼底血流速度測定方法 | |
JP2003019116A (ja) | 眼科測定装置 | |
JP3997030B2 (ja) | 眼底血流計 | |
JPH0731596A (ja) | 眼底血流計 | |
US6302850B1 (en) | Fundus blood flow metering method | |
JP2021191552A (ja) | 眼科検査装置 | |
US6685650B2 (en) | Fundus blood flowmeter | |
JP4194187B2 (ja) | 眼底血流計 | |
JP2001112717A (ja) | 眼底血流計及びその信号処理方法 | |
JP2000287946A (ja) | 眼底血流計 | |
JPH01256924A (ja) | 眼科診断方法及び装置 | |
JP4371496B2 (ja) | 眼底血流計及びその信号処理方法 | |
JP4724320B2 (ja) | 眼血流計 | |
JP2000296108A (ja) | 眼科検査装置 | |
EP4166065A1 (en) | Ophthalmic device, method for controlling ophthalmic device, and program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070326 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070326 |
|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20100218 |
|
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: 20100309 |
|
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: 20100315 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130319 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140319 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |