JP7479030B2 - 血流解析装置、血流解析方法、及びプログラム - Google Patents

血流解析装置、血流解析方法、及びプログラム Download PDF

Info

Publication number
JP7479030B2
JP7479030B2 JP2020101391A JP2020101391A JP7479030B2 JP 7479030 B2 JP7479030 B2 JP 7479030B2 JP 2020101391 A JP2020101391 A JP 2020101391A JP 2020101391 A JP2020101391 A JP 2020101391A JP 7479030 B2 JP7479030 B2 JP 7479030B2
Authority
JP
Japan
Prior art keywords
blood flow
flow velocity
change
unit
cardiac cycle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020101391A
Other languages
English (en)
Other versions
JP2021194146A (ja
Inventor
晃敏 吉田
夏奈 南出
正博 秋葉
潤 酒井
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.)
Asahikawa Medical University
Topcon Corp
Original Assignee
Asahikawa Medical University
Topcon Corp
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 Asahikawa Medical University, Topcon Corp filed Critical Asahikawa Medical University
Priority to JP2020101391A priority Critical patent/JP7479030B2/ja
Priority to US17/319,089 priority patent/US20210386290A1/en
Publication of JP2021194146A publication Critical patent/JP2021194146A/ja
Application granted granted Critical
Publication of JP7479030B2 publication Critical patent/JP7479030B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A61B3/1241—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes specially adapted for observation of ocular blood flow, e.g. by fluorescein angiography
    • 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]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0285—Measuring or recording phase velocity of blood waves

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Ophthalmology & Optometry (AREA)
  • Hematology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Vascular Medicine (AREA)
  • Cardiology (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Eye Examination Apparatus (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Description

特許法第３０条第２項適用 第３６回 日本眼循環学会 講演抄録集 ６４頁 発行日 令和１年６月１４日
本発明は、血流解析装置、血流解析方法、及びプログラムに関する。
動脈硬化などの循環器系の疾患の診断及び予防には、全身循環の動態の評価は有用である。全身循環の動態は、例えば血圧波形を用いて評価される。一方、近年では、網膜の血流量と全身循環の動態との関係性を評価する手法が着目されている。
例えば、非特許文献１には、レーザースペックル血流計を用いて、視神経乳頭から出る血管内を流れる赤血球により生じるスペックルパターンの変化に基づいて血流量を測定し、測定された血流量から年齢依存性を評価する手法が開示されている。
Ｎｉｋｏｌａｕｓ Ｌｕｆｔ ｅｔ ａｌ．，"Ｏｃｕｌａｒ Ｂｌｏｏｄ Ｆｌｏｗ Ｍｅａｓｕｒｅｍｅｎｔｓ ｉｎ Ｈｅａｌｔｈｙ Ｗｈｉｔｅ Ｓｕｂｊｅｃｔｓ Ｕｓｉｎｇ Ｌａｓｅｒ Ｓｐｅｃｋｌｅ Ｆｌｏｗｇｒａｐｈｙ"， ＰＬＯＳ ＯＮＥ，２０１６年１２月１３日，ＤＯＩ：１０．１３７１／ｊｏｕｒｎａｌ．ｐｏｎｅ．０１６８１９０
しかしながら、レーザースペックル血流計で測定される血流は、脈絡膜血流と網膜血流とを含む。それにより、脈絡膜血流又は網膜血流などの眼底における所定の層領域の血流に起因した全身循環の動態を高精度に評価することができない。
本発明は、このような事情を鑑みてなされたものであり、その目的は、眼底における所定の層領域の血流に起因した全身循環の動態を高精度に評価するための新たな技術を提供することにある。
いくつかの実施形態の第１形態は、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得部と、前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出部と、を含む、血流解析装置である。
いくつかの実施形態の第２態様では、第１態様において、前記抽出部は、前記血流情報に基づいて心周期を特定する特定部を含み、前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第３態様では、第２態様において、前記抽出部は、前記血流情報に基づいて前記血流速度の振幅及び前記心周期を正規化する正規化部を含み、前記正規化部により正規化された前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第４態様では、第２態様又は第３態様において、前記抽出部は、前記心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間部を含み、前記補間部により補間された前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第５態様では、第２態様～第４態様のいずれかにおいて、前記１以上のパラメータは、前記心周期の開始タイミングから前記血流速度のピーク到達タイミングまでの間の立ち上がり時間と、前記血流速度の立ち上がり期間に前記血流速度が第１血流速度になるタイミングから前記ピーク到達タイミングまでの間の第１時間と、前記ピーク到達タイミングから前記ピーク到達タイミング以降の前記血流速度の立ち下がり期間に前記血流速度が前記第１血流速度になるタイミングまでの間の第２時間とを含む。
いくつかの実施形態の第６態様では、第１態様～第５態様のいずれかにおいて、前記１以上のパラメータは、拡張期における前記血流速度の経時的変化の波形下面積と、前記拡張期の前半における前記血流速度の経時的変化の波形下面積と、前記拡張期の後半における前記血流速度の経時的変化の波形下面積とを含む。
いくつかの実施形態の第７態様は、第１態様～第６態様のいずれかにおいて、前記１以上のパラメータに基づいて、前記血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価部を含む。
いくつかの実施形態の第８態様では、第１態様～第７態様のいずれかにおいて、前記網膜動脈又は前記網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である。
いくつかの実施形態の第９態様は、第１態様～第８態様のいずれかにおいて、略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて前記血流情報を生成する生成部を含む。
いくつかの実施形態の第１０態様は、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、を含む、血流解析方法である。
いくつかの実施形態の第１１態様では、第１０態様において、前記抽出ステップは、前記血流情報に基づいて心周期を特定する特定ステップを含み、前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第１２態様では、第１１態様において、前記抽出ステップは、前記血流情報に基づいて前記血流速度の振幅及び前記心周期を正規化する正規化ステップを含み、前記正規化ステップにおいて正規化された前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第１３態様では、第１１態様又は第１２態様において、前記抽出ステップは、前記心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間ステップを含み、前記補間ステップにおいて補間された前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する。
いくつかの実施形態の第１４態様では、第１１態様～第１３態様のいずれかにおいて、前記１以上のパラメータは、前記心周期の開始タイミングから前記血流速度のピーク到達タイミングまでの間の立ち上がり時間と、前記血流速度の立ち上がり期間に前記血流速度が第１血流速度になるタイミングから前記ピーク到達タイミングまでの間の第１時間と、前記ピーク到達タイミングから前記ピーク到達タイミング以降の前記血流速度の立ち下がり期間に前記血流速度が前記第１血流速度になるタイミングまでの間の第２時間とを含む。
いくつかの実施形態の第１５態様では、第１０態様～第１４態様のいずれかにおいて、前記１以上のパラメータは、拡張期における前記血流速度の経時的変化の波形下面積と、前記拡張期の前半における前記血流速度の経時的変化の波形下面積と、前記拡張期の後半における前記血流速度の経時的変化の波形下面積とを含む。
いくつかの実施形態の第１６態様は、第１０態様～第１５態様のいずれかにおいて、前記１以上のパラメータに基づいて、前記血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価ステップを含む。
いくつかの実施形態の第１７態様では、第１０態様～第１６態様のいずれかにおいて、前記網膜動脈又は前記網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である。
いくつかの実施形態の第１８態様は、第１０態様～第１７態様のいずれかにおいて、略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて前記血流情報を生成する生成ステップを含む。
いくつかの実施形態の第１９態様は、コンピュータに、第１０態様～第１８態様のいずれかに記載の血流解析方法の各ステップを実行させるプログラムである。
なお、上記した複数の態様に係る構成を任意に組み合わせることが可能である。
本発明によれば、眼底における所定の層領域の血流に起因した全身循環の動態を高精度に評価するための新たな技術を提供することができる。
実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の構成の一例を表す概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作を説明するための概略図である。 実施形態に係る血流解析装置の動作の一例を表すフロー図である。 実施形態に係る血流解析装置の動作の一例を表すフロー図である。 実施形態に係る血流解析装置の動作の一例を表すフロー図である。
この発明に係る血流解析装置、血流解析方法、及びプログラムの実施形態の例について、図面を参照しながら詳細に説明する。なお、この明細書において引用された文献の記載内容や任意の公知技術を、以下の実施形態に援用することが可能である。
実施形態に係る血流解析装置は、被検眼の眼底における所定の層領域における１本の動脈又は静脈の血流速度の経時的変化を表す血流情報を取得し、取得された血流情報から血流速度の変化に対応する１以上のパラメータを抽出する。例えば、抽出された１以上のパラメータに基づいて全身循環の動態が評価される。いくつかの実施形態では、動脈又は静脈は、網膜動脈又は網膜静脈である。網膜動脈又は網膜静脈は、内境界膜からブルッフ膜までの間の層領域に動脈又は静脈である。いくつかの実施形態では、動脈又は静脈は、脈絡膜動脈又は脈絡膜静脈である。いくつかの実施形態では、血流情報は、光コヒーレンストモグラフィ（ＯＣＴ）を用いて取得される。いくつかの実施形態では、血流解析装置は、血流情報を取得するためのＯＣＴを実行する機能を有する。いくつかの実施形態では、血流解析装置は、ＯＣＴを実行することにより得られた血流情報を外部から受信することで血流情報を取得する。
実施形態に係る血流解析方法は、実施形態に係る血流解析装置により実行される。実施形態に係るプログラムは、血流解析方法の各ステップをコンピュータに実行させる。
以下の実施形態では、フーリエドメインタイプのＯＣＴを実行する眼科装置が実施形態に係る血流解析装置の機能を実現する場合について説明する。
また、以下では、スウェプトソースＯＣＴと眼底カメラとを組み合わせた眼科装置について説明するが、実施形態に係る構成はこれに限定されない。例えば、ＯＣＴの種別はスウェプトソースＯＣＴには限定されず、スペクトラルドメインＯＣＴ等であってもよい。ここで、スウェプトソースＯＣＴは、波長可変光源（波長掃引光源）からの光を測定光と参照光とに分割し、被検物からの測定光の戻り光を参照光と重ね合わせて干渉光を生成し、この干渉光をバランスドフォトダイオード等で検出し、波長の掃引及び測定光のスキャンに応じて収集された検出データにフーリエ変換等を施して画像を形成する手法である。一方、スペクトラルドメインＯＣＴは、低コヒーレンス光源からの光を測定光と参照光とに分割し、被検物からの測定光の戻り光を参照光と重ね合わせて干渉光を生成し、この干渉光のスペクトル分布を分光器で検出し、検出されたスペクトル分布にフーリエ変換等を施して画像を形成する手法である。換言すると、スウェプトソースＯＣＴは時分割でスペクトル分布を取得する手法であり、スペクトラルドメインＯＣＴは空間分割でスペクトル分布を取得する手法である。なお、実施形態に適用可能なＯＣＴの手法はこれらに限定されず、他の任意の手法（例えば、タイムドメインＯＣＴ）を適用することが可能である。
実施形態に係る眼科装置は、眼底カメラのような被検眼の写真（デジタル写真）を取得する機能を備えていてもよいし、備えていなくてもよい。また、眼底カメラの代わりに、走査型レーザ検眼鏡（ＳＬＯ）や、スリットランプ顕微鏡や、前眼部撮影カメラや、手術用顕微鏡など、任意のモダリティが設けられてもよい。なお、眼底写真等の正面画像は、眼底の観察やスキャンエリアの設定やトラッキングなどに利用可能である。
＜構成＞
図１に示すように、眼科装置１は、眼底カメラユニット２、ＯＣＴユニット１００、及び演算制御ユニット２００を含む。眼底カメラユニット２には、従来の眼底カメラとほぼ同様の光学系が設けられている。ＯＣＴユニット１００には、ＯＣＴを実行するための光学系や機構が設けられている。演算制御ユニット２００はプロセッサを含む。被検者の顔を支持するための顎受けや額当てが、眼底カメラユニット２に対向する位置に設けられている。
本明細書において「プロセッサ」は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、プログラマブル論理デバイス（例えば、ＳＰＬＤ（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ）、ＣＰＬＤ（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ）、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ））等の回路を意味する。プロセッサは、例えば、記憶回路や記憶装置に格納されているプログラムを読み出し実行することで、実施形態に係る機能を実現する。
＜光学系＞
［眼底カメラユニット２］
眼底カメラユニット２には、被検眼Ｅの眼底Ｅｆを撮影するための光学系や機構が設けられている。眼底Ｅｆを撮影して得られる画像（眼底像、眼底写真等と呼ばれる）には、観察画像や撮影画像がある。観察画像は、例えば、近赤外光を用いた動画撮影により得られる。撮影画像は、例えば、可視フラッシュ光を用いて得られるカラー画像若しくはモノクロ画像、又は近赤外フラッシュ光を用いて得られるモノクロ画像である。眼底カメラユニット２は、フルオレセイン蛍光画像やインドシアニングリーン蛍光画像や自発蛍光画像などを取得可能であってよい。
眼底カメラユニット２は、照明光学系１０と撮影光学系３０とを含む。照明光学系１０は被検眼Ｅに照明光を照射する。撮影光学系３０は、被検眼Ｅからの照明光の戻り光を検出する。ＯＣＴユニット１００からの測定光は、眼底カメラユニット２内の光路を通じて被検眼Ｅに導かれ、その戻り光は、同じ光路を通じてＯＣＴユニット１００に導かれる。
照明光学系１０の観察光源１１は、例えばハロゲンランプ又は発光ダイオード（ＬＥＤ）である。観察光源１１から出力された光（観察照明光）は、曲面状の反射面を有する反射ミラー１２により反射され、集光レンズ１３を経由し、可視カットフィルタ１４を透過して近赤外光となる。更に、観察照明光は、撮影光源１５の近傍にて一旦集束し、ミラー１６により反射され、リレーレンズ１７、１８、絞り１９及びリレーレンズ２０を経由する。そして、観察照明光は、孔開きミラー２１の周辺部（孔部の周囲の領域）にて反射され、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて被検眼Ｅ（特に眼底Ｅｆ）を照明する。
被検眼Ｅからの観察照明光の戻り光は、対物レンズ２２により屈折され、ダイクロイックミラー４６を透過し、孔開きミラー２１の中心領域に形成された孔部を通過する。孔部を通過した戻り光は、ダイクロイックミラー５５を透過し、撮影合焦レンズ３１を経由し、ミラー３２により反射される。更に、この戻り光は、ハーフミラー３３Ａを透過し、ダイクロイックミラー３３により反射され、集光レンズ３４によりＣＣＤイメージセンサ３５の受光面に結像される。ＣＣＤイメージセンサ３５は、例えば所定のフレームレートで戻り光を検出する。なお、撮影光学系３０のピントが眼底Ｅｆに合っている場合には眼底Ｅｆの観察画像が得られ、ピントが前眼部に合っている場合には前眼部の観察画像が得られる。
撮影光源１５は、例えば、キセノンランプ又はＬＥＤを含む可視光源である。撮影光源１５から出力された光（撮影照明光）は、観察照明光と同様の経路を通って眼底Ｅｆに照射される。被検眼Ｅからの撮影照明光の戻り光は、観察照明光の戻り光と同じ経路を通ってダイクロイックミラー３３まで導かれ、ダイクロイックミラー３３を透過し、ミラー３６により反射され、集光レンズ３７によりＣＣＤイメージセンサ３８の受光面に結像される。
液晶ディスプレイ（ＬＣＤ）３９は、被検眼Ｅを固視させるための固視標を表示する。ＬＣＤ３９から出力された光束（固視光束）は、その一部がハーフミラー３３Ａにて反射され、ミラー３２に反射され、撮影合焦レンズ３１及びダイクロイックミラー５５を経由し、孔開きミラー２１の孔部を通過する。孔開きミラー２１の孔部を通過した固視光束は、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて眼底Ｅｆに投射される。ＬＣＤ３９の画面における固視標の表示位置を変更することにより被検眼Ｅの固視位置を変更できる。なお、ＬＣＤ３９の代わりに、複数のＬＥＤが２次元的に配列されたマトリクスＬＥＤや、光源と可変絞り（液晶絞り等）との組み合わせなどを、固視光束生成手段として用いることができる。
眼底カメラユニット２にはアライメント光学系５０とフォーカス光学系６０が設けられている。アライメント光学系５０は、被検眼Ｅに対する光学系のアライメントに用いられるアライメント指標を生成する。フォーカス光学系６０は、被検眼Ｅに対するフォーカス調整に用いられるスプリット指標を生成する。
アライメント光学系５０のＬＥＤ５１から出力されたアライメント光は、絞り５２及び５３並びにリレーレンズ５４を経由し、ダイクロイックミラー５５により反射され、孔開きミラー２１の孔部を通過する。孔開きミラー２１の孔部を通過した光は、ダイクロイックミラー４６を透過し、対物レンズ２２により被検眼Ｅに投射される。
アライメント光の角膜反射光は、対物レンズ２２、ダイクロイックミラー４６及び上記孔部を経由し、その一部がダイクロイックミラー５５を透過し、撮影合焦レンズ３１を通過する。撮影合焦レンズ３１を通過した角膜反射光は、ミラー３２により反射され、ハーフミラー３３Ａを透過し、ダイクロイックミラー３３に反射され、集光レンズ３４によりＣＣＤイメージセンサ３５の受光面に投影される。ＣＣＤイメージセンサ３５による受光像（２つの輝点からなるアライメント指標像）に基づき、従来と同様のマニュアルアライメントやオートアライメントを行うことができる。
フォーカス光学系６０は、撮影光学系３０の光路（撮影光路）に沿った撮影合焦レンズ３１の移動に連動して、照明光学系１０の光路（照明光路）に沿って移動される。反射棒６７は、照明光路に対して挿脱可能である。
フォーカス調整を行う際には、反射棒６７の反射面が照明光路に斜設される。ＬＥＤ６１から出力されたフォーカス光は、リレーレンズ６２を通過し、スプリット視標板６３により２つの光束に分離され、二孔絞り６４を通過し、ミラー６５により反射され、集光レンズ６６により反射棒６７の反射面に一旦結像されて反射される。更に、フォーカス光は、リレーレンズ２０を経由し、孔開きミラー２１に反射され、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて眼底Ｅｆに投射される。
フォーカス光の眼底反射光は、アライメント光の角膜反射光と同じ経路を通ってＣＣＤイメージセンサ３５により検出される。ＣＣＤイメージセンサ３５による受光像（２つの輝線像からなるスプリット指標像）に基づき、従来と同様のマニュアルアライメントやオートアライメントを行うことができる。
ダイクロイックミラー４６は、眼底撮影用の光路とＯＣＴ用の光路とを合成する。ダイクロイックミラー４６は、ＯＣＴに用いられる波長帯の光を反射し、眼底撮影用の光を透過させる。ＯＣＴ用の光路には、ＯＣＴユニット１００側から順に、コリメータレンズユニット４０、光路長変更部４１、光スキャナ４２、ＯＣＴ合焦レンズ４３、ミラー４４、及びリレーレンズ４５が設けられている。
光路長変更部４１は、図１に示す矢印の方向に移動可能とされ、ＯＣＴ用の光路の光路長を変更する。この光路長の変更は、被検眼Ｅの眼軸長に応じた光路長の補正や、干渉状態の調整などに利用される。光路長変更部４１は、例えばコーナーキューブと、これを移動する機構とを含む。
光スキャナ４２は、被検眼Ｅの瞳孔と光学的に共役な位置に配置される。光スキャナ４２は、ＯＣＴ用の光路を通過する測定光ＬＳの進行方向を変更する。それにより、被検眼Ｅが測定光ＬＳでスキャンされる。光スキャナ４２は、ｘｙ平面の任意方向に測定光ＬＳを偏向可能であり、例えば、測定光ＬＳをｘ方向に偏向するガルバノミラーと、ｙ方向に偏向するガルバノミラーとを含む。
［ＯＣＴユニット１００］
図２に例示するように、ＯＣＴユニット１００には、被検眼ＥのＯＣＴを実行するための光学系が設けられている。この光学系の構成は、従来のスウェプトソースＯＣＴと同様である。すなわち、この光学系は、光源からの光を測定光と参照光とに分割し、被検眼Ｅからの測定光の戻り光と参照光路を経由した参照光とを重ね合わせて干渉光を生成し、この干渉光を検出する干渉光学系を含む。干渉光学系により得られる検出結果（検出信号）は、干渉光のスペクトルを示す信号であり、演算制御ユニット２００に送られる。
光源ユニット１０１は、一般的なスウェプトソースＯＣＴと同様に、出射光の波長を高速で変化させることが可能な波長可変光源を含む。波長可変光源は、例えば、近赤外レーザ光源である。
光源ユニット１０１から出力された光Ｌ０は、光ファイバ１０２により偏波コントローラ１０３に導かれてその偏光状態が調整される。更に、光Ｌ０は、光ファイバ１０４によりファイバカプラ１０５に導かれて測定光ＬＳと参照光ＬＲとに分割される。
参照光ＬＲは、光ファイバ１１０によりコリメータ１１１に導かれて平行光束に変換され、光路長補正部材１１２及び分散補償部材１１３を経由し、コーナーキューブ１１４に導かれる。光路長補正部材１１２は、参照光ＬＲの光路長と測定光ＬＳの光路長とを合わせるよう作用する。分散補償部材１１３は、参照光ＬＲと測定光ＬＳとの間の分散特性を合わせるよう作用する。
コーナーキューブ１１４は、入射した参照光ＬＲの進行方向を逆方向に折り返す。コーナーキューブ１１４に対する参照光ＬＲの入射方向と出射方向は互いに平行である。コーナーキューブ１１４は、参照光ＬＲの入射方向に移動可能であり、それにより参照光ＬＲの光路長が変更される。
図１及び図２に示す構成では、測定光ＬＳの光路（測定光路、測定アーム）の長さを変更するための光路長変更部４１と、参照光ＬＲの光路（参照光路、参照アーム）の長さを変更するためのコーナーキューブ１１４の双方が設けられている。いくつかの実施形態では、光路長変更部４１とコーナーキューブ１１４のいずれか一方のみが設けられる。また、これら以外の光学部材を用いて、測定光路長と参照光路長との差を変更することも可能である。
コーナーキューブ１１４を経由した参照光ＬＲは、分散補償部材１１３及び光路長補正部材１１２を経由し、コリメータ１１６によって平行光束から集束光束に変換され、光ファイバ１１７に入射する。光ファイバ１１７に入射した参照光ＬＲは、偏波コントローラ１１８に導かれてその偏光状態が調整され、光ファイバ１１９によりアッテネータ１２０に導かれて光量が調整され、光ファイバ１２１によりファイバカプラ１２２に導かれる。
一方、ファイバカプラ１０５により生成された測定光ＬＳは、光ファイバ１２７により導かれてコリメータレンズユニット４０により平行光束に変換される。平行光束に変換された測定光ＬＳは、光路長変更部４１、光スキャナ４２、ＯＣＴ合焦レンズ４３、ミラー４４及びリレーレンズ４５を経由し、ダイクロイックミラー４６により反射され、対物レンズ２２により屈折されて被検眼Ｅに入射する。測定光ＬＳは、被検眼Ｅの様々な深さ位置において散乱・反射される。被検眼Ｅからの測定光ＬＳの戻り光は、往路と同じ経路を逆向きに進行してファイバカプラ１０５に導かれ、光ファイバ１２８を経由してファイバカプラ１２２に到達する。
ファイバカプラ１２２は、光ファイバ１２８を介して入射された測定光ＬＳと、光ファイバ１２１を介して入射された参照光ＬＲとを重ね合わせて干渉光を生成する。ファイバカプラ１２２は、所定の分岐比（例えば１：１）で干渉光を分岐することにより、一対の干渉光ＬＣを生成する。一対の干渉光ＬＣは、それぞれ光ファイバ１２３及び１２４を通じて検出器１２５に導かれる。
検出器１２５は、例えばバランスドフォトダイオード（Ｂａｌａｎｃｅｄ Ｐｈｏｔｏ Ｄｉｏｄｅ）である。バランスドフォトダイオードは、一対の干渉光ＬＣをそれぞれ検出する一対のフォトディテクタを有し、これらによる検出結果の差分を出力する。検出器１２５は、その検出結果（検出信号）をＤＡＱ（Ｄａｔａ Ａｃｑｕｉｓｉｔｉｏｎ Ｓｙｓｔｅｍ）１３０に送る。
ＤＡＱ１３０には、光源ユニット１０１からクロックＫＣが供給される。クロックＫＣは、光源ユニット１０１において、波長可変光源により所定の波長範囲内で掃引される各波長の出力タイミングに同期して生成される。光源ユニット１０１は、例えば、各出力波長の光Ｌ０を分岐することにより得られた２つの分岐光の一方を光学的に遅延させた後、これらの合成光を検出した結果に基づいてクロックＫＣを生成する。ＤＡＱ１３０は、検出器１２５から入力される検出信号をクロックＫＣに基づきサンプリングする。ＤＡＱ１３０は、検出器１２５からの検出信号のサンプリング結果を演算制御ユニット２００に送る。
［演算制御ユニット２００］
演算制御ユニット２００は、眼底カメラユニット２、表示装置３及びＯＣＴユニット１００の各部を制御する。また、演算制御ユニット２００は、各種の演算処理を実行する。例えば、演算制御ユニット２００は、一連の波長走査毎に（Ａライン毎に）、検出器１２５により得られた検出結果に基づくスペクトル分布にフーリエ変換等の信号処理を施すことにより、各Ａラインにおける反射強度プロファイルを形成する。更に、演算制御ユニット２００は、各Ａラインの反射強度プロファイルを画像化することにより画像データを形成する。そのための演算処理は、従来のスウェプトソースＯＣＴと同様である。
演算制御ユニット２００は、例えば、プロセッサ、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ハードディスクドライブ、通信インターフェイスなどを含む。ハードディスクドライブ等の記憶装置には各種のコンピュータプログラムが格納されている。演算制御ユニット２００は、操作デバイス、入力デバイス、表示デバイスなどを含んでいてもよい。
＜制御系＞
眼科装置１の制御系（処理系）は、演算制御ユニット２００を中心に構成される。
眼科装置１の制御系（処理系）の構成例を図３～図７に示す。図３は、眼科装置１の制御系の機能ブロック図の一例を表す。図４は、図３の画像形成部２２０の機能ブロック図の一例を表す。図５は、図３のデータ処理部２３０の機能ブロック図の一例を表す。図６は、図５の血流情報生成処理部２３１の機能ブロック図の一例を表す。図７は、図５の血流情報解析処理部２３２の機能ブロック図の一例を表す。
演算制御ユニット２００は、図３に示すように、制御部２１０と、画像形成部２２０と、データ処理部２３０とを含む。
（制御部２１０）
制御部２１０は、眼科装置１の各部を制御する。制御部２１０は、プロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブなどを含む。制御部２１０の機能は、回路を含むハードウェアと、制御ソフトウェアとの協働により実現される。制御部２１０は、主制御部２１１と記憶部２１２とを含む。
（主制御部２１１）
主制御部２１１は各種の制御を行う。特に、主制御部２１１は、眼底カメラユニット２の合焦駆動部３１Ａ及び４３Ａ、ＣＣＤイメージセンサ３５及び３８、ＬＣＤ３９、光路長変更部４１、及び光スキャナ４２などを制御する。また、主制御部２１１は、ＯＣＴユニット１００の光源ユニット１０１、参照駆動部１１４Ａ、検出器１２５、及びＤＡＱ１３０などを制御する。更に、主制御部２１１は、図１及び図２に示す光学系を駆動する光学系駆動部（図示せず）を制御する。
合焦駆動部３１Ａは、主制御部２１１からの制御を受け、撮影光学系３０の光軸に沿って撮影合焦レンズ３１を移動させる。合焦駆動部３１Ａには、撮影合焦レンズ３１を保持する保持部材と、この保持部材を移動するための駆動力を発生するアクチュエータと、この駆動力を伝達する伝達機構とが設けられる。アクチュエータは、例えばパルスモータにより構成される。伝達機構は、例えば歯車の組み合わせやラック・アンド・ピニオンなどによって構成される。それにより、主制御部２１１からの制御を受けた合焦駆動部３１Ａが撮影合焦レンズ３１を移動することにより、撮影光学系３０の合焦位置が変更される。なお、手動又はユーザーの操作部２４２に対する操作により合焦駆動部３１Ａが撮影光学系３０の光軸に沿って撮影合焦レンズ３１を移動するようにしてもよい。
合焦駆動部４３Ａは、主制御部２１１からの制御を受け、ＯＣＴユニット１００における干渉光学系の光軸（測定光の光路）に沿ってＯＣＴ合焦レンズ４３を移動させる。合焦駆動部４３Ａには、ＯＣＴ合焦レンズ４３を保持する保持部材と、この保持部材を移動するための駆動力を発生するアクチュエータと、この駆動力を伝達する伝達機構とが設けられる。アクチュエータは、例えばパルスモータにより構成される。伝達機構は、例えば歯車の組み合わせやラック・アンド・ピニオンなどによって構成される。それにより、主制御部２１１からの制御を受けた合焦駆動部４３ＡがＯＣＴ合焦レンズ４３を移動することにより、測定光の合焦位置が変更される。なお、手動又はユーザーの操作部２４２に対する操作により合焦駆動部４３Ａが干渉光学系の光軸に沿ってＯＣＴ合焦レンズ４３を移動するようにしてもよい。
主制御部２１１は、ＣＣＤイメージセンサ３５の露光時間（電荷蓄積時間）、感度、フレームレート等を制御することが可能である。主制御部２１１は、ＣＣＤイメージセンサ３８の露光時間、感度、フレームレート等を制御することが可能である。
主制御部２１１は、ＬＣＤ３９に対して固視標や視力測定用視標の表示制御を行うことが可能である。それにより、被検眼Ｅに呈示される視標の切り替えや視標の種別の変更が可能になる。また、ＬＣＤ３９における視標の表示位置を変更することにより、被検眼Ｅに対する視標呈示位置を変更することが可能である。
主制御部２１１は、光路長変更部４１を制御することにより、参照光ＬＲの光路長と測定光ＬＳの光路長との差を相対的に変更することが可能である。主制御部２１１は、被検眼Ｅの対象部位がＯＣＴ画像のフレーム内における所定の範囲に描出されるように光路長変更部４１を制御する。具体的には、主制御部２１１は、被検眼Ｅの対象部位がＯＣＴ画像のフレーム内における所定のｚ位置（深さ方向の位置）に描出されるように光路長変更部４１を制御することが可能である。
主制御部２１１は、光スキャナ４２を制御することにより被検眼Ｅの眼底Ｅｆ又は前眼部における測定光ＬＳの走査位置を変更することが可能である。光スキャナ４２の制御には、測定光ＬＳをｘ方向に偏向するガルバノミラーによる走査位置、走査範囲又は走査速度の制御、測定光ＬＳをｙ方向に偏向するガルバノミラーによる走査位置、走査範囲又は走査速度の制御などがある。また、主制御部２１１は、光スキャナ４２による測定光ＬＳの走査態様を制御することができる。測定光ＬＳの走査態様としては、例えば、水平スキャン、垂直スキャン、十字スキャン、放射スキャン、円スキャン、同心円スキャン、螺旋スキャンなどがある
主制御部２１１は、光源ユニット１０１を制御することにより、光Ｌ０の点灯と消灯の切り替えや、光Ｌ０の光量の変更などを制御することが可能である。
主制御部２１１は、参照駆動部１１４Ａを制御することにより、参照光ＬＲの光路長と測定光ＬＳの光路長との差を相対的に変更することが可能である。参照駆動部１１４Ａは、参照光路に設けられたコーナーキューブ１１４を移動させる。それにより、参照光路の長さが変更される。主制御部２１１は、被検眼Ｅの対象部位がＯＣＴ画像のフレーム内における所定の範囲に描出されるように参照駆動部１１４Ａを制御する。具体的には、主制御部２１１は、被検眼Ｅの対象部位がＯＣＴ画像のフレーム内における所定のｚ位置に描出されるように参照駆動部１１４Ａを制御することが可能である。主制御部２１１は、光路長変更部４１及び参照駆動部１１４Ａの少なくとも一方を制御することにより、参照光ＬＲの光路長と測定光ＬＳの光路長との差を相対的に変更することが可能である。以下では、主制御部２１１は、光路長変更部４１だけを制御することにより測定光ＬＳと参照光ＬＲとの光路長差調整を行うものとして説明するが、参照駆動部１１４Ａだけを制御することにより参照光ＬＲと測定光ＬＳとの光路長差調整を行ってもよい。
主制御部２１１は、検出器１２５の露光時間（電荷蓄積時間）、感度、フレームレート等を制御することが可能である。また、主制御部２１１は、ＤＡＱ１３０を制御することが可能である。
図示しない光学系駆動部は、眼科装置１に設けられた光学系（図１及び図２に示す光学系）を３次元的に移動する。主制御部２１１は、被検眼Ｅと装置光学系との位置関係を維持するように光学系駆動部を制御することが可能である。この制御は、アライメントやトラッキングにおいて用いられる。トラッキングとは、被検眼Ｅの運動に合わせて装置光学系を移動させるものである。トラッキングを行う場合には、事前にアライメントとピント合わせが実行される。トラッキングは、被検眼Ｅを動画撮影して得られる画像に基づき被検眼Ｅの位置や向きに合わせて装置光学系をリアルタイムで移動させることにより、アライメントとピントが合った好適な位置関係を維持する機能である。
（記憶部２１２）
記憶部２１２は各種のデータを記憶する。記憶部２１２に記憶されるデータとしては、例えば、ＯＣＴ画像の画像データ、眼底像の画像データ、被検眼情報などがある。被検眼情報は、被検者ＩＤや氏名などの被検者情報や、左眼／右眼の識別情報や、電子カルテ情報などを含む。
（画像形成部２２０）
画像形成部２２０は、ＤＡＱ１３０によりサンプリングされた検出器１２５からの検出信号に基づいて、眼底Ｅｆの断層像の画像データと位相画像の画像データとを形成する。画像形成部２２０はプロセッサを含む。なお、この明細書では、「画像データ」と、それに基づく「画像」とを同一視することがある。
図４に示すように、画像形成部２２０は、断層像形成部２２１と、位相画像形成部２２２とを含む。
例えば、特開２０１７－７９８８６号公報に開示されているように、実施形態では、血流計測の対象部位を設定するための走査（予備計測）と、予備計測の結果から設定された部位に関する血流情報を取得するための走査（血流計測）とが実行される。
予備計測では、眼底Ｅｆの複数の断面を測定光ＬＳで反復的に走査する。複数の断面の走査は、例えば、互いに平行な複数の走査線に沿ったラスタースキャンである。ラスタースキャンの走査線の本数は例えば１２８本であり、それにより眼底Ｅｆの３次元領域が走査される（３次元スキャン）。なお、走査パターンや走査線の本数は本例に限定されない。
血流計測では、予備計測の結果から設定された眼底Ｅｆの部位に基づき２種類の走査（第１走査及び第２走査）が実行される。予備計測の結果から設定された部位を注目部位と呼ぶ。注目部位は、注目血管と注目断面とを含む。注目血管は、予備計測の結果から設定された眼底Ｅｆの血管である。注目血管は、脈波が明瞭に観察される血管であってよく、例えば動脈である。注目断面は、注目血管に交差する断面である。更に、注目部位は、注目断面の近傍に位置する１以上の断面を含んでいてよい。この１以上の断面は、注目血管の傾きを求めるためのＯＣＴスキャンの対象となる。
第１走査では、注目血管に交差する２以上の断面が測定光ＬＳで走査される。第１走査により取得されたデータは、注目断面における注目血管の傾き（向き）を求めるために用いられる。一方、第２走査は、注目血管に交差する注目断面が測定光ＬＳで反復的に走査される。第１走査が行われる断面は、注目断面の近傍に配置される。第２走査は、ＯＣＴを用いたドップラー計測である。
第１走査及び第２走査の対象断面は、ｘｙ面において、注目血管の走行方向に対して直交するように向き付けられることが望ましい。図８の眼底像Ｄに示すように、実施形態では、例えば、視神経乳頭Ｄａの近傍に、第１走査が行われる２つの断面Ｃ１１及びＣ１２と、第２走査が行われる注目断面Ｃ２とが注目血管Ｄｂに交差するように設定される。２つの断面Ｃ１１及びＣ１２の一方は注目断面Ｃ２に対して注目血管Ｄｂの上流側に位置し、他方は下流側に位置する。注目断面Ｃ２に対する各断面Ｃ１１及びＣ１２の距離（断面間距離）は、事前に決定される。
第２走査は、被検者の心臓の少なくとも１心周期の間にわたって実行されることが望ましい。それにより、心臓の全ての時相における血流情報が得られる。第２走査の実行時間は、あらかじめ設定された一定の時間（例えば２秒間）であってもよいし、被検者ごとに又は検査毎に設定された時間であってもよい。後者の一例において、心電計を用いて得られる被検者の心拍データを利用することができる。
（断層像形成部２２１）
断層像形成部２２１は、ＤＡＱ１３０による検出信号のサンプリング結果に基づいて、眼底Ｅｆの断層像の画像データを形成する。断層像形成部２２１は、予備計測において取得されたデータに基づいて断層像を形成することができる。また、断層像形成部２２１は、血流計測（第１走査及び／又は第２走査）により取得されたデータに基づいて断層像を形成することができる。例えば、断層像形成部２２１は、断面Ｃ１１及びＣ１２に対する第１走査により得られた干渉光ＬＣの検出結果に基づいて、断面Ｃ１１の形態を表す断層像と、断面Ｃ１２の形態を表す断層像とを形成する。このとき、断面Ｃ１１を１回走査して１枚の断層像を形成し、かつ、断面Ｃ１２を１回走査して１枚の断層像を形成することができる。或いは、断面Ｃ１１を複数回走査して得られた複数の断層像に基づき１枚の断層像を取得し、かつ、断面Ｃ１２を複数回走査して得られた複数の断層像に基づき１枚の断層像を取得することができる。複数の断層像から１枚の断層像を取得する処理の例として、複数の断層像を平均して画質向上を図る処理や、複数の断層像から最適な１枚を選択する処理がある。
また、断層像形成部２２１は、注目断面Ｃ２に対する第２走査により得られた干渉光ＬＣの検出結果に基づいて、注目断面Ｃ２の形態の時系列変化を表す断層像群を形成する。この処理についてより詳しく説明する。第２走査では、上記のように注目断面Ｃ２が繰り返し走査される。断層像形成部２２１には、第２走査に応じて、ＯＣＴユニット１００の検出器１２５から検出信号が逐次入力される。断層像形成部２２１は、注目断面Ｃ２の１回分の走査に対応する検出信号群に基づいて、注目断面Ｃ２の１枚の断層像を形成する。断層像形成部２２１は、この処理を第２走査の反復回数だけ繰り返すことで、時系列に沿った一連の断層像を形成する。ここで、これら断層像を複数の群に分割し、各群の断層像を平均して画質の向上を図ってもよい。
断層像形成部２２１が実行する処理は、従来のスペクトラルドメインＯＣＴと同様に、ノイズ除去（ノイズ低減）、フィルタ処理、ＦＦＴ（Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ）などを含む。他のタイプのＯＣＴが適用される場合、断層像形成部２２１は、そのタイプに応じた公知の処理を実行する。
（位相画像形成部２２２）
位相画像形成部２２２は、眼底Ｅｆの所定部位を繰り返しＯＣＴスキャンすることにより得られたデータ（ＯＣＴデータ）に基づいて位相画像を形成する。位相画像形成部２２２は、予備計測において取得されたデータに基づいて、スキャンされた断面における位相差の時系列変化を表す位相画像を形成する。予備計測では、眼底Ｅｆの複数の断面が繰り返し走査される。位相画像形成部２２２は、予備計測により取得されたデータに基づいて、これら断面のそれぞれにおける位相画像を形成する。
また、位相画像形成部２２２は、血流計測の第２走査により取得されたデータに基づいて位相画像を形成する。例えば、位相画像形成部２２２は、注目断面Ｃ２に対する第２走査により得られた干渉光ＬＳの検出結果に基づいて、注目断面Ｃ２における位相差の時系列変化を表す位相画像を形成する。
位相画像の形成に用いられるデータは、断層像形成部２２１が断層像を形成するために用いられるデータと同じであってよい。この場合、スキャンされた断面の断層像と位相画像とを容易に位置合わせすることができる。つまり、同じデータに基づき形成された断層像と位相画像とについて、断層像の画素と位相画像の画素とを自然に対応付けることが可能である。
位相画像の形成方法の例を説明する。本例の位相画像は、隣り合うＡライン複素信号（隣接する走査点に対応する信号）の位相差を算出することにより得られる。換言すると、本例の位相画像は、スキャンされた断面の断層像の各画素について、その画素の画素値（輝度値）の時系列変化に基づき形成される。任意の画素について、位相画像形成部２２２は、その輝度値の時系列変化のグラフを考慮する。位相画像形成部２２２は、このグラフにおいて所定の時間間隔Δｔだけ離れた２つの時点ｔ１及びｔ２（ｔ２＝ｔ１＋Δｔ）の間における位相差Δφを求める。そして、この位相差Δφを時点ｔ１（より一般に２つの時点ｔ１及びｔ２の間の任意の時点）における位相差Δφ（ｔ１）として定義する。あらかじめ設定された多数の時点のそれぞれについてこの処理を実行することで、当該画素における位相差の時系列変化が得られる。
位相画像は、各画素の各時点における位相差の値を画像として表現したものである。この画像化処理は、例えば、位相差の値を表示色や輝度で表現することで実現できる。このとき、時系列に沿って位相が増加したことを表す色（例えば赤）と、減少したことを表す色（例えば青）とを違えることができる。また、位相の変化量の大きさを表示色の濃さで表現することもできる。このような表現方法を採用することで、血流の向きや大きさを色や濃度で提示することが可能となる。以上の処理を各画素について実行することにより位相画像が形成される。
なお、位相差の時系列変化は、上記の時間間隔Δｔを十分に小さくして位相の相関を確保することにより得られる。このとき、測定光ＬＳの走査において断層像の分解能に相当する時間未満の値に時間間隔Δｔを設定したオーバーサンプリングが実行される。
（データ処理部２３０）
データ処理部２３０は、画像形成部２２０により形成された画像データに対して各種のデータ処理（画像処理）や解析処理を施す。例えば、データ処理部２３０は、画像の輝度補正や分散補正等の補正処理を実行する。また、データ処理部２３０は、ＣＣＤイメージセンサ３５、３８を用い得られた画像（前眼部像等）に対して各種の画像処理や解析処理を施す。
データ処理部２３０は、断層像の間の画素を補間する補間処理などの公知の画像処理を実行することにより、被検眼Ｅのボリュームデータ（ボクセルデータ）を形成することができる。ボリュームデータに基づく画像を表示させる場合、データ処理部２３０は、このボリュームデータに対してレンダリング処理を施して、特定の視線方向から見たときの擬似的な３次元画像を形成する。
データ処理部２３０は、取得されたボリュームデータ（３次元データセット、スタックデータ等）に各種のレンダリングを施すことで、任意断面におけるＢモード画像（縦断面像、軸方向断面像）、任意断面におけるＣモード画像（横断面像、水平断面像）、プロジェクション画像、シャドウグラムなどを形成することができる。Ｂモード画像やＣモード画像のような任意断面の画像は、指定された断面上の画素（ピクセル、ボクセル）を３次元データセットから選択することにより形成される。プロジェクション画像は、３次元データセットを所定方向（ｚ方向、深さ方向、軸方向）に投影することによって形成される。シャドウグラムは、３次元データセットの一部（たとえば特定層に相当する部分データ）を所定方向に投影することによって形成される。Ｃモード画像、プロジェクション画像、シャドウグラムのような、被検眼の正面側を視点とする画像を正面画像（ｅｎ－ｆａｃｅ画像）と呼ぶ。
また、データ処理部２３０は、ＯＣＴにより時系列に収集されたデータ（例えば、Ｂスキャン画像データ）に基づいて、網膜血管や脈絡膜血管が強調されたＢモード画像や正面画像（血管強調画像、アンギオグラム）を構築することができる。例えば、被検眼Ｅの略同一部位を反復的にスキャンすることにより、時系列のＯＣＴデータを収集することができる。
いくつかの実施形態では、データ処理部２３０は、略同一部位に対するＢスキャンにより得られた時系列のＢスキャン画像を比較し、信号強度の変化部分の画素値を変化分に対応した画素値に変換することにより当該変化部分が強調された強調画像を構築する。更に、データ処理部２３０は、構築された複数の強調画像から所望の部位における所定の厚さ分の情報を抽出してｅｎ－ｆａｃｅ画像として構築することでＯＣＴアンギオグラムを形成する。
更に、データ処理部２３０は、画像形成部２２０により形成された断層像において、所定の層領域を特定することが可能である。特定可能な層領域として、内境界膜、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、網膜色素上皮層、脈絡膜、強膜、各層領域の界面などがある。
断層像から所定の層領域を特定する処理は、典型的には、セグメンテーション処理を含む。セグメンテーション処理は、断層像中の部分領域を特定するための公知の処理である。データ処理部２３０は、例えば、断層像における各画素の輝度値に基づきセグメンテーション処理を行う。すなわち、眼底Ｅｆのそれぞれの層領域は特徴的な反射率を有し、これら層領域に相当する画像領域もそれぞれ特徴的な輝度値を有する。データ処理部２３０は、これら特徴的な輝度値に基づきセグメンテーション処理を実行することにより、目的の画像領域（層領域）を特定することができる。データ処理部２３０は、少なくとも網膜動脈又は網膜静脈が存在する層領域（例えば、内境界膜からブルッフ膜までの間の層領域）を特定する。いくつかの実施形態では、データ処理部２３０は、断層像において特定された層領域に対応する位相画像中の画像領域を特定する。
データ処理部２３０は、画像形成部２２０により形成された画像（断層像、位相画像）に基づいて血流情報を生成することが可能である。具体的には、データ処理部２３０は、画像形成部２２０により形成された画像（断層像、位相画像）において上記のように特定された所定の層領域（例えば、内境界膜からブルッフ膜までの間の層領域）の血流情報を生成する。また、データ処理部２３０は、生成された血流情報を解析して１以上のパラメータを抽出し、抽出された１以上のパラメータに基づいて被検者の全身循環の動態を評価することが可能である。このようなデータ処理部２３０は、図５に示すように、血流情報生成処理部２３１と、血流情報解析処理部２３２とを含む。
（血流情報生成処理部２３１）
血流情報生成処理部２３１は、予備計測に基づいて設定された部位の血流計測にて取得されたデータの処理を実行する。本例において、血流情報生成処理部２３１は、１本の血管について、第１走査により取得されたデータに基づいて、注目断面における注目血管の傾きを求める。更に、血流情報生成処理部２３１は、第１走査に基づき求められた注目血管の傾きと、第２走査により取得されたデータとに基づいて、注目血管に関する血流情報を求める。血流情報は、例えば、血流速度や血流量を含む。このような血流情報生成処理部２３１は、図６に示すように、血管領域特定部２３１１と、傾き算出部２３１２と、血流速度算出部２３１３と、血管径算出部２３１４と、血流量算出部２３１５とを含む。
（血管領域特定部２３１１）
血管領域特定部２３１１は、断層像形成部２２１により形成された断層像において、注目血管に対応する血管領域を特定する。更に、血管領域特定部２３１１は、位相画像形成部２２２により形成された位相画像において、注目血管に対応する血管領域を特定する。血管領域の特定は、各画像の画素値を解析することにより行われる（例えば閾値処理）。いくつかの実施形態では、血管領域特定部２３１１は、断層像又は位相画像に対して、閾値処理、エッジ検出、二値化、細線化、領域拡張法（リージョングローイング）などの公知の画像処理を実行することにより血管領域を特定する。
いくつかの実施形態では、位相画像に描出される血管領域の形状が略円形状である場合、血管領域特定部２３１１は、輝度差又は位相差の変化が大きい略円形状の領域を特定することにより血管領域を特定する。
この実施形態では、同じデータから断層像と位相画像とを形成することが可能である。それにより、血管領域特定部２３１１は、断面形態が比較的明瞭に描出される断層像内の血管領域を特定し、特定された血管領域に対応する位相画像内の領域を血管領域として特定することが可能である。また、断層像と位相画像とが異なるデータから作成された場合であっても、それらのレジストレーションが直接的又は間接的に可能である場合には、一方の画像の血管領域の特定結果を他方の画像の血管領域の特定に利用することができる。
いくつかの実施形態では、血管領域特定部２３１１は、ユーザーによる操作部２４２（後述）に対する操作内容に基づいて指定された血管領域を特定する。ユーザーは、主制御部２１１により表示部２４１（後述）に表示された眼底Ｅｆの正面画像又は断面の位相画像を見ながら操作部２４２に対して操作を行うことにより、所望の注目血管に対応する血管領域を指定することができる。
（傾き算出部２３１２）
傾き算出部２３１２は、第１走査により取得されたデータに基づいて注目断面Ｃ２における注目血管Ｄｂの傾きを算出する。このとき、第２走査により得られたデータを更に用いることも可能である。傾き算出部２３１２は、断面間距離と血管領域の特定結果とに基づいて、注目断面Ｃ２における注目血管Ｄｂの傾きを算出する。断面間距離は、断面Ｃ１１と断面Ｃ１２との間の距離を含んでよい。また、断面間距離は、断面Ｃ１１と注目断面Ｃ２との間の距離と、断面Ｃ１２と注目断面Ｃ２との間の距離とを含んでよい。
注目血管Ｄｂの傾きの算出方法の例を、図９を参照しつつ説明する。断層像Ｇ１１及びＧ１２は、それぞれ、第１走査が適用される断面Ｃ１１を表す断層像及び断面Ｃ１２を表す断層像である。また、断層像Ｇ２は、第２走査が適用される注目断面Ｃ２を表す断層像である。符号Ｖ１１、Ｖ１２及びＶ２は、それぞれ、断層像Ｇ１１内の血管領域、断層像Ｇ１２内の血管領域、及び断層像Ｇ２内の血管領域を示す。なお、これら血管領域は注目血管Ｄｂの断面に相当する。図９において、ｚ座標軸は下方向を向いており、これは測定光ＬＳの照射方向（測定光ＬＳの光路の光軸の方向、軸線方向、Ａライン方向）と実質的に一致するものとする。また、隣接する断層像（断面）の間隔をＬとする。
１つの例において、傾き算出部２３１２は、３つの血管領域Ｖ１１、Ｖ１２及びＶ２の位置関係に基づいて、注目断面Ｃ２における注目血管Ｄｂの傾きＵを算出する。この位置関係は、例えば、３つの血管領域Ｖ１１、Ｖ１２及びＶ２を結ぶことによって得られる。具体的には、傾き算出部２３１２は、３つの血管領域Ｖ１１、Ｖ１２及びＶ２のそれぞれの特徴点を特定し、これら特徴点を結ぶ。この特徴点としては、中心位置、重心位置、最上部（ｚ座標値が最小の位置）、最下部（ｚ座標値が最大の位置）などがある。また、これら特徴点の結び方としては、線分で結ぶ方法、近似曲線（スプライン曲線、ベジェ曲線等）で結ぶ方法などがある。
更に、傾き算出部２３１２は、これら特徴点を結ぶ線に基づいて傾きＵを算出する。線分が用いられる場合、例えば、注目断面Ｃ２内の血管領域Ｖ２の特徴点と断面Ｃ１１内の血管領域Ｖ１１の特徴点とを結ぶ第１線分の傾きと、血管領域Ｖ２の当該特徴点と断面Ｃ１２内の血管領域Ｖ１２の特徴点とを結ぶ第２線分の傾きとに基づいて、傾きＵが算出される。この算出処理の例として、２つの線分の傾きの平均値を求めることができる。また、近似曲線で結ぶ場合の例として、近似曲線と注目断面Ｃ２との交差位置における近似曲線の傾きを求めることができる。なお、断面間距離Ｌは、線分や近似曲線を求める処理において、これら断層像Ｇ１１、Ｇ１２及びＧ２をｘｙｚ座標系に埋め込むときに用いられる。
本例では、３つの断面における血管領域を考慮しているが、２つの断面の血管領域を考慮して傾きを求めることも可能である。具体例として、断面Ｃ１１内の血管領域Ｖ１１と断面Ｃ１２内の血管領域Ｖ１２とに基づいて、注目断面Ｃ２における注目血管Ｄｂの傾きＵを求めるよう構成できる。或いは、上記第１線分又は第２線分の傾きを傾きＵとして用いることも可能である。
注目断面における注目血管の傾きを算出する方法は、上記のものには限定されない。例えば、次のような方法を適用することができる。まず、注目断面に交差し、かつ注目血管に沿う断面に対してＯＣＴスキャンを行う。次に、このＯＣＴスキャンにより取得されたデータに基づいて断層像を形成し、この断層像のセグメンテーションを行って所定組織に相当する画像領域（層領域）が特定される。特定される層領域は、例えば、内境界膜に相当する層領域（ＩＬＭ領域）である。内境界膜は、網膜と硝子体との境界を規定する網膜の組織であり、比較的明瞭に描出される。更に、特定された層領域の形状を近似する線分を求め、その傾きを求める。近似線分の傾きは、例えば、ｚ座標軸に対する角度として、又は、ｘｙ面（つまりｚ座標軸に直交する平面）に対する角度として表現される。
（血流速度算出部２３１３）
血流速度算出部２３１３は、位相画像として得られる位相差の時系列変化に基づいて、注目血管Ｄｂ内を流れる血液の注目断面Ｃ２における血流速度を算出する。この算出対象は、或る時点における血流速度でもよいし、この血流速度の時系列変化（血流速度変化情報）でもよい。前者の場合、例えば心電図の所定の時相（例えばＲ波の時相）における血流速度を選択的に取得することが可能である。また、後者における時間の範囲は、注目断面Ｃ２を走査した時間の全体又は任意の一部である。
血流速度変化情報が得られた場合、血流速度算出部２３１３は、当該時間の範囲における血流速度の統計値を算出することができる。この統計値としては、平均値、標準偏差、分散、中央値、最大値、最小値、極大値、極小値などがある。また、血流速度の値についてのヒストグラムを作成することもできる。
血流速度算出部２３１３は、前述のようにドップラーＯＣＴの手法を用いて血流速度を算出する。このとき、傾き算出部２３１２により算出された注目断面Ｃ２における注目血管Ｄｂの傾きＵが考慮される。具体的には、傾き算出部２３１２は次式を用いる。
Figure 0007479030000001
式（１）において、Δｆは、測定光ＬＳの散乱光が受けるドップラーシフトを表し、ｎは、媒質（血液）の屈折率を表し、ｖは、媒質の流速（血流速度）を表し、θは、測定光ＬＳの入射方向と媒質の流れの方向（傾きＵ）とが成す角度を表し、λは、測定光ＬＳの中心波長を表す。
実施形態では、ｎとλは既知であり、Δｆは位相差の時系列変化から得られ、θは傾きＵから得られる（又はθは傾きＵとして得られる）。これらの値を式（１）に代入することにより、血流速度ｖが算出される。
（血管径算出部２３１４）
血管径算出部２３１４は、眼底像又はＯＣＴ画像を解析することにより、注目断面Ｃ２における注目血管Ｄｂの径を算出する。眼底像が用いられる場合、注目断面Ｃ２の位置を含む眼底Ｅｆの部位の撮影が行われ、それにより得られた眼底像（例えば、カラー眼底像、レッドフリー画像等）に基づいて、注目断面Ｃ２における注目血管Ｄｂの径、つまり血管領域Ｖ２の径を算出する。ＯＣＴ画像が用いられる場合、このＯＣＴ画像は、例えば、第２走査に基づき形成された断層像、又は、予備計測に基づき形成された画像である。このような血管径の算出は、従来と同様にして実行される。
（血流量算出部２３１５）
血流量算出部２３１５は、血流速度の算出結果と血管径の算出結果とに基づいて、注目血管Ｄｂ内を流れる血液の流量を算出する。この処理の一例を以下に説明する。
血管内における血流がハーゲン・ポアズイユ流（Ｈａｇｅｎ－Ｐｏｉｓｅｕｉｌｌｅ ｆｌｏｗ）と仮定する。また、血管径をｗとし、血流速度の最大値をＶｍとすると、血流量Ｑは次式で表される。
Figure 0007479030000002
血流量算出部２３１５は、血管径算出部２３１４による血管径の算出結果ｗと、血流速度算出部２３１３による血流速度の算出結果に基づく最大値Ｖｍとを式（２）に代入することにより、血流量Ｑを算出する。
（血流情報解析処理部２３２）
血流情報解析処理部２３２は、血流情報生成処理部２３１により生成された血流情報の解析処理を実行する。本例において、血流情報解析処理部２３２は、血流情報生成処理部２３１により生成された血流情報から、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得し、取得された血流情報から被検眼の血流速度の変化に対応する１以上のパラメータを抽出する。いくつかの実施形態では、血流情報解析処理部２３２は、抽出された１以上のパラメータを評価する。このような血流情報解析処理部２３２は、図７に示すように、抽出部２３２１と、パラメータ評価部２３２２とを含む。
（抽出部２３２１）
抽出部２３２１は、血流情報生成処理部２３１により生成された血流情報からあらかじめ決められた１以上のパラメータを抽出する。抽出部２３２１は、血流情報生成処理部２３１により生成された血流情報から血流速度の経時的変化を表す血流速度変化情報（広義には、血流情報）を特定し、特定された血流速度変化情報から上記の１以上のパラメータを抽出する。このような抽出部２３２１は、図７に示すように、心周期特定部２３２１Ａと、正規化部２３２１Ｂと、補間部２３２１Ｃとを含む。
（心周期特定部２３２１Ａ）
心周期特定部２３２１Ａは、血流速度の経時的変化を表す血流速度変化情報から心周期を特定する。血流速度の経時的変化は、心臓の拍動に起因した変化である。心周期特定部２３２１Ａは、例えば、拍動に対応した血流速度の変化点を特定することにより心周期を特定する。
図１０に、抽出部２３２１の動作説明図を示す。図１０は、血流速度の経時的変化を表す血流速度波形の一例を表す。図１０において、横軸は時間を表し、縦軸は血流速度を表す。
心周期特定部２３２１Ａは、血流速度変化情報から血流速度が所定の値となる２つのタイミングの間の期間を１心拍に要する心周期ＴＨＲとして特定する。所定の値として、極小値、極大値、最小値、最大値などがある。実施形態では、心周期特定部２３２１Ａは、血流速度が最小値となる２つのタイミングの間の期間を心周期ＴＨＲとして特定する。心周期特定部２３２１Ａは、血流速度変化情報から特定された２以上の心周期の統計値を心周期ＴＨＲとして出力してもよい。この統計値としては、平均値、標準偏差、分散、中央値、最大値、最小値、極大値、極小値などがある。
（正規化部２３２１Ｂ）
正規化部２３２１Ｂは、心周期特定部２３２１Ａにより特定された心周期ＴＨＲの間の血流速度変化情報を正規化する。具体的には、正規化部２３２１Ｂは、心周期ＴＨＲに基づいて、時間軸方向（横軸方向）及び血流速度軸方向（縦軸方向）に血流速度変化情報を正規化する。時間軸方向の正規化は、心周期の正規化に相当する。血流速度軸方向の正規化や、血流速度の振幅の正規化に相当する。
正規化部２３２１Ｂは、所定の時間（例えば、１秒間又は２秒間）内に血流速度の計測サンプリング数が所定の計測サンプリング数となるように時間軸方向に正規化する。また、正規化部２３２１Ｂは、図１０に示すように、心周期ＴＨＲ内の血流速度の最大値と最小値との差ＶＨＲを最大値が「１」になり最小値が「０」になるように血流速度軸方向に正規化する。
いくつかの実施形態では、正規化により所定の計測サンプリング数を超える場合、正規化部２３２１Ｂは、１以上の計測サンプリング点を間引く。いくつかの実施形態では、正規化により所定の計測サンプリング数に満たない場合、後述の補間部２３２１Ｃは、１以上のサンプリング点を補間により追加する。
（補間部２３２１Ｃ）
補間部２３２１Ｃは、上記のように正規化により所定の計測サンプリング数に満たない場合、１以上のサンプリング点（血流速度）を補間により追加する。いくつかの実施形態では、スプライン補間等の公知の補間方法を用いて１以上のサンプリング点が追加される。
いくつかの実施形態では、正規化部２３２１Ｂにより正規化された血流速度変化情報に対し、補間部２３２１Ｃによる補間処理が省略される。いくつかの実施形態では、補間部２３２１Ｃが血流速度変化情報を補間し、正規化部２３２１Ｂが、補間された血流速度変化情報を時間軸方向及び血流速度軸方向に正規化する。
このような正規化処理及び補間処理を行うことによって、心拍数（心周期）が異なる被検者について、同一のサンプリング点数の血流速度変化情報からパラメータを抽出することが可能になる。それにより、被検眼から抽出されるパラメータの精度が向上し、パラメータに基づく評価の精度を向上させることが可能になる。
抽出部２３２１は、心周期特定部２３２１Ａ、正規化部２３２１Ｂ及び補間部２３２１Ｃにより生成された血流速度変化情報から１以上のパラメータを抽出する。１以上のパラメータは、心周期ＴＨＲ内で、基準タイミングから血流速度の変化の特徴点までに要する時間又は波形下面積を含む。
図１１Ａ及び図１１Ｂに、実施形態に係る１以上のパラメータの説明図を示す。図１１Ａ及び図１１Ｂにおいて、横軸は時間を表し、縦軸は血流速度を表す。
抽出部２３２１により抽出されるパラメータとして、波形下全面積（Ａｒｅａ ｔｏｔａｌ）Ａｔ、上行脚面積（Ａｒｅａ Ｅｌｅｖａｔｉｏｎ）Ａｅ、下降部面積（Ａｒｅａ Ｄｅｃｌｉｎａｔｉｏｎ）Ａｄがある（図１１Ａ参照）。ここで、Ａｔ＝Ａｅ＋Ａｄである。波形下全面積Ａｔは、心周期ＴＨＲの開始タイミング（時刻）から終了タイミングまでの間の波形下面積である。上行脚面積Ａｅは、心周期ＴＨＲの開始タイミングから血流速度のピーク到達タイミングＴｐまでの間の波形下面積である。下降部面積Ａｄは、血流速度のピーク到達タイミングＴｐから心周期ＴＨＲの終了タイミングまでの間の波形下面積である。
また、抽出部２３２１により抽出されるパラメータとして、収縮期面積（Ａｒｅａ Ｓｙｓｔｌｉｃ）Ｓａ、拡張期面積（Ａｒｅａ Ｄｉａｓｔｏｌｉｃ）Ｓｄ、Ｓ１部、Ｓ２部、Ｄ１部、Ｄ２部がある（図１１Ｂ参照）。ここで、Ａｔ＝Ｓａ＋Ｓｄ、Ｓａ＝Ｓ１＋Ｓ２、Ｓｄ＝Ｄ１＋Ｄ２である。収縮期面積Ｓａは、心周期ＴＨＲ内の収縮期における波形下面積である。実施形態では、心周期ＴＨＲの開始タイミングからタイミングＴｂまでの期間を収縮期とする。タイミングＴｂは、血流速度の立ち下がり期間に血流速度が最大値Ｖ１００（ピーク到達タイミングＴｐ）の半分のＶ５０になるタイミングである。拡張期面積Ｓｄは、心周期ＴＨＲ内の拡張期における波形下面積（拡張期における血流速度の経時的変化の波形下面積）である。実施形態では、タイミングＴｂから心周期ＴＨＲの終了タイミングまでの期間を拡張期とする。Ｓ１部は、心周期ＴＨＲ（収縮期）の開始タイミングからピーク到達タイミングＴｐまでの間の波形下面積である。Ｓ２部は、心周期ＴＨＲ（収縮期）のピーク到達タイミングＴｐからタイミングＴｂまでの間の波形下面積である。Ｄ１部は、心周期ＴＨＲ（拡張期）内のタイミングＴｂからタイミングＴｃまでの間の波形下面積（拡張期の前半における血流速度の経時的変化の波形下面積）である。タイミングＴｂから心周期ＴＨＲの終了タイミングまでの間の期間Ｔｄとしたとき、タイミングＴｃは、タイミングＴｂからＴｄ／２だけ経過したタイミングである（Ｔｃ＝Ｔｂ＋Ｔｄ／２）。Ｄ２部は、心周期ＴＨＲ（拡張期）内のタイミングＴｃから心周期ＴＨＲの終了タイミングまでの間の波形下面積（拡張期の後半における血流速度の経時的変化の波形下面積）である。なお、タイミングＴｂ、Ｔｄの少なくとも一方は任意に変更可能であってよい。
更に、抽出部２３２１により抽出されるパラメータとして、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、立ち上がり時間（Ｕｐｓｔｒｏｋｅ Ｔｉｍｅ）（幅）Ｔｕｐがある。半値幅上行期間Ｔ１は、血流速度の立ち上がり期間（収縮期）に血流速度がＶ５０（第１血流速度）になるタイミングＴａからピーク到達タイミングＴｐまでの間の時間（第１時間）である。半値幅下降期間Ｔ２は、ピーク到達タイミングＴｐからピーク到達タイミング以降の血流速度の立ち下がり期間に血流速度がＶ５０になるタイミングＴｂまでの間の時間（第２時間）である。立ち上がり時間Ｔｕｐは、心周期ＴＨＲの開始タイミングから血流速度のピーク到達タイミングＴｐまでの間の時間である。
いくつかの実施形態では、抽出部２３２１は、収縮期面積Ｓａ、拡張期面積Ｓｄ、Ｓ１部、Ｓ２部、Ｄ１部、及びＤ２部のうち少なくとも１つを用いた新たなパラメータを求める。例えば、抽出部２３２１は、収縮期面積Ｓａに対する収縮期面積Ｓａ又は拡張期面積Ｓｄの比、収縮期面積Ｓａと拡張期面積Ｓｄとの比、収縮期面積Ｓａに対するＳ１部又はＳ２部の比、Ｓ１部とＳ２部との比、拡張期面積Ｓｄに対するＤ１部又はＤ２部の比、Ｄ１部とＤ２部との比などを求める。
いくつかの実施形態では、抽出部２３２１は、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、及び立ち上がり時間Ｔｕｐのうち少なくとも１つを用いた新たなパラメータを求める。例えば、抽出部２３２１は、立ち上がり時間Ｔｕｐに対する半値幅上行期間Ｔ１の比、半値幅上行期間Ｔ１と半値幅下降期間Ｔ２との比、半値幅上行期間Ｔ１と半値幅下降期間Ｔ２との和に対する半値幅上行期間Ｔ１又は半値幅下降期間Ｔ２の比（Ｔ１／（Ｔ１＋Ｔ２）又はＴ２／（Ｔ１＋Ｔ２））などを求める。
実施形態に係る抽出部２３２１によるパラメータの抽出手法は、上記の手法に限定されるものではない。抽出部２３２１は、フーリエ変換、微分法、ピーク検出、面積計算、モーメント計算等による公知の波形処理手法を用いて１以上のパラメータを抽出することが可能である。
（パラメータ評価部２３２２）
パラメータ評価部２３２２は、抽出部２３２１により抽出された１以上のパラメータに基づいて、所定の評価処理を実行する。例えば、複数の被検眼に対するＯＣＴ計測によって得られた複数の血流速度の経時的変化を事前に解析することにより得られた所定の評価項目の時間依存性に基づいて、パラメータ評価部２３２２は、抽出部２３２１により抽出された１以上のパラメータに基づいて、被検者の全身循環の動態に関する評価を行う。
いくつかの実施形態では、パラメータ評価部２３２２は、被検眼の血流速度の経時的変化が加齢に伴う循環動態の変化に起因した加齢性変化であるか否かを評価する。いくつかの実施形態では、パラメータ評価部２３２２は、被検眼の血流速度の変化から全身循環の動態の程度（例えば、動脈硬化指標、心機能に関する指標）を評価する。
複数の血流速度の経時的変化に対する事前の解析には、回帰分析（単回帰分析、重回帰分析）、ランダムフォレストやサポートベクタマシン等の機械学習等の公知の解析手法が用いられる。この場合、パラメータ評価部２３２２は、回帰分析により得られた回帰直線又は回帰曲線を用いて、抽出部２３２１により抽出された１以上のパラメータに対応した評価結果を取得することができる。或いは、パラメータ評価部２３２２は、機械学習により得られた学習済みモデルを用いて、抽出部２３２１により抽出された１以上のパラメータに対応した評価結果を取得することができる。
図１２～図１５に、抽出部２３２１により抽出されるパラメータと年齢との関係を示す。図１２～図１５は、年齢５０．０±２４．３才の３９例３９眼（男女比は２１：１８）を対象とした無散瞳下でのＯＣＴ計測結果から求められたパラメータと年齢との関係を表す。但し、眼疾患又は心疾患の既往歴がある被検者、高血圧（収縮期血圧が１４０以上、及び／又は、拡張期血圧が９０以上）の被検者は除外されている。
図１２は、（ａ）波形下全面積Ａｔと年齢との相関関係、（ｂ）上行脚面積Ａｅと年齢との相関関係、（ｃ）下降部面積Ａｄと年齢との相関関係を表す。図１２の（ａ）～（ｃ）において、横軸は年齢を表し、縦軸は任意単位を表す。
図１２の（ａ）～（ｃ）の相関係数に相当するｒ値及び有意確率に相当するＰ値から、少なくとも上行脚面積Ａｅは年齢と正の相関関係があり、下降部面積Ａｄは年齢と負の相関関係がある。例えば、複数の上行脚面積Ａｅ又は下降部面積Ａｄにおける回帰直線（回帰曲線）を特定することで、被検眼の血流速度の経時的変化に基づいて抽出された上行脚面積Ａｅ又は下降部面積Ａｄから、被検眼の血流速度の変化が加齢性変化であるか否かを評価することが可能である。また、上行脚面積Ａｅ又は下降部面積Ａｄと動脈硬化指標（又は心機能に関する指標）と関連付けられる場合、被検眼の血流速度の経時的変化に基づいて抽出された上行脚面積Ａｅ又は下降部面積Ａｄから全身循環の動態の程度を評価することが可能である。
図１３は、（ａ）収縮期面積Ｓａと年齢との相関関係、（ｂ）Ｓ１部と年齢との相関関係、（ｃ）Ｓ２部と年齢との相関関係を表す。図１３の（ａ）～（ｃ）において、横軸は年齢を表し、縦軸は任意単位を表す。
図１３の（ａ）～（ｃ）のｒ値及びＰ値から、少なくともＳ１部は年齢と正の相関関係がある。例えば、複数のＳ１部における回帰直線（回帰曲線）を特定することで、被検眼の血流速度の経時的変化に基づいて抽出されたＳ１部から、被検眼の血流速度の変化が加齢性変化であるか否かを評価することが可能である。また、Ｓ１部と動脈硬化指標（又は心機能に関する指標）と関連付けられる場合、被検眼の血流速度の経時的変化に基づいて抽出されたＳ１部から全身循環の動態の程度を評価することが可能である。
図１４は、（ａ）拡張期面積Ｓｄと年齢との相関関係、（ｂ）Ｄ１部と年齢との相関関係、（ｃ）Ｄ２部と年齢との相関関係を表す。図１４の（ａ）～（ｃ）において、横軸は年齢を表し、縦軸は任意単位を表す。
図１４の（ａ）～（ｃ）のｒ値及びＰ値から、拡張期面積Ｓｄは年齢と負の相関関係があり、Ｄ１部は年齢と負の相関関係があり、Ｄ２部は年齢と負の相関関係がある。例えば、複数の拡張期面積Ｓｄ、Ｄ１部又はＤ２部における回帰直線（回帰曲線）を特定することで、被検眼の血流速度の経時的変化に基づいて抽出された拡張期面積Ｓｄ、Ｄ１部又はＤ２部から、被検眼の血流速度の変化が加齢性変化であるか否かを評価することが可能である。また、拡張期面積Ｓｄ、Ｄ１部又はＤ２部部と動脈硬化指標（又は心機能に関する指標）と関連付けられる場合、被検眼の血流速度の経時的変化に基づいて抽出された拡張期面積Ｓｄ、Ｄ１部又はＤ２部から全身循環の動態の程度を評価することが可能である。
図１５は、（ａ）立ち上がり時間Ｔｕｐと年齢との相関関係、（ｂ）半値幅上行期間Ｔ１と年齢との相関関係、（ｃ）半値幅下降期間Ｔ２と年齢との相関関係、（ｄ）Ｔ１／（Ｔ１＋Ｔ２）と年齢との相関関係を表す。図１５の（ａ）～（ｄ）において、横軸は年齢を表し、縦軸は任意単位を表す。
図１５の（ａ）～（ｃ）のｒ値及びＰ値から、立ち上がり時間Ｔｕｐは年齢と正の相関関係があり、半値幅上行期間Ｔ１は年齢と正の相関関係があり、半値幅下降期間Ｔ２は年齢と負の相関関係があり、Ｔ１／（Ｔ１＋Ｔ２）は年齢と正の相関関係がある。
例えば、複数の立ち上がり時間Ｔｕｐ、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、又はＴ１／（Ｔ１＋Ｔ２）における回帰直線（回帰曲線）を特定することで、被検眼の血流速度の経時的変化に基づいて抽出された立ち上がり時間Ｔｕｐ、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、又はＴ１／（Ｔ１＋Ｔ２）から、被検眼の血流速度の変化が加齢性変化であるか否かを評価することが可能である。また、立ち上がり時間Ｔｕｐ、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、又はＴ１／（Ｔ１＋Ｔ２）と動脈硬化指標（又は心機能に関する指標）と関連付けられる場合、被検眼の血流速度の経時的変化に基づいて抽出された立ち上がり時間Ｔｕｐ、半値幅上行期間Ｔ１、半値幅下降期間Ｔ２、又はＴ１／（Ｔ１＋Ｔ２）から全身循環の動態の程度を評価することが可能である。
いくつかの実施形態では、パラメータ評価部２３２２は、上記のパラメータの少なくとも２つのパラメータのそれぞれのパラメータに対する評価結果に基づいて、被検者の全身循環の動態に関する評価を行う。例えば、パラメータ評価部２３２２は、第１パラメータ（例えば、Ｓ１部）に対する評価結果と、第２パラメータ（例えば、半値幅上行期間Ｔ１）に対する評価結果とに基づいて、被検眼の血流速度の変化が加齢性変化であるか否かを評価したり、全身循環の動態の程度を評価したりすることが可能である。
（ユーザーインターフェイス２４０）
図３に示すように、眼科装置１には、ユーザーインターフェイス２４０が設けられている。ユーザーインターフェイス２４０は表示部２４１と操作部２４２とを含む。表示部２４１は表示装置３を含む。操作部２４２は各種の操作デバイスや入力デバイスを含む。ユーザーインターフェイス２４０は、例えばタッチパネルのような表示機能と操作機能とが一体となったデバイスを含んでいてもよい。ユーザーインターフェイス２４０の少なくとも一部を含まない実施形態を構築することも可能である。例えば、表示デバイスは、眼科装置に接続された外部装置であってよい。
眼科装置１は、実施形態に係る「血流解析装置」の一例である。血流情報解析処理部２３２において外部から血流情報を受信するインターフェイス処理部（不図示）、又は図１において対物レンズ２２からＯＣＴユニット１００までの光学系と演算制御ユニット２００（画像形成部２２０、データ処理部２３０（血流情報生成処理部２３１））は、実施形態に係る「取得部」の一例である。心周期特定部２３２１Ａは、実施形態に係る「特定部」の一例である。パラメータ評価部２３２２は、実施形態に係る「評価部」の一例である。血流情報生成処理部２３１は、実施形態に係る「生成部」の一例である。
［動作］
実施形態に係る眼科装置１の動作について説明する。
眼科装置１の動作の例を図１６～図１８に示す。記憶部２１２には、図１６～図１８に示す処理を実現するためのコンピュータプログラムが記憶されている。主制御部２１１は、このコンピュータプログラムに従って動作することにより、図１６～図１８に示す処理を実行する。なお、アライメント、フォーカス調整、干渉感度調整、ｚ位置調整など、一般的な準備処理は既に完了しているものとする。
（Ｓ１：血流情報を生成）
まず、制御部２１０（主制御部２１１）は、被検眼の１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を血流情報生成処理部２３１に生成させる。
具体的には、制御部２１０は、ＯＣＴユニット１００等の光学系を制御してＯＣＴ計測を実行させ、取得された計測結果を用いて血流情報を血流情報生成処理部２３１生成させる。ステップＳ１の詳細は後述する。
（Ｓ２：パラメータを抽出）
次に、制御部２１０は、ステップＳ１において得られた血流情報から上記の１以上パラメータを抽出部２３２１に抽出させる。ステップＳ２の詳細は後述する。
（Ｓ３：パラメータを評価）
続いて、制御部２１０は、ステップＳ２において得られたパラメータを用いて所定の評価処理をパラメータ評価部２３２２に実行させる。
図１２～図１５に示すように、所定の評価処理に対応したパラメータと年齢との相関関係について回帰直線（回帰曲線）があらかじめ特定されているものとする。例えば、パラメータ評価部２３２２は、この回帰直線に対するパラメータのずれが所定の範囲内であるときに被検眼の血流速度の変化が加齢性変化であると評価し、当該ずれが所定の範囲外であるときに被検眼の血流速度の変化が加齢性変化ではないと評価する。
以上で、眼科装置１の動作は終了である（エンド）。
図１６のステップＳ１の処理は、図１７に示すように行われる。
（Ｓ１１：計測位置を指定）
まず、ユーザー又は制御部２１０は、予備計測のスキャン範囲を設定する。スキャン範囲の設定においては、例えば、リアルタイムで取得される眼底Ｅｆの赤外観察像、或いは、過去に取得された眼底Ｅｆのデータ（眼底像、ＯＣＴ画像、ＳＬＯ画像、位相画像等）が参照される。
制御部２１０は、ＯＣＴユニット１００等を制御することにより、設定されたスキャン範囲（複数の断面）に対して予備計測を実行させる。一例として、設定されるスキャン範囲は３次元領域であり、この３次元領域のラスタースキャンが繰り返し実行される。位相画像形成部２２２は、予備計測により取得されたデータに基づいて、複数の断面のそれぞれについて位相画像を形成する。
（Ｓ１２：計測位置近傍の断面を設定）
次に、制御部２１０は、眼底Ｅｆの正面画像を表示させる。正面画像は、例えば、リアルタイムで取得される赤外観察像、又は、過去に取得された画像である。ユーザーは、表示された情報を参照し、血流計測の対象となる断面（注目断面Ｃ２）を設定する。例えば、ユーザーは、眼底Ｅｆの正面画像とともに表示された設定支援情報を参照することにより、この正面画像内の所望の動脈領域（注目血管Ｄｂ）を指定する。
本例では注目断面Ｃ２の設定を手動で行っているが、これに限定されない。例えば、眼底Ｅｆの赤外観察像と計測部位情報とに基づいて、注目血管Ｄｂや注目断面Ｃ２を自動で設定することができる。注目血管Ｄｂの設定は、例えば、動脈領域のいずれか（例えば最も太いもの）を選択することにより行われる。また、注目断面Ｃ２は、例えば、注目血管Ｄｂの向きが所定の許容範囲内である位置に、注目血管Ｄｂの走行方向に対して直交するように設定される。
更に、ユーザー又は制御部２１０（又はデータ処理部２３０）は、第１走査の対象となる断面Ｃ１１及びＣ１２の設定を行う。
（Ｓ１３：ＯＣＴ計測）
続いて、制御部２１０は、ステップＳ１２において設定された断面Ｃ１１及びＣ１２のＯＣＴスキャンを実行する（第１走査）。更に、制御部２１０は、ステップＳ１２において設定された注目断面Ｃ２の反復的なＯＣＴスキャンを実行する（第２走査）。
（Ｓ１４：画像を形成）
制御部２１０は、ステップＳ１３における第１走査により取得されたデータに基づいて、断面Ｃ１１及びＣ１２に対応する断層像Ｇ１１及びＧ１２を断層像形成部２２１に形成させる。また、制御部２１０は、ステップＳ１３における第２走査により取得されたデータに基づいて、注目断面Ｃ２の位相画像を位相画像形成部２２２に形成させる。また、断層像形成部２２１は、当該データに基づいて注目断面Ｃ２の断層像を形成する。
（Ｓ１５：血管領域を特定）
制御部２１０は、ステップＳ１４において形成された断層像における注目血管Ｄｂを特定させる。このとき、制御部２１０は、データ処理部２３０によるセグメンテーション処理により特定された所定の層領域（例えば、内境界膜からブルッフ膜までの間の層領域）における注目血管Ｄｂに相当する画像領域を血管領域として特定する。
（Ｓ１６：血管の傾きを算出）
続いて、制御部２１０は、上記のように注目断面Ｃ２における注目血管Ｄｂの傾きＵを傾き算出部２３１２に算出させる。
（Ｓ１７：血流速度を算出）
続いて、制御部２１０は、ステップＳ１３における第１走査に基づき算出された傾きＵと、ステップＳ１３における第２走査により取得された位相画像とに基づいて、注目断面Ｃ２における血流速度を血流速度算出部２３１３に算出させる。
（Ｓ１８：血管径を算出）
次に、制御部２１０は、注目断面Ｃ２における注目血管Ｄｂの径を血管径算出部２３１４に算出させる。
（Ｓ１９：血流量を算出）
続いて、制御部２１０は、ステップＳ１７において算出された血流速度と、ステップＳ１８において算出された血管径とに基づいて、注目血管Ｄｂ内を流れる血液の流量を血流量算出部２３１５に算出させる。
制御部２１０は、少なくとも１心周期より長い期間にわたる時系列の血流速度を含む血流情報を血流情報生成処理部２３１に生成させる。
以上で、図１６のステップＳ１の処理は終了である（エンド）。
図１６のステップＳ２の処理は、図１８に示すように行われる。
（Ｓ２１：心周期を特定）
まず、制御部２１０は、図１６のステップＳ１において得られた血流情報から血流速度の経時的変化を表す血流情報を特定し、特定された血流情報から上記のように心周期を心周期特定部２３２１Ａに特定させる。
（Ｓ２２：正規化）
次に、制御部２１０は、ステップＳ２１において特定された心周期における血流速度の経時的変化を上記のように正規化部２３２１Ｂに正規化させる。
（Ｓ２３：補間）
続いて、制御部２１０は、ステップＳ２２において正規化された血流速度の経時的変化を上記のように補間部２３２１Ｃに補間させる。
（Ｓ２４：パラメータを特定）
制御部２１０は、ステップＳ２３において補間された血流速度の経時的変化から図１１Ａ及び図１１Ｂに示すパラメータを抽出部２３２１に特定させる。抽出部２３２１は、特定されたパラメータをパラメータ評価部２３２２に出力する。
パラメータ評価部２３２２は、ステップＳ３において、ステップＳ２４において抽出されたパラメータを用いて上記の評価処理を実行する。
以上で、図１６のステップＳ２の処理は終了である（エンド）。
［作用・効果］
実施形態に係る血流解析装置、血流解析方法、及びプログラムの作用及び効果について説明する。
いくつかの実施形態に係る血流解析装置（眼科装置１）は、取得部（血流情報解析処理部２３２において外部から血流情報を受信するインターフェイス処理部、又は対物レンズ２２からＯＣＴユニット１００までの光学系と演算制御ユニット２００（画像形成部２２０、データ処理部２３０（血流情報生成処理部２３１））と、抽出部（２３２１）とを含む。取得部は、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報（血流速度変化情報）を取得する。抽出部は、血流情報から血流速度の変化に対応する１以上のパラメータを抽出する。
このような構成によれば、１本の網膜動脈又は網膜静脈の血流速度の経時的変化から血流速度の変化に対応する１以上のパラメータが抽出される。それにより、眼底における網膜動脈又は網膜静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態では、抽出部は、血流情報に基づいて心周期を特定する特定部（心周期特定部２３２１Ａ）を含み、心周期内の血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような構成によれば、心周期内での血流速度の変化に対応する１以上のパラメータを抽出するようにしたので、心臓の拍動又は血管の状態等に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態では、抽出部は、血流情報に基づいて血流速度の振幅及び心周期を正規化する正規化部（２３２１Ｂ）を含み、正規化部により正規化された血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような構成によれば、血流速度及び心周期が異なる複数の被検者にとって有用なパラメータを抽出することができる。それにより、複数の被検者（被検眼）の血流速度の経時的変化を用いて全身循環の動態を高精度に評価することが可能になる。
いくつかの実施形態では、抽出部は、心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間部（２３２１Ｃ）を含み、補間部により補間された血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような構成によれば、血流速度及び心周期が異なる複数の被検者に対して一定の計測サンプリング数の血流情報からパラメータを抽出することができる。それにより、複数の被検者（被検眼）の血流速度の経時的変化を用いて全身循環の動態を高精度に評価することが可能になる。
いくつかの実施形態では、１以上のパラメータは、心周期の開始タイミングから血流速度のピーク到達タイミングまでの間の立ち上がり時間（Ｔｕｐ）と、血流速度の立ち上がり期間に血流速度が第１血流速度（Ｖ５０）になるタイミングからピーク到達タイミング（Ｔｐ）までの間の第１時間（Ｔ１）と、ピーク到達タイミングからピーク到達タイミング以降の血流速度の立ち下がり期間に血流速度が第１血流速度になるタイミングまでの間の第２時間（Ｔ２）とを含む。
このような構成によれば、上記の立ち上がり時間、第１時間、又は第２時間を用いて被検者の全身循環の動態を高精度に評価するための血流解析装置を提供することが可能になる。
いくつかの実施形態では、１以上のパラメータは、拡張期における血流速度の経時的変化の波形下面積（Ｓｄ）と、拡張期の前半における血流速度の経時的変化の波形下面積（Ｄ１部）と、拡張期の後半における血流速度の経時的変化の波形下面積（Ｄ２部）とを含む。
このような構成によれば、波形下面積Ｓｄ、Ｄ１部、又はＤ２部を用いて被検者の全身循環の動態を高精度に評価するための血流解析装置を提供することが可能になる。
いくつかの実施形態は、１以上のパラメータに基づいて、血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価部（パラメータ評価部２３２２）を含む。
このような構成によれば、血流速度の変化が加齢性変化であるか否かを高精度に評価するための血流解析装置を提供することが可能になる。
いくつかの実施形態では、網膜動脈又は網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である。
このような構成によれば、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態は、略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて血流情報を生成する生成部（血流情報生成処理部２３１）を含む。
このような構成によれば、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を生成し、生成された血流情報から血流速度の変化に対応する１以上のパラメータが抽出される。それにより、眼底における網膜動脈又は網膜静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能な血流解析装置を提供することができる。
いくつかの実施形態に係る血流解析方法は、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、を含む。
このような方法によれば、１本の網膜動脈又は網膜静脈の血流速度の経時的変化から血流速度の変化に対応する１以上のパラメータが抽出される。それにより、眼底における網膜動脈又は網膜静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態では、抽出ステップは、血流情報に基づいて心周期を特定する特定ステップを含み、心周期内の血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような方法によれば、心周期内での血流速度の変化に対応する１以上のパラメータを抽出するようにしたので、心臓の拍動又は血管の状態等に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態では、抽出ステップは、血流情報に基づいて血流速度の振幅及び心周期を正規化する正規化ステップを含み、正規化ステップにおいて正規化された血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような方法によれば、血流速度及び心周期が異なる複数の被検者にとって有用なパラメータを抽出することができる。それにより、複数の被検者（被検眼）の血流速度の経時的変化を用いて全身循環の動態を高精度に評価することが可能になる。
いくつかの実施形態では、抽出ステップは、心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間ステップを含み、補間ステップにおいて補間された血流速度の経時的変化に基づいて１以上のパラメータを抽出する。
このような方法によれば、血流速度及び心周期が異なる複数の被検者に対して一定の計測サンプリング数の血流情報からパラメータを抽出することができる。それにより、複数の被検者（被検眼）の血流速度の経時的変化を用いて全身循環の動態を高精度に評価することが可能になる。
いくつかの実施形態では、１以上のパラメータは、心周期の開始タイミングから血流速度のピーク到達タイミングまでの間の立ち上がり時間（Ｔｕｐ）と、血流速度の立ち上がり期間に血流速度が第１血流速度（Ｖ５０）になるタイミングからピーク到達タイミング（Ｔｐ）までの間の第１時間（Ｔ１）と、ピーク到達タイミングからピーク到達タイミング以降の血流速度の立ち下がり期間に血流速度が第１血流速度になるタイミングまでの間の第２時間（Ｔ２）とを含む。
このような方法によれば、上記の立ち上がり時間、第１時間、又は第２時間を用いて被検者の全身循環の動態を高精度に評価するための血流解析方法を提供することが可能になる。
いくつかの実施形態では、１以上のパラメータは、拡張期における血流速度の経時的変化の波形下面積（Ｓｄ）と、拡張期の前半における血流速度の経時的変化の波形下面積（Ｄ１部）と、拡張期の後半における血流速度の経時的変化の波形下面積（Ｄ２部）とを含む。
このような方法によれば、波形下面積Ｓｄ、Ｄ１部、又はＤ２部を用いて被検者の全身循環の動態を高精度に評価するための血流解析方法を提供することが可能になる。
いくつかの実施形態は、１以上のパラメータに基づいて、血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価ステップを含む。
このような方法によれば、血流速度の変化が加齢性変化であるか否かを高精度に評価するための血流解析方法を提供することが可能になる。
いくつかの実施形態では、網膜動脈又は網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である。
このような方法によれば、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いくつかの実施形態は、略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて血流情報を生成する生成ステップを含む。
このような方法によれば、１本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を生成し、生成された血流情報から血流速度の変化に対応する１以上のパラメータが抽出される。それにより、眼底における網膜動脈又は網膜静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能な血流解析方法を提供することができる。
いくつかの実施形態に係るプログラムは、コンピュータに、上記のいずれかに記載の血流解析方法の各ステップを実行させる。
このようなプログラムによれば、１本の網膜動脈又は網膜静脈の血流速度の経時的変化から血流速度の変化に対応する１以上のパラメータが抽出される。それにより、眼底における網膜動脈又は網膜静脈の血流に起因した全身循環の動態を高精度に評価するためのパラメータを抽出することが可能になる。
いずれかの実施形態に係るプログラムを記録したコンピュータ可読な非一時的記録媒体を作成することが可能である。この非一時的記録媒体は任意の形態であってよく、その例として、磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどがある。
以上に説明した実施形態は本発明の一例に過ぎない。本発明を実施しようとする者は、本発明の要旨の範囲内における変形（省略、置換、付加等）を任意に施すことが可能である。
１ 眼科装置
１００ ＯＣＴユニット
２１０ 制御部
２１１ 主制御部
２２０ 画像形成部
２２１ 断層像形成部
２２２ 位相画像形成部
２３０ データ処理部
２３１ 血流情報生成処理部
２３１１ 血管領域特定部
２３１２ 傾き算出部
２３１３ 血流速度算出部
２３１４ 血管径算出部
２３１５ 血流量算出部
２３２ 血流情報解析処理部
２３２１ 抽出部
２３２１Ａ 心周期特定部
２３２１Ｂ 正規化部
２３２１Ｃ 補間部
２３２２ パラメータ評価部
Ｅｆ 眼底

Claims (15)

  1. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得部と、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出部と、
    を含み、
    前記抽出部は、
    前記血流情報に基づいて心周期を特定する特定部と、
    前記血流情報に基づいて前記血流速度の振幅及び前記心周期を正規化する正規化部と、
    を含み、前記正規化部により正規化された前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する、血流解析装置。
  2. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得部と、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出部と、
    を含み、
    前記抽出部は、
    前記血流情報に基づいて心周期を特定する特定部と、
    前記心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間部と、
    を含み、前記補間部により補間された前記血流速度であって前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する、血流解析装置。
  3. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得部と、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出部と、
    を含み、
    前記抽出部は、前記血流情報に基づいて心周期を特定する特定部を含み、前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出し、
    前記１以上のパラメータは、前記心周期の開始タイミングから前記血流速度のピーク到達タイミングまでの間の立ち上がり時間と、前記血流速度の立ち上がり期間に前記血流速度が第１血流速度になるタイミングから前記ピーク到達タイミングまでの間の第１時間と、前記ピーク到達タイミングから前記ピーク到達タイミング以降の前記血流速度の立ち下がり期間に前記血流速度が前記第１血流速度になるタイミングまでの間の第２時間とを含む、血流解析装置。
  4. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得部と、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出部と、
    を含み、
    前記１以上のパラメータは、拡張期における前記血流速度の経時的変化の波形下面積と、前記拡張期の前半における前記血流速度の経時的変化の波形下面積と、前記拡張期の後半における前記血流速度の経時的変化の波形下面積とを含む、血流解析装置。
  5. 前記１以上のパラメータに基づいて、前記血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価部を含む
    ことを特徴とする請求項１～請求項４のいずれか一項に記載の血流解析装置。
  6. 前記網膜動脈又は前記網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である
    ことを特徴とする請求項１～請求項５のいずれか一項に記載の血流解析装置。
  7. 略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて前記血流情報を生成する生成部を含む
    ことを特徴とする請求項１～請求項６のいずれか一項に記載の血流解析装置。
  8. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、
    を含み、
    前記抽出ステップは、
    前記血流情報に基づいて心周期を特定する特定ステップと、
    前記血流情報に基づいて前記血流速度の振幅及び前記心周期を正規化する正規化ステップと、
    を含み、前記正規化ステップにおいて正規化された前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する、血流解析方法。
  9. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、
    を含み、
    前記抽出ステップは、
    前記血流情報に基づいて心周期を特定する特定ステップと、
    前記心周期内の血流速度が所定の計測サンプリング数となるように血流速度を補間する補間ステップと、
    を含み、前記補間ステップにおいて補間された前記血流速度であって前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出する、血流解析方法。
  10. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、
    を含み、
    前記抽出ステップは、前記血流情報に基づいて心周期を特定する特定ステップを含み、前記心周期内の前記血流速度の経時的変化に基づいて前記１以上のパラメータを抽出し、
    前記１以上のパラメータは、前記心周期の開始タイミングから前記血流速度のピーク到達タイミングまでの間の立ち上がり時間と、前記血流速度の立ち上がり期間に前記血流速度が第１血流速度になるタイミングから前記ピーク到達タイミングまでの間の第１時間と、前記ピーク到達タイミングから前記ピーク到達タイミング以降の前記血流速度の立ち下がり期間に前記血流速度が前記第１血流速度になるタイミングまでの間の第２時間とを含む、血流解析方法。
  11. １本の網膜動脈又は網膜静脈の血流速度の経時的変化を表す血流情報を取得する取得ステップと、
    前記血流情報から血流速度の変化に対応する１以上のパラメータを抽出する抽出ステップと、
    を含み、
    前記１以上のパラメータは、拡張期における前記血流速度の経時的変化の波形下面積と、前記拡張期の前半における前記血流速度の経時的変化の波形下面積と、前記拡張期の後半における前記血流速度の経時的変化の波形下面積とを含む、血流解析方法。
  12. 前記１以上のパラメータに基づいて、前記血流速度の変化が加齢とともに変化する加齢性変化であるか否かを評価する評価ステップを含む
    ことを特徴とする請求項８～請求項１１のいずれか一項に記載の血流解析方法。
  13. 前記網膜動脈又は前記網膜静脈は、内境界膜からブルッフ膜までの間の層領域の動脈又は静脈である
    ことを特徴とする請求項８～請求項１２のいずれか一項に記載の血流解析方法。
  14. 略同一の計測部位に対して異なるタイミングで光コヒーレンストモグラフィを実行することにより得られた複数の信号に基づいて前記血流情報を生成する生成ステップを含む
    ことを特徴とする請求項８～請求項１３のいずれか一項に記載の血流解析方法。
  15. コンピュータに、請求項８～請求項１４のいずれか一項に記載の血流解析方法の各ステップを実行させることを特徴とするプログラム。
JP2020101391A 2020-06-11 2020-06-11 血流解析装置、血流解析方法、及びプログラム Active JP7479030B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2020101391A JP7479030B2 (ja) 2020-06-11 2020-06-11 血流解析装置、血流解析方法、及びプログラム
US17/319,089 US20210386290A1 (en) 2020-06-11 2021-05-13 Blood flow analysis apparatus, blood flow analysis method, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020101391A JP7479030B2 (ja) 2020-06-11 2020-06-11 血流解析装置、血流解析方法、及びプログラム

Publications (2)

Publication Number Publication Date
JP2021194146A JP2021194146A (ja) 2021-12-27
JP7479030B2 true JP7479030B2 (ja) 2024-05-08

Family

ID=78826375

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020101391A Active JP7479030B2 (ja) 2020-06-11 2020-06-11 血流解析装置、血流解析方法、及びプログラム

Country Status (2)

Country Link
US (1) US20210386290A1 (ja)
JP (1) JP7479030B2 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007132865A1 (ja) 2006-05-16 2007-11-22 Retinal Information Diagnosis Research Institute Inc. 血管老化の検出システム
WO2008069062A1 (ja) 2006-12-01 2008-06-12 Kyushu Tlo Company, Limited 血流速度画像化装置
WO2010131550A1 (ja) 2009-05-13 2010-11-18 国立大学法人九州工業大学 血流画像診断装置
JP2019042264A (ja) 2017-09-05 2019-03-22 株式会社トプコン 眼科撮影装置及び眼科情報処理装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08215150A (ja) * 1995-02-09 1996-08-27 Canon Inc 眼科診断装置
EP3272276B1 (en) * 2012-03-12 2019-05-15 Kabushiki Kaisha Topcon Image displaying apparatus and image displaying method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007132865A1 (ja) 2006-05-16 2007-11-22 Retinal Information Diagnosis Research Institute Inc. 血管老化の検出システム
WO2008069062A1 (ja) 2006-12-01 2008-06-12 Kyushu Tlo Company, Limited 血流速度画像化装置
WO2010131550A1 (ja) 2009-05-13 2010-11-18 国立大学法人九州工業大学 血流画像診断装置
JP2019042264A (ja) 2017-09-05 2019-03-22 株式会社トプコン 眼科撮影装置及び眼科情報処理装置

Also Published As

Publication number Publication date
US20210386290A1 (en) 2021-12-16
JP2021194146A (ja) 2021-12-27

Similar Documents

Publication Publication Date Title
JP6086345B2 (ja) 眼科装置
JP6580448B2 (ja) 眼科撮影装置及び眼科情報処理装置
JP7348374B2 (ja) 眼科情報処理装置、眼科撮影装置、眼科情報処理方法、及びプログラム
WO2020044712A1 (ja) 眼科装置、及びその制御方法
JP7090438B2 (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JP2024133106A (ja) 眼科撮影装置
JP7096116B2 (ja) 血流計測装置
JP6633468B2 (ja) 血流計測装置
JP2019037646A (ja) 眼科装置
JP6736734B2 (ja) 眼科撮影装置及び眼科情報処理装置
JP7260426B2 (ja) 光コヒーレンストモグラフィ装置、その制御方法、光計測方法、プログラム、及び記憶媒体
JP7141279B2 (ja) 眼科情報処理装置、眼科装置、及び眼科情報処理方法
JP2023038280A (ja) 血流計測装置
JP6646021B2 (ja) 眼科画像処理装置
JP7479030B2 (ja) 血流解析装置、血流解析方法、及びプログラム
JP2019054994A (ja) 眼科撮影装置、眼科情報処理装置、プログラム、及び記録媒体
CN111787843B (zh) 血流测量装置
JP6942627B2 (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
JP7221628B2 (ja) 血流計測装置、情報処理装置、情報処理方法、及びプログラム
JP7216514B2 (ja) 血管解析装置
JP6991075B2 (ja) 血流計測装置
WO2023042577A1 (ja) 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム
JP7068366B2 (ja) 血流計測装置
JP2019054993A (ja) 血流計測装置
JP2022127783A (ja) 血流解析装置、眼科装置、血流解析方法、及びプログラム

Legal Events

Date Code Title Description
A80 Written request to apply exceptions to lack of novelty of invention

Free format text: JAPANESE INTERMEDIATE CODE: A80

Effective date: 20200624

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230518

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231117

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231212

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240111

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240415

R150 Certificate of patent or registration of utility model

Ref document number: 7479030

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150