JP2020156609A5 - - Google Patents
Download PDFInfo
- Publication number
- JP2020156609A5 JP2020156609A5 JP2019057000A JP2019057000A JP2020156609A5 JP 2020156609 A5 JP2020156609 A5 JP 2020156609A5 JP 2019057000 A JP2019057000 A JP 2019057000A JP 2019057000 A JP2019057000 A JP 2019057000A JP 2020156609 A5 JP2020156609 A5 JP 2020156609A5
- Authority
- JP
- Japan
- Prior art keywords
- blood pressure
- waveform
- individual
- related information
- surge
- 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
Links
- 230000036772 blood pressure Effects 0.000 claims description 276
- 238000012545 processing Methods 0.000 claims description 75
- 238000000034 method Methods 0.000 claims description 70
- 206010005746 Blood pressure fluctuation Diseases 0.000 claims description 61
- 230000008569 process Effects 0.000 claims description 31
- 230000010349 pulsation Effects 0.000 claims description 9
- 101150104494 CAV1 gene Proteins 0.000 description 26
- 208000008784 apnea Diseases 0.000 description 21
- 201000002859 sleep apnea Diseases 0.000 description 15
- 230000035488 systolic blood pressure Effects 0.000 description 15
- 238000004891 communication Methods 0.000 description 12
- 238000012935 Averaging Methods 0.000 description 11
- 230000008452 non REM sleep Effects 0.000 description 11
- 230000036385 rapid eye movement (rem) sleep Effects 0.000 description 11
- 210000004204 blood vessel Anatomy 0.000 description 9
- 239000000463 material Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 206010021079 Hypopnoea Diseases 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 210000001367 artery Anatomy 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 208000024172 Cardiovascular disease Diseases 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 5
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 210000000056 organ Anatomy 0.000 description 5
- 230000007958 sleep Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 description 4
- 206010062519 Poor quality sleep Diseases 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000035487 diastolic blood pressure Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 101150047856 Cav2 gene Proteins 0.000 description 2
- 101150110214 Cav3 gene Proteins 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000004424 eye movement Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000004461 rapid eye movement Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 210000000707 wrist Anatomy 0.000 description 2
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 230000037007 arousal Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000005401 electroluminescence Methods 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 210000002321 radial artery Anatomy 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Description
この発明は血圧関連情報表示装置および血圧関連情報表示方法に関し、より詳しくは、被験者の血圧サージに関連する情報を可視化して表示する装置および方法に関する。また、この発明は、そのような血圧関連情報表示方法をコンピュータに実行させるためのプログラムに関する。
睡眠時無呼吸症候群(SAS;Sleep Apnea Syndrome)を罹患している患者においては、無呼吸後の呼吸再開時に血圧が急激に上昇し、その後に下降することが知られている。本明細書では、このような急激な血圧変動を「血圧サージ」(または単に「サージ」)と呼ぶ。患者に発生した血圧サージに関連する情報(例えば、血圧サージの波形)を可視化して表示することは、SASの診断や治療に役立つと考えられる。
従来、例えば特許文献1(WO 2017/082107 A1)の図3では、睡眠時無呼吸症候群(SAS)患者における変動状態(血圧サージ)の波形がグラフとして表示されている。
しかしながら、特許文献1では、血圧サージの波形は、収縮期血圧値(SBP)の時系列データについて時間軸を拡大して、脈拍毎のピークが分かる形態で表示されているに過ぎず、血圧の上昇、下降の傾向を示す曲線(形状)としては表示されていない。
このような状況下で、本出願人は、先に、SAS患者おける血圧サージの波形のパタン(形状)を分類し、可視化して表示する発明を出願した(特願2017-050066号)。これにより、医師は、比較的短い時間内に、各SAS患者おける血圧サージの波形のパタンを把握することができる。さらに、一歩進んで、医師が、各SAS患者における血圧サージの代表的波形と波形ばらつきを重ね合わせて見ることができれば便利であろう。また、SAS患者に限らず、被験者の血圧サージの代表的波形と波形ばらつきを重ね合わせて見ることは、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。
そこで、この発明の課題は、血圧サージの代表的波形と波形ばらつきを重ね合わせて表示できる血圧関連情報表示装置および血圧関連情報表示方法を提供することにある。また、この発明の課題は、そのような血圧関連情報表示方法をコンピュータに実行させるためのプログラムを提供することにある。
上記課題を解決するため、この開示の血圧関連情報表示装置は、
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示装置であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する血圧サージ検出部と、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める個別波形取得部と、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める統計処理部と、
表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う表示処理部とを備え、
上記統計処理部は、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
上記統計処理部は、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする
ことを特徴とする。
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示装置であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する血圧サージ検出部と、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める個別波形取得部と、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める統計処理部と、
表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う表示処理部とを備え、
上記統計処理部は、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
上記統計処理部は、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする
ことを特徴とする。
本明細書で、「予め定められた判定基準」とは、典型的には、睡眠時無呼吸症候群(SAS)患者において血圧サージを検出するための基準を指す。例えば、特願2017-048946号、特願2017-050066号に開示されているような、ピーク検出区間(例えば、15拍分の期間)内にサージ開始点の時刻からサージピーク点の時刻までが含まれるとともに、サージ開始点の時刻における収縮期血圧値とピーク点の時刻における収縮期血圧値との差(血圧変動量)が20mmHg(または15mmHg)以上であること、サージ開始点の時刻とピーク点の時刻との間の期間が5拍分よりも大きいこと、かつ、ピーク点の時刻とサージ終了点の時刻との間の期間が7拍分よりも大きいことを指す。
また、包絡線を作成するに用いられる、血圧サージをなす複数の「拍対応ピーク」とは、連続瞬時血圧波形の中での、収縮期血圧に相当するピークを指す。ただし、拡張期血圧値(DBP)に相当するピークであってもよい。
また、「統計処理」とは、個別波形を平均する処理、中央値を求める処理などを指す。血圧サージの「代表的波形」とは、例えば、複数の個別波形を平均して得られた平均的波形、または、複数の個別波形の中央値に相当する波形などを指す。また、血圧サージの「波形ばらつき」とは、例えば複数の個別波形の分布の幅を指す。
また、「表示画面」は、典型的には表示器の画面を指すが、例えば、プリンタによって出力される紙面であってもよい。
また、表示画面上の例えば代表的波形を示す曲線の「表示」は、典型的には、横方向座標が時間経過(例えば、拍の進行)を表し、縦方向座標が血圧サージに伴う血圧変動量を表す態様で行われる。
本明細書で、「個別波形のピーク」とは、各血圧サージのピークを指し、各血圧サージをなす複数の拍対応ピークのうちの最大の拍対応ピークに相当する。
また、「等価な状態」とは、上記座標平面上で上記複数の個別波形のピークの位置を実際に描いて揃える処理などを要さず、メモリ領域のデータが、上記座標平面上で規定された状態と等価な状態になっていれば足りるという意味である。
本明細書で、「個別波形のピーク」とは、各血圧サージのピークを指し、各血圧サージをなす複数の拍対応ピークのうちの最大の拍対応ピークに相当する。
また、「等価な状態」とは、上記座標平面上で上記複数の個別波形のピークの位置を実際に描いて揃える処理などを要さず、メモリ領域のデータが、上記座標平面上で規定された状態と等価な状態になっていれば足りるという意味である。
この開示の血圧関連情報表示装置では、血圧サージ検出部は、被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する。個別波形取得部は、検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める。統計処理部は、複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める。表示処理部は、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う。これにより、ユーザ(典型的には、医師、看護師などの医療関係者を指す。被験者であってもよい。以下同様。)は、上記表示画面を見ることによって、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。このことは、SASの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。
この血圧関連情報表示装置では、上記統計処理部は、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記複数の個別波形の血圧変動量データを統計処理するので、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記座標平面上で上記複数の個別波形の横方向寸法が互いに幾つかの拍数相当分異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。
一実施形態の血圧関連情報表示装置では、上記複数の個別波形の血圧変動量データは、それぞれ、その個別波形の開始点の血圧値を基準とした変動量であることを特徴とする。
この一実施形態の血圧関連情報表示装置では、上記代表的波形と上記波形ばらつきを把握し易くなる。
一実施形態の血圧関連情報表示装置では、上記複数の個別波形の血圧変動量データは、上記複数の個別波形のピークの高さを揃えるように正規化されていることを特徴とする。
この一実施形態の血圧関連情報表示装置では、上記代表的波形と上記波形ばらつきをさらに把握し易くなる。
一実施形態の血圧関連情報表示装置では、
上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められ、
上記統計処理部は、上記拍番号毎に上記複数の個別波形の血圧変動量データを統計処理することを特徴とする。
上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められ、
上記統計処理部は、上記拍番号毎に上記複数の個別波形の血圧変動量データを統計処理することを特徴とする。
ここで、心臓の拍動によって拍間隔は一定ではないため、或る個別波形をなす血圧サージ発生時の1拍分の期間と、別の個別波形をなす血圧サージ発生時の1拍分の期間とは、互いに相違する場合がある。「上記時間経過を表す横方向座標は…拍番号によって定められ」とは、その相違にかかわらず、上記時間経過が拍の進行によって表される、という意味である。
この一実施形態の血圧関連情報表示装置では、上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められる。したがって、或る個別波形をなす血圧サージ発生時の1拍分の期間と、別の個別波形をなす血圧サージ発生時の1拍分の期間との相違にかかわらず、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記統計処理部は、上記複数の個別波形の血圧変動量データを統計処理する際に、上記拍番号毎に上記血圧変動量データの統計処理を行えば足りる。したがって、横方向に関して連続的に血圧変動量データの統計処理を行う場合に比して、上記統計処理部の演算量が少なくて済む。
一実施形態の血圧関連情報表示装置では、
上記統計処理部が求める上記波形ばらつきは、kを自然数としたとき、上記複数の個別波形の血圧変動量データがなす分布の標準偏差の±k倍として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±k倍に相当する幅をもつ領域を表示することを特徴とする。
上記統計処理部が求める上記波形ばらつきは、kを自然数としたとき、上記複数の個別波形の血圧変動量データがなす分布の標準偏差の±k倍として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±k倍に相当する幅をもつ領域を表示することを特徴とする。
この一実施形態の血圧関連情報表示装置では、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±k倍(ただし、kは自然数である。典型的には、kは1、2または3として設定される。)に相当する幅をもつ領域が表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。この波形ばらつきを示す領域の表示は、例えば個別波形のばらつきが正規分布として取り扱われるのに十分である場合(例えば、個別波形の数が数十個以上の場合)に、特に有益となる。
一実施形態の血圧関連情報表示装置では、
上記統計処理部が求める上記波形ばらつきは、上記複数の個別波形の血圧変動量データがなす分布の四分位範囲として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域を表示することを特徴とする。
上記統計処理部が求める上記波形ばらつきは、上記複数の個別波形の血圧変動量データがなす分布の四分位範囲として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域を表示することを特徴とする。
この一実施形態の血圧関連情報表示装置では、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域が表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。また、上記波形ばらつきを示す領域を標準偏差を用いて定義する場合に比して、多くの計算を必要とせず、計算が簡単で済む。したがって、処理速度の向上やメモリの節約に貢献することができる。
一実施形態の血圧関連情報表示装置では、上記表示処理部は、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線とは区別可能な態様で、上記複数の個別波形を示す曲線を表示する処理を行うことを特徴とする。
この一実施形態の血圧関連情報表示装置では、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線と上記複数の個別波形を示す曲線とが表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。この複数の個別波形を示す曲線の表示は、例えば個別波形の数が数個以下の場合に、ユーザに個別波形のばらつきを把握させるのに有益となる。
一実施形態の血圧関連情報表示装置では、
上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態が特定された身体状態特定期間を表す情報を入力する入力部を備え、
上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求め、
上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うことを特徴とする。
上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態が特定された身体状態特定期間を表す情報を入力する入力部を備え、
上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求め、
上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うことを特徴とする。
本明細書で、「身体状態特定期間」とは、無呼吸期間、レム睡眠期間、ノンレム睡眠期間、覚醒期間、および/または、経皮的動脈血酸素飽和度(SpO2)が低い期間など、被験者が血圧サージの要因となり得る身体状態をとっている期間を指す。
この一実施形態の血圧関連情報表示装置では、入力部は、上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態特定期間を表す情報を入力する。上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める。上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う。これにより、ユーザは、上記身体状態特定期間毎に、言い換えれば、被験者が血圧サージの要因となり得る身体状態をとっている期間毎に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。したがって、ユーザは、例えば最も深刻な血圧サージを引き起こしている要因(身体状態)を把握することができる。
別の局面では、この開示の血圧関連情報表示方法は、
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示方法であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出するステップと、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求めるステップと、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求めるステップとを備え、
このステップでは、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
このステップでは、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにし、
さらに、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うステップを備える
ことを特徴とする。
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示方法であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出するステップと、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求めるステップと、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求めるステップとを備え、
このステップでは、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
このステップでは、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにし、
さらに、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うステップを備える
ことを特徴とする。
この開示の血圧関連情報表示方法によれば、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。このことは、SASの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。また、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記複数の個別波形の血圧変動量データを統計処理するので、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記座標平面上で上記複数の個別波形の横方向寸法が互いに幾つかの拍数相当分異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。
さらに別の局面では、この開示のプログラムは、上記血圧関連情報表示方法をコンピュータに実行させるためのプログラムである。
この開示のプログラムをコンピュータに実行させることによって、上記血圧関連情報表示方法を実施することができる。
以上より明らかなように、この開示の血圧関連情報表示装置および血圧関連情報表示方法によれば、血圧サージの代表的波形と波形ばらつきを重ね合わせて表示できる。また、この開示のプログラムをコンピュータに実行させることによって、上記血圧関連情報表示方法を実施することができる。
以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。
(システムの構成)
図1は、この発明の血圧関連情報表示装置をネットワーク上の一実施形態のシステム(全体を符号100で示す。)として構成した例を示している。このシステム100は、表示画面としての表示器240を有する病院端末200A,200Bと、サーバ300と、トノメトリ方式の血圧計400と、PSG(polysomnography;睡眠ポリグラフ検査)装置500とを含んでいる。これらの病院端末200A,200B、サーバ300、血圧計400およびPSG装置500は、病院内LAN(Local Area Network)であるネットワーク900を介して互いに通信可能になっている。このネットワーク900を介した通信は、無線、有線のいずれでも良い。この実施の形態において、ネットワーク900は、病院内LAN(Local Area Network)であるが、これに限定されず、インターネットのような他の種類のネットワークであってもよいし、USBケーブルなどを用いた1対1の通信であってもよい。なお、病院端末200A,200Bは、この例では2台のみを示しているが、3台以上設けられていてもよい。
図1は、この発明の血圧関連情報表示装置をネットワーク上の一実施形態のシステム(全体を符号100で示す。)として構成した例を示している。このシステム100は、表示画面としての表示器240を有する病院端末200A,200Bと、サーバ300と、トノメトリ方式の血圧計400と、PSG(polysomnography;睡眠ポリグラフ検査)装置500とを含んでいる。これらの病院端末200A,200B、サーバ300、血圧計400およびPSG装置500は、病院内LAN(Local Area Network)であるネットワーク900を介して互いに通信可能になっている。このネットワーク900を介した通信は、無線、有線のいずれでも良い。この実施の形態において、ネットワーク900は、病院内LAN(Local Area Network)であるが、これに限定されず、インターネットのような他の種類のネットワークであってもよいし、USBケーブルなどを用いた1対1の通信であってもよい。なお、病院端末200A,200Bは、この例では2台のみを示しているが、3台以上設けられていてもよい。
図2に示すように、病院端末200Aは、本体200Mと、この本体200Mに搭載された、制御部210と、メモリ220と、操作部230と、表示画面としての表示器240と、通信部290とを含む。この病院端末200Aは、市販のノート型パーソナルコンピュータからなり、後述の処理を行わせるようにアプリケーションソフトウェア(コンピュータプログラム)をインストールしたものであり、サーバ300にアクセスできるようになっている。
制御部210は、CPU(Central Processing Unit;中央演算処理装置)およびその補助回路を含み、病院端末200Aの各部を制御し、メモリ220に記憶されたプログラムおよびデータに従って後述の処理を実行する。すなわち、操作部230、および、通信部290から入力されたデータを処理し、処理したデータを、メモリ220に記憶させたり、表示器240で表示させたり、通信部290から出力させたりする。
メモリ220は、制御部210でプログラムを実行するために必要な作業領域として用いられるRAM(Random Access Memory)と、制御部210で実行するための基本的なプログラムを記憶するためのROM(Read Only Memory)とを含む。また、メモリ220の記憶領域を補助するための補助記憶装置の記憶媒体として、半導体メモリ(メモリカード、SSD(Solid State Drive))などが用いられてもよい。
操作部230は、この例では、キーボードおよびマウスで構成され、典型的にはユーザとしての医師による操作を示す操作信号を制御部210に入力する。また、操作部230は、キーボードおよびマウスに替えて、または、加えて、タッチパネルなどの他の操作デバイスで構成されるようにしてもよい。
表示器240は、表示画面(例えば、LCD(Liquid Crystal Display)またはEL(Electroluminescence)ディスプレイなど)を含む。表示器240は、制御部210によって制御されて、所定の映像を表示画面に表示させる。
通信部290は、制御部210からの情報をネットワーク900を介してサーバ300へ送信する。
簡単のため図示を省略するが、他の病院端末200B,…も病院端末200Aと同様の構成になっている。
図3に示すように、サーバ300は、制御部310と、記憶部320と、操作部330と、表示器340と、通信部390とを含む。このサーバ300は、汎用のコンピュータ装置に、後述の処理を行わせるようにプログラム(ソフトウェア)をインストールしたものである。
制御部310は、CPUおよびその補助回路を含み、サーバ300の各部を制御し、記憶部320に記憶されたプログラムおよびデータに従って所定の処理を実行し、操作部330、および、通信部390から入力されたデータを処理し、処理したデータを、記憶部320に記憶させたり、表示器340で表示させたり、通信部390から出力させたりする。
記憶部320は、制御部310でプログラムを実行するために必要な作業領域として用いられるRAMと、制御部310で実行するための基本的なプログラムを記憶するためのROMとを含む。記憶部320には、多くの被験者から送られてきた血圧測定データを含むデータベース321が設けられている。また、記憶部320の記憶領域を補助するための補助記憶装置の記憶媒体として、磁気ディスク(HD(Hard Disk)、FD(Flexible Disk))、光ディスク(CD(Compact Disc)、DVD(Digital Versatile Disk)、BD(Blu-ray Disc))、光磁気ディスク(MO(Magneto-Optical disk))、または、半導体メモリ(メモリカード、SSD)などが用いられてもよい。
操作部330は、この例では、キーボードおよびマウスで構成され、ユーザによる操作を示す操作信号を制御部310に入力する。また、操作部330は、キーボードおよびマウスに替えて、または、加えて、タッチパネルなどの他の操作デバイスで構成されるようにしてもよい。
表示器340は、表示画面(例えば、LCDまたはELディスプレイなど)を含む。表示器340は、制御部310によって制御されて、所定の映像を表示画面に表示させる。
通信部390は、制御部310からの情報をネットワーク900を介して他の装置(この例では病院端末200A)へ送信するとともに、他の装置からネットワーク900を介して送信されてきた情報を受信して制御部310に受け渡す。
図1中に示す血圧計400は、この例では、特願2017-050066号に開示されているようなトノメトリ方式の血圧計からなっている。この血圧計400は、図4中に示すように、トノメトリ方式で被験者90の被測定部位(例えば、左手首)を通る橈骨動脈の圧脈波を一拍毎に連続的に検出する圧力センサ部402と、この圧力センサ部402が検出した圧力の変化を血圧の時系列データとして出力する本体ユニット401とを備えている。圧力センサ部402と本体ユニット401との間は、信号ケーブル403によって接続されている。トノメトリ法は血管を圧力センサ部402(例えば、圧脈波センサ)で圧扁することにより圧脈波を計測し血圧を決定する手法である。血管の厚さが一様な円管と見なすと、血管内の血液の流れ、拍動の有無に関係なく血管壁を考慮してラプラスの法則に従い、血管の内圧(血圧)と血管の外圧(圧脈波の圧力)との関係式を導くことができる。この関係式で押圧面において血管が圧扁されている条件下では、血管の外壁及び内壁の半径を近似することにより、圧脈波の圧力と血圧とが等しいと近似できる。従って以降は、圧脈波の圧力は血圧と同一値になる。この結果、血圧計400は、被測定部位(例えば、左手首)の血圧値を一心拍ごとに測定して、例えば図6(A)に示すような態様で、測定時刻(時間)と血圧とを対応付けた血圧の時系列データ801を出力するようになっている。一晩の時系列データ801は、約3万拍分の拍対応ピーク(収縮期血圧(SBP)または拡張期血圧(DBP)に相当するピーク)を含んでいる。
また、図1中に示すPSG装置500は、この例では、市販のPSG装置(例えば、日本光電工業株式会社製のNeurofax(登録商標)EEG-9200)からなり、後述の[B.応用的な血圧関連情報表示方法]を実行する際にのみ用いられる。このPSG装置500は、図4中に示すように、センサ群502と、このセンサ群502からの信号を処理して被験者90の身体状態を特定する情報を出力する本体ユニット501とを備えている。この例では、センサ群502は、脳波を検出する脳波検出電極502Aと、眼球運動を検出する眼球運動検出電極502Bと、呼吸に伴う気流を検出する気流センサ502Cと、心電図を得るための心電図電極502Dと、経皮的動脈血酸素飽和度(SpO2)を検出するパルスオキシメータ502Eと、筋電図を得るための筋電図電極502Fとを含んでいる。センサ群502と本体ユニット501との間は、図示しないケーブルボックスを介して図示しない信号ケーブルによって接続されている。この例では、PSG装置500は、被験者90の身体状態が特定された期間(これを「身体状態特定期間」と呼ぶ。)を表す情報として、無呼吸期間(若しくは低呼吸期間)、レム睡眠期間、ノンレム睡眠期間、覚醒期間、および/または、SpO2が低い期間を表す情報を出力することができる。ここで、睡眠中の「無呼吸」とは、「呼吸が10秒以上止まる状態」を指す。また、「低呼吸」とは、呼吸による換気が10秒以上50%以下に低下することを指す。また、「レム睡眠」とは急速眼球運動を伴う睡眠を指し、「ノンレム睡眠」とは急速眼球運動を伴わない睡眠を指し、「覚醒」とは目覚めている状態を指す。SpO2とは、動脈血の中を流れている赤血球に含まれるヘモグロビンの何%に酸素が結合しているかを、皮膚を通して(経皮的に)測定した値を指す。SpO2が「低い」期間とは、この例ではSpO2が90%未満となった期間を指す。
例えば図6(A)の一部を拡大して示す図6(B)中に、「無呼吸」期間が、斜線付きのバー805,805,…で示されている。また、「低呼吸」期間が、白抜きのバー804,804,…で示されている。このような態様で、PSG装置500は、被験者90についての身体状態特定期間を表す情報を複数種類出力することができる。
(血圧関連情報表示方法)
このシステム100は、大別して、サーバ300での画像データの作成と、病院端末(例えば200A)での画像データの表示とを含む、以下のような[A.基本的な血圧関連情報表示方法]と[B.応用的な血圧関連情報表示方法]を実施することができる。
このシステム100は、大別して、サーバ300での画像データの作成と、病院端末(例えば200A)での画像データの表示とを含む、以下のような[A.基本的な血圧関連情報表示方法]と[B.応用的な血圧関連情報表示方法]を実施することができる。
[A.基本的な血圧関連情報表示方法]
この基本的な血圧関連情報表示方法では、PSG装置500を用いることなく、図6(A)に例示したような、血圧計400が出力する測定時刻(時間)と血圧とを対応付けた血圧の時系列データ801のみを用いて、次のように画像データの作成・表示を行う。
この基本的な血圧関連情報表示方法では、PSG装置500を用いることなく、図6(A)に例示したような、血圧計400が出力する測定時刻(時間)と血圧とを対応付けた血圧の時系列データ801のみを用いて、次のように画像データの作成・表示を行う。
(サーバ300での画像データの作成)
i) 図4に示した態様で、被験者90について、一晩(睡眠期間を含む。)にわたって、血圧計400を用いた連続的な血圧測定が行われたものとする(この例では、PSG装置500を用いた身体状態の測定は行われないものとする。)。この例では、血圧計400からの血圧の時系列データが、病院端末200Aによって受信され、メモリ220に一旦記憶される。続いて、例えばユーザとしての医師が病院端末200Aの操作部230を操作して、被験者90について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク900を介してサーバ300に送信する。
i) 図4に示した態様で、被験者90について、一晩(睡眠期間を含む。)にわたって、血圧計400を用いた連続的な血圧測定が行われたものとする(この例では、PSG装置500を用いた身体状態の測定は行われないものとする。)。この例では、血圧計400からの血圧の時系列データが、病院端末200Aによって受信され、メモリ220に一旦記憶される。続いて、例えばユーザとしての医師が病院端末200Aの操作部230を操作して、被験者90について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク900を介してサーバ300に送信する。
ii) サーバ300は、この例では常時動作状態にあり(ただし、メンテナンス期間等を除く。)、図5A(基本的な血圧関連情報表示方法を実行する際のサーバ300の動作フローを示す。)のステップS101に示すように、病院端末200Aからのデータを待っている。サーバ300の制御部310は、ネットワーク900から入力部としての通信部390を介して、病院端末200Aからの上述の血圧の時系列データを受信すると(ステップS101でYES)、受信した血圧の時系列データを記憶部320のデータベース321に記憶する。そして、次に述べるステップS102~S107の処理を実行する。
iii) まず、ステップS102に示すように、制御部310は前処理部として働いて、血圧の時系列データに公知の移動平均などなどを用いた平滑化、ノイズ除去、ローパスフィルタを用いた高周波成分除去等の前処理を施す。
iv) 次に、ステップS103に示すように、制御部310は血圧サージ検出部として働いて、被験者90の血圧の時系列データから、例えば特願2017-048946号、特願2017-050066号に開示されているように、予め定められた判定基準に基づいて血圧サージを検出する。これにより、例えば図6(B)中に、破線の矩形枠803,803,…で示すように、複数の血圧サージを検出する。血圧サージは一晩に数百個も発生することがあると言われている。この例では、収縮期血圧(SBP)に相当するピークに黒丸●が付され、それらが包絡線で結ばれている。また、拡張期血圧(DBP)に相当するピークに白丸○が付され、それらが包絡線で結ばれている。
ここで、血圧サージを検出するための「予め定められた判定基準」とは、例えば図7(血圧サージの個別波形の一例を曲線Cで示す。)に示すように、ピーク検出区間(例えば、15拍分の期間)内にサージ開始点P1の時刻からサージピーク点P2の時刻までが含まれるとともに、サージ開始点P1の時刻における収縮期血圧値(SBP)とピーク点P2の時刻における収縮期血圧値(SBP)との差(血圧変動量)L1が20mmHg(または15mmHg)以上であること、サージ開始点P1の時刻とピーク点P2の時刻との間の期間T1が5拍分よりも大きいこと、かつ、ピーク点P2の時刻とサージ終了点P4の時刻との間の期間T3が7拍分よりも大きいことを指す。この例では、サージ開始点P1は、ピーク点P2よりも前で、収縮期血圧値(SBP)の極小を与える点として規定されている。また、サージ終了点P4は、ピーク点P2よりも後で、ピーク点P2からL1×0.75(=L3)の分だけ血圧が降下した点として規定されている。
v)次に、図5AのステップS104に示すように、制御部310は個別波形取得部として働いて、検出された各血圧サージ803,803,…について、その血圧サージ803をなす複数の拍対応ピーク(この例では、収縮期血圧(SBP)に相当するピーク)を結ぶ包絡線を個別波形として求める。図7に示したように、血圧サージの個別波形は、山状の形状をなす曲線Cとして表される。
vi)次に、図5AのステップS105に示すように、制御部310は統計処理部として働いて、この例では、得られた全ての個別波形を統計処理して、上記時系列データ801内での血圧サージの代表的波形と波形ばらつきを求める。ここで、「統計処理」とは、個別波形を平均する処理を指す。この例では、血圧サージの「代表的波形」とは、複数の個別波形を平均して得られた平均的波形を指す。また、この例では、血圧サージの「波形ばらつき」とは、複数の個別波形の分布の幅を指す。
概念的に説明すると、血圧サージの個別波形を平均する処理とは、次のような処理を意味する。図8中に例示するように、血圧サージとして「サージNo.1」の個別波形C1、「サージNo.2」の個別波形C2が存在するものとする。この場合、血圧サージの個別波形を平均する処理とは、時間経過を表す横方向座標Xと、血圧サージに伴う血圧変動量を表す縦方向座標Yとを含む座標平面Q上で、横方向Xに関して複数の個別波形C1,C2を相対的にスライドさせて複数の個別波形C1,C2のピークの位置を揃えた状態で、横方向座標X毎に複数の個別波形C1,C2の血圧変動量データを平均して、代表的波形としての平均的波形を示す曲線Cavを求める処理に相当する。ここで、図8では、時間経過を表す横方向座標Xは、拍対応ピークを特定する拍番号によって定められている。縦方向座標Yは血圧サージに伴う収縮期血圧値(SBP)の血圧変動量によって定められている。なお、具体的なデータ処理(血圧サージの個別波形の統計処理)については、後に詳述する。
また、この例では、血圧サージの「波形ばらつき」の範囲は、複数の個別波形の血圧変動量データがなす分布(横方向座標X毎の分布)の標準偏差σを用いて、横方向座標X毎に標準偏差σの±k倍(ただし、kは自然数である。)の範囲として定義される。典型的には、kは1、2または3として設定される。以下では、特段断らない限り、血圧サージの「波形ばらつき」の範囲は、±σの範囲とする。この波形ばらつき±σを示す領域(例えば、後述の図9(B)中に示す領域Sd)の表示は、個別波形のばらつきが正規分布として取り扱われるのに十分である場合(例えば、数十個以上の場合)に、特に有益な表示となる。
これにより、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとが得られる。
vii)次に、図5AのステップS106に示すように、制御部310は表示処理部の一部として働いて、血圧サージの波形ばらつきを示す領域(符号Sdで表す。)に、平均的波形を示す曲線Cavを重ね合わせた態様の画像データを作成する。
例えば、図9(A)に示すように、図5AのステップS104で個別波形C11,C12,C13,…が得られたものとする。この場合、図5AのステップS106によって、例えば図9(B)に示すように、波形ばらつきを示す領域Sdに、平均的波形を示す曲線Cavを重ね合わせた態様の画像データImが得られる。この例では、平均的波形を示す曲線Cavは、或る一定の濃度を有する実線で表されている。波形ばらつきを示す領域Sdは、曲線Cavと背景領域(この例では、白色領域)Bgとの間の中間濃度を有する領域として表されている。
viii)この後、図5AのステップS107に示すように、サーバ300の制御部310は、作成した画像データImを、ネットワーク900を介して、この例では、血圧の時系列データの提供元である病院端末200Aへ送信する。
なお、サーバ300の制御部310は、画像データImを、例えば病院端末200Aを操作するユーザの指定に応じて、病院端末200A以外の他の病院端末200B等へ送信してもよい。
(血圧サージの個別波形の統計処理)
図10は、図5AのステップS105で述べた、サーバ300の制御部310による血圧サージの個別波形についての統計処理の詳細なフローを示している。
図10は、図5AのステップS105で述べた、サーバ300の制御部310による血圧サージの個別波形についての統計処理の詳細なフローを示している。
図10のステップS111に示すように、制御部310は、データベース321から、検出された全ての血圧サージについて、個別波形の拍対応ピークの時系列データを取得する。この例では、取得された血圧サージの数はN個であったものとする。
例えば、図11のテーブルMT1に示すように、検出された血圧サージNo.1は、開始時刻(開始点の時刻)22:20:15、ピーク時刻(ピークの時刻)22:20:21、終了時刻(終了点の時刻)22:20:26で特定されている。また、血圧サージNo.2は、開始時刻22:25:35、ピーク時刻22:25:44、終了時刻22:25:52で特定されている。なお、時刻表記の仕方は、「時:分:秒」である(以下同様)。
この場合、例えば血圧サージNo.1について、個別波形の拍対応ピークの時系列データは、例えば図12のテーブルMT2に示すように、開始時刻22:20:15では、拍対応ピークとしてのSBPに相当するピークが130[mmHg]、DBPに相当するピークが78[mmHg]、脈拍数PRが61[拍/分]というように記録されている。次の拍に相当する時刻22:20:16では、SBPに相当するピークが134[mmHg]、DBPに相当するピークが80[mmHg]、脈拍数PRが60[拍/分]というように記録されている。このようにして順次拍対応ピークの時系列データが記録され、終了時刻22:20:26では、SBPに相当するピークが133[mmHg]、DBPに相当するピークが82[mmHg]、脈拍数PRが61[拍/分]というように記録されている。簡単のため図示を省略するが、血圧サージNo.2についても、個別波形の拍対応ピークの時系列データは、同様に記録されている。
次に、図10のステップS112に示すように、制御部310は、N個の個別波形について、それぞれその個別波形の開始点P1を基準(=0[mmHg])としてオフセットする。すなわち、N個の個別波形について、それぞれ、各時刻の血圧値(この例では、SBP値)からその個別波形の開始点P1の血圧値(この例では、SBP値)を差し引く。これにより、時刻毎の血圧変動量データを得る。このようにオフセットを行うことによって、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとを把握し易くなる。
次に、図10のステップS113に示すように、制御部310は、N個の個別波形のうちX方向(横方向)寸法が最長の個別波形を特定する。また、その最長のX方向寸法を拍数M個として求める。
次に、図10のステップS114に示すように、制御部310は、記憶部320に、演算を実行するための作業領域としてN行M列のメモリ領域を確保する。
ここで簡単のため、図8中に例示したように、取得された血圧サージの数N=2個であり、「サージNo.1」の個別波形C1のX方向寸法が拍数12個分であり、「サージNo.2」の個別波形C2のX方向寸法が拍数17個分であったとする。このとき、制御部310は、X方向寸法が最長の個別波形をC2であると特定し、その最長のX方向寸法を拍数M=17個として求める。そして、制御部310は、記憶部320に、図13のテーブルMT3の2行目~3行目に示すように、「サージNo.1」、「サージNo.2」のための演算を実行するための作業領域として2行17列のメモリ領域を確保する。なお、実際に、取得された血圧サージの数Nが2より大きい場合は、2より多い行数のメモリ領域が確保される。テーブルMT3の最上段(1行目)には、X方向座標を表す拍番号が表されている。
次に、図10のステップS115に示すように、制御部310は、N行M列のメモリ領域において横方向Xに関してN個の個別波形を相対的にスライドさせてN個の個別波形のピークの位置(拍番号)を揃える。例えば図13の例では、横方向Xに関して「サージNo.1」、「サージNo.2」の2個の個別波形C1,C2のピークの位置が、拍番号10(破線の丸印で囲んで示す)の所に揃えられている。この結果、テーブルMT3の2行目に、「サージNo.1」の血圧変動量データが、拍番号4~15にわたって拍番号毎に「0」、「4」、「10」、「13」、…、「11」、「3」というように表されている(単位はmmHg。以下同様。)。また、テーブルMT3の3行目に、「サージNo.2」の血圧変動量データが、拍番号1~17にわたって拍番号毎に「0」、「2」、「3」、「5」、…、「4」、「2」というように表されている。このように、N個の個別波形のピークの位置(拍番号)を揃えた状態は、図8中に示した座標平面Q上で、横方向Xに関して複数の個別波形C1,C2を相対的にスライドさせて複数の個別波形C1,C2のピークの位置を揃えた状態と等価になっている。これにより、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとを、さらに把握し易くなる。なお、「サージNo.1」の個別波形C1の開始点(拍番号4)の値、「サージNo.2」の個別波形C2の開始点(拍番号1)の値がいずれも「0」になっているのは、既述のオフセット(図10のステップS112)の結果である。
次に、図10のステップS116に示すように、制御部310は、N個の個別波形の血圧変動量データを、拍番号(列)毎に統計処理(この例では、平均化)して、平均値を求める。例えば図13の例では、テーブルMT3の最下段(4行目)に、上述の「サージNo.1」、「サージNo.2」の個別波形の場合の血圧変動量データについての拍番号毎の平均値が表されている。この拍番号の平均値によって、代表的波形としての平均的波形を示す曲線Cavが定められる。
このとき、制御部310は、N個の個別波形のうちX方向寸法が最長の個別波形(上の例では、「サージNo.2」)よりも短い個別波形(上の例では、「サージNo.1」)について、最長の個別波形(上の例では、「サージNo.2」)よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする。例えば図13の例では、「サージNo.1」の個別波形について、X方向寸法が「サージNo.2」の個別波形よりも不足する部分、すなわち、拍番号1~3の部分、および、拍番号16~17の部分について、上記血圧変動量データの統計処理に対する寄与をゼロにする。これにより、座標平面Q上でN個の個別波形のX方向寸法が互いに異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。
さらに、図10のステップS117に示すように、制御部310は、N個の個別波形の血圧変動量データがなす分布の標準偏差σを、拍番号(列)毎に求める。拍番号毎に±σの領域が、既述の波形ばらつきを示す領域Sdを定める。
このように、このフローでは、N個の個別波形の血圧変動量データを統計処理する際に、拍番号(列)毎に血圧変動量データの統計処理を行えば足りる。したがって、横方向Xに関して連続的に血圧変動量データの統計処理を行う場合に比して、制御部310の演算量が少なくて済む。
なお、上述の図10のステップS112で、N個の個別波形について、それぞれその個別波形の開始点P1を基準(=0[mmHg])としてオフセットするのに加えて、さらに、N個の個別波形のピークの高さを揃えるように正規化してもよい。これにより、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとを、さらに把握し易くなる。
(病院端末での画像データの表示)
病院端末200Aは、図9(B)に示したような画像データImを、ネットワーク900を介して、サーバ300から受信する。
次に、病院端末200Aの操作部230を介したユーザの表示指示があると、病院端末200Aの制御部210は、表示処理部の一部として働いて、表示器240の表示画面上に、その画像データImを表示する。このとき、表示器240の表示画面上に、血圧変動量(図7中に示したサージ開始点P1の時刻における収縮期血圧値(SBP)とピーク点P2の時刻における収縮期血圧値(SBP)との差L1)の平均値、分布の標準偏差σ(または±σ,±2σ,±3σ)の値などを併せて表示してもよい。
ユーザは、表示器240の表示画面を見ることによって、被験者90における血圧サージの代表的波形を示す曲線Cavと波形ばらつきを示す領域Sdを把握することができる。このことは、SASの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。
なお、ユーザは、病院端末200Aの操作部230を介して転送指示を入力して、その画像データImを、病院端末200A以外の他の病院端末200B等へ送信してもよい。
[B.応用的な血圧関連情報表示方法]
この応用的な血圧関連情報表示方法では、血圧計400が出力する測定時刻(時間)と血圧とを対応付けた血圧の時系列データと、PSG装置500が出力する身体状態特定期間を表す情報とを用いて、次のように画像データの作成・表示を行う。
この応用的な血圧関連情報表示方法では、血圧計400が出力する測定時刻(時間)と血圧とを対応付けた血圧の時系列データと、PSG装置500が出力する身体状態特定期間を表す情報とを用いて、次のように画像データの作成・表示を行う。
(サーバ300での画像データの作成)
i) 図4に示した態様で、被験者90について、一晩(睡眠期間を含む。)にわたって、血圧計400を用いた連続的な血圧測定と、PSG装置500を用いた身体状態の測定が行われたものとする。この例では、血圧計400からの血圧の時系列データと、PSG装置500からの身体状態特定期間を表す情報とが、病院端末200Aによって受信され、メモリ220に一旦記憶される。続いて、例えばユーザとしての医師が病院端末200Aの操作部230を操作して、被験者90について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク900を介してサーバ300に送信する。
i) 図4に示した態様で、被験者90について、一晩(睡眠期間を含む。)にわたって、血圧計400を用いた連続的な血圧測定と、PSG装置500を用いた身体状態の測定が行われたものとする。この例では、血圧計400からの血圧の時系列データと、PSG装置500からの身体状態特定期間を表す情報とが、病院端末200Aによって受信され、メモリ220に一旦記憶される。続いて、例えばユーザとしての医師が病院端末200Aの操作部230を操作して、被験者90について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク900を介してサーバ300に送信する。
ii) サーバ300は、この例では常時動作状態にあり(ただし、メンテナンス期間等を除く。)、図5B(応用的な血圧関連情報表示方法を実行する際のサーバ300の動作フローを示す。)のステップS151に示すように、病院端末200Aからのデータを待っている。サーバ300の制御部310は、ネットワーク900から入力部としての通信部390を介して、病院端末200Aからの上述のデータ(血圧の時系列データと身体状態特定期間を表す情報)を受信すると(ステップS151でYES)、受信した血圧の時系列データと身体状態特定期間を表す情報を記憶部320のデータベース321に記憶する。そして、次に述べるステップS152~S157の処理を実行する。
iii) まず、ステップS152に示すように、制御部310は前処理部として働いて、血圧の時系列データに公知の移動平均などなどを用いた平滑化、ノイズ除去、ローパスフィルタを用いた高周波成分除去等の前処理を施す。
iv) 次に、ステップS153に示すように、制御部310は血圧サージ検出部として働いて、被験者90の血圧の時系列データから、例えば特願2017-048946号、特願2017-050066号に開示されているように、予め定められた判定基準に基づいて血圧サージを検出する。これにより、例えば図6(B)中に、破線の矩形枠803,803,…で示すように、複数の血圧サージを検出する。ここで、血圧サージを検出するための「予め定められた判定基準」は、図5AのステップS103に関して図7を用いて説明した判定基準と同様である。
v)次に、図5BのステップS154に示すように、制御部310は個別波形取得部として働いて、検出された各血圧サージ803,803,…について、その血圧サージ803をなす複数の拍対応ピーク(この例では、収縮期血圧(SBP)に相当するピーク)を結ぶ包絡線を個別波形として求める。図7に示したように、血圧サージの個別波形は、山状の形状をなす曲線Cとして表される。
vi)次に、図5BのステップS155に示すように、制御部310は統計処理部として働いて、この例では既述の身体状態特定期間毎に、それぞれ複数得られた個別波形を統計処理して、上記時系列データ801内での血圧サージの代表的波形と波形ばらつきを求める。この例では、或る身体状態特定期間(例えば、無呼吸期間)が血圧サージの開始点P1から終了点P4までの期間と一部でも重なっていれば、その血圧サージはその身体状態特定期間に発生したものとして処理される。ここで、「統計処理」とは、図5AのステップS105に関して図8(および図10~図13)を用いて説明したのと同様に、個別波形を平均する処理を指す。すなわち、血圧サージの個別波形を平均する処理は、時間経過を表す横方向座標Xと、血圧サージに伴う血圧変動量を表す縦方向座標Yとを含む座標平面Q上で、横方向Xに関して複数の個別波形C1,C2を相対的にスライドさせて複数の個別波形C1,C2のピークの位置を揃えた状態で、横方向座標X毎に複数の個別波形C1,C2の血圧変動量データを平均して、代表的波形としての平均的波形を示す曲線Cavを求める処理に相当する。また、この例では、図5AのステップS105に関して説明したのと同様に、血圧サージの「波形ばらつき」の範囲は、複数の個別波形の血圧変動量データがなす分布(横方向座標X毎の分布)の標準偏差σを用いて、横方向座標X毎に±σの範囲として定義される。
これにより、身体状態特定期間毎に、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとが得られる。例えば、一晩の時系列データで特定された身体状態特定期間が無呼吸期間(低呼吸期間を含む。)、レム睡眠期間、ノンレム睡眠期間、覚醒期間の4種類である場合、無呼吸期間について、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとが得られる。同様に、レム睡眠期間、ノンレム睡眠期間、覚醒期間のそれぞれについて、代表的波形としての平均的波形を示す曲線Cavと、横方向座標X毎の波形ばらつき±σとが得られる。
vii)次に、図5BのステップS156に示すように、制御部310は表示処理部の一部として働いて、既述の身体状態特定期間毎に、図9(B)に例示したような、血圧サージの波形ばらつきを示す領域Sdに、平均的波形を示す曲線Cavを重ね合わせた態様の画像データImを作成する。
例えば、後述の図15A中に示すように、無呼吸期間(低呼吸期間を含む。)、レム睡眠期間(REM)、ノンレム睡眠期間(NREM)、覚醒期間のそれぞれについて、波形ばらつきを示す領域Sd1,Sd2,Sd3,Sd4に、それぞれ平均的波形を示す曲線Cav1,Cav2,Cav3,Cav4を重ね合わせた態様の画像データIm1,Im2,Im3,Im4が得られる。
viii)この後、図5BのステップS157に示すように、サーバ300の制御部310は、身体状態特定期間毎に作成した画像データIm1,Im2,Im3,Im4を、ネットワーク900を介して、この例では、血圧の時系列データの提供元である病院端末200Aへ送信する。
なお、サーバ300の制御部310は、画像データIm1,Im2,Im3,Im4を、例えば病院端末200Aを操作するユーザの指定に応じて、病院端末200A以外の他の病院端末200B等へ送信してもよい。
(病院端末での画像データの表示)
図14は、サーバ300から身体状態特定期間毎に作成された画像データを受信した場合の、病院端末200Aでの表示処理のフローを示している。なお、この例では、病院端末200Aについて説明するが、他の病院端末200B等も同様の表示処理が可能になっている。
図14は、サーバ300から身体状態特定期間毎に作成された画像データを受信した場合の、病院端末200Aでの表示処理のフローを示している。なお、この例では、病院端末200Aについて説明するが、他の病院端末200B等も同様の表示処理が可能になっている。
病院端末200Aは、まず図14のステップS201に示すように、身体状態特定期間毎に作成された画像データを、ネットワーク900を介して、サーバ300から受信する。この例では、身体状態特定期間は、無呼吸期間(低呼吸期間を含む。)、レム睡眠期間、ノンレム睡眠期間、覚醒期間の4種類であるものとする。
次に、病院端末200Aの制御部210は、表示処理部の一部として働いて、図15Aに示すように、表示器240の表示画面上に、身体状態特定期間毎に作成された画像データIm1,Im2,Im3,Im4を並べてサムネイル表示する(図14のステップS202)。図15A中に示すように、この例では、表示器240の表示画面には、身体状態特定期間を特定する身体状態特定欄241と、その下方に対応して配置された画像データ欄242とが表示されている。この例では、身体状態特定欄241には、左から右へ向かって、無呼吸期間を表す「無呼吸」欄241aと、レム睡眠期間を表す「REM」欄241bと、ノンレム睡眠期間を表す「NREM」欄241cと、覚醒期間を表す「覚醒反応」欄241dとが含まれている。画像データ欄242には、無呼吸期間について作成された画像データIm1を表示する欄242aと、レム睡眠期間について作成された画像データIm2を表示する欄242bと、ノンレム睡眠期間について作成された画像データIm3を表示する欄242cと、覚醒期間について作成された画像データIm4を表示する欄242dとが含まれている。
この例では、無呼吸期間について作成された画像データIm1は、波形ばらつきを示す領域Sd1に、平均的波形を示す曲線Cav1を重ね合わせた態様になっている。レム睡眠期間について作成された画像データIm2は、波形ばらつきを示す領域Sd2に、平均的波形を示す曲線Cav2を重ね合わせた態様になっている。ノンレム睡眠期間について作成された画像データIm3は、波形ばらつきを示す領域Sd3に、平均的波形を示す曲線Cav3を重ね合わせた態様になっている。また、覚醒期間について作成された画像データIm4は、波形ばらつきを示す領域Sd4に、平均的波形を示す曲線Cav4を重ね合わせた態様になっている。
ユーザは、これらの画像データIm1,Im2,Im3,Im4のサムネイル表示を見ることによって、身体状態特定期間毎に、言い換えれば、被験者90が血圧サージの要因となり得る身体状態をとっている期間毎に、被験者90における血圧サージの代表的波形を示す曲線Cavと波形ばらつきを示す領域Sdを把握することができる。したがって、ユーザは、例えば最も深刻な血圧サージを引き起こしている要因(身体状態)を容易に把握することができる。例えば、画像データIm1,Im2,Im3,Im4のうち、無呼吸期間について作成された画像データIm1の血圧変動量が最も大きいことから、「無呼吸」が最も深刻な血圧サージを引き起こしていると判断してもよい。
次に、病院端末200Aの制御部210は、図14のステップS203に示すように、表示されている画像データIm1,Im2,Im3,Im4が選択されたか否かを判断する。ユーザが操作部230を操作することによって画像データIm1,Im2,Im3,Im4のうちいずれかを選択した場合(図14のステップS203でYES)、制御部210は、図15Bに示すように、表示器240の表示画面上に、その選択された身体状態特定期間についての画像データ(この例では、無呼吸期間について作成された画像データIm1)を拡大表示する(図14のステップS204)。このとき、表示器240の表示画面上に、その選択された身体状態特定期間についての血圧変動量(図7中に示したサージ開始点P1の時刻における収縮期血圧値(SBP)とピーク点P2の時刻における収縮期血圧値(SBP)との差L1)の平均値、分布の標準偏差σ(または±σ,±2σ,±3σ)の値などを併せて表示してもよい。
なお、画像データIm1,Im2,Im3,Im4がいずれも選択されない場合(図14のステップS203でNO)、制御部210は、画像をオフして(図14のステップS205)、この表示処理を終了する。
なお、ユーザは、病院端末200Aの操作部230を介して転送指示を入力して、画像データIm1,Im2,Im3,Im4を、病院端末200A以外の他の病院端末200B等へ送信してもよい。
上述のように、血圧サージの波形ばらつきを示す領域に平均的波形を示す曲線を重ね合わせた態様の画像データIm1,Im2,Im3,Im4を身体状態特定期間毎に表示することは、SASの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。
(変形例1)
上述の実施形態では、サーバ300が画像データIm(または、身体状態特定期間毎の画像データIm1,Im2,Im3,Im4。以下同様。)を作成し、病院端末200A,200B,…の表示器240が画像を表示したが、これに限られるものではない。サーバ300の制御部310は、画像データImを病院端末200A,200B,…へ送信する代わりに、画像作成のためのデータのみを送信し、病院端末200A,200B,…の制御部210が専ら画像データImを作成しても良い。
上述の実施形態では、サーバ300が画像データIm(または、身体状態特定期間毎の画像データIm1,Im2,Im3,Im4。以下同様。)を作成し、病院端末200A,200B,…の表示器240が画像を表示したが、これに限られるものではない。サーバ300の制御部310は、画像データImを病院端末200A,200B,…へ送信する代わりに、画像作成のためのデータのみを送信し、病院端末200A,200B,…の制御部210が専ら画像データImを作成しても良い。
(変形例2)
また、上述の実施形態では、本発明の血圧関連情報表示装置を、病院端末200A,200B,…とサーバ300とを含むネットワーク上のシステム100として構成したが、これに限られるものではない。
また、上述の実施形態では、本発明の血圧関連情報表示装置を、病院端末200A,200B,…とサーバ300とを含むネットワーク上のシステム100として構成したが、これに限られるものではない。
例えば、本発明の血圧関連情報表示装置を、病院端末200A,200B,…のいずれかのみで構成しても良い。つまり、病院端末(例えば200A)に上述の血圧関連情報表示方法の全部(血圧計400からの血圧の時系列データおよびPSG装置500からの身体状態特定期間を表す情報の受信から、表示器240の表示画面上への画像データIm等の表示までを含む。)を実行させてもよい。
その場合、病院端末200Aのメモリ220には、制御部210に、上述の血圧関連情報表示方法を実行させるプログラムをインストールしておく。これにより、本発明の血圧関連情報表示装置を小型かつコンパクトに構成することができる。
また、上述の血圧関連情報表示方法を、ソフトウェア(コンピュータプログラム)として、CD(コンパクトディスク)、DVD(デジタル万能ディスク)、フラッシュメモリなどの非一時的(non-transitory)にデータを記憶可能な記録媒体に記録してもよい。このような記録媒体に記録されたソフトウェアを、パーソナルコンピュータ、PDA(パーソナル・デジタル・アシスタンツ)、スマートフォンなどの実質的なコンピュータ装置にインストールすることによって、それらのコンピュータ装置に、上述の血圧関連情報表示方法を実行させることができる。
(変形例3)
上述の実施形態では、画像データImを構成する波形ばらつきを示す領域Sdは、横方向座標X毎に血圧変動量データがなす分布の標準偏差σの±k倍(典型的には、±σ,±2σまたは±3σ)の範囲として定義された。しかしながら、波形ばらつきを示す領域Sdの範囲は、標準偏差σを用いて定義されたものに限られない。
上述の実施形態では、画像データImを構成する波形ばらつきを示す領域Sdは、横方向座標X毎に血圧変動量データがなす分布の標準偏差σの±k倍(典型的には、±σ,±2σまたは±3σ)の範囲として定義された。しかしながら、波形ばらつきを示す領域Sdの範囲は、標準偏差σを用いて定義されたものに限られない。
例えば、波形ばらつきを示す領域Sdは、横方向座標X毎に血圧変動量データがなす分布の四分位範囲として定義されてもよい。すなわち、処理としては、横方向座標X毎に血圧変動量データを昇順に並べる。そして、第一四分位を領域Sdの下限とし、第三四分位を領域Sdの上限とする。これにより、波形ばらつきを示す領域Sdを求めることができる。
このように、波形ばらつきを示す領域Sdを、横方向座標X毎に血圧変動量データがなす分布の四分位範囲とした場合、標準偏差σを用いて定義する場合に比して、多くの計算を必要とせず、計算が簡単で済む。したがって、処理速度の向上やメモリの節約に貢献することができる。
また、血圧サージの「代表的波形」は、複数の個別波形を平均して得られた平均的波形としたが、これに限られるものではない。血圧サージの「代表的波形」は、横方向座標X毎に血圧変動量データがなす分布の第二四分位(中央値)であってもよい。この場合、処理としては、横方向座標X毎に血圧変動量データを昇順に並べ、横方向に中央値を連ねて、代表的波形を得ることかできる。
(変形例4)
上述の実施形態では、画像データImを構成する波形ばらつきを示す領域Sdは、表示器240の表示画面上で、例えば図9(B)に示したように、代表的波形を示す曲線Cavと背景領域Bgとの間の中間濃度を有する領域として表示された。しかしながら、波形ばらつきを示す領域Sdの表示態様は、中間濃度を有する領域に限られない。
上述の実施形態では、画像データImを構成する波形ばらつきを示す領域Sdは、表示器240の表示画面上で、例えば図9(B)に示したように、代表的波形を示す曲線Cavと背景領域Bgとの間の中間濃度を有する領域として表示された。しかしながら、波形ばらつきを示す領域Sdの表示態様は、中間濃度を有する領域に限られない。
例えば、図16(A)は、図9(A)と同様に、複数の個別波形C11,C12,C13,…を例示している。この場合において、図16(B)に示すように、波形ばらつきを示す領域Sd′をなすように、代表的波形を示す曲線Cav′とは区別可能な態様で、複数の個別波形C11,C12,C13,…を示す曲線Ciを表示してもよい。この例では、曲線Cav′が実線で表示されているのに対して、複数の個別波形C11,C12,C13,…を示す曲線Ciはそれぞれ破線で表示されている。この例では、この波形ばらつきを示す領域Sd′(複数の曲線Ciを含む)に代表的波形を示す曲線Cav′を重ね合わせた態様の画像データが符号Im′で表されている。
これにより、ユーザは、直感的に、被験者90における血圧サージの代表的波形を示す曲線Cav′と波形ばらつきを示す領域Sd′を把握できる。この複数の個別波形を示す曲線Ciの表示は、例えば個別波形の数が数個以下の場合に、ユーザに個別波形C11,C12,C13,…のばらつきを把握させるのに有益となる。
なお、代表的波形を示す曲線Cav′と複数の個別波形C11,C12,C13,…を示す曲線Ciとは区別可能であれば良い。例えば、代表的波形を示す曲線Cav′が赤色の実線で表示される一方、複数の個別波形C11,C12,C13,…を示す曲線Ciが青色の実線で表示されてもよい。
また、身体状態特定期間毎に取得された血圧サージの数Nに対して閾値Nα(例えば、Nα=10個とする。)を設けて、取得された血圧サージの数Nが閾値Nα以上であれば、図9(B)に示したような波形ばらつきを示す領域Sd(縦方向に±σの幅の中間濃度を有する領域)を含む第1の画像データImの作成・表示を行う一方、取得された血圧サージの数Nが閾値Nα未満であれば、図16(B)に示したような波形ばらつきを示す領域Sd′(複数の個別波形C11,C12,C13,…を示す曲線Ci)を含む第2の画像データIm′の作成・表示を行ってもよい。この血圧サージの数Nが閾値Nαよりも大きいか否かの判定は、例えば、サーバ300の制御部310が血圧サージ数判定部として働いて行うことができる。
または、それに代えて、例えばサーバ300の制御部310が、図9(B)に示したような波形ばらつきを示す領域Sdを含む第1の画像データImと、図16(B)に示したような波形ばらつきを示す領域Sd′を含む第2の画像データIm′との両方を作成してもよい。この場合、例えば、サーバ300から病院端末(例えば200A)へ第1の画像データと第2の画像データとの両方を送信し、病院端末(例えば200A)の操作部230を介したユーザの入力に応じて、表示器240の表示画面上で、それらの第1の画像データImの表示と第2の画像データIm′の表示とを切り替えるようにしてもよい。
(変形例5)
また、上述の実施形態では、血圧計400はトノメトリ方式の血圧計であるものとしたが、これに限られるものではない。血圧計400は、被測定部位のうち対応する部分を通る動脈へ向けて光を照射する発光素子と、その光の反射光(または透過光)を受光する受光素子とを備えて、動脈の脈波を容積の変化に基づいて連続的に血圧を検出してもよい(光電方式)。また、血圧計400は、被測定部位に当接された圧電センサを備えて、被測定部位のうち対応する部分を通る動脈の圧力による歪みを電気抵抗の変化として検出し、この電気抵抗の変化に基づいて連続的に血圧を検出してもよい(圧電方式)。さらに、血圧計400は、被測定部位のうち対応する部分を通る動脈へ向けて電波(送信波)を送る送信素子と、その電波の反射波を受信する受信素子とを備えて、動脈の脈波による動脈とセンサとの間の距離の変化を送信波と反射波との間の位相のずれとして検出し、この位相のずれに基づいて連続的に血圧を検出してもよい(電波照射方式)。また、血圧を算出することができる物理量を観測することができれば、これらの以外の方式を適用してもよい。
また、上述の実施形態では、血圧計400はトノメトリ方式の血圧計であるものとしたが、これに限られるものではない。血圧計400は、被測定部位のうち対応する部分を通る動脈へ向けて光を照射する発光素子と、その光の反射光(または透過光)を受光する受光素子とを備えて、動脈の脈波を容積の変化に基づいて連続的に血圧を検出してもよい(光電方式)。また、血圧計400は、被測定部位に当接された圧電センサを備えて、被測定部位のうち対応する部分を通る動脈の圧力による歪みを電気抵抗の変化として検出し、この電気抵抗の変化に基づいて連続的に血圧を検出してもよい(圧電方式)。さらに、血圧計400は、被測定部位のうち対応する部分を通る動脈へ向けて電波(送信波)を送る送信素子と、その電波の反射波を受信する受信素子とを備えて、動脈の脈波による動脈とセンサとの間の距離の変化を送信波と反射波との間の位相のずれとして検出し、この位相のずれに基づいて連続的に血圧を検出してもよい(電波照射方式)。また、血圧を算出することができる物理量を観測することができれば、これらの以外の方式を適用してもよい。
以上の実施形態は例示であり、この発明の範囲から離れることなく様々な変形が可能である。上述した複数の実施の形態は、それぞれ単独で成立し得るものであるが、実施の形態同士の組みあわせも可能である。また、異なる実施の形態の中の種々の特徴も、それぞれ単独で成立し得るものであるが、異なる実施の形態の中の特徴同士の組みあわせも可能である。
100 システム
200A,200B,… 病院端末
240 表示器
300 サーバ
400 血圧計
500 PSG装置
200A,200B,… 病院端末
240 表示器
300 サーバ
400 血圧計
500 PSG装置
Claims (10)
- 複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示装置であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する血圧サージ検出部と、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める個別波形取得部と、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める統計処理部と、
表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う表示処理部とを備え、
上記統計処理部は、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
上記統計処理部は、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする
ことを特徴とする血圧関連情報表示装置。 - 請求項1に記載の血圧関連情報表示装置において、
上記複数の個別波形の血圧変動量データは、それぞれ、その個別波形の開始点の血圧値を基準とした変動量であることを特徴とする血圧関連情報表示装置。 - 請求項1または2に記載の血圧関連情報表示装置において、
上記複数の個別波形の血圧変動量データは、上記複数の個別波形のピークの高さを揃えるように正規化されていることを特徴とする血圧関連情報表示装置。 - 請求項1から3までのいずれか一つに記載の血圧関連情報表示装置において、
上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められ、
上記統計処理部は、上記拍番号毎に上記複数の個別波形の血圧変動量データを統計処理することを特徴とする血圧関連情報表示装置。 - 請求項1から4までのいずれか一つに記載の血圧関連情報表示装置において、
上記統計処理部が求める上記波形ばらつきは、kを自然数としたとき、上記複数の個別波形の血圧変動量データがなす分布の標準偏差の±k倍として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±k倍に相当する幅をもつ領域を表示することを特徴とする血圧関連情報表示装置。 - 請求項1から4までのいずれか一つに記載の血圧関連情報表示装置において、
上記統計処理部が求める上記波形ばらつきは、上記複数の個別波形の血圧変動量データがなす分布の四分位範囲として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域を表示することを特徴とする血圧関連情報表示装置。 - 請求項1から4までのいずれか一つに記載の血圧関連情報表示装置において、
上記表示処理部は、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線とは区別可能な態様で、上記複数の個別波形を示す曲線を表示する処理を行うことを特徴とする血圧関連情報表示装置。 - 請求項1から7までのいずれか一つに記載の血圧関連情報表示装置において、
上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態が特定された身体状態特定期間を表す情報を入力する入力部を備え、
上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求め、
上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うことを特徴とする血圧関連情報表示装置。 - 複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示方法であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出するステップと、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求めるステップと、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求めるステップとを備え、
このステップでは、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
このステップでは、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにし、
さらに、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うステップを備える
ことを特徴とする血圧関連情報表示方法。 - 請求項9に記載の血圧関連情報表示方法をコンピュータに実行させるためのプログラム。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019057000A JP7267055B2 (ja) | 2019-03-25 | 2019-03-25 | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム |
PCT/JP2020/008976 WO2020195615A1 (ja) | 2019-03-25 | 2020-03-03 | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム |
CN202080024289.4A CN113645897A (zh) | 2019-03-25 | 2020-03-03 | 血压关联信息显示装置、血压关联信息显示方法、以及程序 |
DE112020001440.1T DE112020001440T5 (de) | 2019-03-25 | 2020-03-03 | Anzeigevorrichtung für blutdruckbezogene Informationen, Anzeigeverfahren für blutdruckbezogene Informationen und Programm |
US17/484,422 US20220007951A1 (en) | 2019-03-25 | 2021-09-24 | Blood pressure-related information display device, blood pressure-related information display method, and computer-readable recording medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019057000A JP7267055B2 (ja) | 2019-03-25 | 2019-03-25 | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2020156609A JP2020156609A (ja) | 2020-10-01 |
JP2020156609A5 true JP2020156609A5 (ja) | 2022-03-23 |
JP7267055B2 JP7267055B2 (ja) | 2023-05-01 |
Family
ID=72610007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019057000A Active JP7267055B2 (ja) | 2019-03-25 | 2019-03-25 | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220007951A1 (ja) |
JP (1) | JP7267055B2 (ja) |
CN (1) | CN113645897A (ja) |
DE (1) | DE112020001440T5 (ja) |
WO (1) | WO2020195615A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6766708B2 (ja) * | 2017-03-15 | 2020-10-14 | オムロンヘルスケア株式会社 | 情報処理装置、方法及びプログラム |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7527597B2 (en) * | 2001-01-16 | 2009-05-05 | Biomedical Acoustic Research Corporation | Acoustic detection of vascular conditions |
JP3968699B2 (ja) * | 2002-03-29 | 2007-08-29 | 有限会社 ソフトロックス | 生体器官のデータ監視装置及びその監視方法 |
WO2009144723A1 (en) * | 2008-05-28 | 2009-12-03 | Ofer Hornick | Method and apparatus for co2 evaluation |
JP6111741B2 (ja) * | 2013-03-04 | 2017-04-12 | オムロンヘルスケア株式会社 | 電子血圧計 |
CN104398257A (zh) * | 2014-12-08 | 2015-03-11 | 蚌埠医学院 | 一种心电波形周期性统计分析方法 |
JP6508473B2 (ja) | 2015-08-31 | 2019-05-08 | 東芝ライテック株式会社 | 照明器具 |
JP6548522B2 (ja) | 2015-08-31 | 2019-07-24 | カルソニックカンセイ株式会社 | 吸着器の作製方法 |
JP6547838B2 (ja) | 2015-11-10 | 2019-07-24 | オムロンヘルスケア株式会社 | 診断支援装置、診断支援方法、診断支援プログラム |
JP6766708B2 (ja) * | 2017-03-15 | 2020-10-14 | オムロンヘルスケア株式会社 | 情報処理装置、方法及びプログラム |
-
2019
- 2019-03-25 JP JP2019057000A patent/JP7267055B2/ja active Active
-
2020
- 2020-03-03 DE DE112020001440.1T patent/DE112020001440T5/de active Pending
- 2020-03-03 WO PCT/JP2020/008976 patent/WO2020195615A1/ja active Application Filing
- 2020-03-03 CN CN202080024289.4A patent/CN113645897A/zh active Pending
-
2021
- 2021-09-24 US US17/484,422 patent/US20220007951A1/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5584413B2 (ja) | 患者監視システム及び監視する方法 | |
EP2898823B1 (en) | Monitoring apparatus | |
JP2004195204A (ja) | 脈波測定装置 | |
US20220079470A1 (en) | Pattern recognition system for quantifying the likelihood of the contribution of multiple possible forms of chronic disease to patient reported dyspnea | |
JP6162977B2 (ja) | 生体情報監視装置 | |
US20190388037A1 (en) | Information processing apparatus, method, and program | |
CN110051336A (zh) | 处理生理数据的方法、装置和存储介质 | |
JP2016112049A (ja) | 脈波解析装置 | |
WO2020195615A1 (ja) | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム | |
JP6547838B2 (ja) | 診断支援装置、診断支援方法、診断支援プログラム | |
CN105852863B (zh) | 一种呼吸率测量方法及装置 | |
JP2020156609A5 (ja) | ||
JP6762061B1 (ja) | ソフトウェア | |
CN106572813B (zh) | 用于量化慢性疾病的多种可能形式对被报告呼吸困难的患者的贡献的可能性的模式识别系统 | |
US20150221111A1 (en) | System and method for displaying physiological parameters with a sparkline | |
JP5296392B2 (ja) | 自己組織化マップを用いた健康状態判定支援システム及び自己組織化マップ生成システム | |
JP7256049B2 (ja) | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム | |
WO2020170826A1 (ja) | 血圧関連情報処理装置、血圧関連情報処理方法、およびプログラム | |
JP7253419B2 (ja) | 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム | |
US20220039689A1 (en) | Gas exchange systems and methods for calculating ventilatory threshold and evaluating pulmonary arterial hyertension | |
Bickel et al. | The Use of Pneumoperitoneum During Laparoscopic Surgery as a Model to Study Pathophysiologic Phenomena: The Correlation of Cardiac Functionality with Computerized Acoustic Indices—Preliminary Data | |
EP4379734A1 (en) | System and method for patient monitoring | |
JP2022167320A (ja) | 脈波解析装置、脈波解析方法、および脈波解析プログラム | |
WO2024115345A1 (en) | System and method for patient monitoring | |
JP2020156608A5 (ja) |