JP7267055B2 - 血圧関連情報表示装置、血圧関連情報表示方法、およびプログラム - Google Patents

Description

この発明は血圧関連情報表示装置および血圧関連情報表示方法に関し、より詳しくは、被験者の血圧サージに関連する情報を可視化して表示する装置および方法に関する。また、この発明は、そのような血圧関連情報表示方法をコンピュータに実行させるためのプログラムに関する。

睡眠時無呼吸症候群（ＳＡＳ；Sleep Apnea Syndrome）を罹患している患者においては、無呼吸後の呼吸再開時に血圧が急激に上昇し、その後に下降することが知られている。本明細書では、このような急激な血圧変動を「血圧サージ」（または単に「サージ」）と呼ぶ。患者に発生した血圧サージに関連する情報（例えば、血圧サージの波形）を可視化して表示することは、ＳＡＳの診断や治療に役立つと考えられる。

従来、例えば特許文献１（ＷＯ ２０１７/０８２１０７ Ａ１）の図３では、睡眠時無呼吸症候群（ＳＡＳ）患者における変動状態（血圧サージ）の波形がグラフとして表示されている。

ＷＯ ２０１７/０８２１０７ Ａ１

しかしながら、特許文献１では、血圧サージの波形は、収縮期血圧値（ＳＢＰ）の時系列データについて時間軸を拡大して、脈拍毎のピークが分かる形態で表示されているに過ぎず、血圧の上昇、下降の傾向を示す曲線（形状）としては表示されていない。

このような状況下で、本出願人は、先に、ＳＡＳ患者おける血圧サージの波形のパタン（形状）を分類し、可視化して表示する発明を出願した（特願２０１７－０５００６６号）。これにより、医師は、比較的短い時間内に、各ＳＡＳ患者おける血圧サージの波形のパタンを把握することができる。さらに、一歩進んで、医師が、各ＳＡＳ患者における血圧サージの代表的波形と波形ばらつきを重ね合わせて見ることができれば便利であろう。また、ＳＡＳ患者に限らず、被験者の血圧サージの代表的波形と波形ばらつきを重ね合わせて見ることは、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。

そこで、この発明の課題は、血圧サージの代表的波形と波形ばらつきを重ね合わせて表示できる血圧関連情報表示装置および血圧関連情報表示方法を提供することにある。また、この発明の課題は、そのような血圧関連情報表示方法をコンピュータに実行させるためのプログラムを提供することにある。

上記課題を解決するため、この開示の血圧関連情報表示装置は、
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示装置であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する血圧サージ検出部と、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める個別波形取得部と、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める統計処理部と、
表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う表示処理部とを備え、
上記統計処理部は、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
上記統計処理部は、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする
ことを特徴とする。

本明細書で、「予め定められた判定基準」とは、典型的には、睡眠時無呼吸症候群（ＳＡＳ）患者において血圧サージを検出するための基準を指す。例えば、特願２０１７－０４８９４６号、特願２０１７－０５００６６号に開示されているような、ピーク検出区間（例えば、１５拍分の期間）内にサージ開始点の時刻からサージピーク点の時刻までが含まれるとともに、サージ開始点の時刻における収縮期血圧値とピーク点の時刻における収縮期血圧値との差（血圧変動量）が２０ｍｍＨｇ（または１５ｍｍＨｇ）以上であること、サージ開始点の時刻とピーク点の時刻との間の期間が５拍分よりも大きいこと、かつ、ピーク点の時刻とサージ終了点の時刻との間の期間が７拍分よりも大きいことを指す。

また、包絡線を作成するに用いられる、血圧サージをなす複数の「拍対応ピーク」とは、連続瞬時血圧波形の中での、収縮期血圧に相当するピークを指す。ただし、拡張期血圧値（ＤＢＰ）に相当するピークであってもよい。

また、「統計処理」とは、個別波形を平均する処理、中央値を求める処理などを指す。血圧サージの「代表的波形」とは、例えば、複数の個別波形を平均して得られた平均的波形、または、複数の個別波形の中央値に相当する波形などを指す。また、血圧サージの「波形ばらつき」とは、例えば複数の個別波形の分布の幅を指す。

また、「表示画面」は、典型的には表示器の画面を指すが、例えば、プリンタによって出力される紙面であってもよい。

また、表示画面上の例えば代表的波形を示す曲線の「表示」は、典型的には、横方向座標が時間経過（例えば、拍の進行）を表し、縦方向座標が血圧サージに伴う血圧変動量を表す態様で行われる。
本明細書で、「個別波形のピーク」とは、各血圧サージのピークを指し、各血圧サージをなす複数の拍対応ピークのうちの最大の拍対応ピークに相当する。
また、「等価な状態」とは、上記座標平面上で上記複数の個別波形のピークの位置を実際に描いて揃える処理などを要さず、メモリ領域のデータが、上記座標平面上で規定された状態と等価な状態になっていれば足りるという意味である。

この開示の血圧関連情報表示装置では、血圧サージ検出部は、被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する。個別波形取得部は、検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める。統計処理部は、複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める。表示処理部は、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う。これにより、ユーザ（典型的には、医師、看護師などの医療関係者を指す。被験者であってもよい。以下同様。）は、上記表示画面を見ることによって、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。このことは、ＳＡＳの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。

この血圧関連情報表示装置では、上記統計処理部は、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記複数の個別波形の血圧変動量データを統計処理するので、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記座標平面上で上記複数の個別波形の横方向寸法が互いに幾つかの拍数相当分異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。

一実施形態の血圧関連情報表示装置では、上記複数の個別波形の血圧変動量データは、それぞれ、その個別波形の開始点の血圧値を基準とした変動量であることを特徴とする。

この一実施形態の血圧関連情報表示装置では、上記代表的波形と上記波形ばらつきを把握し易くなる。

一実施形態の血圧関連情報表示装置では、上記複数の個別波形の血圧変動量データは、上記複数の個別波形のピークの高さを揃えるように正規化されていることを特徴とする。

この一実施形態の血圧関連情報表示装置では、上記代表的波形と上記波形ばらつきをさらに把握し易くなる。

一実施形態の血圧関連情報表示装置では、
上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められ、
上記統計処理部は、上記拍番号毎に上記複数の個別波形の血圧変動量データを統計処理することを特徴とする。

ここで、心臓の拍動によって拍間隔は一定ではないため、或る個別波形をなす血圧サージ発生時の１拍分の期間と、別の個別波形をなす血圧サージ発生時の１拍分の期間とは、互いに相違する場合がある。「上記時間経過を表す横方向座標は…拍番号によって定められ」とは、その相違にかかわらず、上記時間経過が拍の進行によって表される、という意味である。

この一実施形態の血圧関連情報表示装置では、上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められる。したがって、或る個別波形をなす血圧サージ発生時の１拍分の期間と、別の個別波形をなす血圧サージ発生時の１拍分の期間との相違にかかわらず、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記統計処理部は、上記複数の個別波形の血圧変動量データを統計処理する際に、上記拍番号毎に上記血圧変動量データの統計処理を行えば足りる。したがって、横方向に関して連続的に血圧変動量データの統計処理を行う場合に比して、上記統計処理部の演算量が少なくて済む。

一実施形態の血圧関連情報表示装置では、
上記統計処理部が求める上記波形ばらつきは、ｋを自然数としたとき、上記複数の個別波形の血圧変動量データがなす分布の標準偏差の±ｋ倍として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±ｋ倍に相当する幅をもつ領域を表示することを特徴とする。

この一実施形態の血圧関連情報表示装置では、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±ｋ倍（ただし、ｋは自然数である。典型的には、ｋは１、２または３として設定される。）に相当する幅をもつ領域が表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。この波形ばらつきを示す領域の表示は、例えば個別波形のばらつきが正規分布として取り扱われるのに十分である場合（例えば、個別波形の数が数十個以上の場合）に、特に有益となる。

一実施形態の血圧関連情報表示装置では、
上記統計処理部が求める上記波形ばらつきは、上記複数の個別波形の血圧変動量データがなす分布の四分位範囲として定義され、
上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域を表示することを特徴とする。

この一実施形態の血圧関連情報表示装置では、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域が表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。また、上記波形ばらつきを示す領域を標準偏差を用いて定義する場合に比して、多くの計算を必要とせず、計算が簡単で済む。したがって、処理速度の向上やメモリの節約に貢献することができる。

