JP3854888B2 - Atherosclerosis detection system - Google Patents

Atherosclerosis detection system Download PDF

Info

Publication number
JP3854888B2
JP3854888B2 JP2002106211A JP2002106211A JP3854888B2 JP 3854888 B2 JP3854888 B2 JP 3854888B2 JP 2002106211 A JP2002106211 A JP 2002106211A JP 2002106211 A JP2002106211 A JP 2002106211A JP 3854888 B2 JP3854888 B2 JP 3854888B2
Authority
JP
Japan
Prior art keywords
fundus
arteriosclerosis
detection system
electrocardiogram signal
fundus image
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.)
Expired - Lifetime
Application number
JP2002106211A
Other languages
Japanese (ja)
Other versions
JP2003299621A (en
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.)
BML Inc
Original Assignee
BML Inc
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 BML Inc filed Critical BML Inc
Priority to JP2002106211A priority Critical patent/JP3854888B2/en
Priority to US10/137,283 priority patent/US6999812B2/en
Publication of JP2003299621A publication Critical patent/JP2003299621A/en
Application granted granted Critical
Publication of JP3854888B2 publication Critical patent/JP3854888B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Description

【0001】
【発明の属する技術分野】
本発明は、動脈硬化の検出システムに関する発明である。
【0002】
【従来の技術】
本発明者は、高齢者が自立するためには、本人の精神的な自覚は勿論のこと、種々の疾病の発症を未然に防ぐ、予防医学的なアプローチが不可欠である、との考えのもと、循環器疾患の多くの場合の発症要因である動脈硬化を、簡便、かつ、的確に検出する手段を提供した。
【0003】
すなわち、本発明者は、心電図信号を検出することが可能な心電図信号の検出手段と、この心電図信号の特定のパターンを感知する心電図信号の感知手段と、これにより感知された心電図信号の特定のパターンに同期させて眼底像を検出可能な眼底像の検出手段を備えた、動脈硬化の検出装置を提供した(WO01/30235A1号公報)。この検出装置により、ウインドケッセル現象による、心臓の拡張期と収縮期の時相のズレにかかわらず、眼底写真による眼底動脈と眼底静脈の径の測定を、正確に行うことが可能となり、眼底像を指標として、動脈硬化の進展度合いを把握することも可能となった。
【0004】
【発明が解決しようとする課題】
本発明者は、この動脈硬化の検出装置による動脈硬化の検出が、予防医学上、非常に意義深いものであることを確信し、この装置をさらに汎用化するための手段を講ずることが必要であると考えた。
【0005】
【課題を解決するための手段】
その結果、本発明者は、上記の検出装置において、最も重要かつ困難な過程の一つである、心電図信号と眼底像の同期化を、コンピューター画面上で行うことにより、上記の課題を克服することが可能であることを見出し、本発明を完成した。
【0006】
すなわち、本発明は、心電図信号の検出手段と、これにより検出された心電図信号に同期させた眼底像を検出可能な眼底像の検出手段を備えた動脈硬化の検出システムにおいて、眼底像の検出が、眼底像の動画から、任意の心電図信号と同期させた眼底の静止像をコンピュータの表示画面上において抽出することにより、心電図信号と同期させた眼底像を提供可能なソフトウエアにより行われる、動脈硬化の検出システムであって、当該ソフトウエアが、眼底動脈と眼底静脈の交叉部近傍における眼底静脈の狭窄の度合いについての演算を、異なる心電図信号に対して行い、これらの演算値の統合値と動脈硬化の程度とを関連付けて、動脈硬化を検出可能な手段が施されているソフトウエアである、動脈硬化の検出システム(以下、本検出システムともいう)を提供する発明である。
【0007】
また、本発明は、上記の心電図信号と同期させた眼底像を提供可能なソフトウエア(以下、本ソフトウエアともいう)を提供する発明でもある。
本ソフトウエアは、コンピュータ端末における表示手段上で、眼底像と心電図の動画を同時に表示しつつ、心電図信号の特定のパターンと同期させた眼底像を抽出可能なソフトウエアであることが好適である。
【0008】
本検出システムにおける「心電図信号の検出手段」は、心電図信号を的確に検出可能な手段であれば、特に限定されず、例えば、圧電素子からなる電極センサを被検者の胸部又は他の生体部位に装着し、導出される心電図信号を検出することができる手段を挙げることができる。典型的には、既存の心電図計が具備する機構を、心電図信号の検出手段として用いることができる。
【0009】
検出する心電図信号は、心電図上の確立したパターンとして把握可能な信号であれば、特に限定されず、P波、Q波、R波、S波又はT波のいずれのパターンを選択することが可能であるが、血液を心臓から体内に向けて排出する段階のパターン信号であるR波、または、心室興奮の回復過程を示すT波を選択することが、好適であり、かつ、現実的である。本検出システムにおいては、本ソフトウエアにより、任意の心電図信号と同期させた眼底の静止像を、眼底像の動画から、コンピュータの表示画面上において抽出して、心電図信号と同期させた眼底像が提供され、任意のタイミングで同期された眼底像を抽出可能であるが、R波とT波に対応する心臓の脈波拍動タイミングは、心臓の脈波拍動サイクルの中でも、最も顕著なタイミングである。よって、本検出システムにおいて抽出する任意の心電図信号に、心電図信号のR波および/またはT波を示す信号が含まれていることは、好適な態様の一つである(これについては、さらに後述する)。
【0010】
また、被検者から心電図信号を得るための誘導法は、特に限定されず、いわゆる「標準12誘導」等から選択することが可能である。誘導方法を選択する場合の基準は、上記の選択する特定の心電図信号の種類であることが好ましい。すなわち、特定の心電図信号を検出することが、可能な限り容易な誘導方法を選択することが好ましい。特定の心電図信号として、R波を選択する場合には、被検者の左手と右手の間の電位差を検出する、II誘導、I誘導、aL 誘導、V1 誘導等を選択することが好ましい。
【0011】
「眼底像の検出手段」は、心電図の検出手段により検出された心電図信号に同期させて、眼底像を検出可能な手段である。
本発明において、「同期させる」とは、心電図信号の特定のパターンに対して一定のタイミングで眼底像の検出手段を呼応させることを意味する。例えば、心電図の信号の特定パターンとして、R波を選択する場合には、R波におけるいずれかの時点、例えば、立ち上がり時点から、一定のタイミングで眼底像の検出手段を作動させることを意味する。このタイミングは、一定に保たれており、かつ、心電図信号の同一のパターンが再び発生するタイミング(例えば、R波であれば、次のR波が発生する時点)よりも短いタイミングであれば、特に限定されない。このように、眼底像を心電図信号に同期させて検出することにより、眼底血管についての情報、具体的には、本発明における動脈硬化についての指標を得る上で不可欠な血管径についての情報を的確に得ることができる。すなわち、前述したように、随時任意に撮影された眼底写真では、ウインドケッセル現象により、心臓の拡張と収縮に対応して変化する眼底血管径について、正確に評価することが困難であったが、眼底像を心電図信号に同期させて検出すれば、一定の脈波拍動タイミングにおける眼底血管像を得ることが可能であり、眼底血管径についても、正確な評価を行うことが可能となる。
【0012】
眼底像の検出手段として、眼底を撮影可能な機構を備えるカメラ(具体的には、いわゆる眼底カメラが挙げられる:アナログカメラであってもデジタルカメラであってもよい)が挙げらるが、本検出システムにおいては、眼底像を、デジタルイメージ情報として連続的に得ることができるデジタルビデオカメラで検出することが、コンピュータにおける、眼底像の心電図信号との同期化に適している。
【0013】
デジタルビデオカメラにより撮影された眼底像の動画を、例えば、DV端子(メディアコンバーターも可能)と、IEEE1394カード、EZDV(カノープス社)、DVRapter(カノープス社)、DVRex(カノープス社)等のDVキャプチャカード等を介してデジタル情報としてコンピューターに取込みつつ、A/D変換器等によりデジタル信号に変換された心電図信号を、コンピューターに取込む。次に、取り込んだ眼底像の動画データーと心電図信号のデーターの、並列複合化を行うことで、眼底像の動画データーと心電図信号を、同一のフレームにおいて同期させて、眼底像の動画データーと心電図信号のデジタル同期化データーを得ることができる。このデジタル同期化データーは、本検出システムを行う上で必要な要素を損なわない限りにおいて、圧縮化を行うことが可能である。かかる圧縮化を含めた符号化は、MPEG等の符号化方式に従うことで行うことができる。
【0014】
このようにして得られた、デジタル同期化データは、例えば、磁気テープ、磁気ディスク、CD−ROM、MO、DVD−R等に保存することができる。なお、これらの記憶媒体は、あくまで例示であり、データ記憶技術の変遷や向上に応じたデータ記憶媒体を用いることが可能である。
【0015】
このように得られたデジタル同期化データーにおける、眼底像に関する演算値(後述するVレシオ等)の算出は、かかるデーターを静止画像、すなわち、1フレーム単位のデジタルデータを抽出することにより行われる。
【0016】
後述するVレシオ等の演算値の、単位時間当りの変化量を求める場合も、任意の心電図信号(時刻t)における、眼底像の画像データーを、動画データーから抽出し、さらに、適切な時間を置いた時点(t+Δt)における眼底像の画像データーを動画データーから抽出し、両者の静止画像データを基に算出することができる(後述するように、この場合のtも、眼底像の動画に同期させた心電図信号に依存させて選択することが好適である)。
【0017】
眼底像の動画データーと心電図信号のデータのデジタル同期化データからの、静止画像データーの抽出は、コンピュータ端末における、コンピューターディスプレイ等の表示手段上で、眼底像と心電図の動画を同時に表示しつつ行うことで、抽出作業を視覚化することが可能であり、かつ、好適である。よって、本ソフトウエアには、このコンピュータ端末の表示手段上における視覚化手段が、アルゴリズムとして搭載されていることが好ましい。
【0018】
本ソフトウエアは、一般的なコンピュータプログラム言語により、所望するアルゴリズムを構築して作出することができる。
コンピュータプログラム言語として、例えば、機械語、アセンブラ言語等の低水準言語;Fortran、ALGOL、COBOL、C、BASIC、PL/I、Pascal、LISP、Prolog、APL、Ada、Smalltalk、C++、Java(登録商標)等の高水準言語;第4世代言語、エンドユーザー言語等を選択して用いることが可能である。また、必要に応じて、特殊問題向き言語を用いることもできる。
【0019】
本ソフトウエアを格納可能な電子媒体は、特に限定されず、例えば、磁気テープ、磁気ディスク、CD−ROM、MO、DVD−R等を用いることができる。本発明は、この本ソフトウエアを格納した電子媒体をも提供する(これらの記憶媒体は、あくまで例示であり、データ記憶技術の変遷や向上に応じたデータ記憶媒体を用いることが可能である)。
【0020】
前述したように、本検出システムにより、一定の脈波拍動タイミングにおける眼底血管像を得ることが可能である。この眼底血管像から得られる情報を、本ソフトウエアにおけるアルゴリズムにより、被検者の動脈硬化の進展と関連付けることにより、被検者の動脈硬化を検出することができる。本発明における「動脈硬化の検出」とは、被検者の動脈の血管形態の器質的変化の程度と、伸展可能度や機能的変化の程度、言い換えれば、動脈のしなやかさ(弾力性)の程度の検出を意味するものである。すなわち、本検出システムにより、動脈硬化の進展が認められる場合とは、被検者の動脈の器質的変化に加えて、弾力性が低下して、その伸展可能度が低くなり、動脈のしなやかさが失われつつある場合である。そして、このような場合には、例えば、虚血性疾患のリスクが高くなっていることを意味している。
【0021】
上記の関連付けの態様は、特に限定されないが、眼底血管径についての情報が、本ソフトウエアのアルゴリズムにより、被検者の動脈の血管の器質的な変化の程度や、しなやかさの程度を検出する上で、特に有力な情報となる。具体的には、これらの眼底血管径についての情報を、本ソフトウエアにおいて、従来の動脈硬化の判定基準、例えば、前述した、「1969年の文部省高血圧研究班の条件」、「K−W分類」、「Seheie分類」等との相関を鑑みることにより、被検者の動脈の血管の器質的な変化の程度や、しなやかさの程度を検出することができる。
【0022】
本発明者は、眼底動脈と眼底静脈(本発明において、「眼底動脈」とは、医学用語にいう「網膜動脈」を意味し、「眼底静脈」とは、医学用語にいう「網膜静脈」を意味するものとする)の交叉部近傍における眼底静脈の狭窄の度合いが、動脈硬化の進展と関連付けることが可能な非常に有用な指標となり得ることを見出した。そして、一定の脈波拍動タイミングにおける眼底血管像を得ることが可能な本検出システムは、この指標を得るために非常に有用である。
【0023】
この眼底静脈の狭窄の度合いは、本検出システムにより得られる眼底像を、目視で直接観察することにより検出することも可能である。また、上述の眼底像の検出手段に、上記の眼底静脈の狭窄を検出可能な、眼底静脈の狭窄の検出手段を設けて、この過程を自動化することもできる。この眼底静脈の狭窄の検出手段としては、本ソフトウエアに、例えば、眼底静脈と眼底動脈の交叉部を検出するためのアルゴリズム及び/又は交叉部近傍における眼底静脈の狭窄の度合いを把握して算出するためのアルゴリズム等を搭載することが挙げられる。このような本ソフトウエアで、上述の眼底像のデータを処理することにより、簡便かつ確実に、所望の眼底静脈の狭窄を検出することができる。
【0024】
ここで、本発明者が見出した、動脈硬化を検出する上において有用な、眼底静脈の狭窄の指標について説明する。
眼底動脈と眼底静脈の交叉部において、これらの眼底血管は、外膜同士を共有することが知られている。そして、かかる交叉部において、動脈硬化が認められる場合には、眼底動脈の内部が硬化を起こし、この動脈の内膜の硬化による「引きつれ力」(「ズリ応力」を含む。以下、同様である。)により、眼底静脈が、交叉部近傍において狭窄を起こす。この引きつれ力による狭窄が著しいほど、交叉部近傍において、動脈硬化が進展していることを表している。
【0025】
本発明者は、本来の眼底静脈の外径と、眼底動脈と眼底静脈の交叉部近傍における眼底静脈の外径との比(以下、Vレシオともいう)が、動脈硬化を示す眼底静脈の交叉部近傍における狭窄の指標として、極めて有用であることを見出した。
【0026】
具体的には、Vレシオは、以下のようにして算出する。
第1図は、眼底静脈の引きつれ狭窄が認められる、眼底動脈と眼底静脈の模式図であり、眼底動脈10の下側を眼底静脈20が交叉し、交叉部30近傍において、眼底動脈10の内膜が硬化を起こし(図示せず)、その硬化による引きつれ力により、眼底静脈20の交叉部30近傍における狭窄が認められる。
【0027】
眼底動脈10の外径をΦAとして、交叉部30から、眼底静脈20の末梢側に3×ΦA程度離れた部位を確認して、その位置の眼底静脈20の外径を、交叉部30において引きつれ狭窄が全く認められない場合の眼底静脈の外径とみなして、これをV2として、交叉部30の眼底動脈10の直側部(例えば、交叉部30から1/10×ΦA程度の距離)の眼底静脈20の外径をV1とする。
【0028】
この場合のV1とV2の比(通常は、V1/V2)を、Vレシオとして算出する。Vレシオ(V1/V2)が小さい程、交叉部30近傍における、眼底動脈10の硬化が進展していることを示し、動脈硬化のリスクが高いことを表している。
【0029】
第2図は、眼底動脈11の血管内膜の硬化が、交叉部31より、少しズレた方向に拡がっている場合(硬化部分:111及び112)を示す模式図であり、眼底静脈21において、交叉部31近傍で、引きつれ力が互い違いに働いた結果、眼底静脈21が交叉部31近傍で、位置のズレを起こしていることを示している。
【0030】
このような場合には、上述したVレシオ(V1/V2)を、若干補正する必要がある。具体的には、交叉部31において、眼底動脈11を挟んだ、眼底静脈21の左右の外径の中央を結んだベクトル線(211)が、眼底静脈21の外径の0.5〜1.0倍ズレている場合には、上述したVレシオの値から0.1を引いた補正値を動脈硬化の検出に用いるべきVレシオとする。また、上記の外径のズレが、眼底静脈21の外径の1.0倍を超えてズレている場合には、上述したVレシオの値から0.2を引いた補正値を、動脈硬化の検出に用いるべき、総合Vレシオとする。
【0031】
また、第3図は、第2図に示す眼底静脈のズレとは若干異なり、交差する眼底動脈において、末梢側とその反対側の眼底静脈の長さ方向の中心線の間に、角度のズレがある場合を示している。すなわち、第3図においては、眼底動脈12を交差する眼底静脈22において、眼底静脈12の末梢側の長さ方向の中心線221と、反対側の長さ方向の中心線222とが、第2図のように平行ではなく、ズレ角度が認められる。この場合のズレ角度αを、末梢静脈側の交叉部32と反対側の交叉部33から1/10×ΦA程度の距離の、外径線321と331と、各々の静脈の長さ方向の中心線の交点同士を結んだ場合の、各々の長さ方向の中心線との間に認められる角度と定義する。
【0032】
このような場合にも、上述したVレシオ(V1/V2)を、若干補正する必要がある。具体的には、上記で得られたVレシオから、ズレ角度αに、1/100を乗じた値をS値とすると、このS値を、上述したVレシオから減ずることにより得られる補正値を、動脈硬化の検出に用いるべき、総合Vレシオとする。
【0033】
なお、このような補正処理をすべき、ズレ角度αは、好適には10°(絶対値)以上であり、さらに好適には5°(絶対値)以上である。このように、本ソフトウエアが、眼底動脈と眼底静脈の交叉部近傍における眼底静脈の狭窄の度合いについて、総合Vレシオ(交叉部に、ズレやズレ角度が認められず、上記の2種類の補正が行われていない「Vレシオ」を含む。以下、特に断わらない限り、総合Vレシオを、補正が行われていないVレシオを含むものとする)等を算出する演算を行い、この演算値と動脈硬化の程度とを関連付けて、動脈硬化を検出可能な手段が施されている態様であることは、本検出システムの最も好適な態様の一つである。
【0034】
本検出システムにおいて、複数(2か所以上、通常は、2〜3か所)の、乳頭部より1乳頭径以上離れた眼底動脈と眼底静脈の交叉部における総合Vレシオ等の演算値を、例えば、平均値として統合して、これを被検者の総合Vレシオ等として、動脈硬化の検出指標とすることは、演算値の信頼性を向上させる上で好適である。
【0035】
また、本ソフトウエアのアルゴリズムの実行により得られる演算値が、異なる心電図信号に対して得られる演算値の統合値とすることも、本検出システムにおいては、任意のタイミングで同期させた眼底像のデータを抽出することが可能であるため、容易である。例えば、心電図信号のR波とT波における眼底像における総合Vレシオ等の演算値の平均値として統合して、これを被検者の総合Vレシオ等として、動脈硬化の検出指標とすることも、演算値の信頼性を向上させる上で、好適である。
【0036】
さらに、この総合Vレシオ等により求められる、眼底動脈と眼底静脈の交叉部近傍における眼底静脈の狭窄の度合いを、複数の被検者の年令に応じて標準化処理して得られる標準化情報と、被検者各々における前記の狭窄の度合いとを、相互の値の高低等を比較することにより関連付けるアルゴリズムを、本ソフトウエアに搭載して、被検者各々の動脈における血管の形態の器質的変化と血管拍動を、機能的に評価することで、動脈の老化の程度(いわば、血管年齢)を検出することも可能である〔例えば、被検者の年齢と性別の標準化情報よりも、被検者の総合Vレシオが小さければ、被検者の血管は、標準よりも動脈硬化が進展していると判断され、実際の被検者の年齢よりも、動脈の血管の形態の器質と血管拍動の機能的な評価として高齢であると評価されることとなる(いわば、血管年齢が、実年齢よりも高いとして評価される)。逆に、被検者の年齢と性別の標準化情報よりも、被検者の総合Vレシオが大きければ、被検者の血管は、標準よりも動脈硬化が進展していないと判断され、動脈の血管の形態の器質と血管拍動の機能的な評価として良好であると評価されることとなる(いわば、血管年齢が、実年齢よりも若いとして評価される)。
【0037】
なお、本ソフトウエアに、異なる心電図信号間における演算値の変化率を算出することが可能なアルゴリズムが搭載されていることも好適な態様の一つである。
【0038】
この変化率を求める場合、予め、ターゲット信号とした、異なる2種類の心電図信号に同期させて得た眼底像における、総合Vレシオ等の演算値の差を求めることにより把握することができる。典型的には、心臓の脈波拍動サイクルの中でも、最も顕著なタイミングであるR波のピーク部分と、T波の終了部分に同期させて得た眼底像における、総合Vレシオの差を表す数値を、上記の変化率の指標として挙げることができる。
【0039】
さらに、総合Vレシオの変化率としては、総合Vレシオの単位時間当りの変化量を求めることも効果的である。すなわち、総合Vレシオをyとして、最初の総合Vレシオの測定時刻をt、次の総合Vレシオの測定時刻までの経過時間をΔtとすると、y=f(t)として表される、総合Vレシオyの時刻tに対する関数を、一次関数f(t)=at+bに近似させることができる程度に、経過時間Δtを小さくした場合の、一次関数f(t)=at+bの傾きaの絶対値を表す数値を、総合Vレシオの変化率として用いることも可能である。この場合、変化率が大きければ、血管に弾力性が認められることとなり、小さければ、血管に弾力性がなく、硬直化していることとなる。このような、総合Vレシオの変化率は、上記の血管拍動の機能的な評価を行う一要素として用いることができる。すなわち、このようなアルゴリズムの実行により、血管に弾力性が認められれば、動脈の血管の器質と血管拍動の機能的な評価が高く判断される要素となり、血管の硬直化が認められれば、同評価が低く判断される要素となる。
【0040】
なお、本検出システムにおいては、時刻tとして、任意の時刻を設定することができるが、単位時間毎の総合Vレシオの変化量は、同一被検者においても、ウインドケッセル現象による、血管脈拍の変動周期に依存しており、この変動周期の一周期内において異なっている。すなわち、選択した時刻tを、眼底像と同期させた心電図信号で特定して、眼底像を、血管脈拍の変動周期の基準となる心電図信号に依存させ、血管脈拍の変動周期内における総合Vレシオの変化量の本来的な差異による測定誤差を解消することが、総合Vレシオの変化を、動脈硬化の進展度合いと、より正確に関連付けるために必要となる。よって、時刻tを任意に選択することができるといっても、この時刻tは、眼底像に同期させた心電図信号に依存させて選択することが必要であり、心電図信号の典型的な波信号であるR波やT波を基準とした時刻として設定することが、好適、かつ、現実的である(上述した、R波のピーク部分とT波の終了部分に同期させた総合Vレシオの差異は、時刻tをR波のピーク部分とT波の終了部分として選択し、時間Δtを、R波のピーク部分からT波の終了部分にまでにかかる時間、または、T波の終了部分からR波のピーク部分にまでにかかる時間として選択して算出される値である)。この時間Δtにおける総合Vレシオの相違が、動脈硬化の程度を推測する、血管の拍動性変化率を規定することとなる。
【0041】
【発明の実施の形態】
以下、本発明の実施の形態を説明する。
第4図は、本検出システムの構成を示すブロック図である。
【0042】
本検出システム40は、心電図信号の検出部41、眼底像の検出部43、および、コンピューター44で構成されている。
本検出システム40は、心電図信号と眼底像の同期化等を、コンピュータ44において行う、本発明の動脈硬化の検出システムの実施態様の一つを示した図面である。
【0043】
さらに、心電図信号の検出部41は、電極センサ411、増幅部412、及び、出力部413で構成されている。さらに、眼底像の検出部43は、DV撮像部431及び出力部432で構成されている。
【0044】
心電図信号の検出部41において、電極センサ411は、例えば、圧電素子からなり、被検者の胸部又は他の生体部位に装着され、導出される心電図信号を検出する機構であり、増幅部412は、電極センサ411で検出された心電図信号の増幅処理を行う機構であり、出力部413は、増幅処理された心電図信号の出力を行う機構である。
【0045】
本検出システム40では、同期化処理をコンピューター44において行うために、心電図信号の検出部41の出力部413から、直接、コンピューター44の入力部441に、心電図信号が入力される、この心電図信号は、A/D変換器(414)等により、デジタル化処理が行われていることが好適である。
【0046】
また、眼底像の検出部43では、被検者の眼底像をDV撮像部(デジタルビデオカメラの撮像部)により撮像を行い、これにより得られる眼底像の動画のイメージ信号を抽出し、この動画信号を、出力部432からDV端子を介して、コンピューター44に、入力部442からDVキャプチャカード等を介して入力する。なお、DV撮像部431のデジタルビデオカメラは、眼底静脈径の微妙な変化を測定する必要上、可能な限り高い解像度であることが好適である。具体的には、200万画素以上の解像度を有することが好適である。また、DV撮像部431には、通常の眼底カメラが具備する、被検者の眼底を撮像するための機構、例えば、接眼レンズ、光源、アライメント機構、画角調整機構等を、必要に応じて具備していることは勿論である。
【0047】
コンピューター44に入力された、眼底像の動画のデジタル信号は、コンピューター44の処理装置443において、入力部441から入力された眼底像の動画のイメージ信号の並列複合化を行うことで、眼底像の動画データーと心電図信号を、同一のフレーム毎に同期させて(同期化処理4431)、眼底像の動画データーと心電図信号のデジタル同期化データー(4432)を得ることができる。同期化データー4432は、必要に応じた、圧縮等の処理が行われていてもよい。
【0048】
また、この際、デジタル化された心電図信号及び/又は眼底像の動画のデジタル信号を、逆量子化、補間等の過程を経て、コンピューター44の表示手段において、心電図及び/又は眼底像として表示や出力を行うことが可能とすることが好適である。
【0049】
なお、上述したように、同期化データー4432は、そのまま、その後の眼底静脈径の測定等の工程に用いることが可能であり、一旦、電子媒体に保存することも可能である。
【0050】
基礎データの計測工程4433は、同期化データー4432を基にして、少なくとも1か所の眼底静脈と眼底動脈の交叉部を、ターゲット部位として選択して、かかるターゲット部位における眼底動脈と眼底静脈の血管径等、動脈硬化についての判断を行うための基礎データを計測する工程である。
【0051】
なお、これらの基礎データは、各ターゲット部位において、異なるタイミングで計測することが好適である。この異なるタイミングは、得られるべき血管径等のデータの変化を感知することができることを、最小のタイミングとして、自由に設定することができる。
【0052】
解析工程4434は、基礎データの計測工程4433において計測された、ターゲット部位である眼底動脈と眼底静脈の交叉部の血管径等の基礎データを基に、適切な解析をすることによって、この基礎データをを、さらに有益な情報へと変換することができる。この有益な情報とは、例えば、被検者の、前述した総合Vレシオ、さらには血管年令等が挙げられる。
【0053】
また、異なるタイミングで、基礎データを計測した場合に、各ターゲット部位における各タイミング間の総合Vレシオ等の変化を計測することにより、眼底静脈径の変化を計測することにより、心電図信号に依存した単位時間当りの総合Vレシオ等の変化を算出することが有益である。すなわち、前述したように、この変化を測定することにより、被検者の血管の器質的な変化の程度や、動脈のしなやかさの程度を算出することが可能である。
【0054】
このようにして、本検出システム40では、コンピューター44において得られる、心電図信号に同期された眼底像、すなわち、ウインドケッセル現象に左右されない安定した眼底像を、連続的に得ることが可能であり、これにより、総合Vレシオ等の動脈硬化の指標となるデーターを、多用にサンプリングして、これらの統合値を算出することで、得られる総合Vレシオ等の信頼性を向上させることが可能であり、さらに、被検者の血管の器質的な変化や、動脈のしなやかさに関する指標を、容易に得ることが可能となる。
【0055】
第5図(1)〜(12)は、本検出システム40のコンピューター44の処理装置において用いる本ソフトウエアのフローシートの一実施例(400)を示した図面である。本実施例における、各プロセスにおけるコンピュータ44の操作は、必要に応じてマウス操作とすることが好適である(本例においては、「マウス操作」の表示により、その一例が示してある)。
【0056】
第5図(1)において、0000「スタート」は、コンピューター44を、フローシート400に示す処理を行う本ソフトウエアを実行することが可能な状態にセットアップすることを示している。セットアップ後、初期メニュー画面を呼出し(0001)、映像出力を、コンピューター44の表示手段において行い、例えば、「登録,修正,表示,集計,終了」等の作動メニューを表示する(0002)。
【0057】
次いで、初期メニューの選択を行い(0003)、この初期メニューの選択ステップ(0004)において、例えば、「登録」を選択し、初期データーを入力すると、「初期データ入力画面」が表示され(0005)、「表示」を選択すると、「データ確認表示画面」が表示され(0006)、「集計」を選択すると、「データ集計画面」が表示される(0007)〔「表示」、または、「集計」選択後のプロセス(エ)(オ)については、後述する。また、「修正」を選択すると、データ修正のための修正メニューが表示され、「終了」を選択すると、終了確認プロセス実行後、プログラムの実行が終了する〕。
【0058】
初期データ入力画面の画像表示(0008・0009)後、初期データ(カルテ番号,氏名,性別,年齢,血圧値,コレステロール値,血糖値等)のキーボード等からの入力を行い(0010)、初期データの入力を完了する(0011)。
【0059】
次に、眼底像の撮影準備プロセスについて説明する。眼底像の撮影に際しては、まず、眼底像の検出部443における撮影準備を行い(0100)、入力部441からコンピューター44に入力された心電図信号を、コンピューター44の表示手段における表示確認を行う(0101)。
【0060】
心電図の確認終了(0102)後、撮影条件の選択画面の表示を行い(0103)、眼底像の検出部43における眼底像の撮影画角(50°または25°)の選択を行い、検出部43では、選択された撮影画角に応じた撮影条件が設定される(0105,0106)。
【0061】
ステップ0106で撮影画角50°を選択すると(撮像プロセスA)、眼底像の検出部43において、画角50°で眼底が撮影され(0200)、その画像のデジタルデータが、入力部442を経て、コンピューター44に入力され、その表示手段において、映像出力表示が行われる(0201,0202)。この入力された眼底画像のデジタルデータと心電図信号との並列複合化により、眼底像の動画データーと心電図信号を、同一のフレーム毎に同期させて、眼底像の動画データーと心電図信号のデジタル同期化データーを得ることができる。この撮像プロセスA終了後、再び、撮影条件選択表示プロセス(0103)に戻り、必要に応じて、撮像プロセスAを複数回行うことも可能である。また、この撮像プロセスA完了後、さらに、特定部分の拡大画面を得るために、撮影画角25°における撮像プロセスを含めたプロセス(撮像プロセスB)に移行することもできる。
【0062】
第5図(2)(3)は、この撮像画角25°における撮像プロセスを含めた一部を示している。まず、撮像画角50°における眼底画像から、動脈硬化の候補部分を、例えば、3ポイント選択して(0203)、各ポイントを選択する(例えば、各ポイントを、Aポイント、Bポイント、および、Cポイントとする。このポイントの選択にあたっては、映像出力を行い、例えば、各ポイントを、×印で表示する(0204,0205)。各ポイントを選択後、この撮像画角50°における各ポイントの、特定の心電図信号と同期されている撮影データを保存する(0207,0208)。この保存に際しては、例えば、上述したカルテ番号、血圧値等の初期データを、撮像データと共に保存することが好適である。
【0063】
次に、上述のように選択したAポイントの解析プロセスを行うために、上記の撮像画角50°による眼底画像の表示を行い(0209,0210)、改めて、Aポイント近傍を、画角25°で撮影し(0211)、Aポイントにおける動脈硬化部分を選択するためのマウス操作を行い(0212)、動脈硬化部分と認めた箇所を、改めて撮影画角25°での撮影を行い(0213)、この撮影に係るAポイントの映像出力を行い(0214)、その映像データの保存を行う(0215〜0217)。
【0064】
次に、BポイントとCポイントの画角25°における眼底画像の映像出力〜保存のプロセスを、上記のAポイントの場合と同様の手順で行う(Bポイント:0218〜0227、Cポイント:0228〜0237)。
【0065】
次に、眼底静脈の外径の測定に関するプロセスについて説明する〔第5図(3)〕。まず、ポイント選択画面の表示を行う(0238)。具体的には、ポイント選択画面の映像出力を行い(0239)、本例では、与えられたポイントである、Aポイント、Bポイント、および、Cポイントのうち、測定ポイントを選択する(0240)。Aポイントを選択する場合には(0241)、マウス操作により、測定画面からAポイント部分の表示を行い(0242,0243)、動脈硬化部分をマウス操作により拡大して(0244)、この拡大部分の映像出力を行う(0245〜0247)。
【0066】
次に、眼底動脈の外径の測定ステップに移行する〔第5図(4):具体的には、Aポイントについてのみ記載するが、BポイントとCポイントについても、同様のステップを行うことができる〕。
【0067】
眼底動脈の外径の測定ステップ(0248)は、上記のAポイントの拡大部分の表示を行い、マウスにより、対象となる眼底動脈10の外径をドラックすることにより行われる(0249)。眼底動脈10の径の測定値をA1として(0250)、動脈外径A1の1/10の値を、1/10A1として(0251)、測定値の表示を行う(0252)。再び、当該部分の映像出力を行い(0253)、眼底動脈と眼底静脈の交差部分における、眼底動脈10の眼底静脈の末梢側の外径のクリックを行う(0255)。次に、この部分から、眼底静脈の末梢側へ1/10A1だけ離れた位置と数値の表示を行い(0256)、映像出力を行う(0257)。
【0068】
次に、第1の眼底静脈の外径の測定ステップ(0259)に移行する。映像出力(0257)を行った後、さらに、マウス操作により、前記の末梢側から1/10A1だけ離れた位置における眼底静脈径(V1)を、マウスでドラックを行うことにより測定し、V1の値を表示しつつ、映像出力を行う(0260〜0264)。
【0069】
上記の1/10A1だけ離れた位置よりも、さらに眼底動脈から離れた位置における末梢側眼底静脈径の測定ステップ(0265)に移行する。まず、眼底動脈と眼底静脈の交差部分における、眼底動脈の眼底静脈の末梢側の外径のクリックを行う(0266)。さらに、眼底動脈径A1の3倍の距離を計算し、これをA2とし、眼底静脈の末梢側A2の位置と数字の表示を行い、映像出力を行う(0267〜0270)。
【0070】
次に、第2の眼底静脈の外径の測定プロセス(0271)に移行する。映像出力を行った後、さらに、マウス操作により、前記の末梢側からA2だけ離れた位置における眼底静脈径(V2)を、マウスでドラックを行うことにより測定し〔第5図(5)0272,0273)、V2の値を表示しつつ(0274)、映像出力を行う(0275)。
【0071】
次に、Vレシオの計算プロセス(0300)に移行する。すなわち、上述のようにして得られたV1とV2の値を基に、V1/V2をVレシオとして計算して(0301)、これの映像表示を行う(0302,0303)と共に、Vレシオをデータとして保存する(0303〜0306)。
【0072】
次に、眼底静脈の末梢側の中心線の測定プロセス(0401)に移行する。このステップでは、まず、V1をクリックして、その値を呼び出し(0402,0403)、この値に1/2を乗じた値であるX1を算出する(0405)。次いで、この算出されたX1を基に、V1をドラックして得た線分の中点(X1)の映像出力を行う(0406,0407)。さらに、V2をクリックして、その値を呼び出し(0408,0409)、この値に1/2を乗じた値であるX2を算出する(0410,0411)。次いで、この算出されたX2を基に、V2をドラックして得た線分の中点(X2)の映像出力を行う(0412,0413)。そして、最後に、中心点X1とX2を結んだベクトル線X1X2を表示して(0414)、映像出力を行う(0415)。
【0073】
次に、眼底静脈の末梢側と反対側の中心線の測定プロセス(0416)に移行する。このステップは、上述した眼底静脈の中心線の測定プロセス(0401)における手順と、実質的に同様の手順で行われる。
【0074】
まず、Aポイントの眼底動脈と眼底静脈の交差部分における、眼底動脈の眼底静脈の末梢側の反対側の外径のクリックを行う〔第5図(6)0417〕。次に、この外径から、眼底静脈の末梢側の反対側へ1/10A1だけ離れた位置と数値の表示を行い(0418,0419)、映像出力を行う(0420)。
【0075】
次に、マウス操作により、前記の末梢側の反対側から1/10A1だけ離れた位置における眼底静脈径(V3)を、マウスでドラックを行うことにより測定し(0421,0422)、V3の値を表示しつつ、映像出力を行う(0423,0424)。さらに、V3の中心を測定するために(0425)、V3をクリックして、その値を呼び出し(0426,0427)、この値に1/2を乗じた値であるY1を算出する(0428)。次いで、この算出されたY1を基に、V3をドラックして得た線分の中点(Y1)の映像出力を行う(0429,0430)。
【0076】
次に、マウス操作により、前記の末梢側の反対側からA2だけ離れた位置(A4)を表示し(0431,0432)、このA4における眼底静脈径(V4)を測定する(0435)。すなわち、マウスでドラックを行うことにより、V4を測定し(0436,0437)、V4の値を表示しつつ、映像出力を行う(0438,0439)。さらに、V4をクリックして、その値を呼び出し(0440,0441)、この値に1/2を乗じた値であるY2を算出する(0442)。次いで、この算出されたY2を基に、V4をドラックして得た線分の中点(Y2)の映像出力を行う(0443)。そして、最後に、中心点Y1とY2を結んだベクトル線Y1Y2を表示して、映像出力を行う〔第5図(7)0444,0445〕。
【0077】
以上のようにして得られた、ベクトル線X1X2とY1Y2の記憶を行い、データとして保存する(0446〜0448)。
次に、総合Vレシオを算出するための、眼底静脈のズレの測定プロセスに移行する。
【0078】
まず、ベクトル線X1X2とY1Y2の間の平均距離を算出し、この値をXY1とし(0500)、眼底静脈径(V4とV2の平均値:V5)で除した値、すなわち、XY1/V5=XY1a、を算出し(0501)、記憶を行う(0502,0503)。
【0079】
次に、眼底静脈のズレ角度の測定プロセス(0600)に移行する。メニューから「角度測定」を選んでクリックし(0601)、マウス操作により、X1とY1を結び(0602,0603)、このベクトル線X1Y1の映像出力を行う(0604,0605)。末梢側静脈のズレを判定するために(0606)、再び、メニューから「角度測定」を選んでクリックし(0607)、ベクトル線Y1Y2を、ズレ角度測定の基準線としてクリックし(0608)、マウス操作によりベクトル線X1Y1を指定し(0609,0610)、ベクトル線Y1Y2とX1Y1の角度(絶対値)を測定する(0611)。
【0080】
ステップ612において、このズレ角度(絶対値)が、5°以上の場合(YES)には、ズレ角度を補正要素として、総合Vレシオを判定するプロセス(ズレ角度補正プロセス)に移行する。
【0081】
ズレ角度補正プロセスは、上記の角度の測定プロセスにおいて測定された、ベクトル線Y1Y2とX1Y1の角度を、Vレシオの補正パラメーターとして、適切な数値(S値)に変換して、この数値でVレシオを補正して、総合Vレシオを判定するプロセスである(0700)。まず、測定角度/100を算出する〔第5図(8)0701〕。例えば、測定角度が10°であれば、S値は0.1であり、同15°であれば、0.15である。次に、VレシオからS値を減じた値を、「総合Vレシオ」として算出し(0702)、測定角度とS値と総合Vレシオの映像出力表示を行い(0703,0704)、これらのデータを保存する(0705)ために、マウス操作(0706)を行い、測定角度等のデータの記憶を行う(0707)。
【0082】
ここで、これまでの測定により得られた測定データの一覧表示を行うことができる(0708)。すなわち、マウス操作により、「測定データの表示」をクリックすると(0709,0710)、上述したAポイントのV1値、V2値、Vレシオ、総合Vレシオ、その他の測定データをはじめ、Aポイントと同様の工程により得られたBポイントとCポイントにおける、これらの測定データが一覧表示される(0711,0712)。このようにして、ズレ角度補正プロセスを行うことができる。この、ズレ角度補正プロセスの完了後、後述する「標準化情報と測定された総合Vレシオとの比較判定プロセス」に移行する〔プロセス(イ)〕。
【0083】
次に、ステップ612において、ズレ角度が5°未満の場合(NO)は、プロセス(ア)に移行し、下記の総合Vレシオを判定するプロセス(ズレ補正プロセス)が行われる。
【0084】
ズレ補正プロセスは、ステップ0501で算出されたXY1aをもとに、Vレシオにおける補正値を算出して、Vレシオを、総合Vレシオとして補正する、総合Vレシオの判定プロセス(0800)に移行する。
【0085】
まず、ステップ0801,0802は、XY1aが0.5以上か否かを選別するステップであり、NOの場合、すなわち、XY1aが0.5未満の場合は、補正値は0であり、ステップ0301で算出されたVレシオは、そのまま、データとして保存される〔プロセス(ウ)〕。これに対して、YESの場合、すなわち、XY1aが0.5以上の場合は、さらに、XY1aが1.0以上か否かを判定するステップ(0803,0804)に移行する。このステップがNOの場合、すなわち、XY1aが0.5〜1.0未満の場合は、補正値は−0.1であり、ステップ0301で算出されたVレシオから0.1を減じた値(V11:0808)が、より正確に動脈硬化に進展度に対応する値であり、この値が、総合VレシオV11として記憶され、保存される(V22:0809,0810)。このステップがYESの場合、すなわち、XY1aが1.0以上の場合は、補正値は−0.2であり、ステップ0301で算出されたVレシオから0.2を減じた値(0805)が、より正確に動脈硬化に進展度に対応する値であり、この値が、総合VレシオV22として記憶され、保存される(0806,0807)。このようにして、ズレ補正プロセスを行うことができる。この、ズレ補正プロセスの完了後、後述する「標準化情報と測定総合Vレシオとの比較判定プロセス」に移行する。
【0086】
次に、標準化情報と測定総合Vレシオとの比較判定プロセス〔第5図(9)0900〕について説明する。予め、データとして入力された標準化情報(年齢に対応した、標準的な総合Vレシオの値)(0902)と上記で測定された、測定総合Vレシオ情報(0901)とを比較する(0903)。すなわち、ステップ0904は、標準化情報と測定総合Vレシオ情報の大小比較を行うステップであり、このステップで、NO、すなわち、標準化情報よりも測定総合Vレシオ情報が大きい場合は、被検者の現実の年齢よりも、若い血管年齢が検出され(0907)、これが、血管年齢Bとして記憶される(0908,0909)。これに対して、ステップ0904で、YES、すなわち、標準化情報と測定総合Vレシオ情報が同一か、標準化情報よりも測定総合Vレシオ情報が小さい場合には、被検者の現実の年齢と同一か、老いた血管年齢が検出され(0905)、これが、血管年齢Aとして記憶される(0906,0910)。
【0087】
次に、血管柔軟度の測定プロセス(1000)に移行する。このプロセスにおいては、R波のピーク時に同期させた眼底像から得られる総合Vレシオと、T波の終了時に同期させた眼底像から得られる総合Vレシオを測定して、両者を比較することにより、血管の柔軟度を測定するステップである。
【0088】
まず、心電図信号のR波のピーク時に同期させた眼底像の表示を行う(1001)。すなわち、マウス操作(1002)により、R波のピーク時に同期させた眼底像(撮影画角50°)の映像出力を行い(1003〜1006)、測定画面の選択を行う(1007)。前述したように、ここでは、Aポイントの選択を行う(1008)。すなわち、マウス操作(1009)により、Aポイント部分の表示を行い(1010)、さらに、動脈硬化部分を拡大した画像の映像出力を行う(撮影画角25°)〔第5図(10)1011〜1014〕。次いで、動脈硬化部分の眼底動脈10の外径(AMX2)を、マウスによりドラッグすることにより測定し、その測定値の映像出力を行う(1015〜1019)。
【0089】
さらに、心電図信号のT波の終了時に同期させた眼底像の表示を行う(1020)。すなわち、マウス操作(1021)により、この心電図信号のT波の終了時に同期させた眼底像(撮影画角50°)の映像出力を行い(1022〜1025)、測定画面の選択を行う(1026)。前記と同様、Aポイントの選択を行う(1027)。すなわち、マウス操作(1028)により、Aポイント部分の表示を行い(1029)、さらに、動脈硬化部分を拡大した画像の映像出力を行う(撮影画角25°)(1030〜1033)。次いで、動脈硬化部分の眼底動脈10の外径(AMN2)を、マウスによりドラッグすることにより測定し、その測定値の映像出力を行う〔第5図(11)1034〜1038〕。
【0090】
次に、これらのプロセスにおいて得られた、AMX2とAMN2の差(FXレシオ)を求める(1039)。このFXレシオを、血管の弾力性を示す指標として用いることができる。最後に、FXレシオの数値と、これに伴う画像を保存する(1040,1041)。
【0091】
このようにして、本検出システム40による、被検者の動脈硬化の検出プロセスの一実施態様が終了する(1042)。
データ確認表示プロセス
上述したフローシートの初期メニューの選択ステップ(0004)において、「表示」を選択した場合のデータ確認プロセス(2000)について説明する。
【0092】
初期メニューの選択ステップ(0004)において、「表示」を選択し、プロセス(エ)に移行すると、キーボード入力画面が表れ(2001)、被検者のID番号のキーボード入力を行うと(2002)、対象の被検者の各データが映像出力される(2003)。ここで、対象の被検者の各データの内容を確認して判定するステップに移行する(2004)。判定ステップ2005において、NO、すなわち、対象の被検者の各データの内容に誤りがある等、修正の必要性がある場合には、上記のキーボード入力画面において、修正情報のキーボード入力を行い、上記のプロセスを繰り返す〔繰り返しプロセス(a)〕。これに対して、YES、対象の被検者の各データの内容に修正の必要が認められない場合には、各ポイント(Aポイント、Bポイント、Cポイント)の映像出力がなされ(2006,2007)、各ポイントが、画面上に表示される。これらのポイントの中から、所望するポイント画面を選択してクリックし(2008,2009)、マウス操作(2010)により拡大画面の映像出力表示を行い(2011)、拡大画面と各種データの表示を行い〔第5図(12)2012〕、データ確認表示ステップが終了する。
【0093】
データ一覧表示プロセス
上述したフローシートの初期メニューの選択ステップ(0004)において、「集計」を選択した場合のデータの一覧表示のプロセス(3000)について説明する。
【0094】
初期メニューの選択ステップ(0004)において、「集計」を選択し、プロセス(オ)に移行すると(3001)、データ表示画面が表れると共に、データ集計様式のメニューが表れる。これらのメニューの中から、所望する集計様式を選択すると(3002)、各数値項目により並べ変えられた一覧表示の映像出力が行われる(3003)。このようなデータの一覧表示を行い、データ一覧表示ステップが終了する。
【0095】
【発明の効果】
本発明により、簡便かつ鋭敏な、動脈硬化の検出手段が提供される。
【図面の簡単な説明】
【図1】眼底静脈の引きつれ狭窄が認められる、眼底動脈と眼底静脈の模式図である。
【図2】眼底動脈の血管内膜の硬化が、交叉部より、少しズレた方向に拡がっている場合を示す模式図である。
【図3】眼底動脈と眼底静脈の交叉部において、眼底静脈にズレ角度が認められる場合を示す模式図である。
【図4】本検出システムの構成を示すブロック図である。
【図5】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第1番目である。
【図6】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第2番目である。
【図7】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第3番目である。
【図8】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第4番目である。
【図9】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第5番目である。
【図10】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第6番目である。
【図11】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第7番目である。
【図12】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第8番目である。
【図13】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第9番目である。
【図14】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第10番目である。
【図15】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第11番目である。
【図16】本検出システムのコンピューターの処理装置において用いる本ソフトウエアのフローシートの一実施例を示した図面の第12番目である。
【符号の説明】
40:本検出システム
41:心電図信号の検出部
43:眼底像の検出部
44:コンピュータ
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an arteriosclerosis detection system.
[0002]
[Prior art]
The present inventor believes that in order for the elderly to become independent, a preventive medical approach that prevents the onset of various diseases as well as the mental awareness of the individual is essential. And a means for simply and accurately detecting arteriosclerosis, which is a causative factor in many cases of cardiovascular disease.
[0003]
That is, the inventor of the present invention provides an electrocardiogram signal detection means capable of detecting an electrocardiogram signal, an electrocardiogram signal detection means for sensing a specific pattern of the electrocardiogram signal, and a specific electrocardiogram signal sensed thereby. An arteriosclerosis detection device provided with a fundus image detection means capable of detecting a fundus image in synchronization with a pattern has been provided (WO01 / 30235A1). This detection device makes it possible to accurately measure the diameter of the fundus artery and fundus vein using a fundus photograph, regardless of the time difference between the diastolic and systolic phases of the heart due to the wind kessel phenomenon. It was also possible to grasp the degree of progression of arteriosclerosis using as an index.
[0004]
[Problems to be solved by the invention]
The inventor is convinced that the detection of arteriosclerosis by this arteriosclerosis detection device is very significant in preventive medicine, and it is necessary to take measures to further generalize this device. I thought it was.
[0005]
[Means for Solving the Problems]
As a result, the present inventor overcomes the above problem by performing synchronization of the electrocardiogram signal and the fundus image on the computer screen, which is one of the most important and difficult processes in the above detection device. The present invention has been completed.
[0006]
  That is, the present invention relates to an arteriosclerosis detection system comprising an electrocardiogram signal detection means and a fundus image detection means capable of detecting a fundus image synchronized with the electrocardiogram signal detected thereby. Artery performed by software capable of providing a fundus image synchronized with an electrocardiogram signal by extracting a still image of the fundus synchronized with an arbitrary electrocardiogram signal from a moving image of the fundus image on a computer display screen. A curing detection system,The software calculates the degree of stenosis of the fundus vein near the intersection of the fundus artery and fundus vein for different ECG signals, and associates the integrated value of these calculations with the degree of arteriosclerosis. , Software with a means to detect arteriosclerosisThe present invention provides an arteriosclerosis detection system (hereinafter also referred to as this detection system).
