JP6059726B2 - 歪み低減された信号検出 - Google Patents
歪み低減された信号検出 Download PDFInfo
- Publication number
- JP6059726B2 JP6059726B2 JP2014532513A JP2014532513A JP6059726B2 JP 6059726 B2 JP6059726 B2 JP 6059726B2 JP 2014532513 A JP2014532513 A JP 2014532513A JP 2014532513 A JP2014532513 A JP 2014532513A JP 6059726 B2 JP6059726 B2 JP 6059726B2
- Authority
- JP
- Japan
- Prior art keywords
- reflection
- signal
- index element
- basic
- disturbing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001514 detection method Methods 0.000 title description 18
- 230000002829 reductive effect Effects 0.000 title description 2
- 238000012545 processing Methods 0.000 claims description 44
- 238000000034 method Methods 0.000 claims description 26
- 238000010606 normalization Methods 0.000 claims description 22
- 230000000737 periodic effect Effects 0.000 claims description 19
- 230000000295 complement effect Effects 0.000 claims description 18
- 230000005670 electromagnetic radiation Effects 0.000 claims description 18
- 230000011514 reflex Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 5
- 230000000996 additive effect Effects 0.000 claims description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 229910052760 oxygen Inorganic materials 0.000 claims description 5
- 239000001301 oxygen Substances 0.000 claims description 5
- 230000036962 time dependent Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000010349 pulsation Effects 0.000 claims description 4
- 230000036387 respiratory rate Effects 0.000 claims description 4
- 230000003139 buffering effect Effects 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims 1
- 239000013598 vector Substances 0.000 description 154
- 239000011159 matrix material Substances 0.000 description 33
- 238000013459 approach Methods 0.000 description 31
- 230000005855 radiation Effects 0.000 description 31
- 238000001914 filtration Methods 0.000 description 22
- 238000005286 illumination Methods 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 11
- 230000033001 locomotion Effects 0.000 description 11
- 230000036555 skin type Effects 0.000 description 10
- 238000012544 monitoring process Methods 0.000 description 9
- 230000001419 dependent effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000013186 photoplethysmography Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 239000003086 colorant Substances 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000002526 effect on cardiovascular system Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000004931 aggregating effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001815 facial effect Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 238000002106 pulse oximetry Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 208000018262 Peripheral vascular disease Diseases 0.000 description 1
- 206010039792 Seborrhoea Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000005352 clarification Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 238000002565 electrocardiography Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000037312 oily skin Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 229920001451 polypropylene glycol Polymers 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000001454 recorded image Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000014624 response to red light Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000035900 sweating Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 230000000007 visual effect Effects 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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
- A61B5/14551—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 for measuring blood gases
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/10—Image acquisition
- G06V10/12—Details of acquisition arrangements; Constructional details thereof
- G06V10/14—Optical characteristics of the device performing the acquisition or on the illumination arrangements
- G06V10/143—Sensing or illuminating at different wavelengths
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/7715—Feature extraction, e.g. by transforming the feature space, e.g. multi-dimensional scaling [MDS]; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30076—Plethysmography
Description
− 画像のシーケンスを得るステップと、
− 画像中に表された対象物の状態を分類するためのデータを得るために画像のシーケンスの少なくとも1つに対して視覚ベースの分析を実行するステップと、
− 画像のシーケンスの少なくとも幾つかにおいて表される生物の生理学的パラメータの少なくとも1つの値を決定するステップであって、生理学的パラメータの少なくとも1つの値は、視覚ベースの分析が実行された少なくとも1つの画像が得られた同じ画像のシーケンスからの画像データの分析によって決定される、ステップと、
− 視覚ベースの分析及び生理学的パラメータの少なくとも1つの値によって得られたデータを用いて対象物の状態を分類するステップと
を含む、画像分析のための方法および装置を開示している。
− オブジェクトにより反射された電磁放射から導出可能なデータストリームを受信するためのインタフェースであって、前記データストリームは、連続又は別個の時系列特性信号を有し、前記特性信号は、信号空間に関連付けられた少なくとも2つの主成分を有し、前記信号空間は、前記特性信号を表すための相補的なチャネルを有し、前記特性信号の主成分は、前記信号空間のそれぞれの相補的なチャネルに関連する、インタフェースと、
− 線形代数方程式を特定するように実質的な線形代数信号組成モデルの考察の下で、規定された成分表現に前記特性信号をマッピングするためのコンバータ手段であって、前記特性信号は、生理学的反射成分及び妨害反射成分を少なくとも実質的に含み、前記生理学的反射成分は、観察されるべきオブジェクトの基本的な生理学的特性に対して起因性のある基本的反射部分と少なくとも1つの少なくとも部分的に周期的なバイタル信号に対して起因性のある指示的反射部分とを含む、コンバータ手段とを有し、
− 前記基本的反射部分は、基本的反射インデックス要素及び基本反射係数を有し、
− 前記指示的反射部分は、指示的反射インデックス要素及び時間依存性指示的反射係数を有し、前記指示的反射係数は、関心のありそうな前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を高度に示し、
− 前記妨害反射成分は、妨害反射インデックス要素及び妨害反射係数を有し、
当該システムは、
− 前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を高度に示す表現を導出するために、前記基本的反射インデックス要素、前記指示的反射インデックス要素及び前記妨害反射インデックス要素の向きの少なくとも近似の推定の考察の下で前記線形代数方程式を少なくとも部分的に解くための処理手段であって、前記の高度に示す表現は、前記指示的反射係数を少なくとも有する、処理手段を有する。
− オブジェクトにより反射された電磁放射から導出可能なデータストリームを受信するステップであって、前記データストリームは、連続又は別個の時系列特性信号を有し、前記特性信号は、信号空間に関連付けられた少なくとも2つの主成分を有し、前記信号空間は、前記特性信号を表すための相補的なチャネルを有し、前記特性信号の主成分は、前記信号空間のそれぞれの相補的なチャネルに関連する、ステップと、
− 線形代数方程式を特定するために実質的な線形代数的信号組成モデルの考察の下で規定された成分表現に対して前記特性信号をマッピングするステップであって、前記特性信号は、生理学的反射成分及び妨害反射成分を少なくとも実質的に含み、前記生理学的反射成分は、観察されるべきオブジェクトの基本的生理学的特性に対して起因性のある基本的反射部分と、少なくとも1つの少なくとも部分的に周期的なバイタル信号に対して起因性のある指示的反射部分とを有する、ステップとを有し、
− 前記基本的反射部分は、基本的反射インデックス要素及び基本的反射係数を有し、
− 前記指示的反射部分は、指示的反射インデックス要素及び時間依存性指示的反射係数を有し、前記指示的反射係数は、関心のありそうな少なくとも1つの少なくとも部分的に周期的なバイタル信号を高度に示し、
− 前記妨害反射成分は、妨害反射インデックス要素及び妨害反射係数を有し、
当該方法は、
− 前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を高度に示す表現を導出するために、基本的反射インデックス要素、指示的反射インデックス要素及び妨害反射インデックス要素の向きの少なくとも近似の推定の考察の下で線形代数方程式を少なくとも部分的に解くステップであって、前記の高度に示す表現は、指示的反射係数を少なくとも有する、ステップを有する。
及び
の割合
が解析され得る。
ここで、
及び
は、それぞれ、信号
及び
の平均値に対応する。割合は、LEDの強度のための、及び、皮膚色のための信号を正規化するように見なされる。基本的に、信号
は、心血管パルスを示し、心拍数は、基本的に、生ずる信号
の周波数分析を介して決定され得る。
の平均
が基本的に1に等しいと推定するときに、式(1)が書き直され得る。信号
の(例えば心拍数に関連した)交互の挙動が特に興味あるものであり得るので、平均
(即ち、1)は、適宜減算され得る。
この文書において、赤外線に対する赤色の割合は、所望の信号の抽出に対する例示的なアプローチであることが理解されるべきである。例えば、緑色に対する赤色の割合が利用され得る。規定されたスペクトル部分に関連する信号成分までの関連する他の割合は、同様に見なされ得る。
の表現を単純化するために利用され得る。これは、割合
及び
が両方とも1に近いことに従う仮定を考慮し、対数関数的表現、例えば、
の近似値のための対数演算子のテイラー列を更に利用して、
をもたらす。対数演算子の利点は、基本的に、(例えば、光源(例えばLED)の強度、及び、関心のありそうなオブジェクトが示す皮膚の色に向けられる)輝度正規化は、対数関数的演算子が
をもたらすので、回避され得るという事実において理解され得る。条件
は、関心のありそうな周波数範囲の外側の歪みを低減するような信号のバンドパスフィルタリングを考慮したときに、除去され得る。
ここで、iは、入射角であり、eは、放射率の角度であり、gは、図3に示された位相角である。図3において、参照符号70は、照明(入射光)の方向に対応する。参照符号74,68は、説明の便宜上図3において一体化されたそれぞれの面法線を示す。参照符号76は、関心のありそうなオブジェクトの界面で反射される鏡面反射成分を表す。参照符号84は、表示方向を示す。更に、完全性のために、図3は、鏡面反射76の方向とセンサ手段24の位置及び/又は向きを示す表示方向84との間の角度を示す角度sを示す。
以下において、角度sは、式(5)のスペクトルパワー分布sと混同されるべきでない。
に移され得る。例えば、赤色、緑色及び青色フィルタのようなセンサ手段のフィルタ特性が考慮され得る。全体で、式(5)は、反射された光の全体強度をもたらすように適用される照明のスペクトルに渡って一体化され得る。式(6)におけるこの結果は、関心のありそうな表面上の特定のポイントのための、信号空間(例えばRGB色空間)に関する二色性反射モデルを与える。大きさms及びmbは、スカラーとみなされ得る。ベクトル
は、測定された色、または、特性信号を表す。係数ms及びmbは、問題となっているポイントの反射の大きさと見なされ得る。ベクトル
及び
は、関心のありそうなオブジェクトの(鏡面反射とも呼ばれる)界面反射及び(生理学的反射とも呼ばれる)ボディ反射の色を表す。
(特性信号とも呼ばれる)は、オブジェクト表面上の一組の(ピクセル)ポイントを表し、これは、多少一様に着色され得る(皮膚色)と仮定され得る。基本的に、幾何学に関連した因子は、次から次へと異なる。それ故、スケールファクタms及びmbは、次から次へと変化し得る。しかしながら、界面及びボディ反射の色
及び
は、幾何学によって変化しないs(λ)及びb(λ)のスペクトル投射とみなされ得るので、実質的に関心のありそうな表面上のあらゆる点において不変のままである。換言すれば、ベクトル表現を適用したときに、各ピクセル信号
は、双方の(二色性プレーンとも呼ばれる)色ベクトル
及び
に渡る面にある。従って、色
は、(同じスペクトルであるが異なる強度を有する)照明状態から独立しているものとみなされ得る。
又は生理学的反射成分に関する他の仮定により拡張され得る。上で述べられるように、オブジェクトの皮膚組織の吸収度の僅かな変化は、関心のありそうなバイタル信号を抽出するために検出され得る。換言すれば、(実際の)ボディ反射ベクトル
は、平均ボディ反射ベクトル
の周りで時間的に(僅かに)振動する。便宜のために、これらの振動は、ベクトルにより述べられ得るか、又は、ベクトルにより少なくとも線形的に近付けられ得ると仮定され得る。従って、ベクトルは、心拍ベクトル
として分類され得る。ベクトル
を二色性反射モデルに取り込むことは、信号検出及び処理を促進するために二色性反射モデルを拡張するための適切な手段と見なされ得る。従って、関心のありそうなピクセル(又は、ピクセルパターン)の特性信号
(例えば、RGB値)は、(平均)ボディ反射ベクトル
、界面ベクトル
及び心拍ベクトル
の実質的な線形的組み合わせにより述べられ得る。
便宜のために、平均ボディ反射ベクトル
のハットは、以下のように落とされ得る。
(例えば、心拍ベクトル)は、一定のベクトルではなく、照明の色、カメラの応答性、赤色、緑色及び青色チャネルのフィルタ応答及び対象物の皮膚色に依存することが理解されるべきである。基本的に、同じことがベクトル
及び
に当てはまる。従って、ベクトル
は、センサ手段(カメラ)ごとに、及び、放射線供給源(照明)のタイプごとに決定されなければならない。しかしながら、以下の解明に基づいて、皮膚タイプ(即ち、皮膚色)は役割を果たさないと結論付けられ得る。それ故、一般に、これは、或る皮膚タイプを有する関心のありそうな例示的なオブジェクトのための安定した照明状態の下でセンサ手段のための心拍ベクトル
を推定するのに十分であると見なされ得る。その後、(実際の)心拍ベクトル
は、実際の皮膚タイプのためのボディ反射ベクトル
の推定を用いることにより個別に推定され得る。潜在的に、この洞察は、以下のように利用され得る。関心のありそうな対象物が、時間の短い期間の間(例えば、数秒)、依然として残る場合には、心拍ベクトル
は、時間のこの短い期間において関心のありそうなオブジェクトの皮膚ピクセルの平均の赤色、平均の緑色及び平均の青色値に適用される本質成分分析を介した中間手段において容易に決定され得る。基本的に、それぞれの平均の赤色、平均の緑色及び平均の青色値は、実際のボディ反射ベクトル
を示すものとみなされ得る。
は、それぞれの赤色、緑色及び青色信号の指示的周波数ピークの考察を介して導出され得る。この目的を達成するために、例えば、フーリエ変換は、それぞれの大きさスペクトラムを得るために信号に適用され得る。同じ周波数での赤色、緑色及び青色信号の各々に少なくとも実質的にあるかなりのピークは、心拍ベクトル
を示し得る。従って、それぞれの赤色、緑色及び青色ピークの大きさは、ベクトル
を"構成する"ために利用され得る。それ故、心拍ベクトル
の向きが決定され得る。言うまでもなく、心拍ベクトル
は、ベクトル
が心拍に関連した信号部分
の(想定された)向きを単に示すだけなので、単位ベクトル又は類似のものに更に変換され得る。
は、強度値(例えば、pR,pG及びpB)を示す主成分92a,92b,92cにより表され、特性信号98(又は
)は、それぞれの相補的なチャネル90a,90b,90cに関連するもので構成される。参照符号94は、輝度インデックス要素を示す。輝度インデックス要素94は、入射放射線(例えば光源)の供給源を示す。輝度インデックス要素94は、信号空間88を横断する斜めのベクトルと見なされ得る。これは、とりわけ、放射線源16が基本的に単調な白色光を放射するときに適用する。好ましくは、放射線源16の"色"は、信号空間88の白色ポイント(例えば1,1,1)に等しい。
)は、生理学的反射成分100及び妨害反射部分102(又は
)の加法的な組み合わせと見なされ得る。生理学的な反射成分100は、基本的な反射部分104(又は、
及び、
のそれぞれ)及び指示的反射部分106(又は
)で構成されると更にみなされ得る。換言すれば、生理学的反射成分100は、基本的反射部分104及び指示的反射部分106の線形的組み合わせと見なされ得る。基本的反射部分104は、平均ボディ反射ベクトル
(又は
)及びそれぞれの係数c0を表す。上で概説されるように、(実際の)ボディ反射ベクトル
は、平均ボディ反射ベクトル
(又は
)周辺で時間的に振動する。それ故、完全性のために、図4a及び図4bは、生理学的反射成分100とも呼ばれる、それぞれの実際のボディ反射条件(
)を示す。条件
、又は、妨害反射部分102は、鏡面反射が基本的に入射光を映すので、輝度インデックス要素94に平行と見なされる。それ故、それぞれの信号成分102は、例えば輝度インデックス要素94により示される、光源特性を少なくとも部分的に示す。
及び
に及ぶ。以下において、
,
及び
、は、向き及び方向を示す単位ベクトルと見なされ得る一方で、c0,c1及びc2、は、長さを示すそれぞれの係数と見なされ得る。指示的反射部分106(又は
)を無視するときには、残りの信号成分102,104は、二色性プレーン96にあると見なされ得る。
は、基本的にゼロに等しいと推定され得る。従って、センサ手段を介して得られる特性信号
は、ボディ反射ベクトル
及び心拍ベクトル
(c2=0による式(7)参照)の線形的組み合わせにより述べられ得る。
ここで、ベクトル
は、(以下において主成分とも呼ばれる)それぞれの赤色及び赤外線信号を含み得る。nは、時間に対応する。又は、より明示的には、
である。この文書において、RGB信号もベクトル
の主成分として利用され得ることが理解されるべきである。例えば、赤色及び緑色信号は、関心のありそうな所望のバイタル信号を検出するように解析され得る。
(基本的に心拍信号はゼロ平均を示す)の考察の下で条件を更に正規化したときには、
平均ベクトル
は、
に等しい。
及びボディ反射ベクトル
による除算を有する。ボディ反射ベクトル
は、オブジェクトの皮膚により赤色光及び赤外線の(平均の)反射の結果と見なされ得る。この除算は、正規化されたベクトル
をもたらす。
ここで、
は、常の
標準を示し、対角行列Bは、対角線上にボディ反射ベクトル
を有し、ベクトル
は、値1を含むベクトルであり、ベクトル
は、皮膚タイプから独立している正規化された心拍ベクトルに対応する。
に等しい。代替手段において、1つは、ボディ反射ベクトル
及び正規化された心拍ベクトル
の推定から始まり得る。ここから、心拍ベクトル
は以下のように導出され得る。例えば、RGB信号空間において、
である。この表現は、単位長に対する正規化と組み合わせにおける式(9)に基づく。
は、ボディ反射ベクトル
、心拍ベクトル
及び界面反射ベクトル
の線形的組み合わせとして述べられ得る(式(7)も参照)。
により除算(又はこれらに対して正規化)され得る。しかしながら、これは、オプショナルなステップと見なされ得る。とりわけ、小さな強度変化だけが予期されるときには、強度正規化は省略され得る。他の例示的な正規化ステップは、生ずる信号のボディ反射ベクトル
による除算を有し得る。
及び
において書き直され得る。ここで、c0は、基本的反射係数とも呼ばれ得る。c1は、指示的反射係数とも呼ばれ得る。c2は、妨害反射係数とも呼ばれ得る。
式(12)は、規定された成分表現の考察の下で構築される線形代数方程式と見なされ得る。従って、関心のありそうなバイタル信号u(n)、例えば、心拍信号は、
又は単純に
になる。更に、信号uは、関心のありそうな周波数範囲の外側の周波数コンテンツをフィルタリングするためにバンドパスフィルタによりフィルタリングされ得る。
及び鏡面反射ベクトル
が推定されなければならない。更に、上で示すように、心拍ベクトル
は、式(10)を参照して、ボディ反射ベクトル
から導出されてもよい。鏡面反射条件に関して、皮膚のような高い油分又は水分を有する材料に関して、材料の表面から反射される光は、基本的に、発光物と同じ色を呈するように見えることが考慮され得る。第1の近似値として、照明が白色(例えば、RGB=1,1,1)であることが明確に想定され得る。即ち、鏡面反射ベクトル
は、1のベクトル
に等しいと見なされ得る。更に、ボディ反射ベクトル
は、フレームの獲得されたシーケンスの最初のフレームにおいて関心のありそうなオブジェクトの皮膚ピクセルの平均の赤色、平均の緑色及び平均の青色と同じであると見なされ得る。従って、単位長に対する正規化が適用され得る。
それ故、関心のありそうな信号uは、式(12)及び(13)を適用することにより得られ得る。
及び界面ベクトル
に直交するベクトルを指示的(心拍)ベクトル
として用いることが想定されてもよい。このアプローチは、有利な想定によれば、指示的(心拍)ベクトル
は、
及び
並びにそれ故、基本的に二色性プレーンの外のポイントから独立している(図4a中の参照符号96を参照)という洞察に基づいている。従って、変化する照明状態は、指示的(心拍)ベクトル
に大幅に影響を及ぼさない。
の(主成分とも呼ばれる)平均の赤色、緑色及び青色値は、鏡面反射により依然として影響を受ける。それ故、平均の鏡面反射値は、依然として、平均の赤色、緑色及び青色値にある。
。
以下の条件が適宜生じる。
ここで、対角行列Bの要素は、ボディ反射ベクトル
の要素に等しいものと見なされる。この表現から、鏡面反射を伴わない理想的な事例に関して、ボディ反射ベクトル
及び心拍ベクトル
に及ぶRGB空間は、強度ベクトル
及び正規化された心拍ベクトル
に及ぶlog RGB空間に変換され得ることが結論付けられ得る。
。そして以下の条件が導出され得る。
この表現から、ボディ反射ベクトル
、心拍ベクトル
及び界面反射ベクトル
に及ぶRGB空間は、正規化されたボディ反射ベクトル
、正規化された心拍ベクトル
及び変更された界面反射ベクトル
に及ぶlog RGB空間に変換され得ることが結論付けられ得る。即ち、心拍信号は、先に述べたものと類似の手段で抽出され得る。
。従って、以下の条件が導出され得る。
対角行列Sの要素は、界面反射ベクトル
の要素に等しい。この表現は、ボディ反射ベクトル
、心拍ベクトル
及び界面反射ベクトル
に及ぶRGB空間は、界面反射ベクトル
、変更された心拍ベクトル
及び変更されたボディ反射ベクトル
に及ぶlog RGB空間に変換され得ることに従った結論を可能にする。結果として、既知の心拍ベクトル
を有する代わりに、発光物の色に依存する心拍ベクトルが得られ得る。更に、心拍信号の振幅は、ボディ反射の量
の代わりに、鏡面反射の量
に依存する。従って、この場合において、鏡面反射の量が比較的大きいので、log RGB空間は、心拍信号を抽出するのにあまり適切ではないと見なされ得る。
。結果として、以下の条件が導出され得る。
ここで、
であり、対角行列Gの要素は、ベクトル
の要素に等しい。この表現に基づいて、ボディ反射ベクトル
、心拍ベクトル
及び界面反射ベクトル
に及ぶRGB空間が、ベクトル
及び変更された心拍ベクトル
に及ぶlog RGB空間に変換され得ると結論付けられ得る。
に直交するプレーン上へRGB空間を最初に投射することにより鏡面反射を低減することが考えられ得る。この場合、鏡面反射を適切に低減するように、光源色を決定することが有利であるだろう。しかしながら、光源色が知られている場合、例えば、基本的に白色の発光物を含む最も一般的な事例に対して、このアプローチは適用可能である。
を処理するために設けられ、上で導入された、及び、以下で示される例により利用され得る幾つかのステップを有する。言うまでもなく、ステップの幾つかは、オプションであると見なされ、これは、破線により示される。第1のステップ110において、特性入力信号
がフィルタリングされる。フィルタリングは、バンドパスフィルタリング、ローパスフィルタリング及びハイパスフィルタリングを有し得る。他のフィルタ特性も同様に適用され得る。(入力)フィルタリングは、非指示的周波数部分を減衰させること、及び、特性信号
の潜在的な指示的周波数部分を強化することに向けられ得る。
ベクトル−行列表現において、この式は、条件
により表され得る。
又は、少なくとも
を得るように、ステップ114において少なくとも部分的に解かれ得る。
又は平均係数
により、これを除算することにより、強度に対して正規化され得る。従って、関心のありそうな所望の信号を高度に示す信号u(n)が得られ得る。幾つかのアプリケーションに関して、係数
(又は
)は、その不所望な変換を回避するように、オプショナルなステップ116の前にフィルタリングステップを経るべきでないことに留意されるべきである。おそらく、これらの変換は、ステップ116において計算されるべき割合に更に悪影響を与え得る。
における平均の赤色、平均の緑色及び平均の青色値は、行列乗算によってベクトル
に変換される。そして、オプショナルな手段として、心拍信号c1は強度に対して補正され、心拍信号をもたらす。この心拍信号を更に強化することに関して、心拍信号uは、バンドパスフィルタによりフィルタリングされ得る。典型的には、関心のありそうな周波数範囲、例えば、約35〜210BPM(beats per minute)の範囲を含む通過帯域を示すことが利用され得る。基本的に、以下の例は、ベクトル空間ベースの観点に依存する。
が二色性プレーンに直交すると見なされる内積
は、所望の(心拍)信号を高度に示すと見なされ得る。
は、例えばフレームのシーケンスの第1のフレームにおける皮膚ピクセルの平均の赤色、平均の緑色及び平均の青色値を計算することにより、一回推定され得る。更に、発光物の色が白色(即ち、界面反射ベクトル
)であることが明確に想定され得る。従って、心拍ベクトル
は、式(10)を参照して、ボディ反射ベクトルから直接追跡し得る。変換行列Tは、
に等しく、従って、式(12)及び(13)は、心拍信号uを計算するために用いられ得る。加えて、この心拍信号は、関心のありそうな周波数範囲の外の周波数コンテンツをフィルタリングするために、バンドパスフィルタによりフィルタリングされ得る。
は、二色性のプレーンに直交するベクトル、即ち
に等しいと見なされ得る。このアプローチは、発光物スペクトルが知られていない場合において好ましく、従って、正規化された心拍ベクトル
は決定され得ない。それ故、3つのベクトルが与えられると、変換行列Tは、
と構成され、前と同じように、式(12)及び(13)が心拍信号uを計算するために利用され得る。
は、これを(対角行列Bにより表される)ボディ反射ベクトル
で除算することにより正規化され得る。
好ましくは、正規化は即座に実行される。正規化は、心拍ベクトル
から独立しており、基本的に一定であり、それ故に心拍信号を含まないボディ反射ベクトル
に依存する。
ベクトルにより表され、心拍ベクトルは、正規化された心拍ベクトル
になり、界面(鏡面)反射ベクトルは、正規化された界面反射ベクトル
になる。従って、これらの正規化されたベクトルは、行列Tを構成するために用いられ、
前と同じように、式(12)及び(13)が心拍信号uを計算するために考慮され得る。この例は、主に含まれたボディ反射を示す行列Bがあるので、基本的に前の例に等しいと見なされ得ることに留意されるべきである。
が規定され得る。従って、対角線として
を有する対応する対角行列
が確立され得る。それ故、界面反射の近似値
が導出され得る。
ここで、
は、近付けられたボディ反射ベクトルと見なされ得る。従って、心拍ベクトル
は、
に等しく、行列Tは、
のように構成され、前と同じように、式(12)及び(13)が心拍信号uを計算するために用いられ得る。代わりとして、
は、ボディ反射ベクトル
及び界面反射ベクトル
に直交するベクトルに等しい。
照明の色が時間的に(ゆっくり)変化するものと予期される場合において、ボディ反射ベクトルを連続的に推定及びフィルタリングすることが考慮され得る。フィルタリングは、心拍信号を強化するために、図5中のフィルタリングステップ110を参照して、ローパスフィルタを利用して実行され得る。従って、行列T及びT−1は、連続的に変更及び適合され得る。
に適用される。上記のように、ゼロ値又はより小さい数の対数をとることは回避されるべきである。基本的に、行列Tは、例えば、前の例のうちの1つに従って構成され得る。従って、逆行列T−1は、ベクトル
に向けた信号に適用され得る。
この例によれば、この得られたベクトル
の第3の列は、基本的に鏡面反射に対応するので、破棄され得る。換言すれば、T−1の最初の2列に対応する列を有する2×3行列
が規定され、この行列は、信号pに適用され得る。
を有する強度信号を表す一方で、信号c1は、ゼロに近い平均を有する心拍信号を表す。それ故、この信号
に対数演算子を適用することを可能とするために、基本変化は、条件
に行列Uを適用することにより実行され得る。この変換の目的は、心拍信号がゼロより大幅に大きな平均を有することを確認することである。この行列Uは、基本的に双方の信号がおおよそ
の平均を有する手法において選択され得る。もちろん、例示的な行列Uから離れた他の行列が想定され得る。
を利用することが有利であるだろう。
は、これらの平均の値を含めるために取り込まれ得る(Nは10〜15秒に対応する)。その後、逆行列T−1は、ベクトル
を得るためにこのベクトルに適用され得る。
における信号の各列は、関心のありそうな周波数範囲の外のコンテンツを減衰させるようにフィルタリングされたバンドパスであり得る。バンドパスフィルタリングされた信号c1は、個々の動作があまり激しくないと想定すると、既に、潜在的な指示的心拍信号を含み得る。例えば、これは、信号c1をフーリエ領域に変換して信号対雑音比を検出することにより、2番目に大きなピークに対する最も大きなピークの割合を計算することにより、確認されることができる。この割合が或る閾値を越えた場合、バンドパスフィルタリングされた信号c1は、テンプレートとして用いられるのに十分きれいであると見なされ得る。バンドパスフィルタリングされた信号c0及びc2における心拍振幅a0及びa2は、
及び
により計算され得る。ここで、
は、内積に対応する。基本的には、以下の条件が計算される。
及び
更に、
であるので、a1=1であることに留意されるべきである。従って、この例は、心拍ベクトル
をもたらし、RGB信号空間を参照したときには、心拍ベクトル
は、
になる。
の幾何学的な関連及び信号空間の他の成分を利用して段階的な変換により補正又は置換され得ることが理解されるべきである。換言すれば、行列演算は、段階的なベクトル演算により容易に置換され得る。
Claims (15)
- 検出された特性信号から情報を抽出するためのデバイスであって、
当該デバイスは、
オブジェクトにより反射された電磁放射から導出可能なデータストリームを受信するためのインタフェースであって、前記データストリームは、連続又は別個の時系列特性信号を有し、前記特性信号は、信号空間に関連付けられた少なくとも2つの主成分を有し、前記信号空間は、前記特性信号を表すための相補的なチャネルを有し、前記特性信号の主成分は、前記信号空間のそれぞれの相補的なチャネルに関連する、インタフェースと、
線形代数方程式を特定するように実質的な線形代数信号組成モデルの考察の下で、規定された成分表現に前記特性信号をマッピングするためのコンバータ手段であって、前記特性信号は、生理学的反射成分及び妨害反射成分を少なくとも実質的に含み、前記生理学的反射成分は、観察されるべきオブジェクトの基本的な生理学的特性に対して起因性のある基本的反射部分と少なくとも1つの少なくとも部分的に周期的なバイタル信号に対して起因性のある指示的反射部分とを含む、コンバータ手段とを有し、
前記基本的反射部分は、基本的反射インデックス要素及び基本的反射係数を有し、
前記指示的反射部分は、指示的反射インデックス要素及び時間依存性指示的反射係数を有し、前記指示的反射係数は、関心のありそうな前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を示し、
前記妨害反射成分は、妨害反射インデックス要素及び妨害反射係数を有し、
当該デバイスは、
前記基本的反射インデックス要素、前記指示的反射インデックス要素及び前記妨害反射インデックス要素に基づいて前記線形代数方程式を少なくとも部分的に解くことにより、前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を示す計算された表現を導出するための処理手段であって、前記計算された表現は、前記指示的反射係数を少なくとも有する、処理手段を有する、デバイス。 - 前記少なくとも1つの少なくとも部分的に周期的なバイタル信号は、心拍数、心拍、呼吸数、心拍変動性、Traube -Hering-Mayer波及び酸素飽和度からなるグループから選択され、前記信号空間は、加法的又は減法的な色信号空間であり、前記少なくとも2つの主成分は、前記相補的なチャネルにより示されたそれぞれの異なる色成分を表し、前記相補的なチャネルは、規定されたスペクトル部分に関連する、請求項1に記載のデバイス。
- 前記計算された表現から前記バイタル信号を抽出するための抽出手段を更に有し、前記バイタル信号は、前記指示的反射係数、又は、前記指示的反射係数と前記基本的反射係数との割合の考察の下で抽出される、請求項1に記載のデバイス。
- 前記基本的反射部分は、関心のありそうなオブジェクトが示す平均組織色を少なくとも実質的に示し、前記指示的反射部分は、関心のありそうなオブジェクトの組織色の少なくとも部分的に周期的な脈動を少なくとも実質的に示し、前記妨害反射成分は、妨害鏡面反射を少なくとも実質的に示す、請求項1に記載のデバイス。
- 前記の推定された指示的反射インデックス要素は、前記基本的反射インデックス要素の推定を有する表現及び予め決められた例示的な指示的反射インデックス要素の正規化された表現の考察の下で得られる、請求項1に記載のデバイス。
- 前記の推定された指示的反射インデックス要素は、前記基本的反射インデックス要素及び前記妨害反射インデックス要素に及ぶ信号面に対して直交すると見なされ、前記の推定された妨害反射インデックス要素は、関心のありそうなオブジェクトに作用する電磁放射供給源の特性を示す輝度インデックス要素により少なくとも実質的に配列される、請求項1に記載のデバイス。
- 前記の推定された基本的反射インデックス要素は、関心のありそうなオブジェクトの指示的エリアの特性信号の凝集された主成分から導出され、前記の推定された妨害反射インデックス要素は、関心のありそうなオブジェクトに作用する電磁放射供給源の特性を示す輝度インデックス要素により少なくとも実質的に配列され、当該デバイスは、関心のありそうなオブジェクトを観察するときに、前記基本的反射インデックス要素の最初の決定を適用するように更に適合される、請求項1に記載のデバイス。
- 前記の推定された指示的反射インデックス要素は、前記基本的反射インデックス要素に対して少なくとも実質的に直交し、前記妨害反射インデックス要素に対して少なくとも実質的に直交するもの見なされる、請求項7に記載のデバイス。
- 対数関数的信号表現は、規定された成分表現に適用され、前記基本的反射係数と前記妨害反射係数との間の想定された割合が、前記指示的反射係数の決定のために考慮される、請求項1に記載のデバイス。
- 基本的変化は、前記指示的反射係数の平均値がゼロより大きくなることを保証するように、規定された成分表現に適用され、前記バイタル信号は、前記指示的反射係数の変換された表現から導出され得る、請求項9に記載のデバイス。
- 前記基本的反射インデックス要素の考察の下での前記特性信号の正規化を更に有する、請求項1に記載のデバイス。
- 前記妨害反射インデックス要素は、関心のありそうなオブジェクトに作用する電磁放射供給源の特性を示す輝度インデックス要素に基づいて得られた前記基本的反射インデックス要素の予め決められた正規化された値を有する表現の考察の下で得られる、請求項1に記載のデバイス。
- 前記特性信号並びに規定された期間の間の前記基本的反射係数、前記指示的反射係数及び前記妨害反射係数の生ずる表現をバッファリングするためのバッファ手段を更に有し、前記処理手段は、支配的な周波数成分を決定するために前記のバッファリングされた指示的反射係数に対してフーリエ変換を適用するように更に適合され、前記基本的反射係数、前記指示的反射係数及び前記妨害反射係数は、それぞれの係数の内積を有する変換された信号表現に変換され、前記指示的反射インデックス要素は、前記の変換された信号表現の考察の下で決定される、請求項1に記載のデバイス。
- 検出された特性信号から情報を抽出するための方法であって、
当該方法は、
オブジェクトにより反射された電磁放射から導出可能なデータストリームを受信するステップであって、前記データストリームは、連続又は別個の時系列特性信号を有し、前記特性信号は、信号空間に関連付けられた少なくとも2つの主成分を有し、前記信号空間は、前記特性信号を表すための相補的なチャネルを有し、前記特性信号の主成分は、前記信号空間のそれぞれの相補的なチャネルに関連する、ステップと、
線形代数方程式を特定するために実質的な線形代数的信号組成モデルの考察の下で規定された成分表現に対して前記特性信号をマッピングするステップであって、前記特性信号は、生理学的反射成分及び妨害反射成分を少なくとも実質的に含み、前記生理学的反射成分は、観察されるべきオブジェクトの基本的生理学的特性に対して起因性のある基本的反射部分と、少なくとも1つの少なくとも部分的に周期的なバイタル信号に対して起因性のある指示的反射部分とを有する、ステップとを有し、
前記基本的反射部分は、基本的反射インデックス要素及び基本的反射係数を有し、
前記指示的反射部分は、指示的反射インデックス要素及び時間依存性指示的反射係数を有し、前記指示的反射係数は、関心のありそうな前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を示し、
前記妨害反射成分は、妨害反射インデックス要素及び妨害反射係数を有し、
当該方法は、
前記基本的反射インデックス要素、前記指示的反射インデックス要素及び前記妨害反射インデックス要素に基づいて前記線形代数方程式を少なくとも部分的に解くことにより、前記少なくとも1つの少なくとも部分的に周期的なバイタル信号を示す計算された表現を導出するステップであって、前記計算された表現は、指示的反射係数を少なくとも有する、ステップを有する、方法。 - コンピュータプログラムがコンピュータ上で実行されたときに、前記コンピュータが請求項14に記載の方法のステップを実行することをもたらすためのプログラムコード手段を有する、コンピュータプログラム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161540596P | 2011-09-29 | 2011-09-29 | |
US61/540,596 | 2011-09-29 | ||
PCT/IB2012/054734 WO2013046082A2 (en) | 2011-09-29 | 2012-09-12 | Distortion reduced signal detection |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014531934A JP2014531934A (ja) | 2014-12-04 |
JP6059726B2 true JP6059726B2 (ja) | 2017-01-11 |
Family
ID=47046662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014532513A Active JP6059726B2 (ja) | 2011-09-29 | 2012-09-12 | 歪み低減された信号検出 |
Country Status (7)
Country | Link |
---|---|
US (1) | US9241674B2 (ja) |
EP (1) | EP2745240A2 (ja) |
JP (1) | JP6059726B2 (ja) |
CN (1) | CN103890780B (ja) |
IN (1) | IN2014CN02002A (ja) |
RU (1) | RU2602975C2 (ja) |
WO (1) | WO2013046082A2 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014064575A2 (en) * | 2012-10-23 | 2014-05-01 | Koninklijke Philips N.V. | Device and method for obtaining vital sign information of a living being |
RU2015156254A (ru) | 2014-05-07 | 2017-07-04 | Конинклейке Филипс Н.В. | Устройство, система и способ выделения физиологической информации |
US11553844B2 (en) | 2014-10-14 | 2023-01-17 | East Carolina University | Methods, systems and computer program products for calculating MetaKG signals for regions having multiple sets of optical characteristics |
WO2016061041A1 (en) | 2014-10-14 | 2016-04-21 | East Carolina University | Methods, systems and computer program products for determining hemodynamic status parameters using signals derived from multispectral blood flow and perfusion imaging |
EP3188651A4 (en) | 2014-10-14 | 2018-07-04 | East Carolina University | Methods, systems and computer program products for visualizing anatomical structures and blood flow and perfusion physiology using imaging techniques |
US10531820B2 (en) | 2014-11-24 | 2020-01-14 | Koninklijke Philips N.V. | Device, system and method for determining the concentration of a substance in the blood of a subject |
US20170035333A1 (en) * | 2015-08-04 | 2017-02-09 | Osram Opto Semiconductors Gmbh | Apparatus for Measuring a Vital Function of a Living Creature |
US10595785B2 (en) * | 2015-10-01 | 2020-03-24 | Silicon Laboratories Inc. | Plethysmography heart rate monitor noise reduction using differential sensors |
RU2018116893A (ru) * | 2015-10-06 | 2019-11-07 | Конинклейке Филипс Н.В. | Устройство, система и способ получения информации о показателе жизненно важной функции живого существа |
GB201608170D0 (en) * | 2016-05-10 | 2016-06-22 | Isis Innovation | A method of determining the frequency of a periodic physiological process of a subject, and a device and system for determining the frequency |
US10335045B2 (en) | 2016-06-24 | 2019-07-02 | Universita Degli Studi Di Trento | Self-adaptive matrix completion for heart rate estimation from face videos under realistic conditions |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6701170B2 (en) * | 2001-11-02 | 2004-03-02 | Nellcor Puritan Bennett Incorporated | Blind source separation of pulse oximetry signals |
BR0304716A (pt) * | 2002-05-08 | 2004-07-13 | Yeda Res & Dev | Método para formação de imagens de mri-bold on line, sensibilizada |
JP4196209B2 (ja) * | 2003-06-30 | 2008-12-17 | 日本光電工業株式会社 | 信号処理方法及びそれを適用したパルスフォトメータ |
US7072709B2 (en) * | 2004-04-15 | 2006-07-04 | Ge Medical Information Technologies, Inc. | Method and apparatus for determining alternans data of an ECG signal |
US8154612B2 (en) * | 2005-08-18 | 2012-04-10 | Qualcomm Incorporated | Systems, methods, and apparatus for image processing, for color classification, and for skin color detection |
US8139850B2 (en) | 2008-12-05 | 2012-03-20 | Tandent Vision Science, Inc. | Constraint generation for use in image segregation |
DE202010017895U1 (de) | 2009-03-06 | 2013-01-08 | Koninklijke Philips Electronics N.V. | System zur Verarbeitung von Bildern von mindestens einem Lebewesen |
US8634591B2 (en) | 2009-08-20 | 2014-01-21 | Koninklijke Philips N.V. | Method and system for image analysis |
JP5834011B2 (ja) | 2009-10-06 | 2015-12-16 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 生体中の少なくとも周期的現象を表す成分を含む信号を処理する方法およびシステム |
JP5856960B2 (ja) | 2009-10-06 | 2016-02-10 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 第1の信号の少なくとも一つの周期的成分を特徴付けるための分析のため第1の信号を得るための方法及びシステム |
US9124777B2 (en) | 2011-01-05 | 2015-09-01 | Koninklijke Philips N.V. | Device and method for extracting information from characteristic signals |
-
2012
- 2012-09-12 CN CN201280048064.8A patent/CN103890780B/zh active Active
- 2012-09-12 EP EP12775311.9A patent/EP2745240A2/en not_active Withdrawn
- 2012-09-12 WO PCT/IB2012/054734 patent/WO2013046082A2/en active Application Filing
- 2012-09-12 JP JP2014532513A patent/JP6059726B2/ja active Active
- 2012-09-12 US US14/239,798 patent/US9241674B2/en active Active
- 2012-09-12 RU RU2014117212/08A patent/RU2602975C2/ru not_active IP Right Cessation
-
2014
- 2014-03-14 IN IN2002CHN2014 patent/IN2014CN02002A/en unknown
Also Published As
Publication number | Publication date |
---|---|
RU2602975C2 (ru) | 2016-11-20 |
WO2013046082A2 (en) | 2013-04-04 |
EP2745240A2 (en) | 2014-06-25 |
US20140213861A1 (en) | 2014-07-31 |
US9241674B2 (en) | 2016-01-26 |
CN103890780A (zh) | 2014-06-25 |
JP2014531934A (ja) | 2014-12-04 |
IN2014CN02002A (ja) | 2015-05-29 |
WO2013046082A3 (en) | 2013-06-13 |
RU2014117212A (ru) | 2015-11-10 |
CN103890780B (zh) | 2018-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6059726B2 (ja) | 歪み低減された信号検出 | |
JP6067706B2 (ja) | 歪が低減された信号検出 | |
JP6615176B2 (ja) | 非干渉的な皮膚組織ハイドレーション測定デバイス及び関連方法 | |
CN108471989B (zh) | 用于生成承载对象的生命体征信息的光体积描记图像的设备、系统和方法 | |
RU2589389C2 (ru) | Устройство и способ извлечения информации из характеристических сигналов | |
US10004410B2 (en) | System and methods for measuring physiological parameters | |
EP3664704B1 (en) | Device, system and method for determining a physiological parameter of a subject | |
WO2013038326A1 (en) | Distortion reduced signal detection | |
CN105813564A (zh) | 用于确定对象的生命体征的设备和方法 | |
Fan et al. | Non-contact remote estimation of cardiovascular parameters | |
Karthick et al. | Analysis of vital signs using remote photoplethysmography (RPPG) | |
WO2019145142A1 (en) | Device, system and method for determining at least one vital sign of a subject | |
EP3838128A1 (en) | Device and method for determining a vital sign of a subject | |
Panigrahi et al. | Contactless HR Measurement from Facial Videos Using Alternative Color Spaces with CEEMDAN | |
Chwyl | A Statistical Framework for Non-Contact Heart Rate Estimation via Photoplethysmogram Imaging | |
Balmaekers | Vital Signs Camera |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150820 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160719 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160721 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20161006 |
|
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: 20161117 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161209 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6059726 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |