JP6866167B2 - 情報処理装置、情報処理方法及びプログラム - Google Patents
情報処理装置、情報処理方法及びプログラム Download PDFInfo
- Publication number
- JP6866167B2 JP6866167B2 JP2017006977A JP2017006977A JP6866167B2 JP 6866167 B2 JP6866167 B2 JP 6866167B2 JP 2017006977 A JP2017006977 A JP 2017006977A JP 2017006977 A JP2017006977 A JP 2017006977A JP 6866167 B2 JP6866167 B2 JP 6866167B2
- Authority
- JP
- Japan
- Prior art keywords
- motion contrast
- image
- data
- tomographic
- eye
- 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
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/0016—Operational features thereof
- A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Multimedia (AREA)
- Eye Examination Apparatus (AREA)
Description
被検眼の同一位置を同一方向に複数回走査するように制御された測定光に基づいて得られた前記被検眼の断層を示す複数の断層データを含むデータセットであって、互いに同一の断層データを含まない複数のデータセット間からモーションコントラスト画像を生成することなく、前記被検眼の深さ方向の断面を示す少なくとも1のモーションコントラスト画像であって、前記複数のデータセットのそれぞれから少なくとも1のモーションコントラスト画像を生成することにより、前記複数のデータセットに対応する複数のモーションコントラスト画像を生成する第1生成手段と、
前記複数のモーションコントラスト画像に対応する複数の3次元のモーションコントラスト画像を用いて、1の3次元のモーションコントラスト画像を生成する第2生成手段と、
前記データセットに含まれる複数の断層データのうち少なくとも1つの断層データを用いて得た層境界に関する情報を用いて、前記被検眼の深さ方向における前記1の3次元のモーションコントラスト画像の一部の範囲を特定し、前記特定された一部の範囲に対応するモーションコントラストのEn−face画像を生成する第3生成手段と、を備える。
被検眼の同一位置を同一方向に複数回走査するように制御された測定光に基づいて得られた前記被検眼の断層を示す複数の断層データを含むデータセットであって、互いに同一の断層データを含まない複数のデータセット間からモーションコントラスト画像を生成することなく、前記被検眼の深さ方向の断面を示す少なくとも1のモーションコントラスト画像であって、前記複数のデータセットのそれぞれから少なくとも1のモーションコントラスト画像を生成することにより、前記複数のデータセットに対応する複数のモーションコントラスト画像を生成する第1生成工程と、
前記複数のモーションコントラスト画像に対応する複数の3次元のモーションコントラスト画像を用いて、1の3次元のモーションコントラスト画像を生成する第2生成工程と、
前記データセットに含まれる複数の断層データのうち少なくとも1つの断層データを用いて得た層境界に関する情報を用いて、前記被検眼の深さ方向における前記1の3次元のモーションコントラスト画像の一部の範囲を特定し、前記特定された一部の範囲に対応するモーションコントラストのEn−face画像を生成する第3生成工程と、を含む。
以下、図面を参照しながら、第1実施形態に係る情報処理装置を備える情報処理システムについて説明する。
本実施形態のOCT装置100における測定光学系及び分光器の構成の一例について図2(b)を用いて説明する。
本実施形態の情報処理装置101の構成の一例について図1を用いて説明する。
次に、図3を参照して本実施形態の情報処理装置101の処理手順の一例を示す。図3は、本実施形態における本システム全体の動作処理の流れを示すフローチャートである。
検者は入力部103を操作することにより、情報処理装置101の撮影制御部101−03に対してモーションコントラスト画像の撮影条件の設定および撮影開始の指示を行う。撮影制御部101−03は、入力部103を介して検者の指示を受信する。
撮影制御部101−03はS301で指示された撮影条件に基づいてクラスタ走査を行うようOCT装置100を制御する。すなわち、撮影制御部101−03は、眼底の同一位置におけるBスキャンをm回繰り返すように、OCT装置100を制御する。OCT装置100は、当該制御に基づいて眼底の同一位置におけるBスキャンをm回繰り返す。すなわち、OCT装置100はクラスタ走査を実行する。
移動量計算部101−31は、SLO眼底画像(動画像)の時間的に隣接するフレーム同士を位置合わせすることにより、被検眼200の眼底における移動量を計算する。また、移動量計算部101−31は、計算した移動量を記憶部101−02へ格納する。
画像取得部101−01は被検眼200の眼底の2次元画像を取得する。具体的には、SLO走査手段204を構成するX走査スキャナとY走査スキャナによって眼底上を走査することで得られた戻り光に基づいて、画像取得部101−01はSLO眼底画像を取得する。なお、画像取得部101−01が取得するSLO眼底画像は例えば動画である。すなわち、画像取得部101−01はSLO眼底画像を順次取得する。
移動量計算部101−31は、記憶部101−02に保存された以前のSLO眼底画像と現在のSLO眼底画像の2枚の眼底画像を用いて移動量を算出する。具体的には、移動量計算部101−31は、画像上の着目領域における2次元(X,Y方向)の変位量を検出することによって、眼底の移動量を算出する。さらに、X,Y方向の変位だけでなく回転方向の角度の変化量を算出し、移動量に含めても良い。
撮影制御部101−03は算出された眼底の移動量を記憶部101−02へ格納する。これは後述する走査位置の補正や再走査の際に情報処理装置101によって使用されるためである。
撮影制御部101−03は取得した眼底画像を次の移動量算出で用いるために記憶部101−02へ格納する。なお、ステップS303−04はステップS303−01と同時に実施されることとしても良い。
図3の説明に戻り、ステップS304において、判定部101−32は測定光の走査中に眼底の移動が発生したかどうかを判定する。本実施形態においては、例えば、判定部101−32は、ステップS303において計算・保存した位置ずれ量を記憶部101−02から読み出せた場合は、眼底の移動が発生したものと判定する。また、判定部101−32は、記憶部101−02に位置ずれ量が保存されていない場合は、眼底の移動が発生していないと判定する。眼底の移動があった場合には、処理はステップS305へ、眼底の移動が生じなかった場合には処理はステップS308へ進む。なお、判定部101−32眼底の移動が所定閾値以下の場合には眼底の移動がないと判定することとしてもよい。
判定部101−32は、ステップS303で計算された眼底の移動量が所定の閾値以上か否かを判定する。本実施形態においては、記憶部101−02に格納された移動量が閾値以上であるかどうかを判定する。閾値は再走査を行うかどうかを判断するための閾値である。
移動量が閾値以上であった場合、決定部101−33はその部位を後で再走査するために、再走査すべきBスキャン番号を決定し、記憶部101−02へ格納する。決定部101−33は、再走査すべきBスキャン番号として、閾値以上の移動量が検出されたSLO画像に対応するBスキャン画像を含むクラスタの識別情報を決定し記憶部101−02に格納する。例えば、決定部101−33は、閾値以上の移動量が検出されたSLO眼底画像が取得されたタイミングと略同一タイミングで取得された断層画像を含むクラスタを決定する。なお、Bスキャン番号とは眼底のどの位置におけるどのクラスタのBスキャンであるか特定可能な識別情報であってもよい。決定部101−33は、例えばステップS302にてクラスタスキャンが実行された際のスキャナの制御情報からBスキャン番号を決定することが可能である。
指示部101−34は、眼底の同一部位の走査を継続ため、ステップS303で計算した眼底の移動量をそれまでの走査位置に対してオフセットした位置へXYスキャナ216を移動させるようOCT装置100に指示する。なお、推定する移動量に回転方向の角度も含まれる場合には、回転方向の補正も行う。
撮影制御部101−03は再走査すべきBスキャン番号が記憶部101−02に格納されているかどうかを確認する。もし該Bスキャン番号が格納されていればS309へ処理を進め、格納されていなければS310へ処理を進める。なお再走査が何度も繰り返されると撮影時間が長くなり被検者の負担が増えるため、再走査の回数に上限を設けてもよい。その場合、再走査の回数が上限に達した場合は再走査すべきBスキャン番号が記憶部101−02に格納されていても、撮影制御部101−03はステップS310へ遷移する。
撮影制御部101−03は再走査すべきBスキャン番号を記憶部101−02に設定し、ステップS302に戻る。ステップS302では、指示部101−34が記憶部101−02からBスキャン番号を読出し、当該Bスキャン番号に対応するクラスタ走査を行うようOCT装置100に指示を行う。この指示に基づいてOCT装置100は再度クラスタ走査を実行する。すなわち、再走査すべきと判定されたBスキャンを含むクラスタが再度走査されることとなる。この再走査により得られた信号データを画像取得部101−01が取得する。すなわち、画像取得部101−01は、決定手段により決定されたデータセットに含まれる数の断層データのみを再走査により再取得する再取得手段の一例に相当する。
撮影制御部101−03はXYスキャナ216を次の走査位置へ移動させる。例えば、撮影制御部101−03は異なる副走査方向の位置を測定光で走査するようにXYスキャナ216を制御する。
撮影制御部101−03は眼底の同一位置におけるクラスタ走査の実行回数が所定の回数rに達したか否かを判定する。該クラスタ走査の実行回数がrに達していればS312へ、未達であれば被検者が瞬目可能な所定の時間の経過後にS302へ処理を進める。なお、S302へ処理を進める際の待ち時間の間に表示部104に被検者が瞬目可能な時間であることを示す(所定の色、形状を持つ)マークもしくは文字を表示させてもよい。また瞬目可能な時間を設けることは本発明に必須の処理ではなく、直ちにS302へ処理をすすめてもよい。また、本実施形態ではステップS310の後にステップS311を実行することで、副走査方向の各位置でクラスタ走査を1回行った後、再び副走査方向の各位置でクラスタ走査を行うこととした。しかしながら、ステップS310をステップS311の後に実行することで、副走査方向の1つの位置でクラスタ走査をr回実行した後に、次の副走査方向の所定位置でクラスタ走査をr回実行することとしてもよい。
撮影制御部101−03は全てのBスキャンが終了したかどうかを判定する。未終了のBスキャンがある場合、情報処理装置101の制御はステップS302に戻る。
断層画像生成部101−11は、画像取得部101−01が取得した断層画像の信号データを再構成して断層画像を生成する。なお、断層画像の生成は、ステップS302以降のタイミングであればステップS312以前に実行することとしてもよい。
画像処理部101−04は、副走査方向(Y軸方向)の走査位置yiにおけるk番目のクラスタ(k=1,2,...,r)での信号処理を開始する。
画像処理部101−04は、Y軸方向の走査位置yi(i=1,2,...,n−1,n)における信号処理を開始する。
画像処理部101−04は、Y軸方向の走査位置yiにおけるk番目のクラスタに属するj番目(j=1,2,...,m)のBスキャン画像に対する信号処理を開始する。
断層画像生成部101−11は、断層画像を生成する。具体的には、断層画像生成部101−11は、干渉信号データから固定ノイズを除去する。次に、断層画像生成部101−11は、スペクトラルシェーピング、分散補償を行い、この信号データに離散フーリエ変換を施して得られる複素信号の絶対値を計算することで、深さに対する強度(Intensity)データを得る。断層画像生成部101−11は、例えば、該強度データから任意の領域の切り出し処理を行って断層画像を生成する。なお、ステップS313−04の処理はステップS302以降のタイミングであればステップS312以前に実行することとしてもよい。
画像処理部101−04は、Y軸方向の位置yiにおけるk番目のクラスタでの信号処理の繰り返し回数が所定数(m)に達したかを判定する。すなわち、位置yiでの画像再構成がm回繰り返されたかを判断する。もし位置yiでの画像再構成の回数が所定数(m)に満たない場合はステップS313−03に戻り、同一Y位置での画像再構成を繰り返す。位置yiでの画像再構成の回数が所定数(m)に達した場合は、ステップS313−06へ処理を進める。
位置合わせ部101−41はある位置yiにおけるk番目のクラスタに属するm枚のBスキャンの位置合わせを行う。具体的には、まず位置合わせ部101−41はm枚のうち任意の1枚を基準画像として選択し、他の(m−1)枚のBスキャン画像に対して位置合わせを行う。位置合わせ方法としては、例えば、画像の類似度を表す評価関数を事前に定義しておき、この評価関数の値が最も良くなるように位置合わせ部101−41は画像を変形する。本実施形態では評価関数として相関係数を用いる。これに限らず、任意の公知な類似度もしくは相違度の指標を評価関数として用いてよい。また、本実施形態では画像の変形処理としてアフィン変換を用いるが、これに限らず、任意の公知な位置合わせ手法を用いて位置合わせしてよい。
画像特徴取得部101−44は、各クラスタにおける断層画像の少なくとも1枚を用いて眼底の層境界データを取得する。本実施形態においては、m枚の断層画像を加算平均した重ね合わせ画像に対して可変形状モデルに基づく層境界データの検出処理を行う。本発明はこれに限定されるものではなく、画像特徴取得部101−44は、例えば各クラスタから代表画像を選択して該代表画像から層境界を検出してもよい。代表画像とは例えば、所定フレーム目に撮影した断層画像あるいは最もコントラストの高い断層画像である。また層境界データの検出方法は可変形状モデルに限定されるものではなく、任意の公知のセグメンテーション手法を用いてよい。
第一生成部101−42はY軸方向の位置yiにおけるk番目のクラスタに属するm枚のBスキャンを用いて1以上のモーションコントラスト値を算出する。第一生成部101−42はモーションコントラスト値を算出することでモーションコントラスト画像を生成する。本実施形態では、m枚の断層画像間((m−1)組)において各画素位置で無相関値(Decorrelation)を計算し、該無相関値をモーションコントラスト値として(m−1)枚のモーションコントラスト画像を生成する。さらに、第一生成部101−42は、該(m−1)枚のモーションコントラスト画像を平均化することにより、1クラスタあたり1のモーションコントラスト画像を生成する。なお、本ステップにおいてモーションコントラスト画像を平均化せず後述するステップS313−11で平均化を行うこととしてもよい。すなわち、本ステップにより1クラスタあたり少なくとも1枚のモーションコントラスト画像が得られる。
画像処理部101−04は、Y軸方向の位置yiにおけるクラスタ内のモーションコントラストデータの平均化処理が位置ynに至るまで実行されたかを判定する。クラスタ内のモーションコントラストデータの平均化処理完了したy位置がynに至っていない場合はS313−02に戻り、ynに到達した場合はS313−10へ処理を進める。
画像処理部101−04は、Y軸方向の位置yiでのk番目のクラスタにおける信号処理がr回繰り返されたかを判定する。所定数(r)に満たない場合はS313−01に戻り、所定数(r)に到達した場合はS313−11へ処理を進める。なお、画像処理部101−04は副走査方向の各位置で一度モーションコントラスト画像を生成した後、再び副走査方向の各位置でモーションコントラスト画像を生成することとしたが、本発明はこれに限定されるものではない。例えば、画像処理部101−04は副走査方向の各位置で全てのクラスタに対してモーションコントラスト画像を生成した後に、次の副走査方向における位置でモーションコントラスト画像を生成することとしてもよい。
第二生成部101−43は、各y位置において複数のクラスタから得られた複数のモーションコントラスト画像を合成する。例えば、3つのクラスタが存在する場合には、第二生成部101−43は、それぞれのクラスタから得られた3枚のモーションコントラスト画像を合成する。すなわち、第二生成部101−43は第1生成手段により生成された複数のモーションコントラスト画像を合成することで前記1のモーションコントラスト画像を生成する。
投影部101−45は、S313−07で取得した層境界データに基づいて深度方向の投影範囲を特定した上で、S313−11において生成されたモーションコントラスト値に関する3次元データを深度方向に投影することにより2次元のモーションコントラスト画像(En−face画像)を生成する。本実施形態では投影方法として最大値投影を用いるが、これに限らず任意の公知の投影方法を用いてよい。なお、En−face画像が不要であれば本ステップは省略することとしてもよい。
画像処理部101−04は、S302〜S312で取得した断層画像またはモーションコントラスト画像を検査日時、披検眼を同定する情報と関連付けて記憶部102へ保存する。
図3に戻り、ステップS314において、表示制御部101−05は、ステップS313で作成した断層画像およびモーションコントラスト画像の少なくとも一方を表示部104に表示させる。なお、表示されるモーションコントラスト画像は、En−face画像であってもよいし、3次元画像であってもよい。
上記の実施形態においては複数(r)回のクラスタ走査を眼底の同一位置で行う場合について説明したが、本発明はこれに限定されるものではない。
上記の実施形態では、再走査が行われた場合、再走査が必要と判定された詩信号データの再構成処理は行わず再走査して得られた信号データのみ再構成処理を行う例について説明したが、本発明はこれに限定されない。
第1の実施形態においては、副走査方向の各位置において複数のクラスタ走査を実行する例について説明したが本実施形態に係る情報処理装置は、同一位置で実行するクラスタ走査の回数を眼部の領域ごとに変更する例について説明する。必要な領域のみ他の領域よりも多くクラスタ走査を実行することで、不要なクラスタ走査を削減しつつ必要な領域に関して高コントラストなモーションコントラスト画像を取得できるようになる。また、不要なクラスタ走査を行わないため、眼が安定しなくなる可能性を低減することができる。従って、モーションコントラスト画像を安定して取得することが可能となる。
画像特徴取得部101−44は、画像取得部101−11が取得した眼部の画像から眼部の画像特徴を検出し、眼底の領域ごとに撮影条件を変更(撮影パラメータを重み付け)する。なお、具体的な画像特徴の検出処理および撮影パラメータの重み付けについては図8(b)に示すステップS801−01〜S801−03を用いて説明する。
画像取得部101−01は、OCT装置100から被検眼の画像特徴を取得するために必要な眼部の画像を取得する。本実施形態では眼部の画像としてSLO眼底画像を取得する。なおSLO眼底画像に限らず、画像取得部101−01は例えば眼部の断層画像を取得しても良い。
画像特徴取得部101−44は、ステップS801−01で取得した画像から眼部の画像特徴を取得する。本実施形態では、画像特徴取得部101−44はステップS801−01で取得したSLO眼底画像から眼部の画像特徴として例えば中心窩、網膜血管アーケードおよび視神経乳頭領域の少なくとも1つを取得する。また、画像特徴取得部101−44はSLO眼底画像から網膜動静脈交叉を取得することとしてもよい。眼部の画像特徴の取得方法としては任意の公知の画像処理手法を用いてよく、例えば網膜血管アーケードはSLO眼底画像に対して任意の公知の線強調フィルタを用いることで取得できる。
決定部101−33は、ステップS801−02で取得された眼部の画像特徴に基づいて追加走査方法(追加走査の対象領域及び走査回数)を決定する。
病変の好発部位としては、中心窩無血管領域(FAZ)の境界、中心窩下の網膜外層及び脈絡膜、網膜動静脈交叉部、視神経乳頭部などが挙げられる。また、病変の一例としては毛細血管閉塞または毛細血管瘤等が挙げられる。
眼底形状が所定の条件を満たす領域としては、例えば網膜や脈絡膜を構成する層の厚みが所定の範囲外にある領域が挙げられる。例えば、眼底の所定の層の厚みが正常値未満となっている部位を眼底形状が所定の条件を満たす領域としてもよい。
例えば、1回目のクラスタ走査により得られた断層像から生成したモーションコントラストデータに対して計算した指標値(コントラストや血管密度等)が許容値の範囲外であるような領域である。
なお、ここで示した追加走査に係るm及びrはあくまでも例であり、他の値としてもよい。例えばクラスタ走査間で走査回数(m)の値が異なっていても良い。
撮影制御部101−03はOCT装置100にS801で指示された撮影条件に基づいてクラスタ走査、すなわち眼底の同一位置におけるBスキャンをm回繰り返し実行させる。本実施形態では本ステップでの1回目クラスタ走査の対象領域として図4(c)に示すような視神経乳頭部と黄斑部を含む10mm×10mmの領域が設定されている。そのため、OCT装置100は、10mm×10mmの領域においてm=3のクラスタ走査を行う。
撮影制御部101−03は略同一位置におけるクラスタ走査の実行回数が所定の回数(r)に達したか否かを判定する。本実施形態では、中心窩を含む3mm×3mmの領域に対してのみr=2とする。該クラスタ走査の実行回数が所定の回数(r)に達していればS812へ、未達であれば追加走査を行うためにS802へ処理を進める。なお、10mm×10mmの領域を一度クラスタ走査した後に3mm×3mmの領域をOCT装置100にクラスタ走査させることとしてもよいが、本発明はこれに限定されない。例えば、10mm×10mmの領域をクラスタ走査している際に副走査方向が3mm×3mmの領域に達した時に3mm×3mmの領域に対して2回目のクラスタ走査をOCT装置100に行わせることとしてもよい。
ステップS813の処理はステップS313の処理と基本的には同様であるが、第一実施形態と異なる処理について図7(b)を参照しながら説明する。
第二生成部101−43は、第一生成部101−42が生成したモーションコントラストデータを合成する。具体的には、第一生成部101−42が生成したモーションコントラストデータ同士を位置合わせした上で、各画素位置においてモーションコントラスト値の平均値を持つ画像を生成する。なお、ステップS313−11と同様に平均値以外に中央値などを用いることとしてもよい。
以上、実施形態例を詳述したが、開示の技術は例えば、システム、装置、方法、プログラム若しくは記録媒体(記憶媒体)等としての実施態様をとることが可能である。具体的には、複数の機器(例えば、ホストコンピュータ、インタフェース機器、撮像装置、webアプリケーション等)から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。
101−42 第一生成部
101−43 第二生成部
Claims (20)
- 被検眼の同一位置を同一方向に複数回走査するように制御された測定光に基づいて得られた前記被検眼の断層を示す複数の断層データを含むデータセットであって、互いに同一の断層データを含まない複数のデータセット間からモーションコントラスト画像を生成することなく、前記被検眼の深さ方向の断面を示す少なくとも1のモーションコントラスト画像であって、前記複数のデータセットのそれぞれから少なくとも1のモーションコントラスト画像を生成することにより、前記複数のデータセットに対応する複数のモーションコントラスト画像を生成する第1生成手段と、
前記複数のモーションコントラスト画像に対応する複数の3次元のモーションコントラスト画像を用いて、1の3次元のモーションコントラスト画像を生成する第2生成手段と、
前記データセットに含まれる複数の断層データのうち少なくとも1つの断層データを用いて得た層境界に関する情報を用いて、前記被検眼の深さ方向における前記1の3次元のモーションコントラスト画像の一部の範囲を特定し、前記特定された一部の範囲に対応するモーションコントラストのEn−face画像を生成する第3生成手段と、
を備えることを特徴とする情報処理装置。 - 再走査が必要な断層データを判定する判定手段と、
前記複数のデータセットのうち、前記判定手段により再走査が必要であると判定された断層データが含まれるデータセットを決定する決定手段と、
前記決定手段により決定されたデータセットに含まれる数の断層データのみを再走査により再取得する再取得手段と、
を更に備えることを特徴とする請求項1に記載の情報処理装置。 - 前記断層データの取得の成否を判定する判定手段と、
前記判定手段により取得失敗と判定された断層像データが含まれるデータセットを決定する決定手段と、を更に備え、
前記第2生成手段は、前記複数のデータセットのうち、前記決定手段により決定されたデータセット以外のデータセットから得られたモーションコントラスト画像に基づいて前記被検眼の所定位置における断面を示す1のモーションコントラスト画像を生成することを特徴とする請求項1に記載の情報処理装置。 - 前記複数の3次元のモーションコントラスト画像は、前記同一方向に交差する方向における各位置で1のデータセットを1のクラスタとして取得して得た複数のデータセットと、前記同一方向に交差する方向における各位置で1のデータセットを1のクラスタとして再び取得して得た複数のデータセットとを用いて生成されることを特徴とする請求項1乃至3のいずれか1項に記載の情報処理装置。
- 前記複数のデータセットのうち異なるデータセットに含まれる断層データは、前記被検眼の断層の深さ方向に直交するXY平面におけるX方向またはY方向の少なくともいずれか一方の方向において異なる位置を走査するように制御された測定光に基づいて得られた断層データであることを特徴とする請求項1乃至3のいずれか1項に記載の情報処理装置。
- 前記複数のデータセットそれぞれに含まれる断層データの数は、データセット間で等しいことを特徴とする請求項1乃至5の何れか1項に記載の情報処理装置。
- 前記複数のデータセットそれぞれに含まれる断層データの数は2または3であることを特徴とする請求項1乃至6のいずれか1項に記載の情報処理装置。
- 前記第1生成手段は、前記判定手段により再走査が必要であると判定された断層データを含まないデータセットからモーションコントラスト画像を生成し、且つ、前記再走査により取得された断層データからなるデータセットからモーションコントラスト画像を生成し、
前記第2生成手段は、前記複数のモーションコントラスト画像を合成することで、前記1の3次元のモーションコントラスト画像を生成することを特徴とする請求項2に記載の情報処理装置。 - 前記第1生成手段は、再走査が必要であると判定された断層データが含まれるデータセットをモーションコントラスト画像の生成に用いないことを特徴とする請求項8に記載の情報処理装置。
- 前記第2生成手段は、前記複数のモーションコントラスト画像を合成することで、前記1の3次元のモーションコントラスト画像を生成することを特徴とする請求項1乃至9のいずれか1項に記載の情報処理装置。
- 前記第2生成手段は、前記複数のモーションコントラスト画像の間のモーションコントラストに関する統計値を用いて、前記1の3次元のモーションコントラスト画像を生成することを特徴とする請求項1乃至10のいずれか1項に記載の情報処理装置。
- 操作者からの指示により設定された撮像条件を用いて、前記測定光が前記被検眼の同一位置を同一方向に複数回走査することで1回のクラスタ走査とするクラスタ走査を行うように、前記測定光を制御するための撮影制御手段を更に備えることを特徴とする請求項1乃至11のいずれか1項に記載の情報処理装置。
- 異なる時間に前記被検眼を撮影して得た複数の眼底画像を用いて移動量を算出する移動量算出手段を更に備え、
前記撮影制御手段は、前記算出された移動量を用いて、1のクラスタ走査と次のクラスタ走査との間の走査位置を補正するように、前記測定光を制御することを特徴とする請求項12に記載の情報処理装置。 - 前記データセットに含まれる複数の断層データを用いて得た重ね合わせ画像を用いて、前記層境界に関する情報を取得する取得手段を更に備えることを特徴とする請求項1乃至13のいずれか1項に記載の情報処理装置。
- 前記データセットに含まれる複数の断層データのいずれかに対応する断層画像であって、代表画像として選択された断層画像を用いて、前記層境界に関する情報を取得する取得手段を更に備えることを特徴とする請求項1乃至13のいずれか1項に記載の情報処理装置。
- 前記3次元のモーションコントラスト画像に対応する前記被検眼の第1の領域の一部の領域を、追加走査する第2の領域として決定する決定手段を更に備えることを特徴とする請求項1乃至15のいずれか1項に記載の情報処理装置。
- 前記被検眼において前記測定光を走査する走査手段と、前記測定光を照射した前記被検眼からの戻り光と参照光とを合波して得た干渉光を検出する検出手段と、を備えるOCT装置と、
請求項1乃至16のいずれか1項に記載の情報処理装置と、
を備えることを特徴とするシステム。 - 前記OCT装置は、異なる時間に前記被検眼を撮影して複数の眼底画像を取得するためのSLO光学系の光路と固視灯用の光路とに含まれる一部共通の光路に設けられ、前記一部共通の光路に沿って駆動可能な焦点合わせ用のレンズと、前記一部共通の光路に設けられるSLO走査手段と、を更に備えることを特徴とする請求項17に記載のシステム。
- 被検眼の同一位置を同一方向に複数回走査するように制御された測定光に基づいて得られた前記被検眼の断層を示す複数の断層データを含むデータセットであって、互いに同一の断層データを含まない複数のデータセット間からモーションコントラスト画像を生成することなく、前記被検眼の深さ方向の断面を示す少なくとも1のモーションコントラスト画像であって、前記複数のデータセットのそれぞれから少なくとも1のモーションコントラスト画像を生成することにより、前記複数のデータセットに対応する複数のモーションコントラスト画像を生成する第1生成工程と、
前記複数のモーションコントラスト画像に対応する複数の3次元のモーションコントラスト画像を用いて、1の3次元のモーションコントラスト画像を生成する第2生成工程と、
前記データセットに含まれる複数の断層データのうち少なくとも1つの断層データを用いて得た層境界に関する情報を用いて、前記被検眼の深さ方向における前記1の3次元のモーションコントラスト画像の一部の範囲を特定し、前記特定された一部の範囲に対応するモーションコントラストのEn−face画像を生成する第3生成工程と、
を含むことを特徴とする情報処理方法。 - 請求項19に記載の情報処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017006977A JP6866167B2 (ja) | 2017-01-18 | 2017-01-18 | 情報処理装置、情報処理方法及びプログラム |
US15/872,647 US10980411B2 (en) | 2017-01-18 | 2018-01-16 | Information processing apparatus, information processing method, and storage medium |
EP18152041.2A EP3354191B1 (en) | 2017-01-18 | 2018-01-17 | Information processing apparatus, information processing method, and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017006977A JP6866167B2 (ja) | 2017-01-18 | 2017-01-18 | 情報処理装置、情報処理方法及びプログラム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2018114121A JP2018114121A (ja) | 2018-07-26 |
JP2018114121A5 JP2018114121A5 (ja) | 2020-07-16 |
JP6866167B2 true JP6866167B2 (ja) | 2021-04-28 |
Family
ID=61002882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017006977A Active JP6866167B2 (ja) | 2017-01-18 | 2017-01-18 | 情報処理装置、情報処理方法及びプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US10980411B2 (ja) |
EP (1) | EP3354191B1 (ja) |
JP (1) | JP6866167B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6949628B2 (ja) * | 2017-09-04 | 2021-10-13 | キヤノン株式会社 | 画像処理装置、光干渉断層撮像装置、画像処理方法、及びプログラム |
CN114343566B (zh) * | 2021-11-30 | 2022-11-01 | 图湃(北京)医疗科技有限公司 | 一种辅助oct进行扫描的方法、pc终端、存储介质及系统 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7301644B2 (en) * | 2004-12-02 | 2007-11-27 | University Of Miami | Enhanced optical coherence tomography for anatomical mapping |
WO2012170722A2 (en) * | 2011-06-07 | 2012-12-13 | California Institute Of Technology | Enhanced optical angiography using intensity contrast and phase contrast imaging methods |
US8433393B2 (en) * | 2011-07-07 | 2013-04-30 | Carl Zeiss Meditec, Inc. | Inter-frame complex OCT data analysis techniques |
WO2013004801A1 (en) * | 2011-07-07 | 2013-01-10 | Carl Zeiss Meditec Ag | Improved data acquisition methods for reduced motion artifacts and applications in oct angiography |
CN104411233B (zh) * | 2012-02-03 | 2017-08-29 | 俄勒冈健康科学大学 | 体内光学流成像 |
JP6062688B2 (ja) * | 2012-08-30 | 2017-01-18 | キヤノン株式会社 | 眼科装置、眼科装置の制御方法、およびプログラム |
JP6217185B2 (ja) * | 2013-07-02 | 2017-10-25 | 株式会社ニデック | 眼科撮影装置及び眼科画像処理プログラム |
JP6507615B2 (ja) * | 2013-12-13 | 2019-05-08 | 株式会社ニデック | 光コヒーレンストモグラフィ装置、及びプログラム |
JP6520099B2 (ja) * | 2014-06-30 | 2019-05-29 | 株式会社ニデック | 光コヒーレンストモグラフィ装置、及びデータ処理プログラム |
US9759544B2 (en) * | 2014-08-08 | 2017-09-12 | Carl Zeiss Meditec, Inc. | Methods of reducing motion artifacts for optical coherence tomography angiography |
JP6598503B2 (ja) * | 2015-05-01 | 2019-10-30 | キヤノン株式会社 | 画像生成装置、画像生成方法及びプログラム |
JP6606881B2 (ja) | 2015-06-16 | 2019-11-20 | 株式会社ニデック | Oct信号処理装置、oct信号処理プログラム、およびoct装置 |
JP6544071B2 (ja) * | 2015-06-17 | 2019-07-17 | 株式会社ニデック | 光コヒーレンストモグラフィ装置、および光コヒーレンストモグラフィ制御プログラム |
JP6843125B2 (ja) * | 2015-09-24 | 2021-03-17 | カール ツァイス メディテック インコーポレイテッドCarl Zeiss Meditec Inc. | 高感度フロー視覚化方法 |
US10402965B1 (en) * | 2015-11-12 | 2019-09-03 | Carl Zeiss Meditec, Inc. | Systems and methods for reducing artifacts in OCT angiography images |
JP2017158836A (ja) * | 2016-03-10 | 2017-09-14 | キヤノン株式会社 | 眼科装置および撮像方法 |
JP6703730B2 (ja) * | 2016-06-16 | 2020-06-03 | 株式会社ニデック | 光コヒーレンストモグラフィ装置、および光コヒーレンストモグラフィ制御プログラム |
-
2017
- 2017-01-18 JP JP2017006977A patent/JP6866167B2/ja active Active
-
2018
- 2018-01-16 US US15/872,647 patent/US10980411B2/en active Active
- 2018-01-17 EP EP18152041.2A patent/EP3354191B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
US10980411B2 (en) | 2021-04-20 |
JP2018114121A (ja) | 2018-07-26 |
EP3354191A1 (en) | 2018-08-01 |
EP3354191B1 (en) | 2021-07-28 |
US20180199806A1 (en) | 2018-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7250653B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP6460618B2 (ja) | 光干渉断層撮像装置およびその制御方法 | |
CN104799810B (zh) | 光学相干断层成像设备及其控制方法 | |
JP6025349B2 (ja) | 画像処理装置、光干渉断層撮像装置、画像処理方法および光干渉断層撮像方法 | |
US9918625B2 (en) | Image processing apparatus and control method of image processing apparatus | |
JP6652281B2 (ja) | 光断層撮像装置、その制御方法、及びプログラム | |
US20120253183A1 (en) | Image processing apparatus, ophthalmic apparatus, and ophthalmic system | |
US10123699B2 (en) | Ophthalmologic apparatus and imaging method | |
JP7090438B2 (ja) | 眼科撮影装置、その制御方法、プログラム、及び記録媒体 | |
JP6866167B2 (ja) | 情報処理装置、情報処理方法及びプログラム | |
JP6274728B2 (ja) | 光干渉断層撮像装置およびその制御方法 | |
JP6736734B2 (ja) | 眼科撮影装置及び眼科情報処理装置 | |
JP5975650B2 (ja) | 画像形成方法及び装置 | |
WO2020075719A1 (ja) | 画像処理装置、画像処理方法及びプログラム | |
US11074694B2 (en) | Image processing apparatus, optical coherence tomography apparatus, image processing method, and computer-readable medium | |
JP6976778B2 (ja) | 情報処理装置、情報処理方法及びプログラム | |
JP2017158867A (ja) | 眼科撮影装置 | |
US20210049742A1 (en) | Image processing apparatus, image processing method, and non-transitory computer-readable storage medium | |
JP7246862B2 (ja) | 画像処理装置、画像処理装置の制御方法及びプログラム | |
JP7024240B2 (ja) | 眼科システムおよび眼科システム制御プログラム | |
JP6486427B2 (ja) | 光干渉断層撮像装置およびその制御方法 | |
JP2019195381A (ja) | 眼科撮影装置及び眼科撮影プログラム | |
JP7013134B2 (ja) | 情報処理装置、情報処理方法及びプログラム | |
JP7086683B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP7446730B2 (ja) | 画像処理装置、画像処理方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20191217 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200525 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201020 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20201021 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201106 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210105 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210127 |
|
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: 20210309 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210407 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6866167 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |