WO2014132470A1 - 視野検査支援装置 - Google Patents
視野検査支援装置 Download PDFInfo
- Publication number
- WO2014132470A1 WO2014132470A1 PCT/JP2013/073426 JP2013073426W WO2014132470A1 WO 2014132470 A1 WO2014132470 A1 WO 2014132470A1 JP 2013073426 W JP2013073426 W JP 2013073426W WO 2014132470 A1 WO2014132470 A1 WO 2014132470A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- inspection
- test
- estimated
- visual field
- result
- Prior art date
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/02—Subjective types, i.e. testing apparatus requiring the active assistance of the patient
- A61B3/024—Subjective types, i.e. testing apparatus requiring the active assistance of the patient for determining the visual field, e.g. perimeter types
-
- 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/0091—Fixation targets for viewing direction
-
- 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/02—Subjective types, i.e. testing apparatus requiring the active assistance of the patient
- A61B3/06—Subjective types, i.e. testing apparatus requiring the active assistance of the patient for testing light sensitivity, e.g. adaptation; for testing colour vision
- A61B3/063—Subjective types, i.e. testing apparatus requiring the active assistance of the patient for testing light sensitivity, e.g. adaptation; for testing colour vision for testing light sensitivity, i.e. adaptation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
Definitions
- control unit 11 uses a mixture distribution based on the parameter as a prior distribution, and performs Bayes estimation (specifically, Bayesian estimation) representing an estimated inspection result vector representing an estimation result of visual sensitivity at each inspection point as a result of visual field inspection.
- Bayes estimation specifically, Bayesian estimation
- Information output processing based on the estimated test result vector obtained by the estimation.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Data Mining & Analysis (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Eye Examination Apparatus (AREA)
- Image Processing (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
【課題】各検査点における現在の視野の状況を推定できる視野検査支援装置を提供する。 【解決手段】予め定めた混合分布を事前分布として、過去に得られた複数のケースにおける視野検査結果の時系列情報を参照して混合分布に係るパラメータを学習する。視野検査の結果である各検査点での視覚感度の推定結果を表す推定検査結果ベクトルを、ベイズ推定して演算し、当該推定により得られた推定検査結果ベクトルに基づく情報出力処理を行う視野検査支援装置である。
Description
本発明は、視覚の視野検査を支援する視野検査支援装置に関する。
緑内障等の疾病診断で行われる視野検査では、想定されている視野範囲に複数の検査点を配置し、各検査点において指定明度の光を点灯させて、被験者に見えるか見えないかを伺う検査方法が知られている。こうした検査方法には例えば、ハンフリー視野検査法などが知られている。
ところが、こうした視野検査法では、多数の検査点の各々において、明度を変化させつつ被験者の反応を伺う必要があり、被験者だけでなく検査担当者の負担も大きいのが現状である。
なお、緑内障の患者の視野の変化を推定する技術がViswanathan AC, Fitzke FW, Hitchings RA. Early Detection of Visual Field Progression in Glaucoma: A Comparison of PROGRESSOR and Statpac2. Br J Ophthalmol
1997;81(12):1037-1042.に開示されている。
1997;81(12):1037-1042.に開示されている。
上述のように、視野検査における負担の軽減が望まれており、具体的には各検査点における現在の視野の状況を推定できる視野検査支援装置が望まれているところである。
本発明は上記実情に鑑みて為されたもので、各検査点における現在の視野の状況を推定できる視野検査支援装置を提供することを、その目的の一つとする。
上記従来例の問題点を解決する本発明は、予め想定された視野範囲に配された複数の検査点の各々における視覚感度を検査する視野検査を支援する装置であって、予め定めた混合分布を事前分布として、過去に得られた複数のケースにおける視野検査結果の時系列情報を参照して学習した、前記混合分布に係るパラメータを取得する取得手段と、前記取得したパラメータに基づく混合分布を事前分布として、視野検査の結果である前記各検査点での視覚感度の推定結果を表す推定検査結果ベクトルを、ベイズ推定して演算する推定演算手段と、当該推定により得られた推定検査結果ベクトルに基づく情報出力処理を行う出力手段と、を含み、前記推定演算手段は、視野検査中に、処理の時点までで得られている検査点が少なくとも一つあれば、当該検査点での視覚感度を取得して前記事前分布を更新し、処理の時点では得られていない検査点での視覚感度の推定結果を表す部分推定検査結果ベクトルを、ベイズ推定し、前記出力手段が、当該部分推定検査結果ベクトルに基づく情報出力処理を行うこととしたものである。
本発明によると、現在の視覚感度の推定結果が得られる。
本発明の実施の形態について図面を参照しながら説明する。本発明の実施の形態に係る視野検査支援装置1は、一般的なコンピュータを用いて実現でき、具体的には図1に例示するように、制御部11,記憶部12,操作部13,表示部14,入出力部15を含んで構成される。
ここで制御部11は、CPU等のプログラム制御デバイスであり、記憶部12に格納されたプログラムに従って動作する。本実施の形態においてこの制御部11は、予め想定された視野範囲に配された複数の検査点の各々における視覚感度を検査する視野検査を支援する現在視野の推定処理と、将来視野の推定処理とを行う。このうち現在視野の推定処理では、予め定めた混合分布を事前分布として、過去に得られた複数のケースにおける視野検査結果の時系列情報を参照してこの混合分布に係るパラメータを学習して記憶しておく。そしてこの制御部11は、当該パラメータに基づく混合分布を事前分布として、視野検査の結果である各検査点での視覚感度の推定結果を表す推定検査結果ベクトルを、ベイズ推定(具体的にはベイズ線形回帰による推定)して演算し、当該推定により得られた推定検査結果ベクトルに基づく情報出力処理を行う。
またこの制御部11は、視野検査中に、処理の時点までで得られている検査点における視覚感度の情報が入出力部15を介して入力されると、当該部分的な視覚感度の情報を得て事前分布を更新し、当該処理の時点では得られていない検査点での視覚感度の推定結果を表す部分推定検査結果ベクトルを、ベイズ線形回帰により推定する。そして制御部11は、当該部分推定検査結果ベクトルに基づく情報出力処理を行う。これらの処理については、後に詳しく述べる。
さらに制御部11は、将来視野の推定処理として、検査点の各点における視覚感度の時間変化の情報を取得し、当該取得した検査点ごとの視覚感度の時間変化についての回帰分析処理を行う。ここで制御部11は、当該回帰分析処理において得られたスロープに係るp値が、予め定めた分析しきい値を下回る検査点の個数をカウントし、検査点の個数についての二項検定を行い、当該検定の結果である二項検定p値を演算する。制御部11は、演算の結果得られた二項検定p値について平均化処理を実行し、平均化した二項検定p値を演算し、当該平均化した二項検定p値が予め定めた検定しきい値を下回るか否かにより、視覚感度変化の有無を判定して当該判定の結果を出力する。この制御部11の処理についても後に述べる。
記憶部12は、メモリデバイスやディスクデバイス等であり、制御部11によって実行されるプログラムを保持する。このプログラムは、コンピュータ可読な記録媒体に格納されて提供され、この記憶部12に複写されたものであってもよい。さらにこの記憶部12は、パラメータ値など、制御部11の処理に必要な情報を保持する。また記憶部12は、制御部11のワークメモリとしても動作する。
操作部13は、キーボードやマウス等を含む。この操作部13は、利用者の指示操作を受け入れて、当該指示操作の内容を制御部11に出力する。表示部14は、ディスプレイやプリンタ等であり、制御部11から入力される指示に従って情報を出力する。
入出力部15は、データの入出力を行うシリアル、またはパラレルインタフェースであり、外部の装置から情報を受け入れて制御部11に出力する。またこの入出力部15は制御部11から入力される指示に従い、情報を外部の装置宛に出力する。
ここで制御部11による現在視野の推定処理について説明する。以下の説明では視野の検査結果を表す検査点の数Dt×検査回数Nmの行列Tm(mは検査の対象となる目を識別する識別子で、m=1,2,…,Mとする)を確率変数として、
とモデル化しておく。ここでwmは視野検査の切片及び進行速度を表すベクトルであり、
Lは精度行列(分散共分散行列の逆行列)、λは精度を決定するスカラ値である。またここで、丸囲みの「×」は、クロネッカー積を表しており、
は、多変量正規分布を表す。さらにxnは、初回の検査からの経過時間を表し、IDtを検査点Dt×Dtの単位行列とする。
Lは精度行列(分散共分散行列の逆行列)、λは精度を決定するスカラ値である。またここで、丸囲みの「×」は、クロネッカー積を表しており、
また、視野検査の切片及び進行速度を表すベクトルwmを、混合正規分布として、
とし、精度行列についても所定の混合分布を用いて記述する。一例として精度行列について混合ウィシャート分布を用いる場合、
とする。
このように混合分布としたのは、視野の欠損が必ずしも一ヶ所から発生するものではなく、複数個所にて発生することがあるという経験的知識に基づく。
制御部11では、現在の視野検査結果を推定する処理に先立ち、事前に、この(*)式に、(1)乃至(6)式を代入した式を用いて、パラメータ(ハイパーパラメータ)としての
を推定しておく。この推定は具体的にはEMアルゴリズムを利用することができる。
なお、平均場近似における最適な事後分布を、それぞれの記号の肩にアスタリスク(*)を付して表すものとし、Nm,Dwをそれぞれ、目mについての検査回数、並びにベクトルμmの次元とする。
また、Mステップ(最大化ステップ)におけるハイパーパラメータの再推定方程式が、
について、
となる。また、νhについては、
の解を新たなνh、つまり
とし、aについては
の解を新たなa、つまり
とする。なお、ここで、ψ(・)は、ディガンマ関数を表す。
制御部11は、一例として上述のEステップ、Mステップを交互に行い、最適なハイパーパラメータを演算する。なお、尤度P(T)の下界は、次のように演算できる。
また制御部11は、操作部13から現在の視野検査の結果を推定するべき旨の指示を受けて、次の推定処理を実行する。以下の推定処理では、(1)過去にN回の視野検査を受けている被験者についての推定処理と、(2)初めて視野検査を受ける被験者についての推定処理とで異なるので、これらに分けて述べる。なお、制御部11は、どちらの場合であるかの指示を受け入れて、次のいずれかの処理を実行するものとする。
[過去にN回の視野検査を受けている被験者の場合]
過去にN回の視野検査を受けている被験者の場合、制御部11は、現在視野検査を受けた場合の各検査点での視覚感度の推定結果を表す推定検査結果ベクトルをtN+1、過去N回分の実際の各検査点での視覚感度の測定結果を表す検査結果ベクトルをそれぞれt1,t2…tNとして、この推定検査結果ベクトルtN+1の完全尤度を
とおいてベイズ推定演算する。
過去にN回の視野検査を受けている被験者の場合、制御部11は、現在視野検査を受けた場合の各検査点での視覚感度の推定結果を表す推定検査結果ベクトルをtN+1、過去N回分の実際の各検査点での視覚感度の測定結果を表す検査結果ベクトルをそれぞれt1,t2…tNとして、この推定検査結果ベクトルtN+1の完全尤度を
具体的にはtn(n=1,2,…N)を得た後の事後分布を平均場近似により、
として推定する。最適な変分事後分布は具体的には次のようになる。すなわち、
ただし
ただし
ただし
ただし
ただし
ただし
なお、
である。また、
については、記憶部12に保持されている、事前に演算したパラメータを用いる。
従って制御部11が求めるべきtN+1の推定結果は、
として求められる。制御部11はこの(15)式で表される多変量正規分布により、現在視野検査を受けた場合の各検査点での視覚感度の推定結果を表す推定検査結果ベクトルtN+1を得る。そして制御部11は、この推定結果ベクトルtN+1に基づく情報出力処理を行う。この情報出力処理については、後に述べる。
さらに制御部11は、実際の視野検査中においても、情報出力内容を更新する。すなわち、視野検査により実際の検査結果の一部(一部の検査点における視覚感度の情報)が入出力部15を介して入力されると、制御部11は、当該処理の時点までで得られている検査点での視覚感度の情報を参照して事前分布を更新する。具体的に制御部11は、検査結果ベクトルtN+1を、この処理の時点で得られていない検査点での視覚感度の情報(ベク
トル)taと、処理の時点で得られている検査点での視覚感度の情報(ベクトル)tbとに分割する。この分割された要素を便宜的に(検査点の測定順が予め定まっていれば、ベクトルtの要素を予めこの測定順にしておけばよい)、
と表す。また、対応するμ,Λをそれぞれ
とする。そして、制御部11は、未測定の検査点についての推定結果taを、部分推定検査結果ベクトルとして、
を用いて演算する。これにより複数の未測定の検査点の各々を注目検査点とした、注目検査点における推定の検査結果の情報を得る。
トル)taと、処理の時点で得られている検査点での視覚感度の情報(ベクトル)tbとに分割する。この分割された要素を便宜的に(検査点の測定順が予め定まっていれば、ベクトルtの要素を予めこの測定順にしておけばよい)、
制御部11は、これにより得られた推定結果ta(P(ta|tb))に基づく情報出力処理を行う。なお、情報出力処理については、後に述べる。
[初めて視野検査を受ける被験者の場合]
初めて視野検査を受ける被験者の場合は、過去の情報が存在しないため、制御部11は、初めて視野検査を受ける被験者についての推定を行うべき旨の指示を受けると、各検査点での視覚感度の情報の推定結果たる事後分布を、混合正規分布を用いて近似的に推定する。
初めて視野検査を受ける被験者の場合は、過去の情報が存在しないため、制御部11は、初めて視野検査を受ける被験者についての推定を行うべき旨の指示を受けると、各検査点での視覚感度の情報の推定結果たる事後分布を、混合正規分布を用いて近似的に推定する。
具体的に、時刻0における完全尤度を、N回の視野検査を受けている被験者に倣って
とおき、これを
として近似する。制御部11は、カルバック・ライブラー・ダイバージェンス
を最小化する混合正規分布を求める。
制御部11は、予め定めた混合正規分布関数を用い、そのパラメータの初期値を適宜(例えばランダムに)設定する。そして制御部11は当該混合正規分布関数のカルバック・ライブラー・ダイバージェンスを演算し、これを最小化するようパラメータを更新して、再度カルバック・ライブラー・ダイバージェンスを演算する、という処理を、カルバック・ライブラー・ダイバージェンスが収束するまで繰返し行う。この演算の具体的態様は広く知られている方法が採用できる。
制御部11は、この処理により、最適な分布となったt0について
を得る。さらに制御部11はこれをw,zについて周辺化してt0を得る。具体的に制御部11は、
を演算してt0を得る。なお、制御部11は、これらの演算において記憶部12に格納されているパラメータの値を用いる。またこのときのカルバック・ライブラー・ダイバージェンスは、
と表される。なおここで第3項と第10項とがキャンセルしあい、第1項と第8項との和が
であり、第2項と第9項との和が
と表せる。また、第4項から第7項、第11項、第12項の期待値はそれぞれ
と表すことができる。制御部11は、ここで得たt0に基づいて情報出力処理を行う。なお、情報出力処理については、後に述べる。
さらに、制御部11は、初めて検査を受ける被験者についての実際の視野検査中においても、情報出力内容を更新する。すなわち、視野検査により実際の検査結果の一部(一部の検査点における視覚感度の情報)が入出力部15を介して入力されると、制御部11は、当該処理の時点までで得られている検査点での視覚感度の情報を参照して事前分布を更新する。具体的に制御部11は、検査結果ベクトルtを、この処理の時点で得られていない検査点での視覚感度の情報(ベクトル)taと、処理の時点で得られている検査点での視覚感度の情報(ベクトル)tbとに分割する。この分割された要素についても便宜的に(検査点の測定順が予め定まっていれば、ベクトルtの要素を予めこの測定順にしておけばよい)、
と表す。また、対応するμ,Λをそれぞれ
とする。そして、制御部11は、一部の検査点についての推定結果taを、部分推定検査結果ベクトルとして、
を用いて演算する。ここで
である。これにより複数の未測定の検査点の各々を注目検査点とした、注目検査点における推定の検査結果の情報を得る。
制御部11は、これにより得られた推定結果ta(P(ta|tb))に基づく情報出力処理を行う。なお、情報出力処理については、次に述べる。
[情報出力処理]
制御部11は、以上のようにして、未測定の検査点についての視覚感度の推定結果を含む、推定検査結果ベクトル、または部分推定検査結果ベクトルを得て情報出力処理を行う。この情報出力処理では、制御部11は、例えば各検査点を特定する情報と、推定検査結果ベクトル、または部分推定検査結果ベクトルに含まれる、対応する検査点(未測定の検査点)について、その少なくとも一部についての視覚感度の推定結果を一覧としてそのまま表示してもよい。
制御部11は、以上のようにして、未測定の検査点についての視覚感度の推定結果を含む、推定検査結果ベクトル、または部分推定検査結果ベクトルを得て情報出力処理を行う。この情報出力処理では、制御部11は、例えば各検査点を特定する情報と、推定検査結果ベクトル、または部分推定検査結果ベクトルに含まれる、対応する検査点(未測定の検査点)について、その少なくとも一部についての視覚感度の推定結果を一覧としてそのまま表示してもよい。
これにより検査担当者は、各検査点で、被験者が所定の確率で視認すると推定される明度の情報を得ることができ、当該確からしい情報が表す明度を参考に、各検査点で最初に、被験者に対して提示する明度を選択することで、ある被験者に対して明らかに視認可能な明度を試したりする機会を低減でき、検査点ごとの繰返し回数を低減して、検査担当者、被験者の両者に対する検査の負担を軽減できる。
また制御部11は、推定検査結果ベクトルまたは部分推定検査結果ベクトルに含まれる検査結果から、未測定の検査点の一つを注目検査点として次の処理を行ってもよい。すなわち制御部11は、当該注目検査点近傍にある検査点における視覚感度の推定結果または測定結果(実際の検査の結果)との差を演算し、当該演算により得られた差が予め定めた空間変化しきい値を超える場合に、当該注目検査点の近傍空間を特定する情報を出力する。
一例として、図2に例示するように、未測定の検査点Pi,jと、この検査点Pi,jに隣接して予め設定されている検査点Pi-1,j、Pi,j-1,Pi+1,j,Pi,j+1(いわゆる四近傍の検査点)につき、それぞれ推定された視覚感度の値がV[Pi,j],V[Pi-1,j],V[Pi,j-1],V[Pi+1,j],V[Pi,j+1]であるとき、制御部11は、これらの差Δi-1,j=|Pi-1,j-Pi,j|,Δi,j-1=|V[Pi,j-1]-Pi,j|,Δi+1,j=|V[Pi+1,j]-Pi,j|,Δi,j+1=|V[Pi,j+1]-Pi,j|(ここで|*|は、*の絶対値を意
味する)を演算する。
味する)を演算する。
そして制御部11は、これらの差のいずれかが予め定めたしきい値(空間変化しきい値と呼ぶ)を超える場合に、注目検査点Pi,jの近傍空間を特定する情報を表示部14に出力する。このとき制御部11は、併せて、当該特定した空間に、さらに多くの検査点を設けるよう示唆する表示を表示部14に行ってもよい。ここで注目検査点Pi,jの近傍空間は、具体的には注目検査点Pi,jの位置そのものでよい。また注目検査点Pi,jの近傍空間は、注目検査点Pi,jから予め定めた距離rの範囲内などであってもよい。
検査担当者は、この表示を参照して、必要に応じて当該表示された範囲乃至位置の周囲にさらに多くの検査点を配して、当該検査点における測定を行ってもよい。これにより、より精度の高い検査を行うことが可能となる。
[将来推定]
さらに本実施の形態の視野検査支援装置1は、次のようにして将来の視野検査結果についての予測を行い、その予測結果を出力する処理を行ってもよい。この場合制御部11は、指定された被験者について、当該被験者が過去に受けた検査の結果を受け入れる。この検査の結果には、過去の複数の時点での各検査点における視覚感度の情報が含まれる。
さらに本実施の形態の視野検査支援装置1は、次のようにして将来の視野検査結果についての予測を行い、その予測結果を出力する処理を行ってもよい。この場合制御部11は、指定された被験者について、当該被験者が過去に受けた検査の結果を受け入れる。この検査の結果には、過去の複数の時点での各検査点における視覚感度の情報が含まれる。
そして制御部11は、検査点ごとに、上記取得した情報が表す視覚感度の時間変化についての回帰分析処理を行って、その切片とスロープとを表す値、並びにそれらに係るp値を得る。制御部11は、このうちスロープに係るp値が、予め定めた分析しきい値を下回る(検査結果についての時間変化が進行している)検査点の個数をカウントする。また制御部11は、検査点の個数についての片側二項検定を行い、当該検定の結果である二項検定p値(検査点の個数が有意でない確率)を演算する。なお、進行している検査点の個数については、事前の割合を得ておく。
制御部11は、分析しきい値を予め定めた互いに異なる複数の値に設定しつつ、複数の分析しきい値において上記の処理を繰返し行う。ここで分析しきい値の設定は、例えば0.025,0.5,0.075,0.1の4つとする。制御部11は、それぞれの分析しきい値に対応して得られた二項検定p値について平均化した二項検定p値を演算する(具体的には各二項検定p値の中央値を採用すればよい)。制御部11は、当該平均化した二項検定p値が予め定めた検定しきい値を下回る場合に、視覚感度の時間変化が有意である(検査結果(視覚感度)についての時間変化が進行している)と判断する。そして制御部11は、この視覚感度の時間変化の有無についての判定結果を表示部14に出力する。
また制御部11は、視覚感度の時間変化があると判定した場合は、所定の分析しきい値を下回るp値(回帰分析結果のスロープに関するp値)となっている検査点の、視野空間内の位置を表示する。これにより、視覚感度の時間変化があると推定される検査点の視野内での配置が提示されることとなる。そして制御部11は、当該視覚感度の時間変化があると判断された各検査点についての回帰分析結果から、将来の指定時点における視野の推定結果を得て、当該推定結果を出力する。これは視覚感度の時間変化があると判断された各検査点についての回帰分析結果と、指定時点までの時間とを回帰式に代入し、上記各検査点における視覚感度の情報を得て表示することで行われる。
これにより視覚感度について進行がある部分と、当該部分の将来の視覚感度が提示されることとなり、検査担当者及び被験者は、将来的にどのような視野検査結果が得られることとなるかの推定結果を得られることとなる。この場合、例えば医療担当者が、被験者に対して、当該被験者の生活の状況に鑑みて、どのような状況になるかを具体的に指摘可能となり、治療等の指針決定に資することが可能となる。
[動作]
本実施の形態の視野検査支援装置1は、以上の構成を備えてなり、次のように動作する。すなわちこの視野検査支援装置1は、予め定められた混合分布を事前分布としたベイズ推定により、事後分布としての現在または将来の視野検査の結果を推定する。この際、当該混合分布に関するパラメータを予め、既知の情報である、過去に測定された複数の被験者についての、複数の機会における測定結果(時間をおいて測定された時系列の情報)を用いた学習処理により得ておく。
本実施の形態の視野検査支援装置1は、以上の構成を備えてなり、次のように動作する。すなわちこの視野検査支援装置1は、予め定められた混合分布を事前分布としたベイズ推定により、事後分布としての現在または将来の視野検査の結果を推定する。この際、当該混合分布に関するパラメータを予め、既知の情報である、過去に測定された複数の被験者についての、複数の機会における測定結果(時間をおいて測定された時系列の情報)を用いた学習処理により得ておく。
検査の際には、この視野検査支援装置1は、過去に検査を受けている被験者に対しては、当該過去の検査の結果に係る情報と、先に求めたパラメータを用いた混合分布とを用いたベイズ推定を行って、現在の検査における視野検査の結果(各検査点で所定の確率で被験者が視認可能となると推定される明度の情報)を推定する。そして視野検査支援装置1は当該推定の結果を検査担当者に提示する。検査担当者は、当該明度の情報を参考に、各検査点で、最初に被験者に提示するべき明度を決定する。
また検査中には、検査担当者が、検査により得られた各検査点での明度の情報をこの視野検査支援装置1に入力する(あるいは検査に用いる装置から人為的操作なしに当該情報が入力されるよう設定しておく)。視野検査支援装置1は、検査中に得られた検査点での明度の情報を受け入れると、当該情報で事前分布を更新し、ベイズ推定により、未測定の検査点における視野検査の結果の推定を行う。そして視野検査支援装置1は、当該新たに推定した結果を検査担当者に提示する。検査担当者は、当該新たに推定された結果を参照して、未測定の検査点において最初に被験者に提示するべき明度を決定する。
なお、初めて検査を受ける被験者に対しては、予め定められた混合正規分布により推定を行うことで、同様の処理を実現する。これにより検査担当者は、試行錯誤的に検査結果を得る場合に比べ、所定の確率で確からしい検査結果に基づき、初期条件を設定でき、従って比較的短時間で検査結果を得られることが期待できる。
[変形例]
ここまでの説明における視野検査支援装置1において、ある検査点(注目検査点)近傍にある検査点における視覚感度の推定結果または測定結果(実際の検査の結果)との差を演算し、当該演算により得られた差が予め定めた空間変化しきい値を超える場合に、当該注目検査点の近傍の空間を特定する情報を出力することとしていたが、この空間変化しきい値は、視野空間全体で同じ値である必要はない。
ここまでの説明における視野検査支援装置1において、ある検査点(注目検査点)近傍にある検査点における視覚感度の推定結果または測定結果(実際の検査の結果)との差を演算し、当該演算により得られた差が予め定めた空間変化しきい値を超える場合に、当該注目検査点の近傍の空間を特定する情報を出力することとしていたが、この空間変化しきい値は、視野空間全体で同じ値である必要はない。
例えば、視野空間において上半分にある検査点と下半分にある検査点とでは、各検査点における検査結果の値に相関が小さいことが知られている。そこで視野空間の上下半分をまたぐ二点については、当該二点間の差分Δと比較するべき空間変化しきい値を、他の部分における空間変化しきい値より大きくしておいてもよい。
また、ここまでに説明した視野検査支援装置1では、ベイズ推定に用いるパラメータは一つのセットとなっていたが、人種ごと、性別ごと、あるいは住地域ごと、などといったように複数の層別にわけて、それぞれの層別に属する被験者の検査結果を用いて学習させることで、層別ごとのパラメータを得てもよい。
この場合、視野検査支援装置1は、推定を行おうとする被験者と同じ層別に属する被験者の情報に基づいて学習して得たパラメータを選択して用いて、当該推定を行おうとする被験者についての検査結果を推定する。
Claims (5)
- 予め想定された視野範囲に配された複数の検査点の各々における視覚感度を検査する視野検査を支援する装置であって、
予め定めた混合分布を事前分布として、過去に得られた複数のケースにおける視野検査結果の時系列情報を参照して学習した、前記混合分布に係るパラメータを取得する取得手段と、
前記取得したパラメータに基づく混合分布を事前分布として、視野検査の結果である前記各検査点での視覚感度の推定結果を表す推定検査結果ベクトルを、ベイズ推定して演算する推定演算手段と、
当該推定により得られた推定検査結果ベクトルに基づく情報出力処理を行う出力手段と、
を含み、
前記推定演算手段は、視野検査中に、処理の時点までで得られている検査点が少なくとも一つあれば、当該検査点での視覚感度を取得して前記事前分布を更新し、処理の時点では得られていない検査点での視覚感度の推定結果を表す部分推定検査結果ベクトルを、ベイズ推定し、前記出力手段が、当該部分推定検査結果ベクトルに基づく情報出力処理を行う視野検査支援装置。 - 請求項1記載の視野検査支援装置であって、
前記出力手段は、前記推定演算手段が推定した、推定検査結果ベクトルに基づく情報出力処理、または、前記処理の時点で得られていない検査点での視覚感度の推定結果を表す部分推定検査結果ベクトルに基づく情報出力処理として、
推定検査結果ベクトルまたは部分推定検査結果ベクトルに含まれる検査結果から、未測定の検査点の少なくとも一部を注目検査点として、各注目検査点で被験者が所定の確率で視認すると推定される明度に係る情報を演算し、各注目検査点を特定する情報と、当該特定された注目検査点について推定された情報とを出力する。 - 請求項1記載の視野検査支援装置であって、
前記出力手段は、前記推定演算手段が推定した、推定検査結果ベクトルに基づく情報出力処理、または、前記処理の時点で得られていない検査点での視覚感度の推定結果を表す部分推定検査結果ベクトルに基づく情報出力処理として、
推定検査結果ベクトルまたは部分推定検査結果ベクトルに含まれる検査結果から、未測定の検査点の一つを注目検査点として、当該注目検査点近傍にある検査点における視覚感度の推定結果または測定結果との差を演算し、当該演算により得られた差が予め定めた空間変化しきい値を超える場合に、当該注目検査点の近傍空間を特定する情報を出力する。 - 請求項1に記載の視野検査支援装置であって、
前記検査点の各点における視覚感度の時間変化の情報を取得する手段と、
当該取得した検査点ごとの視覚感度の時間変化についての回帰分析処理を行う手段と、
前記回帰分析処理において得られたスロープに係るp値が、予め定めた分析しきい値を下回る検査点の個数をカウントする手段と、
前記分析しきい値を相異なる複数の値のそれぞれに設定して前記回帰分析処理及びカウントの処理を実行させ、各分析しきい値に対応する検査点の個数についての二項検定を行い、各分析しきい値に対応する検定の結果である二項検定p値を演算する手段と、
前記演算の結果得られた、各分析しきい値に対応する二項検定p値について平均化処理を実行し、平均化した二項検定p値を演算し、当該平均化した二項検定p値が予め定めた検定しきい値を下回るか否かにより、視覚感度変化の有無を判定する判定手段と、
当該判定の結果を出力する出力手段と、をさらに含む。 - 請求項4記載の視野検査支援装置であって、
視覚感度変化ありと判定される検査点の空間配置と、前記回帰分析処理により得られたスロープとに基づいて将来の指定時点における視野の推定結果を出力する手段をさらに含む。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/771,150 US11291360B2 (en) | 2013-02-28 | 2013-08-30 | Visual field test assistance apparatus |
EP13876498.0A EP2984983B1 (en) | 2013-02-28 | 2013-08-30 | Visual-field-test assistance device |
CN201380074073.9A CN105188515B (zh) | 2013-02-28 | 2013-08-30 | 视场检查支援装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013-039796 | 2013-02-28 | ||
JP2013039796A JP6053574B2 (ja) | 2013-02-28 | 2013-02-28 | 視野検査支援装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014132470A1 true WO2014132470A1 (ja) | 2014-09-04 |
Family
ID=51427764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/073426 WO2014132470A1 (ja) | 2013-02-28 | 2013-08-30 | 視野検査支援装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11291360B2 (ja) |
EP (1) | EP2984983B1 (ja) |
JP (1) | JP6053574B2 (ja) |
CN (1) | CN105188515B (ja) |
WO (1) | WO2014132470A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016190200A1 (ja) * | 2015-05-28 | 2016-12-01 | 国立大学法人 東京大学 | 視野検査装置 |
JP7459623B2 (ja) | 2020-03-31 | 2024-04-02 | 株式会社ニコン | 視野検査方法、視野検査装置、およびプログラム |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201613923D0 (en) | 2016-08-15 | 2016-09-28 | Ucl Business Plc | Method and apparatus for performing a visual field test |
EP3704126A4 (en) | 2017-11-03 | 2020-12-09 | Université de Montréal | COMPOUNDS AND THEIR USE IN THE EXPANSION OF STEM CELLS AND / OR PROGENITOR CELLS |
WO2019169166A1 (en) * | 2018-03-01 | 2019-09-06 | The Schepens Eye Research Institute, Inc. | Visual field progression |
JP7082889B2 (ja) * | 2018-03-19 | 2022-06-09 | 株式会社クリュートメディカルシステムズ | 視野検査装置、その制御方法及び視野検査プログラム |
CN209404742U (zh) * | 2018-04-03 | 2019-09-20 | 林臣 | 一种视野计 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012506761A (ja) * | 2008-10-28 | 2012-03-22 | オレゴン・ヘルス・アンド・サイエンス・ユニバーシティ | 視野モニタリングのための方法及び装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4818091A (en) | 1985-12-31 | 1989-04-04 | University Of Southern California | System and method of detecting visual field defects |
US6033076A (en) | 1996-07-31 | 2000-03-07 | Virtual-Eye.Com, Inc. | Visual field testing via telemedicine |
JP4171995B2 (ja) | 1997-10-16 | 2008-10-29 | 潤 福原 | フリッカー感度分布測定装置及びフリッカー感度分布測定プログラムを記録したコンピューター読み取り可能な記録媒体 |
CA2706385C (en) * | 2007-07-17 | 2015-11-17 | Heidelberg Engineering Gmbh | Method and apparatus for evaluating the field of vision |
US7938538B2 (en) * | 2008-05-09 | 2011-05-10 | University Of Southern California | Methods and devices for rapid measurement of visual sensitivity |
JP4187227B1 (ja) | 2008-07-28 | 2008-11-26 | 知博 蔦 | 静的視野スキャン装置 |
JP4888579B2 (ja) | 2010-04-21 | 2012-02-29 | パナソニック電工株式会社 | 視機能検査装置 |
CN101966075A (zh) | 2010-11-16 | 2011-02-09 | 天津滨海华医光电技术有限公司 | 便携式青光眼快速筛选分析仪 |
CN102499623B (zh) | 2011-09-28 | 2013-08-21 | 北京工业大学 | 一种自校正面阵源视野仪 |
WO2013170091A1 (en) * | 2012-05-09 | 2013-11-14 | The Schepens Eye Research Institute, Inc. | Rapid measurement of visual sensitivity |
CN102793530A (zh) | 2012-08-31 | 2012-11-28 | 天津滨海华医光电技术有限公司 | 便携式静态多点同步阈值筛查兼单点阈值检测的视野计 |
-
2013
- 2013-02-28 JP JP2013039796A patent/JP6053574B2/ja active Active
- 2013-08-30 WO PCT/JP2013/073426 patent/WO2014132470A1/ja active Application Filing
- 2013-08-30 CN CN201380074073.9A patent/CN105188515B/zh active Active
- 2013-08-30 US US14/771,150 patent/US11291360B2/en active Active
- 2013-08-30 EP EP13876498.0A patent/EP2984983B1/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012506761A (ja) * | 2008-10-28 | 2012-03-22 | オレゴン・ヘルス・アンド・サイエンス・ユニバーシティ | 視野モニタリングのための方法及び装置 |
Non-Patent Citations (2)
Title |
---|
HIROSHI MURATA ET AL.: "Gankai no Tameno Sentan Iryo Dai 148 Kai Shiya no Atarashii Tokeigakuteki Kaiseki Hoho", JOURNAL OF THE EYE, vol. 30, no. 4, 30 April 2013 (2013-04-30), pages 509 - 511, XP008180485 * |
VISWANATHAN AC; FITZKE FW; HITCHINGS RA: "Early Detection of Visual Field Progression in Glaucoma: A Comparison of PROGRESSOR and Statpac2", BR J OPHTHALMOL, vol. 81, no. 12, 1997, pages 1037 - 1042 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016190200A1 (ja) * | 2015-05-28 | 2016-12-01 | 国立大学法人 東京大学 | 視野検査装置 |
JP2016220832A (ja) * | 2015-05-28 | 2016-12-28 | 国立大学法人 東京大学 | 視野検査装置 |
US10631722B2 (en) | 2015-05-28 | 2020-04-28 | The University Of Tokyo | Visual field testing device |
JP7459623B2 (ja) | 2020-03-31 | 2024-04-02 | 株式会社ニコン | 視野検査方法、視野検査装置、およびプログラム |
Also Published As
Publication number | Publication date |
---|---|
US20160015263A1 (en) | 2016-01-21 |
US11291360B2 (en) | 2022-04-05 |
EP2984983A1 (en) | 2016-02-17 |
JP2014166250A (ja) | 2014-09-11 |
CN105188515A (zh) | 2015-12-23 |
EP2984983B1 (en) | 2020-09-30 |
EP2984983A4 (en) | 2016-11-02 |
JP6053574B2 (ja) | 2016-12-27 |
CN105188515B (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6053574B2 (ja) | 視野検査支援装置 | |
JP5847735B2 (ja) | 視覚機能検査方法、視覚機能検査システム及びソフトウェア | |
Elmahdy | A new approach for Weibull modeling for reliability life data analysis | |
US10133703B2 (en) | Anomaly detection method, program, and system | |
Norton et al. | Suboptimal criterion learning in static and dynamic environments | |
CN108983744B (zh) | 异常诊断设备和异常诊断方法 | |
Li et al. | A nonlinear Wiener process degradation model with autoregressive errors | |
JP2018025932A (ja) | センサと機械学習部を備えた作業管理システム | |
JP6702297B2 (ja) | プロセスの異常状態診断方法および異常状態診断装置 | |
Zhou et al. | EWMA smoothing and Bayesian spatial modeling for health surveillance | |
Lasko | Nonstationary gaussian process regression for evaluating repeated clinical laboratory tests | |
Panteghini et al. | Effects of the strain-hardening law in the numerical simulation of wire drawing processes | |
Hassani et al. | Model validation and error estimation in multi-block partial least squares regression | |
Chen et al. | Voxel-based Bayesian lesion-symptom mapping | |
JP2009230209A (ja) | ソフトセンサの性能評価装置 | |
Zhang et al. | On the time‐varying predictive performance of longitudinal biomarkers: Measure and estimation | |
Lu | Bayesian semiparametric mixed-effects joint models for analysis of longitudinal-competing risks data with skew distribution | |
Taylor et al. | Calibrating longitudinal models to cross-sectional data: the effect of temporal changes in health practices | |
Levinson | Clinical validation of biomarkers for predicting risk | |
Kangan et al. | Prevalence and associated factors of hypertension among primary school children: A cross-sectional study in Kuching, Sarawak | |
WO2023054268A1 (ja) | 認知機能評価方法、認知機能検査方法、システムおよびコンピュータプログラム | |
Paranjape | Advancing Nitinol Implant Design and Simulation Through Data-Driven Methodologies | |
Pfuhl et al. | Identifying the computational parameters gone awry in psychosis | |
Jin et al. | Unidimensionality, agreement and concordance probability | |
Yang | Joint models for longitudinal and survival data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201380074073.9 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13876498 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14771150 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013876498 Country of ref document: EP |