JP4389032B2 - 光コヒーレンストモグラフィーの画像処理装置 - Google Patents

光コヒーレンストモグラフィーの画像処理装置 Download PDF

Info

Publication number
JP4389032B2
JP4389032B2 JP2007009574A JP2007009574A JP4389032B2 JP 4389032 B2 JP4389032 B2 JP 4389032B2 JP 2007009574 A JP2007009574 A JP 2007009574A JP 2007009574 A JP2007009574 A JP 2007009574A JP 4389032 B2 JP4389032 B2 JP 4389032B2
Authority
JP
Japan
Prior art keywords
oct
image
interference signal
scan
envelope
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
JP2007009574A
Other languages
English (en)
Other versions
JP2008175698A (ja
Inventor
嘉晃 安野
豊彦 谷田貝
修一 巻田
瑛周 洪
Original Assignee
国立大学法人 筑波大学
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 国立大学法人 筑波大学 filed Critical 国立大学法人 筑波大学
Priority to JP2007009574A priority Critical patent/JP4389032B2/ja
Publication of JP2008175698A publication Critical patent/JP2008175698A/ja
Application granted granted Critical
Publication of JP4389032B2 publication Critical patent/JP4389032B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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]

Landscapes

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

Description

本発明は、光コヒーレンストモグラフィー（ＯＣＴ）の画像処理装置に関し、特に、ＯＣＴの光源のスペクトルに由来するＯＣＴ干渉信号（スペクトル干渉信号）におけるノイズを補正する画像処理装置に関する。
より詳細には、本発明は、半導体製品の製造等の工業分野、動物、昆虫の生体観察や植物の構造観察のような動植物の研究・観察分野、各種文化財の解析・鑑定技術の分野、ロボット技術分野（ここでは、植物、昆虫、動物、人等の各部器官を観察し、その構造・機能をロボットの機能として応用するための技術）、医療用検査装置等、各種の技術分野で利用される光コヒーレンストモグラフィー（ＯＣＴ）装置で取得した画像データの画像処理装置に関する。
従来、各種被計測物体の非破壊断層計測技術の１つとして、光断層画像装置「光コヒーレンストモグラフィー」（ＯＣＴ）がある（特許文献１参照）。ＯＣＴは、光を計測プローブとして用いるため、被計測物体の屈折率分布、分光情報、偏光情報（複屈折率分布）等が計測できるという利点がある。
基本的なＯＣＴ４３は、マイケルソン干渉計を基本としており、その原理を図６で説明する。光源４４から射出された光は、コリメートレンズ４５で平行化された後に、ビームスプリッター４６により参照光と物体光に分割される。物体光は、物体アーム内の対物レンズ４７によって被計測物体４８に集光され、そこで散乱・反射された後に再び対物レンズ４７、ビームスプリッター４６に戻る。
一方、参照光は参照アーム内の対物レンズ４９を通過した後に参照鏡５０によって反射され、再び対物レンズ４９を通してビームスプリッター４６に戻る。このようにビームスプリッター４６に戻った物体光と参照光は、物体光とともに集光レンズ５１に入射し光検出器５２（フォトダイオード等）に集光される。
ＯＣＴの光源４４は、時間的に低コヒーレンスな光（異なった時刻に光源から出た光同士は極めて干渉しにくい光）の光源を利用する。時間的低コヒーレンス光を光源としたマイケルソン型の干渉計では、参照アームと物体アームの距離がほぼ等しいときにのみ干渉信号が現れる。この結果、参照アームと物体アームの光路長差（τ）を変化させながら、光検出器５２で干渉信号の強度を計測すると、光路長差に対する干渉信号（インターフェログラム）が得られる。
そのインターフェログラムの形状が、被計測物体４８の奥行き方向の反射率分布を示しており、１次元の軸方向走査により被計測物体４８の奥行き方向の構造を得ることができる。このように、ＯＣＴ４３では、光路長走査により、被計測物体４８の奥行き方向の構造を計測できる。
このような軸方向の走査のほかに、横方向の機械的走査を加え、２次元の走査を行うことで被計測物体の２次元断面画像が得られる。この横方向の走査を行う走査装置としては、被計測物体を直接移動させる構成、物体は固定したままで対物レンズをシフトさせる構成、被計測物体も対物レンズも固定したままで、対物レンズの瞳面付近においたガルバノミラーの角度を回転させる構成等が用いられている。
以上の基本的なＯＣＴが発展したものとして、分光器を用いてスペクトル信号を得るスペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）と、光源の波長を走査してスペクトル干渉信号を得る波長走査型ＯＣＴ（Swept Source OCT、略して「ＳＳ−ＯＣＴ」という。）がある。ＳＤ−ＯＣＴには、フーリエドメインＯＣＴ（Fourier Domain OCT、略して「ＦＤ−ＯＣＴ」という。特許文献２参照）、及び偏光感受型ＯＣＴ（Polarization-Sensitive OCT、略して「ＰＳ−ＯＣＴ」という。特許文献３参照）がある。
フーリエドメインＯＣＴは、被計測物体からの反射光の波長スペクトルを、スペクトロメーター（スペクトル分光器）で取得し、このスペクトル強度分布に対してフーリエ変換することで、実空間（ＯＣＴ信号空間）上での信号を取り出すことを特徴とするものであり、このフーリエドメインＯＣＴは、奥行き方向の走査を行う必要がなく、ｘ軸方向の走査を行うことで被計測物体の断面構造を計測可能である。
波長走査型ＯＣＴは、高速波長スキャニングレーザーにより光源の波長を変え、スペクトル信号と同期取得された光源走査信号を用いて干渉信号を最配列し、信号処理を加えることで３次元光断層画像を得るものである。なお、光源の波長を変える手段として、モノクロメーターを利用したものでも、波長走査型ＯＣＴとして利用可能である。
特開２００２−３１０８９７号公報 特開平１１−３２５８４９号公報 特開２００４−０２８９７０号公報
ところで、ＯＣＴの光源のスペクトルは、一般的に、きれいなガウシアンエンベロープ（正規分布）ではなく、形状が不定で、多くのリップルのピークを持っている。これはＯＣＴ画像の深さ方向の分解能の低下をもたらす。
また、光源のスペクトル分布がガウシアンエンベロープの場合は（図７（ａ）参照）、光学系の性能を評価する点応答関数（ＰＳＦ）は対称であり１山である（図７（ｂ）参照）。しかし、光源のスペクトル分布が非対称の場合（図７（ｃ）参照）は、複数の山ができたり（図７（ｄ）参照）、サイドローブ（山の裾野の方の持ち上がり）の原因となり、画像の劣化およびクロストークの原因となる。
本発明は、このような問題点を解決することを目的とするものであり、ＯＣＴ干渉信号からＯＣＴの光源のスペクトル分布及びそのエンベロープを求め、このエンベロープから被計測物体の計測で得られたＯＣＴ干渉信号を補正可能なフィルターを作成し、ＯＣＴ画像の深さ方向の分解能の低下を改善する画像処理装置を実現することを課題とする。
本発明は上記課題を解決するために、光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、光源のスペクトルに由来するＯＣＴ画像の深さ方向の分解能の低下を改善することを特徴とする光コヒーレンストモグラフィーの画像処理装置を提供する。
本発明は上記課題を解決するために、光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、該エンベロープ補正手段で補正して得られたＯＣＴ画像の情報エントロピーを評価関数とし、画像の位相を操作しながら評価関数の値をもとめ、その値を最小化するように、位相を操作するというステップを繰り返し、情報エントロピーの値が最小となった画素それぞれの位相を最適位相とし、ＯＣＴ画像の位相を最適化することによって、ＯＣＴ画像の鮮明度及び分解能を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置を提供する。
本発明は上記課題を解決するために、光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化するＯＣＴ干渉信号平滑化手段と、該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、を備えており、ＯＣＴ干渉信号（スペクトル干渉信号）からスペクトルエンベロープを差し引き、干渉成分のみをぬきだし、干渉成分のみから成るＯＣＴ干渉信号を得て、２次元画像のプレビューの高速化をすることを特徴とする光コヒーレンストモグラフィーの画像処理装置を提供する。
本発明は上記課題を解決するために、光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化し、又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、該エンベロープ補正手段で補正して得られたＯＣＴ画像の、光源のコヒーレンスに起因するスペックルを除去するために、注目している画素における光強度の値を、その画素の周囲の画素における光強度の値のうち、もっとも大きな値で置き換えることにより、画像のスペックルを除去し、ＯＣＴ画像の出力を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置を提供する。
本発明は上記課題を解決するために、光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、該エンベロープ補正手段で補正して得られたＯＣＴ画像を高コントラスト化するために、そのＯＣＴ画像の、位相情報を無視し、強度情報のみを用い、被計測物体の散乱光強度に対応した画素を抽出し、３次元マスクを作成する手段と、を備えており、そのマスクがかかっている部分のＯＣＴ画像を抽出することにより、ＯＣＴ画像から特定の構造物を抽出して、ＯＣＴ画像の位相を最適化することによって、ＯＣＴ画像の鮮明度及び分解能を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置を提供する。
以上のような構成から成る本発明に係る光コヒーレンストモグラフィーの画像処理装置によれば、ＯＣＴ画像の深さ方向の分解能の低下を改善することができ、画像の劣化を防止できるという効果が生じる。
本発明に係る光コヒーレンストモグラフィー（以下、「ＯＣＴ」と言う。）の画像処理装置を実施するための最良の形態を実施例に基づいて図面を参照して、以下に説明する。
まず、本発明の画像処理装置を実施するための前提となるスペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）としてフーリエドメインＯＣＴ（ＦＤ−ＯＣＴ）を概略説明する。また、波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）について概略説明する。
（ＦＤ−ＯＣＴ）
図１は、ＦＤ−ＯＣＴ１の全体構成を示す図である。広帯域光源２、低コヒーレンス干渉計３、及び分光器４（スペクトロメーター）とを備えている。このＦＤ−ＯＣＴ１は、低コヒーレンス干渉の原理を用いて奥行き方向の分解能を得ているため、光源として、ＳＬＤ（スーパールミネツセントダイオード）や超短パルスレーザー等の広帯域光源２が用いられる。
広帯域光源２から出た光は、まずビームスプリッター５で物体光と参照光に分割される。このうち物体光は、レンズ６を通してガルバノミラー７で反射され被計測物体８（眼底などの生体試料）を照射し、そこで反射、散乱された後に分光器４に導かれる。一方、参照光はレンズ９を通して参照鏡１０（平面鏡）で反射された後に物体光と並行に分光器４に導かれる。これらの二つの光は分光器４の回折格子１１によって同時に分光され、スペクトル領域で干渉し、結果、スペクトル干渉縞がＣＣＤ１２によって計測される。
このスペクトル干渉縞に対して適当な信号処理を行うことで、被計測物体８のある点における深さ方向１次元の屈折率分布の微分、つまり、反射率分布を得ることが可能となる。さらに、被計測物体８上の計測点をガルバノミラー７を駆動し１次元走査することにより２次元断層画像（ＦＤ−ＯＣＴ画像）を得ることができる。
通常のＯＣＴでは、２次元断層画像を得るために、深さ（光軸）方向の走査（この走査を「Ａ−スキャン」と言い、この方向を「Ａ−方向」、「Ａスキャン方向」とも言う。）と、縦方向の操作（この走査を「Ｂ−スキャン」と言い、この方向を「Ｂ−方向」、「Ｂスキャン方向」とも言う。）の２次元の機械的走査が必要なのに対して、ＦＤ−ＯＣＴ１では、Ａ−スキャンは不要で一回の測定で深さ方向の後方散乱データを取得することができるから、Ｂ−スキャンの１次元の機械的走査しか必要とされない。
なお、Ａ−方向とＢ−方向で形成される平面に垂直な方向のスキャンを「Ｃ−スキャン」と言い、この方向を「Ｃ−方向」、「Ｃスキャン方向」とも言う。要するに、ＦＤ−ＯＣＴ１では、面内に２次元走査（Ｂ−スキャンおよびＣ−スキャン）をすることにより高速な断層計測が可能で、被計測物体８内部の２次元および３次元情報を得ることができる。
図２は、本発明の実施例３を適用する波長走査型ＯＣＴ（ＳＤ−ＯＣＴ）２４の全体構成を示す図である。波長走査型光源２５から出射された出力光を、ファイバー２６を通してファイバーカップラー２７に送る。この出力光を、ファイバーカップラー２７において、ファイバー２８を通して被計測物体２９への照射する物体光と、ファイバー３０を通して固定参照鏡３１に照射する参照光に分割する。
物体光は、ファイバー２８、レンズ３２、角度が可変な走査鏡３３及びレンズ３４を介して、被計測物体２９に照射、反射され、同じルートでファイバーカップラー２７に戻る。参照光は、ファイバー３０、レンズ３５及びレンズ３６を介して固定参照鏡３１に照射、反射されて同じルートでファイバーカップラー２７に戻る。
そして、これらの物体光と参照光はファイバーカップラー２７で重ねられ、ファイバー３７を通して光検知器３８（ＰＤ（フォトダイオード）等のポイントセンサが使用される。）に送られ、スペクトル干渉信号として検出され、コンピュータ３９に取り込まれる。光検知器３８における検知出力に基づいて、被計測物体２９の奥行き方向（Ａ方向）と走査鏡の走査方向（Ｂ方向）の断面画像が形成される。コンピュータ３９にはディスプレー４０が接続されている。
ここで、波長走査型光源２５は、時間的に波長を変化させて走査する光源であり、即ち波長が時間依存性を有する光源である。これにより、固定参照鏡３１を走査（移動。Ａ−スキャン）することなく、 被計測物体２９の奥行き方向の反射率分布を得て奥行き方向の構造を取得することができ、１次方向の走査（Ｂ−スキャン）をするだけで、二次元の断層画像を形成することができる。
以上のようなＳＤ−ＯＣＴ（上記ＦＤ−ＯＣＴ等）やＳＤ−ＯＣＴによって得られるＯＣＴ画像は、光源のスペクトルにより影響を受ける。即ち、ＯＣＴの光源のスペクトルは、一般的に、きれいなガウシアンエンベロープ（正規分布）ではなく、形状が不定で、多くのリップルのピークを持っている。これはＯＣＴ画像の深さ方向の分解能の低下をもたらす。
本発明は、このような光源のスペクトルに由来するＯＣＴ画像の深さ方向の分解能の低下を補正することで、特に被計測物体の奥行き方向の分解能の低下による画像の劣化を防止することを特徴とする画像処理装置である。
（原理等）
本発明に係る画像処理装置の原理を図３により説明する。光源のスペクトル分布は、一例を図３（ａ）に波長（λ）と光強度（Ｉ）の関係で示すように、ガウシアンエンベロープ（正規分布）ではない。
この光源を使用したＯＣＴで、被計測物体をＢ方向（被計測物体の奥行き方向であるＡスキャン方向に直角な方向）にスキャンしたときのＯＣＴ干渉信号を、図３（ｂ）に示す。図３（ｂ）において、横軸はＡスキャン方向を示し、縦軸は光強度（Ｉ）を示している。なお、ＯＣＴ信号を平均化することで、被測定物体の特性はなくなるので、被計測物体としては、特に校正用試料を用いなくても光源の特性のみ残る。但し、校正用として一定の被計測体を決めておいて、それを標準試料として用いてもよい。
この図３（ｂ）に示すように、ＯＣＴ干渉信号は、主に光源のスペクトル分布、さらにＯＣＴの光学系（ＯＣＴを構成するファイバー、レンズ等の機器）の特性の影響を受け、対称状態（ガウシアン等の正規分布状態）とならない。そこで、本発明では、図３（ｂ）に示すＯＣＴ干渉信号について、エンベロープ（包絡線）信号を検出する。
Ｂスキャン方向にわたった複数の点におけるＯＣＴ干渉信号（スペクトル干渉信号）（図３（ｂ）参照）を合計して平均をとって得られた平滑化されたエンベロープ信号を図３（ｃ）に示す。図３（ｃ）において、横軸は、Ａスキャン方向（奥行き方向）を示しているが、これは、スペクトルの波長方向に対応している。縦軸は光強度（Ｉ）を示している。このエンベロープ信号は、光源のスペクトルだけでなく、ＯＣＴの光学系の特性を含んだスペクトル特性を示すものである。
そして、本発明では、図３（ｃ）で示すエンベロープ信号が、図３（ｄ）に示すようなガウシアン分布となるようなフィルター（このようなフィルターを「逆フィルター」という。）を作成する。この逆フィルター作成は、図３（ｃ）に示すエンベロープ信号と、図３（ｄ）に示すガウシアン分布に基づき、作成する。例えば、逆フィルターは、図３（ｄ）を図３（ｃ）で割り算したもの、あるいはウィナーフィルターを用いる。
このような逆フィルターにより、被計測物体の計測で得られたＯＣＴ干渉信号を補正すれば、Ｂスキャン方向の各点で得られるＯＣＴ干渉信号を合計して平均した平滑されたＯＣＴ干渉信号は、例えば、図３（ｅ）（横軸はＡスキャン方向、縦軸は光強度）に示すようになる。このＯＣＴ干渉信号のエンベロープ信号は、ガウシアン分布となる。
なお、図３（ｅ）に示すエンベロープ信号のガウシアン分布の形状は、理論的、理想的には被測定物体に異存しない。しかし、現実的には、例えば、皮膚と眼底では被測定対象までの光路の光吸収率などが異なる。このような違いをＯＣＴの光学系の特性と見なせば、測定対象に応じた補正が必要となる。
図４は、本発明に係る画像処理装置を実施するための画像処理装置１３の基本構成を示す図である。この画像処理装置１３は、ＯＣＴを構成するＣＣＤ光検出器１４に接続されたコンピュータから成る。
この画像処理装置１３は、ＯＣＴ干渉信号入力部１５、ＯＣＴ干渉信号平滑化手段１６、エンベロープ検出手段１７、逆フィルター作成手段１８、エンベロープ補正手段１９、及びＯＣＴ干渉信号出力部２０（補正されたＯＣＴ画像信号の出力部）を備えている。
この画像処理装置１３では、ＯＣＴ干渉信号入力部で入力されたＯＣＴ干渉信号を利用して、ＯＣＴ干渉信号平滑化手段１６、エンベロープ検出手段１７及び逆フィルター作成手段１８により逆フィルターを作成し、これをエンベロープ補正手段１９に送信する。
この画像処理装置１３で、被計測物体を計測し、そのＯＣＴ画像を補正する場合は、ＯＣＴ干渉信号入力部１５で取得したＯＣＴ干渉信号をエンベロープ補正手段１７に送り、ここで逆フィルターで補正してＯＣＴ干渉信号出力部２０から出力する。なお、通常、光コーヒレントトモグラフィーのＯＣＴ干渉信号を補正することがない場合は、ＯＣＴ干渉信号入力部１５から入力したＯＣＴ干渉信号の生信号をＯＣＴ干渉信号出力部２０から出力してＯＣＴ画像を得る。
画像処理装置１３によって、ＯＣＴ干渉信号を補正する画像処理を、以下、順次説明する。
（スペクトルエンベロープ信号の取得）
本発明に係る画像処理装置は、ＯＣＴ干渉信号から、光源のスペクトルエンベロープ信号（「スペクトル包絡線」又は「分光分布特性」ともいう。）を取得することが基本となる。このスペクトルエンベロープ信号に基づいて、逆フィルターを作成したり、後述するプレビューの高速化処理を行う。
従って、まず光源のスペクトルエンベロープ信号の取得について、ＦＤ−ＯＣＴのようなスペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）の場合と、波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）の場合を例に挙げて説明する。
（１）スペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）の場合
スペクトルドメインＯＣＴの場合は、図４において、図３（ａ）に示すようなスペクトルを有する光源を用い、ＣＣＤ光検出器１４で検出したＡスキャン方向（被計測物体の奥行き方向）のＯＣＴ干渉信号を、画像処理装置１３のＯＣＴ干渉信号入力部に入力する。そして、ＯＣＴ干渉信号平滑化手段において、Ａスキャン方向のＯＣＴ干渉信号をＢスキャン方向で合計して平均をとり（アンサンブル平均をとり）、図３（ｂ）に示すような、平滑化されたＯＣＴ干渉信号を取得する。
即ち、ＯＣＴ干渉信号平滑化手段１６は、Ａスキャン方向のＯＣＴ干渉信号を、Ｂスキャン方向の全ての点について合計し、その平均をとり、Ａスキャン方向のＯＣＴ干渉信号成分のばらつきをなくし、平滑化する処理を行う。
なお、Ｂスキャン方向のかわりに、上記同様にＣスキャン方向（Ａスキャン方向及びＢスキャン方向に直交する方向）でアンサンブル平均をとり、ＯＣＴ干渉信号成分を平滑化し、スペクトルエンベロープ信号を抽出してもよい。また、Ｂスキャン方向およびＣスキャン方向でアンサンブル平均をとり、ＯＣＴ干渉信号成分を平滑化し、スペクトルエンベロープ信号を抽出してもよい。
上記のとおり、ＯＣＴ干渉信号平滑化手段１６で平滑化されたＯＣＴ干渉信号（図３（ｂ）参照）は、エンベロープ検出手段１７において、そのエンベロープ（包絡線）を検出する処理が行われ、図３（ｃ）に示すようなスペクトルエンベロープ信号（「スペクトルエンベロープ情報」とも言う。）が形成される。
（２）波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）の場合
波長走引型ＯＣＴの場合も、基本的には、画像処理装置によってスペクトルドメインＯＣＴの場合と同様に、ＯＣＴ干渉信号からスペクトルエンベロープ信号を検出し、このスペクトルエンベロープ信号に基づいて逆フィルターを作成したり、後述するプレビューの高速化処理を行うものである。しかし、波長走引型ＯＣＴの場合は、下記に説明するように、ＯＣＴ干渉信号を取得及び画像処理について、ＳＳ−ＯＣＴとは若干異なる。
波長走引型ＯＣＴの場合は、ＳＳ−ＯＣＴにおいてＯＣＴ干渉信号を取得するに際して、光源の相対強度雑音（ＲＩＮ）の影響を減らすため、光ファイバー及びカプラで２つに分けられた信号強度をそれぞれＣＣＤで検出して、両者の信号強度差を検出（これを「バランス検出」と言う。）する。
この信号強度差がＯＣＴ干渉信号として検出される。そのため、ＯＣＴ信号のエンベロープは、光源のスペクトルのエンベロープと光ファイバーなどの光学系のスペクトル透過特性などが重ね合わされる。
上記取得したＯＣＴ干渉信号を、図５に示す画像処理装置２１のＯＣＴ干渉信号入力部に入力し、スペクトルエンベロープ信号を取得する処理を行うのであるが、波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）の場合は、スペクトルドメインＯＣＴで説明したような単純な平滑化処理では、スペクトルエンベロープ信号に、ファイバーやカプラ等の光学素子のスペクトル特性が乗せられたものが抽出されてしまい、光源のみのスペクトルエンベロープ信号を求めることはできない。
そこで、図４のＯＣＴ干渉信号平滑化手段１６及びエンベロープ検出手段１７に替えて、図５に示すように、干渉信号変換エンベロープ検出手段２２において、下記のイ．ロ．又はハ．のいずれかを用いて、ＯＣＴ干渉信号からスペクトルエンベロープ信号を取得する。
イ．ヒルベルト変換（Hilbert変換）を用いる場合
ヒルベルト変換を用いる場合を図８を参照して説明する。Ｂ又はＣスキャン方向のある１点において、Ａスキャン方向のＯＣＴ干渉信号Ｓ（ｋ）に注目する（図８（ａ）参照）。ここで、ｋは光源のスペクトルの波数に対応している。ここで、ｋは波数であり、波長λの逆数×２π、即ちｋ＝２π／λで表されるので、Ａスキャン方向の深さ方向の位置に対応している。
Ｂ又はＣスキャン方向の上記ある１点において、ｋについてのアンサンブル平均値＜Ｓ（ｋ）＞（Ａスキャン方向の値を合算の平均値。図８（ｂ）参照）を、それぞれのデータから減じ、Ｓ’（ｋ）＝Ｓ（ｋ）−＜Ｓ（ｋ）＞を求める（図８（ｃ）参照）。
ＯＣＴ干渉信号のスペクトルエンベロープ信号は、次の数式１で表すことができる。ここで、ＮはＡスキャン方向のデータ総数、ｊはＡスキャン方向のデータ番号、Ｈはヒルベルト変換である。このスペクトルエンベロープ信号が求まる処理はエンベロープ検出手段で行われる。
Figure 0004389032
ロ．フーリエ変換（fast Fourier transform: FFT変換）を用いる場合
フーリエ変換を用いる場合を図９を参照して説明する。Ｂ又はＣスキャンのある１点において、Ａスキャン方向のＯＣＴ干渉信号Ｓ（ｋ）に注目する（図９（ａ）参照）。この信号をフーリエ変換し（図９（ｂ）参照）、周波数領域の信号を計算する。
この信号において、エンベロープ信号に対応した周波数成分のみを抽出する（図９（ｃ）参照）ような、ローパスフィルター又はバンドパスフィルターをかけ、再び、逆フーリエ変換を行うことで、ＯＣＴ干渉信号のエンベロープ信号を抽出することができる（図９（ｄ）参照）。
ハ．スペクトルを２乗して平均化する場合
より高速な方法として、スペクトルを２乗して平均化する方法がある。Ｓ（ｋ）は包絡線成分Ｅ（ｋ）とＯＣＴ干渉信号成分ｃｏｓ（αｋ）の成分の積で表される。
即ち、Ｓ（ｋ）＝Ｅ（ｋ）ｃｏｓ（αｋ）である。Ｓ（ｋ）２ のアンサンブル平均をとり、ＯＣＴ干渉信号成分を平均化する。
つまり、
＜Ｓ（ｋ）２ ＞＝＜Ｅ（ｋ）２ ｃｏｓ２（αｋ）＞＝Ｅ（ｋ）２＜ ｃｏｓ２（αｋ）＞
＝１／２Ｅ（ｋ）２
となるため、
Ｅ（ｋ）＝（＜Ｓ（ｋ）２ ＞×２）１／２ で、ＯＣＴ干渉信号のエンベロープ信号を求めることができる。
なお、Ｓ（ｋ）２ を求めるとき、あらかじめ変換テーブル(look-up table)を作成しておき、それを参照すると演算時間を短縮することができる。ここで、ルックアップテーブルとは、コンピュータで大きな負担がかかる処理をコンピュータに行わせる場合、予め先に計算できるデータは計算しておき、その値を配列として保存しておき、この配列から目的のデータを取り出すことによって、計算の負担を軽減し効率アップを図るが、このような配列を言う。
また、整数演算に限定することで、たとえば、１４ビットの整数の二乗であれば、事前に１６３８４回の計算を行い、変換テーブルを作成しておけば、アドレス変換のみで二乗演算を行うことができ、大幅な処理の高速化が実現される。
（ウィナーフィルターの利用による点応答関数形状の改善）
以上の説明のようにして取得したスペクトルエンベロープ信号に基づき、スペクトルエンベロープ信号がガウス分布となるような逆フィルターを作成し、ＯＣＴで取得したＯＣＴ干渉信号をこの逆フィルターにかけて補正（較正）することで、光源のスペクトルに由来するＯＣＴ画像の深さ方向の分解能の低下を改善する。
具体的には、図４又は図５において、逆フィルター作成手段１８において、逆フィルターとして機能するウィナーフィルター（Wiener filter）を作成し、このウィナーフィルターを利用することで、点応答関数（point spread function: ＰＳＦ）形状を改善することができる。
ＰＳＦの形状が広がったり、サイドローブが生じたり、複数のピークを持つと、点状の物体の像の形状が広がったり、サイドローブが生じたり、複数の点として観測されるため、ＯＣＴ画像の画質が低下する。従って光学系などで発生するＰＳＦの形状の劣化を補正することにより画質が向上する。
ウィナーフィルターは、通常の逆フィルターに信号と雑音のパワーを考慮し雑音を含んだ画像の回復に用いられるものである。信号強度Ｓ（ｋ）はスペクトル干渉信号Ｆ（ｋ）と雑音（ノイズ）Ｎ（ｋ）を用いてＳ（ｋ）＝Ｆ（ｋ）Ｅ（ｋ）＋Ｎ（ｋ）となる。ここで、Ｅ（ｋ）はスペクトルエンベロープである。単純な逆フィルターは、Ｗ（ｋ）＝１／Ｅ（ｋ）となるが、ウィナーフィルターは、ノイズのパワースペクトルＰＮ（ｋ）とノイズのパワースペクトルＰＳ（ｋ）を用いてＷ（ｋ）＝１／（Ｅ（ｋ）＋ＰＮ（ｋ）＋ＰＳ（ｋ））
であらわされ、劣化信号の復元に広く用いられている。
以下、上記同様に、スペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）、フーリエドメインＯＣＴ（ＦＤ−ＯＣＴ）の場合、及び波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）の場合について、それぞれ説明する。
観測される生のＯＣＴ干渉信号Ｓ（ｋ）（図１０（ａ）参照）は、光源や光学系のスペクトル特性が理想的であった時のＯＣＴ干渉信号Ｆ（ｋ）（図１０（ｂ）参照）と光源や光学系のスペクトル特性Ｅ（ｋ）（図１０（ｃ）参照）の積で表される。つまり、Ｓ（ｋ）＝Ｆ（ｋ）×Ｅ（ｋ）となる。
光源や光学系のスペクトル特性が理想的であった時のＯＣＴ干渉信号が最も分解能が高く、ＰＳＦの形も理想的で、雑音が最小の画像が得られる。観測される生のＯＣＴ干渉信号から光源や光学系のスペクトル特性が理想的であった時のＯＣＴ干渉信号を求める発明が以下のものである。
（１）スペクトルドメインＯＣＴ（ＳＤ−ＯＣＴ）の場合
（ｉ）図３（ａ）〜（ｃ）により説明した上記の場合と同様に、スペクトルエンベロープ信号を抽出する。
（ii ）抽出されたエンベロープ信号に基づいて、ウィナーフィルター（Wiener filter）などの逆フィルターを作成する。
（iii）この逆フィルターを、スペクトルエンベロープ信号を抽出するために使用したもとのＯＣＴ干渉信号（平滑化される前のＯＣＴ信号）に乗ずることにより、光源のスペクトル分布を補正することができる。これにより、測定光学系の点応答関数（point spread function：ＰＳＦ）の形状を改善することができる。このＰＳＦの形を改善することにより、ＯＣＴ画像のクロストークやボケを改善することができる。
（２）波長走引型ＯＣＴ（ＳＳ−ＯＣＴ）の場合
（ｉ）上記イ．またはロ．の方法で、スペクトルエンベロープ信号を抽出する。
（ii ）抽出されたエンベロープ信号において、スペクトル変換フィルターを設計する。
スペクトル変換フィルターは、上記ウィナーフィルター（Wiener filter）などの逆フィルターを用いる。また、スペクトル変換フィルターは、Ｗ（ｋ）＝＜Ｓ（ｋ）＞／（＜Ｓ（ｋ）２＞＋ｎｃ）×Ｇａｕｓｓ（ｋ）を用いることもできる。ここで、ｎｃはシステム感度に依存する定数でスペクトルの形状をより平坦化するように決める。Ｇａｕｓｓ（ｋ）はガウシアン窓関数である。
スペクトルフィルターはｋを変数とする関数（Ｗ（ｋ））（図１０（ｅ）参照）になる。理想的な（スペクトルエンベロープが補正された）干渉信号Ｆ（ｋ）（図１０（ｆ）参照）は、計測した生のＯＣＴ信号Ｓ（ｋ）（図１０（ｄ）参照）と、ウィナーフィルター（Ｗｉｅｎｅｒ ｆｉｌｔｅｒ）などの逆フィルター、又は上記Ｗ（ｋ）＝＜Ｓ（ｋ）＞／（＜Ｓ（ｋ）２＞＋ｎｃ）×Ｇａｕｓｓ（ｋ）で示されるスペクトル変換フィルターの積で計算される。つまり、図１０（ｄ）〜（ｆ）に示すように、Ｆ（ｋ）＝Ｓ（ｋ）×Ｗ（ｋ）で計算される。
このスペクトル変換フィルターを、スペクトルエンベロープ信号を抽出するために使用したもとのＯＣＴ干渉信号に乗ずることにより、光源のスペクトル分布と光学系のスペクトル透過特性を補正することができる。
ＰＳＦは点物体の結像特性を示すもので、ＰＳＦの広がりが小さいほど、よい光学系とされる。ＯＣＴシステムにおけるＰＳＦは光学系の特性と光源のスペクトルエンベロープの形状のガウシアンからのずれにより劣化する。本発明により光源のスペクトル分布と光学系のスペクトル透過特性を補正することにより、測定光学系のＰＳＦの形状を改善することができる。
そして、ＰＳＦの形を改善することにより、ＯＣＴ画像のクロストークやボケを改善することができる。ＰＳＦの形状が広がったり、サイドローブが生じたり、複数のピークを持つと、点状の物体の像の形状が広がったり、サイドローブが生じたり、複数の点として観測されるため、ＯＣＴ画像の画質が低下する。従って光学系などで発生するＰＳＦの形状の劣化を補正することにより画質が向上する。
（情報エントロピーを利用した分散補正）
従来、ＯＣＴ画像の位相情報の最適化による後処理による画像の改善の方法はいくつか提案されている。その一例を挙げると、次のとおりである。
例えば、マサチューセッツ総合病院（ＭＧＨ）の方法（ｈｔｔｐ：／／ｗｗｗ．ｏｐｔｉｃｓｉｎｆｏｂａｓｅ．ｏｒｇ／ａｂｓｔｒａｃｔ．ｃｆｍ？ｉｄ＝８０１５０）は、眼の網膜における黄斑部のくぼみの形状を最適化するように画像処理を行う。この方法ではこのくぼみがＯＣＴ像に写っていないと適用できない。
また、マサチューセッツ工科大学（ＭＩＴ）の方法（ｈｔｔｐ：／／ｗｗｗ．ｏｐｔｉｃｓｉｎｆｏｂａｓｅ．ｏｒｇ／ａｂｓｔｒａｃｔ．ｃｆｍ？ｉｄ＝８０１４７）では、画像鮮明度（ｓｈａｐｎｅｓｓ）を評価関数とし、それを最大化する。この評価関数は多くのローカルミニマム（局所的な最適解）を持つため最適化が困難である。
そこで、本発明者は、画像の情報エントロピーを評価関数とし、その値を最小化するように画像の位相情報の最適化について発明した。これについて説明する。
（ａ）まず本実施例の上記「スペクトルエンベロープ信号の取得」で説明した場合と同様に光源および光学系のスペクトル特性を補正されたＯＣＴ画像を得る。この画像には、光源や光学系の位相特性が線形でないことに起因する高次の位相分散が含まれている。位相分散はＯＣＴ画像に対して、画質の劣化や分解能の低下などをもたらす。
本発明の画像の位相情報の最適化は、この位相分散を自動的に補償することである。ＯＣＴ画像の位相の２次の分散を補償することによって、ＯＣＴ画像の画質を向上させることができ、分解能も向上する。
（ｂ）ＯＣＴ画像の位相の２次の分散量を変化させると、ＯＣＴ画像の画質を変えることができる。ＯＣＴ画像の情報エントロピー量を、画像の２次の位相量の繰り返し最適化のための評価関数として用いる。情報エントロピーは、数式２で表される。ここで、数式３は、Ａ−，Ｂ−スキャン画像Ｉ（ｘ、ｚ）を規格化したものである。
Figure 0004389032
Figure 0004389032
この情報エントロピーの値は、画像の２次の位相分散量を変化させたときなめらかに、かつ単調に変化するため、上記マサチューセッツ工科大学（ＭＩＴ）の方法のように、他の評価関数などを用いる場合と比較して、ローカルミニマムに陥ることなく、位相の補正係数についてグローバルな最適解を求めることができる。
本発明では、最適値は繰り返し最適化法で求める。２次の位相分散量を微小に変化させ、評価関数を計算し、評価関数の値が小さくなれば、前記２次の位相分散量の変化を容認し、大きくなる場合は、前記２次の位相分散変化量を取り消し逆方向に変化させ、再び評価関数を計算し、その値について同様の手順で、位相分散変化量を繰り返し変化させていき、その値が最小値となる２次の位相分散量を求める（図１１（ａ）参照）。
情報エントロピー関数は、ローカルミニマムを持たないため、２次の位相分散量である変化量の初期値がＡ、Ｂ、Ｃ、Ｄいずれの場合でも、評価関数の値を小さくするように変化量を変えていった場合最適値に到達することができる。
上記他の評価関数を用いる場合に到達するローカルミニマムとは、真の最小値ではなく、周囲より値の小さい局所的な最小値のことであり（図１１（ｂ）参照）、この値を用いても最適な画質改善は得られない。つまり、変化量の初期値がＣあるいはＤであれば、評価関数の値を小さくするように変化量を変えていった場合最適値に到達するが、変化量の初期値がＡあるいはＢの場合、評価関数の値を小さくするように変化量を変えていった場合、最終的にローカルミニマムに落ちいって安定してしまう。
これに対して、情報エントロピー関数は、ローカルミニマムを持たないため、繰り返し回数も少なくてすみ、高速にＯＣＴ画像の最適化を行うことができる。
なお、３次元の体積画像の最適化を行う場合は、上記「位相の補正係数」は、最初のフレーム（Ａ、Ｂスキャン２次元画像）について計算されたものを、すべてのフレーム（Ｃスキャン方向の画像）に適用すればよい。このため、計算量が少なく、高速な画像補整が可能である。
（ＯＣＴの２次元画像のプレビュー）
ところで、ＯＣＴの２次元画像をプレビューする場合は、画像は、ある程度大まかでもよいが、高速にプレビュー画像が得られることが望ましい。
従来、ＯＣＴの２次元画像のプレビュー方法として、ＯＣＴ信号をＡスキャン方向（深さ方向）に積分して、干渉信号を平均化し、強度反射・散乱画像を得るという方法が知られている（ｈｔｔｐ：／／ｗｗｗ．ｏｐｔｉｃｓｉｎｆｏｂａｓｅ．ｏｒｇ／ａｂｓｔｒａｃｔ．ｃｆｍ？ｉｄ＝８２３８１参照）
しかしながら、従来技術では、大まかな２次元画像のプレビューしかできない。そこで、本発明で取得したスペクトルエンベロープ信号を利用すれば、光源のスペクトルに影響されず、ノイズの少ないＯＣＴの２次元画像を高速でプレビューすることが可能となる。
ＯＣＴ信号は、反射・散乱光と参照光の干渉信号に、光源や光学系のスペクトルエンベロープが加算されたものとなっている。従って被測定試料の情報は干渉信号にのみ含まれていると考えられる。この２次元画像のプレビューでは、ＯＣＴで得たＯＣＴ干渉信号（スペクトル干渉縞）から、本発明で取得したスペクトルエンベロープ信号を差し引き、干渉成分のみをぬきだし、干渉成分のみから成るＯＣＴ干渉信号を得る。
ＯＣＴ干渉信号は正負の値を持っているので、干渉信号から反射・散乱強度情報を取り出すためには正の符号のみにする必要がある。また、プレビュー画面では半透明被測定試料の反射・散乱体の情報を奥行き方向（Ａスキャン方向）を圧縮して、ＢおよびＣスキャン２次元画像が確認できればよいので、奥行き方向に積分してり和をとって情報を圧縮する。ＯＣＴ干渉信号を二乗し奥行き方向に対して積分する（和をとる）ことで、光源のスペクトルに影響されずノイズの少ない２次元のプレビュー画像を高速に得ることができる。
なお、２次元画像のプレビューでは、ＳＤ−ＯＣＴの場合は、ＯＣＴ干渉信号（スペクトル干渉縞）からスペクトルエンベロープ信号を差し引けばよいが、ＳＳ−ＯＣＴの場合は、前記したとおりのバランス検出されたＯＣＴ干渉信号から平均化されたＯＣＴ干渉信号を差し引いてもよい。
また、上記処理のうち、二乗処理が一番時間のかかる処理である。そこで、計算時間を短縮することを目的として、次のような処理をしてもよい。
（イ）上記のような二乗処理をすることなく、ＯＣＴ干渉信号の絶対値を積分してもよい。
（ロ）ルックアップテーブル（look-up table）を利用することで上記処理の処理時間を短縮することができる。
さらに、この２次元画像のプレビューにおいて上記信号処理の高速化を図る手段として、次のような手段がある。
（ハ）上記信号処理では、浮動小数点・固定小数点を使用することで行うことができるが、するために、上記信号処理を整数型で行う（浮動小数点演算を使わない）ことで、高速化を図る。
（ニ）上記スペクトルエンベロープ信号を差し引き干渉成分のみから成るＯＣＴ干渉信号の一部分のみを用いることで、高速に２次元画像を表示することができる。この一部分としては、例えば、以下の部分を抽出して表示する。
・ＯＣＴ干渉信号のモジュレーションの大きい部分。モジュレーションとは干渉縞の山の光強度をＩｍａｘ、谷の光強度Ｉｍｉｎをとすると、（Ｉｍａｘ−Ｉｍｉｎ）／（Ｉｍａｘ＋Ｉｍｉｎ）で定義される値である。
・光源のスペクトルの中心付近に対応した画像
・散乱が大きいなど、ＯＣＴ干渉信号の強度が大きい部分
・画像のＢあるいはＣスキャンラインを規則的に間引く
（デスペックルフィルター）
ＯＣＴ画像には、光源のコヒーレンスに起因するスペックルパターンが重畳することが多い。これを低減するために、注目している画素の付近の、異なるＢスキャン画像の、対応する画素を複数個選択し、その中で最大値（もっとも明るいデータ）で注目している画素における光強度の値を置き換える。
このような画素の置き換え処理を、本発明では、前記情報エントロピーを利用した分散補正をした後で行うと、光源のコヒーレンスに起因するスペックルパターンが取り除かれたよりノイズのないＯＣＴ画像を得ることができる。
（本発明の応用例）
本発明の画像処理装置は、各種の工業分野、動植物の研究・観察分野、各種文化財の解析・鑑定技術の分野、ロボット技術分野、医療用検査装置等の各種の技術分野に応用可能である。近年、ロボット技術が急速に発展し、植物、昆虫、動物、人等の器官（一例として、昆虫、動物、人の眼、皮膚、鼻等）の構造及び機能の研究を目的として、人間型ロボットに適用する技術開発も行われてる。例えば、人の眼の構造物の観察に本発明の画像処理装置を適用する例として、眼底血管を造影する場合について本発明の画像処理装置を適用する場合について説明する。
眼底には、網膜、脈絡膜などがあり、それぞれ別個の血管網がある。従来の眼底写真では、それらが重なって撮影されてしまう。そのために、それぞれの血管について蛍光造影剤を注入し造影する方法が行われているが、これは人体への負担が大きい。
このような従来の眼底血管を造影する場合について、本発明はきわめて有用である。ここでは、散乱光ＯＣＴを利用したＯＣＡ（Optical coherence angiography：光コヒーレンス血管造影法）において、本発明の画像処理装置を適用し、眼底の３次元ＯＣＴ画像から各層内の血管網（特定の構造物）の抽出及びそれに必要な装置を説明する。
まず、眼底血管の位置する部位を被計測物体としたＯＣＴ計測を行いＯＣＴ信号を得て、このＯＣＴ信号を利用して、上記本発明の画像処理装置により、スペクトルエンベロープを補正したＯＣＴ画像を得る。
次に、このＯＣＴ画像について、その層構造を自動的に分離する。この分離は、本願出願人の出願した先願発明である特願２００６−４４５６７（特開２００７−２２５３４９号公報参照）に係る発明を利用する。この先願発明の詳細の説明はここでは省略するが、その概要は、光コヒーレンストモグラフィーを用いて、被計測物体の奥行き方向の軸（Ａスキャン軸）に平行な２次元断層画像（Ｂスキャン画像）を、該２次元断層画像に垂直方向（Ｂスキャン方向）に位置をずらしながら複数取得し３次元画像を構成する光コヒーレンストモグラフィーによる光断層画像化を得る技術手段であって、被計測物体の表面の抽出および平面化の技術手段であり、また被計測物体の層構造を分離する技術手段である。
表面の抽出および平面化の技術手段では、上記取得した３次元画像から、被計測物体表面を、深さ方向の画像強度の最初のピーク位置から自動抽出し、抽出された表面形状を、多項式関数でフィッティングを行い、それを平面化することにより、内部の構造を非可燃体表面からの深さ情報３次元画像に自動変換する構成を特徴とする。
また、層構造を分離する技術手段では、その取得画像から、被計測物体内部の表面下の層構造の厚みの分布を、表面による信号のピークから深さ方向に向かって最初のピークまでの極小値の位置から自動抽出する。これにより、表面に近い層とその下の層に分離することができる構成を特徴とする。
上記技術手段、あるいは他の技術手段を用いて、分離した複数の層のそれぞれの層内で層表面からの深さｚ’の面での画素の光強度の平均値μｚ’と標準偏差σｚ’を求める。血管の散乱光は小さく、ＯＣＴ信号は小さいので、Ｉｎ＜Ｉｚ’（ｘ、ｙ）＜μｚ’−１／２σｚ’の光強度を持つ画素Ｉｚ’（ｘ、ｙ）、即ち、被計測物体のＯＣＴ画像の、位相情報を無視し、散乱光強度に対応した画素を抽出する。
ここで、Ｉｎはノイズレベルで決まる定数である。ｚ’はＡスキャン方向の特定の深さ、ｘ、ｙはＢおよびＣスキャン方向である。ｚ’をＡスキャン方向に動かして、３次元画像内のすべての体積中を上記の手順を繰り返し、抽出された画素を求めると、３次元体積内の血管網を構築することができる。
抽出された３次元分布について、３次元のマスクを作成し、そのマスクがかかっている部分のＯＣＴ画像の部分を元のＯＣＴ３次元画像から抽出して、ＯＣＴ画像の位相を最適化することによって、ＯＣＴ画像の鮮明度及び分解能を向上させ、血管の３次元的な分布が高いコントラストで得られる。
以上、本発明を実施するための最良の形態を実施例に基づいて説明したが、本発明はこのような実施例に限定されることなく、特許請求の範囲記載の技術的事項の範囲内で、いろいろな実施例があることは言うまでもない。
本発明に係る光コヒーレンストモグラフィーの画像処理装置は、以上のような構成であるから、被計測物体の奥行方向に高精度な分解能が要求される各種の技術分野、例えば、半導体製品の製造等の工業分野、動物の生体観察や植物の構造観察のような動植物の研究・観察分野、各種文化財の解析・鑑定技術の分野、ロボット技術分野（植物、昆虫、動物、人等の各部器官を観察し、その構造・機能をロボット技術に応用するための技術）、医療用検査装置等の各種の技術分野で利用される被計測物体の光断層画像取得用ＯＣＴに有用である。
本発明の前提となるＦＤ−ＯＣＴの全体構成を説明する図である。 本発明の前提となるＳＳ−ＯＣＴの全体構成を説明する図である。 本発明の原理を説明する図である。 本発明の実施例を説明する図である。 本発明の実施例を説明する図である。 従来の基本的なＯＣＴを説明する図である。 光源スペクトルの点応答関数を説明する図である。 本発明におけるヒルベルト変換（Hilbert変換）を用いる場合を説明する図である。 本発明におけるフーリエ変換を用いる場合を説明する図である。 本発明におけるスペクトル変換フィルターを説明する図である。 本発明における画像の情報エントロピーを評価関数として用い、画像の位相情報を最適化する場合を説明する図である。
符号の説明
１ ＦＤ−ＯＣＴ
２ 広帯域光源
３ 低コヒーレンス干渉計
４ 分光器
５ ビームスプリッター
６、９、３２、３４、３５、３６ レンズ
７ ガルバノミラー
８ 被計測物体
１０ 参照鏡
１１ 回折格子
１２ ＣＣＤ
１３、２１ 画像処理装置
１４ ＣＣＤ光検出器
１５ ＯＣＴ干渉信号入力部
１６ ＯＣＴ干渉信号平滑化手段
１７ エンベロープ検出手段
１８ 逆フィルター作成手段
１９ エンベロープ補正手段
２０ ＯＣＴ干渉信号出力部
２２ 干渉信号変換エンベロープ検出手段
２４ 波長走査型ＯＣＴ
２５ 波長走査型光源
２６、２８、３７ ファイバー
２７ ファイバーカップラー
２９、４８ 被計測物体
３０ ファイバー
３１ 固定参照鏡
３３ 走査鏡
３８ 光検知器
３９ コンピュータ
４０ ディスプレー
４３ ＯＣＴ
４４ 光源
４５ コリメートレンズ
４６ ビームスプリッター
４７ 物体アーム内の対物レンズ
４９ 参照アーム内の対物レンズ
５０ 参照鏡
５１ 集光レンズ
５２ （フォトダイオード等）光検出器

