JP2019510532A - 皮膚検出に対する装置、システム及び方法 - Google Patents

皮膚検出に対する装置、システム及び方法 Download PDF

Info

Publication number
JP2019510532A
JP2019510532A JP2018540408A JP2018540408A JP2019510532A JP 2019510532 A JP2019510532 A JP 2019510532A JP 2018540408 A JP2018540408 A JP 2018540408A JP 2018540408 A JP2018540408 A JP 2018540408A JP 2019510532 A JP2019510532 A JP 2019510532A
Authority
JP
Japan
Prior art keywords
signal
image data
spectral
skin
unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2018540408A
Other languages
English (en)
Other versions
JP6538287B2 (ja
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2019510532A publication Critical patent/JP2019510532A/ja
Application granted granted Critical
Publication of JP6538287B2 publication Critical patent/JP6538287B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0077Devices for viewing the surface of the body, e.g. camera, magnifying lens
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infrared light
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/0295Measuring 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2433Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0075Measuring 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • A61B5/02427Details of sensor
    • A61B5/02433Details of sensor for infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02438Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20216Image averaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30076Plethysmography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30088Skin; Dermal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Cardiology (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Physiology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Hematology (AREA)
  • General Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Image Analysis (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)

Abstract

本発明は、皮膚検出に対する装置、システム及び方法に関する。信頼できる正確な高速の検出を可能にするために、提案された装置は、シーンの画像データを得る入力インタフェース30であって、前記画像データが時系列の画像フレームを有する、当該入力インタフェースと、前記画像データの関心領域からフォトプレチスモグラフィPPG信号を抽出する抽出ユニット31と、前記PPG信号をスペクトル信号に変換する変換ユニット32と、前記スペクトル信号をソートして、記述子を表すソートされたスペクトル信号を得るソートユニット33と、前記記述子に基づいて前記関心領域を生物の皮膚領域として又は非皮膚領域として分類する分類器34とを有する。

Description

本発明は、皮膚検出に対する装置、システム及び方法に関する。
人のバイタルサイン、例えば心拍(HR)、呼吸数(PR)又は動脈血酸素飽和度は、人の現在の健康状態のインジケータとして及び深刻な医療事象の強力な予測因子として機能する。この理由のため、バイタルサインは、入院患者及び外来患者ケア環境、家庭又は他の健康、レジャー及びフィットネス環境において広く監視される。
バイタルサインを測定する1つの方法は、プレチスモグラフィである。プレチスモグラフィは、一般に、器官又は身体部分の体積変化の測定を指し、特に、心拍ごとに対象の身体を通って移動する心臓血管脈波による体積変化の検出を指す。
フォトプレチスモグラフィ(PPG)は、関心領域又は体積の光反射又は透過の時間変化を評価する光学的測定手法である。PPGは、血液が周辺組織より光を吸収し、したがって、心拍ごとの血液体積の変化が、対応して透過又は反射に影響を与えるという原理に基づく。心拍に関する情報のほかに、PPG波形は、呼吸のような他の生理学的現象に起因する情報を有することができる。異なる波長(典型的には赤及び赤外)における透過率及び/又は反射率を評価することにより、血液酸素(又は他の血液気体/物質)飽和度が、決定されることができる。
近年、控えめな測定に対する非接触の遠隔PPG(rPPG)装置(カメラrPPG装置とも称される)が、導入されている。遠隔PPGは、関心対象から離れて配置された光源、又は一般に放射線源を使用する。同様に、検出器、例えばカメラ又は光検出器も、関心対象から離れて配置されることができる。したがって、遠隔フォトプレチスモグラフィシステム及び装置は、控えめであり、医療及び非医療の日常的な使用に対して適切であると見なされる。この手法は、特に、極めて虚弱な皮膚を持つ新生児集中治療室(NICU)の患者、未熟児、又は広範囲熱傷を持つ患者のような、バイタルサイン監視を必要とする極端な皮膚感度を持つ患者に対して明白な利点を持つ。
Verkruysse他, "Remote plethysmographic imaging using ambient light", Optics Express, 16(26), 22 December 2008, pp. 21434-21445は、フォトプレチスモグラフィ信号が、赤色、緑色及び青色チャネルを使用して、従来の消費者レベルのビデオカメラ及び環境光を使用して遠隔で測定されることができることを示している。
完全に非接触である利点から離れて、カメラ(一般に撮像装置と称される)は、マルチスポット及び大面積測定を可能にする2D情報を提供し、しばしば、追加のコンテキスト情報を含む。特定の測定点/領域における正しい配置に依存する接触センサとは異なり、rPPG手法を使用してパルス信号を測定するのに使用される領域は、実際の画像から決定される。したがって、いかなる照明条件下でも信頼できる皮膚領域の正確な検出は、カメラベースのrPPG装置及び方法の処理チェーンにおいて重要な部分になる。
現在、皮膚領域の信頼できる検出及び追跡に対して既知である2つの主要なアプローチが存在する。
1つのアプローチは、皮膚色(RGBベース)検出及びセグメンテーションに基づく。このアプローチによる方法は、皮膚色を持つ領域の検出及び追跡の両方において高速である。しかしながら、これらは、皮膚領域から反射された光の色を変更する環境光色の変化に対してロバストではなく、低照明条件下又は暗闇で皮膚領域を検出することができない。更に、このような方法は、常に、同じ色を持つ他の表面から皮膚を区別することができない。
他のアプローチは、抽出されたPPG信号に基づく(PPGベース)。このアプローチによる方法は、実際の皮膚領域と同じ皮膚色の他の対象の領域とを区別する際に、よりロバストである。このアプローチは、より強力なPPG信号(最も周期的な信号)を持つ皮膚領域をセグメント化するのにも使用されることができる。しかしながら、このアプローチの信頼性は、PPG信号抽出のロバスト性に依存し、したがって、対象の運動及び血流レベルにより影響される。したがって、パルス信号が、周期的ではない又は弱い場合、カメラベースのシステムは、皮膚領域を検出/セグメント化するのが困難である。更に、このアプローチは、計算的に高価である。
皮膚領域の検出が、rPPG技術に基づくバイタルサイン検出の分野のみの関心ではなく、他の技術分野、例えばプレイヤのジェスチャを認識するのにカメラ技術を使用する遠隔ゲーム応用、顔検出、セキュリティ(監視カメラを使用する人のロバストな検出及びマスクを着用している人の検出又はカメラレジストレーションにおける現実的なマスクからの実際の顔の区別)等においても関心があることは、注意されるべきである。
取得された画像における周期的な色変化自体が、生物の皮膚と非皮膚とを区別するのに使用されることができることが、示されている。しかしながら、特に低い光レベルで又は暗闇での応用、例えば夜間患者監視に対して、特に近赤外(NIR)波長範囲における、これらの吸収変化の比較的低い振幅は、依然として問題であり、皮膚及び非皮膚を区別する実際の試みを失敗させる。
WO2015/049150A1は、対象のバイタルサインに関連した入力信号の処理に対する装置に関し、前記装置は、非侵襲的に検出された入力信号を受信するインタフェースと、前記入力信号の瞬間的な周波数及び/又は振幅表現を含む、前記入力信号の少なくとも1つのフィーチャを抽出するフィーチャ抽出モジュールと、前記少なくとも1つの抽出されたフィーチャに基づいて、前記入力信号に含まれる前記対象のバイタルサインに関する情報を示す、前記入力信号に対する信号情報コンテンツパラメータを決定する処理モジュールと、複数の入力信号の前記信号情報コンテンツパラメータに基づいて前記複数の入力信号を結合して、前記対象の前記バイタルサインを特徴づける結合された出力信号とする結合モジュールとを有する。D1は、更に、対応する方法及び対象のバイタルサインを遠隔で監視する監視システムに関する。
US2013/0296660A1は、血液量パルス、心拍、呼吸波、又は呼吸数のような、対象の1以上の生理学的パラメータを遠隔で測定又は監視する方法及びシステムを記載する。前記方法は、前記対象の一連の画像を取得するステップと、関心のある生理学的パラメータを取得するように前記画像を処理するステップとを含む。これらの方法は、アクティブナイトビジョンカメラから取得された信号を含む、単一チャネル信号を分析するのに使用されることができる。結果として、これらの方法は、日光及び微小光条件の両方で1以上の生理学的パラメータを測定又は監視するのに使用されることができる。
本発明の目的は、特に対象のバイタルサインを検出する装置及び方法において使用する、皮膚の信頼できる正確かつ高速な検出を可能にする装置及び対応する方法並びにシステムを提供することである。
本発明の第1の態様において、
‐シーンの画像データを取得する入力インタフェースであって、前記画像データが、時系列の画像フレームを有する、当該入力インタフェースと、
‐前記画像データの関心領域からフォトプレチスモグラフィPPG信号を抽出する抽出ユニットと、
‐前記PPG信号をスペクトル信号に変換する変換ユニットと、
‐記述子を表すソートされたスペクトル信号を得るように前記スペクトル信号をソートするソートユニットと、
‐前記記述子に基づいて生物の皮膚領域として又は非皮膚領域として前記関心領域を分類する分類器と、
を有する皮膚検出に対する装置が、提示される。
本発明の更なる態様において、対応する方法が、提示される。
本発明の更に他の態様において、
‐シーンの画像データを取得する撮像ユニットと、
‐前記シーンの前記取得された画像データに基づいてここに開示される皮膚検出に対する装置と、
を有する皮膚検出に対するシステムが、提示される。
本発明の更に他の態様において、コンピュータにより実行される場合に、ここに開示される方法のステップを前記コンピュータに実行させるプログラムコード手段を有するコンピュータプログラム、及びコンピュータにより実行される場合に、ここに開示された方法が実行されるようにするコンピュータプログラムを記憶する不揮発性コンピュータ可読記録媒体が、存在する。
本発明の好適な実施例は、従属請求項において規定される。請求された方法、システム、コンピュータプログラム及び媒体が、請求された装置及び従属請求項に規定されたものと同様の及び/又は同一の好適な実施例を持つと理解されるべきである。
本発明は、例えば生物の皮膚領域又は完全な生物又は皮膚領域及び非皮膚領域を示す異なる身体部分を潜在的に含むシーンの画像データから得られたソートされた周波数領域信号においてエネルギの分布を評価することにより、周波数領域においてノイズからパルス信号を区別するアイデアに基づく。特に、(ソートされた(又はランク付けされた)周波数領域表現において)パルス信号のエネルギが、ソートされた(又はランク付けされた)ノイズスペクトルのエネルギより速く降下することが、わかっている。好適には、例えばこのようなランク付けされたスペクトルを入力として使用して、適宜に訓練された分類器を用いて、前記画像データ内の特定の関心領域(ROI)、例えば前記画像データの画像フレームの画像セグメントが、皮膚を示すか否かの信頼できる決定が、なされることができる。
ここで、ソート及びランク付けは、同義語として理解されるべきである。これは、周波数サンプルの順序が、減少する又は増加する振幅で順序付けられたように見えるように変更される、すなわち、前記周波数サンプルが、周波数方向において再順序付けされることを意味する。
したがって、提案された装置、システム及び方法は、実質的に3D時系列画像フレームから1D時間信号(ここで「PPG信号」とも称される)を抽出し、前記時間信号を変換された信号に変換し、記述子を得るように前記変換された信号をソートし、前記記述子に基づいて前記3D時系列画像フレームを分類することにより、3D時系列画像フレーム(特にその中の画像セグメント)が、生きている人間組織を含むか否かを決定することを可能にする。
これに関連して、前記画像データのROIから、例えばビデオシーケンス内のROIから、前記PPG信号を抽出する時間において、前記ROIが実際に皮膚を示すかどうかが、既知ではないことに注意すべきである。この不確定性にかかわらず、前記ROIは、皮膚を含んでいたように扱われ、PPG信号は、従来、皮膚領域に対して行われていたようにこの仮定で前記ROIから抽出される。
前記変換に対して、一般的に、FFT、DCT、ウェーブレット変換、ハール変換等のようないかなる種類の周波数領域変換も、使用されうる。
一実施例によると、前記変換ユニットは、前記PPG信号を位相情報なしのスペクトル信号に、特にパワースペクトル又は絶対スペクトルに変換する、すなわち複素スペクトルを実信号に変換するように構成される。フーリエ変換の出力は、複素(実部及び虚部)信号である。前記位相情報に関心が無い場合、複素サンプルのノルム(複素面内のベクトルの長さ)が取られてもよく、これは、しばしば、絶対スペクトル又は振幅スペクトルと称される。これらの数を二乗することにより、パワースペクトルが得られる。これは、(単純に未知の位相及び周波数を持つ正弦曲線でモデル化されることができるが、前記ソートされたスペクトルにおいて常に第1の(少数の)ビンにおける高い振幅として表れる)パルスと(ランダム(振幅及び位相の両方)信号としてモデル化されることができ、前記(ソートされた又はソートされていない)スペクトルにおいて同じ振幅を持つ広帯域の周波数として表れる)ノイズとを区別することが望ましいので、有利である。これら2つのソートされたスペクトルは、分類するのが容易である。
他の実施例において、前記ソートユニットは、前記スペクトル信号を2以上のスペクトルサブ信号に分割し、前記記述子を表すソートされたスペクトルサブ信号を得るように前記サブ信号を別々にソートするように構成される。このようなサブ信号が、例えばノイズからパルス信号を、したがって非皮膚から皮膚を区別することを更に良好に可能にすることがわかっている。
好ましくは、前記ソートユニットは、前記スペクトル信号を、前記スペクトル信号の第1の周波数範囲をカバーする帯域内サブ信号及び前記スペクトル信号の残りの周波数範囲をカバーする帯域外サブ信号に分割し、前記記述子を表すソートされた帯域内サブ信号及びソートされた帯域外サブ信号を得るように前記帯域内サブ信号及び前記帯域外サブ信号を別々にソートするように構成される。これにより、前記ソートユニットは、前記帯域内サブ信号が、前記スペクトル信号の周波数範囲の下側部分又は前記スペクトル信号の最高周波数ピークの周りの周波数範囲の部分をカバーするように前記スペクトル信号を分割するように更に構成されてもよい。例えば、前記帯域内は、単純にオーバサンプリングされたPPG信号の下半分であることができ、前記帯域外は、単純に同じオーバサンプリングされたPPG信号の上半分であることができる。より洗練されたバージョンにおいて、前記帯域内は、例えば合計周波数ビンの半分を持つ、最高周波数ピークの周りの窓として規定されてもよく、前記帯域外は、この場合、残りの周波数ビンにより規定されてもよい。前記周波数ビンを分離する他の選択肢が、例えば見込みのあるパルススペクトル成分(例えば前記周波数領域における20%の最高ピーク)を帯域内として、残りのスペクトル成分を帯域外として選択することにより、同様に、より正確な帯域内を規定することが可能である。
他の実施例において、前記変換ユニットは、前記スペクトルの(正規化の種類に依存する)エネルギ又は振幅が、後の評価に影響を持たないように、前記スペクトル信号を正規化するように構成される。これは、前記スペクトルを、実質的にそれぞれ皮膚及び非皮膚領域におけるパルス及びノイズ振幅から独立したものにする。前記正規化は、例えば、L1又はL2ノルムにより実施されてもよい。実際的な実施例は、位相スペクトルをすべての周波数において一定値にリセットし、ロバストな分類を得るように振幅スペクトルを正規化することを有してもよい。
前記装置は、2以上の反復を実行するように前記変換ユニット及び前記ソートユニットを制御する制御ユニットを更に有してもよく、前記ソートユニットから出力される前記ソートされたスペクトル信号は、次の反復において前記変換ユニットに対する入力PPG信号として使用される。したがって、前記反復は、前の反復の出力に対して作用し、好ましくは、前記信号の長さを半分にした後に、前記分類器により最終的に使用される前記スペクトルのマルチスケール表現を生じる。特に、最初の反復は、ノイズ周波数成分の位相をリセットすることにより非組織を表す比較的先鋭化した信号を生じる。好ましくは、前記スペクトルの長さは、異なる反復において減少され、すなわち前記変換は、異なるスケールで実行され、結果としてマルチスケール表現を生じ、(より少ないサンプルを持つ)より粗いスペクトルが、依然として、同じ信号を記述するが、より粗いスケールである。
前記変換及びソートの繰り返しが、(先鋭化された)パルススペクトルを平坦なスペクトルに変換するが、(比較的平坦な)ノイズスペクトルが、先鋭化されたスペクトルに変換されることがわかっている。繰り返される変換において、前記位相情報、前記ピークの場所、及び前記スペクトルの振幅は、好ましくは、無視される。前記ソートは、前記ピークが同じ場所にあることを確認し、前記分類を改善する。
他の実施例において、前記分類器は、各反復において前記ソートユニットから出力された前記ソートされたスペクトル信号を連結し、前記関心領域を生物の皮膚領域として又は非皮膚領域として分類するのに、前記連結されたソートされたスペクトル信号を記述子として使用するように構成される。前記連結されたソートされたスペクトル信号において、前記先鋭化された部分及び平坦な部分は、交互であり、前記交互のシーケンスは、非皮膚領域からのPPG信号と比較して皮膚領域からのPPG信号に対して反対である。したがって、非常に信頼できる分類は、このような前記連結されたソートされたスペクトル信号に基づいてなされることができる。
更に、一実施例において、前記抽出ユニットは、画像フレームごとに前記画像データの画素のグループの画像データ値を結合し、特に平均し、前記結合された画像データ値から前記PPG信号を得るように構成される。したがって、画素のグループにわたる平均(結合)が、取られてもよく、連結された平均の時間発展が、前記1D時間信号として見なされてもよい。これは、前記分類の信頼性を改善する。基本的に、(皮膚マップの)解像度の損失は、平均化が(前記パルス信号のノイズを少なくすることにより)パルスとノイズとの間の差を増加するので、前記分類の改善された信頼性に対して取引される。
前記抽出ユニットは、波長又は波長範囲において画像フレームごとに前記画像データの画素のグループの画像データ値を結合し、特に平均し、前記結合された画像データ値から前記PPG信号を得るように更に構成されてもよい。更に、前記抽出ユニットは、画素又は画素のグループごとに及び画像フレームごとに、少なくとも2つの異なる波長チャネルの画像データ値を重み付け平均として結合し、前記結合された画像データ値から前記PPG信号を得るように構成されてもよい。したがって、波長(又は色)ごとの画素のグループにわたる平均が、取られてもよく、連結された平均の時間発展は、前記1D時間信号を得るように波長ごとに連結されてもよい。これにより、前記抽出ユニットは、正規化された血液量パルスベクトルシグネチャベースの方法(すなわちPbv法)、クロミナンスベースの方法(すなわちCHROM法)、ブラインド音源分離法(すなわちBSS法)、主成分分析(PCA)又は独立成分分析(ICA)を使用して前記重みを計算するように構成されてもよい。
一般に、PPG信号は、皮膚の血液量の変化から生じる。したがって、前記変化は、反射/透過光の異なるスペクトル成分において見られる場合に特徴的な拍動性「シグネチャ」を与える。この「シグネチャ」は、基本的に、血液の吸収スペクトル及び血液のない皮膚組織の吸収スペクトルのコントラスト(差)として生じる。検出器、例えばカメラ又はセンサが、光スペクトルの特定の部分を各々感知する離散的な数の色チャネルを持つ場合、これらのチャネルにおける相対的な拍動性は、「正規化された血液量ベクトル」Pbvとも称される、「シグネチャベクトル」に配置されることができる。このシグネチャベクトルが既知である場合に、色チャネル及びシグネチャベクトルに基づく運動ロバストパルス信号抽出が可能であることは、参照によりここに組み込まれる、 G. de Haan and A. van Leest, "Improved motion robustness of remote-PPG by using the blood volume pulse signature", Physiol. Meas. 35 1913, 2014に示されている。前記パルス信号の品質に対して、前記シグネチャが正しいにもかかわらず、他の形で、既知の方法が、前記シグネチャベクトルにより示される正規化された色チャネルを用いる前記パルスベクトルの所定の補正を達成するために出力パルス信号内にノイズを混合することは、不可欠である。
Pbv法及び(「基準生理学情報を示す設定された向きを持つ所定のインデックス要素」と称される)正規化された血液量ベクトルの使用の詳細は、その詳細が参照によりここに組み込まれるUS2013/271591A1にも記載されている。
前記分類器は、前記関心領域が生物の皮膚領域である尤度を決定するように構成されてもよく、すなわち前記分類器は、前記関心領域が皮膚領域であるか否かのバイナリ決定のみならず、前記関心領域が皮膚領域である尤度をも発行する。
前記分類器は、管理学習(supervised learning)(例えばAdaBoost、SVM等)から得られてもよく、入力として前記変換された信号のサンプル(例えばランク付けされた正規化された周波数ビン)を取り、画像セグメントが生きている人間組織であるか否かの尤度を識別する信号(硬い(バイナリ)ラベル、又は回帰値)を出力する。更に、前記分類器は、生きている人間組織を含むセグメントを表す変化する振幅、ノイズのレベル、及び心拍帯域の周波数を持つ正弦曲線を含む1D時間信号と、生きている人間組織を含まないセグメントを表すノイズ信号とのデータセットを使用して訓練されてもよい。
前記装置は、前記画像データの画像フレームをセグメント化するセグメンテーションユニットを更に有してもよく、前記抽出ユニットは、別々の後の処理のために2以上の画像フレームセグメントからPPG信号を抽出するように構成される。
本発明は、好ましくは、遠隔PPG技術の使用によるバイタルサイン取得に関連して使用される。この目的で、前記撮像ユニットは、好ましくは、経時的に前記シーンの画像のシーケンスを取得するように構成され、前記装置は、前記画像のシーケンス内の検出された皮膚領域からの画像情報に基づいて前記シーン内の対象のバイタルサインを検出する場タイルサイン検出器を更に有してもよい。したがって、皮膚領域の提案された検出は、バイタルサインの取得中に皮膚領域を検出及び/又は追跡するのに一度又は連続的に使用されてもよい。
本発明のこれら及び他の態様は、以下に記載される実施例を参照して説明され、明らかになる。
本発明によるシステムの第1の実施例の概略図を示す。 本発明による装置の第1の実施例の概略図を示す。 異なる生物に対する皮膚及び非皮膚領域からの典型的なPPG信号を示す図を示す。 異なる生物に対する皮膚及び非皮膚領域からの典型的なPPG信号を示す図を示す。 本発明の一実施例による方法の様々なステップにおける典型的な信号を示す図を示す。 本発明の一実施例による方法の様々なステップにおける典型的な信号を示す図を示す。 本発明の一実施例による方法の様々なステップにおける典型的な信号を示す図を示す。 本発明の一実施例による方法の様々なステップにおける典型的な信号を示す図を示す。 本発明による装置の第2の実施例の概略図を示す。 第2の実施例による処理に関する記述子を示す図を示す。 本発明の第3の実施例による処理に関する記述子を示す図を示す。 本発明の第3の実施例による処理に関する記述子を示す図を示す。
図1は、皮膚検出に対する装置12を含む本発明によるシステム10の概略図を示す。システム10及び装置12は、好ましくは、対象14の時系列の画像フレームを含む画像データから前記対象のバイタルサインを検出する装置及び方法において使用されうる。対象14、本例において患者は、例えば病院又は他のヘルスケア施設内の、ベッド16に横たわるが、例えば保育器内に横たわる新生児若しくは早産児、又はスポーツを行っているアスリートのような、家庭若しくは異なる環境における人であってもよい。
撮像ユニット18は、前記シーンの画像データ(RGB画像とも称され、可視及び/又は赤外光の波長範囲における画像として理解されるべきである)を取得する、特に、好ましくはPPG信号が得られることができる対象14の皮膚領域15を含む、経時的な対象14の画像フレームのシーケンスを取得するカメラ(検出ユニット、又はカメラベース若しくは遠隔PPGセンサとも称される)を含んでもよい。対象14のバイタルサインを得る装置12の使用において、皮膚領域15は、好ましくは、頬又は額のような顔の領域であるが、手又は腕のような可視の皮膚表面を持つ身体の他の領域であってもよい。
前記撮像により獲得された画像フレームは、特に、例えば(デジタル)カメラ内のアナログ又はデジタル光センサを用いて獲得されたビデオシーケンスに対応してもよい。このようなカメラは、通常、特定のスペクトル範囲(可視、nIR)において動作してもよく、又は異なるスペクトル範囲に対する情報を提供し、特にPPG信号の抽出を可能にしてもよい、CMOS又はCCDセンサのような光センサを含む。前記カメラは、アナログ又はデジタル信号を提供してもよい。前記画像フレームは、関連付けられた画素値を持つ複数の画像画素を含む。特に、前記画像フレームは、光センサの異なる光感知素子で獲得された光強度値を表す画素を含む。これらの光感知素子は、特定のスペクトル範囲(すなわち特定の色を表す)において感知してもよい。前記画像フレームは、人の皮膚部分を表す少なくとも一部の画像画素を含む。これにより、画像画素は、光検出器の1つの光感知素子及びその(アナログ又はデジタル)出力に対応してもよく、又は複数の前記光感知素子の(例えばビニングによる)組み合わせに基づいて決定されてもよい。
カメラ18を使用する場合、システム10は、更に、オプションとして、例えば所定の波長範囲又は複数の範囲(例えば赤、緑及び/又は赤外波長範囲)の光で、患者の顔の皮膚(例えば頬又は額の一部)のような、関心領域24を照射する/光を当てる、ランプ又はLEDのような、照射ユニット22(照射源又は光源又は電磁放射体とも称される)を有してもよい。前記照射に応答して関心領域24から反射された光は、カメラ18により検出される。他の実施例において、専用の光源は提供されないが、環境光が、対象14の照射に対して使用される。前記反射光から、所望の波長範囲(例えば緑及び赤又は赤外光、又は少なくとも2つの波長チャネルをカバーする十分に大きな波長範囲の光)の光のみが、検出及び/又は評価されてもよい。
装置12は、更に、前記決定された情報を表示する及び/又は装置12、カメラ18、照射ユニット22及び/又はシステム10の他のパラメータを変更するインタフェースを医療専門家に提供するインタフェース20に接続される。このようなインタフェース20は、異なるディスプレイ、ボタン、タッチスクリーン、キーボード又は他の人間機械インタフェース手段を有してもよい。
図1に示されるシステム10は、例えば、病院、ヘルスケア施設、又は高齢者ケア施設等に配置されうる。患者の監視から離れて、本発明は、新生児監視、一般的な監視応用、セキュリティ監視又はフィットネス機器、ウェアラブル、若しくはスマートフォンのような携帯装置のようないわゆるライフスタイル環境のようなほかの分野において使用されてもよい。装置12、カメラ18及びインタフェース20の間の一方向又は二方向通信は、無線又は有線通信インタフェースを介して機能しうる。本発明の他の実施例は、スタンドアローンで提供されないが、カメラ18又はインタフェース20内に一体化された装置12を含んでもよい。
図2は、図1に示されるシステム10内の装置12として使用されてもよい、本発明による装置12aの第1の実施例の概略図を示す。対象14の1以上のバイタルサインを得るために、前記対象の皮膚領域は、前記画像データ内で見つけられなければならない。この目的で、提案された装置12aは、シーンの画像データ40を得るための入力インタフェース30を有し、前記画像データは、撮像ユニット18により取得された時系列の画像フレームを有する。抽出ユニット31は、前記画像データの関心領域からPPG信号41を抽出し、前記関心領域は、1以上の画像フレームのセグメンテーションから生じる単一の画素又は画素のグループ又は領域でありうる。変換ユニット32は、前記PPG信号41をスペクトル信号42に変換する。ソートユニット33は、記述子を表すソートされたスペクトル信号43を得るように前記スペクトル信号42をソートする。最終的に、分類器34は、前記記述子に基づいて前記関心領域を生物の皮膚領域として又は非皮膚領域として分類し、バイナリ決定(例えば前記関心領域が皮膚であるか否かの標示)又は前記関心領域が皮膚であるか否かの尤度でありうる、対応する分類結果44を発行する。
ユニット30乃至34は、専用のハードウェア要素として構成されなくてもよいが、適宜にプログラムされるプロセッサ又はコンピュータとして構成されてもよい。装置12aは、例えば共通のハウジング(例えば撮像ユニット18の共通のハウジング)内に、全ての要素及びユニットを含む一体化された装置として、又は前記要素及びユニットが分散されている、すなわち異なる位置に配置された別々の要素及びユニットとして実装される、図1に示されるような、分散された装置として構成されてもよい。
図3A及び3Bは、異なる生物に対する皮膚及び非皮膚領域からの典型的なPPG信号41を示す図を示す。特に、図3Aは、第1の対象(例えば大人)からのパルス信号41a及び第2の対象(例えば新生児)からのパルス信号41cを示す。見られることができるように、異なる対象からのパルス信号は、異なる。更に、単一の対象からのパルス信号でさえ、(例えば位相において)時間変化する。前記パルスが、一般に、明確な周期的成分を持つが、振幅、位相及び周波数(例えば心拍可変性)において変化が存在してもよく、典型的には、前記信号は、また、センサノイズを被り、対象運動により歪められうる。図3Bは、前記第1の対象(例えば大人)からのノイズ信号41b及び前記第2の対象(例えば新生児)からのノイズ信号41dを示す。見られることができるように、ノイズ信号は、学習されることができない不規則な/一貫性のない信号であり、パルス信号も、大幅な可変性を示す。したがって、1つのアイデアは、得られたPPG信号を、管理学習を可能にする異なる表現に変換することである。
パルス及びノイズを2つのクラスとして見なして、前記変換された表現(例えば前記記述子)は、PPG信号の3つの特性を除去すべきである。前記記述子は、パルス/ノイズの位相変化、すなわち異なる瞬間におけるパルスに対して不変であるべきである。更に、前記記述子は、パルス/ノイズの振幅に依存すべきではない。更に、前記記述子は、パルス/ノイズの変化する周波数、すなわち異なる対象から独立であるべきである。
上記の要件を仮定して、以下の典型的なアプローチは、PPG信号41に対して使用され、これは、本発明による方法の様々なステップにおける典型的な信号を示す図を示す図4の使用により示される。第1のステップにおいて、スペクトルブーストは、変換ユニット32によりPPG信号41(図示のために、パルス信号41a及びノイズ信号41aが図4Aに示される)に対して使用される。パルスが周期的な信号であるという有効な仮定に基づいて、変換ユニット32は、例えばフーリエ変換(FT)の使用により、分析のためにPPG信号41を時間領域から周波数領域に変換する。図4Bに示される変換されたパルス42aは、周波数スペクトルにおいて著しいピークを示すのに対し、図4Bに示される変換されたノイズ42bは、このようなパターンを示さない不規則な信号である。これにより、フーリエ変換(又は高速フーリエ変換;FFT)は、離散コサイン変換(DCT)、離散サイン変換(DST)、ウェーブレット変換、又はハール変換等により置き換えられることができる。
前記FTを使用することは、
Figure 2019510532
として記載されることができ、
Figure 2019510532
に対して、長さL=2n,{n|n∈Z,n?2}を持つPPG信号を用い、
フーリエ変換(例えばL=64);
Figure 2019510532
は、FT演算を示す。
Figure 2019510532
の実部及び虚部は、振幅又はパワースペクトルを使用するだけで除去されることができる、変化する位相情報を含む。好ましくは、前記パワースペクトルは、一般的に、パルスの周波数ピークをブーストし、前記ノイズを抑制するので、使用される。
Figure 2019510532
は、半分の冗長性を持つミラーされたスペクトルであるので、前記パワースペクトルを得る前に半分にされ、すなわち
Figure 2019510532
であり、ここでconj(.)は、結合(conjugation)を示し、
Figure 2019510532
は、要素ごとの積を示す。SL/2において、前記位相情報は消え、前記パルスの周波数ピークは、図4Bに示されるように、ノイズのものと比べてブースとされる。
次のステップにおいて、スペクトル正規化が、変換ユニット32により実行されうる。この正規化は、前記スペクトルを、それぞれ皮膚/非皮膚領域におけるパルス及びノイズ振幅から実質的に独立にする。前記スペクトル振幅は、依然としてSL/2において可変であり、これらは、
Figure 2019510532
として正規化され、ここで、||.||Pは、Lpノルムを示す。これは、L1ノルム又はL2ノルムのいずれかであることができる。標準偏差の正規化は、絶対エネルギ情報のみが照射されるべきなので、好適ではないが、変化は、パルス及びノイズを区別するために前記スペクトル内に残るべきである。典型的な実施例において、L2ノルムは、合計エネルギに対してノイズを抑制することができるので、使用される。正規化された
Figure 2019510532
は、スペクトル振幅から独立であるのに対し、エントリの相対エネルギ分布は、図4Cに示されるように、正規化されたスペクトルパルス信号42a'及び正規化されたスペクトルノイズ信号42b'に対して残される。
次のステップにおいて、スペクトルソートは、ソートユニット33により実行される。周波数変化のため、
Figure 2019510532
は、まだ分類に対して使用されることができない。しかしながら、異なる個人は、異なる心拍を持つので、これらのパルス周波数は、ほとんど先鋭化され、特定の(低い)帯域、すなわち、[40,240]拍/分(bpm)に集中されるのに対し、背景(非皮膚)信号は、通常、高周波帯域に広がる白色ノイズである。このために、前記
Figure 2019510532
は、再び分割され、例えば(サブ信号とも称される)2つの半分に、好ましくは前記パルス関連特性がここで暗示的に利用される「帯域内」及び[帯域外」周波数を近似するように下部分及び上部分に等分される。
周波数依存性を除去するために、前記分割されたスペクトルは、ソートされ、次いで、
Figure 2019510532
のように連結され、ここで、sort(.)は、例えば振幅/エネルギの降順にスペクトルエントリをソートすることを示す。
Figure 2019510532
において、パルス及びノイズにおける周波数変化は、除去されるが、下帯域及び上帯域における本質的な差は、ソートされたスペクトルパルス信号43a及びソートされたスペクトルノイズ信号43bを示す図4Dに示されるように、保存される。したがって、このステップにおいて、ランク付け(ソート)手順は、本質的に実行され、周波数ビンに対して作用する。
提案されたアプローチと既知のアプローチとの間の本質的な差は、既知のアプローチが、パルス及びノイズを分離するのに単一の値(例えば正規化されたスペクトルピーク)のみを使用するのに対し、提案されたアプローチは、実際には形状記述子である分類に対してソートされたスペクトル内の全てのエントリを利用することである。本質的に、提案されたアプローチによると、(前記ランク付けされた周波数スペクトルにおける)パルス信号のエネルギが、前記ランク付けされたノイズスペクトルのエネルギより速く降下する現象が、利用される。これらのランク付けされたスペクトルを入力として使用する訓練された分類器を用いて、最適な決定が、得られることができる。換言すると、上に示された第1の実施例において、ソートされたスペクトル信号43a、43bは、下の画像データ内のそれぞれの関心領域が生物の皮膚領域であるか又は非皮膚領域であるかを決定するのに分類器34により使用される。
図5は、本発明による装置12bの第2の実施例の概略図を示す。この実施例において、制御ユニット35は、2以上の反復を実行するように前記変換ユニット32及び前記ソートユニット33を制御するために提供され、前記ソートユニット33から出力されたソートされたスペクトル信号43は、次の反復において変換ユニット32に対する入力PPG信号41'として使用される。したがって、本実施例によると、マルチスケール反復が、以下に説明されるように実行されうる。
装置12aの第1の実施例を用いて、変換された信号
Figure 2019510532
は、入力PPG信号
Figure 2019510532
を条件として得られ、ここで、パルス及びノイズは、自己統合(self-unified)されたが、相互に異なる解釈を持つ。パルス及びノイズの記述子が比較される場合、パルス記述子は、目立つフィーチャ(例えば、第1の場所におけるピーク)を持つのに対し、前記ノイズは、そうではない。より良好な分類性能を得るために、異なるクラスからの記述子は、大きなクラス間変化を必要とし、すなわちパルス及びノイズは、容易に区別可能である。これは、手順(ブースト、正規化、ソート)を反復することにより改善されることができる。ここで、比較的平坦なノイズスペクトルが、明白なピークに変換し、先鋭化されたパルススペクトルが、比較的平坦な結果に変換する。結合された2回の反復は、2つのクラスの間の逆位相パターンを提供し、これは、より容易な分離を生じる。
同様に、同じ変換は、更に、前記変換された信号に対して1回以上反復される。パルス及びノイズにおける新しく生成されたパターンは、反対の順序で、すなわち、得られたパルス信号記述子43a'及びノイズ信号記述子43b'に対する図6に示されるように、「ピーク‐フラット‐ピーク」対「フラット‐ピーク‐フラット」で、生じる。図示された例において、5回の反復のソートされたスペクトル信号が、連結される。
このようにして、より長い記述子Xが、異なるスケールで反復的に変換された信号を収集/連結するように作成され、すなわち、
Figure 2019510532
であり、ここで、
Figure 2019510532
は、長さL/(2×i)を持つi番目の反復における変換された信号である。前記反復が終了した場合、完全な記述子は、更に、L2ノルムにより正規化されてもよい。実際に、提案された記述子は、マルチスケール反復が前記記述子の識別性を改善することができるという仮説上で構築される。このような仮説は、実験的に確認されている。
したがって、前の反復の出力に作用する反復は、好ましくは、前記信号の長さを半分にした後に開始される。この場合、前記反復は、前記スペクトルのマルチスケール表現を前記分類に対して利用可能にする。特に、最初の反復は、ノイズ周波数成分における位相の除去のため、非組織を表す比較的先鋭化された信号を生じる。この理由で、少なくとも2つの順次的な変換、すなわち、FFT‐位相削除‐正規化‐ランク付け‐FFT‐位相削除‐正規化‐ランク付けが、実行されてもよく、ここで、効率のため、第2の変換は、第1の反復から得られた半分のスペクトルに作用してもよい。
更に、パルス及びノイズ表現の間の識別性は、更に改善されうる。一実施例において、前記変換された記述子(43a',43b')における平坦な/先鋭化されたパターンは、更に平坦に/先鋭化されることができる。式(2)は、前記パワースペクトルを得るのに単一の信号を使用する。これは、(一方のフレームがシフトされた)2つの時間的に隣接した信号を使用することにより改善されうる。これは、主に、ノイズクラスの利益になり、すなわち、2つのノイズ信号の結合が、前記パワースペクトルの実部において負のエントリを誘発する。これは、ノイズ信号における高周波成分により、すなわち、背景(非皮膚)信号は、ほとんど白色ノイズであり、したがって、ここで利用される。前記スペクトルの最小の負の値を減算することは、図7A(及び図4D)に示される元のノイズ記述子43bと比較して得られたノイズ記述子43cを示す図7Bに示されるように、第1の反復において前記ノイズ記述子を更に平坦にすることができる。したがって、このような実施例において、前記ブーストするステップは、単一の信号の代わりに(1フレームシフトされた)2つの時間的に隣接した信号を結合することにより修正されうる。
上に示されたように、好適な実施例において、前記位相情報を無視するが、前記絶対スペクトル又は前記パワースペクトルを取ることが、提案される。
好ましくは、前記ランク付けは、2回、すなわち、帯域内ランク付け及び帯域外ランク付けが行われる。帯域内は、単純にオーバサンプリングされた信号の下半分であることができ、帯域外は、上半分である。しかしながら、より洗練されたバージョンにおいて、帯域内は、例えば合計ビン数の半分を持つ、最高周波数ピークの周りの窓として規定され、帯域外は、この場合、残りの周波数ビンにより形成される。
前記分類は、管理学習(例えばAdaBoost、SVM等)から得られ、入力として前記変換された信号(例えば位相情報を持たない、ランク付けされた、正規化された周波数ビン)のサンプルを取り、画像セグメントが生きている人間の組織であるか否かの尤度を識別する信号(硬い(バイナリ)ラベル、又は回帰値)を出力する分類器を使用してもよい。前記管理学習は、皮膚及び非皮膚表面から得られた実際のデータを使用しうるが、良好な性能は、変化する振幅、ノイズのレベル、及び生きている人間の組織を含むセグメントを表す脈拍数帯域の周波数を持つ正弦曲線を含む1D時間信号と、生きている人間の組織を含まないセグメントを表すノイズ信号(ゼロ平均ガウス、又は一様等)のデータセットを使用して前記分類器を訓練することにより得られている。
更に、提案された方法は、セグメンテーションから得られた画像領域を分類するのに使用されてもよく、場合により、運動追跡が、連続した画像において経時的に個別のセグメントを追跡するのに使用されてもよい。
本発明は、好ましくは、人のバイタルサインの取得に対するrPPGの分野において使用される。したがって、撮像ユニットにより得られた画像は、上で説明されたように皮膚領域を検出するのに使用されるのみならず、検出された(及び好ましくは本発明の使用により追跡された)皮膚領域から、心拍、SpO2等のような人のバイタルサインを得るのに使用されるPPG信号が、得られる。撮像ユニット18は、少なくとも前記シーンが(環境光により及び/又は照射により)照射される波長又は波長範囲において感知するが、特に所望のバイタルサインを得るのに必要とされる場合に、同様に他の波長に対して感知してもよい。
本発明の他の実施例において、皮膚検出に対する提案された分析は、一般的に既知であるように皮膚検出に対する他の方法、例えば皮膚領域から反射された構造化された光の時間的拍動性又はクロミナンスの分析と組み合わせられることができる。前記方法は、更なるステップを有してもよく、前記装置の様々な実施例に対して上で説明されたように及びここに開示されるように修正されてもよい。
提案された装置及び方法は、PPG関連バイタルサイン(例えば心拍、SpO2、呼吸)の連続的な目立たない監視に対して使用されることができ、NICU、手術室又は一般病棟において使用されることができる。提案された装置及び方法は、個人のヘルス監視に対して使用されることもできる。一般に、本発明は、皮膚が、シーンの画像において検出される必要があり、特に非皮膚から区別される必要がある全ての応用において使用されることができる。
本発明は、図面及び先行する記載において詳細に図示及び記載されているが、このような図示及び記載は、例示的又は典型的であり、限定的ではないと見なされるべきであり、本発明は、開示された実施例に限定されない。開示された実施例に対する他の変形例は、図面、開示及び添付の請求項の検討から、請求された発明を実施する際に当業者により理解及び達成されることができる。
請求項において、単語「有する」は、他の要素又はステップを除外せず、不定冠詞「a」又は「an」は、複数を除外しない。単一の要素又は他のユニットが、請求項に記載された複数のアイテムの機能を満たしてもよい。特定の方策が相互に異なる従属請求項に記載されているという単なる事実は、これらの方策の組み合わせが有利に使用されることができないことを示さない。
コンピュータプログラムは、他のハードウェアの一部として又は一緒に供給される光記憶媒体又は半導体媒体のような、適切な不揮発性媒体に記憶/分配されてもよいが、インターネット又は他の有線若しくは無線電気通信システムのような、他の形で分配されてもよい。
請求項内の参照符号は、その範囲を限定すると解釈されるべきではない。

Claims (15)

  1. シーンの画像データを得る入力インタフェースであって、前記画像データが、時系列の画像フレームを有する、当該入力インタフェースと、
    前記画像データの関心領域からフォトプレチスモグラフィPPG信号を抽出する抽出ユニットと、
    前記PPG信号をスペクトル信号に変換する変換ユニットと、
    前記スペクトル信号をソートして、記述子を表すソートされたスペクトル信号を得るソートユニットと、
    前記記述子に基づいて前記関心領域を生物の皮膚領域として又は非皮膚領域として分類する分類器と、
    を有する皮膚検出に対する装置。
  2. 前記変換ユニットが、前記PPG信号を、位相情報なしのスペクトル信号に、特にパワースペクトル又は絶対スペクトルに変換する、
    請求項1に記載の装置。
  3. 前記ソートユニットが、前記スペクトル信号を、2以上のスペクトルサブ信号に分割し、前記サブ信号を別々にソートして、前記記述子を表すソートされたスペクトルサブ信号を得る、
    請求項1に記載の装置。
  4. 前記ソートユニットが、前記スペクトル信号を、前記スペクトル信号の第1の周波数範囲をカバーする帯域内サブ信号及び前記スペクトル信号の残りの周波数範囲をカバーする帯域外サブ信号に分割し、前記帯域内サブ信号及び前記帯域外サブ信号を別々にソートして、前記記述子を表すソートされた帯域内サブ信号及びソートされた帯域外サブ信号を得る、
    請求項1に記載の装置。
  5. 前記ソートユニットは、前記帯域内サブ信号が、前記スペクトル信号の周波数範囲の下部分又は前記スペクトル信号の最高周波数ピークの周りの周波数範囲の部分をカバーするように、前記スペクトル信号を分割する、
    請求項4に記載の装置。
  6. 前記変換ユニットが、前記スペクトル信号を正規化する、
    請求項1に記載の装置。
  7. 2以上の反復を実行するように前記変換ユニット及び前記ソートユニットを制御する制御ユニットであって、前記ソートユニットから出力された前記ソートされたスペクトル信号が、次の反復において前記変換ユニットに対する入力PPG信号として使用される、当該制御ユニット、
    を更に有する、請求項1に記載の装置。
  8. 前記分類器が、各反復において前記ソートユニットから出力された前記ソートされたスペクトル信号を連結し、前記関心領域を生物の皮膚領域として又は非皮膚領域として分類するのに前記連結されたソートされたスペクトル信号を記述子として使用する、
    請求項1に記載の装置。
  9. 前記抽出ユニットが、画像フレームごとに前記画像データの画素のグループの画像データ値を結合、特に平均し、前記結合された画像データ値から前記PPG信号を得る、
    請求項1に記載の装置。
  10. 前記抽出ユニットが、波長又は波長範囲において画像フレームごとに前記画像データの画素のグループの画像データ値を結合、特に平均し、前記結合された画像データ値から前記PPG信号を得る、
    請求項1に記載の装置。
  11. 前記抽出ユニットが、画素又は画素のグループごとに及び画像フレームごとに、少なくとも2つの異なる波長チャネルの画像データ値を重み付け平均として結合し、前記結合された画像データ値から前記PPG信号を得る、
    請求項1に記載の装置。
  12. 前記抽出ユニットが、正規化された血液量パルスベクトルシグネチャベースの方法、クロミナンスベースの方法、ブラインド音源分離法、主成分分析又は独立成分分析を使用して前記重みを計算する、
    請求項11に記載の装置。
  13. シーンの画像データを取得する撮像ユニットと、
    前記シーンの前記取得された画像データに基づく請求項1に記載の皮膚検出に対する方法と、
    を有する皮膚検出に対するシステム。
  14. シーンの画像データを得るステップであって、前記画像データが、時系列の画像フレームを有する、ステップと、
    前記画像データの関心領域からフォトプレチスモグラフィPPG信号を抽出するステップと、
    前記PPG信号をスペクトル信号に変換するステップと、
    前記スペクトル信号をソートして、記述子を表すソートされたスペクトル信号を得るステップと、
    前記記述子に基づいて前記関心領域を生物の皮膚領域として又は非皮膚領域として分類するステップと、
    を有する皮膚検出に対する方法。
  15. コンピュータにより実行される場合に、前記コンピュータに、請求項14に記載の方法のステップを実行させるプログラムコード手段を有するコンピュータプログラム。
JP2018540408A 2016-02-08 2017-02-08 皮膚検出に対する装置、システム及び方法 Active JP6538287B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16154665.0 2016-02-08
EP16154665 2016-02-08
PCT/EP2017/052748 WO2017137435A1 (en) 2016-02-08 2017-02-08 Device, system and method for skin detection

Publications (2)

Publication Number Publication Date
JP2019510532A true JP2019510532A (ja) 2019-04-18
JP6538287B2 JP6538287B2 (ja) 2019-07-03

Family

ID=55359412

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2018540408A Active JP6538287B2 (ja) 2016-02-08 2017-02-08 皮膚検出に対する装置、システム及び方法
JP2018540804A Pending JP2019508116A (ja) 2016-02-08 2017-02-08 拍動検出のための装置、システム、及び方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2018540804A Pending JP2019508116A (ja) 2016-02-08 2017-02-08 拍動検出のための装置、システム、及び方法

Country Status (6)

Country Link
US (2) US20190050985A1 (ja)
EP (2) EP3414739B1 (ja)
JP (2) JP6538287B2 (ja)
CN (2) CN108701357B (ja)
RU (1) RU2018132127A (ja)
WO (2) WO2017137415A1 (ja)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190050985A1 (en) * 2016-02-08 2019-02-14 Koninklijke Philips N.V. Device, system and method for pulsatility detection
US10667723B2 (en) 2016-02-19 2020-06-02 Covidien Lp Systems and methods for video-based monitoring of vital signs
KR20230116072A (ko) * 2017-02-24 2023-08-03 서니브룩 리서치 인스티튜트 이미징에서의 노이즈 감소를 위한 시스템 및 방법
WO2019094893A1 (en) 2017-11-13 2019-05-16 Covidien Lp Systems and methods for video-based monitoring of a patient
US11311202B2 (en) 2017-11-14 2022-04-26 Arizona Board Of Regents On Behalf Of Arizona State University Robust real-time heart rate monitoring method based on heartbeat harmonics using small-scale radar
US11712176B2 (en) 2018-01-08 2023-08-01 Covidien, LP Systems and methods for video-based non-contact tidal volume monitoring
WO2019240991A1 (en) 2018-06-15 2019-12-19 Covidien Lp Systems and methods for video-based patient monitoring during surgery
US11311252B2 (en) 2018-08-09 2022-04-26 Covidien Lp Video-based patient monitoring systems and associated methods for detecting and monitoring breathing
US11617520B2 (en) 2018-12-14 2023-04-04 Covidien Lp Depth sensing visualization modes for non-contact monitoring
US11564633B2 (en) * 2018-12-21 2023-01-31 Industrial Technology Research Institute State assessment system, diagnosis and treatment system, and method for operating the diagnosis and treatment system
CN109864740B (zh) * 2018-12-25 2022-02-01 北京津发科技股份有限公司 一种运动状态的表面肌电信号采集传感器和设备
US11315275B2 (en) 2019-01-28 2022-04-26 Covidien Lp Edge handling methods for associated depth sensing camera devices, systems, and methods
EP3698704A1 (en) * 2019-02-20 2020-08-26 Koninklijke Philips N.V. Device, system and method for determining physiological information
US11771380B2 (en) 2019-03-19 2023-10-03 Arizona Board Of Regents On Behalf Of Arizona State University Vital sign monitoring system using an optical sensor
US11783483B2 (en) 2019-03-19 2023-10-10 Arizona Board Of Regents On Behalf Of Arizona State University Detecting abnormalities in vital signs of subjects of videos
KR20200133112A (ko) * 2019-05-17 2020-11-26 삼성전자주식회사 생체 신호 측정 장치
CN112037162B (zh) * 2019-05-17 2022-08-02 荣耀终端有限公司 一种面部痤疮的检测方法及设备
CN112336318B (zh) * 2019-08-09 2022-02-18 复旦大学 一种自适应多模态融合的脉搏位置精确定位方法
CN110859615B (zh) * 2019-12-06 2020-07-31 电子科技大学 一种基于幅度排列的生理信号时间不可逆分析方法
CN111227795B (zh) * 2020-01-13 2022-07-22 天津众阳科技有限公司 一种人脸反射区颜色变化信息的提取方法
US11484208B2 (en) 2020-01-31 2022-11-01 Covidien Lp Attached sensor activation of additionally-streamed physiological parameters from non-contact monitoring systems and associated devices, systems, and methods
US11359915B1 (en) * 2020-05-15 2022-06-14 Lockheed Martin Corporation Active optical compressive sensing
CN112017155B (zh) * 2020-07-13 2023-12-26 浙江华锐捷技术有限公司 健康体征数据的测量方法、装置、系统和存储介质
CN112168171A (zh) * 2020-10-16 2021-01-05 北京中科心研科技有限公司 一种基于穿戴式脉搏信号的皮肤肤色检测方法及装置
CN112633093B (zh) * 2020-12-10 2024-04-16 中国人民解放军32181部队 一种装备状态检测方法及系统
US11769313B2 (en) 2021-05-21 2023-09-26 Ford Global Technologies, Llc Counterfeit image detection
US11636700B2 (en) 2021-05-21 2023-04-25 Ford Global Technologies, Llc Camera identification
US11967184B2 (en) 2021-05-21 2024-04-23 Ford Global Technologies, Llc Counterfeit image detection
EP4098183A1 (en) 2021-06-02 2022-12-07 Koninklijke Philips N.V. Blood pressure detection
SE545755C2 (en) * 2022-02-25 2024-01-02 Detectivio Ab Non-contact oxygen saturation estimation
CN115661131B (zh) * 2022-11-17 2023-03-10 菲特(天津)检测技术有限公司 一种图像识别方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130296660A1 (en) * 2012-05-02 2013-11-07 Georgia Health Sciences University Methods and systems for measuring dynamic changes in the physiological parameters of a subject
WO2014095759A1 (en) * 2012-12-21 2014-06-26 Koninklijke Philips N.V. System and method for extracting physiological information from remotely detected electromagnetic radiation
US20150302158A1 (en) * 2014-04-21 2015-10-22 Microsoft Corporation Video-based pulse measurement

Family Cites Families (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2165706A1 (ja) * 1971-11-23 1973-08-10 Commissariat Energie Atomique
US5309917A (en) * 1991-09-12 1994-05-10 Drexel University System and method of impedance cardiography and heartbeat determination
DE102004002918B4 (de) * 2004-01-20 2016-11-10 Siemens Healthcare Gmbh Vorrichtung zur Untersuchung der Haut
US7446316B2 (en) * 2005-10-31 2008-11-04 Honeywell International Inc. Skin detection sensor
US8504146B2 (en) * 2007-06-29 2013-08-06 The Regents Of The University Of California Multi-channel myoelectrical control using single muscle
CN102341828B (zh) * 2009-03-06 2014-03-12 皇家飞利浦电子股份有限公司 处理至少一个活体的图像
CN101716069A (zh) * 2009-11-17 2010-06-02 中国科学院安徽光学精密机械研究所 人体氧化应激无创荧光检测装置与方法
US8520074B2 (en) * 2010-12-14 2013-08-27 Xerox Corporation Determining a total number of people in an IR image obtained via an IR imaging system
JP5984841B2 (ja) 2011-01-05 2016-09-06 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 特性信号から情報を抽出する装置及び方法
WO2012128614A1 (en) * 2011-03-24 2012-09-27 Erasmus University Medical Center Rotterdam Method to determine the absorption coefficient in turbid media
US8838209B2 (en) * 2012-02-21 2014-09-16 Xerox Corporation Deriving arterial pulse transit time from a source video image
US10452894B2 (en) * 2012-06-26 2019-10-22 Qualcomm Incorporated Systems and method for facial verification
CN104853670B (zh) * 2012-11-23 2018-02-16 皇家飞利浦有限公司 用于提取生理信息的设备和方法
CN103149177B (zh) * 2013-01-14 2014-12-17 天津先阳科技发展有限公司 压力调制近红外光谱实现生物组织检测的装置和方法
EP2967376B1 (en) * 2013-03-14 2023-02-15 Koninklijke Philips N.V. Device and method for determining vital signs of a subject
CN104434078A (zh) * 2013-09-13 2015-03-25 施乐公司 用时间序列信号确定基于视频的脉搏传导时间的系统和方法
EP3052008B1 (en) 2013-10-01 2017-08-30 Koninklijke Philips N.V. Improved signal selection for obtaining a remote photoplethysmographic waveform
CN105636506B (zh) * 2013-10-17 2018-11-27 皇家飞利浦有限公司 用于远程光体积描记法的自动相机调节
US20150148632A1 (en) * 2013-11-26 2015-05-28 David Alan Benaron Calorie Monitoring Sensor And Method For Cell Phones, Smart Watches, Occupancy Sensors, And Wearables
JP6608822B2 (ja) * 2013-12-12 2019-11-20 コーニンクレッカ フィリップス エヌ ヴェ 対象のバイタルサインを決定する装置及び方法
EP3094248B1 (en) * 2014-01-16 2020-05-27 Nokia Technologies Oy Method and device for the detection of respiratory rate
US9336594B2 (en) * 2014-03-07 2016-05-10 Xerox Corporation Cardiac pulse rate estimation from source video data
US10028669B2 (en) * 2014-04-02 2018-07-24 Massachusetts Institute Of Technology Methods and apparatus for physiological measurement using color band photoplethysmographic sensor
US9459201B2 (en) * 2014-09-29 2016-10-04 Zyomed Corp. Systems and methods for noninvasive blood glucose and other analyte detection and measurement using collision computing
US10660533B2 (en) * 2014-09-30 2020-05-26 Rapsodo Pte. Ltd. Remote heart rate monitoring based on imaging for moving subjects
US9770213B2 (en) * 2014-10-30 2017-09-26 Koninklijke Philips N.V. Device, system and method for extracting physiological information
US9432133B2 (en) * 2015-02-03 2016-08-30 Keysight Technologies, Inc. Vector signal generator channel characterization employing orthogonal modulation
US10004427B1 (en) * 2015-08-04 2018-06-26 Verily Life Sciences Llc Methods, systems, and devices for determining a respiration rate
US20190050985A1 (en) * 2016-02-08 2019-02-14 Koninklijke Philips N.V. Device, system and method for pulsatility detection

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130296660A1 (en) * 2012-05-02 2013-11-07 Georgia Health Sciences University Methods and systems for measuring dynamic changes in the physiological parameters of a subject
WO2014095759A1 (en) * 2012-12-21 2014-06-26 Koninklijke Philips N.V. System and method for extracting physiological information from remotely detected electromagnetic radiation
US20150302158A1 (en) * 2014-04-21 2015-10-22 Microsoft Corporation Video-based pulse measurement

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
WENJIN WANG ET AL.: ""Exploiting Spatial Redundancy of Image Sensor for Motion Robust rPPG"", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, vol. 62, no. 2, JPN6018048964, 28 February 2015 (2015-02-28), pages 415 - 425, XP011570468, ISSN: 0003936938, DOI: 10.1109/TBME.2014.2356291 *

Also Published As

Publication number Publication date
JP6538287B2 (ja) 2019-07-03
CN108604376B (zh) 2022-04-01
RU2018132127A (ru) 2020-03-11
US20190057502A1 (en) 2019-02-21
CN108604376A (zh) 2018-09-28
EP3414738A1 (en) 2018-12-19
CN108701357A (zh) 2018-10-23
WO2017137415A1 (en) 2017-08-17
US10818010B2 (en) 2020-10-27
JP2019508116A (ja) 2019-03-28
EP3414738B1 (en) 2020-01-22
EP3414739B1 (en) 2020-04-08
EP3414739A1 (en) 2018-12-19
WO2017137435A1 (en) 2017-08-17
US20190050985A1 (en) 2019-02-14
CN108701357B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
JP6538287B2 (ja) 皮膚検出に対する装置、システム及び方法
EP3355787B1 (en) Device, system and method for extracting physiological information
JP6194105B2 (ja) 遠隔フォトプレチスモグラフィ波形を取得するための改良された信号選択
US10524725B2 (en) Device, system and method for detecting apnoea of a subject
CN106999116A (zh) 用于皮肤检测的设备和方法
WO2014140978A1 (en) Device and method for obtaining vital sign information of a subject
JP6970840B2 (ja) 対象を含むシーンの画像の画像セグメンテーションのための装置、システム及び方法
EP3422931B1 (en) Device, system and method for determining a vital sign of a subject
WO2019145142A1 (en) Device, system and method for determining at least one vital sign of a subject

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190306

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190605

R150 Certificate of patent or registration of utility model

Ref document number: 6538287

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