[0007]
The present invention is also an invention that provides software (hereinafter also referred to as this software) that can provide a fundus image synchronized with the above-described electrocardiogram signal.
The software is preferably software that can extract a fundus image synchronized with a specific pattern of an electrocardiogram signal while simultaneously displaying a fundus image and an electrocardiogram video on a display means in a computer terminal. .
[0008]
The “electrocardiogram signal detection means” in the present detection system is not particularly limited as long as it can accurately detect the electrocardiogram signal. For example, an electrode sensor composed of a piezoelectric element is connected to the chest of the subject or other living body part. And a means capable of detecting the derived electrocardiogram signal. Typically, a mechanism provided in an existing electrocardiograph can be used as a means for detecting an electrocardiogram signal.
[0009]
The electrocardiogram signal to be detected is not particularly limited as long as it is a signal that can be grasped as an established pattern on the electrocardiogram, and any pattern of P wave, Q wave, R wave, S wave, or T wave can be selected. However, it is preferable and realistic to select the R wave, which is a pattern signal at the stage of discharging blood from the heart toward the body, or the T wave indicating the recovery process of ventricular excitation. . In this detection system, this software extracts a still image of the fundus synchronized with an arbitrary electrocardiogram signal from a moving image of the fundus image on the display screen of the computer, and the fundus image synchronized with the electrocardiogram signal is obtained. It is possible to extract a fundus image that is provided and synchronized at an arbitrary timing, but the pulse wave pulsation timing of the heart corresponding to the R wave and the T wave is the most prominent timing in the pulse wave pulsation cycle of the heart. It is. Therefore, it is one of preferred embodiments that an arbitrary ECG signal extracted in this detection system includes a signal indicating the R wave and / or T wave of the ECG signal (this will be further described later). To do).
[0010]
In addition, the guidance method for obtaining an electrocardiogram signal from the subject is not particularly limited, and can be selected from so-called “standard 12 leads” or the like. It is preferable that the reference for selecting the guidance method is the type of the specific electrocardiogram signal to be selected. That is, it is preferable to select a guidance method that is as easy as possible to detect a specific electrocardiogram signal. When an R wave is selected as a specific electrocardiogram signal, a potential difference between the left hand and the right hand of the subject is detected, leads II, I,aVLInduction, V1It is preferable to select induction or the like.
[0011]
The “fundus image detection means” is a means capable of detecting a fundus image in synchronization with an electrocardiogram signal detected by an electrocardiogram detection means.
In the present invention, “synchronize” means that the fundus image detection means is made to respond at a certain timing to a specific pattern of the electrocardiogram signal. For example, when an R wave is selected as a specific pattern of an electrocardiogram signal, this means that the fundus image detection means is operated at a certain timing from any time point in the R wave, for example, the rising time point. This timing is kept constant, and if the timing is shorter than the timing at which the same pattern of the electrocardiogram signal is generated again (for example, if the R wave is the next R wave is generated) There is no particular limitation. In this way, by detecting the fundus image in synchronization with the electrocardiogram signal, information on the fundus blood vessel, specifically, information on the blood vessel diameter indispensable for obtaining an index of arteriosclerosis in the present invention can be accurately obtained. Can get to. That is, as described above, in the fundus photograph taken arbitrarily at any time, it was difficult to accurately evaluate the fundus blood vessel diameter that changes in accordance with the expansion and contraction of the heart due to the wind kessel phenomenon. If the fundus image is detected in synchronization with the electrocardiogram signal, it is possible to obtain a fundus blood vessel image at a constant pulse wave pulsation timing, and it is possible to accurately evaluate the fundus blood vessel diameter.
[0012]
As a fundus image detection means, a camera having a mechanism capable of photographing the fundus (specifically, a so-called fundus camera: an analog camera or a digital camera) may be mentioned. In the system, detection of a fundus image with a digital video camera that can be obtained continuously as digital image information is suitable for synchronization with an electrocardiogram signal of the fundus image in a computer.
[0013]
For example, a DV capture card such as a DV terminal (media converter is also possible), IEEE1394 card, EZDV (Canopus), DVRapter (Canopus), DVRex (Canopus), etc. The electrocardiogram signal converted into a digital signal by an A / D converter or the like is taken into the computer while taking it into the computer as digital information through the above. Next, the fundus image moving image data and the electrocardiogram signal are synchronized in the same frame by performing parallel compounding of the captured fundus image moving image data and the electrocardiogram signal data, so that the fundus image moving image data and the electrocardiogram are synchronized. Digitally synchronized data of the signal can be obtained. This digital synchronization data can be compressed as long as the elements necessary for carrying out the present detection system are not impaired. Encoding including such compression can be performed by following an encoding method such as MPEG.
[0014]
The digital synchronization data obtained in this way can be stored, for example, on a magnetic tape, magnetic disk, CD-ROM, MO, DVD-R, or the like. Note that these storage media are merely examples, and it is possible to use data storage media according to changes and improvements in data storage technology.
[0015]
In the digital synchronization data obtained in this way, calculation values (such as a V ratio described later) regarding the fundus image are calculated by extracting still data, that is, digital data in units of one frame.
[0016]
Also when calculating the amount of change per unit time of a calculated value such as a V ratio, which will be described later, image data of the fundus image at an arbitrary electrocardiogram signal (time t) is extracted from the moving image data, and an appropriate time is determined. The image data of the fundus image at the time of placement (t + Δt) can be extracted from the moving image data and calculated based on the still image data of both (as will be described later, t in this case is also synchronized with the moving image of the fundus image) It is preferable to select depending on the electrocardiogram signal.
[0017]
Extraction of still image data from digitally synchronized data of fundus image moving image data and electrocardiogram signal data is performed while simultaneously displaying the fundus image and the electrocardiogram moving image on a display means such as a computer display in a computer terminal. Thus, it is possible and preferable to visualize the extraction work. Therefore, it is preferable that the present software is equipped with the visualization means on the display means of the computer terminal as an algorithm.
[0018]
This software can be created by constructing a desired algorithm in a general computer program language.
As a computer program language, for example, low-level languages such as machine language and assembler language; Fortran, ALGOL, COBOL, C, BASIC, PL / I, Pascal, LISP, Prolog, APL, Ada, Smalltalk, C ++, Java (registered trademark) ) Or the like; a fourth generation language, an end user language, or the like can be selected and used. A special problem language can also be used as necessary.
[0019]
The electronic medium that can store the software is not particularly limited, and for example, a magnetic tape, a magnetic disk, a CD-ROM, an MO, a DVD-R, or the like can be used. The present invention also provides an electronic medium in which the software is stored (these storage media are merely examples, and data storage media can be used in accordance with changes and improvements in data storage technology). .
[0020]
As described above, the present detection system can obtain a fundus blood vessel image at a constant pulse wave pulsation timing. The arteriosclerosis of the subject can be detected by associating the information obtained from the fundus blood vessel image with the progress of the arteriosclerosis of the subject by the algorithm in this software. “Detection of arteriosclerosis” in the present invention refers to the degree of organic change in the blood vessel morphology of the subject's artery and the degree of extension or functional change, in other words, the flexibility of the artery (elasticity). This means detection of the degree. In other words, when this detection system shows the progression of arteriosclerosis, in addition to the organic changes in the subject's arteries, the elasticity decreases and the degree of extension becomes low, and the flexibility of the arteries. Is being lost. In such a case, for example, it means that the risk of ischemic disease is high.
[0021]
The manner of the above association is not particularly limited, but the information about the fundus blood vessel diameter detects the degree of organic changes in the blood vessels of the subject's arteries and the degree of flexibility using the algorithm of this software. Above, it is particularly useful information. Specifically, information about these fundus blood vessel diameters is obtained in this software using the conventional criteria for arteriosclerosis, for example, “the conditions of the 1969 Ministry of Education, Hypertension Research Group”, “K-W classification”. ”,“ Seheie classification ”, and the like, the degree of organic changes in the arterial blood vessels of the subject and the degree of flexibility can be detected.
[0022]
The present inventor has determined that the fundus artery and the fundus vein (in the present invention, “fundus artery” means “retinal artery” in medical terms, and “retinal vein” means “retinal vein” in medical terms). It has been found that the degree of stenosis of the fundus vein in the vicinity of the intersection (which shall mean) can be a very useful indicator that can be correlated with the progression of arteriosclerosis. The present detection system capable of obtaining a fundus blood vessel image at a constant pulse wave pulsation timing is very useful for obtaining this index.
[0023]
The degree of stenosis of the fundus vein can also be detected by directly observing the fundus image obtained by this detection system visually. Further, the fundus image detection means described above can be provided with a fundus vein stenosis detection means capable of detecting the fundus vein stenosis, and this process can be automated. As a means for detecting stenosis of the fundus vein, for example, the software detects, for example, an algorithm for detecting the intersection of the fundus vein and the fundus artery and / or the degree of stenosis of the fundus vein near the intersection. Mounting an algorithm or the like for doing so. By processing the above-described fundus image data with such software, it is possible to detect a desired stenosis of the fundus vein simply and reliably.
[0024]
Here, the index of stenosis of the fundus vein that is useful in detecting arteriosclerosis found by the present inventor will be described.
It is known that these fundus blood vessels share the adventitia at the intersection of the fundus artery and the fundus vein. When arteriosclerosis is observed at such a crossing portion, the inside of the fundus artery is hardened, and includes “tension force” (“slip stress”) due to hardening of the intima of the artery. In other words, the fundus vein causes stenosis in the vicinity of the intersection. As the stenosis due to the pulling force is more remarkable, it indicates that the arteriosclerosis has progressed in the vicinity of the intersection.
[0025]
The inventor found that the ratio of the outer diameter of the original fundus vein to the outer diameter of the fundus vein in the vicinity of the intersection of the fundus artery and the fundus vein (hereinafter also referred to as V ratio) is the intersection of the fundus vein showing arteriosclerosis. It was found to be extremely useful as an index of stenosis in the vicinity of the head.
[0026]
Specifically, the V ratio is calculated as follows.
FIG. 1 is a schematic diagram of the fundus artery and the fundus vein, in which the fundus vein is constricted, and the fundus vein 20 crosses the lower side of the fundus artery 10, and the fundus artery 10 is located near the crossover 30. The intima is cured (not shown), and a stenosis in the vicinity of the intersection 30 of the fundus vein 20 is recognized due to the pulling force due to the curing.
[0027]
  Assuming that the outer diameter of the fundus artery 10 is ΦA, a site separated by about 3 × ΦA from the crossover portion 30 to the distal side of the fundus vein 20 is confirmed, and the outer diameter of the fundus vein 20 at that position is pulled at the crossover portion 30. Considering this as the outer diameter of the fundus vein when no stenosis is observed, this is taken as V2 and the right side of the fundus artery 10 of the intersection 30 (for example, a distance of about 1/10 × ΦA from the intersection 30) Of the fundus vein 20Outer diameterIs V1.