Claims (5)

  1. 光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、
    前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、
    該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、
    該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、
    該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、
    光源のスペクトルに由来するＯＣＴ画像の深さ方向の分解能の低下を改善することを特徴とする光コヒーレンストモグラフィーの画像処理装置。
  2. 光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、
    前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、
    該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、
    該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、
    該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、
    該エンベロープ補正手段で補正して得られたＯＣＴ画像の情報エントロピーを評価関数とし、画像の位相を操作しながら評価関数の値をもとめ、その値を最小化するように、位相を操作するというステップを繰り返し、情報エントロピーの値が最小となった画素それぞれの位相を最適位相とし、ＯＣＴ画像の位相を最適化することによって、ＯＣＴ画像の鮮明度及び分解能を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置。
  3. 光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、
    前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化するＯＣＴ干渉信号平滑化手段と、
    該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、を備えており、
    ＯＣＴ干渉信号（スペクトル干渉信号）からスペクトルエンベロープを差し引き、干渉成分のみをぬきだし、干渉成分のみから成るＯＣＴ干渉信号を得て、２次元画像のプレビューの高速化をすることを特徴とする光コヒーレンストモグラフィーの画像処理装置。
  4. 光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、
    前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化し、又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、
    該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、
    該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、
    該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、を備えており、
    該エンベロープ補正手段で補正して得られたＯＣＴ画像の、光源のコヒーレンスに起因するスペックルを除去するために、注目している画素における光強度の値を、その画素の周囲の画素における光強度の値のうち、もっとも大きな値で置き換えることにより、画像のスペックルを除去し、ＯＣＴ画像の出力を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置。
  5. 光コヒーレンストモグラフィーによって、被計測物体の奥行き方向の軸（Ａスキャン軸）に垂直方向（Ｂスキャン方向）に位置をずらしながら複数の１次元断層画像を取得して被計測物体の奥行き方向の軸に平行な２次元断層画像（Ｂスキャン画像）を取得し、さらに、該２次元断層画像に垂直方向（Ｃスキャン方向）に位置をずらしながら複数２次元断層画像を取得し３次元画像を形成する光コヒーレンストモグラフィーの画像を処理する光コヒーレンストモグラフィー画像処理装置において、
    前記Ａスキャン方向のＯＣＴ干渉信号（スペクトル干渉信号）を、Ｂスキャン方向に合計してその平均をとって平滑化する又はヒルベルト変換やフーリエ変換を用いて平滑化するＯＣＴ干渉信号平滑化手段と、
    該平滑化したＯＣＴ干渉信号についてスペクトルエンベロープ信号を求めるエンベロープ検出手段と、
    該スペクトルエンベロープ信号に基づき、該スペクトルエンベロープ信号がガウス分布となるように補正するための逆フィルターを求める逆フィルター作成手段と、
    該逆フィルターにより、ＯＣＴ干渉信号（スペクトル干渉信号）を補正するエンベロープ補正手段と、
    該エンベロープ補正手段で補正して得られたＯＣＴ画像を高コントラスト化するために、そのＯＣＴ画像の、位相情報を無視し、強度情報のみを用い、被計測物体の散乱光強度に対応した画素を抽出し、３次元マスクを作成する手段と、を備えており、
    そのマスクがかかっている部分のＯＣＴ画像を抽出することにより、ＯＣＴ画像から特定の構造物を抽出して、ＯＣＴ画像の位相を最適化することによって、ＯＣＴ画像の鮮明度及び分解能を向上させることを特徴とする光コヒーレンストモグラフィーの画像処理装置。