一実施形態の血圧関連情報表示装置では、上記表示処理部は、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線とは区別可能な態様で、上記複数の個別波形を示す曲線を表示する処理を行うことを特徴とする。

この一実施形態の血圧関連情報表示装置では、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線と上記複数の個別波形を示す曲線とが表示される。したがって、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握できる。この複数の個別波形を示す曲線の表示は、例えば個別波形の数が数個以下の場合に、ユーザに個別波形のばらつきを把握させるのに有益となる。

一実施形態の血圧関連情報表示装置では、
上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態が特定された身体状態特定期間を表す情報を入力する入力部を備え、
上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求め、
上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うことを特徴とする。

本明細書で、「身体状態特定期間」とは、無呼吸期間、レム睡眠期間、ノンレム睡眠期間、覚醒期間、および／または、経皮的動脈血酸素飽和度（ＳｐＯ_２）が低い期間など、被験者が血圧サージの要因となり得る身体状態をとっている期間を指す。

この一実施形態の血圧関連情報表示装置では、入力部は、上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態特定期間を表す情報を入力する。上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める。上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う。これにより、ユーザは、上記身体状態特定期間毎に、言い換えれば、被験者が血圧サージの要因となり得る身体状態をとっている期間毎に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。したがって、ユーザは、例えば最も深刻な血圧サージを引き起こしている要因（身体状態）を把握することができる。

別の局面では、この開示の血圧関連情報表示方法は、
複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示方法であって、
被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する第１ステップと、
検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める第２ステップと、
複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める第３ステップとを備え、
上記第３ステップでは、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
上記第３ステップでは、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにし、
さらに、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う第４ステップを備える
ことを特徴とする。

この開示の血圧関連情報表示方法によれば、ユーザは、直感的に、上記被験者における血圧サージの代表的波形を示す曲線と波形ばらつきを示す領域を把握することができる。このことは、ＳＡＳの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。また、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記複数の個別波形の血圧変動量データを統計処理するので、上記代表的波形と上記波形ばらつきを把握し易くなる。また、上記座標平面上で上記複数の個別波形の横方向寸法が互いに幾つかの拍数相当分異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。

さらに別の局面では、この開示のプログラムは、上記血圧関連情報表示方法をコンピュータに実行させるプログラムである。

この開示のプログラムをコンピュータに実行させることによって、上記血圧関連情報表示方法を実施することができる。

以上より明らかなように、この開示の血圧関連情報表示装置および血圧関連情報表示方法によれば、血圧サージの代表的波形と波形ばらつきを重ね合わせて表示できる。また、この開示のプログラムをコンピュータに実行させることによって、上記血圧関連情報表示方法を実施することができる。

この発明の血圧関連情報表示装置をネットワーク上のシステムとして構成した一実施形態を示すブロック図である。 上記システムに含まれた病院端末の構成を示すブロック図である。 上記システムに含まれたサーバの構成を示すブロック図である。 上記システムによって、被験者の一晩の血圧の時系列データと血圧サージの要因となり得る身体状態とを併せて測定する際の態様を示す図である。 基本的な血圧関連情報表示方法を実行する際の上記サーバの動作フローを示す図である。 応用的な血圧関連情報表示方法を実行する際の上記サーバの動作フローを示す図である。 図６（Ａ）は被験者の一晩の血圧の時系列データを示す図である。図６（Ｂ）は、図６（Ａ）の一部を拡大して、その時系列データ上で検出された血圧サージと、ＰＳＧ装置で計測された無呼吸・低呼吸イベントを示す図である。 血圧サージの波形（個別波形）を例示して、血圧サージを検出するための判定基準を説明する図である。 血圧サージの個別波形を平均する処理を概念的に説明する図である。 図９（Ａ）は、複数の個別波形を例示する図である。図９（Ｂ）は、波形ばらつきを示す領域に平均的波形を示す曲線を重ね合わせた態様の画像データを例示する図である。 血圧サージの個別波形についての統計処理の詳細なフローを示す図である。 検出された血圧サージを特定する時刻を記録したテーブルを例示する図である。 或る血圧サージについて、個別波形の拍対応ピークの時系列データを記録したテーブルを例示する図である。 血圧変動量データを統計処理するために用いられるテーブルを例示する図である。 上記応用的な血圧関連情報表示方法を実行する際の上記病院端末での表示処理のフローを示す図である。 無呼吸期間、レム睡眠期間、ノンレム睡眠期間、覚醒期間についてそれぞれ作成された、波形ばらつきを示す領域に平均的波形を示す曲線を重ね合わせた態様の画像データを、表示画面上に並べてサムネイル表示した態様を示す図である。 上記無呼吸期間について作成された画像データを、上記表示画面上に拡大表示した態様を示す図である。 図１６（Ａ）は、図９（Ａ）と同様に、無呼吸期間について得られた複数の個別波形を例示する図である。図１６（Ｂ）は、波形ばらつきの表示の仕方の変形例を示す図である。

以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。

（システムの構成）
図１は、この発明の血圧関連情報表示装置をネットワーク上の一実施形態のシステム（全体を符号１００で示す。）として構成した例を示している。このシステム１００は、表示画面としての表示器２４０を有する病院端末２００Ａ，２００Ｂと、サーバ３００と、トノメトリ方式の血圧計４００と、ＰＳＧ（polysomnography；睡眠ポリグラフ検査）装置５００とを含んでいる。これらの病院端末２００Ａ，２００Ｂ、サーバ３００、血圧計４００およびＰＳＧ装置５００は、病院内ＬＡＮ（Local Area Network）であるネットワーク９００を介して互いに通信可能になっている。このネットワーク９００を介した通信は、無線、有線のいずれでも良い。この実施の形態において、ネットワーク９００は、病院内ＬＡＮ（Local Area Network）であるが、これに限定されず、インターネットのような他の種類のネットワークであってもよいし、ＵＳＢケーブルなどを用いた１対１の通信であってもよい。なお、病院端末２００Ａ，２００Ｂは、この例では２台のみを示しているが、３台以上設けられていてもよい。

図２に示すように、病院端末２００Ａは、本体２００Ｍと、この本体２００Ｍに搭載された、制御部２１０と、メモリ２２０と、操作部２３０と、表示画面としての表示器２４０と、通信部２９０とを含む。この病院端末２００Ａは、市販のノート型パーソナルコンピュータからなり、後述の処理を行わせるようにアプリケーションソフトウェア（コンピュータプログラム）をインストールしたものであり、サーバ３００にアクセスできるようになっている。

制御部２１０は、ＣＰＵ（Central Processing Unit；中央演算処理装置）およびその補助回路を含み、病院端末２００Ａの各部を制御し、メモリ２２０に記憶されたプログラムおよびデータに従って後述の処理を実行する。すなわち、操作部２３０、および、通信部２９０から入力されたデータを処理し、処理したデータを、メモリ２２０に記憶させたり、表示器２４０で表示させたり、通信部２９０から出力させたりする。

メモリ２２０は、制御部２１０でプログラムを実行するために必要な作業領域として用いられるＲＡＭ（Random Access Memory）と、制御部２１０で実行するための基本的なプログラムを記憶するためのＲＯＭ（Read Only Memory）とを含む。また、メモリ２２０の記憶領域を補助するための補助記憶装置の記憶媒体として、半導体メモリ（メモリカード、ＳＳＤ（Solid State Drive））などが用いられてもよい。

操作部２３０は、この例では、キーボードおよびマウスで構成され、典型的にはユーザとしての医師による操作を示す操作信号を制御部２１０に入力する。また、操作部２３０は、キーボードおよびマウスに替えて、または、加えて、タッチパネルなどの他の操作デバイスで構成されるようにしてもよい。

表示器２４０は、表示画面（例えば、ＬＣＤ（Liquid Crystal Display）またはＥＬ（Electroluminescence）ディスプレイなど）を含む。表示器２４０は、制御部２１０によって制御されて、所定の映像を表示画面に表示させる。

通信部２９０は、制御部２１０からの情報をネットワーク９００を介してサーバ３００へ送信する。

簡単のため図示を省略するが、他の病院端末２００Ｂ，…も病院端末２００Ａと同様の構成になっている。

