JP2013140077A - 波長走査型光干渉断層計及びその位相安定化プログラム - Google Patents

波長走査型光干渉断層計及びその位相安定化プログラム Download PDF

Info

Publication number
JP2013140077A
JP2013140077A JP2012000313A JP2012000313A JP2013140077A JP 2013140077 A JP2013140077 A JP 2013140077A JP 2012000313 A JP2012000313 A JP 2012000313A JP 2012000313 A JP2012000313 A JP 2012000313A JP 2013140077 A JP2013140077 A JP 2013140077A
Authority
JP
Japan
Prior art keywords
interference signal
spectrum
phase
oct
shift
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.)
Granted
Application number
JP2012000313A
Other languages
English (en)
Other versions
JP5900950B2 (ja
Inventor
Yoshiaki Yasuno
嘉晃 安野
Masahide Ito
雅英 伊藤
Ying Zhou Hong
暎周 洪
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.)
University of Tsukuba NUC
Original Assignee
University of Tsukuba NUC
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 University of Tsukuba NUC filed Critical University of Tsukuba NUC
Priority to JP2012000313A priority Critical patent/JP5900950B2/ja
Priority to PCT/JP2012/082460 priority patent/WO2013103080A1/ja
Priority to US14/368,456 priority patent/US9354038B2/en
Priority to EP12864429.1A priority patent/EP2801814B1/en
Priority to CN201280062822.1A priority patent/CN104053980B/zh
Publication of JP2013140077A publication Critical patent/JP2013140077A/ja
Application granted granted Critical
Publication of JP5900950B2 publication Critical patent/JP5900950B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02055—Reduction or prevention of errors; Testing; Calibration
    • G01B9/02075—Reduction or prevention of errors; Testing; Calibration of particular errors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02083—Interferometers characterised by particular signal processing and presentation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/0209—Low-coherence interferometers
    • G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Ophthalmology & Optometry (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

【課題】ＳＳ−ＯＣＴにおける光源の波長スキャンと、光検知器でスペクトル干渉信号としてデータを収集するタイミングとの間のジッターを高価で複雑なハードウェアを追加的に使用せず、位相を安定化させる。
【解決手段】ＳＳ−ＯＣＴの波長走査型光源２から出射されて分割され、固定参照鏡８で反射された参照光と、被計測物体６で反射された物体光が重ねられて光検知器１５で検出されたスペクトル干渉信号に基づき、断層画像を生成するコンピュータ１６を、第１の補正手段で大まかに補正し、第２の補正手段でさらに詳細に補正させるように機能させて、ＳＳ−ＯＣＴの位相を安定化する。
【選択図】図１

Description

本発明は、光源の波長を走査してスペクトル干渉信号を得る波長走査型光干渉断層計（Swept Source Optical Coherence Tomography、略して「ＳＳ−ＯＣＴ」という。）及びその位相安定化プログラムに関し、特に、追加的なハードウェアを使用せず、一般的に用いられているＳＳ−ＯＣＴと同じ構成のみにより位相を安定化させる技術に関する。
従来、医療分野等で用いられる非破壊断層計測技術の１つとして、時間的に低コヒーレンスな光をプローブ（探針）として用いる光断層画像化法「光コヒーレンストモグラフィー」（ＯＣＴ）がある（特許文献１参照）。ＯＣＴは、光を計測プローブとして用いるため、被計測物体の屈折率分布、分光情報、偏光情報（複屈折率分布）等が計測できるという利点がある。
基本的なＯＣＴ４３は、マイケルソン干渉計を基本としており、その原理を図３で説明する。光源４４から射出された光は、コリメートレンズ４５で平行化された後に、ビームスプリッター４６により参照光と物体光に分割される。物体光は、物体アーム内の対物レンズ４７によって被計測物体４８に集光され、そこで散乱・反射された後に再び対物レンズ４７、ビームスプリッター４６に戻る。
一方、参照光は参照アーム内の対物レンズ４９を通過した後に参照鏡５０によって反射され、再び対物レンズ４９を通してビームスプリッター４６に戻る。このようにビームスプリッター４６に戻った物体光と参照光は、物体光とともに集光レンズ５１に入射し光検出器５２（フォトダイオード等）に集光される。
ＯＣＴの光源４４は、時間的に低コヒーレンスな光（異なった時刻に光源から出た光同士は極めて干渉しにくい光）の光源を利用する。時間的低コヒーレンス光を光源としたマイケルソン型の干渉計では、参照アームと物体アームの距離がほぼ等しいときにのみ干渉信号が現れる。この結果、参照アームと物体アームの光路長差（τ）を変化させながら、光検出器５２で干渉信号の強度を計測すると、光路長差に対する干渉信号（インターフェログラム）が得られる。
そのインターフェログラムの形状が、被計測物体４８の奥行き方向の反射率分布を示しており、１次元の軸方向走査により被計測物体４８の奥行き方向の構造を得ることができる。このように、ＯＣＴ４３では、光路長走査により、被計測物体４８の奥行き方向の構造を計測できる。
このような軸方向の走査のほかに、横方向の機械的走査を加え、２次元の走査を行うことで被計測物体の２次元断面画像が得られる。この横方向の走査を行う走査装置としては、被計測物体を直接移動させる構成、物体は固定したままで対物レンズをシフトさせる構成、被計測物体も対物レンズも固定したままで、対物レンズの瞳面付近においたガルバノミラーの角度を回転させる構成等が用いられている。
以上の基本的なＯＣＴが発展したものとして、分光器を用いてスペクトル信号を得るスペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）と、光源の波長を走査してスペクトル干渉信号を得る波長走査型ＯＣＴ（Swept Source OCT、略して「ＳＳ−ＯＣＴ」という。）がある。ＳＤ−ＯＣＴには、フーリエドメインＯＣＴ（Fourier Domain OCT、略して「ＦＤ−ＯＣＴ」という。特許文献２参照）、及び偏光感受型ＯＣＴ（Polarization-Sensitive OCT、略して「ＰＳ−ＯＣＴ」という。特許文献３参照）がある。
ＦＤ−ＯＣＴは、被計測物体からの反射光の波長スペクトルを、スペクトロメーター（スペクトル分光器）で取得し、このスペクトル強度分布に対してフーリエ変換することで、実空間（ＯＣＴ信号空間）上での信号を取り出すことを特徴とするものであり、このＦＤ−ＯＣＴは、奥行き方向の走査を行う必要がなく、ｘ軸方向の走査を行うことで被計測物体の断面構造を計測可能である。
ＳＳ−ＯＣＴは、高速波長スキャニングレーザーにより光源の波長を変え、スペクトル信号と同期取得された光源走査信号を用いて干渉信号を最配列し、信号処理を加えることで３次元光断層画像を得るものである。なお、光源の波長を変える手段として、モノクロメーターを利用したものでも、ＳＳ−ＯＣＴとして利用可能である。
ドップラー光コヒーレンストモグラフィー（ドップラーＯＣＴ）による、網膜の血流の分布の計測が知られている。これは、上記ＦＤ−ＯＣＴ等を用いて、網膜の血流の分布が計測を行うことのできる手段であり、同様に、スペクトルドメインＯＣＴを使うことによって、横断面網膜血流画像形成が得られ、また次元の網膜の脈管構造も観察することができる。
本発明者らは、従来よりドップラーＯＣＴに着目し、生体内、特に眼底の血流を非侵襲で計測する手法の研究・開発を行って来た。従来本発明者らはＳＤ−ＯＣＴを技術基盤上にドップラーＯＣＴを実装し、眼底血流検査に成功してきたが、近年ＯＣＴの技術基板はＳＤ−ＯＣＴから、次世代であるＳＳ−ＯＣＴに移行しつつある。
特開２００２−３１０８９７号公報 特開平１１−３２５８４９号公報 特開２００４−０２８９７０号公報
B. Vokoc et.al., "Phase-resolved optical frequency domain imaging,"OpticsExpress 13, 5483- (2005).
本発明者らはＳＳ−ＯＣＴによってもドップラー計測の研究開発を行ってきたが、ドップラー計測には、困難な問題がある。これは、ＳＳ−ＯＣＴ特有の位相不安定性に起因するものである。
ＳＳ−ＯＣＴの位相不安定性を補正する方法は Harvard Medical School により提案されているものの（非特許文献１参照）、この技術は、複雑なハードウェアを必要とし製品のコストの増大がさけられない。
本発明は、上記ＳＳ−ＯＣＴの位相不安定性の問題を解決することを目的とするものであり、従来同様の標準的なＳＳ−ＯＣＴの構成のみを前提とし、そこから得られる信号を適切に処理することで高い位相の安定性を得ることを課題とするものであり、これにより、従来のＳＳ−ＯＣＴを利用して、ＳＳ−ＯＣＴによる計測及びドップラー計測を実現しようとするものである。
本発明は上記課題を解決するために、ＳＳ−ＯＣＴにおいて、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光とが重ねられて光検知器で検出されたスペクトル干渉信号に基づき、断層画像を生成するコンピュータを、スペクトル干渉信号の位相データを補正する補正手段として機能させるＳＳ−ＯＣＴのスペクトル干渉信号の位相を安定化するプログラムであって、前記補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（２’）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴの位相を安定化するプログラムを提供する。
Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（２’）
但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、＊ は複素共約を示す。
本発明は上記課題を解決するために、ＳＳ−ＯＣＴにおいて、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光とが重ねられて光検知器で検出されたスペクトル干渉信号に基づき、断層画像を生成するコンピュータを、スペクトル干渉信号の位相データを補正する第１の補正手段と第２の補正手段として順次機能させるＳＳ−ＯＣＴのスペクトル干渉信号の位相を安定化するプログラムであって、前記第１の補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（２’）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行い、前記第２の補正手段は、前記２つのＡ−スキャンで得られたスペクトル干渉信号に係るドップラー位相シフトを示す下記の式（９）において、ドップラー信号ｖｚ（ζ）がない領域における、ζの切片ｂｍ＝４πτｖｂ／λＣ、傾きａｍ＝−２πβ’／Ｎの１次関数について、下記の式（１０）に示す誤差エネルギーを最小とする傾きａｍと切片ｂｍを求めることで、バルクモーション及びジッターによるスペクトルシフトβの残差β’を求めて、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴの位相を安定化するプログラムを提供する。
Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（２’）
但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、＊ は複素共約を示す。
Δφ（ζ）＝（４πτ／λＣ）ｎｖｚ（ζ）＋（４πτ／λＣ）ｖｂ −２πζβ’／Ｎ
・・・（９）
但し、τは２つのＡ−スキャン間の時間、λＣ は光源の中心波長、ｎは屈折率、変数ζの関数ｖｚ（ζ）は光軸方向の速度、ｖｂはバルクモーションの速度、β’はジッターによるスペクトルシフトの残差、Ｎはサンプリング点数、ζの関数ｖｚ（ζ）は光軸方向の速度、ζは被計測物体の深さを示す変数をそれぞれ示し、ζ＝０，１，・・・・Ｎ−１である。
Ｒ２ ２＝ΣζＷｍ（ζ）[Δφ（ζ）−（ａｍζ＋ｂｍ）]２ ・・・（１０）
但し、ｍ＝０，１，２…… であり、ζは被計測物体の深さを示す変数であり、Ｗｍ（ζ）は重みである。
本発明は上記課題を解決するために、波長走査型光源と、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光が重ねられたスペクトル干渉信号を検出する光検知器と、光検知器検出されたスペクトル干渉信号に基づき断層画像を生成するコンピュータと、を備えたＳＳ−ＯＣＴであって、
前記コンピュータは、スペクトル干渉信号の位相データを補正してＳＳ−ＯＣＴの位相を安定化する補正手段を備えており、前記補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（２’）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴを提供する。
Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（２’）
但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、＊ は複素共約を示す。
本発明は上記課題を解決するために、波長走査型光源と、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光が重ねられたスペクトル干渉信号を検出する光検知器と、光検知器検出されたスペクトル干渉信号に基づき断層画像を生成するコンピュータと、を備えたＳＳ−ＯＣＴであって、前記コンピュータは、スペクトル干渉信号の位相データを順次補正してＳＳ−ＯＣＴの位相を安定化する第１の補正手段と第２の補正手段を備えており、前記第１の補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（２’）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行い、前記第２の補正手段は、前記２つのＡ−スキャンで得られたスペクトル干渉信号に係るドップラー位相シフトを示す下記の式（９）において、ドップラー信号ｖｚ（ζ）がない領域における、ζの切片ｂｍ＝４πτｖｂ／λＣ、傾きａｍ＝−２πβ’／Ｎの１次関数について、下記の式（１０）に示す誤差エネルギーを最小とする傾きａｍと切片ｂｍを求めることで、バルクモーション及びジッターによるスペクトルシフトβの残差β’を求めて、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴを提供する。
Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（２’）
但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、＊ は複素共約を示す。
Δφ（ζ）＝（４πτ／λＣ）ｎｖｚ（ζ）＋（４πτ／λＣ）ｖｂ−２πζβ’／Ｎ
・・・（９）
但し、τは２つのＡ−スキャン間の時間、λＣ は光源の中心波長、ｎは屈折率、変数ζの関数ｖｚ（ζ）は光軸方向の速度、ｖｂはバルクモーションの速度、β’はジッターによるスペクトルシフトの残差、Ｎはサンプリング点数、ζの関数ｖｚ（ζ）は光軸方向の速度、ζは被計測物体の深さを示す変数をそれぞれ示し、ζ＝０，１，・・・・Ｎ−１である。
Ｒ２ ２＝ΣζＷｍ（ζ）[Δφ（ζ）−（ａｍζ＋ｂｍ）]２ ・・・（１０）
但し、ｍ＝０，１，２…… であり、ζは被計測物体の深さを示す変数であり、Ｗｍ（ζ）は重みである。
本発明に係る波長走査型光干渉断層計及びその位相安定化プログラによれば、次のような効果が生じる。
（１）位相安定化ミラー、複数チャンネルデジタイザ等の高価で複雑なハードウェアを追加的に使用せず、一般的に用いられているＳＳ−ＯＣＴと同じ構成のみにより位相を安定化させることができる。
（２）ＳＳ−ＯＣＴで取得された複数のＡ−スキャン（被計測物体の深さ方向のスキャン）のスペクトル干渉信号の時間ジッター（光源の波長スキャンと、光検知器でスペクトル干渉信号の収集のタイミングの不一致に基づく時間軸方向での信号波形の揺らぎ）を補正することができる。
（３）ＳＳ−ＯＣＴで生体を撮影した際、生体の動き（バルクモーション）によるドップラー信号のアーチファクト（artifact（人工産物）：偽の所見）を検出し、補正することができる。
（４）ＳＳ−ＯＣＴで生体を撮影した際、そのドップラーシフト情報を用いて、時間ジッターをさらに高精度に補正することができる。
本発明に係るＰＳ−ＯＣＴシステムの全体構成を示す図である。 本発明に係る位相安定化プログラム及び該プログラムを搭載したＳＳ−ＯＣＴの効果を実証するための実験の結果を示すグラフである。 従来の基本的なＯＣＴを説明する図である。
本発明に係る波長走査型光干渉断層計及びその位相安定化プログラムを実施するための実施するための形態を実施例に基づいて図面を参照して、以下に説明する。
まず、本発明におけるＳＳ−ＯＣＴ（波長走査型光干渉断層計）の本体の構成を説明する。
図１は、ＳＳ−ＯＣＴ１の基本構成を示す図である。波長走査型光源２から出射された出力光を、ファイバー３を通してファイバーカップラー４に送る。この出力光を、ファイバーカップラー４において、ファイバー５を通して被計測物体６への照射する物体光と、ファイバー７を通して固定参照鏡８に照射する参照光に分割する。
物体光は、ファイバー５、レンズ９、角度が可変な走査鏡１０及びレンズ１１を介して、被計測物体６に照射、反射され、同じルートでファイバーカップラー４に戻る。参照光は、ファイバー７、レンズ１２及びレンズ１３を介して固定参照鏡８に照射、反射されて同じルートでファイバーカップラー４に戻る。
そして、これらの物体光と参照光はファイバーカップラー４で重ねられ、ファイバー１４を通して光検知器１５（ＰＤ（フォトダイオード）等のポイントセンサが使用される。）に送られ、スペクトル干渉信号として検出され、コンピュータ１６に取り込まれる。
光検知器１５における検知出力に基づいて、被計測物体６の奥行き方向（Ａ方向）と走査鏡の走査方向（Ｂ方向）の断面画像が形成される。コンピュータ１６にはディスプレー１７が接続されている。
波長走査型光源２は、時間的に波長を変化させて走査する光源であり、即ち波長が時間依存性を有する光源である。これにより、Ａ−スキャン（被計測物体の奥行き方向への走査）をするために、固定参照鏡８を走査（光軸方向に移動して行う光軸方向の走査）することなく、被計測物体６の奥行き方向の反射率分布を得て奥行き方向の構造を取得することができ、Ｂ−スキャン（光軸方向と直交する１次方向の走査）をするだけで、二次元の断層画像を形成することができる。
ところで、ＳＳ−ＯＣＴでは、その光源は時間的に波長を変化させて走査するので、この光源の波長スキャン（波長の変化のタイミング）と、光検知器でスペクトル干渉信号としてデータを収集するタイミングとの間の（不一致により生じる）ジッター（Jitter：時間軸方向での信号波形の揺らぎや、その揺らぎにより生じる画像の乱れ）が問題になる。
このジッターは、スペクトルサンプリングのランダムなシフトをもたらし、結果、ＳＳ−ＯＣＴで取得されるスペクトル干渉信号のジッターとなる。
ドップラー計測では、ドップラー信号は２つのＡ−スキャン（Ｂ−スキャンで得られる第１のＡ−スキャンＡ１と第２のＡ−スキャンＡ２）で取得されるスペクトル干渉信号の位相差によって得られるため、上記ＳＳ−ＯＣＴで取得されるスペクトル干渉信号のジッターは、ドップラー計測におけるデータ誤差に直接影響する。
本発明は、ＳＳ−ＯＣＴを構成するコンピュータに搭載され、そのコンピュータを、ＳＳ−ＯＣＴで取得されるスペクトル干渉信号の位相データを安定化する手段として機能させるプログラム、及びこのプログラムを搭載して成るＳＳ−ＯＣＴに関する。
（ＳＳ−ＯＣＴの位相データを安定化するプログラム）
本発明に係るＳＳ−ＯＣＴのスペクトル干渉信号の位相データを安定化するプログラムは、ＳＳ−ＯＣＴにおけるコンピュータを、スペクトル干渉信号の位相データを安定化させるために、その補正処理を、第１の補正処理手段及び第２の補正処理手段という二段階で処理をさせる手段として機能させる構成を有している。
第１の補正処理手段は、ＳＳ−ＯＣＴの参照光のみを数値的に取り出し、そこからＳＳ−ＯＣＴ光スペクトルの時間相関を計算することで、ＳＳ−ＯＣＴのスペクトル干渉信号の大まかなスペクトルシフト（具体的には、ＳＳ−ＯＣＴのスペクトル干渉信号のジッター）を補正する。
第２の補正処理手段は、被計測物体のバルクモーション（例えば、人の目の動き等であり、計測画像のアーチファクト（artifact（人工産物）：偽の所見）の原因となる。）等によるスペクトルシフトを検出して補正し、第１の補正処理手段よりもさらに高性能なＳＳ−ＯＣＴのスペクトル干渉信号の位相データの補正を行う。
本発明に係るＳＳ−ＯＣＴの位相データを安定化するプログラムを、ＳＳ−ＯＣＴを構成するコンピュータに搭載すれば、従来例において必要であった「位相安定化ミラー」や「複数チャンネルデジタイザ」等の高価で複雑な装置を追加して使用する必要がない。
以下、本発明に係るＳＳ−ＯＣＴのスペクトル干渉信号の位相データを安定化するプログラムについて、その原理とともに説明する。
＜第１の補正処理手段＞（大まかなスペクトルシフトの補正手段）：
複数のＡ−スキャン（Ａ−スキャン方向と直交する方向に移動（要するにＢ−スキャン）して行う複数のＡ−スキャン）によりＳＳ−ＯＣＴのスペクトル干渉信号が得られるが、第１のＡ−スキャンＡ１と第２のスキャンＡ２の２つのＡ−スキャンで得られた２つのスペクトル干渉信号をＳ１（ｊ）及びＳ２（ｊ）とする。
ここで、ｊはスペクトルのサンプリング点、即ちＡ−スキャンにおける深さ方向（光軸方向）の座標であり、波数ｋに比例するとする。なお波長λ、周波数ｆ、角速度ωとすると、ｋω＝λｆであるから、波数ｋは、ｋ＝λｆ／ωであり、またｆ／ω＝２πであるから、ｋ＝２πｆで表される。
そして、２つのスペクトル干渉信号Ｓ１（ｊ）及びＳ２（ｊ）間の相対的なスペクトルシフトを考える。スペクトル干渉信号Ｓ２（ｊ）がジッターによるスペクトルシフトβを受け、その結果、Ｓ２（ｊ）がＳ’２（ｊ）となったとする。すると、Ｓ’２（ｊ）は次の式（１）で表される。
Ｓ’２（ｊ）＝Ｓ２（ｊ）＊ δ（ｊ−β） ・・・（１）
ここで、＊はコンボリューション（畳み込み積分）である。
両辺をフーリエ変換すると、次の式（２）となる。
Ｆ[Ｓ’２（ζ）]＝Ｆ[Ｓ２（ζ）]ｅｘｐ（−ｉ２πζβ／Ｎ） ・・・（２）
ここで、Ｆ[ ]は[ ]内をフーリエ変換した式である。Ｎは、サンプリング点数であり、ζは被計測物体の深さを示す変数で、ζ＝０，１，・・・・Ｎ−１である。ｉは虚数単位を示している。
ところで、前記したとおり、ドップラー計測では、ドップラー信号は２つのＡ−スキャンで取得されるスペクトル干渉信号の位相差によって得られるが、仮に、２つのＡ−スキャンのスペクトル干渉信号の形状が一致する場合は、スペクトル干渉信号のジッターは、ドップラー計測におけるデータ誤差（位相誤差）に直接影響する。
２つのＡ−スキャンによるスペクトル干渉信号Ｓ１（ｊ）及びＳ２（ｊ）が一致する場合（Ｓ１（ｊ）＝Ｓ２（ｊ））は、Ｓ１（ｊ）及びＳ２（ｊ）から得られたドップラー信号は次の式（２’）で表される位相誤差をもつことになる。
Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（２’）
ここでの ＊ は複素共約を示し、Angle[]は[]内の位相を表す式である（Angle[]は位相角を表す）。
この式は、ドップラー信号は、スペクトルシフトβに１次関数的に比例するため、この式で表されるドップラー信号のデータを１次関数でフィットすることで、そのフィット関数の傾きとしてジッターによるスペクトルシフトβを求めることができ、さらに、そのフィット値をこのドップラー信号から減算することで、ジッターの効果を除去することができる、ということを示している。
しかしながら現実的には、次のような要因で誤差が生じ、そのために２つのＡ−スキャン間でスペクトルデータが異なる。
（１）２つのＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）、Ｓ２（ｊ）それぞれについて一般的な誤差がある。
（２）血流などによるドップラーシフトがある。
（３）バルクモーション（人の目のような被検体のぶれや動きであり、アーチファクト（artifact（人工産物）：偽の所見）を生じる）等の要因による誤差がある。
そこで、本発明のプログラムでは、コンピュータが、式（２’）に従って、基準となるＳ１（ｊ）信号については、誤差等含まない参照光（参照光アームから戻ってくる信号）のスペクトルデータを使用してスペクトルシフトβを求め、このスペクトルシフトβによってスペクトル干渉信号の位相データの補正を行なわせる手段として機能させることを特徴とする。
以上が、第１の補正処理手段であるが、これをさらに詳細に説明すると、実際に検出されるスペクトル干渉信号Ｓ１（ｊ）は、次の式（３）で表される。
Ｓ１（ｊ）＝│Ｅｒ（ｊ）│２＋│Ｅｐ１（ｊ）│２＋Ｅｒ（ｊ）Ｅｐ１ ＊（ｊ）＋Ｅｒ ＊（ｊ）Ｅｐ１（ｊ）
・・・（３）
ここでの ＊ は、複素共約を示す。また、Ｅｒ（ｊ）及びＥｐ１（ｊ）は、それぞれＡ−スキャンにおける参照光及び物体光の光波の振幅である。
この式において、第３項のＥｒ（ｊ）Ｅｐ１ ＊（ｊ）及び第４項Ｅｒ ＊（ｊ）Ｅｐ１（ｊ）は光学系の調整により比較的高い搬送周波数を持つように設定されているのに対して、第１項の│Ｅｒ（ｊ）│２の搬送周波数はゼロであるので、数値的な周波数フィルタリングにより第３項、第４項を除去することが可能である。
また、通常のＳＳ−ＯＣＴでは、参照光のＥｒ（ｊ）Ｅｐ１ ＊（ｊ）強度は、物体光の強度より十分大きい。即ち、│Ｅｒ（ｊ）│２≫│Ｅｐ１（ｊ）│２である。
従って、スペクトル干渉信号Ｓ１（ｊ）は、Ｓ’１（ｊ）≒│Ｅｒ（ｊ）│２ と書くことができる。そのフーリエ変換は、次の式（４）ようになる。
Ｆ[Ｓ’１（ｊ）]≒Ｆ[│Ｅｒ（ｊ）│２] ・・・（４）
ここで、Ｆ[ ]は[ ]内をフーリエ変換した式である。
Ｓ’２（ｊ）は、実際のＯＣＴ計測では、次の式（５）で表される。
Ｓ’２（ｊ）＝｛│Ｅｒ（ｊ）│２＋│Ｅｐ２（ｊ）│２＋Ｅｒ（ｊ）Ｅｐ２ ＊（ｊ）
＋Ｅｒ ＊（ｊ）Ｅｐ２（ｊ）｝＊δ（ｊ−β） ・・・（５）
ここでの ＊ は複素共約を示す。また、＊はコンボリューション（畳み込み積分）である。
Ｅｐ２（ｊ）はプローブ光の振幅である。
Ｓ’２（ｊ）はＳ’１（ｊ）の場合と同様な近似を用いて、次の式ように書くことができる。
Ｓ''２（ｊ）≒│Ｅｒ（ｊ）│２ ＊ δ（ｊ−β）
これをフーリエ変換すると次の式（６）となる。
Ｆ[Ｓ''２（ζ）]≒Ｆ[│Ｅｒ（ｊ）│２ ]ｅｘｐ（−ｉ２πζβ／Ｎ） ・・・（６）
従って位相差は、次の式（７）で表される。
Ｆ[Ｓ''２（ζ）]Ｆ[Ｓ'１ ＊（ζ）]≒Ｉ（ζ）ｅｘｐ（−ｉ２πζβ／Ｎ） ・・・（７）
ここで、Ｉ（ζ）＝Ｆ[│Ｅｒ（ｊ）│２ ]は、ＯＣＴ強度信号である。
理論的には、式（７）の誤差成分の位相−ｉ２πζβ／Ｎは、ζの１次関数になるはずであるが、実際には線形な傾きからのずれが生じる。この位相の勾配βは、次の式（８）で表される誤差を最小となるように決めるとよい。
Ｒ２ １＝ΣζＩ（ζ）[−２πζβ／Ｎ−φ（ζ）]２ ・・・（８）
ここで、φ（ζ）＝∠Ｆ[Ｓ''２（ζ）]Ｆ[Ｓ'１ ＊（ζ）]である。
強度Ｉ（ζ）がノイズレベルより低い場合は、Ｉ（ζ）＝０とする。これ（前記段落００５３〜００６３の記載内容）によりジッターによって発生する線形な誤差成分を数値的に補正することができる。
発明者らは、本発明に係るプログラムをＳＳ−ＯＣＴのコンピュータに搭載し第１の補正処理手段を使用した場合と、搭載せずに第１の補正処理手段を使用しない場合について、それぞれスペクトル干渉信号のスペクトルシフトを算出し、本発明に係るプログラム及び該プログラムを搭載したＳＳ−ＯＣＴの効果を実証するための実験を行った。
この実験では、互いに隣接するＡ−スキャンによってそれぞれ取得したスペクトル干渉信号間について、異なる深さ方向における位相差の標準偏差を算出した。
この実験の結果を図２に示す。図２において、横軸はＡ−スキャン方向の被計測物体の深さ（Depth、単位はmm）を示し、縦軸はスペクトル干渉信号の位相差の標準偏差（σφ）を示し、単位は度（degree）である。そして、○印は第１の補正処理手段を使用した場合の算出値を示し、×印は第１の補正処理手段を使用しない場合の算出値を示している。
この実験の結果によると、第１の補正処理手段を使用の有無にかかわらず、Ａ−スキャン方向の被計測物体の深さによって位相差の標準偏差が異なり、より深い領域においては、位相差の標準偏差が大きくなる。
そして、第１の補正処理手段を使用した場合は、使用しない場合より位相差の標準偏差は小さい。従って、第１の補正処理手段を使用すると、Ａ−スキャン方向の被計測物体の深さ方向の全領域にわたって、スペクトル干渉信号の位相差の標準偏差が小さくなり位相の安定性が改善されることが実証された。
＜第２の補正処理手段＞（スペクトルシフトの残差とバルクモーションの補正手段）
ドップラーＯＣＴ像である隣接したＡ−スキャンのスペクトル干渉信号間の位相差は、そのＡ−スキャンと隣のＡ−スキャンの複素共役を掛け合わせた信号の位相成分として得られる（前記式（２’）参照）。
参照光と物体光の光路差がゼロ（ゼロディレイ）付近の大きなスペクトルシフトは、第１の補正処理手段で補正可能である。しかし、ゼロディレイから遠い信号はエラー（スペクトルシフトの残差）が残る。また、バルクモーションも補正されていない。
そこで、ＳＳ−ＯＣＴ計測のスペクトル干渉信号は、第１の補正処理手段で補正してからさらに、次の式（９）で表される「ドップラー位相シフトΔφ（ζ）」で補正する。
Δφ（ζ）＝（４πτ／λＣ）ｎｖｚ（ζ）＋（４πτ／λＣ）ｖｂ −２πζβ’／Ｎ
・・・（９）
ここで、τは２つのＡ−スキャン間の時間、λＣ は光源の中心波長、ｎは屈折率、変数ζの関数ｖｚ（ζ）は光軸方向の速度、ｖｂはバルクモーションの速度、β’はジッターによるスペクトルシフトの残差、Ｎはサンプリング点数、ζの関数ｖｚ（ζ）は光軸方向の速度、ζは被計測物体の深さを示す変数をそれぞれ示し、ζ＝０，１，・・・・Ｎ−１である。
この式において、第１項の（４πτ／λＣ）ｎｖｚ（ζ）は、計測したいドップラー信号である。第２項の（４πτ／λＣ）ｖｂは、バルクモーションに起因する信号であり、深さζに依存しない定数である。
第１項の−２πζβ’／Ｎにおける「β’」は、前記のとおりジッターによるスペクトルシフトの残差であるが、第１の補正処理手段による補正しない状態のスペクトルシフトβは、第１の補正処理手段で補正すると、その大部分は、第１の補正処理手段による補正で取り除かれる。そして、その残差β’は小さいので、２πの位相のたたみ込みはここでは生じない。
ドップラー信号ｖｚ（ζ）がない領域（血液等の動きのない領域）では、第１項は無視できるので、式（９）のΔφ（ζ）は、ζの切片が４πτｖｂ／λＣであり、傾きが−２πβ’／Ｎの１次関数となる。
この切片と傾きは、次の式（１０）で表される誤差エネルギーを最小とする傾きａｍと切片ｂｍを反復法によって求めることで決めることができる。
Ｒ２ ２＝ΣζＷｍ（ζ）[Δφ（ζ）−（ａｍζ＋ｂｍ）]２ ・・・（１０）
ここで、ｍ＝０，１，２…… であり、Ｗｍ（ζ）は、反復法で誤差エネルギーを最小化する際の各ζにおける誤差の重要性を表す重み関数である。
反復法の式（１０）において、ｍ＝０、即ち初期値Ｗ０（ζ）は次の式（１１）で示される。
Ｉ（ζ）＞ε２ の場合：Ｗ０（ζ）＝（Ｉ（ζ））１／２
その他の場合 ：Ｗ０（ζ）＝０
・・・ （１１）
ここで、Ｉ（ζ）はOCTの強度信号，εはノイズレベルである。
重みの更新は、次の式（１２）で行う。
｜Δφ（ζ）−（ａｍ−１ζ＋ｂｍ−１）｜≧π／２ｍの場合： Ｗｍ（ζ）＝０
その他の場合： Ｗｍ（ζ）＝Ｗｍ−１
・・・（１２）
ここで、ｍ＝１，２，…である。つまり、誤差エネルギーが所定の閾値、例えば、π／２ｍより小さい場合は、重みをゼロとする。それ以外の場合は重みは更新されない、ということである。閾値は繰り返し回数ｍが大きくなると低くなる。ｍの最大値は（たとえば６などに）予め決めておく。
以上、本発明に係る波長走査型光干渉断層計及びその位相安定化プログラムを実施するための最良の形態を実施例に基づいて説明したが、本発明はこのような実施例に限定されることなく、特許請求の範囲記載の技術的事項の範囲内で、いろいろな実施例があることは言うまでもない。
本発明に係る波長走査型光干渉断層計及びその位相安定化プログラムは、眼底血流の検査に有用であり、これにより、緑内障、糖尿病網膜症の初期検出が可能である。また、腫瘍イメージングにも有用であり、悪性腫瘍における新生血管の三次元可視化が可能である。その他、消化器OCTなどへの応用、及び、動物実験への応用が期待される。
１ 波長走査型ＯＣＴ
２ 波長走査型光源
３、５、７、１４ ファイバー
４ ファイバーカップラー
６ 被計測物体
８ 固定参照鏡
１０ 走査鏡
１１、１２、１３ レンズ
１５ 光検知器
１７ ディスプレー

Claims (4)

  1. ＳＳ−ＯＣＴにおいて、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光とが重ねられて光検知器で検出されたスペクトル干渉信号に基づき、断層画像を生成するコンピュータを、スペクトル干渉信号の位相データを補正する補正手段として機能させるＳＳ−ＯＣＴのスペクトル干渉信号の位相を安定化するプログラムであって、
    前記補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式に従い、かつ基準となるスペクトル干渉信号をＳ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴの位相を安定化するプログラム。

    Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ
    但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、＊ は複素共約を示す。
  2. ＳＳ−ＯＣＴにおいて、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光とが重ねられて光検知器で検出されたスペクトル干渉信号に基づき、断層画像を生成するコンピュータを、スペクトル干渉信号の位相データを補正する第１の補正手段と第２の補正手段として順次機能させるＳＳ−ＯＣＴのスペクトル干渉信号の位相を安定化するプログラムであって、
    前記第１の補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（１）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行い、
    前記第２の補正手段は、前記２つのＡ−スキャンで得られたスペクトル干渉信号に係るドップラー位相シフトを示す下記の式（２）において、ドップラー信号ｖｚ（ζ）がない領域における、ζの切片ｂｍ＝４πτｖｂ／λＣ、傾きａｍ＝−２πβ’／Ｎの１次関数について、下記の式（３）に示す誤差エネルギーを最小とする傾きａｍと切片ｂｍを求めることで、バルクモーション及びジッターによるスペクトルシフトβの残差β’を求めて、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴの位相を安定化するプログラム。

    Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（１）
    但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、 ＊ は複素共約を示す。

    Δφ（ζ）＝（４πτ／λＣ）ｎｖｚ（ζ）＋（４πτ／λＣ）ｖｂ −２πζβ’／Ｎ
    ・・・（２）
    但し、τは２つのＡ−スキャン間の時間、λＣ は光源の中心波長、ｎは屈折率、変数ζの関数ｖｚ（ζ）は光軸方向の速度、ｖｂはバルクモーションの速度、β’はジッターによるスペクトルシフトの残差、Ｎはサンプリング点数、ζの関数ｖｚ（ζ）は光軸方向の速度、ζは被計測物体の深さを示す変数をそれぞれ示し、ζ＝０，１，・・・・Ｎ−１である。

    Ｒ２ ２＝ΣζＷｍ（ζ）[Δφ（ζ）−（ａｍζ＋ｂｍ）]２ ・・・（３）
    但し、ｍ＝０，１，２…… であり、ζは被計測物体の深さを示す変数であり、Ｗｍ（ζ）は重みである。
  3. 波長走査型光源と、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光が重ねられたスペクトル干渉信号を検出する光検知器と、光検知器検出されたスペクトル干渉信号に基づき断層画像を生成するコンピュータと、を備えたＳＳ−ＯＣＴであって、
    前記コンピュータは、スペクトル干渉信号の位相データを補正してＳＳ−ＯＣＴの位相を安定化する補正手段を備えており、
    前記補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴ。

    Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ
    但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、 ＊ は複素共約を示す。
  4. 波長走査型光源と、波長走査型光源から出射されて分割され、固定参照鏡で反射された参照光と、被計測物体で反射された物体光が重ねられたスペクトル干渉信号を検出する光検知器と、光検知器検出されたスペクトル干渉信号に基づき断層画像を生成するコンピュータと、を備えたＳＳ−ＯＣＴであって、
    前記コンピュータは、スペクトル干渉信号の位相データを順次補正してＳＳ−ＯＣＴの位相を安定化する第１の補正手段と第２の補正手段を備えており、
    前記第１の補正手段は、隣接する２つのＡ−スキャンにおいて、基準となる第１のＡ−スキャンのスペクトル干渉信号Ｓ１（ｊ）と、第２のＡ−スキャンのスペクトル干渉信号であってジッターによるスペクトルシフトβを受けたスペクトル干渉信号Ｓ’２（ｊ）とから、下記の式（１）に従い、かつ基準となるスペクトル干渉信号Ｓ１（ｊ）として第１のＡ−スキャンの参照光のスペクトルを使用し、スペクトルシフトβを求め、このスペクトルシフトβを消去することにより、スペクトル干渉信号の位相データの補正を行い、
    前記第２の補正手段は、前記２つのＡ−スキャンで得られたスペクトル干渉信号に係るドップラー位相シフトを示す下記の式（２）において、ドップラー信号ｖｚ（ζ）がない領域における、ζの切片ｂｍ＝４πτｖｂ／λＣ、傾きａｍ＝−２πβ’／Ｎの１次関数について、下記の式（３）に示す誤差エネルギーを最小とする傾きａｍと切片ｂｍを求めることで、バルクモーション及びジッターによるスペクトルシフトβの残差β’を求めて、スペクトル干渉信号の位相データの補正を行うことを特徴とするＳＳ−ＯＣＴ。

    Angle[Ｆ[Ｓ’２（ζ）]Ｆ[Ｓ１ ＊（ζ）]]＝−ｉ２πζβ／Ｎ ・・・（１）
    但し、Ｎはサンプリング点数であり、ζは被計測物体の深さを示す変数でζ＝０，１，・・・・Ｎ−１であり、Ｆ[ ]は[ ]内をフーリエ変換した式であり、Angle[]は[]内の位相を表す式であり、 ＊ は複素共約を示す。

    Δφ（ζ）＝（４πτ／λＣ）ｎｖｚ（ζ）＋（４πτ／λＣ）ｖｂ −２πζβ’／Ｎ
    ・・・（２）
    但し、τは２つのＡ−スキャン間の時間、λＣ は光源の中心波長、ｎは屈折率、変数ζの関数ｖｚ（ζ）は光軸方向の速度、ｖｂはバルクモーションの速度、β’はジッターによるスペクトルシフトの残差、Ｎはサンプリング点数、ζの関数ｖｚ（ζ）は光軸方向の速度、ζは被計測物体の深さを示す変数をそれぞれ示し、ζ＝０，１，・・・・Ｎ−１である。

    Ｒ２ ２＝ΣζＷｍ（ζ）[Δφ（ζ）−（ａｍζ＋ｂｍ）]２ ・・・（３）
    但し、ｍ＝０，１，２…… であり、ζは被計測物体の深さを示す変数であり、Ｗｍ（ζ）は重みである。
JP2012000313A 2012-01-05 2012-01-05 波長走査型光干渉断層計及びその位相安定化プログラム Active JP5900950B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2012000313A JP5900950B2 (ja) 2012-01-05 2012-01-05 波長走査型光干渉断層計及びその位相安定化プログラム
PCT/JP2012/082460 WO2013103080A1 (ja) 2012-01-05 2012-12-14 波長走査型光干渉断層計及びその位相安定化方法
US14/368,456 US9354038B2 (en) 2012-01-05 2012-12-14 Swept source optical coherence tomography and method for stabilizing phase thereof
EP12864429.1A EP2801814B1 (en) 2012-01-05 2012-12-14 Swept source optical coherence tomograph and method for stabilizing phase thereof
CN201280062822.1A CN104053980B (zh) 2012-01-05 2012-12-14 波长扫描型光学相干断层仪及其相位稳定化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012000313A JP5900950B2 (ja) 2012-01-05 2012-01-05 波長走査型光干渉断層計及びその位相安定化プログラム

Publications (2)

Publication Number Publication Date
JP2013140077A true JP2013140077A (ja) 2013-07-18
JP5900950B2 JP5900950B2 (ja) 2016-04-06

Family

ID=48745143

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012000313A Active JP5900950B2 (ja) 2012-01-05 2012-01-05 波長走査型光干渉断層計及びその位相安定化プログラム

Country Status (5)

Country Link
US (1) US9354038B2 (ja)
EP (1) EP2801814B1 (ja)
JP (1) JP5900950B2 (ja)
CN (1) CN104053980B (ja)
WO (1) WO2013103080A1 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015070919A (ja) * 2013-10-02 2015-04-16 株式会社トプコン 光画像計測装置
JP2015180863A (ja) * 2014-01-31 2015-10-15 ザ ジェネラル ホスピタル コーポレイション 光干渉法において振幅および位相を使用した多次元速度測定を実施する方法および装置
JP2016106652A (ja) * 2014-12-02 2016-06-20 株式会社ニデック 光コヒーレンストモグラフィ装置、及びそれに用いる制御プログラム
KR101640348B1 (ko) * 2015-08-17 2016-07-18 세종대학교산학협력단 초정밀 광 스캐닝 장치
JP2017185342A (ja) * 2017-07-18 2017-10-12 株式会社トプコン 光画像計測装置
CN112203579A (zh) * 2018-04-02 2021-01-08 俄勒冈健康与科学大学 基于相位的功能光学相干层析成像中的体运动补偿的系统和方法
WO2021095852A1 (ja) * 2019-11-14 2021-05-20 国立大学法人筑波大学 信号処理装置、信号処理方法および信号処理プログラム

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9207062B2 (en) * 2012-12-10 2015-12-08 The Johns Hopkins University Distortion corrected optical coherence tomography system
CN105939652B (zh) * 2014-02-04 2018-07-24 南加利福尼亚大学 具有相位敏感的b扫描配准的光学相干断层扫描(oct)系统
CN104168823B (zh) * 2014-02-28 2017-01-18 深圳市斯尔顿科技有限公司 一种血液流速的测量装置及方法
WO2016099923A1 (en) 2014-12-17 2016-06-23 Pgs Geophysical As Optical filter
MX2017008068A (es) * 2014-12-17 2017-09-28 Pgs Geophysical As Interferometro insensible a la presion.
CN104622474B (zh) * 2015-03-09 2016-11-23 莆田市荔城区聚慧科技咨询有限公司 一种表面深度测量装置
CN104849221B (zh) * 2015-05-21 2017-04-26 清华大学深圳研究生院 一种基于光学相干层析成像技术的相位定标方法
US10129131B2 (en) * 2015-06-30 2018-11-13 Dell Products, Lp System and method for device optimization in a network of devices with embedded electronics
CN105342558B (zh) * 2015-09-30 2017-11-14 苏州大学 一种基于光学相干断层扫描成像中相位误差的校正方法
EP3558091A4 (en) 2016-12-21 2020-12-02 Acucela, Inc. MINIATURIZED AFFORDABLE OPTICAL COHERENCE TOMOGRAPHY SYSTEM FOR OPHTHALMIC APPLICATIONS IN THE HOME
US20180360655A1 (en) * 2017-06-16 2018-12-20 Michael S. Berlin Methods and systems for oct guided glaucoma surgery
JP7144822B2 (ja) * 2017-12-22 2022-09-30 株式会社トーメーコーポレーション 光断層画像撮影装置
EP3809948A4 (en) 2018-06-20 2022-03-16 Acucela Inc. MINIATURIZED MOBILE, LOW COST OPTICAL COHERENCE TOMOGRAPHY SYSTEM FOR HOME OPHTHALMIC APPLICATIONS
CN108956533A (zh) * 2018-07-13 2018-12-07 福州大学 一种可消除多色误差的oct共轭镜像去除装置及方法
US11730363B2 (en) 2019-12-26 2023-08-22 Acucela Inc. Optical coherence tomography patient alignment system for home based ophthalmic applications
US10959613B1 (en) 2020-08-04 2021-03-30 Acucela Inc. Scan pattern and signal processing for optical coherence tomography
AU2021324968A1 (en) 2020-08-14 2023-03-02 Acucela Inc. System and method for optical coherence tomography a-scan decurving
CN111829954B (zh) * 2020-09-09 2023-07-25 广东工业大学 一种提高全场扫频光学相干层析测量量程的系统及方法
US11393094B2 (en) 2020-09-11 2022-07-19 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
CA3192083A1 (en) 2020-09-30 2022-04-07 Acucela Inc. Myopia prediction, diagnosis, planning, and monitoring device
US11974806B2 (en) * 2021-02-19 2024-05-07 Topcon Corporation Ophthalmologic information processing apparatus, ophthalmologic apparatus, ophthalmologic information processing method, and recording medium
US11497396B2 (en) 2021-03-24 2022-11-15 Acucela Inc. Axial length measurement monitor

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008175698A (ja) * 2007-01-18 2008-07-31 Univ Of Tsukuba 光コヒーレンストモグラフィーの画像処理方法及び画像処理装置
JP2008545970A (ja) * 2005-06-01 2008-12-18 ザ ジェネラル ホスピタル コーポレイション 位相分解光学周波数領域画像化を行うための装置、方法及びシステム
JP2010014514A (ja) * 2008-07-03 2010-01-21 Fujifilm Corp 光断層画像化装置及び光断層画像化装置における干渉信号の処理方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5501226A (en) * 1994-10-19 1996-03-26 Carl Zeiss, Inc. Short coherence length, doppler velocimetry system
DE19814057B4 (de) 1998-03-30 2009-01-02 Carl Zeiss Meditec Ag Anordnung zur optischen Kohärenztomographie und Kohärenztopographie
JP2002310897A (ja) 2001-04-13 2002-10-23 Japan Science & Technology Corp 光コヒーレンストモグラフィーにおける透光体の動きによる高速光遅延発生方法及びその高速光遅延発生装置
JP4045140B2 (ja) 2002-06-21 2008-02-13 国立大学法人 筑波大学 偏光感受型光スペクトル干渉コヒーレンストモグラフィー装置及び該装置による試料内部の偏光情報の測定方法
CN103082996A (zh) * 2003-10-27 2013-05-08 通用医疗公司 用于使用频域干涉测量法进行光学成像的方法和设备
WO2006107755A2 (en) * 2005-04-01 2006-10-12 Visualsonics Inc. System and method for 3-d visualization of vascular structures using ultrasound
CN100512752C (zh) * 2007-01-12 2009-07-15 天津市先石光学技术有限公司 光程可选择的无创人体成分检测方法及检测装置
CN101720205B (zh) * 2007-04-10 2012-11-28 南加利福尼亚大学 使用多普勒光学相干层析成像进行血液流量测量的方法和系统
US8133127B1 (en) * 2008-07-21 2012-03-13 Synder Terrance W Sports training device and methods of use
JP5725697B2 (ja) * 2009-05-11 2015-05-27 キヤノン株式会社 情報処理装置および情報処理方法
US9192294B2 (en) * 2012-05-10 2015-11-24 Carl Zeiss Meditec, Inc. Systems and methods for faster optical coherence tomography acquisition and processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008545970A (ja) * 2005-06-01 2008-12-18 ザ ジェネラル ホスピタル コーポレイション 位相分解光学周波数領域画像化を行うための装置、方法及びシステム
JP2008175698A (ja) * 2007-01-18 2008-07-31 Univ Of Tsukuba 光コヒーレンストモグラフィーの画像処理方法及び画像処理装置
JP2010014514A (ja) * 2008-07-03 2010-01-21 Fujifilm Corp 光断層画像化装置及び光断層画像化装置における干渉信号の処理方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HONG Y J , ET.AL.: "Non-Invasive Three-Dimensional Angiography of Human Eye by Doppler Optical Coherence Tomography", IQEC/CLEO PACIFIC RIM 2011, JPN6015049242, 28 August 2011 (2011-08-28), pages 769 - 771, XP032175620, ISSN: 0003255494, DOI: 10.1109/IQEC-CLEO.2011.6194122 *
HONG,Y.-J. ET.AL.: ""High-penetration swept source Doppler optical coherence angiography by fully numerical phase stabi", OPTICS EXPRESS, vol. Volume 20, Issue 3, JPN6013000879, 23 January 2012 (2012-01-23), pages 2740 - 2760, XP055080126, ISSN: 0003255493, DOI: 10.1364/OE.20.002740 *
TSAI,M.-T. ET.AL.: ""Microvascular Imaging Using Swept-Source Optical Coherence Tomography with Single-Channel Acquisit", APPLIED PHYSICS EXPRESS, vol. Volume 4, Number 9, JPN6013000881, 12 August 2011 (2011-08-12), pages 097001 - 1, ISSN: 0003255492 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015070919A (ja) * 2013-10-02 2015-04-16 株式会社トプコン 光画像計測装置
JP2015180863A (ja) * 2014-01-31 2015-10-15 ザ ジェネラル ホスピタル コーポレイション 光干渉法において振幅および位相を使用した多次元速度測定を実施する方法および装置
JP2016106652A (ja) * 2014-12-02 2016-06-20 株式会社ニデック 光コヒーレンストモグラフィ装置、及びそれに用いる制御プログラム
KR101640348B1 (ko) * 2015-08-17 2016-07-18 세종대학교산학협력단 초정밀 광 스캐닝 장치
JP2017185342A (ja) * 2017-07-18 2017-10-12 株式会社トプコン 光画像計測装置
CN112203579A (zh) * 2018-04-02 2021-01-08 俄勒冈健康与科学大学 基于相位的功能光学相干层析成像中的体运动补偿的系统和方法
JP2021520244A (ja) * 2018-04-02 2021-08-19 オレゴン ヘルス アンド サイエンス ユニバーシティ 位相ベースの機能光コヒーレンストモグラフィにおけるバルク運動補償のシステム及び方法
US11944382B2 (en) 2018-04-02 2024-04-02 Oregon Health & Science University Systems and methods for bulk motion compensation in phase-based functional optical coherence tomograpgy
WO2021095852A1 (ja) * 2019-11-14 2021-05-20 国立大学法人筑波大学 信号処理装置、信号処理方法および信号処理プログラム

Also Published As

Publication number Publication date
WO2013103080A1 (ja) 2013-07-11
EP2801814A1 (en) 2014-11-12
CN104053980B (zh) 2016-05-18
US20150009507A1 (en) 2015-01-08
JP5900950B2 (ja) 2016-04-06
EP2801814A4 (en) 2015-09-02
EP2801814B1 (en) 2018-04-18
CN104053980A (zh) 2014-09-17
US9354038B2 (en) 2016-05-31

Similar Documents

Publication Publication Date Title
JP5900950B2 (ja) 波長走査型光干渉断層計及びその位相安定化プログラム
US9706915B2 (en) Method of motion correction in optical coherence tomography imaging
US7602501B2 (en) Interferometric synthetic aperture microscopy
JP5626687B2 (ja) ２ビーム型光コヒーレンストモグラフィー装置
JP5787255B2 (ja) Ｐｓ−ｏｃｔの計測データを補正するプログラム及び該プログラムを搭載したｐｓ−ｏｃｔシステム
JP6346410B2 (ja) ジョーンズマトリックスｏｃｔシステム及び該ｏｃｔで得られた計測データを画像処理するプログラム
JP5864258B2 (ja) スペクトル光コヒーレンス断層映像法を用いて構造データを収集する方法および装置
JP2008039651A (ja) 光断層画像の処理方法
JP2007127425A (ja) 光断層画像化法における補正方法
WO2012107307A1 (en) Optimized device for swept source optical coherence domain reflectometry and tomography
JP2016209182A (ja) 撮像装置、撮像装置の作動方法、情報処理装置、及び情報処理装置の作動方法
Subhash et al. Microcirculation imaging based on full-range high-speed spectral domain correlation mapping optical coherence tomography
McNamara et al. Development of a first-generation miniature multiple reference optical coherence tomography imaging device
Ksenofontov et al. Numerical method for axial motion artifact correction in retinal spectral-domain optical coherence tomography
JP6606640B2 (ja) 眼科装置及びその制御方法
JP6579718B2 (ja) ジョーンズマトリックスｏｃｔ装置及びプログラム
Gelikonov et al. Numerical method for axial motion correction in optical coherence tomography
JP6779662B2 (ja) 撮像装置、撮像装置の制御方法、及びプログラム
US20230358527A1 (en) Optical interference tomographic imaging device
Romania Optical coherence tomography in oncological imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141222

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151208

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160126

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: 20160223

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160303

R150 Certificate of patent or registration of utility model

Ref document number: 5900950

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250