JP2007009574A 2007-01-18 2007-01-18 光コヒーレンストモグラフィーの画像処理装置 Active JP4389032B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007009574A JP4389032B2 (ja) 2007-01-18 2007-01-18 光コヒーレンストモグラフィーの画像処理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007009574A JP4389032B2 (ja) 2007-01-18 2007-01-18 光コヒーレンストモグラフィーの画像処理装置

Publications (2)

Publication Number Publication Date
JP2008175698A JP2008175698A (ja) 2008-07-31
JP4389032B2 true JP4389032B2 (ja) 2009-12-24

Family

ID=39702826

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007009574A Active JP4389032B2 (ja) 2007-01-18 2007-01-18 光コヒーレンストモグラフィーの画像処理装置

Country Status (1)

Country Link
JP (1) JP4389032B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9625250B2 (en) 2012-04-27 2017-04-18 Samsung Electronics Co., Ltd. Apparatus and method of generating tomographic image
US11231272B2 (en) 2017-11-16 2022-01-25 Nec Corporation Optical measuring apparatus and optical measuring method

Families Citing this family (69)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8062316B2 (en) 2008-04-23 2011-11-22 Avinger, Inc. Catheter system and method for boring through blocked vascular passages
US9125562B2 (en) 2009-07-01 2015-09-08 Avinger, Inc. Catheter-based off-axis optical coherence tomography imaging system
US8696695B2 (en) 2009-04-28 2014-04-15 Avinger, Inc. Guidewire positioning catheter
JP5355994B2 (ja) * 2008-11-05 2013-11-27 株式会社ニデック 眼科撮影装置
JP5242473B2 (ja) * 2009-03-23 2013-07-24 株式会社ニデック 眼科撮影装置、及び眼科撮影装置のキャリブレーション方法
JP5242492B2 (ja) * 2009-04-28 2013-07-24 株式会社トーメーコーポレーション ３次元画像処理装置
JP6101078B2 (ja) 2009-05-28 2017-03-29 アビンガー・インコーポレイテッドＡｖｉｎｇｅｒ， Ｉｎｃ． バイオイメージングのための光コヒーレンストモグラフィ
JP5473429B2 (ja) 2009-06-25 2014-04-16 キヤノン株式会社 眼底撮像装置及びその制御方法
JP5483408B2 (ja) * 2009-06-26 2014-05-07 四国電力株式会社 連続式濃度測定装置および方法
EP2448502B1 (en) 2009-07-01 2022-04-06 Avinger, Inc. Atherectomy catheter with laterally-displaceable tip
JP5879262B2 (ja) * 2009-07-01 2016-03-08 アビンガー・インコーポレイテッドＡｖｉｎｇｅｒ， Ｉｎｃ． カテーテルベースの軸外光コヒーレンストモグラフィイメージングシステム
GB2478590A (en) * 2010-03-12 2011-09-14 Precitec Optronik Gmbh Apparatus and method for monitoring a thickness of a silicon wafer
JP5466061B2 (ja) * 2010-03-26 2014-04-09 株式会社トプコン 眼科画像処理装置
US11382653B2 (en) 2010-07-01 2022-07-12 Avinger, Inc. Atherectomy catheter
US10548478B2 (en) 2010-07-01 2020-02-04 Avinger, Inc. Balloon atherectomy catheters with imaging
WO2014039096A1 (en) 2012-09-06 2014-03-13 Avinger, Inc. Re-entry stylet for catheter
WO2012003430A2 (en) 2010-07-01 2012-01-05 Avinger, Inc. Atherectomy catheters with longitudinally displaceable drive shafts
WO2012145133A2 (en) 2011-03-28 2012-10-26 Avinger, Inc. Occlusion-crossing devices, imaging, and atherectomy devices
US9949754B2 (en) 2011-03-28 2018-04-24 Avinger, Inc. Occlusion-crossing devices
JP5827024B2 (ja) 2011-03-31 2015-12-02 株式会社吉田製作所 光干渉断層画像生成装置の制御装置、制御方法及び制御プログラム
JP5787255B2 (ja) * 2011-07-12 2015-09-30 国立大学法人 筑波大学 Ｐｓ−ｏｃｔの計測データを補正するプログラム及び該プログラムを搭載したｐｓ−ｏｃｔシステム
EP3653151A1 (en) 2011-10-17 2020-05-20 Avinger, Inc. Atherectomy catheters and non-contact actuation mechanism for catheters
US9345406B2 (en) 2011-11-11 2016-05-24 Avinger, Inc. Occlusion-crossing devices, atherectomy devices, and imaging
JP5900950B2 (ja) * 2012-01-05 2016-04-06 国立大学法人 筑波大学 波長走査型光干渉断層計及びその位相安定化プログラム
JP5664564B2 (ja) * 2012-01-24 2015-02-04 住友電気工業株式会社 光断層画像取得方法
US11406412B2 (en) 2012-05-14 2022-08-09 Avinger, Inc. Atherectomy catheters with imaging
EP2849660B1 (en) 2012-05-14 2021-08-25 Avinger, Inc. Atherectomy catheter drive assemblies
WO2013172972A1 (en) 2012-05-14 2013-11-21 Avinger, Inc. Optical coherence tomography with graded index fiber for biological imaging
US9498247B2 (en) 2014-02-06 2016-11-22 Avinger, Inc. Atherectomy catheters and occlusion crossing devices
US11284916B2 (en) 2012-09-06 2022-03-29 Avinger, Inc. Atherectomy catheters and occlusion crossing devices
EP2967507B1 (en) 2013-03-15 2018-09-05 Avinger, Inc. Tissue collection device for catheter
US10932670B2 (en) 2013-03-15 2021-03-02 Avinger, Inc. Optical pressure sensor assembly
WO2014143064A1 (en) 2013-03-15 2014-09-18 Avinger, Inc. Chronic total occlusion crossing devices with imaging
US20140293289A1 (en) * 2013-03-27 2014-10-02 Kabushiki Kaisha Topcon Method for Generating Two-Dimensional Images From Three-Dimensional Optical Coherence Tomography Interferogram Data
JP6278295B2 (ja) * 2013-06-13 2018-02-14 国立大学法人 筑波大学 脈絡膜の血管網を選択的に可視化し解析する光干渉断層計装置及びその画像処理プログラム
JP6517198B2 (ja) 2013-07-08 2019-05-22 アビンガー・インコーポレイテッドＡｖｉｎｇｅｒ， Ｉｎｃ． 介入療法を案内する弾性層の識別
JP2015129730A (ja) * 2014-01-09 2015-07-16 住友電気工業株式会社 光学的測定方法
CA2938972A1 (en) 2014-02-06 2015-08-13 Avinger, Inc. Atherectomy catheters and occlusion crossing devices
JP6261450B2 (ja) * 2014-05-30 2018-01-17 株式会社トーメーコーポレーション 眼科装置
WO2016007652A1 (en) 2014-07-08 2016-01-14 Avinger, Inc. High speed chronic total occlusion crossing devices
JP6437364B2 (ja) * 2015-03-30 2018-12-12 株式会社Ｓｃｒｅｅｎホールディングス 信号処理方法および画像処理装置
JP6896699B2 (ja) 2015-07-13 2021-06-30 アビンガー・インコーポレイテッドＡｖｉｎｇｅｒ， Ｉｎｃ． 画像誘導治療／診断カテーテル用マイクロ成形アナモフィック反射器レンズ
JP6927986B2 (ja) 2016-01-25 2021-09-01 アビンガー・インコーポレイテッドＡｖｉｎｇｅｒ， Ｉｎｃ． 遅延補正を備えたｏｃｔイメージングカテーテル
JP6109369B2 (ja) * 2016-03-16 2017-04-05 キヤノン株式会社 撮像装置、制御装置及び制御方法
EP3435892B1 (en) 2016-04-01 2024-04-03 Avinger, Inc. Atherectomy catheter with serrated cutter
GB2549515A (en) * 2016-04-20 2017-10-25 Michelson Diagnostics Ltd Processing optical coherence tomography scans
US11344327B2 (en) 2016-06-03 2022-05-31 Avinger, Inc. Catheter device with detachable distal end
US11224459B2 (en) 2016-06-30 2022-01-18 Avinger, Inc. Atherectomy catheter with shapeable distal tip
JP6542178B2 (ja) * 2016-08-05 2019-07-10 株式会社トーメーコーポレーション 偏光情報を利用した光断層画像撮影装置
JP6456421B2 (ja) * 2017-03-09 2019-01-23 キヤノン株式会社 撮像装置
US10578722B2 (en) * 2017-03-27 2020-03-03 Honeywell International Inc. System and method for glare suppression and ranging
JP6754492B2 (ja) * 2018-01-18 2020-09-09 小出　珠貴 生体組織分析装置および生体組織分析方法
JP7019128B2 (ja) * 2018-01-22 2022-02-15 株式会社トーメーコーポレーション 光断層画像撮影装置
JP7113637B2 (ja) 2018-03-23 2022-08-05 株式会社Ｓｃｒｅｅｎホールディングス 光干渉断層撮像装置および光干渉断層撮像方法
US12167867B2 (en) 2018-04-19 2024-12-17 Avinger, Inc. Occlusion-crossing devices
JP7381827B2 (ja) 2019-05-30 2023-11-16 株式会社トプコン 光干渉測定装置および光干渉測定方法
JP6999908B2 (ja) * 2019-05-30 2022-01-19 株式会社トプコン 光干渉測定装置および光干渉測定方法
JP7146698B2 (ja) * 2019-06-26 2022-10-04 株式会社日立製作所 撮像装置、撮像方法、撮像システムの送信装置、及び撮像システムの受信装置
US11793400B2 (en) 2019-10-18 2023-10-24 Avinger, Inc. Occlusion-crossing devices
JP7629639B2 (ja) * 2019-11-14 2025-02-14 国立大学法人 筑波大学 画像処理装置、画像処理方法およびプログラム
KR102339811B1 (ko) * 2020-04-21 2021-12-16 주식회사 필로포스 누화, 에코신호 및 비선형 현상을 개선한 소형 oct용 분광기
JP7779065B2 (ja) * 2021-09-29 2025-12-03 株式会社ニデック 眼科撮影装置及び眼科撮影プログラム
CN113729615B (zh) * 2021-10-12 2023-01-13 中山大学中山眼科中心 具有手持式探头的光学相干断层成像设备
CN116067914B (zh) * 2022-12-23 2025-10-28 南开大学 一种材料表面及亚表面三维微观结构成像方法及应用
CN115984406B (zh) * 2023-03-20 2023-06-20 始终(无锡)医疗科技有限公司 一种深度学习和谱域空域联合子采样的ss-oct压缩成像方法
CN117109450B (zh) * 2023-07-25 2024-05-10 广东工业大学 一种自适应参数谱估计的高分辨光谱共焦干涉薄膜测厚法
CN118583967B (zh) * 2024-07-04 2025-09-19 广东贝洛新材料科技有限公司 连接器的检测方法、装置、计算机设备和存储介质
CN118521516B (zh) * 2024-07-19 2024-10-25 深圳明锐理想科技股份有限公司 一种深度图校正方法、电子设备及存储介质
CN119323538B (zh) * 2024-12-19 2025-02-28 灏麟(天津)生物科技有限公司 一种瘢痕诊断的光学相干断层成像数据优化方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9625250B2 (en) 2012-04-27 2017-04-18 Samsung Electronics Co., Ltd. Apparatus and method of generating tomographic image
US11231272B2 (en) 2017-11-16 2022-01-25 Nec Corporation Optical measuring apparatus and optical measuring method

Also Published As

Publication number Publication date
JP2008175698A (ja) 2008-07-31

Similar Documents

Publication Publication Date Title
JP2008175698A (ja) 光コヒーレンストモグラフィーの画像処理方法及び画像処理装置
US20240065552A1 (en) Intraoral oct with color texture
De Boer et al. Twenty-five years of optical coherence tomography: the paradigm shift in sensitivity and speed provided by Fourier domain OCT
JP5166889B2 (ja) 眼底血流量の定量測定装置
JP6195334B2 (ja) 撮像装置、撮像方法およびプログラム
US9243887B2 (en) Lateral distortion corrected optical coherence tomography system
US20070201033A1 (en) Methods and systems for performing angle-resolved fourier-domain optical coherence tomography
JP6765786B2 (ja) 撮像装置、撮像装置の作動方法、情報処理装置、及び情報処理装置の作動方法
JP2011087813A (ja) 光断層画像生成装置及び光断層画像生成方法
JP4461259B2 (ja) 光断層画像の処理方法
JP3688608B2 (ja) 分光機能を備えた光干渉断層画像計測装置
JP5557397B2 (ja) 半透明物質の画像化の方法および装置
JP7058901B1 (ja) ３次元撮像装置
WO2016023502A1 (en) Phase-inverted sidelobe-annihilated optical coherence tomography
JP5602363B2 (ja) 光干渉断層撮像装置
Huang et al. Applying a Pix2Pix generative adversarial network to a Fourier-domain optical coherence tomography system for artifact elimination
JP2007127425A (ja) 光断層画像化法における補正方法
Yang et al. High-speed low-cost line-field spectral-domain optical coherence tomography for industrial applications
JP6648891B2 (ja) 物質含有量を断層可視化する装置および方法
JP2019154764A (ja) 涙液層厚み測定装置及び方法
CN117999456A (zh) 用于样本成像的光学相干断层扫描系统和方法
JP2008151734A (ja) 光断層画像化方法、装置およびプログラムならびに光断層画像化システム
Dsouza Towards low cost multiple reference optical coherence tomography for in vivo and NDT applications
Liu et al. Progress toward inexpensive endoscopic high-resolution common-path OCT
JP5746741B2 (ja) 画像生成装置、画像生成システム及び画像生成方法

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090526

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090609

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090709

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150