JPH05137696A - 眼底画像処理方法 - Google Patents
眼底画像処理方法Info
- Publication number
- JPH05137696A JPH05137696A JP3331404A JP33140491A JPH05137696A JP H05137696 A JPH05137696 A JP H05137696A JP 3331404 A JP3331404 A JP 3331404A JP 33140491 A JP33140491 A JP 33140491A JP H05137696 A JPH05137696 A JP H05137696A
- Authority
- JP
- Japan
- Prior art keywords
- image
- time
- image processing
- fundus
- series
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Eye Examination Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【目的】 蛍光眼底画像解析の視覚的理解度を向上させ
る 【構成】 眼底カメラ1から入力された時系列眼底画像
をコンピュータ3に導き時系列眼底画像の各部分におけ
る輝度変化を理論式に対応させて、理論式から導かれる
パラメータ等に応じた画像を形成しディスプレイ5に出
力する。
る 【構成】 眼底カメラ1から入力された時系列眼底画像
をコンピュータ3に導き時系列眼底画像の各部分におけ
る輝度変化を理論式に対応させて、理論式から導かれる
パラメータ等に応じた画像を形成しディスプレイ5に出
力する。
Description
【0001】
【産業上の利用分野】本発明は、時系列蛍光眼底画像を
用いて網膜循環機能を解析し、診断の支援をする眼底画
像処理方法に関するものである。
用いて網膜循環機能を解析し、診断の支援をする眼底画
像処理方法に関するものである。
【0002】
【従来の技術】従来における画像処理による網膜循環機
能解析としては、「ビデオ蛍光眼底造形と画像解析によ
る網膜循環時間の測定」(日眼会誌91巻10号)に記
載の解析法が知られている。この解析法では眼底カメラ
等から得た網膜画像を処理し、画像上に投影されている
動脈、静脈の上に複数個の測定用画素を指定して、時系
列画像における測定用画素の輝度平均値の時間変化に基
づいて血流の網膜平均循環時間を求めている。また、局
所脳血流解析等で用いられる核医学、CT等の画像処理
では、画像を小矩形に分割し時系列画像における小矩形
の平均画素値の時間変化をガンマ関数等の理論式に対応
させ、算出したパラメータからその小矩形の画素値を決
めてイメージ画像を作成するファンクショナルイメージ
法があり、時系列画像の解析法としては有効とされてい
る。
能解析としては、「ビデオ蛍光眼底造形と画像解析によ
る網膜循環時間の測定」(日眼会誌91巻10号)に記
載の解析法が知られている。この解析法では眼底カメラ
等から得た網膜画像を処理し、画像上に投影されている
動脈、静脈の上に複数個の測定用画素を指定して、時系
列画像における測定用画素の輝度平均値の時間変化に基
づいて血流の網膜平均循環時間を求めている。また、局
所脳血流解析等で用いられる核医学、CT等の画像処理
では、画像を小矩形に分割し時系列画像における小矩形
の平均画素値の時間変化をガンマ関数等の理論式に対応
させ、算出したパラメータからその小矩形の画素値を決
めてイメージ画像を作成するファンクショナルイメージ
法があり、時系列画像の解析法としては有効とされてい
る。
【0003】
【発明が解決しようとする課題】しかしながら上述の解
析法では、或る画素値の時間変化は単純に時間系列画像
の同一画素の画素値を用いるため、患者が動くなどして
位置ずれがあった場合には精度が悪化する。また、網膜
平均循環時間だけでは網膜全体の循環機能の詳細な解析
にはデータが不十分であり、満足のゆく診断が下せない
ことが多い。
析法では、或る画素値の時間変化は単純に時間系列画像
の同一画素の画素値を用いるため、患者が動くなどして
位置ずれがあった場合には精度が悪化する。また、網膜
平均循環時間だけでは網膜全体の循環機能の詳細な解析
にはデータが不十分であり、満足のゆく診断が下せない
ことが多い。
【0004】本発明の目的は、複数枚の時系列の蛍光眼
底画像に正確に対応させ、今まで核医学、CT等におい
てのみ用いられていたファンクショナルイメージが得ら
れる眼底画像処理方法を提供することにある。
底画像に正確に対応させ、今まで核医学、CT等におい
てのみ用いられていたファンクショナルイメージが得ら
れる眼底画像処理方法を提供することにある。
【0005】
【課題を解決するための手段】上述の目的を達成するた
めの本発明に係る眼底画像処理方法は、時系列の蛍光眼
底画像の各画像を複数の領域に分割して、前記領域中で
血管部に対応する部位の画素値の平均値を求め、該平均
値を基に前記領域の画素値を算出する第1の工程と、前
記領域ごとに前記画素値の時系列変化を理論式に対応さ
せてファンクショナルイメージを作成する第2の工程を
有することを特徴とする。
めの本発明に係る眼底画像処理方法は、時系列の蛍光眼
底画像の各画像を複数の領域に分割して、前記領域中で
血管部に対応する部位の画素値の平均値を求め、該平均
値を基に前記領域の画素値を算出する第1の工程と、前
記領域ごとに前記画素値の時系列変化を理論式に対応さ
せてファンクショナルイメージを作成する第2の工程を
有することを特徴とする。
【0006】
【作用】上述の構成を有する眼底画像処理方法は、時系
列の蛍光眼底画像が複数の領域に分割され、それぞれの
領域に対して求めるべき画素値が第1の工程によって決
定される。各時系列画像間で対応する決定された画素値
は、第2の工程においてファンクショナルイメージ画像
となって表示される。
列の蛍光眼底画像が複数の領域に分割され、それぞれの
領域に対して求めるべき画素値が第1の工程によって決
定される。各時系列画像間で対応する決定された画素値
は、第2の工程においてファンクショナルイメージ画像
となって表示される。
【0007】
【実施例】本発明を図示の実施例に基づいて詳細に説明
する。図1は実施例のシステム構成図である。眼底画像
処理は眼底カメラ1、スライドスキャナ2、コンピュー
タ3、メモリ4、ディスプレイ5、プリンタ6、入力デ
バイス7によりなされ、全ての装置はコンピュータ3に
接続され制御されるようになっている。
する。図1は実施例のシステム構成図である。眼底画像
処理は眼底カメラ1、スライドスキャナ2、コンピュー
タ3、メモリ4、ディスプレイ5、プリンタ6、入力デ
バイス7によりなされ、全ての装置はコンピュータ3に
接続され制御されるようになっている。
【0008】患者の眼底は眼底カメラ1に設けられたC
CDによって蛍光血管撮像装置に連動して、図2に示す
ように時系列的な複数枚の画像となって撮像される。撮
像された画像は例えば最も明るい乳頭部が抽出され、次
いで乳頭を中心に放射状に広がる血管部が追跡抽出さ
れ、輝度値による2値化処理等によって各画素に対応し
て血管部、非血管部の2値マスク画像がメモリ4上に形
成される。
CDによって蛍光血管撮像装置に連動して、図2に示す
ように時系列的な複数枚の画像となって撮像される。撮
像された画像は例えば最も明るい乳頭部が抽出され、次
いで乳頭を中心に放射状に広がる血管部が追跡抽出さ
れ、輝度値による2値化処理等によって各画素に対応し
て血管部、非血管部の2値マスク画像がメモリ4上に形
成される。
【0009】各時系列画像は互いに共通する部分がコン
ピュータによって自動的に、或いはマウス等のポインテ
ィングデバイスを用いて人為的にポイント指定され、時
系列の画像の適宜な大きさのウインド領域が指定され
る。時系列画像は撮像光学系の収差や被検眼の屈折の作
用によって歪んでおり、正確に一致しない。真の対応を
求めるためにアフィン変換等の幾何変換を用いて画像の
変形が行われ、画像情報として真の対応がメモリ4上に
保存される。
ピュータによって自動的に、或いはマウス等のポインテ
ィングデバイスを用いて人為的にポイント指定され、時
系列の画像の適宜な大きさのウインド領域が指定され
る。時系列画像は撮像光学系の収差や被検眼の屈折の作
用によって歪んでおり、正確に一致しない。真の対応を
求めるためにアフィン変換等の幾何変換を用いて画像の
変形が行われ、画像情報として真の対応がメモリ4上に
保存される。
【0010】時系列画像は図3に示すように小矩形に分
割され、それぞれの小矩形はファンクショナルイメー
ジ、つまり血液の流れる様子を示す画像を作成するとき
の画素単位とされる。
割され、それぞれの小矩形はファンクショナルイメー
ジ、つまり血液の流れる様子を示す画像を作成するとき
の画素単位とされる。
【0011】図4は小矩形画像単位の模式図である。時
系列画像はnピクセル×mピクセルの小矩形に分割さ
れ、メモリ4上に記憶された二値画像をこの小矩形に対
応させて血管部分の画素値から平均画像単位値を算出す
る。図4では、血管部分がP02、P11、P22、P33であ
るから、(P02+P11+P22+P33)/4をこの小矩形
の画像単位値と決める。非血管部分の画素は平均画像単
位値を決める際に用いられないからコントラストの高い
画像が得られる。
系列画像はnピクセル×mピクセルの小矩形に分割さ
れ、メモリ4上に記憶された二値画像をこの小矩形に対
応させて血管部分の画素値から平均画像単位値を算出す
る。図4では、血管部分がP02、P11、P22、P33であ
るから、(P02+P11+P22+P33)/4をこの小矩形
の画像単位値と決める。非血管部分の画素は平均画像単
位値を決める際に用いられないからコントラストの高い
画像が得られる。
【0012】各小矩形の時系列画像から得られた図5に
示す画像単位値は経時変化を追跡するために、図6に示
すように色素希釈曲線等の理論式に合致するようにパラ
メータが決定される。
示す画像単位値は経時変化を追跡するために、図6に示
すように色素希釈曲線等の理論式に合致するようにパラ
メータが決定される。
【0013】色素希釈曲線は色素の希釈の度合いの時間
変化を表す理論式であり、濃度Iは背景濃淡度k、ピー
ク時の濃淡度Ip、曲線の傾きα、立ち上がり時間to 、
ピークタイムtp として、 I=k+Ip ・exp [ −α{ log( t−to )/(tp −to )}2 ] という変化を示す。各画像単位値からこの式に合致する
ように、各パラメータk、Ip、α、tp 、to が最小2
乗法を用いて決定される。
変化を表す理論式であり、濃度Iは背景濃淡度k、ピー
ク時の濃淡度Ip、曲線の傾きα、立ち上がり時間to 、
ピークタイムtp として、 I=k+Ip ・exp [ −α{ log( t−to )/(tp −to )}2 ] という変化を示す。各画像単位値からこの式に合致する
ように、各パラメータk、Ip、α、tp 、to が最小2
乗法を用いて決定される。
【0014】求められたパラメータ或いはそれから算出
した数値からは、それぞれに対応した例えば平均循環時
間等の図7に示すようなファンクショナルイメージ画像
が作成され、ディスプレイ5に出力表示される。ファン
クショナルイメージ画像は血流状態等の視覚的理解を助
長し、例えば平均循環時間のファンクショナルイメージ
では、上述の色素希釈曲線のパラメータto 、tp 、α
から平均循環時間tmが、 tm =to +(tp −to )exp(0.75α) の理論式から決定され、各画像単位において、各平均循
環時間tm が各画像単位値として代表されて、値の大き
さに応じて色分けするなどしてディスプレイ5に出力さ
れる。出力された画像は構成が眼底画像に類似してお
り、血管部分が平均循環時間に応じて色わけされてお
り、眼底の血液循環の局所的な停留等の判別が一目で確
認でき、網膜全体の循環機能を一目で把握することを可
能とする。
した数値からは、それぞれに対応した例えば平均循環時
間等の図7に示すようなファンクショナルイメージ画像
が作成され、ディスプレイ5に出力表示される。ファン
クショナルイメージ画像は血流状態等の視覚的理解を助
長し、例えば平均循環時間のファンクショナルイメージ
では、上述の色素希釈曲線のパラメータto 、tp 、α
から平均循環時間tmが、 tm =to +(tp −to )exp(0.75α) の理論式から決定され、各画像単位において、各平均循
環時間tm が各画像単位値として代表されて、値の大き
さに応じて色分けするなどしてディスプレイ5に出力さ
れる。出力された画像は構成が眼底画像に類似してお
り、血管部分が平均循環時間に応じて色わけされてお
り、眼底の血液循環の局所的な停留等の判別が一目で確
認でき、網膜全体の循環機能を一目で把握することを可
能とする。
【0015】
【発明の効果】以上説明したように本発明に係る眼底画
像処理方法は、得られた情報を単に数値として表示する
のではなく画像として表示するので、網膜全体の機能が
視覚的に把握可能であり診断上有効である。
像処理方法は、得られた情報を単に数値として表示する
のではなく画像として表示するので、網膜全体の機能が
視覚的に把握可能であり診断上有効である。
【図1】実施例の眼底画像処理システムの構成図であ
る。
る。
【図2】時系列的な複数枚の画像の概念図である。
【図3】小矩形分に分画した状態の画像の説明図であ
る。
る。
【図4】小矩形画素の模式図である。
【図5】画素値の時間変化のグラフ図である。
【図6】色素希釈曲線に適合したグラフ図である。
【図7】ファンクションイメージ画像の説明図である。
1 眼底カメラ 2 スライドスキャナ 3 コンピュータ 4 メモリ 5 ディスプレイ 6 プリンタ 7 入力デバイス
Claims (2)
- 【請求項1】 時系列の蛍光眼底画像の各画像を複数の
領域に分割して、前記領域中で血管部に対応する部位の
画素値の平均値を求め、該平均値を基に前記領域の画素
値を算出する第1の工程と、前記領域ごとに前記画素値
の時系列変化を理論式に対応させてファンクショナルイ
メージを作成する第2の工程を有することを特徴とする
眼底画像処理方法。 - 【請求項2】 前記第1の工程の前に、前記時系列の蛍
光眼底画像の各画像同志の位置合わせを行う請求項1に
記載の眼底画像処理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP33140491A JP3169244B2 (ja) | 1991-11-19 | 1991-11-19 | 眼底画像処理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP33140491A JP3169244B2 (ja) | 1991-11-19 | 1991-11-19 | 眼底画像処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH05137696A true JPH05137696A (ja) | 1993-06-01 |
JP3169244B2 JP3169244B2 (ja) | 2001-05-21 |
Family
ID=18243313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP33140491A Expired - Fee Related JP3169244B2 (ja) | 1991-11-19 | 1991-11-19 | 眼底画像処理方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3169244B2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007068928A (ja) * | 2005-09-02 | 2007-03-22 | Tomey Corporation | 眼球表面解析システム |
JP2007252402A (ja) * | 2006-03-20 | 2007-10-04 | Topcon Corp | 眼科測定装置 |
US7306336B2 (en) | 2003-12-26 | 2007-12-11 | Nidek Co., Ltd. | Fundus observation apparatus |
US7347550B2 (en) | 2003-08-13 | 2008-03-25 | Canon Kabushiki Kaisha | Ophthalmologic image taking apparatus |
JP2013500101A (ja) * | 2009-07-28 | 2013-01-07 | エフ・ホフマン−ラ・ロシュ・アクチェンゲゼルシャフト | 非侵襲性インビボ光学イメージング方法 |
-
1991
- 1991-11-19 JP JP33140491A patent/JP3169244B2/ja not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7347550B2 (en) | 2003-08-13 | 2008-03-25 | Canon Kabushiki Kaisha | Ophthalmologic image taking apparatus |
US7566130B2 (en) | 2003-08-13 | 2009-07-28 | Canon Kabushiki Kaisha | Ophthalmologic image taking apparatus |
US7306336B2 (en) | 2003-12-26 | 2007-12-11 | Nidek Co., Ltd. | Fundus observation apparatus |
JP2007068928A (ja) * | 2005-09-02 | 2007-03-22 | Tomey Corporation | 眼球表面解析システム |
JP2007252402A (ja) * | 2006-03-20 | 2007-10-04 | Topcon Corp | 眼科測定装置 |
JP2013500101A (ja) * | 2009-07-28 | 2013-01-07 | エフ・ホフマン−ラ・ロシュ・アクチェンゲゼルシャフト | 非侵襲性インビボ光学イメージング方法 |
Also Published As
Publication number | Publication date |
---|---|
JP3169244B2 (ja) | 2001-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5123056A (en) | Whole-leg x-ray image processing and display techniques | |
US7826884B2 (en) | Live fluoroscopic roadmapping including targeted automatic pixel shift for misregistration correction | |
CA1193761A (en) | Digital fluorographic method and system | |
CN108899075A (zh) | 一种基于深度学习的dsa图像检测方法、装置及设备 | |
US20180350080A1 (en) | Analysis method and system of digital subtraction angiographic images | |
JP4958254B2 (ja) | 画像解析システム、及び画像解析プログラム | |
US8411925B2 (en) | X-ray diagnosis apparatus and image processing apparatus | |
CN107564048A (zh) | 基于分叉点特征配准方法 | |
US4802093A (en) | X-ray image-processing apparatus utilizing grayscale transformation | |
Sandor et al. | High precision quantitative angiography | |
CN106651771A (zh) | 一种数字x射线图像拼接方法及系统 | |
JP2006500099A (ja) | 画像データ情報の表示 | |
JP2018114031A (ja) | 眼底画像処理装置 | |
US20200069276A1 (en) | X-ray imaging apparatus and x-ray image processing method | |
JPH05137696A (ja) | 眼底画像処理方法 | |
JPH07239933A (ja) | 画像処理装置 | |
US8023707B2 (en) | Evaluation method for mapping the myocardium of a patient | |
CN114612484B (zh) | 基于无监督学习的视网膜oct图像分割方法 | |
JP2006075390A (ja) | 医用画像表示装置 | |
JP3407761B2 (ja) | 眼底画像解析処理装置 | |
JPS61199389A (ja) | ディジタルサブトラクション装置 | |
US20220313092A1 (en) | Method and device for monitoring the fluorescence emitted at the surface of a biological tissue | |
JPS5980085A (ja) | デイジタルサブトラクシヨンシステム | |
JPH03222588A (ja) | 画像強調方法 | |
Martinez-Perez et al. | Automatic Detection of Optic Disc and Fovea from Colour Fundus Images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |