JP2014110883A - 画像処理装置及び画像処理方法 - Google Patents

画像処理装置及び画像処理方法 Download PDF

Info

Publication number
JP2014110883A
JP2014110883A JP2013159180A JP2013159180A JP2014110883A JP 2014110883 A JP2014110883 A JP 2014110883A JP 2013159180 A JP2013159180 A JP 2013159180A JP 2013159180 A JP2013159180 A JP 2013159180A JP 2014110883 A JP2014110883 A JP 2014110883A
Authority
JP
Japan
Prior art keywords
image
display
image processing
discontinuous
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.)
Pending
Application number
JP2013159180A
Other languages
English (en)
Other versions
JP2014110883A5 (ja
Inventor
Yoshihiko Iwase
好彦 岩瀬
Makoto Sato
眞 佐藤
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP2013159180A priority Critical patent/JP2014110883A/ja
Priority to US14/065,759 priority patent/US20140121506A1/en
Priority to EP13190930.1A priority patent/EP2728549B1/en
Publication of JP2014110883A publication Critical patent/JP2014110883A/ja
Publication of JP2014110883A5 publication Critical patent/JP2014110883A5/ja
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/102Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02027Two or more interferometric channels or interferometers
    • G01B9/02028Two or more reference or object arms in one interferometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02029Combination with non-interferometric systems, i.e. for measuring the object
    • G01B9/0203With imaging systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time
    • G01B9/02067Active error reduction, i.e. varying with time by electronic control systems, i.e. using feedback acting on optics or light
    • G01B9/02068Auto-alignment of optical elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02083Interferometers characterised by particular signal processing and presentation
    • G01B9/02087Combining two or more images of the same region
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • G01B9/02091Tomographic interferometers, e.g. based on optical coherence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/70Using polarization in the interferometer
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10101Optical tomography; Optical coherence tomography [OCT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Pathology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Automation & Control Theory (AREA)
  • Optics & Photonics (AREA)
  • Signal Processing (AREA)
  • Eye Examination Apparatus (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

【課題】 偏光が解消された領域における萎縮箇所の状態によって、診断に有効な情報をユーザに提示すること。
【解決手段】 画像処理装置が、被検体の偏光断層画像から所定の領域を抽出する抽出手段と、前記抽出された所定の領域の不連続な箇所の状態に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させる表示制御手段と、を有する。
【選択図】 図8

Description

本発明は、被検体の画像を処理する画像処理装置及び画像処理方法に関する。
多波長光波干渉を利用した光コヒーレンストモグラフィ(Optical Coherence Tomography:以下、OCT)は、試料(特に眼底)の断層画像を高分解能に得ることができる。
近年、眼科用OCT装置において、眼底組織の形状をイメージングする通常のOCT画像に加えて、眼底組織の光学特性の一つである偏光パラメータ(リターデーションとオリエンテーション)を用いてイメージングする偏光OCT画像が取得されている。
偏光OCTは、偏光パラメータを利用して、偏光OCT画像を構成し、眼底組織の区別やセグメンテーションを行うことができる。偏光OCTは、試料を観察する測定光に円偏光に変調した光を用い、干渉光を2つの直交する直線偏光として分割して検出し、偏光OCT画像を生成する(特許文献1参照)。
WO2010/122118A1
ここで、一般的に、疾患がある被検眼の偏光OCT画像からRPE層の萎縮箇所(所定の層の不連続な箇所)を認識することができる。本発明の目的は、この萎縮箇所の状態によって、診断に有効な情報をユーザに提示することである。
本発明に係る画像処理装置は、
被検体の偏光断層画像を取得する断層画像取得手段と、
前記取得された偏光断層画像から所定の領域を抽出する抽出手段と、
前記抽出された所定の領域の不連続な箇所の状態に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させる表示制御手段と、を有する。
また、本発明に係る画像処理方法は、
被検体の偏光断層画像を取得する工程と、
前記取得された偏光断層画像から所定の領域を抽出する工程と、
前記抽出された所定の領域の不連続な箇所の状態に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させる工程と、を有する。
本発明によれば、被検体の偏光断層画像から抽出された偏光が解消された領域の不連続な箇所の状態(連続性等)に基づいて、不連続な箇所の種類を示す表示形態を表示手段に表示させることができる。これにより、偏光が解消された領域における萎縮箇所の状態によって、診断に有効な情報をユーザに提示することができる。
本実施形態における画像処理装置の全体構成の概略図。 信号処理部190で生成される画像の例。 本実施形態における処理フロー。 本実施形態における画像処理装置の表示部の表示画面における表示例。 本実施形態における断層画像と偏光解消領域とを説明するための図。 本実施形態における病変の断層画像と偏光解消領域とを説明するための図。 本実施形態における病変の眼底画像を説明するための図。 本実施形態における画像処理装置の表示部の表示画面における表示例。 本実施形態における病変の断層画像と偏光解消領域とを説明するための図。
本発明に係る撮影装置は、被検眼、皮膚、内臓等の被検体に適用することができる。また、本発明に係る撮影装置としては、例えば、眼科装置や内視鏡等である。以下、本発明の一例として、本実施形態に係る眼科装置について、図面を用いて詳細に説明する。
[装置の全体構成]
図1は、本実施形態における撮影装置の一例である「眼科装置」の全体構成の概略図である。なお、後述する信号処理部190の少なくとも一部を「画像処理装置」とみなすことができ、また、この場合、「眼科装置」全体を「眼科システム」、あるいは「撮影装置」全体を「撮影システム」とみなすこともできる。
本装置は、偏光OCT(Polarization Sensitive OCT;以下、PS−OCT)100、偏光を利用した走査型検眼鏡(Polarization Sensitive Scanning Laser Ophothalmoscope:以下、PS−SLO)140、前眼部撮像部160、内部固視灯170、制御部200から構成される。
内部固視灯170を点灯して被検眼に注視させた状態で、前眼部観察部160により観察される被検体の前眼部の画像を用いて、装置のアライメントが行われる。アライメント完了後に、PS−OCT100とPS−SLO140による眼底の撮像が行われる。
<PS−OCT100の構成>
PS−OCT100の構成について説明する。
光源101は、低コヒーレント光源であるSLD光源(Super Luminescent Diode)であり、例えば、中心波長850nm、バンド幅50nmの光を出射する。光源101としてSLDを用いたが、ASE光源(Amplified Spontaneous Emission)等、低コヒーレント光が出射できる光源であれば何れでも良い。
光源101から出射された光は、PM(Polarization Maintaining)ファイバ102、偏光コントローラ103を介して、偏光保持機能を有したファイバカップラ104に導かれ、測定光(以下、「断層画像用の測定光」や「OCT測定光」ともいう)と、測定光に対応する参照光とに分割される。
偏光コントローラ103は、光源101から出射された光の偏光の状態を調整するものであり、直線偏光に調整される。ファイバカップラ104の分岐比は、90(参照光):10(測定光)である。
測定光は、PMファイバ105を介してコリメータ106から平行光として出射される。出射された測定光は、眼底Erにおいて測定光を水平方向にスキャンするガルバノミラーから構成されるXスキャナ107、レンズ108、109、眼底Erにおいて測定光を垂直方向にスキャンするガルバノミラーから構成されるYスキャナ110を介し、ダイクロイックミラー111に到達する。Xスキャナ107、Yスキャナ110は、駆動制御部180により制御され、眼底Erの所望の範囲で測定光を走査することができる。なお、測定光が走査される眼底上の範囲は、断層画像(Bスキャンともいう)の取得範囲、断層画像の取得位置、測定光の照射位置としてみなすことができる。また、Xスキャナ107、Yスキャナ110は、PS−OCT用の走査手段の一例であり、共通のXYスキャナとして構成しても良い。ダイクロイックミラー111は、800nm〜900nmの光を反射し、それ以外の光を透過する特性を有する。
ダイクロイックミラー111により反射された測定光は、レンズ112を介し、光軸を回転軸としてP偏光からS偏光に対して45°傾けて設置されたλ/4偏光板113を通過することにより、位相が90°ずれ、円偏光の光に偏光制御される。なお、λ/4偏光板113は、測定光の偏光状態を調整する測定光用の偏光調整部材の一例である。ここで、後述するPS−SLO光学系を適用する場合、λ/4偏光板113をPS−OCT光学系の一部とPS−SLO光学系の一部との共通光路に設けることができる。これにより、PS−SLO光学系で取得した画像と、PS−OCT光学系で取得した画像とに生じる偏光状態のばらつきを比較的に抑制することができる。このとき、PS−SLO用の走査手段と、PS−OCT用の走査手段とは、互いに共役な位置に設けられ、被検眼の瞳と共役な位置に設けることができる。なお、λ/4偏光板113の傾きは、λ/4偏光板113の状態の一例であり、例えば、偏光ビームスプリッタを内蔵したファイバカップラ123の偏光分割面の光軸を回転軸とした所定の位置からの角度である。
また、λ/4偏光板113を光路に対して挿脱可能に構成することができる。例えば、光軸あるいは光軸に平行な軸を回転軸としてλ/4偏光板113を回転する機械的な構成が考えられる。これにより、SLO光学系とPS−SLO光学系とを簡単に切り換え可能な小型な装置を実現することができる。また、OCT光学系とPS−OCT光学系とを簡単に切り換え可能な小型な装置を実現することができる。
ここで、被検眼に入射される光は、λ/4偏光板を45°傾けて設置することで円偏光の光に偏光制御されるが、被検眼の特性により眼底Erにおいて円偏光とならない場合がある。そのため、駆動制御部180の制御により、λ/4偏光板の傾きを微調整できるように構成されている。
円偏光に偏光制御された測定光は、ステージ116上に乗ったフォーカスレンズ114により、被検体である眼の前眼部Eaを介し、眼底Erの網膜層にフォーカスされる。眼底Erを照射した測定光は各網膜層で反射・散乱し、上述の光学経路をファイバカップラ104に戻る。
一方、ファイバカプラ104で分岐された参照光は、PMファイバ117を介してコリメータ118から平行光として出射される。出射された参照光は測定光と同様に、光軸を回転軸としてP偏光からS偏光に対して22.5°傾けて設置されたλ/4偏光板119で偏光制御される。なお、λ/4偏光板119は、参照光の偏光状態を調整する参照光用の偏光調整部材の一例である。参照光は分散補償ガラス120介し、コヒーレンスゲートステージ121上のミラー122で反射され、ファイバカップラ104に戻る。参照光は、λ/4偏光板119を二度通過する事で直線偏光の光がファイバカップラ104に戻ることになる。
コヒーレンスゲートステージ121は、被検者の眼軸長の相違等に対応する為、駆動制御部180で制御される。なお、コヒーレンスゲートとは、測定光の光路における参照光の光路長に対応する位置のことである。本実施形態では、参照光の光路長を変更しているが、測定光の光路と参照光の光路との光路長差を変更できれば良い。
ファイバカップラ104に戻った戻り光と参照光とは合波されて干渉光(以下、「合波光」ともいう)となり、偏光ビームスプリッタを内蔵したファイバカップラ123に入射され、異なる偏光方向の光であるP偏光の光とS偏光の光とに分岐比50:50で分割される。
P偏光の光は、PMファイバ124、コリメータ130を介し、グレーティング131により分光されレンズ132、ラインカメラ133で受光される。同様に、S偏光の光は、PMファイバ125、コリメータ126を介し、グレーティング127により分光されレンズ128、ラインカメラ129で受光される。なお、グレーティング127、131、ラインカメラ129、133は、各偏光の方向に合わせて配置されている。
ラインカメラ129、133でそれぞれ受光した光は、光の強度に応じた電気信号として出力され、断層画像生成部の一例である信号処理部190で受ける。
λ/4偏光板113、119の傾きに関して、偏光ビームスプリッタの偏光分割面の傾きを基準に自動的に調整することができるが、眼底の視神経乳頭中心と黄斑中心を結んだ直線に対して自動的に調整しても良い。このとき、λ/4偏光板113、119の傾きを検知する傾き検知部(不図示)を有することが好ましい。この傾き検知部により、現在の傾きと所定の傾きになったことを検知することができる。もちろん、受光した光の強度に基づいて、λ/4偏光板113、119の傾き具合を検知し、所定の強度になるように傾きを調整しても良い。なお、後述するように、GUI上に傾きを示すオブジェクトを表示して、ユーザがマウスを用いて調整しても良い。また、偏光基準として鉛直方向を基準にして偏光ビームスプリッタ、λ/4偏光板113、119を調整しても同様の効果が得られる。
<PS−SLO140の構成>
PS−SLO140の構成について説明する。
光源141は、半導体レーザであり、本実施例では、例えば、中心波長780nmの光を出射する。光源141から出射された測定光(以下、「眼底画像用の測定光」や「SLO測定光」ともいう)は、PMファイバ142を介し、偏光コントローラ145で直線偏光になるよう偏光制御され、コリメータ143から平行光として出射される。出射された測定光は穴あきミラー144の穴あき部を通過し、レンズ155を介し、眼底Erにおいて測定光を水平方向にスキャンするガルバノミラーから構成されるXスキャナ146、レンズ147、148、眼底Erにおいて測定光を垂直方向にスキャンするガルバノミラーから構成されるYスキャナ149を介し、ダイクロイックミラー154に到達する。Xスキャナ146、Yスキャナ149は駆動制御部180により制御され、眼底上で所望の範囲を測定光で走査できる。なお、Xスキャナ146、Yスキャナ149は、PS−SLO用の走査手段の一例であり、共通のXYスキャナとして構成しても良い。ダイクロイックミラー154は、760nm〜800nmを反射し、それ以外の光を透過する特性を有する。
ダイクロイックミラー154にて反射された直線偏光の測定光は、PS−OCT100と同様の光路を経由し、眼底Erに到達する。
眼底Erを照射した測定光は、眼底Erで反射・散乱され、上述の光学経路をたどり穴あきミラー144に達する。穴あきミラー144で反射された光が、レンズ150を介し、偏光ビームスプリッタ151にて異なる偏光方向の光(本実施形態では、P偏光の光とS偏光の光)に分割され、アバランシェフォトダイオード(APD)152、153で受光され、電気信号に変換されて、眼底画像生成部の一例でもある信号処理部190で受ける。
ここで、穴あきミラー144の位置は、被検眼の瞳孔位置と共役となっており、眼底Erに照射された測定光が反射・散乱された光のうち、瞳孔周辺部を通った光が、穴あきミラー144によって反射される。
本実施例では、PS−OCT、PS−SLOともにPMファイバを用いたが、シングルモードファイバー(SMF)でも偏光コントローラを用い偏光を制御する事で同様の構成と効果が得られる。
<前眼部撮像部160>
前眼部撮像部160について説明する。
前眼部撮像部160は、波長1000nmの照明光を発するLED115−a、115−bから成る照明光源115により前眼部Eaを照射する。前眼部Eaで反射され光は、レンズ114、偏光板113、レンズ112、ダイクロイックミラー111、154を介し、ダイクロイックミラー161に達する。ダイクロイックミラー161は、980nm〜1100nmの光を反射し、それ以外の光を透過する特性を有する。ダイクロイックミラー161で反射された光は、レンズ162、163、164を介し、前眼部カメラ165で受光される。前眼部カメラ165で受光された光は、電気信号に変換され、信号処理部190で受ける。
<内部固視灯170>
内部固視灯170について説明する。
内部固視灯170は、内部固視灯用表示部171、レンズ172で構成される。内部固視灯用表示部171として複数の発光ダイオード(LD)がマトリックス状に配置されたものを用いる。発光ダイオードの点灯位置は、駆動制御部180の制御により撮像したい部位に合わせて変更される。内部固視灯用表示部171からの光は、レンズ172を介し、被検眼に導かれる。内部固視灯用表示部171から出射される光は520nmで、制御部180により所望のパターンが表示される。
<制御部200>
本装置全体を制御するための制御部200について説明する。
制御部200は、駆動制御部180、信号処理部190、表示制御部191、表示部192から構成される。
駆動制御部180は、上述の通り各部を制御する。
信号処理部190は、画像生成部193と画像解析部194から構成される。信号処理部190は、ラインカメラ129及び133、APD152及び153、前眼部カメラ165からそれぞれ出力される信号に基づき、画像の生成、生成された画像の解析、解析結果の可視化情報の生成を行う。なお、画像の生成、解析などの詳細については後述する。
表示制御部191は、眼底画像取得部(不図示)と断層画像取得部(不図示)により、断層画像生成部と眼底画像生成部とでそれぞれ生成された画像を取得した画像等を表示部192の表示画面に表示させる。ここで、表示部192は、例えば、液晶等のディスプレイである。なお、信号処理部190で生成された画像データは、表示制御部191に有線で送信されても良いし、無線で送信されても良い。この場合、表示制御部191を画像処理装置とみなすことができる。なお、撮影システムとして、眼底画像取得部がSLO光学系を含み、断層画像取得部がOCT光学系を含むように構成しても良い。なお、本明細書において、被検眼以外の被検体の場合、「眼底画像(眼底輝度画像)」を「平面画像(平面輝度画像)」と換言することができ、また、「眼底画像取得部」を「平面画像取得部」と換言することができる。
表示部192は、表示制御部191の制御の下、後述するように種々の情報を示す表示形態を表示する。なお、表示制御部191からの画像データは、表示部192に有線で送信されても良いし、無線で送信されても良い。また、表示部192等は、制御部200に含まれているが、本発明はこれに限らず、制御部200とは別に設けられても良い。また、表示制御部191と表示部192とを一体的に構成した、ユーザが持ち運び可能な装置の一例であるタブレットでも良い。この場合、表示部にタッチパネル機能を搭載させ、タッチパネル上で画像の表示位置の移動、拡大縮小、表示される画像の変更等の操作可能に構成することが好ましい。
[画像処理]
次に、信号処理部190を構成する画像生成部193における画像生成について説明する。
画像生成部193は、ラインカメラ129、133から出力されたそれぞれの干渉信号に対して、一般的なSD−OCT(Spectral Domain OCT)に用いられる再構成処理を行うことで、各偏光成分に基づいた2つの断層画像である第一の偏光に対応する断層画像と、第二の偏光に対応する断層画像とを生成する。
まず、画像生成部193は、干渉信号から固定パターンノイズ除去を行う。固定パターンノイズ除去は検出した複数のAスキャン信号を平均することで固定パターンノイズを抽出し、これを入力した干渉信号から減算することで行われる。
次に、画像生成部193は、干渉信号を波長から波数に変換し、フーリエ変換を行うことによって、偏光状態を示す断層信号を生成する。
以上の処理を2つの偏光成分の干渉信号に対して行うことにより、2つの断層画像が生成される。
また、画像生成部193は、APD152、153から出力された信号を、Xスキャナ146、Yスキャナ149の駆動に同期して整列させることにより、各偏光成分に基づいた2つの眼底画像である第一の偏光に対応する眼底画像と、第二の偏光に対応する眼底画像とを生成する。
<断層輝度画像あるいは眼底輝度画像の生成>
画像生成部193は、前述した2つの断層信号から断層輝度画像を生成する。
断層輝度画像は、従来のOCTにおける断層画像と基本的に同じもので、その画素値rは各ラインセンサ129、133から得られた断層信号AおよびAから(式1)によって計算される。
また、同様に、2つの眼底画像から眼底輝度画像を生成する。
図2(a)に視神経乳頭部の輝度画像の例を示す。
なお、表示制御部191は、λ/4偏光板113を光路から外している場合に、従来のOCTの手法により取得した断層輝度画像を表示部192に表示させても良いし、従来のSLOの手法により取得した眼底輝度画像を表示部192に表示させても良い。
<リターデーション画像の生成>
画像生成部193は、互いに直行する偏光成分の断層画像からリターデーション画像を生成する。
リターデーション画像の各画素の値δは、断層画像を構成する各画素の位置において、垂直偏光成分と水平偏光成分とが被検眼で受ける影響の比を示す値であり、各断層信号AおよびAから(式2)によって計算される。
図2(b)は、このように生成された視神経乳頭部のリターデーション画像の例を示したものであり、各Bスキャン画像に対して(式2)を計算することによって得ることができる。ここで、上述した通り、リターデーション画像は、2つの偏光が被検眼で受ける影響の違いを示す断層画像のことである。図2(b)は、上記比を示す値を断層画像としてカラーで表示しており、濃淡の濃い場所は上記比を示す値が小さく、濃淡の淡い場所は上記比を示す値が大きいことを表している。そのため、リターデーション画像を生成することにより、複屈折性のある層を把握することが可能となる。なお、詳細は、「E. Gotzinger et al., Opt. Express 13, 10217, 2005」に記載されている通りである。
また、同様に、信号処理部190は、APD152及び153からの出力に基づいて眼底の平面方向のリターデーション画像を生成することもできる。
<リターデーションマップの生成>
画像生成部193は、複数のBスキャン像に対して得たリターデーション(Retardation)画像からリターデーションマップを生成する。
まず、画像生成部193は、各Bスキャン画像において、網膜色素上皮(以下、「RPE」ともいう)を検出する。RPEは偏光を解消する性質を持っているため、各Aスキャンを深度方向に沿って内境界膜(以下、「ILM」ともいう)からRPEを含まない範囲でリターデーションの分布を調べ、その最大値を当該Aスキャンにおけるリターデーションの代表値とする。
画像生成部193は、以上の処理を全てのリターデーション画像に対して行うことにより、リターデーションマップを生成する。
図2(c)に視神経乳頭部のリターデーションマップの例を示す。濃淡の濃い場所は上記比を示す値が小さく、濃淡の淡い場所は上記比を示す値が大きいことを表している。視神経乳頭部において、複屈折性を持つ層としては網膜神経線維層(以下、「RNFL」ともいう)であり、リターデーションマップは、2つの偏光がRNFLの複屈折性とRNFLの厚みとで受ける影響の違いを示す画像である。そのため、RNFLが厚い箇所では上記比を示す値が大きくなり、RNFLが薄い箇所では上記比を示す値が小さくなる。したがって、リターデーションマップにより、眼底全体のRNFLの厚みを把握することが出来、緑内障の診断に用いることが出来る。
<複屈折マップの生成>
画像生成部193は、先に生成されたリターデーション画像の各Aスキャン画像において、ILMからRNFLの範囲でリターデーションδの値を線形近似し、その傾きを当該Aスキャン画像の網膜上の位置における複屈折として決定する。すなわち、リターデーションはRNFLにおける距離と複屈折と積であるため、各Aスキャン画像において深さとリターデーションの値をプロットすると線形の関係が得られる。したがって、このプロットに対して最小二乗法等により線形近似を行い、その傾きを求めればそれが当該Aスキャン画像におけるRNFLの複屈折の値となる。この処理を取得した全てのリターデーション画像に対して行うことで、複屈折を表すマップを生成する。
図2(d)に視神経乳頭部の複屈折マップの例を示す。複屈折マップは、複屈折の値を直接マップ化するため、RNFLの厚さが変化しない場合であっても、その繊維構造が変化した場合に、複屈折の変化として描出することができる。
<DOPU画像の生成>
画像生成部193は、取得した断層信号AH、とそれらの間の位相差ΔΦから、各画素毎にストークスベクトルSを(式3)により計算する。
ただし、ΔΦは2つの断層画像を計算する際に得られる各信号の位相ΦとΦからΔΦ=Φ−Φとして計算する。
次に画像生成部193は、各Bスキャン画像を概ね計測光の主走査方向に70μm、深度方向に18μm程度の大きさのウィンドウを設定し、各ウィンドウ内において数Cで画素毎に計算されたストークスベクトルの各要素を平均し、(式4)により当該ウィンドウ内の偏光の均一性DOPU(Degree Of Polarization Uniformity)を(式4)により計算する。
ただし、Q、U、Vは各ウィンドウ内のストークスベクトルの要素Q,U,Vを平均した値である。この処理をBスキャン画像内の全てのウィンドウに対して行うことで、図2(e)に示す視神経乳頭部のDOPU画像が生成される。ここで、上述した通り、DOPU画像は、2つの偏光の均一度を示す断層画像のことである。
DOPUは偏光の均一性を表す数値であり、偏光が保たれている箇所においては1に近い数値となり、偏光が解消された保たれない箇所においては1よりも小さい数値となるものである。網膜内の構造においては、RPEが偏光状態を解消する性質があるため、DOPU画像においてRPEに対応する部分は、他の領域に対してその値が小さくなる。図において、濃淡が淡い場所210がRPEを示しており、濃淡が濃い場所220が偏光が保たれている網膜層領域を示している。DOPU画像は、RPE等の偏光を解消する層を画像化しているので、病気などによりRPEが変形している場合においても、輝度の変化よりも確実にRPEを画像化出来る。
また、同様に、信号処理部190は、APD152及び153からの出力に基づいて眼底の平面方向のDOPU画像を生成することもできる。
なお、本明細書において、上述した第一及び第二の偏光に対応する断層画像、リターデーション画像、DOPU画像等を、偏光状態を示す断層画像や偏光断層画像とも言うことにする。また、本明細書において、上述したリターデーションマップや複屈折マップ等を、偏光状態を示す眼底画像や偏光眼底画像とも言うことにする。
[処理動作]
次に本画像処理装置による処理動作について説明する。
図3は、本画像処理装置の処理動作を示すフローチャートである。
<調整>
まず、ステップS101において、被検眼を本装置に配置した状態で、本装置と被検眼のアライメントを行う。アライメントの説明に関して、本実施形態に特有な処理について説明し、ワーキングディスタンス等のXYZ方向のアライメント、フォーカス、コヒーレンスゲートの調整等は一般的であるのでその説明は省略する。
(PS−OCT撮像位置の調整)
図4は、調整時に表示部192に表示されるウィンドウ400を示している。第一の表示領域の一例である表示領域410には、PS−SLO140で撮像され、信号処理部190で生成された眼底画像411が表示され、眼底画像411上に、PS−OCT100の撮像範囲を示す枠412が重畳表示されている。
操作者がマウス等の指示装置(不図示)を用いて、ウィンドウ400に表示されるカーソルで指定し、クリック操作やドラッグ操作等により指示することにより、駆動制御部180の制御の下、撮像範囲の設定が行われる。即ち、カーソルで枠412を指定し、ドラッグ操作することにより、枠412を移動することができる。これにより、駆動制御部180がスキャナの駆動角度を制御することにより撮像範囲を設定する。なお、本実施形態のマウスには、例えば、ユーザの手によってマウス本体が2次元的に移動させたときの移動信号を検出するセンサと、ユーザの手によって押圧されたことを検知するための左右2つのマウスボタンと、左右2つのマウスボタンの間に前後左右に回転可能なホイール機構と、が設けられている。また、指示装置は、表示部にタッチパネル機能を搭載させ、タッチパネル上で取得位置を指定しても良い。
(λ/4偏光板の調整)
λ/4偏光板113の調整について説明する。
図4において、指示部413、414は、λ/4偏光板113の角度を調整するための表示であり、操作者が指示装置を用いて指示することにより、駆動制御部180の制御の下、λ/4偏光板113の角度が調整される。指示部413は反時計回りの調整を、指示部414は時計回りの調整を指示するための表示である。指示部413、414の横に表示されている数値は、現在のλ/4偏光板113の角度を表している。なお、表示制御部191は、λ/4偏光板119の角度を調整する指示部を、指示部413と並べて表示部192に表示させても良いし、指示部413の代わりに表示させても良い。
操作者は、第三の表示領域の一例である表示領域430と、第四の表示領域の一例である表示領域440にそれぞれ表示された各偏光の断層画像の輝度が同じになるように、マウスを用いてカーソルで指示する。なお、各偏光の断層画像431、441と共にピーク輝度値を表示し、あるいは、それぞれの干渉信号の波形そのものを表示し、それを見ながら調整を行う構成でも良い。ここで、各偏光の断層画像431、441は、第一の偏光に対応する断層画像、第二の偏光に対応する断層画像の一例である。なお、各偏光の断層画像431、441(あるいは後述する断層画像531、541)には、それぞれの画像の種類を示す表示形態、例えば、P偏光を示す「P」の文字や、S偏光を示す「S」の文字を画像に重ねて表示させることが好ましい。これにより、ユーザが画像を誤って認識することを防ぐことができる。もちろん、画像に重ねて表示させずに、画像の上側や横側に表示させても良く、画像と対応させるように表示させれば良い。
また、第二の表示領域の一例である表示領域420には、この段階では何も表示させなくても良いし、オート調整等の場合には現在の調整状態を示す表示形態、例えば、「λ/4偏光板の調整中」等のメッセージを表示させても良い。また、ウィンドウ400には、被検眼の左右眼等の患者情報を示す表示形態や、撮影モード等の撮影情報を示す表示形態を表示させても良い。なお、眼底輝度画像と偏光状態を示す断層画像とを交互に取得するように、光路に対するλ/4偏光板113の挿脱を繰り返すことが望ましい。これにより、できるだけ小型な眼科装置において、表示制御部191は、例えば、眼底輝度画像を表示領域410に表示させ、偏光状態を示す断層画像を表示領域420に表示させることができる。
ここで、調整の順番は、前眼部画像や角膜輝点を用いたアライメント調整、偏光状態を示す眼底画像を用いたフォーカス調整、偏光状態を示す断層画像を用いたコヒーレンスゲート調整、λ/4偏光板の調整の順番が好ましい。なお、偏光状態を示す断層画像の取得位置の決定は、偏光状態を示す断層画像を用いたコヒーレンスゲート調整前が好ましいが、偏光状態を示す眼底画像の中心領域を取得するように初期設定で決めるようにしても良い。これにより、偏光状態を示す眼底画像よりも精細で狭い範囲を対象にする偏光状態を示す断層画像を精度良く取得可能に簡単に調整することができる。このとき、コヒーレンスゲート調整の完了に応じてλ/4偏光板を自動的に調整しても良いし、偏光状態を示す画像を取得するための信号の入力に応じてλ/4偏光板を自動的に調整しても良い。もちろん、眼科装置の起動時に初期設定画面等でλ/4偏光板を予め調整しておき、撮影毎に調整しないように構成しても良い。
また、λ/4偏光板を光路に対して挿脱可能に構成している場合、調整の順番は、前眼部画像や角膜輝点を用いたアライメント調整、SLO眼底画像を用いたフォーカス調整、OCT断層画像を用いたコヒーレンスゲート調整、λ/4偏光板を光路に挿入、λ/4偏光板の調整の順番が好ましい。これにより、偏光状態を示す画像の取得前の調整を、ユーザが直感的に慣れている通常のSLO眼底画像やOCT断層画像を用いて行うことができる。ただし、フォーカス調整の後に、λ/4偏光板を挿入してからPS−OCTの偏光状態を示す断層画像を用いたコヒーレンスゲート調整を行っても良い。このとき、コヒーレンスゲート調整の完了あるいはフォーカス調整の完了に応じてλ/4偏光板を自動的に光路に挿入しても良いし、偏光状態を示す画像を取得するための信号の入力に応じてλ/4偏光板を自動的に光路に挿入しても良い。
なお、フォーカス調整は、SLO眼底画像を用いた粗フォーカス調整の後、OCT断層画像を用いた微フォーカス調整を行っても良い。
また、これらの調整は、上記順番で全て自動的に調整しても良いし、表示部に表示された各調整に対応したスライダにカーソルを合わせてドラッグ操作等を行うようにしても良い。また、λ/4偏光板を挿脱する場合、λ/4偏光板を光路に挿入あるいは光路から離脱を指示するためのアイコンを表示部に表示させてもよい。
<撮像>〜<画像生成>
ステップS102〜S103において、光源101、光源141からそれぞれ測定光を出射して、網膜Erからの戻り光を、ラインカメラ129、133、APD152、153で受光して、画像生成部193で、前述の通り各画像を生成する。なお、本実施の形態においては、Xスキャナ110およびYスキャナ107を制御することによって、M本のAスキャンから構成されるN枚のBスキャンが生成されることとする。MおよびNは撮影に要する時間および診断上必要となる領域の大きさから設定できるが、例えば黄斑を中心とした8mm×6mm程度の領域に対してM=1024、N=250程度の値に設定すればよい。
<解析>
ステップS104において、抽出手段の一例である画像解析部194は、(式4)のDOPUを計算し、網膜層における偏光が解消される領域(以下、「偏光解消領域」とも言う)を抽出する。なお、偏光が解消された領域とは、例えば、2つの偏光が被検眼で受ける影響の違いが比較的大きな領域のことである。また、偏光が解消される領域を所定の領域とも呼ぶ。RPEが萎縮などにより乱れた場合、乱れた箇所はDOPUではRPEとして検出されなくなるため、RPEの萎縮箇所を検出することができるというメリットがある。なお、RPEの萎縮箇所は、RPEが不連続な箇所あるいは欠陥箇所(RPEにおける異常箇所)を示している。そして、抽出された偏光が解消された領域の不連続な箇所の状態(連続性等)に基づいて、不連続な箇所の種類を示す表示形態を表示手段に表示させることができる。これにより、偏光が解消された領域における萎縮箇所の状態によって、診断に有効な情報をユーザに提示することができる。
図5に、健常者の黄斑部を撮影した断層画像500、RPE501と、その断層画像上にDOPU計算によって求めた偏光解消領域505を重畳表示した例を示す。図5に示すように、網膜層内においてはRPE501が偏光が解消される領域であるため、RPE501に偏光解消領域505を重畳することになる。
図6に、病変が存在する眼の場合の断層画像と偏光解消領域を重畳した例を示す。図6(a)、(b)は、シュタルガルト病の場合の断層画像600、RPE601、脈絡膜602と、その断層画像上にDOPU計算によって求めた偏光解消領域605を重畳表示した例である。図6(c)、(d)は、加齢黄斑変性で地図上萎縮の断層画像610、RPE611、脈絡膜612と、その断層画像上にDOPU計算によって求めた偏光解消領域615を重畳表示した例である。図6(a)、(c)に示すように、輝度画像である断層画像600と610を比較するだけでは、違いが分かりにくい場合があるが、図6(b)、(d)に示すように、偏光解消領域を検出する事により、それらの違いをみる事が出来る。また、図7には、眼底自発蛍光撮影した場合の画像と、偏光解消領域の厚みを表現したマップの例を示す。偏光解消領域厚みマップは、撮影した複数の断層画像において、図6に示すような偏光解消領域605の深さ方向の厚みを計測して、3次元のデータを2次元のマップとして表現したものである。図7(a)はシュタルガルト病の眼底自発蛍光画像、(b)はその偏光解消領域厚みマップの例である。図7(c)は加齢黄斑変性の眼底自発蛍光画像、(d)はその偏光解消領域厚みマップの例である。図7(a)の701、(c)の711は、RPEが消失、あるいは非薄化してしまっている領域である。図7(b)、(d)に示す偏光解消領域厚みマップは、図7(e)のカラーバー750に示すように、厚みが厚い個所の色を薄く、厚みが薄い個所は色を濃く表示している。図7(a)、(c)に示すように、眼底自発蛍光画像の比較では、違いが分かりにくい場合があるが、図7(b)、(d)に示すように、偏光解消領域を厚みマップとして表現する事でその違いをみる事が出来る。
決定手段の一例である画像解析部194は、上述した偏光解消領域の不連続な状態(連続性等)に基づいて、不連続な箇所の種類を決定する。例えば、偏光解消領域の厚みマップにおいて、偏光解消領域の厚みが一定値以上(50μm以上と仮定する)となる領域の分布の幅が閾値より大きい(10%以上と仮定する)場合、シュタルガルト病の可能性がある。また、偏光解消領域の厚みが一定値以下(5μm以下と仮定する)となる領域の分布の幅が閾値より小さい場合、或いは、偏光解消領域が存在しない領域の分布の幅が閾値より大きい(10%以上と仮定する)場合、加齢黄斑変性の地図上萎縮の可能性がある。このように、深さ方向に対して交差する方向における不連続な箇所の分布の幅の閾値に対する大小に応じて、不連続な箇所における病変の種類を不連続な箇所の種類として決定することができる。あるいは、画像解析部194は、Support vector MachineやAdaboostといった識別器を持ち、この識別器に断層画像、偏光解消領域の特徴を用いて病気の種類を決定してもよい。
ここで、シュタルガルト病 (遺伝子変異による網膜変性)や、加齢黄斑変性の地図状萎縮等の異なる病気において、通常のOCT画像(断層輝度画像)では病気の違いが良く分からない場合がある。
一方、偏光断層画像では、病気の違いが分かる場合がある。そこで、偏光断層画像を用いて自動的に病気を判断する。
具体的には、偏光断層画像の所定の領域(RPE層等の偏光が解消される領域)における不連続な箇所の状態に基づいて、病気を自動識別し、病気の種類(病名)をユーザに提示する。このように、偏光断層画像を用いることにより、ユーザは、病気によってRPE層等に生じる微小な変化を早期に発見し、その病気の種類を早期に知ることができる。
なお、不連続な箇所の状態とは、例えば、RPE層の厚さ、幅、大きさ、面積、隣り合う画素同士の連結性等であり、不連続性を示すパラメータであれば何でも良い。また、不連続な箇所の状態とは、3次元の偏光断層画像におけるRPE層において、不連続な箇所の体積でも良い。そして、不連続性を示すパラメータと閾値とを比較することにより、比較結果に応じて病気の種類を決定することができる。また、RPE層における不連続な箇所の厚さや幅等の大小を比較する基準を断層輝度画像から抽出しても良い。例えば、断層輝度画像から輝度値が大きい領域であるが、偏光が解消しない領域(例えば、ブルッフ膜)を抽出して、該抽出された領域を該基準としても良い。
<出力>
次に、生成した各画像及び解析した結果の出力処理ステップS105について説明する。本実施形態おける出力処理は、ステップS104で求めた偏光解消領域の特徴を効果的に表示する。
信号処理部190内の画像生成部193と画像解析部194において、各画像の生成及び解析が終了すると、その結果に基づき、表示制御部191は、出力情報を生成し、表示部192に出力して表示を行う。
図8は、本実施形態における表示部192における表示例である。図において、800は表示部192に表示されるウィンドウであり、表示領域810、820、830、840を有する。
第一の表示領域である表示領域810には、眼底画像811が表示され、断層画像の撮影位置を示す矩形の枠812と断層画像600の位置を示す線813が重畳されている。ここで眼底画像811は、PS−SLOで撮影されたものである。眼底画像811の付近には、ステップS104で求めた偏光解消領域の特徴を用いて推定をした病気の種類を表示する表示領域815がある。本実施形態においては、シュタルガルト病(Stargardt)の場合を例に示す。図の例では、推定した病変名を表示するのみであるが、これに限らず、確率の高い病変の種類を確率順に3種類程度表示するようにしてもよい。その場合は、病変名だけではなく、その確率も表示する事が望ましい。そのような表示を行う場合、識別器を用いて病変名を推定した結果を表示する。
第二の表示領域である表示領域820には、断層画像600が表示される。更に、表示領域820には、表示される断層画像の種類を選択するための選択部の一例であるボタン822〜825が表示される。なお、ボタン822〜825の代わりにメニューから断層画像の種類を選択するようにしても良い。図8においては、ボタン825が選択された状態を示しており、断層輝度画像の上にセグメンテーション結果を重畳表示し、偏光解消領域605が強調表示されている例を示している。その他のボタン822〜824とその表示について説明をする。操作者が、ボタン822を選択すると、表示領域820に断層輝度画像を表示する。ボタン823を選択すると、表示領域820にリターデーション画像を表示する。ボタン824を選択すると、表示領域820にDOPU画像を表示する。
第三の表示領域の一例である表示領域830には、偏光解消領域の厚みマップ705を表示する。偏光解消領域の厚みマップは図7で説明したものである。そして、偏光解消領域の厚みマップ705上に断層画像600の位置を示す線831を重畳している。これは、813とも対応している。また、第四の表示領域の一例である表示領域840には、偏光解消領域の厚みを示すグラフ841を表示する。これにより、ユーザは、PS−OCTにおける構造変化と網膜層の構造とを対応付けて確認することができるので、被検眼の診断を効果的に行うことができる。
なお、表示制御部191は、上述した各種画像の代わりに、リターデーションマップや複屈折マップを表示部192の各表示領域のいずれかに表示させても良い。また、表示制御部191は、リターデーションマップや複屈折マップを眼底輝度画像811に重ねて表示させても良い。このとき、リターデーションマップや複屈折マップを枠812で示される領域に重ねて表示させることが好ましい。以上説明のように本実施形態によれば、偏光OCT画像から得た偏光成分の情報を用いて、偏光解消物質の情報を表示するとともに、病変の可能性を推定した情報を提供することが出来る。
なお、これらの画像を表示させる表示領域の位置は、本実施形態に限定されず、例えば、眼底画像を表示画面の左側の表示領域に表示させても良い。また、表示させる画像の数も本実施形態に限定されず、例えば、調整時には眼底画像と断層画像との合計2つを表示画面に並べて表示し、撮影後には表示方式を変更し、眼底画像の他に、それぞれ異なる偏光状態を示す複数の断層画像を表示画面に並べて表示させても良い。また、ボタン822〜825が並ぶ順番や位置等も本実施形態に限らない。
(特定の病気において分類された種類の決定)
次に、不連続な箇所について、特定の病気、例えば、シュタルガルト病について、さらに細かく分類された種類を決定する例について、図9を用いて説明する。なお、本例も、不連続な箇所の種類を決定する一例である。
まず、図9(a)、(b)は、シュタルガルト病で視細胞内節外節接合部(以下IS/OS)が消失しRPEにも萎縮または欠損が生じている場合の断層画像900、RPE901、IS/OS902と、その断層画像上にDOPU計算によって求めた偏光解消領域905を重畳表示した例である。また、破線の円領域903は、偏光解消領域の不連続な箇所を示す表示形態の一例であり、IS/OSが消失してしまっている領域を示している。図9(b)に示すように、IS/OSが消失してしまっている領域においては、偏光解消領域905は、RPEと思われる領域には存在せずに、RPEが萎縮してしまっている事が把握出来る。すなわち、輝度画像である断層画像900からでは判断しにくいRPEの萎縮が偏光解消領域を求めることで把握しやすくなる。
次に、図9(c)、(d)は、シュタルガルト病でRPEが内網膜層にまで向かって浸潤している場合の断層画像910、RPE911、IS/OS912と、その断層画像上にDOPU計算によって求めた偏光解消領域915を重畳表示した例である。913および914は、RPEに対して病変が存在するために、偏光解消領域のうち内網膜層に浸潤している箇所である。図9(d)において、偏光解消領域のうちIS/OSまで達している箇所913と、IS/OSの下に留まっている箇所914が示されている。IS/OSは視細胞の接合部であり、そこに病変部が存在すると視覚障害に繋がる可能性が考えられる。
決定手段の一例である画像解析部194は、上述した偏光解消領域の不連続な箇所の状態に基づいて、不連続な箇所の病気の種類を決定する。例えば、画像解析部194は、網膜層境界を検出し、境界の位置と偏光解消領域との位置関係、あるいは、偏光解消領域の連続性を見ることで病変の種類を推定する。網膜層境界に関して、RPEやIS/OSの境界を輝度画像である断層画像900、910から検出をする。境界の検出方法としては、例えば、処理の対象とする断層画像に対して、平滑化の一種としてメディアンフィルタと、エッジ検出の一種としてSobelフィルタとをそれぞれ適用して、それぞれの画像(以下、「メディアン画像」、「Sobel画像」ともいう)を作成する。次に、作成したメディアン画像とSobel画像から、Aスキャン毎にプロファイルを作成する。メディアン画像では輝度値のプロファイル、Sobel画像では勾配のプロファイルとなる。そして、Sobel画像から作成したプロファイル内のピークを検出する。検出したピークの前後やピーク間に対応するメディアン画像のプロファイルを参照しそれらの境界を検出する。検出した境界はそのまま用いてもよいし、検出した境界を基にN次曲線をあてはめて、基準線を推定してもよい。本実施形態においては、画像解析部194で検出したRPE境界に対して2次曲線をあてはめて基準線とする。2次曲線の推定には、M推定やLMedSなどのロバスト推定を用いる。これらの境界の位置と偏光解消領域との位置関係を利用する。
画像解析部194は、RPE境界あるいは基準線に対して、偏光解消領域が連続しているか、不連続な範囲において、偏光解消領域は、RPE境界あるいは基準線に対してどの位置に存在しているかという特徴を用いて推定をする。シュタルガルト病でIS/OSが消失しRPEにも萎縮または欠損が生じている場合は、RPEが萎縮している範囲において、偏光解消領域の不連続性と、RPE境界の下方(深度方向)に偏光解消領域が存在している。また、シュタルガルト病でRPEが内網膜層にまで向かって浸潤している場合は、偏光解消領域の厚みが一部厚い個所が存在し、その領域は、IS/OS境界にまで達するか、IS/OS境界付近にまで達している。さらに、上述したシュタルガルト病で脈絡膜に偏光解消物質が拡散している場合、RPE境界の下方(深度方向)に偏光解消領域の厚みが一定値以上(50μm以上と仮定する)となる領域が広範囲に存在することとなる。従って、画像解析部194は、これらの特徴を用いて種類の違いを推定する事が出来る。
(その他の実施例)
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。

Claims (17)

  1. 被検体の偏光断層画像を取得する断層画像取得手段と、
    前記取得された偏光断層画像から所定の領域を抽出する抽出手段と、
    前記抽出された所定の領域の不連続な箇所の状態に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させる表示制御手段と、
    を有することを特徴とする画像処理装置。
  2. 前記不連続な箇所の連続性に基づいて、前記不連続な箇所の種類を決定する決定手段を有し、
    前記表示制御手段が、前記決定された種類を示す表示形態を前記表示手段に表示させることを特徴とする請求項1に記載の画像処理装置。
  3. 前記決定手段が、深さ方向に対して交差する方向における前記不連続な箇所の分布の幅の閾値に対する大小に応じて、前記不連続な箇所における病変の種類を前記不連続な箇所の種類として決定することを特徴とする請求項2に記載の画像処理装置。
  4. 前記表示制御手段は、前記不連続な箇所の不連続性を示すパラメータに基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させることを特徴とする請求項1乃至3のいずれか1項に記載の画像処理装置。
  5. 前記表示制御手段は、前記不連続な箇所の厚さと幅とのうち少なくとも一方に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させることを特徴とする請求項1乃至4のいずれか1項に記載の画像処理装置。
  6. 前記種類は、病名であることを特徴であることを特徴とする請求項1乃至5のいずれか1項に記載の画像処理装置。
  7. 前記所定の領域は、前記偏光が解消された領域であり、
    前記表示制御手段が、前記偏光が解消された領域を前記被検体の断層輝度画像に重ねて前記表示手段に表示させ、且つ前記不連続な箇所を示す表示形態を前記断層輝度画像に重ねて前記表示手段に表示させることを特徴とする請求項1乃至6のいずれか1項に記載の画像処理装置。
  8. 前記偏光断層画像は、測定光を照射した前記被検体からの戻り光と、該測定光に対応する参照光とを合波した光を分割して得た互いに異なる偏光の光に基づいて、生成された断層画像であることを特徴とする請求項1乃至7のいずれか1項に記載の画像処理装置。
  9. 被検体の偏光断層画像を取得する工程と、
    前記取得された偏光断層画像から所定の領域を抽出する工程と、
    前記抽出された所定の領域の不連続な箇所の状態に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させる工程と、
    を有することを特徴とする画像処理方法。
  10. 前記不連続な箇所の連続性に基づいて、前記不連続な箇所の種類を決定する工程を有し、
    前記表示させる工程では、前記決定された種類を示す表示形態を前記表示手段に表示させることを特徴とする請求項9に記載の画像処理方法。
  11. 前記決定する工程では、深さ方向に対して交差する方向における前記不連続な箇所の分布の幅の閾値に対する大小に応じて、前記不連続な箇所における病変の種類を前記不連続な箇所の種類として決定することを特徴とする請求項10に記載の画像処理方法。
  12. 前記表示させる工程では、前記不連続な箇所の不連続性を示すパラメータに基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させることを特徴とする請求項9乃至11のいずれか1項に記載の画像処理方法。
  13. 前記表示させる工程では、前記不連続な箇所の厚さと幅とのうち少なくとも一方に基づいて、前記不連続な箇所の種類を示す表示形態を表示手段に表示させることを特徴とする請求項9乃至12のいずれか1項に記載の画像処理方法。
  14. 前記種類は、病名であることを特徴であることを特徴とする請求項9乃至13のいずれか1項に記載の画像処理方法。
  15. 前記所定の領域は、前記偏光が解消された領域であり、
    前記表示させる工程では、前記偏光が解消された領域を前記被検体の断層輝度画像に重ねて前記表示手段に表示させ、且つ前記不連続な箇所を示す表示形態を前記断層輝度画像に重ねて前記表示手段に表示させることを特徴とする請求項9乃至14のいずれか1項に記載の画像処理方法。
  16. 前記偏光断層画像は、測定光を照射した前記被検体からの戻り光と、該測定光に対応する参照光とを合波した光を分割して得た互いに異なる偏光の光に基づいて、生成された断層画像であることを特徴とする請求項9乃至15のいずれか1項に記載の画像処理方法。
  17. 請求項9乃至16のいずれか1項に記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
JP2013159180A 2012-10-30 2013-07-31 画像処理装置及び画像処理方法 Pending JP2014110883A (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2013159180A JP2014110883A (ja) 2012-10-30 2013-07-31 画像処理装置及び画像処理方法
US14/065,759 US20140121506A1 (en) 2012-10-30 2013-10-29 Image processing apparatus and image processing method
EP13190930.1A EP2728549B1 (en) 2012-10-30 2013-10-30 Image processing apparatus and image processing method

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012239273 2012-10-30
JP2012239273 2012-10-30
JP2013159180A JP2014110883A (ja) 2012-10-30 2013-07-31 画像処理装置及び画像処理方法

Publications (2)

Publication Number Publication Date
JP2014110883A true JP2014110883A (ja) 2014-06-19
JP2014110883A5 JP2014110883A5 (ja) 2016-09-08

Family

ID=49554012

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013159180A Pending JP2014110883A (ja) 2012-10-30 2013-07-31 画像処理装置及び画像処理方法

Country Status (3)

Country Link
US (1) US20140121506A1 (ja)
EP (1) EP2728549B1 (ja)
JP (1) JP2014110883A (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016107148A (ja) * 2014-01-21 2016-06-20 株式会社トプコン 眼科画像処理装置
JP2016129664A (ja) * 2015-01-07 2016-07-21 キヤノン株式会社 画像処理装置及び画像処理方法
JP2016131881A (ja) * 2015-01-15 2016-07-25 株式会社トプコン 地図状萎縮の同定および測定
WO2020116351A1 (ja) * 2018-12-04 2020-06-11 株式会社ニデック 診断支援装置、および診断支援プログラム
JP2020137103A (ja) * 2019-02-19 2020-08-31 株式会社シンクアウト 情報処理システムおよび情報処理プログラム

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6146951B2 (ja) 2012-01-20 2017-06-14 キヤノン株式会社 画像処理装置、画像処理方法、撮影装置及び撮影方法
JP6061554B2 (ja) * 2012-01-20 2017-01-18 キヤノン株式会社 画像処理装置及び画像処理方法
JP6184114B2 (ja) * 2013-01-31 2017-08-23 キヤノン株式会社 光干渉断層撮像装置およびその制御方法
CA2882784C (en) 2014-11-05 2021-02-02 National Taiwan University Three-dimensional optical coherence tomography apparatus and its application
TWI553294B (zh) * 2014-11-05 2016-10-11 Univ Nat Taiwan 干涉式光學成像裝置、其應用之系統及方法
JP2016093240A (ja) 2014-11-12 2016-05-26 キヤノン株式会社 画像処理装置及び画像処理方法
WO2016110917A1 (en) * 2015-01-07 2016-07-14 Canon Kabushiki Kaisha Image processing apparatus and image processing method for polarization-sensitive optical coherence tomography
JP2017080344A (ja) * 2015-10-30 2017-05-18 キヤノン株式会社 画像処理装置、画像処理方法及び光干渉断層撮影装置
JP6748434B2 (ja) * 2016-01-18 2020-09-02 キヤノン株式会社 画像処理装置、推定方法、システム及びプログラム
JP2022047278A (ja) * 2020-09-11 2022-03-24 株式会社トーメーコーポレーション 光断層画像撮影装置及び断層画像を表示するためのコンピュータプログラム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006077885A1 (ja) * 2005-01-19 2006-07-27 Hitachi Medical Corporation 画像診断支援装置及び画像診断支援プログラム
JP2010094381A (ja) * 2008-10-17 2010-04-30 Canon Inc 画像処理装置、画像処理方法
EP2243420A1 (en) * 2009-04-24 2010-10-27 Schmidt-Erfurth, Ursula Method for determining exudates in the retina
JP2012016525A (ja) * 2010-07-09 2012-01-26 Canon Inc 光断層画像撮像装置及びその撮像方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6522407B2 (en) * 1999-01-22 2003-02-18 The Regents Of The University Of California Optical detection dental disease using polarized light
US7612880B2 (en) * 2003-08-06 2009-11-03 Arizona Board Of Regents On Behalf Of The University Of Arizona Advanced polarization imaging method, apparatus, and computer program product for retinal imaging, liquid crystal testing, active remote sensing, and other applications
US7668342B2 (en) * 2005-09-09 2010-02-23 Carl Zeiss Meditec, Inc. Method of bioimage data processing for revealing more meaningful anatomic features of diseased tissues
US7593559B2 (en) * 2005-11-18 2009-09-22 Duke University Method and system of coregistrating optical coherence tomography (OCT) with other clinical tests
JP4822454B2 (ja) * 2005-11-22 2011-11-24 株式会社松風 歯科用光コヒーレンストモグラフィー装置
US8223143B2 (en) * 2006-10-27 2012-07-17 Carl Zeiss Meditec, Inc. User interface for efficiently displaying relevant OCT imaging data
JP5448353B2 (ja) * 2007-05-02 2014-03-19 キヤノン株式会社 光干渉断層計を用いた画像形成方法、及び光干渉断層装置
US8208996B2 (en) * 2008-03-24 2012-06-26 Carl Zeiss Meditec, Inc. Imaging of polarization scrambling tissue
JP4909378B2 (ja) * 2009-06-02 2012-04-04 キヤノン株式会社 画像処理装置及びその制御方法、コンピュータプログラム
JP5701024B2 (ja) * 2010-11-26 2015-04-15 キヤノン株式会社 画像処理装置及び方法
CN103477205B (zh) * 2011-01-31 2016-08-10 斯凯德医疗有限公司 光学偏振成像
EP2574273B1 (en) * 2011-06-23 2014-09-24 Nidek Co., Ltd. Optical coherence tomography apparatus
JP5988772B2 (ja) * 2012-01-20 2016-09-07 キヤノン株式会社 画像処理装置及び画像処理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006077885A1 (ja) * 2005-01-19 2006-07-27 Hitachi Medical Corporation 画像診断支援装置及び画像診断支援プログラム
JP2010094381A (ja) * 2008-10-17 2010-04-30 Canon Inc 画像処理装置、画像処理方法
EP2243420A1 (en) * 2009-04-24 2010-10-27 Schmidt-Erfurth, Ursula Method for determining exudates in the retina
JP2012016525A (ja) * 2010-07-09 2012-01-26 Canon Inc 光断層画像撮像装置及びその撮像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHRISTIAN AHLERS, ET AL.: ""Imaging of the Retinal Pigment Epithelium in Age-Related Macular Degeneration Using Polarization-S", INVESTIGATIVE OPHTHALMOLOGY AND VISUAL SCIENCE, vol. 51, no. 4, JPN6017020327, April 2010 (2010-04-01), pages 2149 - 2157, ISSN: 0003717173 *
MICHAEL PIRCHER, ET AL.: ""Polarization sensitive optical coherence tomography in the human eye"", PROGRESS IN RETINAL AND EYE RESEARCH, vol. Vol.30 Issue.6, JPN6017020326, November 2011 (2011-11-01), pages 431 - 451, ISSN: 0003871002 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016107148A (ja) * 2014-01-21 2016-06-20 株式会社トプコン 眼科画像処理装置
JP2016129664A (ja) * 2015-01-07 2016-07-21 キヤノン株式会社 画像処理装置及び画像処理方法
JP2016129663A (ja) * 2015-01-07 2016-07-21 キヤノン株式会社 画像処理装置及び画像処理方法
JP2016131881A (ja) * 2015-01-15 2016-07-25 株式会社トプコン 地図状萎縮の同定および測定
WO2020116351A1 (ja) * 2018-12-04 2020-06-11 株式会社ニデック 診断支援装置、および診断支援プログラム
JPWO2020116351A1 (ja) * 2018-12-04 2021-10-28 株式会社ニデック 診断支援装置、および診断支援プログラム
JP2020137103A (ja) * 2019-02-19 2020-08-31 株式会社シンクアウト 情報処理システムおよび情報処理プログラム
JP7304508B2 (ja) 2019-02-19 2023-07-07 株式会社シンクアウト 情報処理システムおよび情報処理プログラム

Also Published As

Publication number Publication date
EP2728549B1 (en) 2018-12-12
EP2728549A1 (en) 2014-05-07
US20140121506A1 (en) 2014-05-01

Similar Documents

Publication Publication Date Title
JP5988772B2 (ja) 画像処理装置及び画像処理方法
JP6071331B2 (ja) 画像処理装置及び画像処理方法
JP6061554B2 (ja) 画像処理装置及び画像処理方法
JP2014110883A (ja) 画像処理装置及び画像処理方法
JP6184232B2 (ja) 画像処理装置及び画像処理方法
JP6061555B2 (ja) 画像処理装置及び画像処理方法
JP6184231B2 (ja) 画像処理装置および画像処理方法
JP6202924B2 (ja) 撮影装置及び撮影方法
JP2014110884A (ja) 画像処理装置及び画像処理方法
JP2016093240A (ja) 画像処理装置及び画像処理方法
JP2017131550A (ja) 画像処理装置及び画像処理方法
EP2617352B1 (en) Imaging apparatus, imaging method, and program
JP2014083285A (ja) 画像処理装置及び画像処理方法
JP6146951B2 (ja) 画像処理装置、画像処理方法、撮影装置及び撮影方法
JP5988883B2 (ja) 画像処理装置及び画像処理方法
JP6381622B2 (ja) 画像処理装置及び画像処理方法
JP6505072B2 (ja) 画像処理装置及び画像処理方法
JP6437055B2 (ja) 画像処理装置及び画像処理方法
JP2013148509A (ja) 画像処理装置及び画像処理方法
JP6039185B2 (ja) 撮影装置
JP5936368B2 (ja) 光干渉断層撮影装置及びその作動方法
JP2019111396A (ja) 画像処理装置及び画像処理方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160721

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160721

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170606

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170807

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180116

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180315

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20180904