[0028]
The ratio of V1 and V2 in this case (usually V1 / V2) is calculated as the V ratio. The smaller the V ratio (V1 / V2), the more the hardening of the fundus artery 10 in the vicinity of the intersection 30 is shown, indicating that the risk of arteriosclerosis is high.
[0029]
FIG. 2 is a schematic diagram showing a case where the hardening of the intima of the fundus artery 11 spreads in a direction slightly shifted from the crossing portion 31 (hardening portions: 111 and 112). As a result of the pulling force acting alternately in the vicinity of the crossing portion 31, it is shown that the fundus vein 21 is displaced in the vicinity of the crossing portion 31.
[0030]
In such a case, it is necessary to slightly correct the V ratio (V1 / V2) described above. Specifically, a vector line (211) connecting the centers of the left and right outer diameters of the fundus vein 21 across the fundus artery 11 at the intersection 31 is 0.5 to 1. When the deviation is 0 times, a correction value obtained by subtracting 0.1 from the above-described V ratio value is set as a V ratio to be used for detection of arteriosclerosis. Further, when the deviation of the outer diameter is more than 1.0 times the outer diameter of the fundus vein 21, the correction value obtained by subtracting 0.2 from the value of the V ratio is used as the arteriosclerosis. The total V ratio that should be used for detection of.
[0031]
3 is slightly different from the deviation of the fundus vein shown in FIG. 2, and in the intersecting fundus artery, there is an angle deviation between the distal center side and the center line in the length direction of the opposite fundus vein. Shows when there is. That is, in FIG. 3, in the fundus vein 22 crossing the fundus artery 12, the longitudinal center line 221 on the distal side of the fundus vein 12 and the longitudinal center line 222 on the opposite side are As shown in the figure, a shift angle is recognized rather than parallel. In this case, the deviation angle α is set to a distance of about 1/10 × ΦA from the crossing portion 33 opposite to the crossing portion 32 on the peripheral vein side, and the center in the length direction of each vein. It is defined as an angle that is recognized between the center lines in the length direction when the intersections of the lines are connected.
[0032]
Even in such a case, it is necessary to slightly correct the V ratio (V1 / V2) described above. Specifically, when a value obtained by multiplying the deviation angle α by 1/100 from the V ratio obtained above is an S value, a correction value obtained by subtracting the S value from the V ratio described above is obtained. The total V ratio to be used for detecting arteriosclerosis.
[0033]
  It should be noted that the deviation angle α to be subjected to such correction processing is preferably 10 ° (absolute value).more thanMore preferably, 5 ° (absolute value)more thanIt is. In this way, the software can determine the degree of stenosis of the fundus vein in the vicinity of the intersection of the fundus artery and fundus vein. In the following, unless otherwise specified, the total V ratio is assumed to include the uncorrected V ratio), and the calculated value and arteriosclerosis are calculated. It is one of the most preferable aspects of the present detection system that the means capable of detecting arteriosclerosis is provided in association with the degree of the above.
[0034]
In this detection system, a plurality of (two or more, usually two to three) computed values such as the total V ratio at the intersection of the fundus artery and the fundus vein that are separated from the nipple by one nipple diameter, For example, integrating as an average value and using this as a subject's total V ratio or the like as a detection index of arteriosclerosis is suitable for improving the reliability of the calculated value.
[0035]
In addition, the detection value obtained by executing the software algorithm may be an integrated value of calculation values obtained for different ECG signals. Since data can be extracted, it is easy. For example, the average value of the calculated values such as the total V ratio in the fundus image of the R wave and T wave of the electrocardiogram signal may be integrated, and this may be used as the detection index of arteriosclerosis as the total V ratio of the subject. This is suitable for improving the reliability of the operation value.
[0036]
  Furthermore, the standardization information obtained by standardizing the degree of stenosis of the fundus vein in the vicinity of the intersection of the fundus artery and the fundus vein, which is obtained by the total V ratio and the like, according to the age of a plurality of subjects, An algorithm that correlates the degree of stenosis in each subject by comparing the level of each other's value, etc. is installed in this software, and the organic changes in the morphology of blood vessels in the arteries of each subject It is also possible to detect the degree of arterial aging (in other words, blood vessel age) by functionally evaluating the blood vessel pulsation [for example, rather than the standardized information on the age and sex of the subject. The examiner's overall V ratio isIf smallThe subject's blood vessels are judged to have developed arteriosclerosis more than the standard, and the elderly as a functional assessment of the arterial blood vessel morphology and vascular pulsation, rather than the actual subject's age (In other words, the blood vessel age is evaluated as being higher than the actual age). Conversely, rather than the standardized information of the subject's age and gender, the subject's total V ratio isIf bigThe subject's blood vessels are judged to have not developed arteriosclerosis more than the standard, and are evaluated as being good as a functional evaluation of the arterial blood vessel morphology and blood vessel pulsation ( In other words, blood vessel age is evaluated as younger than actual age).
[0037]
In addition, it is one of the preferable aspects that this software is equipped with an algorithm capable of calculating the change rate of the operation value between different ECG signals.
[0038]
When this rate of change is obtained, it can be grasped in advance by obtaining the difference between the calculated values such as the total V ratio in the fundus image obtained in synchronization with two different types of electrocardiogram signals as target signals. Typically, it represents the difference between the total V ratio in the fundus image obtained in synchronization with the peak portion of the R wave, which is the most prominent timing, and the end portion of the T wave in the pulse wave cycle of the heart. A numerical value can be cited as an index of the change rate.
[0039]
Furthermore, as the rate of change of the overall V ratio, it is also effective to obtain the amount of change per unit time of the overall V ratio. That is, if the total V ratio is y, the first total V ratio measurement time is t, and the elapsed time until the next total V ratio measurement time is Δt, the total V ratio is expressed as y = f (t). The absolute value of the slope a of the linear function f (t) = at + b when the elapsed time Δt is reduced to such an extent that the function of the ratio y with respect to time t can be approximated to the linear function f (t) = at + b. It is also possible to use the numerical value to represent as the rate of change of the overall V ratio. In this case, if the rate of change is large, the blood vessel has elasticity, and if it is small, the blood vessel has no elasticity and is rigid. Such a change rate of the total V ratio can be used as one element for performing the functional evaluation of the vascular pulsation. In other words, if the elasticity of the blood vessel is recognized by the execution of such an algorithm, the functional evaluation of the arterial blood vessel structure and blood vessel pulsation is highly judged, and if the rigidity of the blood vessel is recognized, It becomes an element that the evaluation is judged low.
[0040]
In the present detection system, an arbitrary time can be set as the time t. However, the amount of change in the total V ratio per unit time is the same as that of the blood vessel pulse due to the wind Kessel phenomenon. It depends on the fluctuation period and is different within one period of this fluctuation period. That is, the selected time t is specified by an electrocardiogram signal synchronized with the fundus image, and the fundus image is made to depend on the electrocardiogram signal which is a reference of the fluctuation cycle of the vascular pulse, and the total V ratio in the fluctuation cycle of the vascular pulse. It is necessary to eliminate the measurement error due to the inherent difference in the amount of change in order to more accurately associate the change in the overall V ratio with the progress of arteriosclerosis. Therefore, even if the time t can be arbitrarily selected, it is necessary to select the time t depending on the electrocardiogram signal synchronized with the fundus image, and a typical wave signal of the electrocardiogram signal. It is preferable and realistic to set the time based on the R wave and T wave (the difference between the total V ratio synchronized with the peak portion of the R wave and the end portion of the T wave described above. Selects time t as the peak portion of the R wave and the end portion of the T wave, and sets the time Δt to the time taken from the peak portion of the R wave to the end portion of the T wave, or from the end portion of the T wave to R This is a value that is selected and calculated as the time it takes to reach the peak of the wave). The difference in the overall V ratio at this time Δt defines the rate of change in the pulsatility of the blood vessel that estimates the degree of arteriosclerosis.
[0041]
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the present invention will be described below.
FIG. 4 is a block diagram showing the configuration of the present detection system.
[0042]
The detection system 40 includes an electrocardiogram signal detection unit 41, a fundus image detection unit 43, and a computer 44.
This detection system 40 is a drawing showing one embodiment of an arteriosclerosis detection system of the present invention in which an electrocardiogram signal and a fundus image are synchronized in a computer 44.
[0043]
The electrocardiogram signal detection unit 41 includes an electrode sensor 411, an amplification unit 412, and an output unit 413. Further, the fundus image detection unit 43 includes a DV imaging unit 431 and an output unit 432.
[0044]
In the electrocardiogram signal detection unit 41, the electrode sensor 411 is formed of, for example, a piezoelectric element, and is a mechanism that detects an electrocardiogram signal that is attached to and is derived from the chest or other living body of the subject. The amplification unit 412 The electrocardiogram signal detected by the electrode sensor 411 is amplified, and the output unit 413 is a mechanism that outputs the amplified electrocardiogram signal.
[0045]
In this detection system 40, in order to perform synchronization processing in the computer 44, an electrocardiogram signal is directly input from the output unit 413 of the electrocardiogram signal detection unit 41 to the input unit 441 of the computer 44. It is preferable that the digitization processing is performed by an A / D converter (414) or the like.
[0046]
The fundus image detection unit 43 captures the fundus image of the subject using a DV imaging unit (imaging unit of a digital video camera), extracts a video signal of the fundus image obtained thereby, A signal is input from the output unit 432 to the computer 44 via the DV terminal and from the input unit 442 via a DV capture card or the like. Note that the digital video camera of the DV imaging unit 431 preferably has a resolution as high as possible in order to measure a subtle change in the fundus vein diameter. Specifically, it is preferable to have a resolution of 2 million pixels or more. Further, the DV imaging unit 431 includes a mechanism for imaging the fundus of the subject, such as an eyepiece lens, a light source, an alignment mechanism, and an angle-of-view adjustment mechanism, which is included in a normal fundus camera, as necessary. Of course, it has.
[0047]
The digital image signal of the fundus image input to the computer 44 is combined with the image signal of the video image of the fundus image input from the input unit 441 in the processing unit 443 of the computer 44, so that the fundus image The moving image data and the electrocardiogram signal are synchronized every same frame (synchronization processing 4431), and the fundus image moving image data and the electrocardiogram signal digital synchronization data (4432) can be obtained. The synchronization data 4432 may be subjected to processing such as compression as necessary.
[0048]
At this time, the digitized electrocardiogram signal and / or the digital image signal of the fundus image are displayed as an electrocardiogram and / or fundus image on the display means of the computer 44 through processes such as inverse quantization and interpolation. It is preferable that output can be performed.
[0049]
As described above, the synchronization data 4432 can be used as it is in subsequent steps such as measurement of the fundus venous diameter, and can be temporarily stored in an electronic medium.
[0050]
In the basic data measurement step 4433, based on the synchronization data 4432, at least one intersection of the fundus vein and the fundus artery is selected as a target site, and the blood vessels of the fundus artery and the fundus vein in the target site are selected. This is a step of measuring basic data for determining arteriosclerosis such as diameter.
[0051]
These basic data are preferably measured at different timings in each target region. This different timing can be freely set as a minimum timing that a change in data such as a blood vessel diameter to be obtained can be detected.
[0052]
In the analysis step 4434, the basic data is obtained by performing appropriate analysis based on the basic data such as the blood vessel diameter of the intersection of the fundus artery and the fundus vein, which is the target site, measured in the basic data measurement step 4433. Can be converted into more useful information. This useful information includes, for example, the above-described total V ratio, blood vessel age, and the like of the subject.
[0053]
In addition, when basic data was measured at different timings, it was dependent on the electrocardiogram signal by measuring changes in the fundus vein diameter by measuring changes in the overall V ratio between timings at each target site. It is useful to calculate changes such as the total V ratio per unit time. That is, as described above, by measuring this change, it is possible to calculate the degree of organic change in the blood vessel of the subject and the degree of flexibility of the artery.
[0054]
In this manner, in the present detection system 40, it is possible to continuously obtain a fundus image synchronized with an electrocardiogram signal obtained by the computer 44, that is, a stable fundus image not affected by the wind Kessel phenomenon, As a result, it is possible to improve the reliability of the obtained overall V ratio, etc. by sampling a large amount of data that is an index of arteriosclerosis, such as the overall V ratio, and calculating these integrated values. Furthermore, it is possible to easily obtain an index relating to the organic changes in the blood vessels of the subject and the flexibility of the artery.
[0055]
FIGS. 5 (1) to (12) are drawings showing an embodiment (400) of the flow sheet of the software used in the processing device of the computer 44 of the detection system 40. FIG. In this embodiment, the operation of the computer 44 in each process is preferably a mouse operation as necessary (in this example, an example is shown by displaying “mouse operation”).
[0056]
In FIG. 5 (1), 0000 “start” indicates that the computer 44 is set up in a state where the software for performing the processing shown in the flow sheet 400 can be executed. After the setup, the initial menu screen is called (0001), and video output is performed on the display means of the computer 44, for example, an operation menu such as “registration, correction, display, aggregation, end” is displayed (0002).
[0057]
Next, an initial menu is selected (0003). In this initial menu selection step (0004), for example, when “registration” is selected and initial data is input, an “initial data input screen” is displayed (0005). When “Display” is selected, “Data Confirmation Display Screen” is displayed (0006), and when “Total” is selected, “Data Total Screen” is displayed (0007) [“Display” or “Total” Process (d) (e) after selection will be described later. Further, when “correction” is selected, a correction menu for data correction is displayed, and when “end” is selected, execution of the program ends after execution of the end confirmation process.
[0058]
After displaying the image on the initial data input screen (0008/0009), input the initial data (medical chart number, name, gender, age, blood pressure level, cholesterol level, blood glucose level, etc.) from the keyboard, etc. (0010). Is completed (0011).
[0059]
Next, the fundus image photographing preparation process will be described. When photographing a fundus image, first, a preparation for photographing is performed in the fundus image detection unit 443 (0100), and an electrocardiogram signal input from the input unit 441 to the computer 44 is confirmed on the display means of the computer 44 (0101). ).
[0060]
After the end of the electrocardiogram confirmation (0102), an imaging condition selection screen is displayed (0103), the fundus image detection field angle (50 ° or 25 °) is selected by the fundus image detection unit 43, and the detection unit 43 is detected. Then, shooting conditions corresponding to the selected shooting angle of view are set (0105, 0106).
[0061]
When a shooting field angle of 50 ° is selected in step 0106 (imaging process A), the fundus image detection unit 43 captures the fundus at an angle of view of 50 ° (0200), and digital data of the image passes through the input unit 442. The image is input to the computer 44, and the video output display is performed on the display means (0201, 0202). By digitally combining the input digital data of the fundus image and the electrocardiogram signal, the moving image data of the fundus image and the electrocardiogram signal are synchronized every same frame, so that the moving image data of the fundus image and the electrocardiogram signal are digitally synchronized. Data can be obtained. After the completion of the imaging process A, it is possible to return to the imaging condition selection display process (0103) again, and the imaging process A can be performed a plurality of times as necessary. Further, after the completion of the imaging process A, it is possible to shift to a process (imaging process B) including an imaging process at an imaging angle of view of 25 ° in order to obtain an enlarged screen of a specific part.
[0062]
FIGS. 5 (2) and 5 (3) show a part including the imaging process at the imaging angle of view of 25 °. First, from a fundus image at an imaging angle of view of 50 °, for example, three points are selected for candidate portions of arteriosclerosis (0203), and each point is selected (for example, each point is designated as A point, B point, and When selecting this point, video is output and, for example, each point is displayed with a cross (0204, 0205) After each point is selected, each point at this imaging angle of view of 50 ° is displayed. The imaging data synchronized with a specific electrocardiogram signal is saved (0207, 0208) In this saving, for example, it is preferable to save the initial data such as the above-described chart number and blood pressure value together with the imaging data. is there.
[0063]
Next, in order to perform the analysis process of the A point selected as described above, the fundus image is displayed at the above-mentioned imaging angle of view of 50 ° (0209, 0210), and the vicinity of the A point is again displayed at the angle of view of 25 °. (0211), a mouse operation for selecting an arteriosclerotic portion at point A is performed (0212), and a portion recognized as an arteriosclerotic portion is again imaged at an imaging angle of view of 25 ° (0213). The video of the A point related to this shooting is output (0214), and the video data is saved (0215 to 0217).
[0064]
Next, the process of outputting and storing the fundus image at the angle of view of 25 points of the B point and the C point is performed in the same procedure as the case of the A point (B point: 0218 to 0227, C point: 0228 to 0237).
[0065]
Next, a process related to the measurement of the outer diameter of the fundus vein will be described [FIG. 5 (3)]. First, the point selection screen is displayed (0238). Specifically, video output of the point selection screen is performed (0239), and in this example, a measurement point is selected from the given points A point, B point, and C point (0240). When the A point is selected (0241), the A point portion is displayed from the measurement screen by operating the mouse (0242, 0243), and the arteriosclerotic portion is enlarged by operating the mouse (0244). Video output is performed (0245 to 0247).
[0066]
Next, the process proceeds to a step of measuring the outer diameter of the fundus artery [FIG. 5 (4): Specifically, only the A point is described, but the same step can be performed for the B point and the C point. it can〕.
[0067]
The outer diameter measurement step (0248) of the fundus artery is performed by displaying the enlarged portion of the A point and dragging the outer diameter of the target fundus artery 10 with the mouse (0249). The measured value is displayed by setting the measured value of the diameter of the fundus artery 10 as A1 (0250), 1/10 of the outer diameter of the artery A1 as 1 / 10A1 (0251). Again, the image of the part is output (0253), and the outer diameter of the fundus artery 10 on the peripheral side of the fundus artery 10 is clicked at the intersection of the fundus artery and the fundus vein (0255). Next, from this portion, a position and a numerical value that are separated by 1 / 10A1 to the distal side of the fundus vein are displayed (0256), and video output is performed (0257).
[0068]
Next, the process proceeds to the outer diameter measuring step (0259) of the first fundus vein. After the video output (0257), the fundus vein diameter (V1) at a position 1 / 10A1 away from the peripheral side is further measured by dragging with the mouse, and the value of V1 is measured. Is output while video is displayed (0260 to 0264).
[0069]
The process proceeds to the peripheral fundus vein diameter measurement step (0265) at a position further away from the fundus artery than the position separated by 1 / 10A1. First, the outer diameter on the distal side of the fundus vein of the fundus artery at the intersection of the fundus artery and fundus vein is clicked (0266). Further, a distance three times as large as the fundus artery diameter A1 is calculated, and this is set to A2, and the position and number of the distal side A2 of the fundus vein are displayed, and video output is performed (0267 to 0270).
[0070]
Next, the process proceeds to a process (0271) for measuring the outer diameter of the second fundus vein. After the video output, the diameter of the fundus vein (V2) at a position away from the peripheral side by A2 is further measured by dragging with the mouse [FIG. 5 (5) 0272, 0273), while displaying the value of V2 (0274), video output is performed (0275).
[0071]
Next, the process proceeds to the V ratio calculation process (0300). That is, based on the values of V1 and V2 obtained as described above, V1 / V2 is calculated as a V ratio (0301), and this video is displayed (0302, 0303), and the V ratio is represented by data. (0303 to 0306).
[0072]
Next, the process proceeds to the measurement process (0401) of the center line on the peripheral side of the fundus vein. In this step, first, V1 is clicked to call that value (0402,0403), and X1 which is a value obtained by multiplying this value by 1/2 is calculated (0405). Next, based on the calculated X1, video output of the midpoint (X1) of the line segment obtained by dragging V1 is performed (0406, 0407). Further, V2 is clicked to call the value (0408, 0409), and X2 that is a value obtained by multiplying this value by 1/2 is calculated (0410, 0411). Next, based on the calculated X2, video output of the midpoint (X2) of the line segment obtained by dragging V2 is performed (0412, 0413). Finally, the vector line X1X2 connecting the center points X1 and X2 is displayed (0414), and video output is performed (0415).
[0073]
Next, the process proceeds to the measurement process (0416) of the center line on the opposite side to the peripheral side of the fundus vein. This step is performed in a procedure substantially similar to the procedure in the above-described fundus vein centerline measurement process (0401).
[0074]
First, at the intersection of the fundus artery and the fundus vein at point A, the outer diameter of the fundus artery opposite to the distal side of the fundus vein is clicked (FIG. 5 (6) 0417). Next, from this outer diameter, a position and a numerical value that are separated by 1 / 10A1 from the distal side of the fundus vein to the distal side are displayed (0418, 0419), and video output is performed (0420).
[0075]
Next, by operating the mouse, the fundus vein diameter (V3) at a position 1 / 10A1 away from the opposite side of the peripheral side is measured by dragging with the mouse (0421, 0422), and the value of V3 is determined. The video is output while displaying (0423, 0424). Further, in order to measure the center of V3 (0425), V3 is clicked to call the value (0426, 0427), and Y1 which is a value obtained by multiplying this value by 1/2 is calculated (0428). Next, based on the calculated Y1, video output of the midpoint (Y1) of the line segment obtained by dragging V3 is performed (0429, 0430).
[0076]
Next, a position (A4) separated from the peripheral side by A2 is displayed by mouse operation (A4, 0432), and the fundus vein diameter (V4) at A4 is measured (0435). That is, by dragging with the mouse, V4 is measured (0436, 0437), and video output is performed while displaying the value of V4 (0438, 0439). Further, V4 is clicked to call the value (0440, 0441), and Y2 which is a value obtained by multiplying this value by 1/2 is calculated (0442). Next, based on the calculated Y2, the video output of the midpoint (Y2) of the line segment obtained by dragging V4 is performed (0443). Finally, the vector line Y1Y2 connecting the center points Y1 and Y2 is displayed and video output is performed [FIG. 5 (7) 0444, 0445].
[0077]
The vector lines X1X2 and Y1Y2 obtained as described above are stored and stored as data (0446-0448).
Next, the process proceeds to a fundus vein deviation measurement process for calculating the total V ratio.
[0078]
First, an average distance between the vector lines X1X2 and Y1Y2 is calculated, and this value is set as XY1 (0500) and divided by the fundus vein diameter (average value of V4 and V2: V5), that is, XY1 / V5 = XY1a. Are calculated (0501) and stored (0502, 0503).
[0079]
Next, the process proceeds to a measurement process (0600) of the deviation angle of the fundus vein. “Angle measurement” is selected from the menu and clicked (0601), and by operating the mouse, X1 and Y1 are connected (0602, 0603), and video output of this vector line X1Y1 is performed (0604, 0605). In order to determine the displacement of the peripheral vein (0606), again select and click "angle measurement" from the menu (0607), click the vector line Y1Y2 as the reference line for the displacement angle measurement (0608), mouse The vector line X1Y1 is designated by the operation (0609, 0610), and the angle (absolute value) between the vector lines Y1Y2 and X1Y1 is measured (0611).
[0080]
In step 612, when the deviation angle (absolute value) is 5 ° or more (YES), the process proceeds to a process of determining the total V ratio (deviation angle correction process) using the deviation angle as a correction factor.
[0081]
In the deviation angle correction process, the angle of the vector lines Y1Y2 and X1Y1 measured in the above-described angle measurement process is converted into an appropriate numerical value (S value) as a V ratio correction parameter, and this value is used as the V ratio. Is a process for determining the total V ratio (0700). First, the measurement angle / 100 is calculated [FIG. 5 (8) 0701]. For example, if the measurement angle is 10 °, the S value is 0.1, and if it is 15 °, it is 0.15. Next, a value obtained by subtracting the S value from the V ratio is calculated as an “total V ratio” (0702), and a video output display of the measurement angle, the S value, and the total V ratio is performed (0703, 0704). Is stored (0705), a mouse operation (0706) is performed and data such as a measurement angle is stored (0707).
[0082]
Here, it is possible to display a list of measurement data obtained by the previous measurements (0708). In other words, when "Display measurement data" is clicked by operating the mouse (0709, 0710), it is the same as the A point including the V1 value, V2 value, V ratio, total V ratio, and other measurement data of the above A point. These measurement data at the B point and C point obtained by the above process are displayed in a list (0711-0712). In this way, the shift angle correction process can be performed. After the deviation angle correction process is completed, the process proceeds to a “comparison / determination process between standardized information and the measured total V ratio” described later [Process (A)].
[0083]
Next, in step 612, when the misalignment angle is less than 5 ° (NO), the process proceeds to process (A), and the following process for determining the overall V ratio (deviation correction process) is performed.
[0084]
The shift correction process shifts to a comprehensive V ratio determination process (0800) in which a correction value in the V ratio is calculated based on XY1a calculated in step 0501 and the V ratio is corrected as the total V ratio. .
[0085]
First, steps 0801 and 0802 are steps for selecting whether or not XY1a is 0.5 or more. If NO, that is, if XY1a is less than 0.5, the correction value is 0, and in step 0301 The calculated V ratio is stored as data [process (c)]. On the other hand, in the case of YES, that is, when XY1a is 0.5 or more, the process proceeds to a step (0803, 0804) for determining whether or not XY1a is 1.0 or more. When this step is NO, that is, when XY1a is less than 0.5 to 1.0, the correction value is −0.1, and a value obtained by subtracting 0.1 from the V ratio calculated in step 0301 ( V11: 0808) is a value that more accurately corresponds to the degree of progression to arteriosclerosis, and this value is stored and stored as the total V ratio V11 (V22: 0809, 0810). When this step is YES, that is, when XY1a is 1.0 or more, the correction value is −0.2, and the value (0805) obtained by subtracting 0.2 from the V ratio calculated in step 0301 is It is a value that more accurately corresponds to the degree of progress in arteriosclerosis, and this value is stored and stored as the overall V ratio V22 (0806, 0807). In this way, the deviation correction process can be performed. After the deviation correction process is completed, the process proceeds to a “comparison determination process between standardized information and measurement total V ratio” described later.
[0086]
  Next, the comparison judgment process [FIG. 5 (9) 0900] between the standardized information and the measurement total V ratio will be described. The standardized information (standard total V ratio value corresponding to age) (0902) input as data in advance (0902) is compared with the measured total V ratio information (0901) measured above (0903). That is, step 0904 is a step of comparing the size of the standardized information and the measured total V ratio information. In this step, NO, that is, the measured total V ratio information is more than the standardized information.largeIn this case, the blood vessel age younger than the actual age of the subject is detected (0907), and this is stored as the blood vessel age B (0908, 0909). On the other hand, in step 0904, YES, that is, whether the standardized information and the measured total V ratio information are the same, or the measured total V ratio information is more than the standardized information.smallIn this case, the age of the blood vessel that is the same as or older than the actual age of the subject is detected (0905), and this is stored as the blood vessel age A (0906, 0910).
[0087]
Next, the process proceeds to the blood vessel flexibility measurement process (1000). In this process, the total V ratio obtained from the fundus image synchronized at the peak of the R wave and the total V ratio obtained from the fundus image synchronized at the end of the T wave are measured and compared. This is a step of measuring the flexibility of the blood vessel.
[0088]
  First, the fundus image synchronized with the peak of the R wave of the electrocardiogram signal is displayed (1001). That is, a mouse operation (1002) outputs a fundus image (shooting angle of view 50 °) synchronized at the peak of the R wave (1003 to 1006), and selects a measurement screen (1007). As described above, here, the A point is selected (1008). That is, the A point portion is displayed by the mouse operation (1009) (1010), and an image of the enlarged arteriosclerotic portion is output (imaging angle of view 25 °) [FIG. 5 (10) 1011 1014]. Then, atherosclerotic partFundus arteryThe outer diameter (AMX2) of 10 is measured by dragging with the mouse, and the measured value is output as an image (1015 to 1019).
[0089]
  Further, the fundus image synchronized with the end of the T wave of the electrocardiogram signal is displayed (1020). That is, the mouse operation (1021) outputs the image of the fundus image (shooting angle of view 50 °) synchronized at the end of the T wave of the electrocardiogram signal (1022 to 1025), and selects the measurement screen (1026). . As described above, the A point is selected (1027). That is, the A point portion is displayed by the mouse operation (1028) (1029), and further, the video output of the enlarged image of the arteriosclerotic portion is performed (imaging angle of view 25 °) (1030 to 1033). Then, atherosclerotic partFundus arteryAn outer diameter (AMN2) of 10 is measured by dragging with a mouse, and an image of the measured value is output [FIG. 5 (11) 1034 to 1038].
[0090]
Next, the difference (FX ratio) between AMX2 and AMN2 obtained in these processes is obtained (1039). This FX ratio can be used as an index indicating the elasticity of blood vessels. Finally, the numerical value of the FX ratio and the image associated therewith are stored (1040, 1041).
[0091]
Thus, one embodiment of the subject's arteriosclerosis detection process by the detection system 40 ends (1042).
Data confirmation display process
The data confirmation process (2000) when “display” is selected in the above-described flow sheet initial menu selection step (0004) will be described.
[0092]
In the initial menu selection step (0004), when “display” is selected and the process (D) is entered, a keyboard input screen appears (2001), and when the subject ID number is entered by the keyboard (2002), Each data of the subject is output as a video (2003). Here, the process proceeds to a step of confirming and determining the contents of each data of the subject subject (2004). In the determination step 2005, NO, that is, if there is an error in the contents of each data of the subject subject, there is a need for correction, the keyboard input of the correction information is performed on the keyboard input screen, The above process is repeated [repeated process (a)]. On the other hand, if YES, the content of each data of the subject is not required to be corrected, video output of each point (A point, B point, C point) is made (2006, 2007) ) Each point is displayed on the screen. From these points, a desired point screen is selected and clicked (2008, 2009), and the video output display of the enlarged screen is performed by the mouse operation (2011) (2011), and the enlarged screen and various data are displayed. [FIG. 5 (12) 2012], the data confirmation display step ends.
[0093]
Data list display process
The process (3000) of displaying a list of data when “total” is selected in the above-described initial menu selection step (0004) of the flow sheet will be described.
[0094]
In the initial menu selection step (0004), when “total” is selected and the process proceeds to process (e) (3001), a data display screen appears and a menu for data totalization appears. When a desired tabulation form is selected from these menus (3002), video output of a list display rearranged by each numerical item is performed (3003). Such data list display is performed, and the data list display step is completed.
[0095]
【The invention's effect】
The present invention provides a simple and sensitive means for detecting arteriosclerosis.
[Brief description of the drawings]
FIG. 1 is a schematic diagram of a fundus artery and a fundus vein in which a stenosis of the fundus vein is observed.
FIG. 2 is a schematic diagram showing a case where the hardening of the vascular intima of the fundus artery spreads in a direction slightly shifted from the intersection.
FIG. 3 is a schematic diagram showing a case where a shift angle is recognized in the fundus vein at the intersection of the fundus artery and the fundus vein.
FIG. 4 is a block diagram showing a configuration of the present detection system.
FIG. 5 is the first of the drawings showing an embodiment of the flow sheet of the software used in the computer processing apparatus of the detection system.
FIG. 6 is the second of the drawings showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system.
FIG. 7 is the third drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system.
FIG. 8 is the fourth drawing showing an embodiment of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 9 is a fifth drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 10 is the sixth drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 11 is the seventh drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 12 is the eighth drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 13 is the ninth drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system;
FIG. 14 is the tenth drawing showing an example of the flow sheet of the software used in the computer processing apparatus of the detection system.
FIG. 15 is the eleventh drawing showing an embodiment of the software flow sheet used in the computer processor of the detection system.
FIG. 16 is the twelfth figure showing an embodiment of the flow sheet of the software used in the computer processing apparatus of the detection system.
[Explanation of symbols]
40: This detection system
41: ECG signal detector
43: Fundus image detection unit
44: Computer

Claims (8)

心電図信号の検出手段と、これにより検出された心電図信号に同期させた眼底像を検出可能な眼底像の検出手段を備えた動脈硬化の検出システムにおいて、眼底像の検出が、眼底像の動画から、任意の心電図信号と同期させた眼底の静止像をコンピュータの表示画面上において抽出することにより、心電図信号と同期させた眼底像を提供可能なソフトウエアにより行われる、動脈硬化の検出システムであって、当該ソフトウエアが、眼底動脈と眼底静脈の交叉部近傍における眼底静脈の狭窄の度合いについての演算を、異なる心電図信号に対して行い、これらの演算値の統合値と動脈硬化の程度とを関連付けて、動脈硬化を検出可能な手段が施されているソフトウエアである、動脈硬化の検出システム。In an arteriosclerosis detection system comprising an electrocardiogram signal detection means and a fundus image detection means capable of detecting a fundus image synchronized with the electrocardiogram signal detected thereby, the fundus image is detected from a moving image of the fundus image. This is a detection system for arteriosclerosis performed by software that can provide a fundus image synchronized with an electrocardiogram signal by extracting a still image of the fundus synchronized with an arbitrary electrocardiogram signal on a computer display screen. The software calculates the degree of stenosis of the fundus vein in the vicinity of the intersection of the fundus artery and fundus vein with respect to different ECG signals, and calculates the integrated value of these calculated values and the degree of arteriosclerosis. A system for detecting arteriosclerosis, which is software associated with a means capable of detecting arteriosclerosis. 前記動脈硬化の検出システムにおける、心電図信号と同期させた眼底像を提供可能なソフトウエアにより得られる演算値が、2か所以上の眼底動脈と眼底静脈の交叉部近傍における眼底画像から得られる演算値の統合値である、請求項1記載の動脈硬化の検出システム。In the arteriosclerosis detection system, a calculation value obtained by software capable of providing a fundus image synchronized with an electrocardiogram signal is obtained from a fundus image near the intersection of two or more fundus arteries and fundus veins. The arteriosclerosis detection system according to claim 1, which is an integrated value. 前記動脈硬化の検出システムにおける、心電図信号と同期させた眼底像を提供可能なソフトウエアに、異なる心電図信号間における演算値の変化率を算出することが可能な手段が施されている、請求項1又は2記載の動脈硬化の検出システム。The software capable of providing a fundus image synchronized with an electrocardiogram signal in the arteriosclerosis detection system is provided with means capable of calculating a rate of change of a calculation value between different electrocardiogram signals. The detection system for arteriosclerosis according to 1 or 2. 前記動脈硬化の検出システムにおける、心電図信号と同期させた眼底像を提供可能なソフトウエアに、演算値の変化率と、動脈のしなやかさの程度を関連付けて、動脈硬化を検出可能な手段が施されている、請求項3記載の動脈硬化の検出システム。In the arteriosclerosis detection system, software capable of providing a fundus image synchronized with an electrocardiogram signal is associated with a means capable of detecting arteriosclerosis by associating the change rate of the calculated value with the degree of arterial flexibility. The detection system for arteriosclerosis according to claim 3. 前記動脈硬化の検出システムにおける、心電図信号と同期させた眼底像を提供可能なソフトウエアが、コンピュータ端末における表示手段上で、眼底像と心電図の動画を同時に表示しつつ、任意の心電図信号と同期させた眼底像を抽出可能なソフトウエアである、請求項1〜4のいずれかに記載の動脈硬化の検出システム。In the arteriosclerosis detection system, software capable of providing a fundus image synchronized with an electrocardiogram signal is synchronized with an arbitrary ECG signal while simultaneously displaying a fundus image and a moving image of the electrocardiogram on the display means in the computer terminal. The system for detecting arteriosclerosis according to any one of claims 1 to 4, which is software capable of extracting a fundus image that has been obtained. 前記動脈硬化の検出システムにおいて、任意の心電図信号に、心電図信号のR波および/またはT波を示す信号が含まれている、請求項1〜5のいずれかに記載の動脈硬化の検出システム。The arteriosclerosis detection system according to any one of claims 1 to 5, wherein in the arteriosclerosis detection system, an arbitrary electrocardiogram signal includes a signal indicating an R wave and / or a T wave of the electrocardiogram signal. 請求項1〜6のいずれかに記載の動脈硬化の検出システムを行うために用いるソフトウエアが格納された電子媒体。The electronic medium in which the software used in order to perform the detection system of arteriosclerosis in any one of Claims 1-6 was stored. 請求項1〜6のいずれかに記載の動脈硬化の検出システムを行うために用いるソフトウエア。Software used to perform the arteriosclerosis detection system according to any one of claims 1 to 6.
JP2002106211A 1999-10-22 2002-04-09 Atherosclerosis detection system Expired - Lifetime JP3854888B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2002106211A JP3854888B2 (en) 2002-04-09 2002-04-09 Atherosclerosis detection system
US10/137,283 US6999812B2 (en) 1999-10-22 2002-05-03 Arteriosclerosis detection system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002106211A JP3854888B2 (en) 2002-04-09 2002-04-09 Atherosclerosis detection system

Publications (2)

Publication Number Publication Date
JP2003299621A JP2003299621A (en) 2003-10-21
JP3854888B2 true JP3854888B2 (en) 2006-12-06

Family

ID=29390599

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002106211A Expired - Lifetime JP3854888B2 (en) 1999-10-22 2002-04-09 Atherosclerosis detection system

Country Status (1)

Country Link
JP (1) JP3854888B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2005261423B2 (en) * 2004-07-08 2011-03-03 Costruzioni Strumenti Oftalmici C.S.O. S.R.L. Reflection microscope for examination of the corneal endothelium and method of operating same
DE102006018445B4 (en) * 2006-04-18 2008-04-24 Imedos Gmbh Apparatus and method for determining arterio-venous ratio values by quantitative analysis of retinal vessels
US20100234744A1 (en) * 2006-05-16 2010-09-16 Retinal Information Diagnosis Research Institute Inc. Blood vessel senescence detection system
JP2007319403A (en) * 2006-05-31 2007-12-13 Topcon Corp Medical support system, apparatus, and program
JP6346656B2 (en) * 2012-02-20 2018-06-20 キヤノン株式会社 Image processing apparatus and image processing method
WO2019125026A1 (en) * 2017-12-20 2019-06-27 주식회사 메디웨일 Method and apparatus for assisting in diagnosis of cardiovascular disease
CN111712186A (en) 2017-12-20 2020-09-25 医鲸股份有限公司 Method and device for assisting in the diagnosis of cardiovascular diseases

Also Published As

Publication number Publication date
JP2003299621A (en) 2003-10-21

Similar Documents

Publication Publication Date Title
JP6003487B2 (en) Blood pressure measuring device, blood pressure measuring method, blood pressure measuring program
RU2334455C2 (en) Device and method for outputting bioinformation and bioinformation report
TWI395572B (en) Blood pressure measuring apparatus and control method therefor
JP3882084B2 (en) Arteriosclerosis analysis system, arteriosclerosis analysis method, and arteriosclerosis analysis program
JP2005087753A (en) Beat onset detector
CN110709006B (en) Non-invasive blood pressure measurement
US20060074327A1 (en) Pulse wave information display apparatus, program product for controlling pulse wave information display apparatus, and method of displaying pulse wave information
JP6766710B2 (en) Blood pressure measuring device, method and program
JP5535477B2 (en) Vascular aging detection system
WO2006024871A1 (en) Methods and apparatus for the measurement of blood pressure
US6999812B2 (en) Arteriosclerosis detection system
US5365933A (en) Apparatus and method for distinguishing heart beats from intra-aortic balloon pump beats
JP3854888B2 (en) Atherosclerosis detection system
JP2008061824A (en) Medical measuring instrument, biosignal waveform extraction method and biosignal waveform extraction program
JP4813815B2 (en) Vascular sclerosis calculating device and vascular sclerosis calculating program
JP4417250B2 (en) Pulse wave detection system
CN107106046B (en) Blood vessel elastic modulus evaluation device
KR101002079B1 (en) Method of measuring blood vessel by blood vessel measurement apparatus
Freithaler et al. Smartphone-Based Blood Pressure Monitoring via the Oscillometric Finger Pressing Method: Analysis of Oscillation Width Variations Can Improve Diastolic Pressure Computation
JP5006509B2 (en) Pulse wave velocity measurement method for measuring pulse wave velocity in a pulse wave velocity measuring device
JP2002045361A (en) Measuring instrument for modulus of longitudinal elasticity of vital blood vessel tissue
WO2001030235A1 (en) A device for detecting arteriosclerosis, and method of use thereof
JP4606836B2 (en) Vascular sclerosis calculating device and vascular sclerosis calculating program
Angyal et al. Algorithm of oscillometric blood pressure measurement applied on microcontroller
JP2009273716A (en) Electronic sphygmomanometer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050404

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060323

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060516

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060718

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060911

R150 Certificate of patent or registration of utility model

Ref document number: 3854888

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090915

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090915

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120915

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130915

Year of fee payment: 7

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term