JP2016521162A - 循環監視システム - Google Patents
循環監視システム Download PDFInfo
- Publication number
- JP2016521162A JP2016521162A JP2016508967A JP2016508967A JP2016521162A JP 2016521162 A JP2016521162 A JP 2016521162A JP 2016508967 A JP2016508967 A JP 2016508967A JP 2016508967 A JP2016508967 A JP 2016508967A JP 2016521162 A JP2016521162 A JP 2016521162A
- Authority
- JP
- Japan
- Prior art keywords
- value
- sensor
- host computer
- signal
- calculated
- 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.)
- Pending
Links
- 230000004087 circulation Effects 0.000 title claims description 11
- 238000012544 monitoring process Methods 0.000 title description 4
- 238000013186 photoplethysmography Methods 0.000 claims abstract description 71
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000003745 diagnosis Methods 0.000 claims abstract description 16
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims abstract description 16
- 230000002093 peripheral effect Effects 0.000 claims abstract description 13
- 230000003595 spectral effect Effects 0.000 claims description 60
- 230000017531 blood circulation Effects 0.000 claims description 7
- 210000000056 organ Anatomy 0.000 claims description 7
- 238000004393 prognosis Methods 0.000 claims description 6
- 238000007477 logistic regression Methods 0.000 claims description 4
- 230000036558 skin tension Effects 0.000 claims description 2
- 206010009192 Circulatory collapse Diseases 0.000 claims 1
- 206010040560 shock Diseases 0.000 claims 1
- 208000030831 Peripheral arterial occlusive disease Diseases 0.000 abstract description 20
- 208000005764 Peripheral Arterial Disease Diseases 0.000 abstract description 19
- 238000004422 calculation algorithm Methods 0.000 abstract description 14
- 210000003371 toe Anatomy 0.000 abstract description 14
- 238000000034 method Methods 0.000 description 62
- 210000003414 extremity Anatomy 0.000 description 33
- 238000005259 measurement Methods 0.000 description 15
- 238000012545 processing Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 230000008859 change Effects 0.000 description 10
- 210000002683 foot Anatomy 0.000 description 9
- 238000010606 normalization Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 201000010099 disease Diseases 0.000 description 8
- 210000003491 skin Anatomy 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 208000035475 disorder Diseases 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 5
- 230000010412 perfusion Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 208000029078 coronary artery disease Diseases 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 description 4
- 208000024891 symptom Diseases 0.000 description 4
- 230000002792 vascular Effects 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 210000002615 epidermis Anatomy 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 201000001320 Atherosclerosis Diseases 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 2
- 208000015592 Involuntary movements Diseases 0.000 description 2
- 206010033425 Pain in extremity Diseases 0.000 description 2
- 208000006011 Stroke Diseases 0.000 description 2
- 208000005475 Vascular calcification Diseases 0.000 description 2
- 238000002266 amputation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000009556 duplex ultrasonography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000003811 finger Anatomy 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 210000004247 hand Anatomy 0.000 description 2
- 208000030175 lameness Diseases 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000004599 local-density approximation Methods 0.000 description 2
- 210000003141 lower extremity Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000017311 musculoskeletal movement, spinal reflex action Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004956 CI calculation Methods 0.000 description 1
- 208000017667 Chronic Disease Diseases 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 206010020565 Hyperaemia Diseases 0.000 description 1
- 208000033892 Hyperhomocysteinemia Diseases 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002730 additional effect Effects 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 210000002565 arteriole Anatomy 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 210000002302 brachial artery Anatomy 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 208000037998 chronic venous disease Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 210000002414 leg Anatomy 0.000 description 1
- 201000002818 limb ischemia Diseases 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 210000004279 orbit Anatomy 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000002106 pulse oximetry Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 208000022064 reactive hyperemia Diseases 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000002460 vibrational spectroscopy Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6825—Hand
- A61B5/6826—Finger
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6829—Foot or ankle
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
- G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Abstract
診断末梢動脈疾患に関する予知診断を提供するための末梢動脈流動障害検出システム。本システムは、ホストコンピュータと、たとえばPPGセンサを用いたフォトプレチスモグラフィにより信号を測定するなど、被験者の手指または足趾から生理信号を検出し測定するために用いられるセンサとを含む。センサデータは、信号波形に対応するいくつかの時間領域および周波数領域計算値を計算するために用いられる前に、処理されフィルタリングされる。この計算値は、予知診断を提供する多面的アルゴリズムを用いた予測モデルで用いられ、予知診断はモニタなどの表示器に表示される。
Description
本出願は、2007年12月11日に出願され、CIRCULATION MONITORING SYSTEM AND METHODと題され、2009年12月8日に米国特許第7,628,760号として現在発行されている、米国出願番号第12/001,505号から引き出された概念に依存しており、その内容および開示は、引用によりその全体が本明細書に明確に組み込まれている。
本発明は一般に、医療監視の分野に関する。より詳細には、本発明は、流動障害を特徴とする末梢動脈疾患に対する被験者の感受性を示すための循環監視および信号処理に関する。
末梢動脈疾患(PAD:Peripheral artery disease)、ならびに関連する冠動脈性心疾患(CHD:coronary heart disease)および頸動脈疾患(CVD:carotid vascular disease)は、死に至る可能性がある。
米国では、推計1000万人がPADを患っている。ほぼ同数が、症状がないため、および診断機器が比較的アクセスしにくいために、未診断であると考えられる。PADの症状のエンドポイントは深刻であり、すなわち身体障害、手足の切断、および死である。PAD、CAD、CVD、および他の合併症を患っている患者を疾患の早期段階で一次診療医が識別することを補助する、より簡単かつよりアクセスしやすいツールを提供することにより、疾患のより深刻な結果の多くを早期に介入および回避できるようにすることが望ましい。
PADは一般的に下肢粥状動脈硬化症に関連するが、CHD、CVD、心臓発作、脳卒中、および切断のリスク上昇に関連する。PADの患者の約75%は、CHDまたはCVDも患っている。脳卒中のリスクは、PADの患者では、その疾患がない患者の場合より3倍高い。PADは、通常は下肢の動脈の狭窄または閉塞となって現れ、粥状動脈硬化症、血栓症、動脈石灰化、糖尿病、およびホモシステイン血症を含むいくつかの要因により引き起こされる。PADは、ふくらはぎの痛みや障害、特に跛行、および重度肢虚血による制限歩行を特徴とする進行性の慢性疾患である。しかしながら、PADの全患者の約半数が、診断時に無症状であったことに留意されたい。
現在の診断方法は通常、跛行または安静時の下肢痛の症状が見られる患者に適用される。一般的な診断経路は、安静時または運動後の足関節上腕血圧比(ABI:Ankle−Brachial Index)の使用、反応性充血、フォトプレチスモグラフィ(PPG:photoplethysmography)、分節血圧分析、容積脈波記録(pulse volume recording)、複式超音波、および末梢血管造影法を含む。
ABIは、典型的には最初に行われる検査であり、通常は診療室または病院の血管検査室で行われる。ABIは、血圧計とドップラー超音波を用いて上腕動脈および足首で取得された収縮期血圧の観測値から計算される。ABIはPADの非侵襲的診断の標準と考えられているが、時間がかかり、導入しにくく、主観的であり、技術依存である。一貫した信頼性のある結果を得るためには、施術者は経験豊富で、専門的な訓練を受けていなければならない。さらにABIは、PADの危険性がある患者において一般的に遭遇する疾患である動脈石灰化を患っている患者には有益な診断ではない。この理由は、ABIが硬直した石灰化動脈を圧迫することに依存しているためである。そのような疾患は、偽陰性診断となることが多い。
従来のフォトプレチスモグラフィシステムは、被験者の組織の領域における血液の心臓リズム量を測定する。従来のパルスオキシメータは、被験者の組織の領域内の赤血球中のヘモグロビンと結合している酸素量を測定する。フォトプレチスモグラフィもパルスオキシメトリも、血流または血循環の質に対する直接的な相関関係を与えない。
Fischer Mら、「Simultaneous measurement of digital artery and skin perfusion pressure by the laser Doppler technique in healthy controls and patients with peripheral arterial occlusive disease」、Eur J Vasc Endovasc Surg.、1995年8月;10(2);231−6頁
Karel J.Zuzakら、「Visible spectroscopic imaging studies of normal and ischemic dermal tissue」、Proc.SPIE 3918、Biomedical Spectroscopy:Vibrational Spectroscopy and Other Novel Techniques,17(2000年5月8日)
Biomedix社(http://www.biomedix.com/products/PADnet_plus.asp)
Jane C GoldenおよびDaniel S Miles、「Assessment of Peripheral Hemodynamics Using Impedance Plethysmography」、PHYS THER.1986年;66:1544−1547頁
Myers,Kenneth、「The Investigation of Peripheral Arterial Disease By Strain Gauge Plethysmography」、ANGIOLOGY、1964年7月、15:293−304頁
本発明は一般的に、閉塞よりも末梢側の微小血管灌流の検知、および後続の信号処理に関する。より詳細には、ホストコンピューティングプラットフォームに動作可能に接続された、一体型の発光および検知部品を有する小型の手指または足趾の形のセンサに関する。ホストコンピュータは、プロセッサおよびメモリを含むものであるが、身体四肢の末梢位置における血流の質を評価する信号処理アルゴリズムを実施する。
検出システムは、被験者の手指または足趾にうまくなじむように輪郭形成され構成された筐体を含む。筐体内にセンサがあり、センサは生理信号を検出および測定することができる。生理信号の測定は、たとえば、プレチスモグラフィセンサを用いた動脈流の検出を含むことができる。測定される他の生理信号は、赤外光、器官容積、皮膚温度、皮膚インピーダンス、微小血管血流速度、または皮膚張力を含むことができる。これらの信号は、いくつかの異なる種類のセンサ、たとえばフォトダイオード、電荷結合素子、圧力カフ、電極または歪みゲージで測定することができる。
筐体およびセンサは、ホストコンピュータに動作可能に接続される。ホストコンピュータは、センサから波形としてデータを受信し、信号処理技法を用いて処理するように構成される。信号処理技法の例には、信号アーティファクト検出、正規化、信号フィルタリング、および他の時間領域および周波数領域技法が含まれる。
センサ信号波形から、信号波形を特徴付けるいくつかの時間領域特徴値および周波数領域特徴値が、ホストコンピュータにより計算される。計算される値の例には、循環指数、高調波の傾き、高調波切片、収縮期上昇、スペクトル信号、スペクトルノイズ、およびスペクトル信号対雑音比が含まれる。測定値は、被験者の各肢から、または被験者の他の様々な位置で取得される。
計算された値を用いて、ホストコンピュータは、計算された値を入力として用いて予知診断を生成する。予知診断は、使用されるセンサの種類および収集される信号の種類に特有の予測モデル式を用いて計算される。予測モデル式において用いられる特定の変数および係数は、既知の診断を有する被験者に対して行われた測定の標本群から得られたセンサデータにロジスティック回帰を行うことで生成される。
末梢動脈流検出システムの一実施形態において、このシステムは、ホストコンピュータと、被験者または患者から生理信号を収集するためのセンサとを含む。他の実施形態では、ホストコンピュータは代わりに、ノートブックコンピュータ、タブレット、モバイルスマートフォン、またはリモートサーバの形をとることもできる。
ホストコンピュータは、データを表示しユーザにシステムとインタラクトさせるソフトウェアアプリケーションを実行するように構成される。加えてシステムは、ホストコンピュータとインターフェースをとるように構成されるセンサを含む。インターフェースは、一実施形態では、USBケーブルを介した有線接続の形をとることができる。
他の実施形態では、センサは、任意数の知られている無線プロトコル、たとえばBlueeTooth(登録商標)またはIEEE802.11によって、ホストコンピュータとインターフェースをとることができる。ソフトウェアアプリケーションが、ホストコンピュータ上で、センサ自体の中で、またはインターネットを介したリモートサーバ上で実行できることを当業者であれば認めよう。
一実施形態では、生理信号は、センサを介して被験者または患者から収集される。一実施形態では、生理信号は、フォトプレチスモグラフィ(PPG)波形であって、PPGセンサを用いて収集されたものである。PPG波形はパルス波形の形であるが、他の生理信号を収集するための異なる種類のセンサを用いた場合に、必ずしもパルス波形となるわけではないことは理解されよう。
センサからのデータは、ホストコンピュータに転送されて処理される。ホストコンピュータは、センサから収集されたデータを表示する。ホストコンピュータは、PPG波形に基づいて、いくつかの値を計算することもできる。ホストコンピュータは、PPG波形に特有の値を用いて、予知診断を計算しモニタまたはディスプレイに結果を表示するように、ソフトウェアアプリケーションなどで構成される。
加えて、種々の異なるセンサ型を用いて種々の生理信号を収集でき、これには、限定はされないが、フォトプレチスモグラフィ(PPG)、レーザードップラー速度計測法(LDV)、赤外線電荷結合素子(CCD)、圧力変換器で構成された圧力カフ、皮膚インピーダンス電極、および皮膚歪みゲージが含まれることを当業者であれば理解しよう。
図1を参照すると、一実施形態では、ホストコンピュータ110はUSBケーブル125により筐体120に接続される。ホストコンピュータ110は、筐体120から収集されたデータを表示するソフトウェアアプリケーションを実行する。
筐体は、末梢肢の一部、たとえば手指または足趾を受けるように輪郭形成されている。筐体120内に配置されているのは、生理信号を検出し測定するためのセンサである。一実施形態では、センサは、動脈流を検出および測定し、パルス波形の形のフォトプレチスモグラフィ信号を生成するPPGセンサである。
前述のように、PPGセンサが、種々の生理信号を収集するための他種のセンサと代替可能であることは理解されよう。USBケーブル125として示されたインターフェースが機能を損なうことなく他の有線または無線データインターフェースと代替可能であることもまた理解されよう。
図1において筐体120は左手135の人差し指130に取り付けられているように図示されているが、筐体120は、他の身体付属肢に取り付けることができる。たとえば、図2を参照すると、筐体120は被験者の右足145の第二趾140に取り付けることができる。本明細書でさらに説明されるように、他の付属肢、器官、または身体上の位置、たとえば手指、足趾、足、手、脚、耳、額、眼窩上、および表皮などからのセンサ読取値を取得するのが望ましいことがある。
一実施形態では、ネットブックコンピュータを、ホストコンピュータとして使用することができる。たとえばホストコンピュータは、Microsoft Windows(登録商標)ベースのオペレーティングシステムなどを実行するように構成することができる。典型的にはそのようなコンピュータは、中央処理装置、たとえばAtom N455 CPUを含み、約1GBのランダムアクセスメモリ(RAM)を含むであろう。加えてそのようなコンピュータは典型的には、ネットワークカード、BlueToothカード、およびユニバーサルシリアルバス(USB)ポートを含む、外部デバイスに接続するための手段を含むであろう。さらにそのようなPCは、小型のLCDディスプレイを含むことがある。
代替のコンピューティングホストは、AppleのiOSを実行しているタブレット、Androidタブレット、またはWindowsタブレットを含むことができる。ホストコンピュータがより小型であれば、医療関係者が患者間で持ち運ぶのがより簡単なので、可搬性が促進される。
一実施形態では、ソフトウェアアプリケーションはC#で書くことができる。代わりに、いくつかの異なる言語が適することがあり、これには:Java(登録商標)、MatLab、C、C++、AJAX、Javascript、Perl、Ruby、Python、VBなどが含まれる。
どの生理信号が収集されるかに応じて、異なるセンサを使用できることは理解されよう。一実施形態では、センサは、手指または足趾を通る相対的な血流を測定するためのプレチスモグラフィセンサである。プレチスモグラフィセンサは、たとえばフォトプレチスモグラフィ(PPG)センサ、誘導性プレチスモグラフィセンサ、歪みゲージ式プレチスモグラフィセンサ、および容積プレチスモグラフィセンサを含む、プレチスモグラフィ用のセンサのクラスを指すことは理解されよう。一実施形態では、PPGセンサは、940nmで、同様の波長に感応するフォトダイオードとペアで動作する赤外発光ダイオード(IR LED)を含む。他の波長、たとえば630nmおよび535nmも使用可能であることは理解されよう。PPGセンサは、本明細書に記載の動作の多くを実施するようにプログラム可能なマイクロコントローラユニット(MCU)をさらに含むことができる。他の実施形態では、PPGセンサのMCUは、ホストコンピュータ自体の機能を提供することができる。適切なMCUの例は、テキサスインスツルメンツにより製造されたものであり、たとえばMSP430シリーズMCUである。他のコンポーネントをPPGセンサに含めることもでき、これには、リニアテクノロジ社のLTC2451機種などのアナログ/デジタル変換器(ADC)またはリニアテクノロジ社のLTC1669機種などのデジタル/アナログ変換器(DAC)が含まれる。
前述のPPGセンサを用いて、電圧出力、電流、またはADCビットカウントなどのデータ信号を提供することができる。そしてこのデータ信号は、ホストコンピュータによりさらに処理することができる。
一実施形態がPPGセンサを用いて記載されているが、血管生理信号を得るために他のセンサタイプを使用できることは当業者には理解されよう。たとえば他の研究者は、以下の文献に記載のように、血管灌流を監視するために様々なセンサタイプが有用であることを実証しており、それらの開示は、引用によりそれらの全体が本明細書に明確に組み込まれている:
レーザードップラー速度計測法:Fischer Mら、「Simultaneous measurement of digital artery and skin perfusion pressure by the laser Doppler technique in healthy controls and patients with peripheral arterial occlusive disease」、Eur J Vasc Endovasc Surg.、1995年8月;10(2):231−6頁。そのようなシステムは、レーザーを用いて表皮の局所領域を照射し、そして後方散乱レーザー光を測定する。この後方散乱光は、その速度に比例した周波数シフトを受ける。したがって、位相シフトに関連する生レーザー光、電圧、電流、またはADCカウントに関して、注目信号を生成し記録することができる。
赤外線CCD:Karel J.Zuzakら、「Visible spectroscopic imaging studies of normal and ischemic dermal tissue」、Proc.SPIE 3918、Biomedical Spectroscopy:Vibrational Spectroscopy and Other Novel Techniques,17(2000年5月8日)。赤外線CCDの場合、波長、電圧、電流、またはADCカウントの形で、信号を生成し記録することができる。
圧力カフおよび変換器:Biomedix社(http://www.biomedix.com/products/PADnet_plus.asp)。そのような圧力変換器は、圧力または圧力変化に関連する電流または電圧を生成することができる。これらの生の信号は、さらなる処理を可能にするためにさらにデジタル化することができる。
インピーダンスプレチスモグラフィ:Jane C GoldenおよびDaniel S Miles、「Assessment of Peripheral Hemodynamics Using Impedance Plethysmography」、PHYS THER.1986年;66:1544−1547頁。名称が示唆するように、皮膚インピーダンスの変化が、皮膚電極で測定される。出力信号は通常は電流であり、電流をADCによりデジタル化しホストコンピュータへ中継することができる。
歪みゲージプレチスモグラフィ:Myers,Kenneth、「The Investigation of Peripheral Arterial Disease By Strain Gauge Plethysmography」、ANGIOLOGY、1964年7月、15:293−304頁。このモダリティは、一般的には測定される器官を中心とした円周状に、表皮に貼付られた歪みゲージを介して実現される。したがって、動脈灌流は器官における容積変化をもたらし、これは歪みゲージの抵抗を変化させる傾向がある。このアナログ電圧を、デジタル化しホストコンピュータへ中継して処理することができる。
上述の技術的モダリティの各々において、動脈灌流は検知され、さらなる処理のための信号としてホストコンピュータへ中継される。
多面的アルゴリズム
多面的アルゴリズム(MFA:multi−faceted algorithm)の動作を以下説明することにする。被験者または患者の異なる肢からPPGセンサから収集されたPPG波形との関連で説明を行うが、前述のような他種のセンサを用いて様々な体の位置において収集される生理信号データを用いて予知診断を分析し導出するために、MFAの方法を使用できることは理解されよう。
フォトプレチスモグラフィ(PPG)のようなセンサから得られる信号トレースを分析して、形態特徴を得ることができる。PPG波形の形態特徴は、肢の疾患状態を示す。PPG波形の形態特徴は、時間領域特徴、たとえば収縮期上昇および正規化された信号を含むことができる。PPG波形の形態特徴は、周波数領域特徴、たとえば高調波の傾き、高調波切片、スペクトル信号、およびスペクトル信号対雑音比をさらに含むことができる。様々な指標をPPG波形から計算することもでき、たとえば、米国特許第7,628,760号(その開示は、引用によりその全体が組み込まれている)に記載の循環指数(CI:Circulation Index)、および本明細書に記載の拍動指数(PI:Pulsatility Index)がある。さらなる計算値を、本明細書に記載のように、これらの特徴の各々から計算することができる。それらの計算の結果は、PPG波形を特徴付けるものであるが、MFAにおける入力成分として用いて、動脈流動障害の確率を予測することができる。
MFAは、PPG波形のスペクトル平坦性だけでなく、PPG波形の他の周波数領域および時間領域特徴も考慮に入れる。PPG波形の形態特徴から生成された計算値は、MFAへの入力として用いられ、この入力を用いて、CIまたはPI単独よりも信頼性高く流動障害の確率を予測することができる。
PPG波形の形態特徴に基づいて生成可能な計算値の種類の例は、以下で説明される。
肢/指ごとの測定値および計算値:
各肢または各指からセンサを介してデータを収集することができる。データは、センサから中継され、測定終了時またはリアルタイムに分析される。一実施形態では、血液量測定値が、約50回毎秒の速度でサンプリングされる。一実施形態では、指の測定は、有限時間、たとえば15秒とすることができ、または、ユーザが終了を決定する連続測定とすることができる。
任意選択で、リアルタイムに、周期的に(たとえば毎秒1回)、または指測定の終了時に、測定値は、信号アーティファクトの存在に関して評価される。様々なアーティファクト検出方法を用いることができる。そのような方法には、バイスペクトル/バイコヒーレンス分析、DC不連続性評価、尖度、歪度が含まれる。一実施形態では、標本尖度法が用いられる。標本尖度法は、式1:
で表される。
信号アーティファクトを有さない測定値は、低尖度尺度、たとえば1.0未満で記録される。一実施形態では、1.0の尖度閾値により、臨床的に意義のあるアーティファクトの検出において97%の感度および97%の特異度を達成することができる。
センサデータは吸収光ではなく反射光を表すので、減算は、吸収光を表すように信号データを変換する単純な手段であり、吸収光は、撮像された細動脈床の容積変化に直接的に関連する。加えて、さらなる時間領域計算は、PPG波形の正常な時間的な定位(temporal orientation)に依存する。加えて、血液量変化の直接的尺度である信号を生成する正規化ステップを行うことができる。以下の形(式2)のような単一の数式で信号を減算および正規化することは、計算機的に効率的かつ簡便である:
ここでrawDataはセンサからのデータであり、initialMedianは初期の生データの、通常は3点の中央値である。初期点ではなく初期中央値を用いて正規化し、それによって初期の過渡状態または他の異常を許容する。他の正規化技法、たとえば(最初のn点の)初期平均、m点(たとえば1秒)の移動平均、全体平均、または固定公称値、たとえば公称レベルに関するデータ(たとえば、16ビットADCからの32768ADC点)を使用できることは認識されたい。
加えて他の減算技法、たとえば(最初のn点の)初期中央値、(最初のm点の)初期平均、移動平均(たとえば1秒)、全体平均、固定公称値(たとえば16ビットADCからの32768ADC点)または固定最大値(たとえば16ビットADCからの65536ADC点)を使用することができる。
平均信号レベル(DC)を表す量により信号を正規化することで、ピーク対トラフ(peak−to−trough)(AC/DC成分)が生成され、これは臨床的に意義があり得る。DC値により別途正規化されなければ、照明源からの平均光量を二倍することは、生のAC成分を二倍することを伴うことになり、これは、初期中央値により提供されるような動的正規化技法の適用が必要なことを示していることを理解されたい。
正規化および減算の後、信号の正規化されたスペクトル密度が、ウェルチ(Welch)法の修正版を用いて推定される。具体的には、各スペクトル密度の正規化は、スペクトル密度推定のために各周波数のパワーを最大値で割ることで行われる。この技法はノイズの影響を低減し、ノイズの影響は、生理信号を絶対値パワーにおいて圧倒する傾向があるが、一過性のために、ウェルチ法に関する平均化およびこの方法への正規化の追加により緩和される。
たとえば、生理信号が、1Hzおよび5000パワー単位で発生し、ウェルチ法の各高速フーリエ変換(FFT)を通して比較的一貫して現れると仮定する。しかしながら、多数のFFTのうちのほんの一部に0.5Hzの周波数において10000単位のパワーを有するノイズが存在する場合、正規化によりそのパワーが1.0単位に変換される。加えて、ノイズは、ウェルチ法で用いられるFFTのごく一部のみに存在することがある。ノイズ軽減の性質のために、スペクトル密度を正規化することで、基本周波数(F0)のより信頼性の高い推定が可能となる。
F0を推定するためのいくつかの方法が存在し、これには:突出(prominence)、単純な最大値、高調波積スペクトル(HPS)およびその変形が含まれる。ノイズを持たない強力な信号の場合、スペクトル密度に関する最大パワーを単に選択するなどの単純な技法が可能になる。そのような選択はさらに強化することができ、周波数範囲をたとえば0.6から2.5Hzまでに狭めることで計算効率が向上する。
生理信号がノイズと混同される場合、F0をより良く推定するためには、より高度な技法が必要となる。HPSは1つのそのような方法であり、基本的に、所与の周波数およびその関連する高調波のパワーの積に関連する。たとえば、1.2Hzのパワーに、2.4Hz、3.6Hz、4.8Hzのパワーを乗算するなどである。この方法の変形例は、積の代わりに総和を含む。
特に有用な他の方法は、周囲の点に対するスペクトルのピークの突出、すなわち相対的な高さの使用を含む。たとえば、1.2Hzの突出は、0.8Hzおよび1.6Hzにおけるパワーの平均に対する1.2Hzにおけるパワーの比率として計算することができる。突出は、低周波においてパワーが高いが周波数が増加するにつれて減衰する低周波ノイズを含むスペクトルにおいて特に有利である。そのような場合、HPSは、突出法と同様に局所ピークを区別しない。隣接比較の幅(中心周波数からのスペクトル「距離」)は、突出法において重要なパラメータである。
分母比較の1つの選択肢はF0/2であり、これは理論的には、F0ピークに対するトラフに対応するはずである。ノイズが多い臨床データを用いた実証試験は、幅35%が最適であることを示していた。たとえば、推定されたF0が1.2Hzの場合、隣接比較周波数は0.78Hzおよび1.62Hzであった。以下の式3は、突出の計算を示す:
ここでP(f)は所与の周波数におけるパワーであり、λはオフセット幅である。F0の選択は、最大の突出に関連する周波数に基づく。
突出およびHPSの組み合わせが、本発明の範囲内であると見なされるべきであることを理解されたい。
F0の推定が行われると、信号をフィルタリングして、注目信号の変更を最小限にしながら低周波成分を除去することができる。すなわち、F0を知ることで、F0を最小限に減衰させるがF0未満の周波数を最大限に減衰させるコーナー周波数を有するハイパスフィルタを実施することができる。
ハイパスフィルタを実施することで、基本的に、信号がトレンド除去され、注目信号に対するより信頼性の高い時間領域および周波数領域処理が可能となる。いくつかの異なる方法が、時間領域および周波数領域の両方におけるハイパスフィルタリングのために使用可能であり、これには:バターワース、チェビシェフ、移動平均、および畳み込みが含まれる。
たとえば、4次の双方向バターワースフィルタにより、急峻な減衰、零位相シフトが可能となり、合理的に計算上効率的かつ安定である。加えて、帯域通過範囲において平坦な周波数応答を提供し、これは、比較的純粋な信号に依存するさらなる周波数領域処理のために特に重要である。F0の約40%のコーナー周波数を選択すると、5%前後にF0の減衰が制限される傾向があるが、たとえば:不随意運動、周囲光の変化、および呼吸などに起因する成分が効果的に減衰する。
米国特許第7,628,760号に記載のように、循環指数(CI)は、信号の準周期成分の強度の尺度であり、基本的にスペクトル平坦性尺度(SFM:Spectral Flatness Measure)を補うものである。SFMのように、CIは0から1までの大きさであるが、パーセンテージ(たとえば0%から100%まで)として表されることもある。SFMは音声処理および他のアプリケーションで使用されているが、米国特許第7,628,760号に記載のような限定された状況を除いて、SFMは心臓血管信号アプリケーションで使用されていない。
CIの計算は、様々なトレンド除去方法により信号の低周波成分を除去することで改善することができ、これらの方法には:ハイパスフィルタリング、線形回帰、多項式フィッティング、キュービックスプラインフィッティング、または平滑化事前分布が含まれる。通常は生理学的プロセスおよび/または環境条件に関連する低周波成分の除去は、細動脈循環の容積変化に関連しない。すなわち、呼吸、不随意運動、および環境光変化は、前処理ステップとして信号から基本的に除去することができ、CIのより臨床的に重要な推定値を生成することができる。
スペクトル平坦性尺度(SFM)は通常、PSDの全周波数範囲にわたって定義されるが、任意の周波数帯域に適用することもできる。具体的には、一実施形態は、0.6Hzから8.0Hzまでの周波数範囲を含み、この範囲にはほとんどの患者の基本周波数および対応する3つの高調波が含まれる。
一代替実施形態は、CIの代わりに、またはCIに加えて、拍動指数(PI)を用いる。前述および米国特許第7,628,760号に記載のように、CIは、陽性PAD診断の確率を予測するための独立メトリックに関連付けられ、独立メトリックとして使用することができる。
一部の患者において、PPG波形のスペクトル平坦性単独では、PAD診断の信頼性の高い予測を与えることができない。近位流動障害の患者では、PPG波形は、高レベルの低周波ノイズと結合された低生理信号振幅を示すことがある。患者の呼吸、不整脈、低周波の血流変化、または他の周囲光および運動アーティファクトもまた、PPG波形のスペクトル平坦性に影響することがある。
被験者の左足から取得されたそのようなPPG波形の一例が図3に示される。PPG波形210が、時間にわたる振幅としてプロットされている。観察され得るように、PPG波形210は、高レベルの低周波ノイズに結合された低生理信号振幅を示している。このことは、PPG波形210に関する対応するパワースペクトル密度図を見るとより明白となる。
図4の対応するスペクトル密度図から分かるように、波形は、低周波ピーク225の突出を特徴とする。この低周波パワーは、高レベルのノイズが存在することを示す。これらの種類の患者では、米国特許第7,628,760号に記載の循環指数(CI)などの方法を用いてスペクトル平坦性尺度単体からのみ計算される予測は、歪んだ結果を与えることがある。この理由は、低周波ノイズによりスペクトル平坦性尺度が低下し、ひいてはCIが異常に高くなるためである。
流動障害の唯一の予測因子としてのスペクトル平坦性に基づくCIのそのような制限を克服するために、代替のアルゴリズムを使用することができる。この代替のアルゴリズムは、拍動指数(PI)と呼ばれるが、以下説明することにする。
拍動指数(PI)は、基本周波数およびその高調波におけるスペクトルパワーを、それらの間のトラフにおけるパワーとの比較で考慮して、非生理信号ノイズを無視することで、この問題に対処する。
図5に、拍動指数を決定するための1つのそのようなアルゴリズムの主な処理ステップを示す。310で示される信号x(n)で開始して、拍動指数を決定する際の最初のステップは、ステップ320に示されるように信号をトレンド除去することである。これは、波形のトレンドを推定し除去することで実現することができる。トレンドは任意数の方法を用いて推定することができ、これには、限定はされないが:平滑化事前分布法、線形回帰、多項式フィッティング、キュービックスプラインの使用、またはハイパスまたはローパスフィルタの適用が含まれる。一実施形態では、波形のトレンドは、式4に示すように、トレンドを推定するための有限インパルス応答(FIR)ローパスフィルタを用いて推定される:
ここでMはインパルス応答h(n)の長さであり、nは離散時間インデックスであり、x(k)は入力信号であり、xt(n)は出力である。一実施形態では、インパルス応答は等しい重みを有する。他の実施形態では、Epanechnikovカーネルが、インパルス応答として用いられる。
カーネルの幅は、カットオフ周波数の見地から指定される。一実施形態では、アルゴリズムは、0.5Hzのカットオフ周波数を用いる。0.5Hzで、カットオフ周波数は、トレンドの推定において心臓周期を含まないようにするのに十分低いが、低周波ドリフトに影響する呼吸および他の要素の付加的な効果を含むのに十分高い。他の実施形態は、F0の推定値に基づいて動的なコーナー周波数を用いることができる。高速フーリエ変換(FFT)アルゴリズムを用いて、アルゴリズムの計算コストを減らすことができる。一実装では、最初および最後の部分(たとえば0.5秒)は、フィルタのエッジ効果を除去するために、切り捨てられる。
信号がトレンド除去された後、ステップ330に示されるように信号のスペクトル密度が推定される。任意数の方法を用いてスペクトル密度を推定することができ、これには、限定はされないが:パワースペクトル密度(PSD)、ウェルチ法、およびバーグ(Burg)法が含まれる。一実施形態では、スペクトル密度推定は、パワースペクトル密度(PSD)である。PSDは、パラメトリック法またはノンパラメトリック法を用いて推定することができる。一実装では、PSDは、ブラックマン−テューキースペクトル推定法を用いて推定される。信号の自己相関が、式7に示されるように、標準バイアス推定量(standard,biased estimator)を用いて推定される:
ここでNはトレンド除去されたセグメントの長さであり、セグメントのインデックスは0からN−1までであり、lは推定値が計算される遅延を示すインデックスである。
そして、自己相関にウィンドウwa(l)が乗算され、ウィンドウ化された自己相関関数の離散時間フーリエ変換が計算される。任意数のウィンドウ化技法を用いることができ、これには、限定はされないが、ブラックマン、ハミング、ハニングなどが含まれる。一実施形態では、式8に示すように、ブラックマンが用いられる:
ここでωは周波数インデックスである。心臓周期によるPSD内のピークを推定するのに十分な持続時間が存在するが、推定値の分散を大幅に低減しスプリアスピークを排除するのに十分な平滑化も与えられるように、ウィンドウが選択されるべきである。一実装では、ブラックマン自己相関ウィンドウは、8秒の持続時間を有する。
PSDが推定された後、基準を最適化する心拍数をあり得る心拍数周波数の範囲にわたって検索することで、ステップ340に示されるように心拍数を推定することができる。一実施形態では、0.7Hzから2.5Hzまでの周波数範囲で検索が行われ、この範囲は、42から150拍毎分に相当する。基準は、高調波の間のトラフにおけるPSDのパワーに対する候補心拍数の高調波におけるパワーの総和に基づく。この方法は、前述の突出法の変形である。高調波のパワーは、式9に示されるように、ピークの整数倍におけるPSDの値として算出される:
ここでNhは使用される高調波の数である。一実施形態は、基本波(F0)および3つの高調波(F1、F2、F3)を用いる。高調波信号のパワーは、他の方法、たとえば、近接する周囲の隣接周波数を有する各ピークの面積を計算することなどで推定することができる。
トラフのパワーは、高調波の間の周波数におけるPSDの値として計算される。一実装では、高調波の各ペアの間の3つの周波数が用いられ、具体的には、式10に示されるように、中点(50%ピーク間点)、40%ピーク間点、および60%ピーク間点である:
心拍数が推定された後、拍動指数を、式13に示されるようにステップ350において計算することができる:
ここでρ(ωhr)は、推定心拍数で評価された、より前段で記述した比である。拍動指数は、図5の参照文字b360で表される。
ピーク周波数における高いパワーに対応する小さい比は、大きい拍動推定値となる。トラフにおけるより大きいパワー、およびピークにおけるより小さいパワーに対応する大きい比は、小さい拍動推定値に対応する結果を返す。
具体的に、MFAに戻ると、トレンド除去により低周波成分が除去された後、低周波ノイズに支配されない従来のスペクトル密度推定を行うことができる。また、以前に論じられたような様々な方法を用いることができ、これには従来のウェルチ法が含まれる。このスペクトル密度推定は、F0をより正確に推定するのに特に有益であり得る。典型的には、0.6Hzから2.5Hzまでのような特定の周波数範囲の間のピークパワーは、F0に対応し得る。また、F0を推定する様々な方法を用いることができる。一実施形態では、突出法が用いられる。
スペクトル密度推定により、スペクトル信号、スペクトルノイズ、スペクトル信号対雑音比(SNR)などの変数の計算のベースが与えられる。図6に、PPG波形405のスペクトル密度推定のスペクトルパワー図400を示す。スペクトルパワー図400から、基本周波数410(F0)、ならびに第1高調波415(F1)、第2高調波420(F2)、および第3高調波425(F3)を、PPG波形405の推定スペクトル密度から識別することができる。
スペクトル信号(「spec.sig」)およびスペクトルノイズ(「spec.noise」)は以下から計算することができる:
ここでPiはF0、F1、F2、...Fnにおけるパワーであり、PjはF0/2、3F0/2、5F0/2、...(N+1)F0/2におけるパワーであり、これは通常は図6の「x」マーク412に対応する。
図7を参照すると、PPG波形455のスペクトル密度推定のスペクトルパワー図450の他の例が示されている。また、基本周波数460(F0)、ならびに第1高調波465(F1)、第2高調波470(F2)、および第3高調波475(F3)は、PPG波形455の推定スペクトル密度から識別することができる。
基本周波数および高調波周波数を識別した場合、高調波減衰(HD:Harmonic Decay)を、式17を用いて定義することができる:
ここでAは高調波切片であり、bは高調波の傾きであり、fは高調波番号(0、1、2、3など)または実際の対応周波数(1.2、2.4、3.6、4.8Hzなど)である。高調波番号または実際の周波数を、高調波減衰変数の計算で用いることができる。いくつかの例では、高調波番号は実際に、絶対周波数よりも統計的に重要な変数を与える。
高調波の傾き(「harm.slope」)の計算は、標準的な最小二乗法を用いるか、または市販のソフトウェアパッケージであるMicrosoft Excelで使用されるLINEST関数およびINDEX関数を用いて、以下の公式で行うことができる:
ここでPF0は基本周波数のパワーであり、PFiはF0の第i高調波のパワー、すなわち第1高調波(PF1)、第2高調波(PF2)などである。任意数の高調波を用いることができるが、十分な臨床的有用性は通常、第3高調波(i=3)までを用いて達成することができる。上記の関数は、高調波のパワーがF0のパワーで正規化されること(たとえばPFi/PF0)を示す。一代替実施形態は、この正規化を除外する。
また、典型的な最小二乗関数を、上記のMicrosoft Excel関数の代わりに用いることができる。同様に、高調波切片を、F0のパワーで正規化してもよいし、しなくてもよい。
正確な高調波減衰変数に加えて、トレンド除去された信号もまた、F0の正確な推定値のベースとなり、これはさらなる時間領域および周波数領域計算において重要である。また、いくつかの方法を用いて、前述のようにF0を推定することができる。一実装は、式3に記載の突出法を用いる。
信頼性の高いF0が推定されると、信号をさらにフィルタリングして、注目信号より高い側にある高周波成分を除去することができる。また、いくつかの方法を用いて、時間領域および周波数領域の両方において高周波成分を除去することができ、これには:バターワース、チェビシェフ、移動平均、および畳み込みが含まれる。一実装では、たとえば、F0の4.5倍のコーナー周波数を有する4次の双方向バターワースローパスフィルタを用いる。たとえば、F0=1.2Hzの場合、コーナー周波数は5.4Hzに設定されることになる。このコーナー周波数の選択を選択することで、F0と第3高調波に続く「トラフ」との間で信号が変わらない。一般的には、コーナー周波数の選択は、以下の式に基づき得る:
ここでHNは高調波番号(たとえば、3)である。
より高い周波数成分の減衰により、注目する周波数成分のみを含む信号が得られる。高周波ノイズを除去することで、たとえば、ピークおよびトラフ検出、または単に「ピーク検出」のより信頼性の高い手段が可能となる。
周波数領域特徴に加えて、PPG波形は、いくつかの時間領域特徴を有することになる。
頑健なピーク検出は、収縮期上昇(SR:Systolic Rise)変数を信頼性高く計算するために重要である。ノイズの多い信号には積極的なフィルタリングが有利であるが、フィルタ仕様は、信号の形状を著しく変更しないように注意して選択すべきである。したがって、F0の0.40倍から4.5倍までの周波数範囲を選択することで、合理的にフィルタリングされた信号が与えられ、これにより注目信号を歪ませることのない頑健なピーク検出が可能になる。
PPG波形のSR特徴を、図8を参照して論じることにする。波形510のSRは、周期520(P)のパーセンテージとして表すことができ、次のように計算される:
ここでSは波形がトラフからピークへ上昇するのに要する時間525であり、Pは信号の周期である。あるいは、SRは絶対的に表してもよい。
収縮期上昇は測定肢におけるPADと相関があり、疾患肢はSR周期がより長い。SRの適切な決定は、波形内の適切なピークおよびトラフの決定に左右される。一実装では、ピークの検出は、まず、最初の候補SRを最大疑似傾き
に基づいて推定することで行われ、ここでdyはパルス振幅の変化でありdxは時間の変化であり、60ミリ秒から推定周期(P0)までの様々なラン距離(run distance)dxを通して行われ、推定周期は、推定されたF0から導出される。そして疑似傾き配列が生成され、疑似傾き配列は、候補SRに基づく固定ラン(fixed run)(dx)を有する
を表す。疑似傾き配列は、PPGのSR特徴の中心を推定するために用いられる。ピークおよびトラフを発見するのに用いられるのはこの特徴(収縮期の傾き)であり、単一点ではない。
そして、疑似傾き配列の最大値として定義された中心インデックスを用いてPPG波形を見直すことで、ピークおよびトラフが位置特定される。まず、
の第1の極大値が位置特定され、これは第1の収縮期傾き特徴に対応する。そして、dxのインデックスを中心インデックスから開始して、引数
を最大化することで、トラフが識別される。最小値が発見され確認されると、対応する最大値が探索される。最小値および最大値のペアは、トラフおよびピークにそれぞれ対応し、ピーク配列に追加される。そして、次の中心インデックスが、1周期分だけ増加される。そして、その中心インデックスが、±P0/2内(半周期内)で傾き配列の最大傾きと一致するように調整される。この処理は、PPG波形全体にわたって反復されて、PPG波形内の周期ごとにトラフとピークのペアを識別するピーク配列が生成される。
まとめると、ピークは以下の基準に従って識別することができる:
・ペアで発生する
・P0の距離で名目上発生する
・信号の最も急峻で最長の上昇部を囲む(収縮期上昇周期525、M2=dy2/dx)
・トラフは:M2の中心からのmaxargx(dy2/dx)により識別される
・トラフは信号配列の端には無い
・ピークはトラフの後に発生する
・ピーク>周囲の2点
・トラフおよびピークの間の距離(d)は:15%P0≦d≦50%P0
・ペアで発生する
・P0の距離で名目上発生する
・信号の最も急峻で最長の上昇部を囲む(収縮期上昇周期525、M2=dy2/dx)
・トラフは:M2の中心からのmaxargx(dy2/dx)により識別される
・トラフは信号配列の端には無い
・ピークはトラフの後に発生する
・ピーク>周囲の2点
・トラフおよびピークの間の距離(d)は:15%P0≦d≦50%P0
比較計算
さらなる診断精度は、後続の肢または指の測定および変数計算により達成することができる。たとえば、足趾から手指へ、または左足の足趾から右足の足趾へ行われた同じ計算を比較することで、疾患を臨床的に示す相対変数を生成することができる。一例として、足趾の循環指数(CI)を、様々な方法で手と比べることができ、たとえば比(たとえばCItoe/CIfinger)または差(たとえばCItoe−CIfinger)などである。さらに、足趾対手指変数は、種々のベース(たとえば除数または減数)を用いることができ、たとえば:平均、「最良(best)」、最大、または最小(たとえばCIleft−toe/CIAVG(fingers)またはCIleft−toe/CIMAX(fingers))などである。加えて、所与の足趾を、反対側の足趾、または足趾の平均と比較することができる(たとえばCIleft−toe/CIAVG(toes))。本明細書に記載の各改変例が、各肢ごとの変数(たとえばCI、PI、高調波の傾きなど)に適用可能であることは認識されたい。
「最良」ベースは単に、手、器官の位置、または一般の肢であり、これは特定の質に基づいて選択される。たとえば、最良の手は、2つの手のうちの最大のSNRまたはCIに基づいて選択されることがある。したがって、その手は、各比較計算のベースとして用いられる。
したがって、いくつかの変数を用いて、相対肢変数(たとえばCIleft−toe/CIAVG(fingers))に加えて絶対肢変数(たとえばCIleft−toe)を用いることで、所与の肢の流動障害の確率を推定することができる。
予測モデル
流動障害に対する予測モデルは、絶対変数および相対変数を用いて生成することができる。一実装は、以下の形のロジスティック関数を用いる:
ここでP(Dx)は流動障害の確率であり、v1、v2、...vnは上記で議論された計算に基づく絶対変数および相対変数であり、C0、C1、C2、...Cnは各変数に対応する係数である。
センサの種類または機種に特有の予測モデルを決定するために、既知の診断を有する被験者の標本集団から、センサデータがたとえば各肢から採取される。上述の技法を用いて、各肢からの信号が処理されて、上述のセンサデータを特徴付ける肢ごとの変数および比較変数が計算される。通常肢および異常肢測定値からの集約データに対して行われるロジスティック回帰モデルを用いて、各変数に対する係数が決定される。
いくつかの変数に対する係数が零であることがあり、これは変数が統計的に有意でなく省略可能であることを意味することに留意されたい。式22は0.0から1.0までの範囲の値を生成し、特定の閾値(たとえば0.5)未満の値は、肢が流動障害であることを示唆する。
上記の閾値は、感度、特異度、および精度を考慮して、感度分析に基づいて決定することができる。一実施形態では、感度が、特異度より有利であることがある。他の実施形態では、閾値を、精度の最大化に基づいて設定することができる。
係数および閾値の値が、当業者であれば上述の方法を用いて容易に決定できることは理解されよう。
一実施形態では、センサはメモリを含むことができ、関連する係数がメモリへロードされ、ホストコンピュータへダウンロードされるようにメモリは構成される。予測モデルに対して計算される関連する係数および変数は、ホストコンピュータに接続されたセンサの種類または機種に応じて異なる場合がある。他の実施形態では、関連する係数はメモリへロードされ、ホストコンピュータに挿入されたメモリカードを介してホストコンピュータにダウンロードされる。さらに他の実施形態では、関連する係数はメモリにロードされ、USBドライブを介してホストコンピュータにダウンロードされる。他の実施形態では、係数は、リモートサーバからホストコンピュータへロードされる。
図9Aおよび図9Bに、4つの機能肢、すなわち右腕、左腕、右足、および左足を有する患者に対する一例示的実装における多面的アルゴリズムの動作ステップを説明する処理フロー図を示す。この例では、ホストコンピュータに接続されたPPGセンサは、PPGセンサにデータを取得させるために、各肢に取り付けられる。
各肢に対して、PPGセンサはステップ710のようにデータ信号を取得する。次いで、ホストコンピュータ、あるいはセンサ自体は、ステップ712に示されるように、任意の信号アーティファクトを識別する。次いで、ステップ714に示されるように、信号は減算され正規化される。ステップ716に示されるように、正規化されたスペクトル密度が、取得されたセンサ信号に対して推定され、基本周波数および高調波がステップ718において識別される。ステップ720において信号はハイパスフィルタで処理されて、低周波ノイズ成分が除去される。
低周波ノイズが除去された後、異なる変数を計算するために、信号をいくつかの異なる方法で処理することができる。1つの経路では、循環指数(CI)がステップ730で計算されて、ステップ732で評価される特定の肢に対するCIが生成される。
図9Aおよび図9Bには明示されないが、CIを計算する代わりに、またはこれに加えて、PIを各々に対して計算できることは理解されよう。
ステップ740のように、信号を処理して、PPG波形の絶対スペクトル密度を推定することもできる。このステップから、スペクトル信号、スペクトルノイズ、およびスペクトル信号対雑音比を、ステップ742に示されるように計算することができる。同じ情報から、高調波の傾きおよび高調波切片を、ステップ744に示されるように計算することもできる。
ステップ740から、スペクトル密度の基本周波数を、ステップ750において決定することができる。そして、ステップ760において信号をローパスフィルタを通過させることができ、これから信号のAC振幅をステップ762において計算することができる。ローパスフィルタステップは、信号のより高い周波数成分を減衰させる働きをして、ステップ770の信頼性の高いピークおよびトラフ検出を容易にする。
ピークおよびトラフがステップ770において検出されると、収縮期上昇周期をステップ772において前述のように計算することができる。
この処理は、患者の四肢の各々に対して反復される。四肢全てからのデータが収集されると、ホストコンピュータは、ステップ780、782、784、786、788、および790に示されたような比較変数を計算する。
比較変数の値を計算した場合、ホストコンピュータは、計算された比較変数および所定の係数の組を用いて、予知診断を計算する。予知診断は、ホストコンピュータに接続または統合されたモニタなどの表示器(indicator)上に表示することができる。
本発明の趣旨および範囲内で、特定の例示の機能ブロックが省略、並べ替え、結合、または分離可能であることを当業者であれば理解しよう。同様に、本発明の趣旨および範囲内で、特定の例示のソフトウェアステップが省略、並べ替え、結合、または分離可能であることを当業者であれば理解しよう。図9Aおよび図9Bに示された処理フローに対する全てのそのような適切なトポロジー的および論理的に適切な代替物は、本発明の趣旨および範囲内であると考えられる。さらに、そのようなステップは、ソフトウェアにより実装されているように示されているが、ハードウェアおよび/またはファームウェアにより適切に実装することができる。
上記のステップはPPGセンサに特有の用語で記載されているが、異なるセンサ技術を使用できることも当業者であれば理解しよう。加えて、上記の例は手および足、または手指および足趾の使用を示しているが、他の身体器官を同様に調べることができる。
実験
PPGセンサを用いた一実験では、実験の目的は、PPG検知デバイスからのデータを用いて末梢動脈流動障害(FO)を予測するアルゴリズムを開発することであった。
PPGセンサからの生データが、70人の患者から収集された。55人の患者は、末梢動脈疾患の疑いがある兆候または症状と共に、5つの血管センターを訪れた。各患者は、各人差し指および各第二趾の両方においてPPGデバイスで検査され、続いて、実績のある画像モダリティ(たとえば複式超音波、血管造影法など)で評価された。加えて、15人の正常対照被験者も同様に、ただし追加の画像モダリティなしで測定された。
次いで、PPGセンサからの生データファイルがソフトウェアアプリケーションで分析され、このソフトウェアアプリケーションは:循環指数(CI)、収縮期上昇(SR)、高調波の傾き(HS)、高調波切片(HI)、スペクトル信号(SS)、スペクトルノイズ(SN)、スペクトル信号対雑音比(SNR)を含むいくつかの異なる変数の計算を行った。加えて相対変数が、足趾に対して、手指と比較して計算された。FOを患っている肢は1の診断値が割り当てられ、FOが無い肢は0が割り当てられた。次いでこのデータは、XLStat(バージョン2012.5.02)により、モンテカルロシミュレーションにおいてロジスティック回帰分析を用いて評価された。すなわち、107肢がランダムに選択され、ロジスティックモデルが作成され、そして他の30肢に対して検証された。
70人の被験者および137肢から、74肢は流動障害があり、63肢は流動障害が存在しなかった。重要な変数には、式24にリスト化されたもの:CI(足/手最大値)、高調波の傾き(足/手最大値)、高調波切片、高調波切片(足−手最大値)、スペクトル信号(足/手最大値)、および収縮期上昇が含まれていた。FOの予測における中央値精度は、89.7%であった。
本発明が、本明細書で説明および例示された方法、または、構造、製造、材料、用途もしくは使用の詳細に限定されないことは理解されよう。実際に、製造、使用または用途の任意の適切な変形例は、一代替実施形態と考えられ、したがって、本発明の趣旨および範囲内である。
使用もしくは動作の用途もしくは方法、構成、製造の方法、形状、大きさ、または原料の任意の変更によって得られる本発明の任意の他の実施形態は、本明細書に含まれる詳細な記載説明または例示では示されないが当業者には理解されるはずのものであり、本発明の範囲内であることがさらに意図される。
最後に、本明細書で説明され例示された発明の方法、システムおよび装置が、ソフトウェア、ファームウェアまたはハードウェア、あるいはそれらの任意の適切な組み合わせで実装できることを当業者であれば理解しよう。低コストおよび柔軟性の目的で、この3つの組み合わせで、方法および装置を実装することが好ましい。したがって、本発明の方法およびシステムの実施形態が、命令が内部で実行されるコンピュータまたはマイクロプロセッサ処理により実装可能であり、その命令がコンピュータ可読媒体上に実行用に格納され、任意の適切な命令プロセッサにより実行されることを当業者であれば理解しよう。
したがって、本発明は、発明された装置の前述の実施形態を参照して示され説明されたが、添付の特許請求の範囲に定義された本発明の趣旨および範囲から逸脱することなく形式および細部の他の変更がなされ得ることは当業者にとって明らかであろう。
Claims (20)
- 末梢動脈流動障害検出システムであって、
末梢肢の一部を受けるように輪郭形成された筐体と、
筐体内にあって、検出された動脈流に基づいてパルス波形出力を生成するプレチスモグラフィ信号センサと、
前記センサと動作可能に接続されたホストコンピュータと
を備え、
ホストコンピュータは、時間領域特徴計算値および周波数領域特徴計算値からなるグループから選択された少なくとも2つの計算値を含む、パルス波形に基づく値を計算するように構成され、
nは、ホストコンピュータにより行われる時間領域特徴計算および周波数領域特徴計算の数の総和に等しい値を有し、
ホストコンピュータは、式
P(Dx)は流動障害の確率であり、
v1からvnは、時間領域特徴計算値および周波数領域特徴計算値を含むグループから選択された少なくとも2つの計算値を含み、
係数c0からcnは所定の係数であり、
計算された予知診断を表示するように構成される表示器
を備える、システム。 - センサがフォトダイオードを備える、請求項1に記載のシステム。
- センサが電荷結合素子(CCD)を備える、請求項1に記載のシステム。
- 生理信号が赤外光を含む、請求項1に記載のシステム。
- 係数c0からc6がロジスティック回帰により決定される、請求項1に記載のシステム。
- プレチスモグラフィ信号センサがフォトプレチスモグラフィセンサである、請求項1に記載のシステム。
- プレチスモグラフィ信号センサが容積プレチスモグラフィセンサである、請求項1に記載のシステム。
- 循環障害検出システムであって、
筐体と、
循環を示す生理信号を受けるように配置された筐体内のセンサと、
生理信号を受信するように前記センサと動作可能に接続され、生理信号を分析し、a)時間領域特徴計算値およびb)周波数領域特徴計算値からなるグループから選択された少なくとも2つの計算値を計算するように構成されるホストコンピュータと
を備え、
ホストコンピュータは、循環障害の確率を予測するロジスティック関数において、前記少なくとも2つの計算値を用いて予知診断を計算するようにさらに構成され、
計算された予知診断を表示するように構成される表示器
を備える、システム。 - センサがフォトダイオードを備える、請求項8に記載のシステム。
- センサが電荷結合素子(CCD)を備える、請求項8に記載のシステム。
- センサが、圧力変換器で構成される圧力カフを備える、請求項8に記載のシステム。
- センサが歪みゲージを備える、請求項8に記載のシステム。
- 生理信号が赤外光を含む、請求項8に記載のシステム。
- 生理信号が器官容積を含む、請求項8に記載のシステム。
- 生理信号が皮膚温度を含む、請求項8に記載のシステム。
- 生理信号が皮膚張力を含む、請求項8に記載のシステム。
- 生理信号が血流速度を含む、請求項8に記載のシステム。
- 末梢動脈流動障害検出システムであって、
末梢肢の一部を受けるように輪郭形成された筐体であって、末梢肢の一部からフォトプレチスモグラフィ信号を検出可能なセンサをさらに備え、センサがパルス波形を生成する、筐体と、
前記センサと動作可能に接続されたホストコンピュータであって、循環指数、高調波の傾き、高調波切片、スペクトル信号、および収縮期上昇を含む、パルス波形に基づく値を得るように構成される、ホストコンピュータと
を備え、
ホストコンピュータは、式
P(Dx)は流動障害の確率であり、
v1は循環指数(足/手最大値)であり、
v2は高調波の傾き(足/手最大値)であり、
v3は高調波切片であり、
v4は高調波切片(足−手最大値)であり、
v5はスペクトル信号(足/手最大値)であり、
v6は収縮期上昇であり、
係数c0からc6は所定の係数であり、
計算された予知診断を表示するように構成される表示器
を備える、システム。 - c0が15.99から20.11までの範囲の値を有し、
c1が−33.94から−38.76までの範囲の値を有し、
c2が3.80から4.38までの範囲の値を有し、
c3が−7.16から−8.46までの範囲の値を有し、
c4が5.12から6.52までの範囲の値を有し、
c5が−2.14から3.28までの範囲の値を有し、
c6が41.81から45.61までの範囲の値を有する、
請求項18に記載のシステム。 - c0が18.05の値を有し、
c1が−36.35の値を有し、
c2が4.09の値を有し、
c3が−7.81の値を有し、
c4が5.82の値を有し、
c5が2.71の値を有し、
c6が43.71の値を有する、
請求項19に記載のシステム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/866,731 US20140316292A1 (en) | 2013-04-19 | 2013-04-19 | Circulation Monitoring System |
US13/866,731 | 2013-04-19 | ||
PCT/US2014/033615 WO2014172177A1 (en) | 2013-04-19 | 2014-04-10 | Circulation monitoring system |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2016521162A true JP2016521162A (ja) | 2016-07-21 |
Family
ID=51729544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016508967A Pending JP2016521162A (ja) | 2013-04-19 | 2014-04-10 | 循環監視システム |
Country Status (10)
Country | Link |
---|---|
US (1) | US20140316292A1 (ja) |
JP (1) | JP2016521162A (ja) |
CN (1) | CN105358058A (ja) |
AU (1) | AU2014254248A1 (ja) |
CA (1) | CA2908656A1 (ja) |
DE (1) | DE112014002005T5 (ja) |
GB (1) | GB2526503A (ja) |
MX (1) | MX2015014609A (ja) |
TW (1) | TW201500032A (ja) |
WO (1) | WO2014172177A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20220049060A (ko) * | 2020-10-13 | 2022-04-21 | 재단법인 아산사회복지재단 | 피부 장력을 이용한 림프부종 진단 장치 및 이를 이용한 림프부종 진단 방법 |
KR20220165509A (ko) * | 2021-06-08 | 2022-12-15 | (주)엘탑 | 생체신호 측정 장치 및 방법 |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105796051B (zh) * | 2014-12-31 | 2019-11-19 | 原相科技股份有限公司 | 三维生理检测系统及其操作方法 |
US9986923B2 (en) * | 2015-01-09 | 2018-06-05 | Xerox Corporation | Selecting a region of interest for extracting physiological parameters from a video of a subject |
EP3282932A4 (en) * | 2015-04-13 | 2018-12-05 | Harbor Vascular, Inc. | Methods and system for assessment of peripheral perfusion |
WO2017048858A1 (en) * | 2015-09-14 | 2017-03-23 | Whoop, Inc. | Physiological signal monitoring |
US10159444B2 (en) | 2015-12-21 | 2018-12-25 | Industrial Technology Research Institute | Method and system for anaerobic threshold heart rate detection |
TWI555508B (zh) * | 2015-12-21 | 2016-11-01 | 財團法人工業技術研究院 | 無氧閾心率偵測方法及系統 |
DE102016115496A1 (de) * | 2016-08-22 | 2018-02-22 | B. Braun Avitum Ag | Verfahren und System zur post-dialytischen Bestimmung des Trockengewichts |
US10743821B2 (en) | 2016-10-21 | 2020-08-18 | Tata Consultancy Services Limited | Anomaly detection by self-learning of sensor signals |
JP2018068428A (ja) * | 2016-10-25 | 2018-05-10 | ソニー株式会社 | 情報処理装置、情報処理方法及びプログラム |
EP3363351B1 (en) * | 2017-02-16 | 2023-08-16 | Tata Consultancy Services Limited | System for detection of coronary artery disease in a person using a fusion approach |
EP3387990B1 (en) * | 2017-03-28 | 2020-04-29 | IMEC vzw | System and method for heart rate detection with motion artifact reduction |
JP6965066B2 (ja) * | 2017-09-12 | 2021-11-10 | オムロン株式会社 | 脈波測定装置、血圧測定装置、機器、脈波測定方法、および血圧測定方法 |
KR102510174B1 (ko) | 2017-09-13 | 2023-03-14 | 삼성전자주식회사 | 질병 예측 모델 생성 장치, 질병 예측 장치 및 방법 |
JP7195323B2 (ja) * | 2017-12-20 | 2022-12-23 | エドワーズ ライフサイエンシーズ コーポレイション | 組織オキシメトリと血圧を用いた自己調節能システムと方法 |
RU2688386C1 (ru) * | 2017-12-25 | 2019-05-22 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Юго-Западный государственный университет" (ЮЗГУ) | Способ мониторинга степени тяжести ишемических поражений нижних конечностей и устройство для его осуществления |
KR20200142265A (ko) | 2019-06-12 | 2020-12-22 | 삼성전자주식회사 | 생체정보 추정 장치 및 방법 |
US11723599B2 (en) * | 2019-08-13 | 2023-08-15 | Tawaun Bell | Lower extremity diagnostic device |
CN110801215B (zh) * | 2019-11-22 | 2022-08-09 | 常州市第二人民医院 | 一种智能检测血量仪 |
CN113100735B (zh) * | 2021-03-30 | 2022-08-05 | 桂林电子科技大学 | 一种基于静脉血流图的血栓检测装置及方法 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1121798A (zh) * | 1994-08-16 | 1996-05-08 | 北京工业大学 | 心血管功能动态参数检测分析方法及其装置 |
IL120881A (en) * | 1996-07-30 | 2002-09-12 | It M R Medic L Cm 1997 Ltd | Method and device for continuous and non-invasive monitoring of peripheral arterial tone |
US6592518B2 (en) * | 2001-04-05 | 2003-07-15 | Kenergy, Inc. | Cardiac monitoring system and method with multiple implanted transponders |
US20030032871A1 (en) * | 2001-07-18 | 2003-02-13 | New England Medical Center Hospitals, Inc. | Adjustable coefficients to customize predictive instruments |
US7628760B2 (en) * | 2007-02-28 | 2009-12-08 | Semler Scientific, Inc. | Circulation monitoring system and method |
WO2004105596A1 (en) * | 2003-06-03 | 2004-12-09 | Orsense Ltd. | Method and system for use in non-invasive optical measurements of blood parameters |
GB2413078C (en) * | 2004-01-08 | 2012-08-15 | Dialog Devices Ltd | A system or method for assessing a subject's pedalblood circulation. |
AU2005286769A1 (en) * | 2004-09-21 | 2006-03-30 | Adidas Ag | Inductive plethysmographic sensors, monitors, and apparel |
CA2675826A1 (en) * | 2006-01-30 | 2007-08-02 | Cardiosense Ltd. | Apparatus, system and method for determining cardio-respiratory state |
GB0603006D0 (en) * | 2006-02-15 | 2006-03-29 | Dialog Devices Ltd | Assessing blood supply to a peripheral portion of an animal |
US20090281413A1 (en) * | 2007-12-18 | 2009-11-12 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Systems, devices, and methods for detecting occlusions in a biological subject |
US20090281412A1 (en) * | 2007-12-18 | 2009-11-12 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | System, devices, and methods for detecting occlusions in a biological subject |
US9173593B2 (en) * | 2010-04-19 | 2015-11-03 | Sotera Wireless, Inc. | Body-worn monitor for measuring respiratory rate |
JP5955341B2 (ja) * | 2011-01-27 | 2016-07-20 | ザ ボード オブ トラスティーズ オブ ザ レランド スタンフォード ジュニア ユニバーシティー | 循環系を観察するためのシステム及び方法 |
-
2013
- 2013-04-19 US US13/866,731 patent/US20140316292A1/en not_active Abandoned
-
2014
- 2014-04-10 AU AU2014254248A patent/AU2014254248A1/en not_active Abandoned
- 2014-04-10 CA CA2908656A patent/CA2908656A1/en not_active Abandoned
- 2014-04-10 MX MX2015014609A patent/MX2015014609A/es unknown
- 2014-04-10 CN CN201480022161.9A patent/CN105358058A/zh active Pending
- 2014-04-10 DE DE112014002005.2T patent/DE112014002005T5/de not_active Withdrawn
- 2014-04-10 JP JP2016508967A patent/JP2016521162A/ja active Pending
- 2014-04-10 WO PCT/US2014/033615 patent/WO2014172177A1/en active Application Filing
- 2014-04-10 GB GB1517436.0A patent/GB2526503A/en not_active Withdrawn
- 2014-04-18 TW TW103114274A patent/TW201500032A/zh unknown
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20220049060A (ko) * | 2020-10-13 | 2022-04-21 | 재단법인 아산사회복지재단 | 피부 장력을 이용한 림프부종 진단 장치 및 이를 이용한 림프부종 진단 방법 |
KR102479289B1 (ko) | 2020-10-13 | 2022-12-21 | 재단법인 아산사회복지재단 | 피부 장력을 이용한 림프부종 진단 장치 및 이를 이용한 림프부종 진단 방법 |
KR20230018353A (ko) * | 2020-10-13 | 2023-02-07 | 재단법인 아산사회복지재단 | 피부 장력을 이용한 림프부종 진행 상태 예측 장치 및 방법 |
KR102507388B1 (ko) | 2020-10-13 | 2023-03-08 | 재단법인 아산사회복지재단 | 피부 장력을 이용한 림프부종 진행 상태 예측 장치 및 방법 |
KR20220165509A (ko) * | 2021-06-08 | 2022-12-15 | (주)엘탑 | 생체신호 측정 장치 및 방법 |
KR102544049B1 (ko) | 2021-06-08 | 2023-06-16 | (주)엘탑 | 생체신호 측정 장치 및 방법 |
Also Published As
Publication number | Publication date |
---|---|
WO2014172177A1 (en) | 2014-10-23 |
CA2908656A1 (en) | 2014-10-23 |
AU2014254248A1 (en) | 2015-10-29 |
DE112014002005T5 (de) | 2016-01-21 |
MX2015014609A (es) | 2016-06-02 |
US20140316292A1 (en) | 2014-10-23 |
TW201500032A (zh) | 2015-01-01 |
GB2526503A (en) | 2015-11-25 |
CN105358058A (zh) | 2016-02-24 |
GB201517436D0 (en) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2016521162A (ja) | 循環監視システム | |
El-Hajj et al. | A review of machine learning techniques in photoplethysmography for the non-invasive cuff-less measurement of blood pressure | |
Charlton et al. | Breathing rate estimation from the electrocardiogram and photoplethysmogram: A review | |
US20210030372A1 (en) | Methods to estimate the blood pressure and the arterial stiffness based on photoplethysmographic (ppg) signals | |
JP6130474B2 (ja) | 体重計装置及び脈波速度取得方法 | |
Sugita et al. | Contactless technique for measuring blood-pressure variability from one region in video plethysmography | |
CN110944575B (zh) | 通过使用相干光源进行生物组织询问的非侵入式血液动力学评估 | |
Natarajan et al. | Photoplethysmography fast upstroke time intervals can be useful features for cuff-less measurement of blood pressure changes in humans | |
JP6407757B2 (ja) | 指標算出装置、指標算出システム、指標算出方法、および制御プログラム | |
Al–Fahoum et al. | A multiple signal classification approach for photoplethysmography signals in healthy and athletic subjects | |
EP2265170A1 (en) | Circulation monitoring system and method | |
Tanaka et al. | Accuracy assessment of a noninvasive device for monitoring beat-by-beat blood pressure in the radial artery using the volume-compensation method | |
Alqudah et al. | Multiple time and spectral analysis techniques for comparing the PhotoPlethysmography to PiezoelectricPlethysmography with electrocardiography | |
Suganthi et al. | Morphological analysis of peripheral arterial signals in Takayasu’s arteritis | |
WO2021249850A1 (en) | Wearable device | |
Johnson et al. | A Review of Photoplethysmography-based Physiological Measurement and Estimation, Part 1: Single Input Methods | |
More et al. | Development of non-invasive diagnostic tool for diseases using Photo Plethysmography | |
Li et al. | The establishment of a non-invasive continuous blood pressure measure system based on pulse transit time | |
Jegan et al. | Methodological role of mathematics to estimate human blood pressure through biosensors | |
Chaurasia et al. | Development of a Low Cost Heart Rate Monitoring and Transmission System using PPG Signal Processing for Wearable Devices | |
Ramsrinivas et al. | Virtual Instrumentation based Non-Invasive Assessment of Arterial Stiffness Using Finger Photoplethysmographic Signal | |
He | Signal Enhancement Applied to Pulse Transit Time Measurement | |
WO2023108177A2 (en) | Pulsewave velocity detection device and the hemodynamic determination of hba1c, arterial age and calcium score | |
El-Hajj et al. | ORCID: 0000-0002-2868-485X (2020). A review of machine learning techniques in photoplethysmography for the non-invasive cuff-less measurement of blood pressure | |
Annapoorani | Ankle Brachial Index (ABI) Based Leg Attack Diagnosis using LabVIEW |