(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2024-09-05
(54)【発明の名称】生理学的システムを特徴付ける際に使用するために生物物理学的信号から視覚的特徴を工学設計するための方法及びシステム
(51)【国際特許分類】
A61B 5/02 20060101AFI20240829BHJP
A61B 5/33 20210101ALI20240829BHJP
【FI】
A61B5/02 310A
A61B5/33 200
【審査請求】未請求
【予備審査請求】未請求
(21)【出願番号】P 2024512055
(86)(22)【出願日】2022-08-19
(85)【翻訳文提出日】2024-04-17
(86)【国際出願番号】 IB2022057794
(87)【国際公開番号】W WO2023026151
(87)【国際公開日】2023-03-02
(32)【優先日】2021-08-23
(33)【優先権主張国・地域又は機関】US
(81)【指定国・地域】
(71)【出願人】
【識別番号】517099409
【氏名又は名称】アナリティクス フォー ライフ インコーポレイテッド
(74)【代理人】
【識別番号】100078282
【氏名又は名称】山本 秀策
(74)【代理人】
【識別番号】100113413
【氏名又は名称】森下 夏樹
(74)【代理人】
【識別番号】100181674
【氏名又は名称】飯田 貴敏
(74)【代理人】
【識別番号】100181641
【氏名又は名称】石川 大輔
(74)【代理人】
【識別番号】230113332
【氏名又は名称】山本 健策
(72)【発明者】
【氏名】ドゥームラ, アビナブ
【テーマコード(参考)】
4C017
4C127
【Fターム(参考)】
4C017AA09
4C017AA19
4C017AB03
4C017AC16
4C017AC28
4C017BC11
4C017BD06
4C127AA02
4C127GG02
4C127GG10
4C127GG16
(57)【要約】
好ましい実施形態では、患者が安静にしている間に患者に配置された表面センサから非侵襲的に取得される心臓/生体電位信号及び/又はフォトプレチスモグラフィ信号などの生物物理学的信号から決定された1つ以上の視覚的特徴又はパラメータの使用を容易にする臨床評価システム及び方法が開示される。視覚的特徴又はパラメータは、モデル又は分類器(例えば、機械学習分類器)において使用され、疾患の存在若しくは不在、医学的コンディション、又はいずれかの兆候を含む、患者の生理学的状態に関連付けられたメトリックを推定することができる。推定されたメトリックを使用して、医師又は他のヘルスケア提供者が、疾患若しくはコンディションの存在又は不在及び/又は重症度及び/又は局在化を診断すること、あるいは当該疾患又はコンディションの治療を支援することができる。
【選択図】
図1
【特許請求の範囲】
【請求項1】
対象の疾患状態又は異常コンディションに関連付けられたメトリックを非侵襲的に査定するための方法であって、前記方法は、
1つ以上のプロセッサによって、前記対象の生物物理学的信号データセットを得ることと、
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットの1つ以上の視覚予測器関連特性の値を決定することと、
前記1つ以上のプロセッサによって、前記1つ以上の視覚予測器関連特性の前記決定された値の、前記メトリックの推定モデルへの適用に部分的に基づいて、前記疾患状態又は異常コンディションに関連付けられたメトリックの存在についての推定値を決定することと、を含み、
前記メトリックの前記存在についての前記推定値は、予想される疾患状態若しくはコンディションの診断において使用するため、又は前記予想される疾患状態若しくはコンディションの治療を方向付けるために、前記予想される疾患状態又はコンディションの存在を非侵襲的に推定するために前記推定モデルにおいて使用される、方法。
【請求項2】
前記生物物理学的信号データセットは、3つの測定チャネルについて取得された生体電位信号を含む、請求項1に記載の方法。
【請求項3】
前記生物物理学的信号データセットは、光学センサから取得されたフォトプレチスモグラフィック信号を含む、請求項1に記載の方法。
【請求項4】
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットの位相空間モデルを生成することと、
前記1つ以上のプロセッサによって、前記位相空間モデルから抽出された1つ以上の特徴の1つ以上の値を決定することと、を更に含み、前記1つ以上の特徴は、
前記位相空間モデルの三次元外周に関連付けられた特徴と、
前記位相空間モデルにおいて決定された最大適合平面に囲まれた面積に関連付けられた特徴と、
前記位相空間モデルにおいて決定された三次元平均曲率又は三次元最大曲率に関連付けられた特徴と、
前記位相空間モデル内のループを定義する点の回転の巨視的尺度に関連付けられた特徴と、
前記位相空間モデルにおいて決定されたカールベクトルの大きさの平均に関連付けられた特徴と、
前記位相空間モデルにおいて定義されたループ内の連続する点の間のユークリッド距離に関連付けられた特徴と、からなる群から選択される、請求項2又は3に記載の方法。
【請求項5】
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットの位相空間モデルの正射影を生成することと、
前記1つ以上のプロセッサによって、前記正射影から抽出された1つ以上の特徴の1つ以上の値を決定することと、を更に含み、前記1つ以上の特徴は、
前記正射影において定義されたループの二次元外周に関連付けられた特徴と、
前記正射影において定義された前記ループの最大二次元外周を有する象限を定義する特徴と、
前記正射影において定義された前記ループの二次元面積に関連付けられた特徴と、
前記正射影において定義された前記ループの最大二次元面積を有する象限を定義する特徴と、
前記正射影において定義された前記ループの前記偏心に関連付けられた特徴と、
前記正射影において決定された二次元平均曲率又は二次元最大曲率に関連付けられた特徴と、
前記正射影における前記ループを定義する点の前記回転の巨視的尺度に関連付けられた特徴と、からなる群から選択される、請求項1~4のいずれか一項に記載の方法。
【請求項6】
前記1つ以上のプロセッサによって、フォトプレチスモグラフィック信号又はその派生物を含む前記生物物理学的信号データセットから抽出された特徴の1つ以上の値を決定することを更に含み、前記1つ以上の特徴は、
前記フォトプレチスモグラフィック信号内の事前定義された起点目印を前記フォトプレチスモグラフィック信号内で決定されたピーク位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記速度プレチスモグラフィック信号内で決定されたピーク位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記速度プレチスモグラフィック信号内で決定された最小位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記速度プレチスモグラフィック信号内で決定されたベース位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された加速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記加速度プレチスモグラフィック信号内で決定されたピーク位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記加速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記加速度プレチスモグラフィック信号内で決定された最小位置に結合するベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記加速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記加速度プレチスモグラフィック信号内で決定されたベース位置に結合するベクトルによって定義される特徴と、からなる群から選択される、請求項1~5のいずれか一項に記載の方法。
【請求項7】
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットの前記位相空間モデルを生成することと、
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットから抽出された特徴の1つ以上の値を決定することと、を更に含み、前記1つ以上の特徴は、
前記心臓信号内の前記心臓信号ピーク目印内の事前定義された起点目印を結合する三次元ベクトルによって定義される特徴と、
前記フォトプレチスモグラフィック信号内の事前定義された起点目印を、前記フォトプレチスモグラフィック信号内で決定されたピーク、最小値、又はベース位置に結合する三次元ベクトルによって定義される特徴と、
前記心臓信号内で前記事前定義された起点目印を前記ピーク目印に結合する前記三次元ベクトルを定義する仰角によって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記速度プレチスモグラフィック信号内で決定されたピーク、最小値、又はベース位置に結合する三次元ベクトルによって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された加速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記加速度プレチスモグラフィック信号内で決定されたピーク、最小値、又はベース位置に結合する三次元ベクトルによって定義される特徴と、
前記フォトプレチスモグラフィック信号内の前記事前定義された起点目印を前記フォトプレチスモグラフィック信号内で決定された前記ピーク、前記最小値、又は前記ベース位置に結合する前記三次元ベクトルを定義する仰角によって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記速度プレチスモグラフィック信号内で決定された前記ピーク、前記最小値、又は前記ベース位置に結合する前記三次元ベクトルを定義する仰角によって定義される特徴と、
i)前記フォトプレチスモグラフィック信号から導出された前記加速度プレチスモグラフィック信号内の前記事前定義された起点目印を、ii)前記加速度プレチスモグラフィック信号内で決定された前記ピーク、前記最小値、又は前記ベース位置に結合する前記三次元ベクトルを定義する仰角によって定義される特徴と、からなる群から選択される、請求項1~6のいずれか一項に記載の方法。
【請求項8】
前記1つ以上のプロセッサによって、前記位相空間モデルを生成することと、
前記1つ以上のプロセッサによって、フォトプレチスモグラフィック信号又はその派生物を含む前記生物物理学的信号データセットから抽出された特徴の1つ以上の値を決定することを更に含み、前記1つ以上の特徴は、
i)前記位相空間モデル内で定義される直交平面のうちの1つ以上の投影内の原点位置と、ii)前記直交平面内のピーク、最小値、又はベース位置との間で定義されたベクトルの二次元大きさによって定義される特徴と、
i)前記位相空間モデル内で定義された前記直交平面のうちの前記1つ以上の前記投影内の前記原点位置と、ii)前記直交平面内の前記ピーク、前記最小値、又は前記ベース位置との間で定義された前記ベクトルの角度によって定義される特徴と、からなる群から選択される、請求項1~7のいずれか一項に記載の方法。
【請求項9】
前記1つ以上のプロセッサによって、前記生物物理学的信号データセットの前記位相空間モデルを生成することと、
前記1つ以上のプロセッサによって、フォトプレチスモグラフィック信号又はその派生物を含む前記生物物理学的信号データセットから抽出された特徴の1つ以上の値を決定することと、を更に含み、前記1つ以上の特徴は、前記位相空間モデル内の定義された幾何学的形状の表面積又は体積パラメータによって定義される特徴を含む、請求項1~8のいずれか一項に記載の方法。
【請求項10】
前記定義された幾何学的形状は、アルファ包形状又は凸包形状を含む、請求項9に記載の方法。
【請求項11】
前記1つ以上のプロセッサによって、前記疾患状態又は異常コンディションの前記存在についての前記推定値の視覚化の生成を引き起こすことを更に含み、前記生成された視覚化は、コンピューティングデバイスのディスプレイにおいてレンダリングかつ表示され、及び/又は報告において提示される、請求項1~10のいずれか一項に記載の方法。
【請求項12】
1つ以上の視覚予測器関連特性の前記値は、線形モデル、決定木モデル、ランダムフォレストモデル、サポートベクトルマシンモデル、及びニューラルネットワークモデルからなる群から選択されるモデルにおいて使用される、請求項1~11のいずれか一項に記載の方法。
【請求項13】
前記モデルは、
1つ以上の脱分極又は再分極波伝播関連特徴と、
1つ以上の脱分極波伝播偏差関連特徴と、
1つ以上のサイクル変動性関連特徴と、
1つ以上の動的システム関連特徴と、
1つ以上の心臓波形トポロジック及び変動関連特徴と、
1つ以上のPPG波形トポロジック及び変動関連特徴と、
1つ以上の心臓又はPPG信号パワースペクトル密度関連特徴と、
1つ以上の心臓又はPPG信号視覚関連特徴と、
1つ以上の予測可能性特徴と、からなる群から選択される特徴を更に含む、請求項12に記載の方法。
【請求項14】
前記疾患状態又は異常コンディションは、冠動脈疾患、肺高血圧、肺動脈高血圧、左心疾患に起因する肺高血圧、肺高血圧につながる稀な障害、左心室心不全又は左側心不全、右心室心不全又は右側心不全、収縮期心不全、拡張期心不全、虚血性心疾患、及び不整脈からなる群から選択される、請求項1~13のいずれか一項に記載の方法。
【請求項15】
測定システムの1つ以上の取得回路によって、前記1つ以上のチャネルにわたって電圧勾配信号を取得することであって、前記電圧勾配信号は、約1kHzを超える周波数において取得される、取得することと、
前記1つ以上の取得回路によって、前記取得された電圧勾配信号から前記得られた生物物理学的データセットを生成することと、を更に含む、請求項1~14のいずれか一項に記載の方法。
【請求項16】
測定システムの1つ以上の取得回路によって、1つ以上のフォトプレチスモグラフィック信号を取得することと、
前記1つ以上の取得回路によって、前記取得された電圧勾配信号から前記得られた生物物理学的データセットを生成することと、を更に含む、請求項1~14のいずれか一項に記載の方法。
【請求項17】
前記1つ以上のプロセッサは、クラウドプラットフォーム内に位置する、請求項1~16のいずれか一項に記載の方法。
【請求項18】
前記1つ以上のプロセッサは、ローカルコンピューティングデバイス内に位置する、請求項1~16のいずれか一項に記載の方法。
【請求項19】
システムであって、
プロセッサと、
命令が記憶されたメモリと、を備え、前記プロセッサによる前記命令の実行は、前記プロセッサに、請求項1~18に記載の方法のいずれか1つを実施させる、システム。
【請求項20】
命令を記憶した非一時的コンピュータ可読媒体であって、前記プロセッサによる前記命令の実行は、前記プロセッサに、請求項1~18に記載の方法のいずれか1つを実施させる、非一時的コンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【0001】
関連出願
本PCT出願は、2021年8月23日に出願された「Methods and Systems for Engineering Visual Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/236,072号の優先権及び利益を主張し、その全体が参照により本明細書に組み込まれる。
【0002】
本開示は、概して、診断用途で使用するために、生物物理学的信号から特徴又はパラメータを工学設計するための方法及びシステムに関し、具体的には、1つ以上の生理学的システム並びにそれらの関連付けられた機能、活動、及び異常を特徴付ける際に使用するための視覚的特徴の工学設計及び使用に関する。特徴又はパラメータはまた、監視若しくは追跡、医療機器の制御のために、又は疾患、医学的コンディション、若しくはいずれかの兆候の治療を誘導するために使用することができる。
【背景技術】
【0003】
ヘルスケア専門家が疾患を診断するのを支援するための多くの方法及びシステムがある。これらのいくつかは、侵襲的若しくは最小侵襲的な技術、放射線、運動若しくはストレス、又は薬理学的薬剤の使用を伴い、時にはそれらの付随するリスク及び他の不利益と組み合わせて使用する。
【0004】
拡張期心不全は、罹患率及び死亡率の主な原因であり、左心室機能が維持されている患者における心不全の症状として定義される。これは、コンプライアンスが低下し、弛緩が損なわれ、左心臓カテーテル法によって測定される左心室の拡張末期圧が増加する、硬い左心室を特徴とする。肺高血圧(PH)及び肺動脈高血圧(PAH)の診断のための現在の臨床標準的なケアは、特に、肺動脈内の圧力を直接測定する心臓の右側の心臓カテーテルを含む。冠動脈造影は、治療医によって説明される冠動脈病変を介して決定される、冠動脈疾患(CAD)を査定するために使用される現在の標準的なケアである。磁気共鳴イメージング及びコンピュータ断層撮影などの非侵襲的イメージングシステムは、放射線科医によってレビューされる患者の血流及び動脈閉塞の画像を取得するために専門的な施設を必要とする。
【0005】
前述の欠点なしに、心臓疾患及び様々な他の疾患及びコンディションの診断においてヘルスケア専門家を支援することができるシステムを有することが望ましい。
【発明の概要】
【課題を解決するための手段】
【0006】
好ましい実施形態では、患者が安静にしている間に患者に配置された表面センサから非侵襲的に取得される心臓/生体電位信号及び/又はフォトプレチスモグラフィ信号などの生物物理学的信号から決定された1つ以上の視覚的特徴又はパラメータの使用を容易にする臨床評価システム及び方法が開示される。視覚的特徴又はパラメータは、モデル又は分類器(例えば、機械学習分類器)において使用され、疾患の存在若しくは不在、医学的コンディション、又はいずれかの兆候を含む、患者の生理学的状態に関連付けられたメトリックを推定することができる。推定されたメトリックを使用して、医師又は他のヘルスケア提供者が、疾患若しくはコンディションの存在又は不在及び/又は重症度及び/又は局在化を診断すること、あるいは当該疾患又はコンディションの治療を支援することができる。
【0007】
疾患、コンディション、又はいずれかの兆候の存在又は不在の推定又は決定された可能性は、疾患又は医学的コンディションの査定のための他の評価又は測定モダリティを置き換える、増強する、又は置き換えることができる。場合によっては、決定は、数値スコア及び関連情報の形態をとることができる。
【0008】
心臓信号の場合、心房脱分極(AD)に関連付けられた波形、心房再分極(AR)に関連付けられた波形、心室脱分極(VD)に関連付けられた波形、及び心室再分極(VR)に関連付けられた波形の特性を含む、心臓電気生理学的信号の反復サイクルの様々な特性を特徴付けるために、心臓又は生体電位信号の位相空間モデルにおいて視覚的特徴が決定される。いくつかの実施形態では、視覚的特徴は、3D位相空間内及び3D空間のオクタント内の心臓ループ(AD、VD、VR)の特性を含む。視覚的特徴と称されるが、ARはVD中に発生し、位相空間ビューでは見えないことがある。いくつかの実施形態では、視覚的特徴は、3D位相空間内及び3D空間のオクタント(又はサブ領域)内の心臓ベクトル(例えば、最大心房脱分極ベクトル(MADV)、心房再分極ベクトル(ARV)、最大心室脱分極ベクトル(MVDV)、初期心室脱分極ベクトル(IVDV)、終末心室脱分極ベクトル(TVDV)、最大心室再分極ベクトル(MVRV))の特性を含む。いくつかの実施形態では、視覚的特徴は、3つの直交面上及びそれぞれの直交面の象限内への心臓ループ(AD、VD、VR)の投影の特性を含む。いくつかの実施形態では、視覚的特徴は、3つの直交平面上への、及びそれぞれの直交平面の象限内の心臓ベクトルの投影の特性を含む。
【0009】
フォトプレチスモグラフィック関連信号に関して、視覚的特徴は、フォトプレチスモグラフィック信号と、フォトプレチスモグラフィック信号から導出された速度プレチスモグラフィック及び加速度プレチスモグラフィック信号との位相空間モデルにおいて決定され、これらの信号の反復サイクルの様々な特性を特徴付ける。「フォトプレチスモグラフィック関連」という用語は、フォトプレチスモグラフィ及び関連する速度プレチスモグラフィック及び加速度プレチスモグラフィック信号を指す。
【0010】
本明細書で使用されるように、用語「特徴」(機械学習及びパターン認識の文脈において、並びに本明細書で使用されるように)は、概して、観察されている現象の個々の計測可能な特性又は特徴を指す。特徴は、分析によって定義され、共通のモデル又は分析フレームワークからの他の特徴と組み合わせてグループで決定され得る。
【0011】
本明細書で使用される場合、「メトリック」は、生理学的システム又はシステム内のいずれかの1つ以上の疾患、コンディション、又は兆候(複数可)の存在、不在、重症度、及び/又は局在化(該当する場合)の推定又は可能性を指す。注目すべきことに、例示される方法及びシステムは、本明細書に記載されるある特定の実施形態において、生物物理学的信号を取得する、及び/又はそうでなければ患者からデータを収集するために、並びに1つ以上のメトリックを介して他の評価モダリティを置き換える、増強する、又は置き換えることができる1つの疾患、コンディション、又はインジケーターについて評価するために、信号処理及び分類器動作におけるそれらの信号及び/又はデータを評価するために使用され得る。場合によっては、メトリックは、数値スコア及び関連情報の形態をとることができる。
【0012】
心臓血管系及び呼吸器系の文脈では、そのようなメトリックが関連することができる疾患及びコンディションの例としては、例えば、(i)心不全(例えば、左側又は右側の心不全;駆出率(HFpEF)が保持された心不全)、(ii)冠動脈疾患(CAD)、(iii)肺動脈高血圧(PAH)を含むがこれらに限定されない様々な形態の肺高血圧(PH)、(iv)左心室駆出率(LVEF)の異常、及び様々な他の疾患又はコンディションが挙げられる。ある特定の形態の心不全の例示的なインジケーターは、上昇した又は異常な左心室拡張末期圧(LVEDP)の存在又は不在である。肺高血圧のある特定の形態の例示的なインジケーターは、上昇した又は異常な平均肺動脈圧(mPAP)の存在又は不在である。
【0013】
本明細書の一部に組み込まれ、それを構成する添付図面は、実施形態を例解し、本明細書と一緒に、方法及びシステムの原理を説明する役割を果たす。
【0014】
本発明の実施形態は、添付の図面と併せて読むと、以下の詳細な説明からよりよく理解され得る。例解のみを目的とするかかる実施形態は、本発明の新規かつ非自明の態様を描画する。図面は、以下の図を含む。
【図面の簡単な説明】
【0015】
【
図1】例示的な実施形態による、視覚的特徴又はパラメータを非侵襲的に計算して、患者の生理学的状態に関連付けられた1つ以上のメトリックを生成するように構成された例示的なモジュール又は構成要素の概略図である。
【
図2】例示的な実施形態による、例示的な生物物理学的信号キャプチャシステム又は構成要素、及び臨床設定における患者の生物物理学的信号を非侵襲的に収集することにおけるその使用を示す。
【
図3A】各々、診断、治療、監視、又は追跡のための実際の用途において、視覚的特徴/パラメータ又はそれらの中間データを使用するための例示的な方法を示す。
【
図3B】各々、診断、治療、監視、又は追跡のための実際の用途において、視覚的特徴/パラメータ又はそれらの中間データを使用するための例示的な方法を示す。
【
図4】例示的な実施形態による、視覚的特徴計算モジュールの一例を例解する。
【
図5】例示的な実施形態による、視覚的特徴計算モジュールの一例を例解する。
【
図6】例示的な実施形態による、
図4又は
図5の視覚的ループ又はベクトル特徴計算モジュールの例示的な動作方法を示す。
【
図7A】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図7B】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図7C】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図7D】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図7E】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図8A】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図8B】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図9A】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図9B】各々、例示的な実施形態による、
図4の視覚的特徴ループ計算モジュールによって実行される例示的な動作を示す。
【
図10A】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10B】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10C】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10D】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10E】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10F】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図10G】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図11A】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図11B】各々、例示的な実施形態による、
図5の視覚ベクトル計算モジュールによって実行される例示的な動作を示す。
【
図12A】例示的な実施形態による、患者の生理的状態に関連付けられた1つ以上のメトリックを生成するために、他の計算された特徴の中でも視覚的特徴を使用するように構成された例示的な臨床評価システムの概略図を示す。
【
図12B】例示的な実施形態による、
図12Aの例示的な臨床評価システムの動作の概略図を示す。
【発明を実施するための形態】
【0016】
本明細書に記載のありとあらゆる特徴、並びにそのような特徴のうちの2つ以上のありとあらゆる組み合わせは、そのような組み合わせに含まれる特徴が相互に矛盾しないことを条件として、本発明の範囲内に含まれる。
【0017】
本開示は、心臓関連病態及びコンディションの診断、追跡、及び治療における生物物理学的信号、例えば、生又は事前処理されたフォトプレチスモグラフィック信号、生体電位/心臓信号などの実用的な査定を対象としているが、このような査定は、生物物理学的信号が生体の任意の関連システムに関わる任意の病態又はコンディションの診断、追跡、及び治療(外科的、最小侵襲的、生活様式、栄養学的、及び/又は薬理学的治療などを含むがこれに限定されない)に適用され得る。査定は、医療機器若しくはウェアラブルデバイスの制御において、又はモニタリングアプリケーションにおいて使用され得る。
【0018】
本明細書で使用される「対象」及び「患者」という用語は、一般的に、例示的なシステム及び方法によって実施される分析を受ける人々を指すために互換的に使用される。
【0019】
本明細書で使用される「心臓信号」という用語は、例えば、心筋の収縮を引き起こす、その信号の電気/電気化学的伝導の態様を含む、心血管系の構造、機能、及び/又は活動に直接的若しくは間接的に関連する1つ以上の信号を指す。心臓信号は、いくつかの実施形態では、生体電位信号若しくは心電図信号、例えば、心電図(ECG)、心臓及びフォトプレチスモグラフィック波形、又は本明細書で後に説明される信号キャプチャ若しくは記録器具、又は他のモダリティを介して取得されるものを含み得る。
【0020】
本明細書で使用される「生物物理学的信号」という用語は、これらに限定されないが、1つ以上の心臓信号、神経信号、バリストカルジオグラフィック信号、及び/又はフォトプレチスモグラフィック信号を含むが、それはまた、情報が得られ得る任意の生理学的信号をより広く包含する。例によって限定されることを意図しないが、生物物理学的信号を、例えば、電気(例えば、時間及び/又は周波数などの様々なドメインにおける電圧/電位(例えば、生体電位)、インピーダンス、抵抗率、伝導性、電流などの測定などの技術によって観察、識別、及び/又は定量化され得る、ある特定の心臓及び神経系関連信号)、磁気、電磁気、光学(例えば、反射率、干渉法、分光法、吸光度、透過率、視覚的観察、フォトプレチスモグラフィ、及び同等物などの技術によって観察、識別、及び/又は定量化され得る信号)、音響、化学、機械(例えば、流体流動、圧力、運動、振動、変位、歪みに関連する信号)、熱、及び電気化学(例えば、グルコースなどのある特定の分析物の存在に相関され得る信号)を含むことができる、タイプ又はカテゴリに分類してもよい。場合によっては、生物物理学的信号は、生理学的系(例えば、呼吸器系、循環器系(循環器系、肺系)、神経系、リンパ系、内分泌系、消化器系、排泄物、筋肉系、骨格系、腎臓系/尿路系/排泄物系、免疫系、外皮系/外分泌系、及び生殖系)、1つ以上の臓器系(例えば、心臓及び肺が一緒に働くときに固有のものであり得る信号)の文脈において、又は組織(例えば、筋肉、脂肪、神経、結合組織、骨)、細胞、器官、分子(例えば、水、タンパク質、脂肪、炭水化物、ガス、遊離ラジカル、有機イオン、鉱物、酸、及び他の化合物)、元素、並びにそれらの亜原子成分の文脈において記載され得る。特に明記しない限り、「生物物理学的信号取得」という用語は、一般に、哺乳類又は非哺乳類生物などの生理学的系から生物物理学的信号を取得する任意の受動的又は能動的手段を指す。受動的及び能動的生物物理学的信号取得は、一般に、身体組織の自然又は誘導電気、磁気、光学、及び/又は音響放射線の観察を指す。受動的及び能動的生物物理学的信号取得手段の非限定的な例としては、例えば、体組織の自然放射線を観察する電圧/電位、電流、磁気、光学、音響、及び他の非能動的な方式が挙げられ、いくつかの例では、そのような放射線を誘導する。受動的及び能動的生物物理学的信号取得手段の非限定的な例としては、例えば、超音波、電波、マイクロ波、赤外線及び/又は可視光(例えば、パルスオキシメトリー又はフォトプレチスモグラフィでの使用のための)、可視光、紫外線、及び電離エネルギー又は放射線(例えば、X線)を含まない身体組織を能動的に調べる他の方式が挙げられる。能動的生物物理学的信号取得は、励起放出分光法(例えば、励起放出蛍光を含む)を含み得る。能動的生物物理学的信号取得はまた、イオン化エネルギー又は放射線(例えば、X線)(「イオン化生物物理学的信号」とも称される)を体組織に伝達することを含み得る。受動的及び能動的生物物理学的信号取得手段は、侵襲的手順(例えば、手術又は侵襲的放射線介入プロトコルを介して)又は非侵襲的(例えば、撮像、アブレーション、心臓収縮調節(例えば、ペースメーカーを介して)、カテーテル留置などを介して)と併せて実行され得る。
【0021】
本明細書で使用される「フォトプレチスモグラフィック信号」という用語は、赤色及び赤外線スペクトル内の波長を有する光などの酸素化及び脱酸素化ヘモグロビンによる光吸収の測定された変化に対応する、光学センサから取得された1つ以上の信号又は波形を指す。フォトプレチスモグラフィック信号(複数可)は、いくつかの実施形態において、パルスオキシメータ又はフォトプレスモグラム(PPG)を介して取得された生の信号(複数可)を含む。いくつかの実施形態において、フォトプレチスモグラフィック信号(複数可)は、健康を監視すること、及び/又は疾患若しくは異常コンディションを診断することを目的として、そのような信号波形を取得するように構成される、市販の、カスタムの、及び/又は専用の機器若しくは回路から取得される。フォトプレチスモグラフィック信号(複数可)は、典型的には、赤色のフォトプレチスモグラフィック信号(例えば、最も優勢に約625~740ナノメートルの波長を有する可視光スペクトル内の電磁信号)及び赤外線のフォトプレチスモグラフィック信号(例えば、可視スペクトルの公称赤縁から最大約1mmまで延在する電磁信号)を含むが、近赤外線、青、及び緑などの他のスペクトルは、用いられているPPGのタイプ及び/又はモードに応じて、異なる組み合わせで使用されてもよい。
【0022】
本明細書で使用される場合、「バリストカルジオグラフィック信号(ballistocardiographic signal)」という用語は、振動、音響、移動、又は向きを通して観察され得る全身を通る血流を概して反映する信号又は信号の群を指す。いくつかの実施形態では、バリストカルジオグラフィック信号は、心臓に近接して搭載されるセンサによって記録されるような身体の振動又は配向を測定することができる、振動、音響、移動、又は配向ベースのサイズモカルジオグラム(SCG)センサなどのウェアラブルデバイスによって取得される。サイズモカルジオグラムセンサは、概して、本明細書において「バリストカルジオグラム」という用語と互換的に使用される「サイズモカルジオグラム」を取得するために使用される。他の実施形態において、バリストカルジオグラフィック信号は、血液が頭部と足との間の長手方向に前後に移動するときの体重の変化などの現象を測定する外部機器、例えば、ベッド又は表面ベースの機器によって取得され得る。そのような実施形態において、各場所における血液量は、動力学的に変化し得、ベッド上の各場所で測定された重量並びにその重量の変化速度に反映され得る。
【0023】
加えて、本明細書の様々な実施形態に記載される方法及びシステムは、それほど限定されず、生体の別の生理学的システム、又はシステム、臓器、組織、細胞などの任意の文脈で利用され得る。例としてのみ、心血管の文脈において有用であり得る2つの生物物理学的信号タイプは、従来の心電図(ECG/EKG)機器を介して取得され得る心臓/生体電位信号、本明細書に記載されるものなどの他の機器から取得され得るバイポーラ広帯域生体電位(心臓)信号、及び、例えば、フォトプレチスモグラフィなどの様々な心電図法によって取得され得る信号を含む。別の例では、2つの生物物理学的信号タイプは、バリストカルジオグラフィック技術によって更に増強され得る。
【0024】
図1は、例示的な実施形態による、視覚的特徴又はパラメータを非侵襲的に計算して、分類器(例えば、機械学習分類器)を介して、患者の生理学的状態に関連付けられた1つ以上のメトリックを生成するように構成された例示的なモジュール又は構成要素の概略図である。モジュール又は構成要素は、製造用途、又は視覚的特徴及び他のクラスの特徴の開発に使用され得る。
【0025】
本明細書に記載の例示的な分析及び分類器は、心臓及び心臓肺関連の病態及び医学的コンディション、又はそのインジケーターの診断及び/又は治療においてヘルスケア提供者を支援するために使用され得る。例としては、本明細書に開示される様々な他の疾患及びコンディションの中で、有意な冠動脈疾患(CAD)、例えば、駆出率が保持された心不全(HFpEF)、うっ血性心不全、様々な形態の不整脈、弁不全、様々な形態の肺高血圧などの1つ以上の形態の心不全が挙げられる。
【0026】
加えて、いくつかの形態の心不全に関連する上昇した又は異常な左心室拡張末期圧(LVEDP)値、いくつかの形態の心不全に関連する異常な左心室駆出率(LVEF)値、又は肺高血圧及び/又は肺動脈高血圧に関連する上昇した平均肺動脈圧(mPAP)値などの疾患又はコンディションの可能性のインジケーターが存在する。本明細書に記載の例示的な分析及び分類器によって提供されるものなど、そのようなインジケーターが異常/上昇又は正常である可能性のインジケーターは、ヘルスケア提供者が、患者が所与の疾患又はコンディションを有するか、又は有しないかを査定又は診断するのを助けることができる。コンディションの疾患状態に関連付けられたこれらのメトリックに加えて、身体検査及び/又は他のテストの結果、患者の病歴、現在の薬物などの他の測定値及び要因の診断を行う際にヘルスケア専門家によって採用することができる。疾患状態又は医学的コンディションの存在又は不在の決定は、そのような疾患の兆候(又は診断に使用される測定のメトリック)を含むことができる。
【0027】
図1では、構成要素は、少なくとも1つの非侵襲的な生物物理学的信号レコーダ又はキャプチャシステム102、及び、例えばクラウド若しくはリモートインフラストラクチャ、又はローカルシステムに位置された査定システム103を含む。生物物理学的信号キャプチャシステム102(生物物理学的信号レコーダシステムとも称される)は、この実施形態では、例えば、同期的に取得された患者の電気信号及び血行動態信号を1つ以上のタイプの生物物理学的信号104として取得し、処理し、記憶し、送信するように構成される。
図1の例では、生物物理学的信号キャプチャシステム102は、測定プローブ106(例えば、血行動態信号104aのための血行動態センサを備える、例えば、プローブ106a及び106b、並びに電気/心臓信号104bのための導線を備えるプローブ106c~106hとして示される)から取得される、第1の生物物理学的信号104a(例えば、他の第1の生物物理学的信号に同期して取得される)及び第2の生物物理学的信号104b(例えば、他の生物物理学的信号に同期して取得される)として示される2つのタイプの生物物理学的信号を同期してキャプチャするように構成される。プローブ106a~hは、例えば、患者108の表面組織(患者場所108a及び108bに示される)に接着される、又はそれに隣接して配置されることによって、その上に配置される。患者は、好ましくは、ヒト患者であるが、任意の哺乳類患者であり得る。取得された生の生物物理学的信号(例えば、106a及び106b)は、生物物理学的信号データセット110(それぞれ、第1の生物物理学的信号データセット110a及び第2の生物物理学的信号データセット110bとして
図1に示される)を一緒に形成し、これは、例えば、好ましくは、記録/信号キャプチャ番号によって、及び/又は患者の名前及び医療記録番号によって識別可能な単一のファイルとして記憶され得る。
【0028】
図1の実施形態では、第1の生物物理学的信号データセット110aは、場所108aでの患者からの酸素化ヘモグロビン及び/又は脱酸素化ヘモグロビンの光吸収の測定された変化に関連付けられた、生のフォトプレチスモグラフィック信号又は血行動態信号(複数可)のセットを含み、第2の生物物理学的信号データセット110bは、心臓の電気信号に関連する生の心臓信号又は生体電位信号(複数可)のセットを含む。
図1において、生のフォトプレチスモグラフィック信号又は血行動態信号(複数可)は、患者の指で取得されるように示されているが、信号は、代替的に、患者のつま先、手首、額、耳たぶ、首などで取得されてもよい。同様に、心臓又は生体電位信号(複数可)は、直交誘導の3つのセットを介して取得されるように示されているが、他の誘導構成(例えば、11誘導構成、12誘導構成など)が使用されてもよい。
【0029】
プロット110a’及び110b’は、それぞれ、第1の生物物理学的信号データセット110a及び第2の生物物理学的信号データセット110aの例を示す。具体的には、プロット110a’は、取得されたフォトプレスモグラフィック又は血行動態信号の例を示す。プロット110a’において、フォトプレチスモグラフィック信号は、2つの光源(例えば、赤外線及び赤光源)から取得されたときの時間の関数としての信号電位を有する時系列信号である。プロット110b’は、3チャネル電位時系列プロットを含む例示的な心臓信号を示す。いくつかの実施形態において、生物物理学的信号キャプチャシステム102は、好ましくは、非侵襲的手段又は構成要素(複数可)を介して生物物理学的信号を取得する。代替の実施形態において、侵襲的又は最小侵襲的手段又は構成要素(複数可)は、非侵襲的手段(例えば、埋め込まれた圧力センサ、化学センサ、加速度計など)を補完するために、又はその代わりとして使用され得る。なおも更に代替の実施形態では、生物物理学的信号を収集することができる非侵襲的及び非接触プローブ又はセンサは、任意の組み合わせ(例えば、本明細書で考察されるように、受動温度計、スキャナ、カメラ、X線、磁気、又は非接触若しくは接触エネルギーデータ収集システムの他の手段)で、非侵襲的及び/又は侵襲的/最小侵襲的手段を補完するために、又はその代わりとして使用され得る。信号の取得及び記録に続いて、生物物理学的信号キャプチャシステム102はその後、例えば、無線又は有線通信システム及び/又はネットワークを送信することによって、取得された生物物理学的信号データセット110(又はそこから導出又は処理されたデータセット、例えば、フィルタリング又は事前処理されたデータ)を査定システム103のデータリポジトリ112(例えば、クラウドベースのストレージエリアネットワーク)に提供する。いくつかの実施形態において、取得された生物物理学的信号データセット110は、分析のために査定システム103に直接送信されるか、又は安全な臨床医のポータルを介してデータリポジトリ112にアップロードされる。
【0030】
生物物理学的信号キャプチャシステム102は、いくつかの実施形態において、回路及びコンピューティングハードウェア、ソフトウェア、ファームウェア、ミドルウェアなどで構成され、キャプチャされた生物物理学的信号の両方を取得、記憶、送信、及び任意選択的に処理して生物物理学的信号データセット110を生成する。例示的な生物物理学的信号キャプチャシステム102及び取得された生物物理学的信号セットデータ110は、「Method and Apparatus for Wide-Band Phase Gradient Signal Acquisition」と題された米国特許第10,542,898号、又は「Method and Apparatus for Wide-Band Phase Gradient Signal Acquisition」と題された米国特許公開第2018/0249960号に記載されており、これらの各々は、参照によりその全体が本明細書に組み込まれる。
【0031】
いくつかの実施形態において、生物物理学的信号キャプチャシステム102は、第1の生物物理学的信号(例えば、フォトプレチスモグラフィック信号)を取得するための第1の信号取得構成要素(図示せず)を含む2つ以上の信号取得構成要素を含み、第2の生物物理学的信号(例えば、心臓信号)を取得するための第2の信号取得構成要素(図示せず)を含む。いくつかの実施形態において、電気信号は、数分間でマルチキロヘルツの速度、例えば、1kHz~10kHzで取得される。他の実施形態において、電気信号は、10kHz~100kHzで取得される。血行動態信号は、例えば、100Hz~1kHzで取得されてもよい。
【0032】
生物物理学的信号キャプチャシステム102は、信号を取得するための1つ以上の他の信号取得構成要素(例えば、機械音響、バリストグラフィック、バリストカルジオグラフィックなどのセンサ)を含み得る。信号キャプチャシステム102の他の実施形態では、信号取得構成要素は、従来の心電図(ECG/EKG)機器(例えば、ホルターデバイス、12誘導ECGなど)を含む。
【0033】
査定システム103は、いくつかの実施形態では、データリポジトリ112と、分析エンジン又は分析器(図示せず、
図12A及び12B参照)とを備える。査定システム103は、特徴モジュール114及び分類器モジュール116(例えば、ML分類器モジュール)を含み得る。
図1において、査定システム103は、取得された生物物理学的信号データセット110を、例えば、データリポジトリ112から取り出し、それを特徴モジュール114において使用するように構成され、これは、視覚的特徴モジュール120及び他のモジュール122(本明細書において後に説明される)を含むように
図1に示される。特徴モジュール114は、分類器モジュール116に提供するために、視覚的特徴の値を含む特徴又はパラメータの値を計算し、分類器モジュールは、患者の生理学的状態(例えば、疾患状態、医学的コンディションの存在又は不在の兆候、あるいはいずれかの兆候)に関連付けられたメトリックの出力118、例えば、出力スコアを計算する。いくつかの実施形態において、出力118は、その後ヘルスケア医師ポータル(図示せず、
図12A及び
図12B参照)に提示されて、病態又は医学的コンディションの診断及び治療のためにヘルスケア専門家によって使用される。いくつかの実施形態において、ポータルは、例えば、患者、介護者、研究者などによるアクセスのために構成されてもよく(例えば、カスタマイズされてもよい)、出力118は、ポータルの意図された視聴者のために構成されてもよい。他のデータ及び情報はまた、出力118(例えば、取得された生物物理学的信号又は他の患者の情報及び病歴)の一部であり得る。
【0034】
分類器モジュール116(例えば、ML分類器モジュール)は、限定はしないが、決定木、ランダムフォレスト、ニューラルネットワーク、線形モデル、ガウス過程、最近傍、SVM、Naive Bayesなどのアルゴリズムに基づいて開発された伝達関数、ルックアップテーブル、モデル、又は演算子を含み得る。いくつかの実施形態では、分類器モジュール116は、2021年8月23日に出願された「Method and System to Non-Invasively Assess Elevated Left Ventricular End-Diastolic Pressure」と題する米国仮特許出願第63/235,960号、「Discovering Novel Features to Use in Machine Learning Techniques,such as Machine Learning Techniques for Diagnosing Medical Conditions」と題する米国特許出願公開第2019/0026430号又は「Discovering Genomes to Use in Machine Learning Techniques」と題する米国特許出願公開第2019/0026431号に記載されているML技術に基づいて開発されたモデルを含むことができ、これらの各々は、参照によりその全体が本明細書に組み込まれる。
【0035】
生物物理学的信号取得の例。
図2は、例示的な実施形態による、生物物理学的信号キャプチャシステム102(102aとして示される)と、臨床設定における患者の生物物理学的信号を非侵襲的に収集することにおけるその使用とを示す。
図2において、生物物理学的信号キャプチャシステム102aは、患者が安静している間に、患者108からの2つのタイプの生物物理学的信号をキャプチャするように構成される。生物物理学的信号キャプチャシステム102aは、患者の、(i)直交して配置されたセンサ(106c~106h、106iは7番目のコモンモード基準リードである)を使用して、胴体からの電気信号(例えば、第2の生物物理学的信号データセット110bに対応する心臓信号)、及び(ii)フォトプレチスモグラフィックセンサ(例えば、信号106a、106bを収集する)を使用して、指からの血行動態信号(例えば、第1の生物物理学的信号データセット110aに対応するPPG信号)を同期的に取得する。
【0036】
図2に示すように、電気信号及び血行動態信号(例えば、104a、104b)は、患者の皮膚に適用された市販のセンサを介して受動的に収集される。信号は、患者が電離放射線又は放射線造影剤に曝露することなく、かつ患者の運動又は薬理学的ストレス要因を使用することなく、有益に取得することができる。生物物理学的信号キャプチャシステム102aは、技術者又は看護師などのヘルスケア専門家が必要なデータを取得することにつながる、及びセルラー信号又はWi-Fi接続を確立することができる、任意の設定で使用され得る。
【0037】
電気信号(例えば、第2の生物物理学的信号データセット110bに対応する)は、基準リードと共に、患者の胸部及び背部にわたって配置された3つの直交ペアの表面電極を使用して収集される。いくつかの実施形態において、電気信号は、数分間(例えば、215秒間)、マルチキロヘルツ速度(例えば、6つのチャネルの各々について、1秒当たり8000サンプル)でローパスアンチエイリアシングフィルタ(例えば、約2kHz)を使用して取得される。代替の実施形態において、生物物理学的信号は、監視のために連続的/間欠的に取得されてもよく、取得された信号の部分が分析のために使用される。血行動態信号(例えば、第1の生物物理学的信号データセット110aに対応する)は、指に配置されたフォトプレチスモグラフィックセンサを使用して収集される。赤色光(例えば、600~750nmの任意の波長)及び赤外線光(例えば、850~950nmの任意の波長)の光吸収は、いくつかの実施形態において、同じ期間にわたって1秒当たり500サンプルの速度で記録される。生物物理学的信号キャプチャシステム102aは、信号中のコモンモード環境ノイズを低減するコモンモード駆動を含み得る。フォトプレチスモグラフィック信号及び心臓信号は、各患者について同時に取得された。データにおけるジッタ(モダリティ間ジッタ)は、約10マイクロ秒(μs)未満であり得る。心臓信号チャネル間のジッタは、10マイクロ秒未満、例えば、約10フェムト秒(fs)であり得る。
【0038】
患者メタデータ及び信号データを含む信号データパッケージは、信号取得手順の完了時にコンパイルされ得る。このデータパッケージは、生物物理学的信号キャプチャシステム102aがパッケージをデータリポジトリ112に転送する前に暗号化され得る。いくつかの実施形態において、データパッケージが、査定システムに転送される(例えば、103)。転送は、いくつかの実施形態において、ユーザの介入なしに信号取得手順の完了後に開始される。データリポジトリ112は、いくつかの実施形態において、患者のデータパッケージ、例えば、Amazon Simple Storage Service(すなわち、「Amazon S3」)に安全な、冗長な、クラウドベースのストレージを提供することができるクラウドストレージサービス上でホストされる。生物物理学的信号キャプチャシステム102aはまた、施術者が不適切な信号取得に対する通知を受信して、施術者に患者から直ちに追加のデータを取得するように警告するためのインターフェースを提供する。
【0039】
例示的な演算方法
図3A~
図3Bは各々、診断、治療、監視、又は追跡のための実際の用途において視覚的特徴又はそれらの中間出力を使用するための例示的な方法を示す。
【0040】
疾患状態又は兆候コンディションの存在の推定。
図3Aは視覚的パラメータ又は特徴を採用し、例えば、診断、追跡、又は治療を支援するために、疾患状態、医学的コンディション、又はいずれかの兆候の存在の推定値を決定する、方法300aを示す。方法300aは、例えば、
図1及び2若しくは本明細書に記載の他の例に関連して説明されるように、患者から生物物理学的信号(例えば、心臓信号、フォトプレチスモグラフィック信号、バリストカルジオグラフィック信号)を取得するステップ(302)を含む。いくつかの実施形態において、取得された生物物理学的信号は、リモートストレージ及び分析のために送信される。他の実施形態において、取得された生物物理学的信号は、ローカルに記憶され、分析される。
【0041】
上述のように、心臓の文脈における一例は、異常な左心室拡張末期圧(LVEDP)又は平均肺動脈圧(mPAP)、有意な冠動脈疾患(CAD)、異常な左心室駆出率(LVEF)、及び肺動脈高血圧(PAH)などの1つ以上の形態の肺高血圧(PH)の存在の推定である。推定され得る他の病態又は兆候コンディションとしては、例えば、本明細書に開示される様々な他の疾患及び医学的コンディションの中で、例えば、駆出率が保持された心不全(HFpEF)、不整脈、うっ血性心不全、弁不全などの1つ以上の形態の心不全が挙げられる。
【0042】
方法300aは、データセットを取り出し、視覚的特徴の値を決定するステップ(304)を更に含む。視覚的特徴の値を決定するための例示的な動作は、本明細書で後に考察される
図4~
図11に関連して提供される。方法300aは、決定された視覚的特徴の推定モデル(例えば、MLモデル)への適用に基づいて、疾患状態、医学的コンディション、又はいずれかの兆候の存在に関する推定値を決定するステップ(306)を更に含む。例示的な実装形態が、
図12A及び
図12Bに関して提供される。
【0043】
方法300aは、例えば、
図1、
図12A、及び
図12B、並びに本明細書に説明される他の例に関連して説明されるように、報告(例えば、疾患状態、医学的コンディション、又はいずれかの兆候の診断又は治療に使用される)内の疾患状態又は異常コンディションの存在に関する推定値(複数可)を出力するステップ(308)を更に含む。
【0044】
診断又はコンディション監視又は追跡。
図3Bは、医療機器又は健康監視デバイスの健康又は制御を監視するための視覚的パラメータ又は特徴を用いる方法300bを示す。方法300bは、患者から生物物理学的信号(例えば、心臓信号、フォトプレチスモグラフィック信号、バリストカルジオグラフィック信号など)を得るステップ(302)を含む。動作は、例えば、報告のための出力を提供するために、又は医療機器若しくは健康監視デバイスのための制御として、連続的又は断続的に実施され得る。
【0045】
方法300bは、例えば、
図4~
図11に関連して説明したように、取得された生物物理学的データセットから視覚的特徴値(複数可)を決定すること(310)を更に含む。
【0046】
方法300bは、視覚的特徴又はパラメータの値(複数可)を(例えば、診断で使用するための報告において、又は制御のための信号として)出力すること(312)を更に含む。監視及び追跡するために、出力は、ウェアラブルデバイス、ハンドヘルドデバイス、又は医療診断機器(例えば、パルス酸素濃度計システム、ウェアラブル健康監視システム)を介して、健康に関連付けられた拡張データを提供することができる。いくつかの実施形態では、出力は、蘇生システム、心臓又は肺ストレス試験機器、ペースメーカーにおいて使用されてもよい。
【0047】
視覚的予測器特徴
図4及び
図5は各々、例示的な実施形態による、1つ以上の視覚予測器関連特性の値を決定するように構成された、合計で2つの例示的なモジュールに対する、例示的な視覚的特徴査定モジュールを示す。特に、
図4の視覚的ループ特徴査定モジュール400は、研究中の波形領域の形態及び変動性を査定するために、各々位相空間モデルで表される、(i)心臓信号内の関心のある波形領域(心房脱分極(AD)波形領域、心室再分極(VR)波形領域、心室脱分極(VR)波形領域として示される)及び(ii)PPG信号(PPG信号#1及び#2として示される)に対する視覚的ループ特徴及びそれらの正射影から生成された特徴の値を計算することができる。ループ特徴は、点群モデルの外周、表面積、体積、曲率、及び渦度、又は関心のある波形の各々について点群から生成されたソリッドモデルに関連付けられた特徴を含むことができる。心臓信号の場合、これらは、AD、VD、及びVR波形領域などの特定の波形領域を含むことができ、PPG信号の場合、これらは、信号全体を含むことができる。
【0048】
図5の視覚ベクトル特徴査定モジュール500は、研究中の波形領域の形態及び変動性を査定するために、各々位相空間モデルで表される、(i)心臓信号内の関心のある波形領域(AD、VD、及びVR波形領域として示される)及び(ii)PPG信号(PPG信号#1及び#2として示される)に対する視覚ベクトル特徴及びそれらの正射影から生成された特徴の値を計算することができる。心臓信号の場合、これらは、AD、VD、及びVR波形領域などの特定の波形領域を含むことができ、PPG信号の場合、これらは、信号全体を含むことができる。
【0049】
実施例#1-視覚的ループ特徴
図4は、2つの例示的な視覚的特徴カテゴリのうちの第1のカテゴリとして、点群の外周、表面積、体積、曲率、及び/又は渦度、又は関心のある波形から生成された点群から生成されたソリッドモデルなどの視覚的特性を特徴付ける視覚的特徴又はパラメータの出力値を決定するように構成された例示的な視覚的ループ特徴査定モジュール400を例解する。2Dソリッドモデル又は投影の象限に対して追加の特徴を決定することができる。心臓信号の場合、心房脱分極(AD)波形領域、心室再分極(VR)波形領域、及び心室脱分極(VD)波形領域などの、関心のある波形としての特定の波形領域の視覚的特性を計算することができる。PPG信号の場合、取得された波形について視覚的特性を計算することができる。点群又はソリッドモデルは、(i)取得された心臓信号、又は(ii)PPG信号及びそのPPG信号から生成された関連するVPG及びAPG信号の3つのチャネルにおけるデータ点から生成することができる。
【0050】
表1は、生物物理学的信号又はその中の波形領域から抽出することができる12個の抽出された視覚的ループ特徴の例示的なセットを示す。一実施形態では、12個の視覚的特徴は、心臓又は生体電位信号の3つの波形領域のうちの一部又は全てに対して、及び2つのPPG信号に対して判定され、最大合計75個の特徴を提供することができる。表1において、記号「*」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、上昇LVEDPの存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。表1において、記号「**」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、冠状動脈疾患の存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。異常な又は上昇したLVEDPの存在又は不在及び有意なCADの存在又は不在の査定において有意な有用性を有すると判定された特定の特徴のリストを、それぞれ、表10A及び表10B、並びに表11A及び表11Bに提供する。
【表1-1】
【表1-2】
【0051】
表2は、表1の表面積及び体積ループ特徴の生物物理学的信号又はその中の波形領域から抽出され得る追加の4つのタイプの抽出された視覚的ループ特徴の例示的なセットを示す。表1の特徴で提供されるような3D位相空間モデル全体の特徴を決定するのではなく、表2の特徴は、3D位相空間モデルのサブ領域(オクタントと称される)の特徴を決定する。一実施形態では、表1の4つの視覚的特徴は、心臓又は生体電位信号の3つの波形領域のうちの一部又は全てに対する8つのオクタント領域に対して、及び2つのPPG信号に対して判定され、最大合計235個の特徴を提供することができる。表2において、記号「*」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、上昇LVEDPの存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。記号「**」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、冠状動脈疾患の存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。異常な又は上昇したLVEDPの存在又は不在の査定において有意な有用性を有すると判定された特定の特徴のリストは、それぞれ、表10A及び表10B並びに表11A及び表11Bに提供される。
【表2】
【0052】
表3は、外周、面積、曲率、及び循環を含む、表1の特徴の投影又はかかる投影の象限分析に基づいて、生物物理学的信号又はその中の波形領域から抽出され得る、追加の8つのタイプの抽出された視覚的ループ特徴の例示的なセットを示す。表3は、追加の特徴偏心を示す。一実施形態では、3つの投影が、表3の8つの特徴タイプの各々に対して判定されることができ、各投影に対して、5つの特徴(全体的値及び4象限内のその細分)が、合計135までの特徴を提供するように査定されることができる。表3において、記号「*」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、上昇LVEDPの存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。表3において、記号「**」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、冠状動脈疾患の存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。異常な又は上昇したLVEDPの存在又は不在及び有意なCADの存在又は不在の査定において有意な有用性を有すると判定された特定の特徴のリストを、それぞれ、表10A及び表10B並びに表11A及び表11Bに提供する。
【表3】
【0053】
図6は、例示的な実施形態による視覚的ループ特徴計算モジュール400の実装形態(600)を示し、これは、視覚的ループ特徴又はパラメータを生成するために全体的に又は部分的に使用することができ、それらの出力は、検討中の患者の生理系に関連付けられたメトリックを決定するために機械学習分類器において使用される。表1の例示的な特徴を決定するために、方法600は、いくつかの実施形態では、(i)信号準備(602)と、(ii)関心のある波形領域を分離すること(604)と、(iii)分離された信号の点群位相空間モデルを生成すること(606)と、(iv)視覚的特徴を決定することと、を含む。
図6は、関心のある3つのセグメント領域、すなわち、AD波形領域612、VD波形領域614、及びVR波形領域616について識別された複数の心周期データを含む例示的な心臓信号位相空間モデル610を示す。位相空間モデル610は、これらの3つの領域の一部ではない領域618も示す。
図6は、PPG信号622及びその関連するVPG信号624及びAPG信号626について生成されたPPG信号位相空間モデル620の例を示す。
図6はまた、収縮期PPGピーク628、拡張期PPGピーク630、ピークVPG632、最小VPG634、ピークAPG636、最小APG638、ベースAPG640、及び原点640を含む関心のある基準点を有するPPG信号位相空間モデル620を示す。追加の説明を以下に提供する。
【0054】
外周特徴。
図7A及び
図7Bは、心臓信号及びPPG信号の外周特徴を計算するための例示的な動作を示す。
図7Aに示される例では、心臓信号データに関して、モジュール400は、点群データのセットを生成することができ、そのセットの各々は、単一の心周期に関するデータと、関心のある波形領域、例えば、AD波形領域612、VD波形領域614、及びVR波形領域616に関するデータとを含む(
図7Aでは、VD波形領域614(614aとして示される)のみに関する例が示される)。次いで、所与の波形領域の外周値を、その波形領域の点群データ内の連続する点の各対の間のユークリッド距離の和として決定することができる。複数の心周期が取得されるため、モジュール400は、各周期の外周特徴値を決定し、次いで、表1に関連して図示及び説明されるように、そのセットのモード、平均、又は中央値、例えば、最大3つの特徴を選択することができる。LVEDP査定のために使用された例示的特徴では、モードが使用された。
【0055】
図7Bに示される例では、PPG信号に対して、モジュール400は、各々、単一PPGサイクルに対する点群データ620(620aとして示される)のセットを生成することができる。次いで、所与の波形領域の外周値は、表1に関連して示され、説明されるように、そのセットのモード、平均、又は中央値が、例えば、追加の2つの特徴まで選択され得る、その点群データ内の連続点の各対間のユークリッド距離の和として判定されることができる。
【0056】
加えて、モジュール400は、3つの直交平面(XY、XZ、及びYZ)に基づいて、点群データ620の投影の外周を計算することができる。心臓信号位相空間の場合、直交平面は、チャネルX及びY、チャネルX及びZ、並びにチャネルY及びZ(本明細書では「ORTH12」、「ORTH13」、及び/又は「ORTH23」平面とも称される)として表すことができる。PPG信号位相空間に関して、直交平面は、PPG及びVPG、PPG及びAPG、並びにVPG及びAPG(本明細書では「PPG_VPG」、「PPG_APG」、及び「PPG_APG」平面とも称される)として表されることができる。
【0057】
図7C及び
図7Dは、心臓信号データ又はPPG信号内の波形領域(例えば、VD波形領域、VR波形領域、AD波形領域)の点群データの投影を3つの直交平面(ORTH12(702)、ORTH13(704)、及びORTH23(706))上に投影するための例示的な動作を示す。
図7Cに示される例では、3つの正射影ORTH12(708)、ORTH13(710)、及びORTH23(712)は、VD波形領域(614a)の3D位相空間データから生成される(VR波形領域及びAD波形領域に対して同様の動作を実行することができる)。
図7Dに示される例では、最大3つ、例えば、ORTH12(708)、ORTH13(710)、及びORTH23(712)の正射影が、PPG信号の3D位相空間データから生成されることができる(類似動作が、他のPPG信号に対して行われることができる)。モジュール400は、3つの正射影の各々に対して、心臓信号データ又は2PPG信号データの3つの波形領域の各サイクルに対する2dperimeter特徴値を決定し、次いで、表3に関して説明されるように、例えば、追加の15個の特徴まで、そのセットのモード、平均、又は中央値を選択することができる。
【0058】
加えて、モジュール400は、3つの直交平面投影の特定の象限の外周を計算することができる。
図7Eは、心臓信号データ又はPPG信号内の波形領域(例えば、VD波形領域、VR波形領域、AD波形領域)の正射影を、3つの直交平面(ORTH12(702)の象限[1..4](716参照)、ORTH13(704)の象限[1..4](716参照)、及びORTH23(706)の象限[1..4](716参照))の各々に対する象限領域716上にセグメント化するための例示的な動作を示す。
図7Eに示される例では、VD波形領域710の正射影(710aとして示される)は、象限718、720、722、及び724に分割されて示される(同様の動作が、VD波形領域の正射影の各々に対して、並びにVR及びAD波形領域並びに2つのPPG信号に対して行われ得る)。したがって、モジュール400は、3つの正射影の象限(4)の各々に対して、心臓信号データ又は2PPG信号データの3つの波形領域の各サイクルに対する2dperimeter特徴値を決定し、次いで、表3に関して説明されるように、例えば、追加の60個の特徴まで、そのセットのモード、平均、又は中央値を選択することができる。
【0059】
表面積、体積、面積投影特徴。
図8A及び
図8Bは、心臓信号及びPPG信号の表面積及び体積特徴を計算するための例示的な動作を示す。
図8Aに示される例では、心臓信号データに対して、モジュール400は、(i)複数の心周期にわたって、関心のある波形領域、例えば、AD波形領域612、VD波形領域614、及びVR波形領域616に対する点群データ802(
図8Aでは、例は、VD波形領域614のみに対して示される)、及び/又は(ii)信号、例えば、PPG信号#1及びPPG信号#2に対する点群データ808を生成することができる。次いで、モジュール400は、点群データ(802、808)からアルファ形状オブジェクト(804、810)又は凸包オブジェクト(806、812)を生成することができる。
【0060】
心臓ループ及びPPG信号ループは、三次元空間において異なる軌道を有することが観察され、軌道は、例えば、固有の生物学的変動に起因して、拍動間又は周期間でわずかに変化する。したがって、心臓及びPPGループは、データ点が集中し、三次元ループを形成する位相空間内の領域を作成することができる。
図8Aは、信号データの32番目のウィンドウに対するVDループ802を示す。モジュール400は、アルファ形状オブジェクト804(例えば、全ての点がカバーされることを保証する最小値に自動的に設定されるアルファ半径を有する)として、及び凸包オブジェクト806として、ループデータ(例えば、802)をラップすることができる。
図8Bは、取得されたPPG信号のPPGループ808を示す。モジュール400は、アルファ形状オブジェクト810(例えば、全ての点がカバーされることを保証する最小値に自動的に設定されるアルファ半径を有する)として、及び凸包オブジェクト812として、ループデータ(例えば、808)をラップすることができる。
【0061】
次いで、モジュール400は、3D AlphaShape及びConvexHull形状の表面積及び体積を特徴値として計算することができ、例えば、表1に関連して図示及び説明されるように、心臓信号データ及び2PPG信号データの3つの波形領域に対して最大20の特徴を計算することができる。
【0062】
加えて、モジュール400は、表2に関連して図示及び説明されるように、特徴値として、八分円領域(例えば、8つの領域-表4参照)の各々に対する3D AlphaShape及びConvexHull形状の表面積及び体積、例えば、心臓信号データの3つの波形領域及び2つのPPG信号データに対する追加の160個までの特徴を計算することができる。オクタント領域の例示的なレイアウトを
図10Bに示す。表4は、オクタント識別領域を示す。
【表4】
【0063】
加えて、モジュール400は、3つの直交平面投影の象限(4)の面積を計算することができる。
図7Eでは、心臓信号データ又は2PPG信号データにおける、波形領域(例えば、VD波形領域、VR波形領域、AD波形領域)の正射影の各象限の面積726、728、730、732、及び象限の総面積は、表3に関連して説明されるように、例えば、追加の75個の特徴まで決定され得る。各象限内の面積を計算するために、TrapZ関数など(例えば、Matlab(登録商標)によって製造される)を使用して、面積を計算することができる。総面積は、ポリエリアMatlab関数を使用して計算することができる。
【0064】
加えて、モジュール400は、表3に関連して説明されるように、例えば、追加の15個までの特徴を提供するために、最大面積を含有する3つの正射影データの中の特定の象限を判定することができる。
【0065】
最大適合平面面積特徴。
図9A及び
図9Bは、心臓信号及びPPG信号に対する最大適合平面面積特徴を計算するための例示的な動作を示す。
図9A及び
図9Bに示される例では、モジュール400は、(i)複数の心周期にわたって、関心のある波形領域、例えば、AD波形領域612、VD波形領域614、及びVR波形領域616の最大適合平面902(
図9Aでは、例は、VD波形領域614aのみについて示される)と、(ii)信号、例えば、PPG信号#1及びPPG信号#2の最大適合平面904とを生成することができる。所与の心臓ループ又はPPGループに関して、ループを包含する点は、必ずしも単一の3D平面内にあるとは限らない。最大適合平面は、所与の心臓又はPPGループの点を通る最良適合平面である。すなわち、最良適合平面は、ループが最小の曲率を有するか又は最も平坦である単一平面である。モジュール400は、次いで、表1に関して図示及び説明されるように、心臓又はPPGループの点(902、904)を最良適合平面上に、例えば、心臓信号データ及び2つのPPG信号データの3つの波形領域のための最大5つの特徴に投影することができる。
【0066】
曲率、渦度、循環、伝播速度、偏心の特徴。心臓ループ又はPPGループであるループ内の単一の点の曲率は、特定の点、所与の点の5ms前の点、及び所与の点の5ms前の点を通過する円の半径の逆数として定義することができる。これらの3つの点は、同一線上にない場合、円を一意に識別する。
【0067】
モジュール400は、表1に関連して示され、説明されるように、心臓信号データ及び2つのPPG信号データの3つの波形領域に対して、最大曲率及び平均曲率の両方、例えば、追加の10個までの特徴を計算することができる。
【0068】
加えて、3つの直交平面における2D投影の平均及び最大曲率が、例えば、心臓信号データ及び2PPG信号データの3つの波形領域に対して、追加の30までの特徴を提供するように決定され得る。
【0069】
渦度は、位相空間におけるループを構成する点の回転を表す。モジュール400は、点を原点に結合するベクトルのカールとして、ループ内の各点について渦度を計算することができる。その点のカールはベクトルとして表すことができ、ベクトルの大きさは、ループ内のその点における回転の大きさを表す。ループ内の各点におけるカールは、式2によって計算することができる。
【数1】
【0070】
式2において、
【数2】
は、原点を心臓ループ内の所与の点に接続するベクトルであり、それぞれ、
【数3】
は、心臓信号又はPPG信号のORTH1/PPG(x軸)、ORTH2/VPG(y軸)、及びORTH3/APG(z軸)の単位ベクトルである。モジュール400は、表1に関連して図示及び説明されるように、ループ内の全ての点にわたる平均カール、例えば、心臓信号データ及び2つのPPG信号データの3つの波形領域に対する最大5つの特徴を計算することができる。
【0071】
循環は、ループを構成する点の回転の巨視的尺度である(一方、渦度は、ループ内の点の局所回転の尺度であり、ループ内の各点で計算される)。言い換えると、循環は、全体としてのループの回転の表現である。これは、式3によって計算することができる。
【数4】
【0072】
式3において、<x.y>は、ベクトルxとyとの間のドット積を表し、Vは、ループ内の2つの連続する点の間の速度ベクトル
【数5】
を表す。加えて、dsは、
【数6】
であるループ内の2つの連続点間の方向距離ベクトルを表し、Mは、心臓ループ又はPPGループの単一サイクルのためのデータ点の数である。モジュール400は、表1に関連して示され、説明されるように、心臓信号データ及び2つのPPG信号データの3つの波形に対して、例えば、追加の5つの特徴まで、循環を計算することができる。
【0073】
加えて、3つの直交平面の各々における2D投影の循環が、例えば、心臓信号データ及び2PPG信号データの3つの波形領域に対して、最大で追加の15個の特徴を提供するように決定され得る。
【0074】
伝播速度は、(i)心臓ループの波形領域、又は(ii)PPGループの任意の2つの連続する点の間のユークリッド距離を、2つの点のサンプリング間の時間差で割ったものとして計算することができる。サンプリングレートは、全ての患者にわたって一定に保たれるため、2つの連続する点の間の距離は、伝播速度のインジケーターとして使用され得る。例示的な心臓信号の場合、時間差は、(1000Hzのダウンサンプリングされた信号の場合)1/1000の一定値である。PPG信号の場合、時間差は1/250の一定値である(250Hzのダウンサンプリングされた信号の場合)。モジュール400は、表1に関連して図示及び説明されるように、最大速度値及び平均速度値、例えば、心臓信号データ及び2つのPPG信号データの3つの波形領域に対する最大10個の特徴を計算することができる。
【0075】
偏心は、式4に従って、ループ上に投影された垂直長軸及び短軸の長さを使用して、円形からのループの偏差の程度を近似することができる。
e=sqrt(1-b2/a2)(式4)
【0076】
式4において、bはループの短軸(736)の長さであり、aはループの長軸(738)の長さである。
図7Eに示される例では、投影の偏心734が、心臓信号データのVD波形領域710aに対して示される。モジュール400は、例えば、最大15の特徴を提供するために、心臓信号データ又は2つのPPG信号の3つの波形領域の各々の3つの正射影の各々に対して偏心を計算することができる。
【0077】
実施例#2-視覚ベクトル特徴
図5は、2つの例示的な視覚的特徴カテゴリのうちの第2のカテゴリとして、ベクトル振幅、ベクトル方位角及び仰角関連特徴、特定のベクトルのオクタント位置、並びに3D位相空間モデルにおいて定義されたベクトルのセットの渦度などの3D位相空間における視覚的特性を特徴付ける視覚的特徴又はパラメータの出力値を決定するように構成された例示的な視覚的ループ特徴査定モジュール500を例解する。
図10A~
図10Fは、これらのベクトルを、心臓信号に対する最大心房脱分極ベクトル(MADV)(1002)、心房再分極ベクトル(ARV)(1004)、最大心室脱分極ベクトル(MVDV)(1006)、初期心室脱分極ベクトル(IVDV)(1008)、終末心室脱分極ベクトル(TVDV)(1010)、及び最大心室再分極ベクトル(MVRV)(1012)として示す。PPG信号データに関して、
図10Gは、収縮期、拡張期、パルスベースPPG目印ベクトル(それぞれ、1014、1016、1018)、VPGピーク、ベース、最小目印ベクトル(それぞれ、1020、1022、1024)、及びAPGピーク、ベース、最小目印ベクトル(それぞれ、1026、1028、1030)を含む、PPG、VPG、及びAPG波形において定義された基準点によって定義され得るベクトルを示す。これらのベクトルの定義を表6に提供する。3Dベクトルの2D投影に対して追加の特徴を決定することができる。
【0078】
表5は、心臓信号又はPPG信号などの生物物理学的信号から抽出されることができる、5つの抽出された視覚ベクトル特徴の例示的なセットを示す。一実施形態では、5つの視覚ベクトル特徴が、心臓又は生体電位信号の一部又は全部に対して、及び2つのPPG信号に対して決定され、最大合計82個の特徴を提供することができる。表5において、記号「*」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、上昇LVEDPの存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。表5において、記号「**」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、冠状動脈疾患の存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。異常な又は上昇したLVEDPの存在又は不在及び有意なCADの存在又は不在の査定において有意な有用性を有すると判定された特定の特徴のリストを、それぞれ、表10A及び表10B並びに表11A及び表11Bに提供する。
【表5】
【表6】
【0079】
表7は、直交平面への振幅の投影又は直交平面に対する角度に基づいて、心臓信号又はPPG信号などの生物物理学的信号から抽出されることができる、追加の6タイプの抽出された視覚ベクトル特徴の例示的なセットを示す。一実施形態では、6つの視覚ベクトル特徴が、心臓又は生体電位信号の一部又は全部に対して、及び2つのPPG信号に対して決定され、最大合計150個の特徴を提供することができる。表7において、記号「*」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、上昇LVEDPの存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。表7において、記号「**」で示される特徴は、少なくとも1つの心臓疾患、医学的コンディション、又はいずれかの兆候の存在又は不在の査定(例えば、冠状動脈疾患の存在又は不在の決定)において有意な有用性を有することが実験的に決定されている。異常な又は上昇したLVEDPの存在又は不在及び有意なCADの存在又は不在の査定において有意な有用性を有すると判定された特定の特徴のリストを、それぞれ、表10A及び表10B並びに表11A及び表11Bに提供する。
【表7】
【0080】
ベクトル増幅特徴。
図10A~
図10Gは、心臓信号及びPPG信号の振幅特徴を計算するための例示的な動作を示す。各ベクトル特徴は、心臓又はPPG信号のループ内の開始位置(例えば、心臓信号の場合、x
start、y
start、z
start)及び終了位置(例えば、x
end、y
end、z
end)によって指定されることができる。モジュール500は、それぞれ、ベクトルの二乗されたORTH1、ORTH2、及びORTH3成分(心臓の場合)の和の平方根、又はベクトルの二乗されたPPG、VPG、及びAPG成分(PPGの場合)の和の平方根を計算することによって、三次元ベクトルの大きさとしてベクトル振幅を決定することができる。モジュール500は、例えば、心臓信号データ及び2つのPPG信号データに対して最大22個の特徴を提供するように、表6の6つの心臓ベクトル及び8つのPPGベクトルに対するベクトル振幅値を決定することができる。
【0081】
加えて、モジュール500は、例えば、投影の直交成分の合計の平方根として、3つの直交平面内の所与のベクトルの投影の大きさを計算することができる。したがって、モジュール500は、例えば、心臓信号データ及び2つのPPG信号データに対して最大66個の特徴を提供するように、3つの直交平面に対する表6の6つの心臓ベクトル及び8つのPPGベクトルに対する2dAmplitude値を決定することができる。
【0082】
Vector Azimuth、Elevation Angle、2dAngle、VDVR_Angle、3dAngle_VectorPair、2dAngle_VectorPair特徴。ベクトル方位角は、心臓信号及びPPG信号の各ベクトルの事前定義された軸に対する3D位相空間内の3Dベクトルの2D投影の角度である。仰角は、所与のベクトルと、事前定義された直交平面(例えば、ORTH12平面)におけるベクトルの投影との角度である。特徴は、度(又はラジアン)で計算され得る。
図11は、ベクトル方位角1102及び仰角1104の図を示す。モジュール500は、例えば、心臓信号データ及び2つのPPG信号データに対して最大28個の特徴を提供するように、(i)6つの心臓ベクトルに対するベクトル方位角値及び仰角値と、(ii)表6の8つのPPGベクトルの仰角とを決定することができる。
【0083】
加えて、モジュール500は、3つの直交平面のうちの所与の直交平面について、投影ベクトルが水平軸となす角度を計算することができる。角度(シータ)の大きさは、度で計算することができる。したがって、モジュール500は、例えば、心臓信号データ及び2つのPPG信号データに対して最大66個の特徴を提供するように、3つの直交平面に対する表6の6つの心臓ベクトル及び8つのPPGベクトルに対する2dAngle値を決定することができる。
【0084】
加えて、モジュール500は、表6の特徴として、心室脱分極平面と心室再分極平面との間の角度(VDVR_Angle)を計算することができる。心室脱分極平面は、心室脱分極ループ内の全ての点について最小二乗残差を有する平面として定義することができる。
【0085】
加えて、モジュール500は、心臓信号のベクトル対間の3D角度(3dAngle_VectorPairs)を計算することができる。ベクトル対の例には、(i)MVDVとMVRVのベクトル対、(ii)IVDVとMVRVのベクトル対、(iii)TVDVとMVRVのベクトル対、(iv)IVDVとMVDVのベクトル対、及び(v)TVDVとMVDVのベクトル対が含まれる。モジュール500はまた、(i)ORTH12平面における2つのベクトルの投影間、(ii)ORTH13平面における2つのベクトルの投影間、及び(iii)ORTH23平面における2つのベクトルの投影間の3つの2D角度を計算することができる。
【0086】
したがって、モジュール500は、例えば、心臓信号データについて最大20個の特徴を提供するために、5つの心臓ベクトル対について3dAngle値及び3つの2dAngle値を決定することができる。
【0087】
オクタント位置及び渦度ベクトル特徴。オクタント位置は、心臓信号データについてのオクタント識別子(表4参照)による所与のベクトルの終点の位置である。モジュール500は、例えば、心臓信号データに対して最大6つの特徴を提供するために、表6の6つの心臓ベクトルに対するオクタント数を決定することができる。
【0088】
渦度は、式2に関して説明したように、位相空間におけるループを構成する点の回転を表す。モジュール500は、表6の5つのMADV、ARV、MVDV、IVDV、TVDVベクトルの渦度、例えば、心臓信号データの最大5つの特徴を決定することができる。
【0089】
2dAngle_Plane特徴。2dAngleは、所与の直交平面について投影ベクトルが水平軸となす角度であり、例えば度である(
図11B参照)。モジュール500は、例えば、心臓信号データに対して最大18の特徴を提供するために、3つの直交平面における表6の6つの心臓ベクトルに対する2dAngleを決定することができる。
【0090】
実施例#3-視覚的特徴
本明細書で、例えば、実施例#1及び実施例#2において考察される特徴に加えて、波形の境界における不連続性を除去して新しいクラスの特徴を生成するために、cos/sin符号化が適用され得る。位相空間におけるデータの「円形度」のために、境界において不連続性が形成されることがある。例えば、10.0度と359.9度は、境界上の遷移を考慮すると10.1度しか離れていない場合があるが、349.0度離れているように見え得る(359.9-10.0)。cos/sin符号化の使用は、1つの特徴から2つの特徴に拡張することを犠牲にして、その循環の明示的な符号化を容易にすることができる。
【0091】
例示的な計算方法
上述したように、
図6は、例示的な実施形態による視覚的ループ特徴計算モジュール400又は500の実装形態(600)を示し、これは、検討中の患者の生理系に関連付けられたメトリックを決定するために機械学習分類器において使用される視覚的ループ特徴又はパラメータ及びそれらの出力を生成するために、全体的に又は部分的に使用することができる。表1の例示的な特徴を決定するために、方法600は、いくつかの実施形態では、(i)信号準備(602)と、(ii)関心のある波形領域を分離すること(604)と、(iii)分離された信号の点群位相空間モデルを生成すること(606)と、(iv)視覚的特徴を決定することと、を含む。
図6は、関心のある3つのセグメント領域、すなわち、AD波形領域612、VD波形領域614、及びVR波形領域616について識別された複数の心周期データを含む例示的な心臓信号位相空間モデル610を示す。位相空間モデル610は、これらの3つの領域の一部ではない領域618も示す。
図6は、PPG信号622及びその関連するVPG信号624及びAPG信号626について生成されたPPG信号位相空間モデル620の例を示す。
図6はまた、収縮期PPGピーク628、拡張期PPGピーク630、ピークVPG632、最小VPG634、ピークAPG636、最小APG638、ベースAPG640、及び原点640を含む関心のある基準点を有するPPG信号位相空間モデル620を示す。
【0092】
心臓信号準備(602)。3つの直交信号を有する例示的な心臓信号の場合、モジュール400又は500は、(i)取得された信号を1kHzにダウンサンプリングし(例えば、デシメーションを介して)、(ii)信号の過渡部分(例えば、取得された信号からの最初の31)を除去し、(iii)ベースラインワンダを除去し(例えば、2番目の次の順方向-逆方向0.67Hzハイパスフィルタを介して)、(iv)信号をウィンドウ処理してサブ信号を生成する(例えば、連続するサブ信号間に15秒の重複を有する30秒の持続時間の8つのサブ信号に分割する)ことができる。
【0093】
心臓信号描出(604)。モジュール400又は500は、前処理された心臓信号を、心房脱分極に関連付けられた波形領域(「P波」とも称される)、心室脱分極に関連付けられた波形領域(QRS波又はピークとも称される)、及び心室再分極に関連付けられた波形領域(「T波」とも称される)に分離することができる。
図12Aは、前処理された心臓信号を関心のある波形領域に分離するためのモジュール400又は500の例示的な実装形態を示す。
図12Aに示す例では、モジュール400又は500は、(i)前処理された心臓信号を単一の時系列
【数7】
に変換し、(ii)変換された時系列
【数8】
の各サイクルにおけるピーク(例えば、心室脱分極)を検出し、(iii)変換された時系列
【数9】
の対応する心室脱分極(VD)オンセット及び心室脱分極(VD)オフセット領域を検出し、(iv)VD
off,iと次の連続するサイクルVD
onset,i+1との間の信号をセグメント化し、セグメントの2/3を心室再分極領域として割り当て、残りの1/3を心房脱分極(P波)領域として割り当てることができる。例えば、D.B.Dubin,Rapid Interpretation of EKG’s:An interactive course,6th ed.tampa:Cover Pub,2000で説明されている他のセグメント化を使用することもできる。
【0094】
心室脱分極は、極大値を検出するように構成された任意の数のピーク検出演算子(例えば、Mathworksによって製造されたMatlab(登録商標)の「findpeak」関数)を使用して、ステップ(ii)ごとに決定することができる。心室脱分極(VD)開始及び心室脱分極(VD)オフセット領域を決定するために、ステップ(iii)に従って、モジュール400又は500は、(a)ウェーブレット変換を使用して(例えば、連続1-Dモーレットウェーブレット又はガウス、メキシカンハット、スプライン、及びメイヤーウェーブレット(例えば、マザーウェーブレットとして)を使用して)チャネルごとに前処理信号を分解し、(b)ウェーブレットスペクトルをフィルタリングし(例えば、40~60Hzの帯域通過フィルタを介して)、(c)帯域通過フィルタリングされたウェーブレットスペクトルから時系列信号を生成し、(d)動的閾値を下回る累積パワーを有する再構成された信号内のインデックスを決定することができる(例えば、QRS開始及びオフセットは、各検出されたQRSピークの前及び後の第1の時間に、累積パワーが動的閾値、例えば、累積パワーの上位25パーセンタイルを下回るとき、時間インデックスに割り当てられる)。
【0095】
PPG信号準備(602)。フォトプレチスモグラフィック信号のセットに対して、モジュール400又は500は、(i)取得された信号を250Hzにダウンサンプリングし、(ii)信号の過渡部分(例えば、両方の信号からの第1の10秒及び最後の2秒)を除去し、(iii)PPG信号から速度プレシズモグラム(VPG)及び加速度プレシズモグラム(APG)を生成し(例えば、PPG信号内の連続点の振幅の差又はVPG信号内の連続点の振幅の差として)、(iv)2つの取得されたPPG信号間のPPG、VPG、及びAPG信号の信号値を正規化する(例えば、zスコア変換を使用して)ことができる。
【0096】
PPG信号描出(604)。モジュール400又は500は、前処理されたPPG、VPG、及びAPG信号を描出して、収縮期、拡張期、パルスベースPPG目印ベクトル(それぞれ1014、1016、1018)、VPGピーク、ベース、最小目印ベクトル(それぞれ1020、1022、1024)、及びAPGピーク、ベース、最小目印ベクトル(それぞれ1026、1028、1030)を識別することができる。例示的な方法が、PPG目印、VPG目印、及びAPG目印について、それぞれ表8A、8B、及び8Cに要約されている。
【表8】
【表9】
【表10】
【0097】
実験結果及び例
特徴セットを開発するためにいくつかの開発研究が行われており、更に、疾患、医学的コンディション、若しくはいずれかの兆候の存在若しくは不在、重症度、又は局在化を推定するために使用され得るアルゴリズムが行われている。ある研究では、異常又は上昇したLVEDPの非侵襲的査定のためのアルゴリズムが開発された。上述のように、異常又は上昇したLVEDPは、様々な形態の心不全のインジケーターである。別の開発研究では、冠動脈疾患の非侵襲的査定のためのアルゴリズム及び特徴が開発された。
【0098】
これら2つの開発研究の一部として、生物物理学的信号キャプチャシステムを使用し、
図2に関連して説明されるプロトコルに従って、成人ヒト患者から臨床データを収集した。対象は、信号取得後に心臓カテーテル留置(CAD及び異常LVEDP評価のための現在の「ゴールドスタンダード」試験)を受け、カテーテル留置結果をCADラベル及びLVEDP値の上昇について評価した。収集されたデータは、特徴/アルゴリズム開発のためのものと、それらの検証のためのものとの別々のコホートに層別化された。
【0099】
特徴開発フェーズ内では、心臓血管系の特性を表すことを意図された生体電位信号(本明細書で考察される心臓信号の例として)及び光吸収信号(本明細書で考察される血行動態又はフォトプレチスモグラフィックの例として)から分析フレームワーク内の特色を抽出するために、視覚的特徴を含む特徴が開発された。対応する分類器も、分類器モデル、線形モデル(例えば、弾性ネット)、決定木モデル(XGB分類器、ランダムフォレストモデルなど)、サポートベクトルマシンモデル、及びニューラルネットワークモデルを使用して開発され、上昇した又は異常なLVEDPの存在を非侵襲的に推定した。単変量特徴選択査定及び交差検証演算を実行して、関心のある特定の疾患兆候のための機械学習モデル(例えば、分類器)で使用する特徴を識別した。機械学習訓練及び査定の更なる説明は、2021年8月23日に出願された「Method and System to Non-Invasively Assess Elevated Left Ventricular End-Diastolic Pressure」と題された米国仮特許出願第63/235,960号に記載されており、その全体が参照により本明細書に組み込まれる。
【0100】
単変量特徴選択査定は、各々がt検定、相互情報量、及びAUC-ROC評価を使用して陰性及び陽性データセット対によって定義される多くのシナリオを査定した。t検定は、未知の分散を有する2つの集団からの2つのサンプル平均の間に差があるかどうかを決定することができる統計的検定である。ここで、t検定は、これらの群、例えば、正常LVEDP対上昇(LVEDPアルゴリズム開発のため)、CAD-対CAD+(CADアルゴリズム開発のため)における特徴の平均間に差異がないという帰無仮説に対して行われた。小さいp値(例えば、≦0.05)は、帰無仮説に対する強い証拠を示す。
【0101】
相互情報(MI)操作を行って、上昇した若しくは異常なLVEDP又は有意な冠動脈疾患のある特定の特徴への依存性を査定した。1より大きいMIスコアは、評価されている変数間のより高い依存性を示す。1未満のMIスコアは、そのような変数のより低い依存性を示し、ゼロのMIスコアは、そのような依存性がないことを示す。
【0102】
受信者動作特性曲線、すなわちROC曲線は、その識別閾値を変化させたときのバイナリ分類器システムの診断能力を例解する。ROC曲線は、様々な閾値設定で偽陽性率(FPR)に対して真陽性率(TPR)をプロットすることによって作成することができる。AUC-ROCは、受信者動作特性(ROC)曲線下の面積を定量化し、この面積が大きいほど、モデルは診断上有用である。ROC値及びAUC-ROC値は、95%信頼区間の下端が0.50より大きい場合に統計的に有意であると考えられる。
【0103】
表6は、単変量特徴選択査定において使用される陰性及び陽性データセット対の例示的なリストを示す。具体的には、表6は、20mmHg又は25mmHgを上回るLVEDP測定値を有するものとして定義される陽性データセットと、12mmHg未満のLVEDP測定値を有するか、又は正常LVEDP読取値を有すると判定された対象群に属するものとして定義された陰性データセットとを示す。
【表11】
【0104】
表10A及び表10Bはそれぞれ、臨床評価システムにおいて実行されるアルゴリズムにおいて、上昇LVEDPの存在及び不在を推定する際に有用性を有すると判定された視覚的特徴のリストを示す。表10A及び表10Bの特徴並びに対応する分類器は、上昇LVEDPを測定するためのゴールドスタンダードの侵襲的方法に匹敵する臨床的性能を有することが検証されている。
【表12】
【表13-1】
【表13-2】
【0105】
表11A及び表11Bはそれぞれ、臨床評価システムにおいて実行されるアルゴリズムにおいて、有意なCADの存在及び不在を推定する際に有用性を有すると判定された視覚的特徴のリストを示す。表11A及び表11Bの特徴並びに対応する分類器は、CADを測定するためのゴールドスタンダードの侵襲的方法に匹敵する臨床的性能を有することが検証されている。表11A及び表11B(並びに10A及び10B)内の特徴の各々は、概して、ベクトル(例えば、「AD」、「IVDV」、「MVDV」、「TVDV」、「VD」、「baseVPG」、「diasVPG」など)、特徴タイプ(例えば、「2Dangle」、「Azimuth」)、及び関連チャネルによって命名される。適用可能な場合、象限情報(例えば、「Quad1」、「Quad2」など)が、特徴タイプと関連付けられたチャネルとの間に提供される。
【表14】
【表15-1】
【表15-2】
【表15-3】
【0106】
ある特定の視覚的特徴が、上昇LVEDPの存在及び不在、又は有意なCADの存在及び不在を推定する際に臨床的有用性を有するという判定は、他の疾患、医学的コンディション、又は特に、限定ではないが、本明細書に説明される心臓疾患若しくはコンディション、又はいずれかの兆候の存在若しくは不在及び/又は重症度及び/又は局在化を推定する際に、これらの視覚的特徴又はパラメータ、並びに本明細書に説明される他の特徴の使用のための基礎を提供する。
【0107】
例示的な臨床評価システム
図12Aは、一実施形態による、分類器(例えば、機械学習分類器)を介して、患者又は対象の生理学的状態に関連付けられた1つ以上のメトリックを生成するように、他の特徴又はパラメータと共に視覚的特徴又はパラメータを非侵襲的に計算するように
図1のモジュールを実装する、例示的臨床評価システム1200(臨床及び診断システムとも称される)を示す。実際に、(例えば、
図1、
図4、及び
図5の)特徴モジュールは、一般に、システム(例えば、臨床評価システム1200)の一部とみなすことができ、システムでは、任意の数及び/又はタイプの特徴が、例えば、異なる構成の特徴モジュールを有する異なる実施形態を用いて、関心のある疾患状態、医学的コンディション、いずれかの兆候、又はそれらの組み合わせに対して利用され得る。これは、
図12Aに更に例解され、臨床評価システム1200は、疾患特異的アドオンモジュール1202(例えば、上昇LVEDP又はmPAP、CAD、PH/PAH、異常なLVEF、HFpEF、及び本明細書に記載される他のものを査定するための)が、システム1200の完全な動作を実現するために、単独で、又は複数の事例において、単一のプラットフォーム(すなわち、ベースシステム1204)と統合されることが可能である、モジュール設計である。モジュール性は、臨床評価システム1200が、同じ同期的に取得された生物物理学的信号及びデータセット、並びにベースプラットフォームを活用して、そのような疾患特有のアルゴリズムが開発されるときに、いくつかの異なる疾患の存在について査定するように設計されることを可能にし、それによって、試験及び認証の時間及びコストを低減する。
【0108】
様々な実施形態では、臨床評価システム1200の異なるバージョンは、関心のある所与の疾患状態(複数可)、医学的コンディション(複数可)、又は対象の兆候コンディション(複数可)のために構成され得る、異なる特徴計算モジュールを含有することを含むことによって、査定システム103(
図1)を実装してもよい。別の実施形態では、臨床評価システム1200は、2つ以上の査定システム103を含むことができ、そのエンジン103の分類器116に固有の異なるスコアを生成するために選択的に利用することができる。このようにして、より一般的な意味での
図1及び
図12のモジュールは、異なる及び/又は複数の対応する分類器116を有する異なる及び/又は複数のエンジン103が、所望のモジュールの構成に応じて使用され得る、モジュール式システムの1つの構成とみなすことができる。したがって、視覚的特定の特徴(複数可)を伴う、又は伴わない、
図1のモジュールの任意の数の実施形態が存在してもよい。
【0109】
図12Aでは、システム1200は、機械学習された疾患特異的アルゴリズムを使用して、1つ以上の生物物理学的信号データセット(例えば、110)を分析し、一例として、病態又は異常状態の上昇LVEDPの可能性について査定することができる。システム1200は、アルゴリズムを使用して推定スコアの分析及び提示を容易にし、医師がそのスコアを使用して、例えば、疾患状態、医学的コンディション、又はいずれかの兆候の存在又は不在について査定することを可能にするように、組み合わせて協働するように設計されるハードウェア及びソフトウェア構成要素を含む。
【0110】
ベースシステム1204は、各アドオンモジュール1202(疾患特有のアルゴリズムを含む)が、次いで、病態又は兆候コンディションを査定するためにインターフェースをとる、機能及び命令の基盤を提供することができる。
図12Aの例に示されるようなベースシステム1204は、ベース分析エンジン又は分析器1206、ウェブサービスデータ転送API1208(「DTAPI」1208として示される)、報告データベース1210、ウェブポータルサービスモジュール1213、及びデータリポジトリ111(112aとして示される)を含む。
【0111】
クラウドベースであり得るデータリポジトリ112aは、信号キャプチャシステム102(102bとして示される)からのデータを記憶する。生物物理学的信号キャプチャシステム102bは、いくつかの実施形態では、7チャネルリードセット及びフォトプレチスモグラム(PPG)センサが確実に取り付けられた(すなわち、取り外し可能ではない)単一ユニットとして設計される、再使用可能デバイスである。信号キャプチャシステム102bは、そのハードウェア、ファームウェア、及びソフトウェアと一緒に、ユーザインターフェースを提供し、その中に入力された患者特有のメタデータ(例えば、名前、性別、生年月日、医療記録番号、身長、及び体重など)を収集し、患者の電気及び血行動態信号を同期して取得する。信号キャプチャシステム102bは、メタデータ及び信号データを単一のデータパッケージとしてクラウドベースのデータリポジトリに直接安全に送信することができる。データリポジトリ112aは、いくつかの実施形態では、患者特有のデータパッケージを受け入れ、記憶し、分析エンジン又は分析器1206若しくは1214によるその読み出しを可能にするように構成された、セキュアなクラウドベースのデータベースである。
【0112】
ベース分析エンジン又は分析器1206は、(「SQA」モジュール1216を介して実行された)取得された信号の品質査定を実施することができる安全なクラウドベースの処理ツールであり、その結果は、ポイントオブケアにおいてユーザに通信することができる。ベース分析エンジン又は分析器1206はまた、取得された生物物理学的信号(例えば、110(
図1参照))の前処理(前処理モジュール1218を介して示される)を実施してもよい。ウェブポータル1213は、ヘルスケア提供者に彼らの患者の報告へのアクセスを提供するように設計された安全なウェブベースのポータルである。ウェブポータル1213の例示的な出力は、視覚化1236によって示される。報告データベース(RD)1212は、セキュアなデータベースであり、出力されたスコア(複数可)(例えば、118)及び関連する情報が患者の一般的な健康記録に統合され、保存されるように、病院内の又は医師がホストする、遠隔でホストする、又は遠隔の電子健康記録システム(例えば、Epic、Cerner、Allscrips、CureMD、Kareoなど)などの他のシステムとセキュアにインターフェースし、通信してもよい。いくつかの実施形態では、ウェブポータル1213は、電話を介して出力臨床情報を提供するためにコールセンターによってアクセスされる。データベース1212は、郵便、宅配便、手交などを介して配達される報告を生成することができる他のシステムによってアクセスされてもよい。
【0113】
アドオンモジュール1202は、ベース分析エンジン(AE)又は分析器1206と共に動作する第2の部分1214(本明細書では分析エンジン(AE)又は分析器1214とも称され、「AEアドオンモジュール」1214として示される)を含む。分析エンジン(AE)又は分析器1214は、所与の疾患特有のアルゴリズムの主要な機能ループ、例えば、特徴計算モジュール1220、分類器モデル1224(「アンサンブル」モジュール1224として示される)、及び外れ値査定及び拒否モジュール1224(「外れ値検出」モジュール1224として示される)を含み得る。ある特定のモジュール式構成では、分析エンジン又は分析器(例えば、1206及び1214)は、単一分析エンジンモジュール内に実装されてもよい。
【0114】
主機能ループは、(i)全ての必要な環境変数値が存在することを確実にするために実行環境を検証し、(ii)取得された生物物理学的信号を含む新しい信号キャプチャデータファイルを分析する分析パイプラインを実行して、疾患特有のアルゴリズムを使用して患者のスコアを計算するための命令を含み得る。分析パイプラインを実行するために、AEアドオンモジュール1214は、
図1に関連して説明されるような様々な特徴モジュール114及び分類器モジュール116のための命令を含み、実行し、患者の生理学的状態に関連付けられたメトリックの出力スコア(例えば、118)を判定することができる。AEアドオンモジュール1214内の分析パイプラインは、特徴又はパラメータを計算することができ(「特徴計算」1220として示される)、特徴に基づいて、外れ値対非外れ値の信号レベル応答に対する外れ値検出リターンを提供することによって、計算された特徴が外れ値であるかどうかを識別する(「外れ値検出」1222として示される)。外れ値は、(モジュール116の)分類器を確立するために使用される訓練データセットに関して査定され得る。AEアドオンモジュール1214は、特徴及び分類器モデルの計算された値を使用して、(例えば、分類器モジュール1224を介して)患者の出力スコア(例えば、118)を生成し得る。上昇LVEDPの推定のための評価アルゴリズムの例では、出力スコア(例えば、118)は、LVEDPスコアである。CADの推定では、出力スコア(例えば、118)はCADスコアである。
【0115】
臨床評価システム1200は、ウェブサービスDTAPI1208(いくつかの実施形態では、HCPPウェブサービスとも称され得る)を使用して、構成要素内及び構成要素にわたってデータを管理することができる。DTAPI1208を使用して、取得された生物物理学的データセットをデータリポジトリ112aから取り出し、信号品質分析結果をデータリポジトリ112aに格納することができる。DTAPI1208はまた、記憶された生物物理学的データファイルを取り出し、分析エンジン又は分析器(例えば、1206、1214)に提供するために呼び出されてもよく、患者信号の分析エンジンの分析の結果は、DTAPI1208を使用して、報告データベース1210に転送されてもよい。DTAPI1208はまた、ヘルスケア専門家による要求に応じて、ウェブポータルモジュール1213に対して所与の患者データセットを読み出すために使用されてもよく、これは、セキュアなウェブアクセス可能インターフェースにおける精査及び解釈のために、報告をヘルスケア施術者に提示してもよい。
【0116】
臨床評価システム1200は、視覚的特徴120及び特徴モジュール122の様々な他の特徴を記憶する、1つ以上の特徴ライブラリ1226を含む。特徴ライブラリ1226は、アドオンモジュール1202(
図12Aに示す)又はベースシステム1204(図示せず)の一部であってもよく、いくつかの実施形態では、AEアドオンモジュール1214によってアクセスされる。
【0117】
モジュールのモジュール性及び様々な構成の更なる詳細は、2021年8月19日に出願された「Modular Disease Assessment System」と題された米国仮特許出願第63/235,960号に提供されており、その全体が参照により本明細書に組み込まれる。
【0118】
モジュール式臨床評価システムの例示的な動作
図12Bは、例示的な実施形態による、
図12Aの臨床評価システム1200の分析エンジン又は分析器(例えば、1206及び1214)の動作及びワークフローの概略図を示す。
【0119】
信号品質査定/拒否(1230)。
図12Bを参照すると、基礎分析エンジン又は分析器1206は、分析パイプラインが実行されている間に、SQAモジュール1216を介して、取得された生物物理学的信号データセットの品質を査定する(1230)。査定の結果(例えば、合格/不合格)は、ユーザによる読み取りのために、信号キャプチャシステムのユーザインターフェースに直ちに返される。信号品質要件を満たす取得された信号データは、容認可能(すなわち、「合格」)とみなされ、更に処理され、AEアドオンモジュール1214によって、病態又は兆候コンディション(例えば、上昇LVEDP又はmPAP、CAD、PH/PAH、異常LVEF、HFpEF)に関連付けられるメトリックの存在についての分析を受ける。許容できないとみなされた取得信号は拒否され(例えば、「不合格」)、患者から追加の信号を直ちに得るようにユーザに通知するために、通知が直ちにユーザに送信される(
図2を参照)。
【0120】
基本分析エンジン又は分析器1206は、信号品質についての2組の査定、すなわち、電気信号についての1組の査定及び血行動態信号についての1組の査定を実行する。電気信号査定(1230)は、電気信号が十分な長さであること、高周波ノイズ(例えば、170Hz超)がないこと、及び環境からの電力線ノイズがないことを確認する。血行動態信号査定(1230)は、血行動態データセット内の外れ値の割合が事前定義された閾値を下回ること、及び血行動態データセットの信号がレール又は飽和される割合及び最大持続時間が事前定義された閾値を下回ることを確認する。
【0121】
特徴値計算(1232)。AEアドオンモジュール1214は、特徴抽出及び計算を実行して、特徴出力値を計算する。LVEDPアルゴリズムの例では、AEアドオンモジュール1214は、いくつかの実施形態では、視覚的特徴(例えば、モジュール120で生成される)を含む、18個の異なる特徴ファミリー(例えば、モジュール120及び122で生成される)に属する合計446個の特徴出力を決定する。CADアルゴリズムの場合、AEアドオンモジュール1214の例示的な実装形態は、同じ18個の特徴ファミリーに対応する456個の特徴を含む特徴のセットを決定する。
【0122】
LVEDPアルゴリズムで使用されるものを含む、様々な特徴並びに他の特徴及びそれらの特徴ファミリーの追加の説明は、2021年8月23日に出願された「Method and System to Non-Invasively Assess Elevated Left Ventricular End-Diastolic Pressure」と題された米国仮特許出願第63/235,960号、2021年8月23日に出願された「Methods and Systems for Engineering Visual Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/236,072号、2021年8月23日に出願された「Methods and Systems for Engineering Power Spectral Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/235,963号、2021年8月23日に出願された「Method and System for Engineering Rate-Related Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/235,966号、2021年8月23日に出願された「Methods and Systems for Engineering Wavelet-Based Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/235,968号、「Method and System to Assess Disease Using Cycle Variability Analysis of Cardiac and Photoplethysmographic Signals」と題された米国仮特許出願第63/130,324号、2021年8月23日に出願された「Methods and Systems for Engineering photoplethysmographic Waveform Features for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/235,971号、2021年8月23日に出願され、「Methods and Systems for Engineering Cardiac Waveform Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/236,193号、2021年8月23日に出願され、「Methods and Systems for Engineering Conduction Deviation Features From Biophysical Signals for Use in Characterizing Physiological Systems」と題された米国仮特許出願第63/235,974号、に記載されており、それらの各々は参照によりその全体が本明細書に組み込まれる、
【0123】
分類器出力計算(1234)。AEアドオンモジュール1214は、次いで、分類器モデル(例えば、機械学習分類器モデル)中の計算された特徴出力を使用して、モデルスコアのセットを生成する。AEアドオンモジュール1214は、モデルスコアのセットを構成要素モデルのアンサンブルに結合し、これは、いくつかの実施形態では、LVEDPアルゴリズムの例における方程式6に示されるように、分類器モデルの出力を平均する。
【数10】
【0124】
いくつかの実施形態では、分類器モデルは、「Discovering Novel Features to Use in Machine Learning Techniques, such as Machine Learning Techniques for Diagnosing Medical Conditions」と題された米国特許出願公開第2019/0026430号、又は「Discovering Genomes to Use in Machine Learning Techniques」と題された米国特許出願公開第2019/0026431号に記載されているML技術に基づいて開発されたモデルを含むことができ、これらの各々は、その全体が参照により本明細書に組み込まれる。
【0125】
LVEDPアルゴリズムの例では、13個の機械学習分類器モデルが各々、計算された特徴出力を使用して計算される。13の分類器モデルは、4つのElasticNet機械学習分類器モデル[9]と、4つのRandomForestClassifier 機械学習分類器モデル[10]と、5つの極端勾配ブースティング(XGB)分類器モデル[11]とを含む。いくつかの実施形態では、年齢、性別、及びBMI値などの患者のメタデータ情報が、使用されてもよい。アンサンブル推定の出力は、連続スコアであってもよい。スコアは、ウェブポータル内の提示のための閾値を減算することによって、ゼロの閾値にシフトされてもよい。閾値は、感度と特異度との間のトレードオフとして選択され得る。閾値は、アルゴリズム内で定義され、試験陽性(例えば、「上昇LVEDPの可能性が高い」)及び試験陰性(例えば、「上昇LVEDPの可能性が低い」)コンディションの判定点として使用されてもよい。
【0126】
いくつかの実施形態では、分析エンジン又は分析器は、モデルスコアのセットを、ボディマス指数ベースの調整又は年齢若しくは性別に基づく調整と融合することができる。例えば、分析エンジン又は分析器は、
【数11】
の形態を有する患者BMIのシグモイド関数を用いてモデル推定を平均化することができる。
【0127】
医師ポータル視覚化(1236)。患者の報告は、取得された患者データ及び信号並びに疾患分析の結果の視覚化1236を含み得る。分析は、いくつかの実施形態では、報告内の複数のビューで提示される。
図12Bに示される例では、視覚化1236は、スコア要約セクション1240(「患者LVEDPスコア要約」セクション1240として示される)と、閾値セクション1242(「LVEDP閾値統計」セクション1242として示される)と、頻度分布セクション1244(「頻度分布」セクション1208として示される)とを含む。ヘルスケア提供者、例えば医師は、報告をレビューし、それを解釈して、疾患の診断を提供するか、又は治療計画を生成することができる。
【0128】
ヘルスケアポータルは、所与の患者の取得された信号データセットが信号品質基準を満たす場合、患者についての報告をリスト化することができる。報告は、信号分析が実行され得る場合に利用可能である疾患特有の結果(例えば、上昇LVEDP)を示すことができる。疾患特異的分析のための患者の推定スコア(視覚的要素118a、118b、118cを介して示される)は、確立された閾値に対して解釈され得る。
【0129】
図12Bの例に示されるスコア要約セクション1240では、患者のスコア118a及び関連付けられた閾値は、2トーンカラーバー(例えば、セクション1240に示される)上に重畳され、閾値は、バーの中心に位置し、定義された値「0」は、試験陽性と試験陰性との間の描出を表す。閾値の左側は、明るい陰影付きの光であってもよく、陰性試験結果(例えば、「上昇LVEDPの可能性が低い」)を示し、一方、閾値の右側は、陽性試験結果(例えば、「上昇LVEDPの可能性が高い」)を示すために暗い陰影付きであってもよい。
【0130】
閾値セクション1242は、患者スコア(例えば、118)の推定のための感度及び特異度を定義する検証集団に提供される閾値の報告された統計を示す。閾値は、個々の患者のスコア(例えば、118)にかかわらず、全ての試験について同じであり、これは、全てのスコア(陽性又は陰性)が、提供された感度及び特異性の情報を考慮して、正確さについて解釈され得ることを意味する。スコアは、臨床評価の更新と同様に、所与の疾患特異的分析について変化し得る。
【0131】
頻度分布セクション1244は、2つの検証集団(例えば、(i)偽陽性推定の可能性を示す上昇していない集団、及び(ii)偽陰性推定の可能性を示す上昇した集団)における全ての患者の分布を例解する。グラフ(1246、1248)は、試験性能検証集団患者に対する患者のスコア118(例えば、118b、118c)を解釈するためのコンテキストを提供するために、滑らかなヒストグラムとして提示される。
【0132】
頻度分布セクション1240は、その疾患、コンディション、又は兆候が存在しない検証集団の分布内の、疾患、コンディション、又は兆候が存在しない可能性を示すスコア(118b)を示す第1のグラフ1246(「非上昇LVEDP集団」1246として示される)と、その疾患、コンディション、又は兆候が存在する検証集団の分布内の、疾患、コンディション、又は兆候が存在する可能性を示すスコア(118c)を示す第2のグラフ1248(「上昇LVEDP集団」1248として示される)とを含む。上昇LVDEPの査定の例では、第1のグラフ1246は、真陰性(TN)及び偽陽性(FP)面積を識別する検証集団の非上昇LVEDP分布を示す。第2のグラフ1248は、偽陰性(TN)及び真陽性(FP)領域を識別する検証集団の上昇LVEDP分布を示す。
【0133】
頻度分布セクション1240はまた、検証集団群内の他の患者に対する患者のスコアの解釈テキストを(割合として)含む。この例では、患者は、LVEDP閾値の左側に位置する-0.08のLVEDPスコアを有し、患者が「上昇LVEDPの可能性が低い」を有することを示す。
【0134】
報告は、例えば、医師又はヘルスケア提供者が左心不全の兆候について診断する際に使用するために、ヘルスケアポータルに提示することができる。兆候は、いくつかの実施形態では、疾患、医学的コンディション、又はいずれかの兆候の存在についての確率又は重症度スコアを含む。
【0135】
外れ値査定及び拒否検出(1238)。AEアドオンモジュール1214が(プロセス1232において)特徴値出力を計算した後、(プロセス1234において)それらを分類器モデルに適用する前に、AEアドオンモジュール1214は、いくつかの実施形態では、特徴値出力の(プロセス1238に示される)外れ値分析を実施するように構成される。外れ値分析評価プロセス1238は、いくつかの実施形態では、機械学習外れ値検出モジュール(ODM)を実行して、検証及び訓練データから生成された特徴値を参照して異常な特徴出力値を識別及び除外することによって、異常な取得された生物物理学的信号を識別及び除外する。外れ値検出モジュールは、観察の残りからの分布外である隔離された領域における疎なクラスタ内に存在する外れ値について査定する。プロセス1238は、外れ値信号が分類器モデルに不適切に適用され、患者又はヘルスケア提供者によって見られる不正確な評価を生成するリスクを低減することができる。外れ値モジュールの精度は、ODMが許容可能な外れ値検出率(ODR)一般化を用いてテストセット内の全てのラベル付けされた外れ値を識別することができるホールドアウト検証セットを使用して検証されている。
【0136】
方法及びシステムは、ある特定の実施形態及び特定の例に関連して説明されてきたが、本明細書の実施形態は、制限的ではなく例解的であることが全ての態様で意図されるため、範囲は、示される特定の実施形態に限定されることは意図されない。本明細書で考察される視覚的特徴は、最終的に、類似又は他の開発アプローチを使用して、例えば、冠動脈疾患、肺高血圧、及び本明細書に説明されるような他の病態などの他の疾患、医学的コンディション、又はいずれかの兆候の存在若しくは不在及び/又は重症度の非侵襲的診断若しくは判定を行うために、あるいはそれを行う際に医師若しくは他のヘルスケア提供者を支援するために採用されてもよい。加えて、視覚的特徴を含む、例示的分析は、他の心臓関連病態及び兆候コンディション並びに神経関連病態及び兆候コンディションの診断及び治療において使用されることができ、そのような査定は、生物物理学的信号が生体の任意の関連システムに関与する、任意の病態又は兆候コンディションの診断及び治療(外科手術、低侵襲、及び/又は薬理学的治療を含む)に適用されることができる。心臓の状況における一例は、CAD、及び本明細書に開示される他の疾患、医学的コンディション、又は兆候コンディションの診断、並びに冠状動脈におけるステントの配置、アテレクトミーの実施、血管形成術、薬物療法の処方、及び/又は運動、栄養及び他の生活様式の変更の処方などのような、任意の数の治療単独若しくは組み合わせによるその処置である。診断され得る他の心臓関連の病態又は兆候コンディションとしては、例えば、不整脈、うっ血性心不全、弁不全、肺高血圧症(例えば、肺動脈高血圧、左心疾患に起因する肺高血圧、肺疾患に起因する肺高血圧、慢性血餅に起因する肺高血圧、及び他の疾患(例えば、血液又は他の障害)に起因する肺高血圧症)、他の心臓関連病態、兆候コンディション及び/又は疾患が挙げられる。診断され得る神経関連疾患、病態、又は兆候コンディションの非限定的な例としては、例えば、てんかん、統合失調症、パーキンソン病、アルツハイマー病(及び他の全ての形態の認知症)、自閉症スペクトラム(アスペルガー症候群を含む)、注意欠陥多動性障害、ハンチントン病、筋ジストロフィー、うつ病、双極性障害、脳/脊髄腫瘍(悪性及び良性)、移動障害、認知機能障害、発話障害、様々な精神病、脳/脊髄/神経損傷、慢性外傷性脳症、クラスタ頭痛、偏頭痛、神経障害(その様々な形態には、末梢神経性疾患を含む)、幻肢痛、慢性疲労症候群、急性及び/又は慢性疼痛(背痛、脊椎手術後痛症候群などを含む)、ジスキネジア、不安症、感染症又は外来因子によって引き起こされる兆候コンディション(例えば、ライム病、脳炎、狂犬病)、ナルコレプシー及び他の睡眠障害、外傷後ストレス障害、脳卒中に関連する神経学的兆候コンディション/影響、動脈瘤、出血性損傷など、耳鳴及び他の聴覚関連の疾患/兆候コンディション、並びに視覚関連の疾患/兆候コンディションなどが挙げられる。
【0137】
加えて、本明細書に記載の臨床評価システムは、とりわけ、心電図(ECG)、脳波(EEG)、ガンマ同期性、呼吸機能信号、パルスオキシメトリー信号、灌流データ信号、準周期的生物学的信号、胎児ECG信号、血圧信号、心臓磁場信号、心拍数信号などの生物物理学的信号を分析するように構成され得る。
【0138】
例示される方法及びシステムと共に使用され得る処理の更なる例は、米国特許第9,289,150号、同第9,655,536号、同第9,968,275号、同第8,923,958号、同第9,408,543号、同第9,955,883号、同第9,737,229号、同第10,039,468号、同第9,597,021号、同第9,968,265号、同第9,910,964号、同第10,672,518号、同第10,566,091号、同第10,566,092号、同第10,542,897号、同第10,362,950号、同第10,292,596号、同第10,806,349号、米国特許公開第2020/0335217号、同第2020/0229724号、同第2019/0214137号、同第2018/0249960号、同第2019/0200893号、同第2019/0384757号、同第2020/0211713号、同第2019/0365265号、同第2020/0205739号、同第2020/0205745号、同第2019/0026430号、同第2019/0026431号、PCT国際公開第WO2017/033164号、同第WO2017/221221号、同第WO2019/130272号、同第WO2018/158749号、同第WO2019/077414号、同第WO2019/130273号、同第WO2019/244043号、同第WO2020/136569号、同第WO2019/234587号、同第WO2020/136570号、同第WO2020/136571号、米国特許出願第16/831,264号、同第16/831,380号、同第17/132869号、PCT国際出願第PCT/IB2020/052889号、同第PCT/IB2020/052890号、に記載されており、これらの各々は、参照によりその全体が本明細書に組み込まれる。
【国際調査報告】