図３に示すように、サーバ３００は、制御部３１０と、記憶部３２０と、操作部３３０と、表示器３４０と、通信部３９０とを含む。このサーバ３００は、汎用のコンピュータ装置に、後述の処理を行わせるようにプログラム（ソフトウェア）をインストールしたものである。

制御部３１０は、ＣＰＵおよびその補助回路を含み、サーバ３００の各部を制御し、記憶部３２０に記憶されたプログラムおよびデータに従って所定の処理を実行し、操作部３３０、および、通信部３９０から入力されたデータを処理し、処理したデータを、記憶部３２０に記憶させたり、表示器３４０で表示させたり、通信部３９０から出力させたりする。

記憶部３２０は、制御部３１０でプログラムを実行するために必要な作業領域として用いられるＲＡＭと、制御部３１０で実行するための基本的なプログラムを記憶するためのＲＯＭとを含む。記憶部３２０には、多くの被験者から送られてきた血圧測定データを含むデータベース３２１が設けられている。また、記憶部３２０の記憶領域を補助するための補助記憶装置の記憶媒体として、磁気ディスク（ＨＤ（Hard Disk）、ＦＤ（Flexible Disk））、光ディスク（ＣＤ（Compact Disc）、ＤＶＤ（Digital Versatile Disk）、ＢＤ（Blu-ray Disc））、光磁気ディスク（ＭＯ（Magneto-Optical disk））、または、半導体メモリ（メモリカード、ＳＳＤ）などが用いられてもよい。

操作部３３０は、この例では、キーボードおよびマウスで構成され、ユーザによる操作を示す操作信号を制御部３１０に入力する。また、操作部３３０は、キーボードおよびマウスに替えて、または、加えて、タッチパネルなどの他の操作デバイスで構成されるようにしてもよい。

表示器３４０は、表示画面（例えば、ＬＣＤまたはＥＬディスプレイなど）を含む。表示器３４０は、制御部３１０によって制御されて、所定の映像を表示画面に表示させる。

通信部３９０は、制御部３１０からの情報をネットワーク９００を介して他の装置（この例では病院端末２００Ａ）へ送信するとともに、他の装置からネットワーク９００を介して送信されてきた情報を受信して制御部３１０に受け渡す。

図１中に示す血圧計４００は、この例では、特願２０１７－０５００６６号に開示されているようなトノメトリ方式の血圧計からなっている。この血圧計４００は、図４中に示すように、トノメトリ方式で被験者９０の被測定部位（例えば、左手首）を通る橈骨動脈の圧脈波を一拍毎に連続的に検出する圧力センサ部４０２と、この圧力センサ部４０２が検出した圧力の変化を血圧の時系列データとして出力する本体ユニット４０１とを備えている。圧力センサ部４０２と本体ユニット４０１との間は、信号ケーブル４０３によって接続されている。トノメトリ法は血管を圧力センサ部４０２（例えば、圧脈波センサ）で圧扁することにより圧脈波を計測し血圧を決定する手法である。血管の厚さが一様な円管と見なすと、血管内の血液の流れ、拍動の有無に関係なく血管壁を考慮してラプラスの法則に従い、血管の内圧（血圧）と血管の外圧（圧脈波の圧力）との関係式を導くことができる。この関係式で押圧面において血管が圧扁されている条件下では、血管の外壁及び内壁の半径を近似することにより、圧脈波の圧力と血圧とが等しいと近似できる。従って以降は、圧脈波の圧力は血圧と同一値になる。この結果、血圧計４００は、被測定部位（例えば、左手首）の血圧値を一心拍ごとに測定して、例えば図６（Ａ）に示すような態様で、測定時刻（時間）と血圧とを対応付けた血圧の時系列データ８０１を出力するようになっている。一晩の時系列データ８０１は、約３万拍分の拍対応ピーク（収縮期血圧（ＳＢＰ）または拡張期血圧（ＤＢＰ）に相当するピーク）を含んでいる。

また、図１中に示すＰＳＧ装置５００は、この例では、市販のＰＳＧ装置（例えば、日本光電工業株式会社製のＮｅｕｒｏｆａｘ（登録商標）ＥＥＧ－９２００）からなり、後述の［Ｂ．応用的な血圧関連情報表示方法］を実行する際にのみ用いられる。このＰＳＧ装置５００は、図４中に示すように、センサ群５０２と、このセンサ群５０２からの信号を処理して被験者９０の身体状態を特定する情報を出力する本体ユニット５０１とを備えている。この例では、センサ群５０２は、脳波を検出する脳波検出電極５０２Ａと、眼球運動を検出する眼球運動検出電極５０２Ｂと、呼吸に伴う気流を検出する気流センサ５０２Ｃと、心電図を得るための心電図電極５０２Ｄと、経皮的動脈血酸素飽和度（ＳｐＯ_２）を検出するパルスオキシメータ５０２Ｅと、筋電図を得るための筋電図電極５０２Ｆとを含んでいる。センサ群５０２と本体ユニット５０１との間は、図示しないケーブルボックスを介して図示しない信号ケーブルによって接続されている。この例では、ＰＳＧ装置５００は、被験者９０の身体状態が特定された期間（これを「身体状態特定期間」と呼ぶ。）を表す情報として、無呼吸期間（若しくは低呼吸期間）、レム睡眠期間、ノンレム睡眠期間、覚醒期間、および／または、ＳｐＯ_２が低い期間を表す情報を出力することができる。ここで、睡眠中の「無呼吸」とは、「呼吸が１０秒以上止まる状態」を指す。また、「低呼吸」とは、呼吸による換気が１０秒以上５０％以下に低下することを指す。また、「レム睡眠」とは急速眼球運動を伴う睡眠を指し、「ノンレム睡眠」とは急速眼球運動を伴わない睡眠を指し、「覚醒」とは目覚めている状態を指す。ＳｐＯ_２とは、動脈血の中を流れている赤血球に含まれるヘモグロビンの何％に酸素が結合しているかを、皮膚を通して（経皮的に）測定した値を指す。ＳｐＯ_２が「低い」期間とは、この例ではＳｐＯ_２が９０％未満となった期間を指す。

例えば図６（Ａ）の一部を拡大して示す図６（Ｂ）中に、「無呼吸」期間が、斜線付きのバー８０５，８０５，…で示されている。また、「低呼吸」期間が、白抜きのバー８０４，８０４，…で示されている。このような態様で、ＰＳＧ装置５００は、被験者９０についての身体状態特定期間を表す情報を複数種類出力することができる。

（血圧関連情報表示方法）
このシステム１００は、大別して、サーバ３００での画像データの作成と、病院端末（例えば２００Ａ）での画像データの表示とを含む、以下のような［Ａ．基本的な血圧関連情報表示方法］と［Ｂ．応用的な血圧関連情報表示方法］を実施することができる。

［Ａ．基本的な血圧関連情報表示方法］
この基本的な血圧関連情報表示方法では、ＰＳＧ装置５００を用いることなく、図６（Ａ）に例示したような、血圧計４００が出力する測定時刻（時間）と血圧とを対応付けた血圧の時系列データ８０１のみを用いて、次のように画像データの作成・表示を行う。

（サーバ３００での画像データの作成）
ｉ） 図４に示した態様で、被験者９０について、一晩（睡眠期間を含む。）にわたって、血圧計４００を用いた連続的な血圧測定が行われたものとする（この例では、ＰＳＧ装置５００を用いた身体状態の測定は行われないものとする。）。この例では、血圧計４００からの血圧の時系列データが、病院端末２００Ａによって受信され、メモリ２２０に一旦記憶される。続いて、例えばユーザとしての医師が病院端末２００Ａの操作部２３０を操作して、被験者９０について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク９００を介してサーバ３００に送信する。

ii） サーバ３００は、この例では常時動作状態にあり（ただし、メンテナンス期間等を除く。）、図５Ａ（基本的な血圧関連情報表示方法を実行する際のサーバ３００の動作フローを示す。）のステップＳ１０１に示すように、病院端末２００Ａからのデータを待っている。サーバ３００の制御部３１０は、ネットワーク９００から入力部としての通信部３９０を介して、病院端末２００Ａからの上述の血圧の時系列データを受信すると（ステップＳ１０１でＹＥＳ）、受信した血圧の時系列データを記憶部３２０のデータベース３２１に記憶する。そして、次に述べるステップＳ１０２～Ｓ１０７の処理を実行する。

