JP2013531262A - 色反射率離散化分析を用いて光コヒーレンストモグラフィー(oct)の結果を分析する方法およびシステム - Google Patents

色反射率離散化分析を用いて光コヒーレンストモグラフィー(oct)の結果を分析する方法およびシステム Download PDF

Info

Publication number
JP2013531262A
JP2013531262A JP2013520870A JP2013520870A JP2013531262A JP 2013531262 A JP2013531262 A JP 2013531262A JP 2013520870 A JP2013520870 A JP 2013520870A JP 2013520870 A JP2013520870 A JP 2013520870A JP 2013531262 A JP2013531262 A JP 2013531262A
Authority
JP
Japan
Prior art keywords
points
oct image
group
reflective layer
oct
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
JP2013520870A
Other languages
English (en)
Other versions
JP2013531262A5 (ja
JP5884195B2 (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 ディオプシス、インコーポレーテッド
Publication of JP2013531262A publication Critical patent/JP2013531262A/ja
Publication of JP2013531262A5 publication Critical patent/JP2013531262A5/ja
Application granted granted Critical
Publication of JP5884195B2 publication Critical patent/JP5884195B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2433Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • 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]
    • 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/10024Color image
    • 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
    • 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/30242Counting objects in image

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

【解決手段】 さまざまな色の点を複数含むOCT画像を分析する方法を記載する。システムは、さまざまな色の点の数を検出して、点を数値化して、統計学的分析を実行して、疾患が発生している可能性を求める。さまざまな色の点は色毎に、異なる種類の網膜細胞および網膜構造に対応している。
【選択図】 図1

Description

本発明は、光コヒーレンストモグラフィー(OCT)に関する。特に、OCT画像の分析および解釈を行う新型システムに関する。
<関連出願>
本国際出願は、米国仮特許出願第61/366,485号(出願日:2010年7月21日)による恩恵を主張する。当該米国特許出願の内容は、参照により本願に組み込まれる。
OCTは、干渉および分光分析を利用した方法であり、生物組織のマイクロメートルレベルの高分解能画像を取得する方法として医療分野で一般的に利用されている。OCT機器は、組織に光の焦点を合わせ、反射光を検出して分析する。組織からの反射光および基準光源を組み合わせることで、干渉およびスペクトルのパターンが得られる。このような処理を用いることで、組織の2次元画像および3次元画像の両方が得られる。
OCTは、どのような生物組織または組織標本でも撮像するために利用され得る。眼科では、OCTは、緑内障が存在するか否かを評価することを目的として網膜を撮像するために用いられるのが一般的である。緑内障を検出するべく、神経節細胞、グリア細胞および他の支持細胞の軸索を含む網膜神経線維束(RNFL)を分離して測定する。RNFL損傷は通常、視野欠損より先に見られるので、RNFL損傷を検出することで、緑内障またはその他の神経系の疾患の初期徴候を発見できる場合がある。神経細胞が失われると、RNFLが薄くなり、その結果としてRNFLが薄くなった1以上の領域が緑内障が発症していることを示すとしてよい。
現在の技術水準では、RNFLの厚みが、緑内障を評価する臨床医が観察する唯一の特性である。しかし、RNFLの内部構造を観察または考慮することなくRNFLの厚みを検査することは、誤った結果を得ることになったり、および/または、緑内障の発見精度が低くなったりする原因となり得る。診断で用いる唯一の測定値が厚みである場合、厚み検出アルゴリズムで誤りがあれば、結果も誤りを含むことになる。また、局所的に薄くなっている領域は、全体的なRNFL厚み測定では平均されてしまうので、見落とされる可能性がある。さらに、軸索の束の局所欠陥は、疾患の初期段階、当該欠陥によりRNFLが薄くなったことが検出可能になる前では見落とされる可能性がある。
本発明では、従来のOCT画像を利用する。従来のOCT画像は他にも、これまで医者によって認識または利用されてこなかったが、緑内障の検出精度を向上させるべく、従来のOCTとは別に、または、従来のOCTを補佐する形で、利用され得る診断上重要な情報を含む。
本発明は、網膜の画像全体に分散した点を数値化する特別なソフトウェアを備えたコンピュータを利用する新型システムに関する。点毎に表す生体細胞の種類が異なり、点の絶対数および割合を数値化することによって、OCT画像から多くの診断情報が追加で得られる。
また、本発明に係るシステムでは、径方向スキャンで点の離散化および数値化を行うことで、視神経(篩板)の分析を実行する。軸索細胞および支持結合組織を表す点をカウントして、緑内障に一致する状態として反射率の低下が見られるか否かを判断する。
通常のOCT画像は、上側のRNFLおよび下型のRNFLに分散している少なくとも4つの異なる色の点を持つことを特徴とする。赤色、黄色、緑色、および青色の点がある。赤色および黄色の点(これに代えて、本明細書では「ホット点」と呼ぶ)は、軸索および神経節細胞を表し、緑色および青色の点(これに代えて、本明細書では「コールド点」と呼ぶ)はグリア細胞および他の支持細胞を表す。
健常なRNFLでは、支持細胞よりも神経細胞の方が多く(このため、コールド点よりもホット点の方が多い)、軸索は、光源に対して垂直に配向されているので、他の種類の細胞よりも反射率が高い。しかし、緑内障という疾患の進行では、軸索が失われ、OCT画像でのホット点が減少することになり、ホット点およびコールド点の差が小さくなる。緑内障を患っている目では、ホット点とコールド点との比率が1に近くなる場合がある。
本発明に係るシステムは、各色の点の絶対数を数値化する特別なソフトウェアを利用して、ホット点とコールド点との比率を算出し、当該比率を分析してRNFLが健常であるか否かの指標を生成する。この比率が1に近いほど、検査対象の目が緑内障を患っている可能性が高くなる。
この比率を決定する処理は、OCT画像全体でも行われるし、セクタ毎でも実行される。
比率分析に加えて、各色の点の絶対数を、データベースに格納されている「通常」絶対数と比較して、罹患可能性を判断する。「通常」数は、健常な患者のOCT画像を検査して、通常時の点の数および分布を判断することによって導出される。絶対数の比較もまた、OCT画像全体でも、セクタ毎でも実行される。
同様の分析は、従来のOCT撮像技術を利用して得られる視神経の径方向スキャン結果および平行スキャン結果に対して行われる。ホット点およびコールド点をカウントして、ホット点およびコールド点の絶対数およびホット点/コールド点の比率を求める。このような分析は、径方向および平行方向のスライス毎に、そして、全部をまとめると視神経を表す複数のスライスを組み合わせたものに対して実行される。全体的な分析およびスライス毎の分析は、点の絶対数と基準データベースとの比較、および、ホット点/コールド点の比率の分析を含む。
本発明の実施形態に係る、局所欠陥を含むRNFLの一部のOCT画像を示す概略図である。 局所欠陥を強調している、図1の概略図である。 本発明の実施形態に係る、OCT画像の上側反射領域および下側反射領域を示す概略図である。 本発明の実施形態に係る、上側反射領域を示す抽出図である。 本発明の実施形態に係る、RNFLの一部における局所欠陥を示す概略拡大図である。 本発明の実施形態に係るOCT画像の点分析を実行するためのステップを示すフローチャートである。 本発明の実施形態に係るOCT画像の点分析を実行するためのステップを示すフローチャートである。
上述した図面を参照しつつ本発明の実施形態を以下で説明する。しかし、本発明の本明細書の図面および説明は、本発明の範囲を限定することを意図したものではない。本発明について本明細書で記載した内容は、本発明の意図から逸脱することなく、さまざまな点で変形することが可能であると理解されたい。また、本明細書で説明する特徴を省略するとしてもよい、特徴を追加するとしてもよいし、および/または、本明細書で説明する特徴を、本明細書で記載した特定の組み合わせとは異なる方法で組み合わせるとしてもよいし、これらは全て本発明の意図から逸脱しないものとする。
上述したように、本発明は、RNFLの内部構造を評価し、数値化分析を実行するための方法および新型ソフトウェアプログラムに関する。この目的は、従来の未加工のOCTデータおよび/またはOCT画像に対してコンピュータで実行する新規のプロセスおよび統計学的分析によって実現する。
図1は、OCT検査結果12を示す概略図である。上述したように、現在の技術水準では、RNFLの厚みを測定して、疾患が発生している可能性を判断する。通常のOCT分析は、上側反射層の上側境界と下側境界との間の距離を測定してRNFLが薄くなっているか否かを検出するために行われる。
しかし、このような分析では、内部欠陥を検査または検出できない。実際、内部欠陥はRNFLが薄くなるよりも前に発生する場合があるので、このような欠陥を検出することで、眼疾患をより早期に検出することが可能になる。
図2は、局所欠陥の一部を強調している図1のOCT画像である。図示しているように、RNFLには色の無い領域14が数多く点在している。これらは、反射率がゼロまたは低くなった領域である。上述したように、反射率が低くなるということは、網膜血管または網膜細胞が欠損していることを意味するとしてよい。このため、緑内障等の眼疾患をできる限り早い段階で検出するのであれば、このような局所欠陥を観察することが重要になる。また、本発明に係るシステムおよび方法では、同様に疾患の早期徴候を検出するべく、OCT画像内の特定の細胞の種類および量を分析して解釈する。
本発明に係る分析の最初のステップは、RNFLを残りのOCT画像から分離して、RNFLに対して分析を実行できるようにする。RNFL領域がこのように分離されると、本発明に係るシステムは、以下で説明するように、構造の分析を実行する。
本発明の一実施形態によると、RNFLは以下のようにして分離する。未加工データ、JPEG、TIFF等の任意のデジタル形式で格納されているOCT画像12から、本来の画像の色を、グレースケール画像に変換する。ソフトウェアは、グレースケール画像を構成しているピクセル値を測定して、平均グレースケール階調を決定する。平均階調をしきい値に設定して、画像のピクセルを二値化して、黒または白とする。しきい値よりも高い値の画像ピクセルは白に変換して、しきい値よりも低い値の画像ピクセルは、黒に変換する。コンピュータは、縦列毎にOCT画像全体に対して二値化プロセスを実行する。このような二値化を行うと、画像の略中点に細い黒の横列16が生じる。この細い黒の横列のすぐ上方にある画像部分18を、RNFLの分析で利用する。
一例として、<100111100011100010>が1つの縦列を表すと仮定し、0が黒い点を表し、1が白い点を表す。この縦列は、<1、−2、4、−3、3、−3、1、−1>とグループ分け(グループ分けA)する。負の数は連続した黒の点の量を表し、正の数は、連続した白の点の量を表す。このグループを降順に並び変えて、<4、3、1、1、−1、−2、−3、−3>とする(グループ分けB)。最初の数字および2番目の数字(4および3)は、連続する白の横列から成るグループのうち最大のグループの2つを表している。ソフトウェアは、グループ分けBの最初の2つの数字(連続した白の横列から成るグループのうち最大のグループ2つを表している、本例では4および3)を取り出して、グループ分けAでのそれぞれの位置を特定するように設定されている。グループ分けAで4および3が特定されると、ソフトウェアは、4と3との間にある数字(つまり、グループ分けAにおいて、グループ分けBで特定された最も大きい2つの正の数の間にあった、数字)を特定する。この例では、4と3とにはさまれているのは−3であることが分かる。−3で表されている3つの連続した黒の横列のグループが、中央の黒の点と判断され、この縦列を分割するものとして算出される。この点を「中点」と呼ぶ。
ソフトウェアは、中点16よりも下方の下側画像部分20は無視して、RNFL分析については中点16よりも上方の部分のみ分析するように構成されている。上述した例のグループ分けAについて、−3は、中点を表し、<1、−2、4>は、RNFL分析を行う領域を表す(これらの数字は、縦列の上端から始まるものと仮定される)。<3、−3、1、−1>と表されるグループが、中点より下方であり、無視される。
本発明の別の実施形態によると、画像認識技術を利用して、分析のためにRNFL領域を分離する。本実施形態によると、RNFL領域を抽出する演算は、例えば、連結成分抽出(connected component labeling)、ブロブ(blob)抽出、領域ラベリング等の演算を用いて実現する。このような演算子を適用する前に、1以上のデジタルフィルタリング処理、例えば、メディアンフィルタ等を利用して、ノイズを低減するとしてよい。エロージョン(erosion)およびダイレーション(dilation)といった他の演算を利用するとしてもよい。
好ましい実施形態では、連結成分抽出フィルタ(CCL)を利用して、OCT画像をスキャンして、ピクセルをグループ分けして、ピクセルの連結性に基づいて複数の成分に分割する。全てのグループ分けを決定すると、各ピクセルを、割り当てられている成分(例えば、ブロブ)に応じてラベリングする。
図3を参照すると、OCT画像には、上側反射層22および下側反射層24を規定している2つの主要ブロブが見られる。CCLフィルタは、n個のブロブを生成し、システムでは、最も大きいk個のブロブを利用する。尚、k<<nであり、k>2である。k個のブロブは、最大値から最小値へと並び替えられる。
上側反射層、つまり、RNFL分析を行う領域を分離する場合、最小サイズの最初のブロブ(上から下へとスキャンする)が上側層22を規定する。ブロブを上側層のフットプリントとして用いることによって、システムは、ブロブを横切るように横方向に移動して、全てのギャップを埋める。
下側反射層(RNFL分析では利用されない)を決定するべく、システムは、最大ブロブ22よりも緯度が低い、残っている2番目に大きいブロブ24を認識する。上側層と同様に、ブロブを横切るように横方向に移動してギャップを埋める。
両方のブロブ(つまり、上側反射領域22および下側反射領域24)をこのように特定した後、上側ブロブ22の上部境界26および下側ブロブ24の上部境界28を規定する点のリストを生成する。この点のリストは、上側反射層を規定し、構造分析の領域として指定される。上側反射層より上方の部分は無視する。図4は、上側反射層を示す抽出画像である。
この点のリストは、最小二乗法を利用して、最良近似多項式を算出することによって平滑化されると理解されたい。さらに、多項式はわずかに増加させて(所定の定数アルゴリズムパラメータを用いて)上側層の一部を含まないようにするとしてよい。
OCT画像の上側反射領域は、複数の点で構成される(図5の拡大図が最も分かり易い)。システムは、OCT画像のうちRNFL分析を行う領域(またはその一部)である上側反射領域を特定すると、元のOCT画像(デジタルメモリに格納されており、JPEG、TIFF、GIFF、PDF等任意のファイル形式で、または、OCTデバイスから得られたままの未加工のデータとして、保存されている)へと戻り、上側反射層として分離された領域に対して点分析を実行する。カラー画像を利用して、ソフトウェアは、RNFLの画像全体に分散しているさまざまな色の点を検出、分類およびカウントする前に、離散化処理を実行する。離散化処理は、各色について、取り得るピクセル色値の範囲を決定することを意味する。好ましくは、複数の同様の色/階調から成る大きいスペクトルを含むように、広い範囲の値とする。例えば、明るいピンク色、暗いピンク色、明るい赤色、暗い赤色、えび茶色を含む色の範囲は、「赤」と分類する。検出すべき他の色(例えば、黄色、緑および青色)についても、同様に色の範囲を設定する。このような色範囲のルールを利用して、ソフトウェアはこの後、各点を分析して、4つ(以上)の色のうち1つに分類する。
例えば、検出された点毎に、システムはピクセル値の範囲を決定する。ピクセル値範囲に基づき、点には色が割り当てられる。
本発明の実施形態によると、ソフトウェアは、離散化によってOCT画像を4つの色(黄色、赤色、緑色および青色)に分類する。このため、各色は<R,G,B>の3色から成る組み合わせと見なされる。尚、Rは赤色、Gは緑色、および、Bは青色である。4つのクラス(つまり、代表点)は、黄色=<0,1,1>、赤色=<1,0,0>、緑色=<0,1,0>、および青色=<0,0,1>と定義される。所与の点は、所与のクラスに対する距離に応じて当該クラスに属するものとする。2つの3D点の間の距離は、以下の式で算出される。
Figure 2013531262
本発明の一実施形態によると、グレースケールのOCT画像を点分析に用いる。本実施形態によると、ピクセル値の所定のしきい値を、さまざまな階調に対応付ける。ソフトウェアは、OCT画像を、離散化により4つの階調に分類する。ソフトウェアは、所与の点と対応する(近い)階調との間の距離の同じ原則を利用する。
画像を分析すると、コンピュータは、各分類の点の数を集計する。コンピュータは、それぞれの点の集計結果を、基準データベースの同様の集計結果と比較して、基準からの逸脱度に基づきスコアを導き出す。本発明の実施形態によると、この比較は、分析された特定のOCTの一部について実行され、OCT画像全体について実行される。
これに加えて、または、これに代えて、コンピュータは、各分類の点の集計結果に基づき、点同士の比率を求めるべく比率分析を実行する。最も好ましいのは、コンピュータが、ホット点の数に対するコールド点の数を検査することである。
一の好ましい実施形態によると、コンピュータは、ホット点とコールド点との間の関係を説明するべく数値指標を割り当てる。この指標は、0と1との間の数字で表現され、少ない方の点の数を多い方の点の数で除算することで得られる。この数値が1に近くなるほど、緑内障の症状が見られる可能性が高くなる。これは、健常な目では、支持細胞よりも軸索の方が多いので、コールド点よりもホット点の方が多くなるためである。しかし、軸索損失が発生すると、軸索と支持細胞との比率は小さくなり、結果的には、コールド点に対するホット点の数が少なくなる。緑内障を患った目では、この比率は、上述した指標で言うと、「1」という高い値になる場合がある。本発明の一実施形態によると、得られた指標を通常指標と比較して、得られた指標と通常指標との差分の程度に基づいてスコアを生成する。スコアが高い場合には、得られた指標と通常指標との差分の程度が大きく、逆もまた然りである。
別の好ましい実施形態によると、比率分析は、全てのホット点(つまり、赤色の点および黄色の点)の集計結果を合計し、コールド点(つまり、緑色の点および青色の点)の集計結果を合計し、ホット点対コールド点の総比率を求めることによって行う。
実施形態によると、比率分析はセクタ毎に行う。この実施形態によると、OCT画像全体を複数のセクタに分割する。各セクタは、個別に分析する。好ましい実施形態によると、OCT画像全体を12個から16個のセクタに分割する。本発明の一実施形態によると、各セクタをさらに、上側部分および下側部分に分割する。絶対数比較および比率決定等の分析は、セクタ毎に、上側部分および下側部分のそれぞれについて実行する。
別の実施形態によると、上述した比率決定に加えて、または、それを補完する分析として、各分類の点の絶対数を、通常の健常なOCT結果を格納しているデータベースに保存されている対応する点の数と比較する。つまり、コンピュータには、健常な症状を見せる人で見られる各分類の点の「通常」分布の範囲がプログラミングされている。OCT画像の分析の一環として、コンピュータは、各分類の点の数を、この特定の分類の点の「通常」の数と比較する。分析されたOCT画像に含まれる点の絶対数が「通常」の範囲内に無い場合、緑内障を患っている可能性がある。本発明の好ましい実施形態によると、「通常」の範囲、「異常」の範囲、および、その間にある「境界線」の範囲がある。つまり、好ましい実施形態では、通常と異常との間には明確な分断が無い。逆に、絶対数が通常の範囲内であれば、システムは、RNFLが「通常」であると判断する。数字が通常範囲から十分に離れていれば(つまり、軸索細胞に大きな減少が見られれば)、システムは、RNFLが異常であると見なす。数字が極値の間にあれば、システムは、RNFLが「境界線」であると判断する。
好ましい実施形態によると、システムは、上述したように(つまり、比率分析または絶対数比較)比較することで、結果または指標を生成する(例えば、表示スクリーン上に測定値を表示したり、印刷したり等する)。好ましい実施形態によると、コンピュータは、「通常」、「異常」および「境界線」のうち1つを生成する。これに代えて、コンピュータは、点分析に基づき、画像がどの程度「通常」または「異常」に近いかを示す指標を生成する。(尚、「点分析」および「構造分析」といった用語は、本明細書で用いる場合、OCT画像に分散した点を計数して、本明細書で説明したように、点の絶対数を比較するか、または、点の絶対数の比率を算出することで、その数を分析する方法を広く意味する。)
臨床医は、「境界線」という結果が出た患者を「フォロー」して、疾患の進行が検出され得るか否かを判断することを希望する場合がある。重要な点は、本発明によれば、臨床医は、次回の測定値と比較する絶対数を持つという利点を得ることにある。さらに、臨床医は、全体的に測定値を比較することに加えて、セグメント単位での比較を何度も繰り返すことができる。これによって、眼疾患について局所的で非常に高精度の検査を実行することが可能となる。
図5は、その一部において、局所欠陥を示すOCT画像である。局所欠陥の領域は、拡大して円32で囲っている。図示しているように、円で囲った領域内では、残りの反射層(円32の外側の領域)に比べて、ドットの密度が少なくなっている。同図から分かるように、RNFLのうち罹患していない領域では、点34が密につまっており互いに接している。しかし、円32で囲った領域には、点がない空間36がある(または、OCT画像の「黒点」)。このような欠陥は、RNFLのセクタ化を利用する本発明に係るシステムによって、つまり、上述した方法のうち1つで、特定され得る。例えば、セクタ毎の分析を利用して、円で囲った領域は、点が無い領域(または、「黒点」)があることから、このセクタ内の点の絶対数は少なくなっているので、異常であると判断される。このように点の数が少なくなっていることは、基準セクタと比較して算出されるとしてよく、または、同じOCT画像内の隣接するセクタと比較して算出されるとしてよい。
これに加えて、または、これに代えて、図5に示す欠陥は、比率分析を行った結果検出されるとしてもよい。つまり、点が無い領域36が一の特定の種類の細胞(例えば、ホット点)が無くなることで発生している場合、円で囲ったセクタ内のホット点とコールド点との比率は、上述したように、異常を示す。
図5に示す欠陥は画像の全体分析に基づいても検出され得ると理解されたい。つまり、通常のホット点対コールド点の比率が全体的に逆転していることによって、または、上述したように、点の絶対数が減っていることによって、欠陥が発生していると判断されるとしてよい。
本発明の別の実施形態によると、上述した比率分析および/または絶対数分析は、視神経(篩板)の画像に対して実行される。視神経のOCT画像は通常、組み合わせると視神経を表す(1枚のピザが複数のスライスから成る円であるのと略同様)複数(例えば、6枚)のスライスを含む。
緑内障を患った目では、視神経による光の反射が、健常な目に比べると、弱い。このように反射率が低下すると、OCT画像でホット点が減る。反射率が低下すると、ホット点はさらに少なくなっていく。上述したのと同様の方法で、システムは、存在するホット点の数を数値化して、緑内障が存在する可能性を判断する。
上述したように、ホット点の絶対数および相対数を求めて、所定の「通常」の範囲と比較する。この比較は、スライス毎に実行され、視神経の画像全体に対しても実行される。
また、比率は、上述した方法に従って求められる。つまり、コールド点の絶対数に対するホット点の絶対数に基づき、数値指標またはスコアを求める。スコアは、ホット点とコールド点との差異の程度を表している。スコアは、緑内障が発生している可能性を判断するために利用される。
尚、本明細書で説明するシステムおよび方法は、網膜画像で利用することに限定されるものではなく、本発明に係るシステムは人体のさまざまな箇所の組織を分析するために用いられ得ることに留意されたい。本開示で説明する原理は、広義には、軸索またはその他の種類の細胞に対応する識別子を数値化することによる組織画像の分析に関する。
本明細書で説明したOCT画像の分析は、さまざまな演算を実行する特別なソフトウェアを持つコンピュータによって実行される。演算を制御して、ルーチンを実行し、データを格納するために用いられるコンピュータは、少なくとも1以上のプロセッサおよびメモリストレージデバイスを備えるとしてよい。コンピュータはさらに、外部と情報をやり取りするべく、多くの入出力を受信するとしてよい。
コンピュータは、オペレーティングシステムおよびソフトウェアアプリケーション、本明細書で説明したルーチンおよびシステムを実行するコンポーネントおよびプログラムの制御下で動作する。概して、本発明の実施形態を実装するべく実行されるルーチンは、オペレーティングシステムまたは特定のアプリケーション、コンポーネント、プログラム、オブジェクト、モジュールまたは命令列の一部として実装されるか否かに関わらず、本明細書において、「システム」または「ソフトウェア」と呼ぶ。
図6は、本発明の実施形態に係るソフトウェアによって実行されるステップの一部を説明するためのフローチャートである。ステップ40において(「S40」と図示)、関連技術分野で公知の方法でOCT画像を取得する。OCT画像はこの後、ソフトウェアによってグレースケールに変換され(S42)、二値化される(S44)。ソフトウェアはこの後、例えば、上述した方法を用いて中点の位置を特定し(S46)、中点より上方の領域を分析のために分離する(S48)。対象領域を分離して画定すると、システムは、分離された領域をマップとして利用する。つまり、システムは、元のカラーOCT画像に対して、しかし、上述したように、分離した領域に対してのみ、点分析を実行する(S50)。これに代えて、システムは、グレースケールOCT画像に対して点分析を実行する(S52)。
図7は、本発明の別の実施形態における、ソフトウェアによって実行されるステップを示すフローチャートである。S54において、OCT画像を取得する。この後、システムは、CCL処理を実行してブロブを検出する(S56)。S56で検出した一群のブロブから、最大のブロブを選択する(S58)。同様に、S56で検出した一群のブロブから、二番目に大きいブロブをこの後検出する(S60)。最初のブロブは、上側反射層として分離する(S62)。この後、OCT画像のうち最初のブロブで画定される領域に対して、点分析を実行する(S64)。
本明細書で説明するシステムは、人体のさまざまな箇所から採取したさまざまな種類の組織のOCT画像を分析するためにも利用され得る。例えば、動脈壁を可視化するべくOCTを利用することは、関連技術分野で公知である。しかし、本発明に係るシステムによれば、動脈を可視化することに加えて、臨床医は、心臓病に関連するさまざまな種類の細胞および/または血小板を数値化することができるようになる。神経学の分野では、本発明に係るシステムは、OCTと組み合わせて、複数の硬化症または任意のその他の神経−眼疾患に関連付けられている神経細胞の喪失を数値化するために用いることができる。歯科学の分野では、本発明は、OCT画像のエナメル構造またはその他の組織構造を数値化するために用いられ得る。
当業者であれば、本発明は医療分野/生物化学分野に限定されないものと理解されたい。逆に、本発明は、任意のOCT画像を構成している複数の層の内容を数値化するために用いられ得る。例えば、OCTは、絵画に含まれる複数の顔料層を個別に正確に検出することが分かっている。本発明は、特定された層に含まれるさまざまな塗料/顔料を数値化するために利用することができる。
具体的な実施形態を参照しつつ本発明を説明してきたが、さまざまな点でさらに変形および変更が可能であることが当業者には明らかであり、上記の説明は本発明を限定することを意味するものではないと理解されたい。本願はこのような変形例および変更例もすべて含むものとする。

Claims (18)

  1. 異常を検出するべくOCT画像を分析する方法であって、
    上側反射層および下側反射層を少なくとも含み、さらに複数の点を含むOCT画像を取得して、デジタルメモリに前記OCT画像を格納する段階と、
    点分析対象領域として、前記OCT画像のうち前記上側反射層の少なくとも一部を分離する段階と、
    それぞれ異なる色に対応する複数の異なるピクセル値範囲を決定する離散化処理を実行する段階と、
    一の決定したピクセル値範囲に基づいて、前記OCT画像の前記上側反射層の複数の点に一の色を割り当てる段階と、
    一の色が割り当てられた前記複数の点を計数する段階と
    を備える方法。
  2. 前記上側反射層の少なくとも一部を分離する段階は、
    前記OCT画像をグレースケール画像に変換する段階と、
    前記グレースケール画像を二値化する段階と、
    前記グレースケール画像において中点を定義する段階と、
    前記中点より上方の領域を分離する段階と
    を有する請求項1に記載の方法。
  3. 前記上側反射層の少なくとも一部を分離する段階は、
    前記OCT画像に対してブロブ抽出処理を実行して一群のブロブを生成する段階と、
    前記一群のブロブから最大ブロブを選択する段階と、
    前記最大ブロブの領域をマッピングする段階と、
    前記最大ブロブの前記領域のマッピングに基づいて、前記上側反射層の領域において点分析を実行する段階と
    を有する請求項1に記載の方法。
  4. 4つの異なるピクセル値範囲を決定する段階を備え、
    前記4つの異なるピクセル値範囲は、4つの異なる色に対応する請求項1に記載の方法。
  5. 前記4つの異なる色は、赤色、黄色、緑色および青色である請求項4に記載の方法。
  6. 前記複数の点を計数して、前記4つの異なる色の各色の点の絶対数を得る段階をさらに備える請求項5に記載の方法。
  7. 通常のOCT画像について決定された点の絶対数に対して、前記4つの異なる色の各色の点の絶対数を比較する段階をさらに備え
    前記通常のOCT画像について決定された点の絶対数は、データベースに格納されている請求項6に記載の方法。
  8. 前記比較する段階に基づき、異常の可能性を含む結果を生成する段階をさらに備える請求項7に記載の方法。
  9. 前記赤色の点および前記黄色の点を第1のグループとし、前記緑色の点および前記青色の点を第2のグループとするグループ分けを行う段階と、
    前記第1のグループの前記赤色の点および前記黄色の点の数を計数し、前記第2のグループの前記緑色の点および前記青色の点の数を計数する段階と、
    前記第2のグループに含まれる点に対する前記第1のグループに含まれる点の比率を求める段階と
    をさらに備える請求項5に記載の方法。
  10. 前記第2のグループに含まれる点に対する前記第1のグループに含まれる点の前記比率に基づき、指標を生成する段階をさらに備える請求項9に記載の方法。
  11. 前記上側反射層は複数のセクタを含み、前記4つの異なる色の各色の点の絶対数を比較する段階は、前記セクタ毎に実行される請求項7に記載の方法。
  12. 前記複数のセクタのうち第1のセクタに含まれる点の数を、前記複数のセクタのうち第2のセクタに含まれる点の数と比較する請求項11に記載の方法。
  13. 前記上側反射層は複数のセクタを含み、前記比率を求める段階は、前記複数のセクタに対して実行される請求項9に記載の方法。
  14. 異常を検出するべく視神経を分析する方法であって、
    上側反射層および下側反射層を少なくとも含み、複数の点をさらに含む、視神経のOCT画像を取得して、前記OCT画像をデジタルメモリに格納する段階と、
    前記視神経の前記OCT画像を複数のスライスに分割する段階と、
    それぞれ異なる色に対応する複数の異なるピクセル値範囲を決定する離散化処理を実行する段階と、
    一の決定されたピクセル値範囲に基づき、前記OCT画像の前記上側反射層に含まれる複数の点に一の色を割り当てる段階と、
    一の色が割り当てられた前記複数の点を計数する段階と
    を備える方法。
  15. 4つの異なるピクセル値範囲を決定する段階をさらに備え、前記4つの異なるピクセル値範囲は、4つの異なる色に対応する
    請求項14に記載の方法。
  16. 前記4つの異なる色は、赤色、黄色、緑色および青色である請求項15に記載の方法。
  17. 前記赤色の点および前記黄色の点を第1のグループとし、前記緑色の点および前記青色の点を第2のグループとするグループ分けを行う段階と、
    前記第1のグループの前記赤色の点および前記黄色の点の数を計数し、前記第2のグループの前記緑色の点および前記青色の点の数を計数する段階と、
    前記第2のグループに含まれる点に対する前記第1のグループに含まれる点の比率を決定する段階と
    をさらに備える請求項16に記載の方法。
  18. OCT画像を分析する方法であって、
    複数の点を含むOCT画像を取得する段階と、
    前記OCT画像の上側反射層の一部を分離する段階と、
    離散化処理を実行する段階と、
    前記OCT画像において複数の点を検出する段階と、
    前記OCT画像において検出された前記複数の点を分類する段階と、
    前記OCT画像で分類された前記複数の点を計数する段階とを備える方法。
JP2013520870A 2010-07-21 2011-07-21 Oct画像を分析する方法および視神経を分析する方法 Expired - Fee Related JP5884195B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US36648510P 2010-07-21 2010-07-21
US61/366,485 2010-07-21
PCT/US2011/044870 WO2012012646A1 (en) 2010-07-21 2011-07-21 Method and system for analyzing optical coherence tomography (oct) results using color reflectivity discrelization analysis

Publications (3)

Publication Number Publication Date
JP2013531262A true JP2013531262A (ja) 2013-08-01
JP2013531262A5 JP2013531262A5 (ja) 2014-08-28
JP5884195B2 JP5884195B2 (ja) 2016-03-15

Family

ID=45493642

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013520870A Expired - Fee Related JP5884195B2 (ja) 2010-07-21 2011-07-21 Oct画像を分析する方法および視神経を分析する方法

Country Status (5)

Country Link
US (2) US8787645B2 (ja)
EP (1) EP2596456B1 (ja)
JP (1) JP5884195B2 (ja)
CA (1) CA2805256C (ja)
WO (1) WO2012012646A1 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5884195B2 (ja) * 2010-07-21 2016-03-15 ディオプシス、インコーポレーテッド Oct画像を分析する方法および視神経を分析する方法
US8632180B2 (en) * 2011-04-25 2014-01-21 Carl Zeiss Meditec, Inc. Automated detection of uveitis using optical coherence tomography
CA2844433A1 (en) * 2011-08-09 2013-02-14 Optovue, Inc. Motion correction and normalization of features in optical coherence tomography
JP2014083285A (ja) * 2012-10-25 2014-05-12 Canon Inc 画像処理装置及び画像処理方法
CN104434027B (zh) * 2013-09-12 2016-03-02 广东福地新视野光电技术有限公司 通过眼科光学相干断层扫描图像获取视盘面积的方法和装置
EP3052002A4 (en) * 2013-09-30 2017-06-28 Barbara L. Dunning Detection of nutrient deficiencies influencing ocular health
CN103854284B (zh) * 2014-02-26 2016-08-31 苏州大学 基于三维图搜索浆液性色素上皮层脱离的视网膜分割方法
US10149610B2 (en) * 2014-04-25 2018-12-11 Carl Zeiss Meditec, Inc. Methods and systems for automatic detection and classification of ocular inflammation
US9918627B2 (en) 2014-05-06 2018-03-20 David Huang Aqueous cell differentiation in anterior uveitis using optical coherence tomography
US9610010B2 (en) 2014-09-03 2017-04-04 The Cleveland Clinic Foundation Systems and methods for analysis of ocular inflammation
CN116453117B (zh) * 2023-02-14 2024-01-26 珠江水利委员会珠江水利科学研究院 一种图像解译图斑定位与精度评估方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006022045A1 (ja) * 2004-08-26 2006-03-02 National University Corporation Nagoya University 光干渉断層装置
JP2008275529A (ja) * 2007-05-02 2008-11-13 Fujifilm Corp 断層画像処理方法および装置ならびにプログラム
JP2009128260A (ja) * 2007-11-27 2009-06-11 Seiko Instruments Inc 検査装置、検査方法、転がり軸受の製造装置、及び転がり軸受の製造方法
JP2009279031A (ja) * 2008-05-19 2009-12-03 Topcon Corp 眼底観察装置
JP2010094381A (ja) * 2008-10-17 2010-04-30 Canon Inc 画像処理装置、画像処理方法
JP2010151685A (ja) * 2008-12-25 2010-07-08 Univ Of Tsukuba 光画像計測装置における画像情報の解析及び組織判別可視化装置

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5421339A (en) * 1993-05-12 1995-06-06 Board Of Regents, The University Of Texas System Diagnosis of dysplasia using laser induced fluoroescence
CA2236268A1 (en) 1995-11-30 1997-06-05 Chromavision Medical Systems, Inc. Method and apparatus for automated image analysis of biological specimens
GB0005792D0 (en) * 2000-03-11 2000-05-03 Shelton Vision Systems Ltd Automatic inspection method
GR1004180B (el) * 2000-03-28 2003-03-11 ����������� ����� ��������� (����) Μεθοδος και συστημα χαρακτηρισμου και χαρτογραφησης αλλοιωσεων των ιστων
EP1754016B1 (en) 2004-05-29 2016-05-18 The General Hospital Corporation Process, system and software arrangement for a chromatic dispersion compensation using reflective layers in optical coherence tomography (oct) imaging
JP3944205B2 (ja) 2004-10-05 2007-07-11 キヤノン株式会社 眼科撮影装置及び眼科画像処理システム
JP2007079577A (ja) * 2005-09-12 2007-03-29 Samsung Electro Mech Co Ltd 両方向スキャン方式のカラーディスプレイ装置
WO2007061769A2 (en) * 2005-11-18 2007-05-31 Duke University Method and system of coregistrating optical coherence tomography (oct) with other clinical tests
US8160676B2 (en) * 2006-09-08 2012-04-17 Medtronic, Inc. Method for planning a surgical procedure
US8160677B2 (en) * 2006-09-08 2012-04-17 Medtronic, Inc. Method for identification of anatomical landmarks
WO2008137077A1 (en) * 2007-05-04 2008-11-13 Bioptigen, Inc. Methods, systems and computer program products for mixed-density optical coherence tomography (oct) imaging
US20090137988A1 (en) * 2007-11-02 2009-05-28 Lensx Lasers, Inc Methods And Apparatus For Improved Post-Operative Ocular Optical Performance
WO2010009447A2 (en) 2008-07-18 2010-01-21 Doheny Eye Institute Optical coherence tomography - based ophthalmic testing methods, devices and systems
CA3060821C (en) 2008-12-19 2023-09-12 University Of Miami System and method for early detection of diabetic retinopathy using optical coherence tomography
US8744159B2 (en) * 2010-03-05 2014-06-03 Bioptigen, Inc. Methods, systems and computer program products for collapsing volume data to lower dimensional representations thereof using histogram projection
JP5884195B2 (ja) * 2010-07-21 2016-03-15 ディオプシス、インコーポレーテッド Oct画像を分析する方法および視神経を分析する方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006022045A1 (ja) * 2004-08-26 2006-03-02 National University Corporation Nagoya University 光干渉断層装置
JP2008275529A (ja) * 2007-05-02 2008-11-13 Fujifilm Corp 断層画像処理方法および装置ならびにプログラム
JP2009128260A (ja) * 2007-11-27 2009-06-11 Seiko Instruments Inc 検査装置、検査方法、転がり軸受の製造装置、及び転がり軸受の製造方法
JP2009279031A (ja) * 2008-05-19 2009-12-03 Topcon Corp 眼底観察装置
JP2010094381A (ja) * 2008-10-17 2010-04-30 Canon Inc 画像処理装置、画像処理方法
JP2010151685A (ja) * 2008-12-25 2010-07-08 Univ Of Tsukuba 光画像計測装置における画像情報の解析及び組織判別可視化装置

Also Published As

Publication number Publication date
EP2596456B1 (en) 2019-09-25
US9367915B2 (en) 2016-06-14
US20140334704A1 (en) 2014-11-13
CA2805256A1 (en) 2012-01-26
US20120020539A1 (en) 2012-01-26
US8787645B2 (en) 2014-07-22
CA2805256C (en) 2019-03-05
JP5884195B2 (ja) 2016-03-15
EP2596456A1 (en) 2013-05-29
EP2596456A4 (en) 2017-11-08
WO2012012646A1 (en) 2012-01-26

Similar Documents

Publication Publication Date Title
JP5884195B2 (ja) Oct画像を分析する方法および視神経を分析する方法
US11972568B2 (en) Automated assessment of glaucoma loss from optical coherence tomography
US12035971B2 (en) Automated determination of arteriovenous ratio in images of blood vessels
Amin et al. A method for the detection and classification of diabetic retinopathy using structural predictors of bright lesions
US11779205B2 (en) Image classification method, device, and program
KR20190087272A (ko) 안저영상을 이용한 녹내장 진단 방법 및 이를 위한 장치
CN103348359A (zh) 使用光学相干断层照相的3d视网膜破裂检测
JP7467205B2 (ja) バイオマーカを光学的に検出するための方法
Hassan et al. Structure tensor graph searches based fully automated grading and 3D profiling of maculopathy from retinal OCT images
Soltanian-Zadeh et al. Weakly supervised individual ganglion cell segmentation from adaptive optics OCT images for glaucomatous damage assessment
Zang et al. A diabetic retinopathy classification framework based on deep-learning analysis of OCT angiography
Gao et al. A deep learning network for classifying arteries and veins in montaged widefield OCT angiograms
Odstrčilík et al. Retinal nerve fiber layer analysis via markov random fields texture modelling
Bahr et al. Deep Learning and Machine Learning Algorithms for Retinal Image Analysis in Neurodegenerative Disease: Systematic Review of Datasets and Models
US20220061920A1 (en) Systems and methods for measuring the apposition and coverage status of coronary stents
Vázquez et al. Evaluation of SIRIUS retinal vessel width measurement in REVIEW dataset
EP3655921B1 (en) Method and computer program for estimating shape parameters of the fovea by optical coherence tomography
Umamageswari et al. Identifying Diabetics Retinopathy using Deep Learning based Classification
Mani et al. An automated hybrid decoupled convolutional network for laceration segmentation and grading of retinal diseases using optical coherence tomography (OCT) images
Sakthi Karthi Durai et al. An Effective Approach Based on Improved Convolutionary Neural Network Classifier for the Detection of Diabetic Retinopathy
Yalvac Artificial Intelligence and Glaucoma
US20230380760A1 (en) Systems and methods for classifying ophthalmic disease severity
Strouthidis et al. Optic Nerve: Heidelberg Retinal Tomography
Marrapu Detection of Glaucoma Using Deep Learning Techniques: Literature Survey
Himbitskaya et al. Generation of Description for Eye Fundus Desease

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140704

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140704

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150224

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150227

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150521

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150622

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150724

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

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20160113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160118

R150 Certificate of patent or registration of utility model

Ref document number: 5884195

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

LAPS Cancellation because of no payment of annual fees