iii） まず、ステップＳ１０２に示すように、制御部３１０は前処理部として働いて、血圧の時系列データに公知の移動平均などなどを用いた平滑化、ノイズ除去、ローパスフィルタを用いた高周波成分除去等の前処理を施す。

iv） 次に、ステップＳ１０３に示すように、制御部３１０は血圧サージ検出部として働いて、被験者９０の血圧の時系列データから、例えば特願２０１７－０４８９４６号、特願２０１７－０５００６６号に開示されているように、予め定められた判定基準に基づいて血圧サージを検出する。これにより、例えば図６（Ｂ）中に、破線の矩形枠８０３，８０３，…で示すように、複数の血圧サージを検出する。血圧サージは一晩に数百個も発生することがあると言われている。この例では、収縮期血圧（ＳＢＰ）に相当するピークに黒丸●が付され、それらが包絡線で結ばれている。また、拡張期血圧（ＤＢＰ）に相当するピークに白丸○が付され、それらが包絡線で結ばれている。

ここで、血圧サージを検出するための「予め定められた判定基準」とは、例えば図７（血圧サージの個別波形の一例を曲線Ｃで示す。）に示すように、ピーク検出区間（例えば、１５拍分の期間）内にサージ開始点Ｐ１の時刻からサージピーク点Ｐ２の時刻までが含まれるとともに、サージ開始点Ｐ１の時刻における収縮期血圧値（ＳＢＰ）とピーク点Ｐ２の時刻における収縮期血圧値（ＳＢＰ）との差（血圧変動量）Ｌ１が２０ｍｍＨｇ（または１５ｍｍＨｇ）以上であること、サージ開始点Ｐ１の時刻とピーク点Ｐ２の時刻との間の期間Ｔ１が５拍分よりも大きいこと、かつ、ピーク点Ｐ２の時刻とサージ終了点Ｐ４の時刻との間の期間Ｔ３が７拍分よりも大きいことを指す。この例では、サージ開始点Ｐ１は、ピーク点Ｐ２よりも前で、収縮期血圧値（ＳＢＰ）の極小を与える点として規定されている。また、サージ終了点Ｐ４は、ピーク点Ｐ２よりも後で、ピーク点Ｐ２からＬ１×０．７５（＝Ｌ３）の分だけ血圧が降下した点として規定されている。

ｖ）次に、図５ＡのステップＳ１０４に示すように、制御部３１０は個別波形取得部として働いて、検出された各血圧サージ８０３，８０３，…について、その血圧サージ８０３をなす複数の拍対応ピーク（この例では、収縮期血圧（ＳＢＰ）に相当するピーク）を結ぶ包絡線を個別波形として求める。図７に示したように、血圧サージの個別波形は、山状の形状をなす曲線Ｃとして表される。

vi）次に、図５ＡのステップＳ１０５に示すように、制御部３１０は統計処理部として働いて、この例では、得られた全ての個別波形を統計処理して、上記時系列データ８０１内での血圧サージの代表的波形と波形ばらつきを求める。ここで、「統計処理」とは、個別波形を平均する処理を指す。この例では、血圧サージの「代表的波形」とは、複数の個別波形を平均して得られた平均的波形を指す。また、この例では、血圧サージの「波形ばらつき」とは、複数の個別波形の分布の幅を指す。

概念的に説明すると、血圧サージの個別波形を平均する処理とは、次のような処理を意味する。図８中に例示するように、血圧サージとして「サージＮｏ．１」の個別波形Ｃ１、「サージＮｏ．２」の個別波形Ｃ２が存在するものとする。この場合、血圧サージの個別波形を平均する処理とは、時間経過を表す横方向座標Ｘと、血圧サージに伴う血圧変動量を表す縦方向座標Ｙとを含む座標平面Ｑ上で、横方向Ｘに関して複数の個別波形Ｃ１，Ｃ２を相対的にスライドさせて複数の個別波形Ｃ１，Ｃ２のピークの位置を揃えた状態で、横方向座標Ｘ毎に複数の個別波形Ｃ１，Ｃ２の血圧変動量データを平均して、代表的波形としての平均的波形を示す曲線Ｃａｖを求める処理に相当する。ここで、図８では、時間経過を表す横方向座標Ｘは、拍対応ピークを特定する拍番号によって定められている。縦方向座標Ｙは血圧サージに伴う収縮期血圧値（ＳＢＰ）の血圧変動量によって定められている。なお、具体的なデータ処理（血圧サージの個別波形の統計処理）については、後に詳述する。

また、この例では、血圧サージの「波形ばらつき」の範囲は、複数の個別波形の血圧変動量データがなす分布（横方向座標Ｘ毎の分布）の標準偏差σを用いて、横方向座標Ｘ毎に標準偏差σの±ｋ倍（ただし、ｋは自然数である。）の範囲として定義される。典型的には、ｋは１、２または３として設定される。以下では、特段断らない限り、血圧サージの「波形ばらつき」の範囲は、±σの範囲とする。この波形ばらつき±σを示す領域（例えば、後述の図９（Ｂ）中に示す領域Ｓｄ）の表示は、個別波形のばらつきが正規分布として取り扱われるのに十分である場合（例えば、数十個以上の場合）に、特に有益な表示となる。

これにより、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとが得られる。

vii）次に、図５ＡのステップＳ１０６に示すように、制御部３１０は表示処理部の一部として働いて、血圧サージの波形ばらつきを示す領域（符号Ｓｄで表す。）に、平均的波形を示す曲線Ｃａｖを重ね合わせた態様の画像データを作成する。

例えば、図９（Ａ）に示すように、図５ＡのステップＳ１０４で個別波形Ｃ１１，Ｃ１２，Ｃ１３，…が得られたものとする。この場合、図５ＡのステップＳ１０６によって、例えば図９（Ｂ）に示すように、波形ばらつきを示す領域Ｓｄに、平均的波形を示す曲線Ｃａｖを重ね合わせた態様の画像データＩｍが得られる。この例では、平均的波形を示す曲線Ｃａｖは、或る一定の濃度を有する実線で表されている。波形ばらつきを示す領域Ｓｄは、曲線Ｃａｖと背景領域（この例では、白色領域）Ｂｇとの間の中間濃度を有する領域として表されている。

viii）この後、図５ＡのステップＳ１０７に示すように、サーバ３００の制御部３１０は、作成した画像データＩｍを、ネットワーク９００を介して、この例では、血圧の時系列データの提供元である病院端末２００Ａへ送信する。

なお、サーバ３００の制御部３１０は、画像データＩｍを、例えば病院端末２００Ａを操作するユーザの指定に応じて、病院端末２００Ａ以外の他の病院端末２００Ｂ等へ送信してもよい。

（血圧サージの個別波形の統計処理）
図１０は、図５ＡのステップＳ１０５で述べた、サーバ３００の制御部３１０による血圧サージの個別波形についての統計処理の詳細なフローを示している。

図１０のステップＳ１１１に示すように、制御部３１０は、データベース３２１から、検出された全ての血圧サージについて、個別波形の拍対応ピークの時系列データを取得する。この例では、取得された血圧サージの数はＮ個であったものとする。

例えば、図１１のテーブルＭＴ１に示すように、検出された血圧サージＮｏ．１は、開始時刻（開始点の時刻）２２：２０：１５、ピーク時刻（ピークの時刻）２２：２０：２１、終了時刻（終了点の時刻）２２：２０：２６で特定されている。また、血圧サージＮｏ．２は、開始時刻２２：２５：３５、ピーク時刻２２：２５：４４、終了時刻２２：２５：５２で特定されている。なお、時刻表記の仕方は、「時：分：秒」である（以下同様）。

この場合、例えば血圧サージＮｏ．１について、個別波形の拍対応ピークの時系列データは、例えば図１２のテーブルＭＴ２に示すように、開始時刻２２：２０：１５では、拍対応ピークとしてのＳＢＰに相当するピークが１３０［ｍｍＨｇ］、ＤＢＰに相当するピークが７８［ｍｍＨｇ］、脈拍数ＰＲが６１［拍／分］というように記録されている。次の拍に相当する時刻２２：２０：１６では、ＳＢＰに相当するピークが１３４［ｍｍＨｇ］、ＤＢＰに相当するピークが８０［ｍｍＨｇ］、脈拍数ＰＲが６０［拍／分］というように記録されている。このようにして順次拍対応ピークの時系列データが記録され、終了時刻２２：２０：２６では、ＳＢＰに相当するピークが１３３［ｍｍＨｇ］、ＤＢＰに相当するピークが８２［ｍｍＨｇ］、脈拍数ＰＲが６１［拍／分］というように記録されている。簡単のため図示を省略するが、血圧サージＮｏ．２についても、個別波形の拍対応ピークの時系列データは、同様に記録されている。

次に、図１０のステップＳ１１２に示すように、制御部３１０は、Ｎ個の個別波形について、それぞれその個別波形の開始点Ｐ１を基準（＝０［ｍｍＨｇ］）としてオフセットする。すなわち、Ｎ個の個別波形について、それぞれ、各時刻の血圧値（この例では、ＳＢＰ値）からその個別波形の開始点Ｐ１の血圧値（この例では、ＳＢＰ値）を差し引く。これにより、時刻毎の血圧変動量データを得る。このようにオフセットを行うことによって、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとを把握し易くなる。

次に、図１０のステップＳ１１３に示すように、制御部３１０は、Ｎ個の個別波形のうちＸ方向（横方向）寸法が最長の個別波形を特定する。また、その最長のＸ方向寸法を拍数Ｍ個として求める。

次に、図１０のステップＳ１１４に示すように、制御部３１０は、記憶部３２０に、演算を実行するための作業領域としてＮ行Ｍ列のメモリ領域を確保する。

ここで簡単のため、図８中に例示したように、取得された血圧サージの数Ｎ＝２個であり、「サージＮｏ．１」の個別波形Ｃ１のＸ方向寸法が拍数１２個分であり、「サージＮｏ．２」の個別波形Ｃ２のＸ方向寸法が拍数１７個分であったとする。このとき、制御部３１０は、Ｘ方向寸法が最長の個別波形をＣ２であると特定し、その最長のＸ方向寸法を拍数Ｍ＝１７個として求める。そして、制御部３１０は、記憶部３２０に、図１３のテーブルＭＴ３の２行目～３行目に示すように、「サージＮｏ．１」、「サージＮｏ．２」のための演算を実行するための作業領域として２行１７列のメモリ領域を確保する。なお、実際に、取得された血圧サージの数Ｎが２より大きい場合は、２より多い行数のメモリ領域が確保される。テーブルＭＴ３の最上段（１行目）には、Ｘ方向座標を表す拍番号が表されている。

次に、図１０のステップＳ１１５に示すように、制御部３１０は、Ｎ行Ｍ列のメモリ領域において横方向Ｘに関してＮ個の個別波形を相対的にスライドさせてＮ個の個別波形のピークの位置（拍番号）を揃える。例えば図１３の例では、横方向Ｘに関して「サージＮｏ．１」、「サージＮｏ．２」の２個の個別波形Ｃ１，Ｃ２のピークの位置が、拍番号１０（破線の丸印で囲んで示す）の所に揃えられている。この結果、テーブルＭＴ３の２行目に、「サージＮｏ．１」の血圧変動量データが、拍番号４～１５にわたって拍番号毎に「０」、「４」、「１０」、「１３」、…、「１１」、「３」というように表されている（単位はｍｍＨｇ。以下同様。）。また、テーブルＭＴ３の３行目に、「サージＮｏ．２」の血圧変動量データが、拍番号１～１７にわたって拍番号毎に「０」、「２」、「３」、「５」、…、「４」、「２」というように表されている。このように、Ｎ個の個別波形のピークの位置（拍番号）を揃えた状態は、図８中に示した座標平面Ｑ上で、横方向Ｘに関して複数の個別波形Ｃ１，Ｃ２を相対的にスライドさせて複数の個別波形Ｃ１，Ｃ２のピークの位置を揃えた状態と等価になっている。これにより、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとを、さらに把握し易くなる。なお、「サージＮｏ．１」の個別波形Ｃ１の開始点（拍番号４）の値、「サージＮｏ．２」の個別波形Ｃ２の開始点（拍番号１）の値がいずれも「０」になっているのは、既述のオフセット（図１０のステップＳ１１２）の結果である。

次に、図１０のステップＳ１１６に示すように、制御部３１０は、Ｎ個の個別波形の血圧変動量データを、拍番号（列）毎に統計処理（この例では、平均化）して、平均値を求める。例えば図１３の例では、テーブルＭＴ３の最下段（４行目）に、上述の「サージＮｏ．１」、「サージＮｏ．２」の個別波形の場合の血圧変動量データについての拍番号毎の平均値が表されている。この拍番号の平均値によって、代表的波形としての平均的波形を示す曲線Ｃａｖが定められる。

このとき、制御部３１０は、Ｎ個の個別波形のうちＸ方向寸法が最長の個別波形（上の例では、「サージＮｏ．２」）よりも短い個別波形（上の例では、「サージＮｏ．１」）について、最長の個別波形（上の例では、「サージＮｏ．２」）よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする。例えば図１３の例では、「サージＮｏ．１」の個別波形について、Ｘ方向寸法が「サージＮｏ．２」の個別波形よりも不足する部分、すなわち、拍番号１～３の部分、および、拍番号１６～１７の部分について、上記血圧変動量データの統計処理に対する寄与をゼロにする。これにより、座標平面Ｑ上でＮ個の個別波形のＸ方向寸法が互いに異なる場合であっても、上記血圧変動量データの統計処理を支障無く行うことができる。

さらに、図１０のステップＳ１１７に示すように、制御部３１０は、Ｎ個の個別波形の血圧変動量データがなす分布の標準偏差σを、拍番号（列）毎に求める。拍番号毎に±σの領域が、既述の波形ばらつきを示す領域Ｓｄを定める。

このように、このフローでは、Ｎ個の個別波形の血圧変動量データを統計処理する際に、拍番号（列）毎に血圧変動量データの統計処理を行えば足りる。したがって、横方向Ｘに関して連続的に血圧変動量データの統計処理を行う場合に比して、制御部３１０の演算量が少なくて済む。

なお、上述の図１０のステップＳ１１２で、Ｎ個の個別波形について、それぞれその個別波形の開始点Ｐ１を基準（＝０［ｍｍＨｇ］）としてオフセットするのに加えて、さらに、Ｎ個の個別波形のピークの高さを揃えるように正規化してもよい。これにより、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとを、さらに把握し易くなる。

（病院端末での画像データの表示）

病院端末２００Ａは、図９（Ｂ）に示したような画像データＩｍを、ネットワーク９００を介して、サーバ３００から受信する。

次に、病院端末２００Ａの操作部２３０を介したユーザの表示指示があると、病院端末２００Ａの制御部２１０は、表示処理部の一部として働いて、表示器２４０の表示画面上に、その画像データＩｍを表示する。このとき、表示器２４０の表示画面上に、血圧変動量（図７中に示したサージ開始点Ｐ１の時刻における収縮期血圧値（ＳＢＰ）とピーク点Ｐ２の時刻における収縮期血圧値（ＳＢＰ）との差Ｌ１）の平均値、分布の標準偏差σ（または±σ，±２σ，±３σ）の値などを併せて表示してもよい。

ユーザは、表示器２４０の表示画面を見ることによって、被験者９０における血圧サージの代表的波形を示す曲線Ｃａｖと波形ばらつきを示す領域Ｓｄを把握することができる。このことは、ＳＡＳの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。

なお、ユーザは、病院端末２００Ａの操作部２３０を介して転送指示を入力して、その画像データＩｍを、病院端末２００Ａ以外の他の病院端末２００Ｂ等へ送信してもよい。

［Ｂ．応用的な血圧関連情報表示方法］
この応用的な血圧関連情報表示方法では、血圧計４００が出力する測定時刻（時間）と血圧とを対応付けた血圧の時系列データと、ＰＳＧ装置５００が出力する身体状態特定期間を表す情報とを用いて、次のように画像データの作成・表示を行う。

（サーバ３００での画像データの作成）
ｉ） 図４に示した態様で、被験者９０について、一晩（睡眠期間を含む。）にわたって、血圧計４００を用いた連続的な血圧測定と、ＰＳＧ装置５００を用いた身体状態の測定が行われたものとする。この例では、血圧計４００からの血圧の時系列データと、ＰＳＧ装置５００からの身体状態特定期間を表す情報とが、病院端末２００Ａによって受信され、メモリ２２０に一旦記憶される。続いて、例えばユーザとしての医師が病院端末２００Ａの操作部２３０を操作して、被験者９０について測定された血圧の時系列データと身体状態特定期間を表す情報とを併せて、ネットワーク９００を介してサーバ３００に送信する。

ii） サーバ３００は、この例では常時動作状態にあり（ただし、メンテナンス期間等を除く。）、図５Ｂ（応用的な血圧関連情報表示方法を実行する際のサーバ３００の動作フローを示す。）のステップＳ１５１に示すように、病院端末２００Ａからのデータを待っている。サーバ３００の制御部３１０は、ネットワーク９００から入力部としての通信部３９０を介して、病院端末２００Ａからの上述のデータ（血圧の時系列データと身体状態特定期間を表す情報）を受信すると（ステップＳ１５１でＹＥＳ）、受信した血圧の時系列データと身体状態特定期間を表す情報を記憶部３２０のデータベース３２１に記憶する。そして、次に述べるステップＳ１５２～Ｓ１５７の処理を実行する。

iii） まず、ステップＳ１５２に示すように、制御部３１０は前処理部として働いて、血圧の時系列データに公知の移動平均などなどを用いた平滑化、ノイズ除去、ローパスフィルタを用いた高周波成分除去等の前処理を施す。

iv） 次に、ステップＳ１５３に示すように、制御部３１０は血圧サージ検出部として働いて、被験者９０の血圧の時系列データから、例えば特願２０１７－０４８９４６号、特願２０１７－０５００６６号に開示されているように、予め定められた判定基準に基づいて血圧サージを検出する。これにより、例えば図６（Ｂ）中に、破線の矩形枠８０３，８０３，…で示すように、複数の血圧サージを検出する。ここで、血圧サージを検出するための「予め定められた判定基準」は、図５ＡのステップＳ１０３に関して図７を用いて説明した判定基準と同様である。

ｖ）次に、図５ＢのステップＳ１５４に示すように、制御部３１０は個別波形取得部として働いて、検出された各血圧サージ８０３，８０３，…について、その血圧サージ８０３をなす複数の拍対応ピーク（この例では、収縮期血圧（ＳＢＰ）に相当するピーク）を結ぶ包絡線を個別波形として求める。図７に示したように、血圧サージの個別波形は、山状の形状をなす曲線Ｃとして表される。

vi）次に、図５ＢのステップＳ１５５に示すように、制御部３１０は統計処理部として働いて、この例では既述の身体状態特定期間毎に、それぞれ複数得られた個別波形を統計処理して、上記時系列データ８０１内での血圧サージの代表的波形と波形ばらつきを求める。この例では、或る身体状態特定期間（例えば、無呼吸期間）が血圧サージの開始点Ｐ１から終了点Ｐ４までの期間と一部でも重なっていれば、その血圧サージはその身体状態特定期間に発生したものとして処理される。ここで、「統計処理」とは、図５ＡのステップＳ１０５に関して図８（および図１０～図１３）を用いて説明したのと同様に、個別波形を平均する処理を指す。すなわち、血圧サージの個別波形を平均する処理は、時間経過を表す横方向座標Ｘと、血圧サージに伴う血圧変動量を表す縦方向座標Ｙとを含む座標平面Ｑ上で、横方向Ｘに関して複数の個別波形Ｃ１，Ｃ２を相対的にスライドさせて複数の個別波形Ｃ１，Ｃ２のピークの位置を揃えた状態で、横方向座標Ｘ毎に複数の個別波形Ｃ１，Ｃ２の血圧変動量データを平均して、代表的波形としての平均的波形を示す曲線Ｃａｖを求める処理に相当する。また、この例では、図５ＡのステップＳ１０５に関して説明したのと同様に、血圧サージの「波形ばらつき」の範囲は、複数の個別波形の血圧変動量データがなす分布（横方向座標Ｘ毎の分布）の標準偏差σを用いて、横方向座標Ｘ毎に±σの範囲として定義される。

これにより、身体状態特定期間毎に、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとが得られる。例えば、一晩の時系列データで特定された身体状態特定期間が無呼吸期間（低呼吸期間を含む。）、レム睡眠期間、ノンレム睡眠期間、覚醒期間の４種類である場合、無呼吸期間について、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとが得られる。同様に、レム睡眠期間、ノンレム睡眠期間、覚醒期間のそれぞれについて、代表的波形としての平均的波形を示す曲線Ｃａｖと、横方向座標Ｘ毎の波形ばらつき±σとが得られる。

vii）次に、図５ＢのステップＳ１５６に示すように、制御部３１０は表示処理部の一部として働いて、既述の身体状態特定期間毎に、図９（Ｂ）に例示したような、血圧サージの波形ばらつきを示す領域Ｓｄに、平均的波形を示す曲線Ｃａｖを重ね合わせた態様の画像データＩｍを作成する。

例えば、後述の図１５Ａ中に示すように、無呼吸期間（低呼吸期間を含む。）、レム睡眠期間（ＲＥＭ）、ノンレム睡眠期間（ＮＲＥＭ）、覚醒期間のそれぞれについて、波形ばらつきを示す領域Ｓｄ１，Ｓｄ２，Ｓｄ３，Ｓｄ４に、それぞれ平均的波形を示す曲線Ｃａｖ１，Ｃａｖ２，Ｃａｖ３，Ｃａｖ４を重ね合わせた態様の画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４が得られる。

viii）この後、図５ＢのステップＳ１５７に示すように、サーバ３００の制御部３１０は、身体状態特定期間毎に作成した画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４を、ネットワーク９００を介して、この例では、血圧の時系列データの提供元である病院端末２００Ａへ送信する。

なお、サーバ３００の制御部３１０は、画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４を、例えば病院端末２００Ａを操作するユーザの指定に応じて、病院端末２００Ａ以外の他の病院端末２００Ｂ等へ送信してもよい。

（病院端末での画像データの表示）
図１４は、サーバ３００から身体状態特定期間毎に作成された画像データを受信した場合の、病院端末２００Ａでの表示処理のフローを示している。なお、この例では、病院端末２００Ａについて説明するが、他の病院端末２００Ｂ等も同様の表示処理が可能になっている。

病院端末２００Ａは、まず図１４のステップＳ２０１に示すように、身体状態特定期間毎に作成された画像データを、ネットワーク９００を介して、サーバ３００から受信する。この例では、身体状態特定期間は、無呼吸期間（低呼吸期間を含む。）、レム睡眠期間、ノンレム睡眠期間、覚醒期間の４種類であるものとする。

次に、病院端末２００Ａの制御部２１０は、表示処理部の一部として働いて、図１５Ａに示すように、表示器２４０の表示画面上に、身体状態特定期間毎に作成された画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４を並べてサムネイル表示する（図１４のステップＳ２０２）。図１５Ａ中に示すように、この例では、表示器２４０の表示画面には、身体状態特定期間を特定する身体状態特定欄２４１と、その下方に対応して配置された画像データ欄２４２とが表示されている。この例では、身体状態特定欄２４１には、左から右へ向かって、無呼吸期間を表す「無呼吸」欄２４１ａと、レム睡眠期間を表す「ＲＥＭ」欄２４１ｂと、ノンレム睡眠期間を表す「ＮＲＥＭ」欄２４１ｃと、覚醒期間を表す「覚醒反応」欄２４１ｄとが含まれている。画像データ欄２４２には、無呼吸期間について作成された画像データＩｍ１を表示する欄２４２ａと、レム睡眠期間について作成された画像データＩｍ２を表示する欄２４２ｂと、ノンレム睡眠期間について作成された画像データＩｍ３を表示する欄２４２ｃと、覚醒期間について作成された画像データＩｍ４を表示する欄２４２ｄとが含まれている。

この例では、無呼吸期間について作成された画像データＩｍ１は、波形ばらつきを示す領域Ｓｄ１に、平均的波形を示す曲線Ｃａｖ１を重ね合わせた態様になっている。レム睡眠期間について作成された画像データＩｍ２は、波形ばらつきを示す領域Ｓｄ２に、平均的波形を示す曲線Ｃａｖ２を重ね合わせた態様になっている。ノンレム睡眠期間について作成された画像データＩｍ３は、波形ばらつきを示す領域Ｓｄ３に、平均的波形を示す曲線Ｃａｖ３を重ね合わせた態様になっている。また、覚醒期間について作成された画像データＩｍ４は、波形ばらつきを示す領域Ｓｄ４に、平均的波形を示す曲線Ｃａｖ４を重ね合わせた態様になっている。

ユーザは、これらの画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４のサムネイル表示を見ることによって、身体状態特定期間毎に、言い換えれば、被験者９０が血圧サージの要因となり得る身体状態をとっている期間毎に、被験者９０における血圧サージの代表的波形を示す曲線Ｃａｖと波形ばらつきを示す領域Ｓｄを把握することができる。したがって、ユーザは、例えば最も深刻な血圧サージを引き起こしている要因（身体状態）を容易に把握することができる。例えば、画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４のうち、無呼吸期間について作成された画像データＩｍ１の血圧変動量が最も大きいことから、「無呼吸」が最も深刻な血圧サージを引き起こしていると判断してもよい。

次に、病院端末２００Ａの制御部２１０は、図１４のステップＳ２０３に示すように、表示されている画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４が選択されたか否かを判断する。ユーザが操作部２３０を操作することによって画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４のうちいずれかを選択した場合（図１４のステップＳ２０３でＹＥＳ）、制御部２１０は、図１５Ｂに示すように、表示器２４０の表示画面上に、その選択された身体状態特定期間についての画像データ（この例では、無呼吸期間について作成された画像データＩｍ１）を拡大表示する（図１４のステップＳ２０４）。このとき、表示器２４０の表示画面上に、その選択された身体状態特定期間についての血圧変動量（図７中に示したサージ開始点Ｐ１の時刻における収縮期血圧値（ＳＢＰ）とピーク点Ｐ２の時刻における収縮期血圧値（ＳＢＰ）との差Ｌ１）の平均値、分布の標準偏差σ（または±σ，±２σ，±３σ）の値などを併せて表示してもよい。

なお、画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４がいずれも選択されない場合（図１４のステップＳ２０３でＮＯ）、制御部２１０は、画像をオフして（図１４のステップＳ２０５）、この表示処理を終了する。

なお、ユーザは、病院端末２００Ａの操作部２３０を介して転送指示を入力して、画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４を、病院端末２００Ａ以外の他の病院端末２００Ｂ等へ送信してもよい。

上述のように、血圧サージの波形ばらつきを示す領域に平均的波形を示す曲線を重ね合わせた態様の画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４を身体状態特定期間毎に表示することは、ＳＡＳの診断や治療のほか、心血管の疾病リスクを評価する材料や、特定の臓器の疾病リスクを評価する材料として役立つと考えられる。

（変形例１）
上述の実施形態では、サーバ３００が画像データＩｍ（または、身体状態特定期間毎の画像データＩｍ１，Ｉｍ２，Ｉｍ３，Ｉｍ４。以下同様。）を作成し、病院端末２００Ａ，２００Ｂ，…の表示器２４０が画像を表示したが、これに限られるものではない。サーバ３００の制御部３１０は、画像データＩｍを病院端末２００Ａ，２００Ｂ，…へ送信する代わりに、画像作成のためのデータのみを送信し、病院端末２００Ａ，２００Ｂ，…の制御部２１０が専ら画像データＩｍを作成しても良い。

（変形例２）
また、上述の実施形態では、本発明の血圧関連情報表示装置を、病院端末２００Ａ，２００Ｂ，…とサーバ３００とを含むネットワーク上のシステム１００として構成したが、これに限られるものではない。

例えば、本発明の血圧関連情報表示装置を、病院端末２００Ａ，２００Ｂ，…のいずれかのみで構成しても良い。つまり、病院端末（例えば２００Ａ）に上述の血圧関連情報表示方法の全部（血圧計４００からの血圧の時系列データおよびＰＳＧ装置５００からの身体状態特定期間を表す情報の受信から、表示器２４０の表示画面上への画像データＩｍ等の表示までを含む。）を実行させてもよい。

その場合、病院端末２００Ａのメモリ２２０には、制御部２１０に、上述の血圧関連情報表示方法を実行させるプログラムをインストールしておく。これにより、本発明の血圧関連情報表示装置を小型かつコンパクトに構成することができる。

また、上述の血圧関連情報表示方法を、ソフトウェア（コンピュータプログラム）として、ＣＤ（コンパクトディスク）、ＤＶＤ（デジタル万能ディスク）、フラッシュメモリなどの非一時的（non-transitory）にデータを記憶可能な記録媒体に記録してもよい。このような記録媒体に記録されたソフトウェアを、パーソナルコンピュータ、ＰＤＡ（パーソナル・デジタル・アシスタンツ）、スマートフォンなどの実質的なコンピュータ装置にインストールすることによって、それらのコンピュータ装置に、上述の血圧関連情報表示方法を実行させることができる。

（変形例３）
上述の実施形態では、画像データＩｍを構成する波形ばらつきを示す領域Ｓｄは、横方向座標Ｘ毎に血圧変動量データがなす分布の標準偏差σの±ｋ倍（典型的には、±σ，±２σまたは±３σ）の範囲として定義された。しかしながら、波形ばらつきを示す領域Ｓｄの範囲は、標準偏差σを用いて定義されたものに限られない。

例えば、波形ばらつきを示す領域Ｓｄは、横方向座標Ｘ毎に血圧変動量データがなす分布の四分位範囲として定義されてもよい。すなわち、処理としては、横方向座標Ｘ毎に血圧変動量データを昇順に並べる。そして、第一四分位を領域Ｓｄの下限とし、第三四分位を領域Ｓｄの上限とする。これにより、波形ばらつきを示す領域Ｓｄを求めることができる。

このように、波形ばらつきを示す領域Ｓｄを、横方向座標Ｘ毎に血圧変動量データがなす分布の四分位範囲とした場合、標準偏差σを用いて定義する場合に比して、多くの計算を必要とせず、計算が簡単で済む。したがって、処理速度の向上やメモリの節約に貢献することができる。

また、血圧サージの「代表的波形」は、複数の個別波形を平均して得られた平均的波形としたが、これに限られるものではない。血圧サージの「代表的波形」は、横方向座標Ｘ毎に血圧変動量データがなす分布の第二四分位（中央値）であってもよい。この場合、処理としては、横方向座標Ｘ毎に血圧変動量データを昇順に並べ、横方向に中央値を連ねて、代表的波形を得ることかできる。

（変形例４）
上述の実施形態では、画像データＩｍを構成する波形ばらつきを示す領域Ｓｄは、表示器２４０の表示画面上で、例えば図９（Ｂ）に示したように、代表的波形を示す曲線Ｃａｖと背景領域Ｂｇとの間の中間濃度を有する領域として表示された。しかしながら、波形ばらつきを示す領域Ｓｄの表示態様は、中間濃度を有する領域に限られない。

例えば、図１６（Ａ）は、図９（Ａ）と同様に、複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を例示している。この場合において、図１６（Ｂ）に示すように、波形ばらつきを示す領域Ｓｄ′をなすように、代表的波形を示す曲線Ｃａｖ′とは区別可能な態様で、複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を示す曲線Ｃｉを表示してもよい。この例では、曲線Ｃａｖ′が実線で表示されているのに対して、複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を示す曲線Ｃｉはそれぞれ破線で表示されている。この例では、この波形ばらつきを示す領域Ｓｄ′（複数の曲線Ｃｉを含む）に代表的波形を示す曲線Ｃａｖ′を重ね合わせた態様の画像データが符号Ｉｍ′で表されている。

これにより、ユーザは、直感的に、被験者９０における血圧サージの代表的波形を示す曲線Ｃａｖ′と波形ばらつきを示す領域Ｓｄ′を把握できる。この複数の個別波形を示す曲線Ｃｉの表示は、例えば個別波形の数が数個以下の場合に、ユーザに個別波形Ｃ１１，Ｃ１２，Ｃ１３，…のばらつきを把握させるのに有益となる。

なお、代表的波形を示す曲線Ｃａｖ′と複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を示す曲線Ｃｉとは区別可能であれば良い。例えば、代表的波形を示す曲線Ｃａｖ′が赤色の実線で表示される一方、複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を示す曲線Ｃｉが青色の実線で表示されてもよい。

また、身体状態特定期間毎に取得された血圧サージの数Ｎに対して閾値Ｎα（例えば、Ｎα＝１０個とする。）を設けて、取得された血圧サージの数Ｎが閾値Ｎα以上であれば、図９（Ｂ）に示したような波形ばらつきを示す領域Ｓｄ（縦方向に±σの幅の中間濃度を有する領域）を含む第１の画像データＩｍの作成・表示を行う一方、取得された血圧サージの数Ｎが閾値Ｎα未満であれば、図１６（Ｂ）に示したような波形ばらつきを示す領域Ｓｄ′（複数の個別波形Ｃ１１，Ｃ１２，Ｃ１３，…を示す曲線Ｃｉ）を含む第２の画像データＩｍ′の作成・表示を行ってもよい。この血圧サージの数Ｎが閾値Ｎαよりも大きいか否かの判定は、例えば、サーバ３００の制御部３１０が血圧サージ数判定部として働いて行うことができる。

または、それに代えて、例えばサーバ３００の制御部３１０が、図９（Ｂ）に示したような波形ばらつきを示す領域Ｓｄを含む第１の画像データＩｍと、図１６（Ｂ）に示したような波形ばらつきを示す領域Ｓｄ′を含む第２の画像データＩｍ′との両方を作成してもよい。この場合、例えば、サーバ３００から病院端末（例えば２００Ａ）へ第１の画像データと第２の画像データとの両方を送信し、病院端末（例えば２００Ａ）の操作部２３０を介したユーザの入力に応じて、表示器２４０の表示画面上で、それらの第１の画像データＩｍの表示と第２の画像データＩｍ′の表示とを切り替えるようにしてもよい。

（変形例５）
また、上述の実施形態では、血圧計４００はトノメトリ方式の血圧計であるものとしたが、これに限られるものではない。血圧計４００は、被測定部位のうち対応する部分を通る動脈へ向けて光を照射する発光素子と、その光の反射光（または透過光）を受光する受光素子とを備えて、動脈の脈波を容積の変化に基づいて連続的に血圧を検出してもよい（光電方式）。また、血圧計４００は、被測定部位に当接された圧電センサを備えて、被測定部位のうち対応する部分を通る動脈の圧力による歪みを電気抵抗の変化として検出し、この電気抵抗の変化に基づいて連続的に血圧を検出してもよい（圧電方式）。さらに、血圧計４００は、被測定部位のうち対応する部分を通る動脈へ向けて電波（送信波）を送る送信素子と、その電波の反射波を受信する受信素子とを備えて、動脈の脈波による動脈とセンサとの間の距離の変化を送信波と反射波との間の位相のずれとして検出し、この位相のずれに基づいて連続的に血圧を検出してもよい（電波照射方式）。また、血圧を算出することができる物理量を観測することができれば、これらの以外の方式を適用してもよい。

以上の実施形態は例示であり、この発明の範囲から離れることなく様々な変形が可能である。上述した複数の実施の形態は、それぞれ単独で成立し得るものであるが、実施の形態同士の組みあわせも可能である。また、異なる実施の形態の中の種々の特徴も、それぞれ単独で成立し得るものであるが、異なる実施の形態の中の特徴同士の組みあわせも可能である。

１００ システム
２００Ａ，２００Ｂ，… 病院端末
２４０ 表示器
３００ サーバ
４００ 血圧計
５００ ＰＳＧ装置

Claims (10)

1. 複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示装置であって、
   被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する血圧サージ検出部と、
   検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める個別波形取得部と、
   複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める統計処理部と、
   表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う表示処理部とを備え、
   上記統計処理部は、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
   上記統計処理部は、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにする
   ことを特徴とする血圧関連情報表示装置。
2. 請求項１に記載の血圧関連情報表示装置において、
   上記複数の個別波形の血圧変動量データは、それぞれ、その個別波形の開始点の血圧値を基準とした変動量であることを特徴とする血圧関連情報表示装置。
3. 請求項１または２に記載の血圧関連情報表示装置において、
   上記複数の個別波形の血圧変動量データは、上記複数の個別波形のピークの高さを揃えるように正規化されていることを特徴とする血圧関連情報表示装置。
4. 請求項１から３までのいずれか一つに記載の血圧関連情報表示装置において、
   上記時間経過を表す横方向座標は、上記拍対応ピークを特定する拍番号によって定められ、
   上記統計処理部は、上記拍番号毎に上記複数の個別波形の血圧変動量データを統計処理することを特徴とする血圧関連情報表示装置。
5. 請求項１から４までのいずれか一つに記載の血圧関連情報表示装置において、
   上記統計処理部が求める上記波形ばらつきは、ｋを自然数としたとき、上記複数の個別波形の血圧変動量データがなす分布の標準偏差の±ｋ倍として定義され、
   上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記標準偏差の±ｋ倍に相当する幅をもつ領域を表示することを特徴とする血圧関連情報表示装置。
6. 請求項１から４までのいずれか一つに記載の血圧関連情報表示装置において、
   上記統計処理部が求める上記波形ばらつきは、上記複数の個別波形の血圧変動量データがなす分布の四分位範囲として定義され、
   上記表示処理部は、上記表示画面上に、上記波形ばらつきを示す領域として、縦方向に関して上記四分位範囲に相当する幅をもつ領域を表示することを特徴とする血圧関連情報表示装置。
7. 請求項１から４までのいずれか一つに記載の血圧関連情報表示装置において、
   上記表示処理部は、上記表示画面上で上記波形ばらつきを示す領域をなすように、上記代表的波形を示す曲線とは区別可能な態様で、上記複数の個別波形を示す曲線を表示する処理を行うことを特徴とする血圧関連情報表示装置。
8. 請求項１から７までのいずれか一つに記載の血圧関連情報表示装置において、
   上記被験者の脈動に連動して変化する血圧の時系列データに伴って、上記被験者の身体状態が特定された身体状態特定期間を表す情報を入力する入力部を備え、
   上記統計処理部は、上記身体状態特定期間毎に、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求め、
   上記表示処理部は、上記表示画面上で、上記身体状態特定期間毎に、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行うことを特徴とする血圧関連情報表示装置。
9. 複数拍にわたって血圧が上昇してピークを示し、そのピークに続いて複数拍にわたって血圧が下降する現象である血圧サージに関連する情報を可視化して表示する血圧関連情報表示方法であって、
   被験者の脈動に連動して変化する血圧の時系列データから、予め定められた判定基準に基づいて血圧サージを検出する第１ステップと、
   検出された各血圧サージについて、その血圧サージをなす複数の拍対応ピークを結ぶ包絡線を個別波形として求める第２ステップと、
   複数得られた上記個別波形を統計処理して、上記時系列データ内での血圧サージの代表的波形と波形ばらつきを求める第３ステップとを備え、
   上記第３ステップでは、時間経過を表す横方向座標と、血圧サージに伴う血圧変動量を表す縦方向座標とを含む座標平面上で、横方向に関して上記複数の個別波形を相対的にスライドさせて上記複数の個別波形のピークの位置を揃えたのと等価な状態をメモリ領域に作成し、この状態で上記横方向座標毎に上記複数の個別波形の血圧変動量データを統計処理して、上記代表的波形と上記波形ばらつきを求めるとともに、
   上記第３ステップでは、上記複数の個別波形のうち横方向寸法が最長の個別波形よりも短い個別波形について、上記最長の個別波形よりも不足する拍数相当分の、上記血圧変動量データの統計処理に対する寄与をゼロにし、
   さらに、表示画面上で、上記血圧サージの上記波形ばらつきを示す領域に上記代表的波形を示す曲線を重ね合わせて表示する処理を行う第４ステップを備える
   ことを特徴とする血圧関連情報表示方法。
10. 請求項９に記載の血圧関連情報表示方法をコンピュータに実行させるプログラム。

Images