▶ 株式会社リコーの特許一覧 ▶ 国立大学法人 東京医科歯科大学の特許一覧
特開2022-101483生体データ処理装置、生体データ計測システム、及びプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022101483
(43)【公開日】2022-07-06
(54)【発明の名称】生体データ処理装置、生体データ計測システム、及びプログラム
(51)【国際特許分類】
A61B 5/248 20210101AFI20220629BHJP
【ＦＩ】
A61B5/248
【審査請求】未請求
【請求項の数】17
【出願形態】ＯＬ
(21)【出願番号】P 2021192184
(22)【出願日】2021-11-26
(31)【優先権主張番号】P 2020214933
(32)【優先日】2020-12-24
(33)【優先権主張国・地域又は機関】JP
(31)【優先権主張番号】P 2021169911
(32)【優先日】2021-10-15
(33)【優先権主張国・地域又は機関】JP
(71)【出願人】
【識別番号】000006747
【氏名又は名称】株式会社リコー
(71)【出願人】
【識別番号】504179255
【氏名又は名称】国立大学法人 東京医科歯科大学
(74)【代理人】
【識別番号】100107766
【弁理士】
【氏名又は名称】伊東 忠重
(74)【代理人】
【識別番号】100070150
【弁理士】
【氏名又は名称】伊東 忠彦
(72)【発明者】
【氏名】星 文和
(72)【発明者】
【氏名】宮野 由貴
(72)【発明者】
【氏名】山形 秀明
(72)【発明者】
【氏名】川端 茂▲徳▼
(72)【発明者】
【氏名】長岡 信頼
(72)【発明者】
【氏名】渡部 泰士
【テーマコード（参考）】
4C127
【Ｆターム（参考）】
4C127AA10
4C127DD03
4C127FF07
4C127GG13
(57)【要約】
【課題】複数の部位に対する刺激に紐づいたトリガー信号に応じて計測される複数の生体データを加算平均処理する場合に、加算平均処理の効果を部位ごとで等しく得ること。
【解決手段】本発明の一態様に係る生体データ処理装置は、少なくとも１つの部位に対する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する加算平均処理部と、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存する保存部と、指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う生体データ処理部と、を有する。
【選択図】図７
【解決手段】本発明の一態様に係る生体データ処理装置は、少なくとも１つの部位に対する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する加算平均処理部と、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存する保存部と、指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う生体データ処理部と、を有する。
【選択図】図７
【特許請求の範囲】
【請求項1】
少なくとも１つの部位に対する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する加算平均処理部と、
少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存する保存部と、
指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う生体データ処理部と、を有する生体データ処理装置。
少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存する保存部と、
指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う生体データ処理部と、を有する生体データ処理装置。
【請求項2】
前記生体データは、生体が発生する磁場データであり、
前記生体データ処理部は、前記磁場データに基づき、前記生体における少なくとも１つの刺激部位ごとでの活動電流の強さを推定する処理を行う請求項１に記載の生体データ処理装置。
前記生体データ処理部は、前記磁場データに基づき、前記生体における少なくとも１つの刺激部位ごとでの活動電流の強さを推定する処理を行う請求項１に記載の生体データ処理装置。
【請求項3】
前記加算回数の指定を受け付ける第１指定受付部を有する請求項１又は２に記載の生体データ処理装置。
【請求項4】
前記所定回数の指定を受け付ける第２指定受付部をさらに有する請求項３に記載の生体データ処理装置。
【請求項5】
更新間隔回数の指定を受け付ける第３指定受付部をさらに有し、
前記所定回数は、前記加算回数が前記所定回数に達するたびに、前記第３指定受付部により受け付けられた前記更新間隔回数が前記所定回数に加算された回数に自動更新される請求項３又は４に記載の生体データ処理装置。
前記所定回数は、前記加算回数が前記所定回数に達するたびに、前記第３指定受付部により受け付けられた前記更新間隔回数が前記所定回数に加算された回数に自動更新される請求項３又は４に記載の生体データ処理装置。
【請求項6】
前記加算回数に応じて前記生体データ処理部に処理方法を切り替えさせる切替部を有する請求項１乃至５の何れか１項に記載の生体データ処理装置。
【請求項7】
前記切替部は、前記加算回数が総加算回数に対応する加算回数である場合にのみ前記処理方法を切り替え、
前記生体データ処理部は、前記加算回数が総加算回数に対応する加算回数以外である場合には、総加算回数に対応する加算回数の場合と比較して高速な処理を行い、前記加算回数が総加算回数に対応する加算回数である場合には、前記加算回数が総加算回数に対応する加算回数以外である場合と比較して高精度な処理を行う請求項６に記載の生体データ処理装置。
前記生体データ処理部は、前記加算回数が総加算回数に対応する加算回数以外である場合には、総加算回数に対応する加算回数の場合と比較して高速な処理を行い、前記加算回数が総加算回数に対応する加算回数である場合には、前記加算回数が総加算回数に対応する加算回数以外である場合と比較して高精度な処理を行う請求項６に記載の生体データ処理装置。
【請求項8】
請求項１乃至７の何れか１項に記載の生体データ処理装置と、
生体を計測する計測装置と、を有する生体データ計測システム。
生体を計測する計測装置と、を有する生体データ計測システム。
【請求項9】
前記生体データ処理部による処理結果に基づき、前記計測装置による計測の中断又は延長を制御する計測制御部をさらに有する請求項８に記載の生体データ計測システム。
【請求項10】
前記生体データ処理部による前記処理結果を表示する表示部と、
前記計測装置による計測の中断又は延長の何れか一方の指示を受け付ける指示受付部と、を有し、
前記計測制御部は、前記指示に応じて前記計測装置による計測を中断又は延長させる請求項９に記載の生体データ計測システム。
前記計測装置による計測の中断又は延長の何れか一方の指示を受け付ける指示受付部と、を有し、
前記計測制御部は、前記指示に応じて前記計測装置による計測を中断又は延長させる請求項９に記載の生体データ計測システム。
【請求項11】
少なくとも１つの部位に対する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理し、
少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存部が保存し、
指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う
処理をコンピュータに実行させるプログラム。
少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存部が保存し、
指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う
処理をコンピュータに実行させるプログラム。
【請求項12】
トリガー信号に応じた生体データの加算回数が所定回数に達するたびに磁場データをセグメント加算平均処理する加算平均処理部と、
複数のセグメント加算平均処理結果であるセグメント加算データと、前記セグメント加算データにおける加算回数と、を対応付けて保存する保存部と、
前記セグメント加算データを表示させる表示部と、
前記保存部を参照して取得される前記セグメント加算データのうち、任意の前記セグメント加算データを加算平均処理して加算平均データを取得する加算平均処理部と、
を有する生体データ処理装置。
複数のセグメント加算平均処理結果であるセグメント加算データと、前記セグメント加算データにおける加算回数と、を対応付けて保存する保存部と、
前記セグメント加算データを表示させる表示部と、
前記保存部を参照して取得される前記セグメント加算データのうち、任意の前記セグメント加算データを加算平均処理して加算平均データを取得する加算平均処理部と、
を有する生体データ処理装置。
【請求項13】
前記セグメント加算データ及び前記加算平均データの少なくとも一方を対象として、前記生体データに基づく処理を行う生体データ処理部を有する請求項１２に記載の生体データ処理装置。
【請求項14】
表示部を有し、
前記生体データは、生体が発生する磁場データであり、
前記生体データ処理部は、前記磁場データに基づく体内の活動電流の強さを推定する処理及び前記磁場データに周波数解析を行う処理の少なくとも一つを行い、その結果を前記表示部に表示させる請求項１３に記載の生体データ処理装置。
前記生体データは、生体が発生する磁場データであり、
前記生体データ処理部は、前記磁場データに基づく体内の活動電流の強さを推定する処理及び前記磁場データに周波数解析を行う処理の少なくとも一つを行い、その結果を前記表示部に表示させる請求項１３に記載の生体データ処理装置。
【請求項15】
前記表示部は、複数の前記セグメント加算データ及び複数の前記加算平均データを並べて表示する請求項１４に記載の生体データ処理装置。
【請求項16】
トリガー信号に応じた生体データの加算回数が所定回数に達するたびに磁場データをセグメント加算平均処理する加算平均処理部と、
セグメント加算データを表示させる表示部と、
トリガー信号に対応付けられた全エポックデータを保存する保存部と、
前記保存部を参照して取得される全エポックデータのうち、任意のエポックデータを選択して加算平均処理し、加算平均データを取得する加算平均処理部と、を有する生体データ処理装置。
セグメント加算データを表示させる表示部と、
トリガー信号に対応付けられた全エポックデータを保存する保存部と、
前記保存部を参照して取得される全エポックデータのうち、任意のエポックデータを選択して加算平均処理し、加算平均データを取得する加算平均処理部と、を有する生体データ処理装置。
【請求項17】
前記所定回数の設定とセグメント幅の設定は、別々に行われる請求項１６に記載の生体データ処理装置。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、生体データ処理装置、生体データ計測システム、及びプログラムに関する。
【背景技術】
【0002】
従来、電気刺激等の刺激に対して発生する生体の磁場データ等の生体データを計測する生体データ計測システムが知られている。このような生体データ計測システムでは、微弱な生体データに含まれるノイズを低減するために、定期的に発生するトリガー信号に応じて計測された生体データが加算平均処理される場合がある。
【0003】
また、加算平均処理を伴う心電図検査の利便性を向上させるために、所定のテンプレートに合致する心電図信号を加算対象とし、合致しないものを加算対象外とするように加算平均処理を実行する構成が開示されている（例えば、特許文献１参照）。
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、複数のトリガー信号に応じて計測される生体データを加算平均処理する場合には、前記トリガー信号に対応する刺激部位ごとでの加算回数の違いにより、加算平均処理の効果を刺激部位ごとで等しく得ることができない場合がある。特許文献１の構成は１つのトリガー信号に応じた生体データの加算平均処理を開示するものであるため、このような課題を解決することはできない。
【0005】
本発明は、複数の部位に対する刺激に紐づいたトリガー信号に応じて計測される複数の生体データを加算平均処理する場合に、加算平均処理の効果を刺激部位ごとで等しく得ることを課題とする。
【課題を解決するための手段】
【0006】
本発明の一態様に係る生体データ処理装置は、少なくとも１つの部位に対する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する加算平均処理部と、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データにおける前記加算回数と、を対応付けて保存する保存部と、指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う生体データ処理部と、を有する。
【発明の効果】
【0007】
本発明によれば、複数の部位に対する刺激に紐づいたトリガー信号に応じて計測される複数の生体データを加算平均処理する場合に、加算平均処理の効果を前記トリガー信号に対応する刺激部位ごとで等しく得ることができる。
【図面の簡単な説明】
【0008】
【図1】実施形態に係る生体データ計測システムの全体構成例を示す図である。
【図2】実施形態に係る計測装置の構成例を説明する図である。
【図3】実施形態に係るコンピュータのハードウェア構成例のブロック図である。
【図4】第１実施形態に係る計測ＷＳの機能構成例のブロック図である。
【図5】第１実施形態に係るデータ保存サーバの機能構成例のブロック図である。
【図6】第１実施形態に係る解析ＷＳの機能構成例のブロック図である。
【図7】第１実施形態に係る計測ＷＳの動作例のフローチャートである。
【図8】第１実施形態に係る解析ＷＳの動作例のフローチャートである。
【図9】所定回数及び総加算回数の指定画面例を示す図である。
【図10】計測画面及び操作画面例を示す図である。
【図11】加算平均データリスト例を示す図である。
【図12】活動電流推定中の表示画面例を示す図である。
【図13】適切でない活動電流推定結果の表示画面例を示す図である。
【図14】適切な活動電流の推定結果の表示画面例を示す図である。
【図15】第２実施形態に係る解析ＷＳの機能構成例のブロック図である。
【図16】第３実施形態に係る計測ＷＳの機能構成例のブロック図である。
【図17】第３実施形態に係る所定回数及び総加算回数の指定画面の第１例の図である。
【図18】第３実施形態に係る所定回数及び総加算回数の指定画面の第２例の図である。
【図19】第４～６実施形態に係る生体データ計測システムの全体構成例の図である。
【図20】第４実施形態に係る計測ＷＳの機能構成図のブロック図である。
【図21】第４実施形態に係るデータ保存サーバの機能構成図のブロック図である。
【図22】第４実施形態に係る解析ＷＳの機能構成図のブロック図である。
【図23】第４実施形態に係る生体データ処理装置の動作例のフローチャートである。
【図24】第４実施形態に係る所定回数、セグメント幅及び総加算回数の指定画面の第1例の図である。
【図25】第４実施形態に係る所定回数、セグメント幅及び総加算回数の指定画面の第２例の図である。
【図26】計測画面または解析画面上に表示された波形データ例の図である。
【図27】計測画面または解析画面上に表示された周波数スペクトル例の図である。
【図28】セグメント加算データリスト例を示す図である。
【図29】第５実施形態に係る計測ＷＳの機能構成図のブロック図である。
【図30】第５実施形態に係るデータ保存サーバの機能構成図のブロック図である。
【図31】第５実施形態に係る解析ＷＳの機能構成図のブロック図である。
【図32】第５実施形態に係る生体データ処理装置の動作例のフローチャートである。
【図33】第６実施形態に係る所定回数、セグメント幅及び総加算回数の指定画面を示す図である。
【図34】第６実施形態に係る生体データ処理装置の動作例のフローチャートである。
【発明を実施するための形態】
【0009】
以下、図面を参照して発明を実施するための形態について説明する。各図面において、同一の構成部分には同一符号を付し、重複した説明を省略する場合がある。
【0010】
また以下に示す実施形態は、本発明の技術思想を具体化するための生体データ処理装置を例示するものであって、本発明を以下に示す実施形態に限定するものではない。以下に記載されている構成部品の形状、その相対的配置、パラメータの値等は特定的な記載がない限り、本発明の範囲をそれのみに限定する趣旨ではなく、例示することを意図したものである。また図面が示す部材の大きさや位置関係等は、説明を明確にするため、誇張している場合がある。
【0011】
実施形態に係る生体データ処理装置は、計測装置で計測された生体データに基づく処理を行う装置である。生体データに基づく処理は、例えば生体の活動電流を推定する処理である。
【0012】
計測装置は、電気刺激等の刺激に対して発生する生体の磁場データを生体データとして計測する脊磁計等の計測装置である。また脊磁計とは、脊髄の神経活動によって生じるわずかな磁場を計測し、身体を傷つけることなく神経活動を可視化する計測装置をいう（例えば、「https://reader.elsevier.com/reader/sd/pii/S1388245718313415?token=C588DF7ADF1D16EAB0CF51EEF358B69D0570B5948A51FB4CFE60D5406A1D8A4DF9D54238F814CBD7DF87EFFF38BBDBC0」参照）。
【0013】
生体の活動電流とは、生体の細胞や組織が刺激を受けたときに発生する活動電位により、刺激を受けた部分が他の部分に対して負の電位をもつことで、電位差が生じて流れる微弱な電流をいう。
【0014】
実施形態では、複数のトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理された結果である加算平均データと、該加算平均データにおける加算回数とを対応付けて保存する。
【0015】
そして、指定された加算回数に基づき保存部を参照して取得されるトリガー信号に対応した加算平均データを用いて、生体データに基づく処理を行う。生体データに基づく処理は、例えば少なくとも１つの刺激部位ごとでの活動電流の強さを推定する処理である。
【0016】
実施形態では、少なくとも１つの刺激部位ごとで同じ加算回数による加算平均データを用いることで、刺激部位ごとでの加算回数の違いをなくし、刺激部位ごとで加算平均処理の効果を等しく得ることを可能にする。
【0017】
以下では、生体の磁場を計測する計測装置と、該計測装置により計測された生体データに対して、該生体データから生体の活動電流を推定する生体データ処理装置とを有する生体データ計測システムを一例として実施形態を説明する。また、実施形態では、生体に電気刺激を与えることで、生体における脊髄内の神経を流れる電流の強さを推定する例について説明する。
【0018】
なお、上記及び以下に示す「ユーザ」は、生体データ計測システムを使用するユーザである。より具体的には、生体データ計測システムを使用して生体データを取得する技師、或いは診察又は診断を行う医師である。
【0019】
【0020】
図１は、生体データ計測システム１の全体構成の一例を説明する図である。図１に示すように生体データ計測システム１は、計測装置２と、計測ＷＳ（Work Station）３と、解析ＷＳ４と、データ保存サーバ５とを有する。これらは有線又は無線によって、相互に通信可能に接続されている。これらのうち、計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５は、生体データ処理装置１０を構成する。
【0021】
計測装置２は、複数のトリガー信号のそれぞれに応じた電気刺激等の刺激に対して生体で発生する磁場データを計測する脊磁計である。磁場データは計測データの一例である。計測装置２は複数の部位に対する刺激に紐づいたトリガー信号に対する計測結果である磁場データを、複数のトリガー信号とともに計測ＷＳ３に送信する。
【0022】
計測ＷＳ３は、計測装置２から受信した複数のトリガー信号をカウントして複数の刺激部位ごとでの加算回数を取得し、複数の刺激部位ごとでの加算回数が所定回数に達するたびに磁場データを加算平均処理する。そして、加算平均処理結果である加算平均データと、該加算平均データにおける加算回数の情報とを対応付けてデータ保存サーバ５に送信する。
【0023】
データ保存サーバ５は、計測ＷＳ３から受信した加算平均データと加算回数とを対応付けて保存する。
【0024】
解析ＷＳ４は、ユーザが指定した加算回数に基づいて、データ保存サーバ５を参照して複数の刺激部位ごとでの加算平均データを取得し、取得された加算平均データを用いて、複数の刺激部位ごとでの活動電流の強さを推定する。解析ＷＳ４は、解析ＷＳ４が有するディスプレイに推定結果を表示したり、データ保存サーバ５に送信して保存させたり、外部サーバ等の外部装置に送信したりすることができる。
【0025】
なお、本実施形態では、計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５の３台の装置により生体データ処理装置１０を構成する構成を例示するが、これに限定されるものではない。計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５の機能を集約した１台の装置で生体データ処理装置１０を構成してもよいし、計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５の機能を分散させた４台以上の装置で生体データ処理装置１０を構成してもよい。
【0026】
また生体データ計測システム１は、計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５以外の装置を通信可能に備えてもよいし、計測装置２以外の他の生体データの計測装置を通信可能に備えてもよい。
【0027】
＜計測装置２の構成例＞
次に図２を参照して、計測装置２の構成について説明する。
次に図２を参照して、計測装置２の構成について説明する。
【0028】
【0029】
磁気センサアレイ２００は、複数の磁気センサ２０１がアレイ状に配置された生体センサであり、被検者１００の首の背後に設置されている。ここで、被検者１００は、「生体」の一例である。
【0030】
複数の磁気センサ２０１のそれぞれは、図２に矢印で示したｘ軸、ｙ軸、及びｚ軸の各方向における生体の磁場を計測し、磁場データを出力する。図２の例では、磁気センサアレイ２００は、７×５個の磁気センサを含んでおり、複数の磁気センサ２０１のそれぞれが計測した磁場データは、生体データ処理装置１０に出力される。なお、被検者１００に対して磁気センサアレイ２００を設置する位置は、マーカーコイル等を用いて予め調整されている。
【0031】
また、デュワー２１０の内部は、液体ヘリウムが充填されており、磁気センサアレイ２００を極低温で動作させるための冷却が行われている。
【0032】
なお、実施形態では、磁気センサアレイ２００上の点２４０の位置を、ｘ軸、ｙ軸及びｚ軸の原点とする。磁気センサアレイ２００上の点２４０の位置をｘ軸、ｙ軸、及びｚ軸の原点とすることで、磁気センサアレイ２００における複数の磁気センサ２０１間の相対的な位置関係を、全てｘ座標、ｙ座標及びｚ座標により表すことができる。
【0033】
また、計測装置２による磁場の計測方法には、特開２０１８―０８９１０４号公報等に記載された公知技術を適用できるため、ここではさらに詳細な説明を省略する。
【0034】
さらに、図２では、磁気センサアレイ２００が被検者１００の首の背後に設置された例を示したが、以下では、磁気センサアレイ２００が被検者１００の腰の背後に設置され、腰付近の生体内部の磁場データが計測されて、脊髄内の神経を流れる電流の強さを推定する例を説明する。但し、推定対象は、脊髄内の神経を流れる電流の強さに限定されるものではない。例えば腕や足等の四肢の末梢神経を流れる電流の強さを推定対象にすることもできる。
【0035】
＜コンピュータのハードウェア構成例＞
本実施形態における計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５はそれぞれコンピュータにより構築可能である。図３を参照して、このコンピュータのハードウェア構成について説明する。
本実施形態における計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５はそれぞれコンピュータにより構築可能である。図３を参照して、このコンピュータのハードウェア構成について説明する。
【0036】
図３は、コンピュータのハードウェア構成の一例を説明するブロック図である。なお、図３では計測ＷＳ３を構築するコンピュータのハードウェア構成を例示するが、解析ＷＳ４及びデータ保存サーバ５のそれぞれを構築するコンピュータのハードウェア構成も図３と同様である。
【0037】
図３に示すように、計測ＷＳ３は、ＣＰＵ（Central Processing Unit）５０１と、ＲＯＭ（Read Only Memory）５０２と、ＲＡＭ（Random Access Memory）５０３と、ＨＤ（Hard Disk）５０４と、ＨＤＤ（Hard Disk Drive）コントローラ５０５と、ディスプレイ５０６と、外部機器接続Ｉ／Ｆ（Interface）５０８と、ネットワークＩ／Ｆ５０９とを有する。
【0038】
また、計測ＷＳ３は、データバス５１０と、キーボード５１１と、ポインティングデバイス５１２と、ＤＶＤ－ＲＷ（Digital Versatile Disk Rewritable）ドライブ５１４と、メディアＩ／Ｆ５１６とを有する。
【0039】
これらのうち、ＣＰＵ５０１は、計測ＷＳ３全体の動作を制御する。ＲＯＭ５０２は、ＩＰＬ（Initial Program Loader）等のＣＰＵ５０１の駆動に用いられるプログラムを記憶する。
【0040】
ＲＡＭ５０３は、ＣＰＵ５０１のワークエリアとして使用される。ＨＤ５０４は、プログラム等の各種データを記憶する。ＨＤＤコントローラ５０５は、ＣＰＵ５０１の制御に従って、ＨＤ５０４に対する各種データの読み出し又は書き込みを制御する。
【0041】
ディスプレイ５０６は、カーソル、メニュー、ウィンドウ、文字、又は画像などの各種情報を表示する。外部機器接続Ｉ／Ｆ５０８は、各種の外部機器を接続するためのインターフェースである。この場合の外部機器は、例えば、ＵＳＢ（Universal Serial Bus）メモリやプリンタ等である。
【0042】
ネットワークＩ／Ｆ５０９は、ネットワークを利用してデータ通信をするためのインターフェースである。データバス５１０は、ＣＰＵ５０１等の各構成要素を電気的に接続するためのアドレスバスやデータバス等である。
【0043】
キーボード５１１は、文字、数値、各種指示などの入力のための複数のキーを備えた入力手段の一種である。ポインティングデバイス５１２は、各種指示の選択や実行、処理対象の選択、カーソルの移動などを行う入力手段の一種である。
【0044】
ＤＶＤ－ＲＷドライブ５１４は、着脱可能な記録媒体の一例としてのＤＶＤ－ＲＷ５１３に対する各種データの読み出し又は書き込みを制御する。なお、ＤＶＤ－ＲＷに限らず、ＤＶＤ－Ｒ等であってもよい。メディアＩ／Ｆ５１６は、フラッシュメモリ等の記録メディア５１５に対するデータの読み出し又は書き込み（記憶）を制御する。
【0045】
［第１実施形態］
＜生体データ処理装置１０の機能構成例＞
次に図４乃至図６を参照して、生体データ処理装置１０を構成する計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５のそれぞれの機能構成について説明する。
＜生体データ処理装置１０の機能構成例＞
次に図４乃至図６を参照して、生体データ処理装置１０を構成する計測ＷＳ３、解析ＷＳ４及びデータ保存サーバ５のそれぞれの機能構成について説明する。
【0046】
【0047】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することで実現される機能又は機能する手段である。なお図４は、計測ＷＳ３が有する主な構成を示すが、計測ＷＳ３はこれら以外の構成を有してもよい。例えば計測ＷＳ３は計測装置２から受信した磁場データを示す波形データを表示する表示部を備えてもよい。
【0048】
通信部３１は、計測装置２、解析ＷＳ４及びデータ保存サーバ５のそれぞれとの間でデータ及び信号の送受信を行う。
【0049】
加算平均処理部３２は、ユーザがキーボード５１１（図３参照）等を用いて入力した所定回数及び総加算回数の情報を取得する。総加算回数は、加算平均処理部３２が磁場データの加算を行う総回数である。なお、加算平均処理部３２は、ＨＤ５０４等に予め保存された所定回数及び総加算回数の情報を、ＨＤ５０４等を参照して取得してもよい。
【0050】
加算平均処理部３２は、計測装置２から通信部３１を介して、複数のトリガー信号を受信する。また加算平均処理部３２は、複数のトリガー信号のそれぞれに対して、磁場計測が開始された時刻から所定の周期で計測装置２が計測した磁場データを、通信部３１を介して受信して加算処理する。なお、所定の周期はトリガー信号ごとで異なる場合がある。
【0051】
加算平均処理部３２は、受信した複数のトリガー信号をカウントして各トリガー信号に対応する刺激部位ごとでの加算回数を取得し、加算回数が所定回数に達するたびに、各トリガー信号に対応した磁場データを加算平均処理する。各トリガー信号に対応する刺激部位ごとで加算回数が異なる場合には、同じ時刻で加算平均処理する刺激部位と加算平均処理しない刺激部位が生じる。
【0052】
加算平均処理は、計測装置２が計測した磁場データを逐次加算処理した値を加算回数で除算して平均値を算出する処理である。加算平均処理部３２は、処理結果である加算平均データと、加算平均データにおける加算回数の情報とを対応付けて、通信部３１を介してデータ保存サーバ５に送信する。
【0053】
計測制御部３３は、解析ＷＳ４が有する推定部（後述）による活動電流の強さの推定結果に基づく指示を、通信部３１を介して受信し、この推定結果に基づき、計測装置２に計測を中断させ、又は延長させることができる。計測制御部３３は、解析ＷＳ４からの計測の中断又延長の指示を、指示があった任意のタイミングで割込みデータとして受信可能である。
【0054】
【0055】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することで実現される機能又は機能する手段である。なお図５は、データ保存サーバ５が有する主な構成を示すが、データ保存サーバ５はこれら以外の構成を有してもよい。
【0056】
通信部５１は、計測ＷＳ３及び解析ＷＳ４のそれぞれとの間でデータ及び信号の送受信を行う。
【0057】
保存部５２は、通信部５１を介して計測ＷＳ３から受信した加算平均データ５２２と加算平均で用いた加算回数５２１とを対応付けて保存する。なお、加算平均データ５２２は、複数の加算平均データの総称表記であり、加算回数５２１は、複数の加算回数の総称表記であり、パラメータ５２３は、複数のパラメータの総称表記である。
【0058】
（解析ＷＳ４の機能構成例）
次に図６は、解析ＷＳ４の機能構成の一例を説明するブロック図である。図６に示すように、解析ＷＳ４は、通信部４１と、推定部４２と、第１指定受付部４３と、指示受付部４４と、表示部４５と、判断部４６とを有する。
次に図６は、解析ＷＳ４の機能構成の一例を説明するブロック図である。図６に示すように、解析ＷＳ４は、通信部４１と、推定部４２と、第１指定受付部４３と、指示受付部４４と、表示部４５と、判断部４６とを有する。
【0059】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することで実現される機能又は機能する手段である。なお図６は、解析ＷＳ４が有する主な構成を示すが、解析ＷＳ４はこれら以外の構成を有してもよい。
【0060】
通信部４１は、計測ＷＳ３及びデータ保存サーバ５のそれぞれとの間でデータ及び信号の送受信を行う。
【0061】
第１指定受付部４３は、ユーザがキーボード５１１等を用いて指定した加算回数の情報を受け付ける。例えば、第１指定受付部４３は、通信部４１を介して保存部５２が保存する加算平均データのリストを取得し、取得した加算平均データのリストを表示部４５によりディスプレイ５０６等に表示させる。第１指定受付部４３は、ユーザが加算平均データのリストを視認して選択した結果から加算回数の情報を受け付けることができる。
【0062】
推定部４２は、生体データに基づく処理を行う生体データ処理部の一例である。推定部４２は、計測装置２が計測した磁場データに基づき、生体の活動電流の強さを推定する処理を行う。
【0063】
具体的には、推定部４２は、指定された加算回数に基づき、通信部４１を介してデータ保存サーバ５の保存部５２を参照して複数の刺激部位ごとでの加算平均データを取得する。保存部５２は、生体の複数の刺激部位ごとで、加算回数と加算平均データを対応付けて保存しているため、加算回数を指定することで、加算回数に対応する加算平均データを複数の刺激部位ごとで取得できる。
【0064】
推定部４２は、取得した加算平均データを用いて、複数の刺激部位ごとでの活動電流の強さを推定する。この推定アルゴリズムには、「Array-Gain Constraint Minimum-Norm Spatial Filter With Recursively Updated Gram Matrix」（例えば、「https://ieeexplore.ieee.org/document/5415622」参照）等を用いることができる。
【0065】
ここで、計測装置２による磁場データの計測の周期は、生体の複数の刺激部位ごとで異なる場合がある。あるいは、特定の刺激部位のトリガーにのみ、トリガーの発生または測定に不具合が生じ、一部のデータが欠落する場合もある。そのため、例えば計測時刻に基づき刺激部位ごとでの加算平均データを取得すると、刺激部位ごとで加算回数が異なる場合がある。加算回数が異なると、加算平均処理の効果を刺激部位ごとで等しく得ることができず、活動電流の強さの推定精度が刺激部位ごとで異なるものとなる。
【0066】
これに対し、本実施形態では、刺激部位ごとで同じ加算回数に対応する加算平均データを用いて活動電流の強さの推定するため、加算平均処理の効果を刺激部位ごとで等しく得ることができる。これにより、活動電流の強さの推定精度が刺激部位ごとで等しくなる。
【0067】
表示部４５は、推定部４２による活動電流の強さの推定結果を表示する。例えば表示部４５は、ディスプレイ５０６に推定結果を表示し、ユーザに視認させることができる。また、表示部４５は、計測装置２が計測した磁場データを示す波形データを、計測ＷＳ３を介して受信し、これを表示することもできる。
【0068】
推定部４２は、測定途中であっても、すでに保存されている（途中の加算回数の）データについて活動電流の強さを推定し、ユーザは、表示部４５が表示した推定部４２による推定結果を視認できる。
【0069】
ユーザは、ディスプレイ５０６に表示された活動電流の強さの推定結果を視認し、刺激部位ごとでの推定結果が妥当であるか否か、或いは十分な数量の磁場データを計測できたか否か等を判断できる。
【0070】
ユーザは活動電流の推定結果が妥当でないと判断した場合には、キーボード５１１等を用いて計測の中断を指示する。指示受付部４４は、受け付けた中断の指示を、通信部４１を介して計測ＷＳ３に送信する。計測ＷＳ３の計測制御部３３は、指示に応答して計測装置２に計測を中断させることができる。
【0071】
また、ユーザは磁場データの数量が不足していると判断した場合には、キーボード５１１等を用いて計測の延長を指示する。指示受付部４４は、受け付けた延長の指示を、通信部４１を介して計測ＷＳ３に送信する。計測ＷＳ３の計測制御部３３は、指示に応答して計測装置２に計測を延長させることができる。
【0072】
判断部４６は、推定部４２による推定結果が総加算回数に対応するものであるか否かを判断し、総加算回数に対応する加算平均データがあるか否かを判断する。また、計測を延長するか否かを判断する。そして、判断部４６はこれらの判断結果に応じて、計測の中断、又は延長の指示を示す信号を指示受付部４４に出力する。
【0073】
指示受付部４４は、ユーザの判断による計測の中断、又は延長の何れか一方の指示だけでなく、判断部４６による計測の中断、又は延長の何れか一方の指示を受け付け、通信部４１を介して計測ＷＳ３に送信することができる。
【0074】
【0075】
【0076】
まず、ステップＳ７１において、加算平均処理部３２は、ユーザがキーボード５１１等を用いて入力した所定回数及び総加算回数の情報を取得する。なお、加算平均処理部３２は、ＨＤ５０４等に予め保存された所定回数及び総加算回数の情報をＨＤ５０４等から取得してもよい。
【0077】
続いて、ステップＳ７２において、加算平均処理部３２は、複数のトリガー信号と、複数のトリガー信号のそれぞれに対応する複数の刺激部位ごとで計測された磁場データを、計測装置２から通信部３１を介して受信して加算処理する。
【0078】
続いて、ステップＳ７３において、加算平均処理部３２は、受信した複数のトリガー信号をカウントして複数の刺激部位ごとでの加算回数を取得し、複数の刺激部位ごとでの加算回数が所定回数に達したか否かを判定する。
【0079】
ステップＳ７３で達したと判定された場合には（ステップＳ７３、Ｙｅｓ）、ステップＳ７４において、加算平均処理部３２は、複数の刺激部位ごとの磁場データを加算平均処理する。一方、達していないと判定された場合には（ステップＳ７３、Ｎｏ）、ステップＳ７２以降の動作が再度行われる。
【0080】
続いて、ステップＳ７４において、加算平均処理部３２は、加算平均処理結果である加算平均データと、加算平均データにおける加算回数の情報とを対応付けて、通信部３１を介してデータ保存サーバ５に送信する。データ保存サーバ５は、受信した加算平均データと、加算回数の情報とを対応付けて保存できる。
【0081】
続いて、ステップＳ７６において、計測制御部３３は、計測の中断の指示があったか否かを判定する。
【0082】
ステップＳ７６で、指示があったと判定された場合には（ステップＳ７６、Ｙｅｓ）、動作はステップＳ７９に移行する。一方、指示がなかったと判定された場合には（ステップＳ７６、Ｎｏ）、動作はステップＳ７７に移行する。
【0083】
なお、ステップＳ７６の動作は、解析ＷＳ４からの割込みデータに基づくものであり、任意のタイミングで行われる。従って、ステップＳ７１乃至Ｓ７９のどの順番でステップＳ７６の動作が行われてもよい。
【0084】
続いて、ステップＳ７７において、計測制御部３３は、加算回数は総加算回数に達したか否かを判定する。なお、この判定は、計測制御部３３でなく、加算平均処理部３２が行ってもよい。
【0085】
ステップＳ７７で、達していないと判定された場合には（ステップＳ７７、Ｎｏ）、ステップＳ７２以降の動作が再度行われる。一方、ステップＳ７７で、達したと判定された場合には（ステップＳ７７、Ｙｅｓ）、ステップＳ７８において、計測制御部３３は、計測の延長の指示があったか否かを判定する。
【0086】
ステップＳ７８で、指示があったと判定された場合には（ステップＳ７８、Ｙｅｓ）、ステップＳ７２以降の動作が再度行われる。一方、指示がなかったと判定された場合には（ステップＳ７８、Ｎｏ）、ステップＳ７９において、計測制御部３３は、計測装置２に計測を終了させる。
【0087】
なお、ステップＳ７８の動作は、解析ＷＳ４からの割込みデータに基づくものであり、任意のタイミングで行われる。従って、ステップＳ７１乃至Ｓ７９のどの順番でステップＳ７８の動作が行われてもよい。
【0088】
このようにして、計測ＷＳ３は加算平均処理を実行し、また計測の中断又は延長の指示に応じた計測装置２の制御を行うことができる。
【0089】
（解析ＷＳ４の動作例）
次に図８は、解析ＷＳ４の動作の一例を示すフローチャートである。図８は、ユーザが解析ＷＳ４を操作して加算平均データを用いて生体の活動電流の推定処理を開始した時点をトリガーにした解析ＷＳ４の動作を示している。
次に図８は、解析ＷＳ４の動作の一例を示すフローチャートである。図８は、ユーザが解析ＷＳ４を操作して加算平均データを用いて生体の活動電流の推定処理を開始した時点をトリガーにした解析ＷＳ４の動作を示している。
【0090】
まず、ステップＳ８１において、第１指定受付部４３は、通信部４１を介して保存部５２が保存する加算平均データのリストを取得し、取得した加算平均データのリストを表示部４５によりディスプレイ５０６等に表示させる。ユーザは、表示された加算平均データのリストを視認し、キーボード５１１等を用いて加算平均データを選択することができる。
【0091】
続いて、ステップＳ８２において、第１指定受付部４３はユーザが選択した結果から加算回数の情報を受け付ける。
【0092】
続いて、ステップＳ８３において、推定部４２は、第１指定受付部４３が受け付けた加算回数に基づき、通信部４１を介してデータ保存サーバ５の保存部５２を参照して複数の刺激部位ごとでの加算平均データを取得する。
【0093】
続いて、ステップＳ８４において、推定部４２は、取得した加算平均データ用いて、複数の刺激部位ごとでの活動電流の強さを推定する。
【0094】
続いて、ステップＳ８５において、表示部４５は、推定部４２による活動電流の強さの推定結果を表示する。例えば表示部４５は、ディスプレイ５０６に推定結果を表示し、ユーザに視認させる。推定部４２は、測定途中であっても、すでに保存されている（途中の加算回数の）データについて活動電流の強さを推定し、ユーザは、表示部４５が表示した推定部４２による推定結果を視認できる。
【0095】
続いて、ステップ８６において、指示受付部４４は、活動電流の強さの推定結果を視認したユーザによる、推定結果が妥当であるか否かの判断結果を受け付ける。
【0096】
ステップＳ８６で、妥当であるという判断結果を受け付けた場合には（ステップＳ８６、Ｙｅｓ）、解析ＷＳ４は動作を終了する。一方、妥当でないという判断結果を受け付けた場合には（ステップＳ８６、Ｎｏ）、ステップＳ８７において、判断部４６は、推定結果は総加算回数に対応する加算平均データによるものであるか否かを判断する。ここで、妥当であるか否かの判断には、測定途中のデータを用いて推定された活動電流の強さを視認して判断する場合には、測定途中のデータとして妥当であるか否かの判断も含まれる。
【0097】
ステップＳ８７において、総加算回数に対応する加算平均データによるものでないと判断された場合には（ステップＳ８７、Ｎｏ）、ステップＳ８８において、判断部４６は、総加算回数に対応する加算平均データがあるか否かを判断する。
【0098】
ステップＳ８８で、総加算回数に対応する加算平均データがあると判断された場合には（ステップＳ８８、Ｙｅｓ）、解析ＷＳ４は動作を終了する。一方、総加算回数に対応する加算平均データがないと判断された場合には（ステップＳ８８、Ｎｏ）、ステップＳ８９において、指示受付部４４は判断部４６から計測の中断の指示を受け付け、通信部４１を介して計測ＷＳ３に送信する。その後、解析ＷＳ４は動作を終了する。
【0099】
一方で、ステップＳ８７で、総加算回数に対応する加算平均データによるものであると判断された場合には（ステップＳ８７、Ｙｅｓ）、ステップＳ９０において、判断部４６は、計測を延長するか否かを判断する。
【0100】
ステップＳ９０で、計測を延長しないと判断された場合には（ステップＳ９０、Ｎｏ）、解析ＷＳ４は動作を終了する。一方、計測を延長すると判断された場合には（ステップＳ９０、Ｎｏ）、ステップＳ９１において、指示受付部４４は判断部４６から計測の延長の指示を受け付け、通信部４１を介して計測ＷＳ３に送信する。その後、解析ＷＳ４は動作を終了する。
【0101】
このようにして、解析ＷＳ４は、生体の活動電流の推定処理を実行し、また推定結果に基づく計測の中断又は延長の指示を計測ＷＳ３に対して行うことができる。
【0102】
＜各種表示画面例＞
次に、生体データ計測システム１が表示する各種表示画面について説明する。
次に、生体データ計測システム１が表示する各種表示画面について説明する。
【0103】
（所定回数及び総加算回数の指定画面例）
まず図９は、計測ＷＳ３が表示する所定回数及び総加算回数の指定画面の一例を示す図である。計測ＷＳ３は、計測装置２による磁場データの計測を開始する際に、図９の画面をディスプレイに表示させる。
まず図９は、計測ＷＳ３が表示する所定回数及び総加算回数の指定画面の一例を示す図である。計測ＷＳ３は、計測装置２による磁場データの計測を開始する際に、図９の画面をディスプレイに表示させる。
【0104】
図９に示す例では、２０００回、２５００回、３０００回及び３５００回のそれぞれが所定回数に該当する。また４０００回は総加算回数に該当する。
【0105】
生体データ計測システム１では、加算回数が所定回数に達するたびに加算平均処理を行って、加算平均データと加算平均データにおける加算回数を対応付けてデータ保存サーバ５に保存する。また加算回数が所定回数に達するたびに、加算平均処理した加算平均データに基づいて、生体の複数の刺激部位ごとでの活動電流の強さを推定し、推定結果を表示して、推定結果が妥当であるか否か等を途中確認させることもできる。
【0106】
【0107】
ここで、図１０は、計測装置２による計測中に表示される計測画面と操作画面の一例を説明する図である。図１１は、加算平均データリストの一例を説明する図である。図１２は、活動電流推定中の表示画面の一例を説明する図である。図１３は、適切でない活動電流推定結果の表示画面の一例を説明する図である。図１４は、適切な活動電流の推定結果の表示画面の一例を説明する図である。
【0108】
図１０に示すように、計測画面６０は、操作画面６１と、計測データ画面６２とを含む。これらのうち、操作画面６１は、計測の開始や計測の終了、計測データの表示方法の変更等をユーザが指示するために操作する画面である。
【0109】
また、計測データ画面６２は、計測装置２が計測した磁場計測データを表示する画面である。計測データ画面６２は、ｘ計測データ画面６２１と、ｙ計測データ画面６２２と、ｚ計測データ画面６２３とを含む。
【0110】
ｘ計測データ画面６２１は、図２のｘ軸方向における磁場データを表示する。また、ｙ計測データ画面６２２は、図２のｙ軸方向における磁場データを表示し、ｚ計測データ画面６２３は、図２のｚ軸方向における磁場データを表示する。
【0111】
計測データ画面６２に表示された波形データ６３は、磁気センサアレイ２００に含まれる１つの磁気センサによる磁場データを表示するものである。波形データ６３の横軸は時間を示し、縦軸は磁場強度を示している。波形データ６３は、磁気センサアレイ２００に含まれる複数の磁気センサのそれぞれによる磁場データをリアルタイムで表示する。
【0112】
ｘ計測データ画面６２１、ｙ計測データ画面６２２及びｚ計測データ画面６２３のそれぞれに含まれる波形データ６３の数は、磁気センサアレイ２００に含まれる磁気センサの数に対応している。
【0113】
ここで、加算平均処理部３２は、時系列に取得される波形データ６３を加算平均処理する。具体的には、加算平均処理部３２は、波形データ６３における時間ごとでの磁場データを加算し、時間ごとでの磁場データの加算結果を加算回数で除算することで加算平均データを取得する。
【0114】
波形データ６３内に含まれるデータのうち、ノイズに該当するデータは、時間的にランダムに発生するため、加算平均処理することで打ち消される。これに対し、波形データ６３内に含まれるデータのうち、磁場データは、加算すると積み上げられる。これにより加算平均処理は、ノイズに該当するデータと比較して磁場データを増幅させることができる。
【0115】
波形データ６３は、生体の磁場データの一例であり、また生体データの一例である。複数の波形データ６３が加算平均されることで生成された波形データが加算平均データに対応する。
【0116】
また、図１０において、破線の四角で示した開始ボタン６４は、ユーザが、推定部４２による活動電流の強さの推定処理の開始指示を行うためのボタンである。ユーザが、図３におけるポインティングデバイス５１２のカーソルを用いて開始ボタン６４を押すと、解析ＷＳ４における指示受付部４４は、データ保存サーバ５に加算平均データのリストを要求して取得する。
【0117】
表示部４５は、取得された加算平均データのリストをディスプレイ５０６に表示させる。図１１は、加算平均データリストの一例を示している。図１１における計測Ａ及び計測Ｂはそれぞれ生体としての被検者を示す情報を示している。トリガーα及びβは、生体の複数の刺激部位のそれぞれに与える電気刺激を発生させるためのトリガー信号を示し、複数のトリガー信号の一例である。最終フラグは、総加算回数を示す情報である。
【0118】
表示された加算平均データリストの中からユーザが加算回数を選択すると、推定部４２は通信部４１及び通信部５１を介して、保存部５２を参照して加算回数に対応する加算平均データを取得する。推定部４２は、ユーザが選択した加算回数に対応する加算平均データを用いて推定処理を行うことができる。
【0119】
ここで、開始ボタン６４が押された後、図１２に示すような脊髄位置指定画面７０を表示させてもよい。脊髄位置指定画面７０は、被検者毎に異なる脊髄位置を指定するために用いられる画面である。図１２に示すように、脊髄位置指定画面７０は、レントゲン撮影画像画面７０１と、推定開始指示受付画面７０２とを含む。
【0120】
これらのうち、レントゲン撮影画像画面７０１は、被検者１００（図２参照）の側方からのレントゲン撮影画像画面７０１を表示する。なお、レントゲン撮影画像は、通信部４１を介して外部装置から入力された画像である。
【0121】
ユーザは、レントゲン撮影画像画面７０１を視認し、図３のポインティングデバイス５１２のカーソルを用いて画面上の点を指定することで、生体の内部で活動電流の強さを推定する位置を指定することができる。指定された位置を含む領域内で活動電流の強さが推定される。
【0122】
図１２において、レントゲン撮影画像画面７０１に含まれる曲線７０１１は、ユーザがレントゲン撮影画像画面７０１上で指定した点が含まれるように自動的に描画された曲線であり、側方から見た被検者１００の内部における脊髄位置に対応している。
【0123】
ユーザは、曲線７０１１上の点を指定後、推定開始指示受付画面７０２における開始ボタン７０２１を図３のポインティングデバイス５１２のカーソルを用いて押すことで、指定した点に対応する被検者１００内部の位置での活動電流の強さを推定する処理を開始させることができる。
【0124】
次に図１３に示すように、推定結果表示画面８０は、レントゲン撮影画像画面７０１と、活動電流の強さの分布図８０１を含んでいる。分布図８０１は、磁気センサアレイ２００に含まれる各磁気センサによって計測されたデータに基づき、推定された活動電流の強さの二次元分布を、色に置き換えて表示した図である。
【0125】
【0126】
図１３の場合、分布図８０１内では、本来生体内部の活動電流が存在しないはずの場所、すなわち人体の外に強い活動電流が推定されていることが分かる。これは、活動電流の強さの推定結果にノイズの影響が大きく生じていることを示している。換言すると、図１３は、加算平均処理部３２による処理に用いられた計測データ数が十分でなかったために、ノイズが多く含まれた磁場の計測データに基づいて活動電流が推定された、適切でない場合の推定結果を示している。
【0127】
一方、図１４の場合、分布図９０１内では、本来生体内部の活動電流が存在するはずの場所、すなわち人体の脊髄上に強い活動電流が推定されていることが分かる。これは、活動電流の強さの推定結果で、ノイズの影響が抑制されていることを示している。
【0128】
換言すると、図１４は、加算平均処理部３２による処理に用いられた計測データ数が十分であったために、ノイズが抑制された磁場の計測データに基づいて活動電流が推定された、適切である場合の推定結果を示している。
【0129】
ユーザは、図１３や図１４で例示される分布図を視認して、推定結果が適切であるか否かを判定できる。適切でないと判定された場合は、ユーザが開始ボタン６４（図１０参照）を操作することで、再度、推定処理の開始指示が行われる。
【0130】
また、開始指示を受け付ける代わりに、図１３に一点鎖線の四角で示したエディットボックス８１に計測データ数を数値で入力し、入力した数値に該当する数の計測データ数が取得された時に、推定部４２による推定処理が開始されるようにしてもよい。
【0131】
＜生体データ処理装置１０の作用効果＞
以上説明したように、実施形態では、加算平均処理部３２は、少なくとも１つの部位に対応する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する。保存部５２は、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、該加算平均データにおける加算回数とを対応付けて保存する。推定部４２は、指定された加算回数に基づき保存部５２を参照して取得される少なくとも１つの刺激部位ごとでの加算平均データを用いて、生体データに基づく処理を行う。例えば、推定部４２は、少なくとも１つの刺激部位ごとでの活動電流の強さを推定することにより生体データに基づく処理を行う。
以上説明したように、実施形態では、加算平均処理部３２は、少なくとも１つの部位に対応する刺激に紐づいたトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理する。保存部５２は、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、該加算平均データにおける加算回数とを対応付けて保存する。推定部４２は、指定された加算回数に基づき保存部５２を参照して取得される少なくとも１つの刺激部位ごとでの加算平均データを用いて、生体データに基づく処理を行う。例えば、推定部４２は、少なくとも１つの刺激部位ごとでの活動電流の強さを推定することにより生体データに基づく処理を行う。
【0132】
このように、少なくとも１つの刺激部位ごとで同じ加算回数による加算平均データを用いることで、刺激部位ごとでの加算回数の違いをなくし、刺激部位ごとで加算平均処理の効果を等しく得ることができる。
【0133】
また本実施形態では、推定部４２による活動電流の強さの推定結果に基づき、計測装置２による計測の中断又は延長を制御する計測制御部３３を有する。また推定部４２による活動電流の強さの推定結果を表示する表示部４５と、計測装置２による計測の中断又は延長の何れか一方の指示を受け付ける指示受付部４４とを有し、計測制御部３３は、指示受付部４４が受け付けた指示に応じて計測装置２による計測を中段又は延長させる。
【0134】
これにより、計測の途中に推定結果が妥当であるか否か、又は計測を延長する必要があるか否かの判断をユーザに判断させ、判断結果に応じて計測装置２による計測を中断又は延長させることができる。
【0135】
換言すると、計測の途中に推定結果が妥当であるか否か、又は計測を延長する必要があるか否かの判断をユーザに途中確認させ、確認結果に応じて計測装置２による計測を中断又は延長させることができる。
【0136】
これにより、無駄な計測を抑制するとともに、必要に応じて計測を延長することができる。
【0137】
［第２実施形態］
次に第２実施形態に係る生体データ計測システム１ａについて説明する。なお、第１実施形態で説明したものと同じ構成部には、同じ符号を付し、重複する説明を適宜省略する。この点は、以降に示す実施形態においても同様とする。
次に第２実施形態に係る生体データ計測システム１ａについて説明する。なお、第１実施形態で説明したものと同じ構成部には、同じ符号を付し、重複する説明を適宜省略する。この点は、以降に示す実施形態においても同様とする。
【0138】
本実施形態では、加算回数に応じて、推定部に推定方法を切り替えさせることで、推定精度を高く確保する。具体的には、推定部は、加算回数が総加算回数に対応する加算回数以外である場合には、総加算回数に対応する加算回数の場合と比較して高速な処理を行い、加算回数が総加算回数に対応する加算回数である場合には、加算回数が総加算回数に対応する加算回数以外である場合と比較して高精度な処理を行う。
【0139】
【0140】
切替部４７は、第１指定受付部４３が受け付けた加算回数の指定回数に応じて、推定部４２による生体の活動電流の推定処理方法を切り替えさせる。
【0141】
推定方法には、例えば上述した「Array-Gain Constraint Minimum-Norm Spatial Filter With Recursively Updated Gram Matrix」を適用でき、推定処理方法の切替では、この推定処理の繰り返し回数の切替が行われる。推定処理方法は、生体データ処理部による処理方法の一例である。指定回数と繰り返し回数とを予め対応付けておくことで、切替部４７は、指定回数に応じて繰り返し回数を切り替えることができる。
【0142】
このように、本実施形態では、加算回数に応じて推定部４２に推定処理方法を切り替えさせることで、推定精度を高く確保することができる。また高精度化を行わない場合には、推定部に推定方法を切り替えさせることで、推定を高速に行うことができる。
【0143】
［第３実施形態］
次に第３実施形態に係る生体データ計測システム１ｂについて説明する。
次に第３実施形態に係る生体データ計測システム１ｂについて説明する。
【0144】
図１６は、生体データ計測システム１ｂが有する解析ＷＳ４ｂの機能構成の一例を説明するブロック図である。図１６に示すように、解析ＷＳ４ｂは、所定回数の指定を受け付ける第２指定受付部４８と、加算回数を自動更新するために使用される更新間隔回数の指定を受け付ける第３指定受付部４９と、を有する。
【0145】
解析ＷＳ４ｂは、第２指定受付部４８により、所定回数として次回の加算回数を受け付ける。解析ＷＳ４ｂは、加算回数が所定回数に達するたびに、第２指定受付部４８により受け付けた次回の加算回数により所定回数を更新できる。また解析ＷＳ４ｂは、加算回数が所定回数に達するたびに、第３指定受付部４９により受け付けられた更新間隔回数が所定回数に加算された回数に、所定回数を自動更新する。
【0146】
図１７及び図１８は、本実施形態に係る所定回数及び総加算回数の指定画面１７０を例示する図であり、図１７は第１例の図、図１８は第２例の図である。指定画面１７０は、ディスプレイ５０６等に表示され、キーボード５１１又はポインティングデバイス５１２等を用いて操作入力されるＧＵＩ（グラフィックユーザインターフェース）である。
【0147】
【0148】
スイッチ１７１は、加算回数の自動更新を行う（ＯＮ）か、又は行わないか（ＯＦＦ）の設定を切り替えるスイッチである。
【0149】
第１入力ボックス１７２は、更新間隔回数を意味する「加算回数の自動更新間隔」を入力するためのボックスである。図１６における第３指定受付部４９は、第１入力ボックス１７２を介して更新間隔回数を受け付ける。
【0150】
第２入力ボックス１７３は、次回の加算回数を入力するためのボックスである。図１６における第２指定受付部４８は、第２入力ボックス１７３を介して、所定回数として次回の加算回数を受け付ける。
【0151】
第３入力ボックス１７４は、総加算回数を入力するためのボックスである。
【0152】
図１７では、スイッチ１７１は、加算回数の自動更新を行わない（ＯＦＦ）ように設定されている。自動更新は行われないため、第１入力ボックス１７２は無効化されている。なお、図１７の第１入力ボックス１７２におけるドットハッチングは、無効化されていることを意味する。
【0153】
第２入力ボックス１７３には、次回の加算回数として「１０００回」が入力されている。解析ＷＳ４ｂは、加算回数が１０００回に達した後に次回の加算回数として例えば１５００回の入力を受け付ける等して、所定回数を任意の回数に指定できる。
【0154】
図１８では、スイッチ１７１は、加算回数の自動更新を行う（ＯＮ）ように設定されている。自動更新が行われるため、第１入力ボックス１７２には、加算回数の自動更新間隔（更新間隔回数）として「１０００回」が入力されており、第２入力ボックス１７３は、無効化されている。なお、図１８の第２入力ボックス１７３におけるドットハッチングは無効化されていることを意味する。
【0155】
解析ＷＳ４ｂは、加算回数が所定回数に達するたびに、所定回数に更新間隔回数の「１０００回」を加算して所定回数を更新する。図１１では、総加算回数として４０００回が指定されているため、所定回数は、加算回数が所定回数に達するたびに、自動更新により１０００回、２０００回、３０００回、４０００回と変化する。
【0156】
このようにして、生体データ計測システム１ｂが第２指定受付部４８と、第３指定受付部４９と、を有することにより、生体データ計測システム１ｂのユーザは、生体データ計測システム１ｂによる計測の状況に応じて、所定回数の指定方法を適宜選択可能になる。その結果、生体データ計測システム１ｂは利便性を向上させることができる。
【0157】
［第４実施形態］
次に、第４実施形態に係る生体データ計測システム１８１について説明する。
次に、第４実施形態に係る生体データ計測システム１８１について説明する。
【0158】
生体データ計測システムでは、計測中に一時的にでも生体信号が低下すること、またはノイズが混入することを検出することが好ましい。
【0159】
例えば、特開２０１９―１６２２５３号公報等に記載の方法では、突発的な大きなノイズを検出して除外することは可能であるが、加算平均処理を行わないければ明らかにならないような小さいノイズを検出して除外することはできない。また生体信号の低下を検出して除外することはできない。例えば計測中に、何らかの異常が発生した場合に、この異常を精度よく検出することは難しく、最終的な加算平均データの品質を保証するためには十分ではない。
【0160】
本実施形態では、生体データ計測システム１８１による計測中に異常を検出し、異常が発生した区間のデータは排除することにより、十分な信号品質を保証しつつ、計測の信頼性を保証する。
【0161】
＜生体データ計測システム１８１の全体構成＞
図１９は、生体データ計測システム１８１の全体構成の一例を説明する図である。生体データ計測システム１８１は、計測装置１８２と、計測ＷＳ１８３と、解析ＷＳ１８４と、データ保存サーバ１８５とを有する。これらは有線又は無線によって、相互に通信可能に接続されている。これらのうち、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５は、生体データ処理装置１８６を構成する。
図１９は、生体データ計測システム１８１の全体構成の一例を説明する図である。生体データ計測システム１８１は、計測装置１８２と、計測ＷＳ１８３と、解析ＷＳ１８４と、データ保存サーバ１８５とを有する。これらは有線又は無線によって、相互に通信可能に接続されている。これらのうち、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５は、生体データ処理装置１８６を構成する。
【0162】
計測装置１８２は、複数のトリガー信号のそれぞれに応じた電気刺激等の刺激に対して生体の複数の部位で発生する磁場データを計測する脊磁計である。磁場データは計測データの一例である。計測装置１８２は、複数の部位ごとでの計測結果である磁場データを、複数の部位のそれぞれに対応した複数のトリガー信号とともに計測ＷＳ１８３に送信する。
【0163】
計測ＷＳ１８３は、計測装置１８２から受信したトリガー信号をカウントして加算回数を取得し、加算回数が所定回数に達するたびに磁場データの加算平均処理を行う。加算平均処理は、所定回数から別途指定していた回数分遡った磁場データのみを対象とする。別途指定していた回数をセグメント幅と呼称し、得られたデータをセグメント加算データと呼称する。
【0164】
計測装置１８２は、セグメント加算データを、セグメント幅と、セグメント幅の始点の加算平均回数と、終点の加算平均回数を対応付けてデータ保存サーバ５に送信する。以降では、上記のセグメント加算に係る3種の情報のうち、少なくとも一つの情報を含む情報をセグメント加算に係る情報と記述する。
【0165】
データ保存サーバ１８５は、計測ＷＳ１８３から受信した加算平均データと、セグメント加算に係る情報と、を対応付けて保存する。
【0166】
解析ＷＳ１８４は、データ保存サーバ１８５を参照してすべてのセグメント加算データを取得し、その中から不正なセグメント加算データを除外したうえですべてのセグメント加算データを加算平均し、加算平均データを得る。解析ＷＳ１８４は、得られた加算平均データに基づいて、活動電流の強さを推定する。解析ＷＳ１８４は、解析ＷＳ１８４が有するディスプレイに推定結果を表示したり、データ保存サーバ１８５に送信して保存させたり、外部サーバ等の外部装置に送信したりすることができる。
【0167】
なお、本実施形態では、生体データ処理装置１８６が、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５の３台の装置により構成される例を示すが、これに限定されるものではない。生体データ処理装置１８６は、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５の機能を集約した１台の装置を有しもよいし、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５の機能を分散させた４台以上の装置を有してもよい。
【0168】
生体データ計測システム１８１は、計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５以外の装置を通信可能に有してもよいし、計測装置１８２以外の他の生体データの計測装置を通信可能に有してもよい。
【0169】
計測装置１８２の構成は、第１実施形態において図２を用いて説明した計測装置２の構成と同様である。コンピュータのハードウェア構成は、第１実施形態において図３を用いて説明したコンピュータのハードウェア構成と同様である。
【0170】
＜生体データ処理装置１８６の機能構成例＞
図２０、図２１及び図２２を参照して、生体データ処理装置１８６を構成する計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５それぞれの機能構成について説明する。
図２０、図２１及び図２２を参照して、生体データ処理装置１８６を構成する計測ＷＳ１８３、解析ＷＳ１８４及びデータ保存サーバ１８５それぞれの機能構成について説明する。
【0171】
（計測ＷＳ１８３の機能構成例）
図２０は、計測ＷＳ１８３の機能構成の一例を説明するブロック図である。計測ＷＳ１８３は、通信部１８７と、加算平均処理部１８８と、計測制御部１８９と、表示部１９０とを有する。
図２０は、計測ＷＳ１８３の機能構成の一例を説明するブロック図である。計測ＷＳ１８３は、通信部１８７と、加算平均処理部１８８と、計測制御部１８９と、表示部１９０とを有する。
【0172】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図２０は、計測ＷＳ１８３が有する主な構成を示すが、計測ＷＳ１８３はこれら以外の構成を有してもよい。
【0173】
通信部１８７は、計測装置１８２、解析ＷＳ１８４及びデータ保存サーバ１８５のそれぞれとの間でデータ及び信号の送受信を行う。
【0174】
加算平均処理部１８８は、ユーザがキーボード５１１（図３参照）等を用いて入力した所定回数、セグメント幅、及び総加算回数の情報を取得する。総加算回数は、加算平均処理部１８８が磁場データの加算を行う総回数である。なお、加算平均処理部１８８は、ＨＤ５０４等に予め保存された所定回数、セグメント幅及び総加算回数の情報を、ＨＤ５０４等を参照して取得してもよい。
【0175】
加算平均処理部１８８は、計測装置１８２から通信部１８７を介して、少なくとも１つのトリガー信号を受信する。また加算平均処理部１８８は、少なくとも１つのトリガー信号に対応し、所定回数からセグメント幅分だけ遡って計測装置１８２が計測した磁場データを、通信部１８７を介して受信して加算平均処理する。
【0176】
加算平均処理は、計測装置１８２が計測した磁場データを所定回数からセグメント幅分だけ遡って加算処理した値を加算回数で除算して平均値を算出する処理である。加算平均処理部１８８は、処理結果であるセグメント加算データを表示部に送る。さらに加算平均処理部１８８は、セグメント加算データをセグメント加算に係る情報に対応付け、通信部１８７を介してデータ保存サーバ１８５に送信する。
【0177】
表示部１９０は、セグメント加算データを波形データとしてディスプレイに表示させる。表示部１９０は、セグメント加算データを表示させる際には、複数の所定回数に紐づいたセグメント加算データを波形データとして並列表示させてもよい。また表示部１９０は、セグメント加算データを何らかの手段により解析した結果である周波数スペクトル等を表示してもよい。計測ＷＳ１８３は、表示された波形データを評価し、セグメントにおいて計測が正常に行われた否かを判断し、記録する。
【0178】
計測ＷＳ１８３は、所定のアルゴリズムに従って、セグメントにおいて計測が正常に行われた否かを判断する。但し、ユーザが目視によりセグメントにおいて計測が正常に行われた否かを判断してもよい。例えばユーザは、周波数スペクトルを目視し、生体信号以外の周波数の箇所にピークが立っている場合には異常と判断したり、波形データから得られる生体信号の振幅が他のセグメント加算データの振幅と比較して明らかに低下していれば異常と判断したりすることができる。「明らかに低下」とは、例えばセグメント加算データに含まれるノイズレベル以上に低下していることをいう。
【0179】
計測ＷＳ１８３は、セグメント加算データが異常なデータであると判断された場合には、加算平均データを生成する際にセグメント加算データを除外するために、セグメント幅分だけ加算回数を増加させる。例えば、計測ＷＳ１８３は、キーボード５１１等を用いたユーザによる計測の延長指示を受け付け、セグメント幅分だけ加算回数を増加させる。計測制御部１８９は、延長指示に応答して計測装置１８２に計測を延長させることができる。
【0180】
計測制御部１８９は、表示部１９０に表示された波形データを評価した結果に基づく指示を、通信部１８７を介して受信し、この評価結果に基づき、計測装置１８２に計測を中断させ、又は延長させる。計測制御部１８９は、計測の中断又延長の指示を、指示があった任意のタイミングにおいて、割込みデータとして受信可能である。
【0181】
（データ保存サーバ１８５の機能構成例）
図２１は、データ保存サーバ１８５の機能構成の一例を説明するブロック図である。データ保存サーバ１８５は、通信部１９１と、保存部１９２とを有する。
図２１は、データ保存サーバ１８５の機能構成の一例を説明するブロック図である。データ保存サーバ１８５は、通信部１９１と、保存部１９２とを有する。
【0182】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図２１は、データ保存サーバ１８５が有する主な構成を示すが、データ保存サーバ１８５はこれら以外の構成を有してもよい。
【0183】
通信部１９１は、計測ＷＳ１８３及び解析ＷＳ１８４のそれぞれとの間でデータ及び信号の送受信を行う。
【0184】
保存部１９２は、通信部１９１を介して計測ＷＳ１８３から受信したセグメント加算データ５２４と、セグメント加算に係る情報５２５と、を対応付けて保存する。なお、セグメント加算データ５２４は、複数のセグメント加算データの総称表記であり、セグメント加算に係る情報５２５は、複数のセグメント加算に係る情報の総称表記である。
【0185】
（解析ＷＳ１８４の機能構成例）
次に図２２は、解析ＷＳ１８４の機能構成の一例を説明するブロック図である。解析ＷＳ１８４は、通信部１９３と、加算平均処理部１９４と、推定部１９５と、第４指定受付部１９６と、表示部１９７とを有する。
次に図２２は、解析ＷＳ１８４の機能構成の一例を説明するブロック図である。解析ＷＳ１８４は、通信部１９３と、加算平均処理部１９４と、推定部１９５と、第４指定受付部１９６と、表示部１９７とを有する。
【0186】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図２２は、解析ＷＳ１８４が有する主な構成を示すが、解析ＷＳ１８４はこれら以外の構成を有してもよい。
【0187】
通信部１９３は、計測ＷＳ１８３及びデータ保存サーバ１８５のそれぞれとの間でデータ及び信号の送受信を行う。
【0188】
加算平均処理部１９４は、通信部１９３を介して保存部１９２が保存するセグメント加算データのリストを取得し、取得したセグメント加算データのリストを表示部１９７によりディスプレイ５０６等に表示させる。第４指定受付部１９６は、ユーザがセグメント加算データのリストを視認して選択した結果から加算回数の情報を受け付けることができる。選択された少なくとも一つのセグメント加算データは、それぞれのデータの加算回数に基づいて加算平均処理が行われる。解析ＷＳ１８４は、提示されたセグメント加算データのリストから異常なセグメント加算データを除外し、加算平均データを取得できる。
【0189】
推定部１９５は、生体データに基づく処理を行う生体データ処理部の一例である。推定部１９５は、加算平均処理部１９４が生成した磁場データに基づき、生体の活動電流の強さを推定する処理を行う。推定部１９５は、取得した加算平均データを用いて、複数の部位ごとでの活動電流の強さを推定する。この推定アルゴリズムには、上述した「Array-Gain Constraint Minimum-Norm Spatial Filter With Recursively Updated Gram Matrix」（例えば、「https://ieeexplore.ieee.org/document/5415622」参照）等を用いることができる。
【0190】
表示部１９７は、推定部１９５による活動電流の強さの推定結果を表示する。例えば表示部１９７は、ディスプレイ５０６に推定結果を表示し、ユーザに視認させることができる。また表示部１９７は、取得した加算平均データやセグメント加算データを示す波形データを表示することもできる。
【0191】
ユーザは、ディスプレイ５０６に表示された活動電流の強さの推定結果を視認し、セグメントにおける計測が正常に行われているかを確認できる。解析ＷＳ１８４は、セグメント加算データが異常なデータであるとユーザにより判断された場合には、加算平均データを生成する際にセグメント加算データを除外する。
【0192】
＜生体データ処理装置１８６の動作例＞
次に図２３を参照して、生体データ処理装置１８６の動作について説明する。図２３は、生体データ処理装置１８６の動作の一例を示すフローチャートである。生体データ処理装置１８６は、ユーザによる生体データ処理の開始指示を受け付けること等により図２３の動作を開始する。
次に図２３を参照して、生体データ処理装置１８６の動作について説明する。図２３は、生体データ処理装置１８６の動作の一例を示すフローチャートである。生体データ処理装置１８６は、ユーザによる生体データ処理の開始指示を受け付けること等により図２３の動作を開始する。
【0193】
まず、ステップＳ１０１において、生体データ処理装置１８６は、加算平均処理部１８８により、ユーザがキーボード５１１等を用いて入力した所定回数とセグメント幅及び総加算回数の情報を取得する。なお、加算平均処理部１８８は、ＨＤ５０４等に予め保存された所定回数とセグメント幅及び総加算回数の情報をＨＤ５０４等から取得してもよい。
【0194】
続いて、ステップＳ１０２において、生体データ処理装置１８６は、加算平均処理部１８８により、少なくとも一つのトリガー信号と、このトリガー信号に対応して計測された磁場データと、を計測装置１８２から通信部１８７を介して受信して蓄積する。なお、本実施形態の説明では、一つのトリガー信号に対応した一つの磁場データをエポックデータと呼称する。
【0195】
続いて、ステップＳ１０３において、生体データ処理装置１８６は、加算平均処理部１８８により、受信したトリガー信号をカウントしてトリガーごとに加算回数を取得し、加算回数が所定回数に達したか否かを判定する。
【0196】
ステップＳ１０３で達したと判定された場合には（ステップＳ１０３、Ｙｅｓ）、ステップＳ１０４において、生体データ処理装置１８６は、加算平均処理部１８８により、所定回数からセグメント幅分だけ遡って磁場データの加算平均処理を実行する。一方、達していないと判定された場合には（ステップＳ１０３、Ｎｏ）、生体データ処理装置１８６は、ステップＳ１０２の動作を再度行う。
【0197】
続いて、ステップＳ１０５において、生体データ処理装置１８６は、加算平均処理部１８８により、セグメント加算平均処理結果であるセグメント加算データと、セグメント加算データにおける加算回数の情報と、を対応付け、通信部１８７を介してデータ保存サーバ１８５に送信する。データ保存サーバ１８５は、受信したセグメント加算データと、セグメント加算に係る情報と、を対応付けて保存する。また生体データ処理装置１８６は、セグメント加算データを計測ＷＳ１８３の表示部１９０に表示させる。
【0198】
続いて、ステップＳ１０６において、生体データ処理装置１８６は、計測ＷＳ１８３により、表示部１９０に表示されたセグメント加算データを確認し、正常であるか否かを判断する。例えば計測ＷＳ１８３は、キーボード５１１等を用いたユーザの操作入力を受け付けることにより、正常であるか否かを判断できる。
【0199】
ステップＳ１０６において、正常ではないと判断された場合（ステップＳ１０６、Ｎｏ）、ステップＳ１０７において、生体データ処理装置１８６は、計測ＷＳ１８３により、このセグメントを異常セグメントとして記録する。
【0200】
続いて、ステップＳ１０８において、生体データ処理装置１８６は、計測ＷＳ１８３により、最終加算回数の増加、すなわち加算平均回数を追加して即時に最終加算回数を更新し、または計測を中断する。その後、生体データ処理装置１８６は、ステップＳ１０２に動作を移行する。
【0201】
なお、生体データ処理装置１８６は、ステップＳ１０１乃至Ｓ１０７における任意の順番においてステップＳ１０８の動作を行ってもよい。また生体データ処理装置１８６は、計測ＷＳ１８３により、セグメント、所定回数についても同様に任意のタイミングにおいて指示を受けて設定を更新してもよい。
【0202】
続いて、ステップＳ１０９において、生体データ処理装置１８６は、計測制御部１８９により、加算回数が最終加算回数（総加算回数）に到達したか否かを判定する。なお、この判定は、計測制御部１８９でなく、加算平均処理部１８８が行ってもよい。
【0203】
ステップＳ１０９において、達していないと判定された場合には（ステップＳ１０９、Ｎｏ）、生体データ処理装置１８６は、ステップＳ１０２以降の動作を再度行う。一方、ステップＳ１０９において、到達したと判定された場合には（ステップＳ１０９、Ｙｅｓ）、生体データ処理装置１８６は、計測装置１８２に計測を終了させる。
【0204】
このようにして、生体データ処理装置１８６は、計測ＷＳ１８３により、加算平均処理を実行し、また計測の中断又は延長の指示に応じた計測装置２の制御を行うことができる。
【0205】
続いて、ステップＳ１１０において、生体データ処理装置１８６は、解析ＷＳ１８４の第４指定受付部１９６により、通信部１９３を介して保存部１９２が保存するセグメント加算データのリストを取得し、取得したセグメント加算データのリストを表示部１９７によりディスプレイ５０６等に表示させる。ユーザは、表示された加算平均データのリストを視認し、キーボード５１１等を用いて加算平均データを選択できる。このとき、生体データ処理装置１８６は、異常セグメントとして記録されているセグメント加算データを除外する。
【0206】
続いて、ステップＳ１１１において、生体データ処理装置１８６は、選択された少なくとも一つのセグメント加算データを、セグメント幅に応じて加算平均処理する。この処理により、生体データ処理装置１８６は、少なくとも一つのトリガーに対応した加算平均データを得る。
【0207】
このようにして、生体データ処理装置１８６は、解析ＷＳ４により、生体の活動電流の推定処理に必要な加算平均データを得ることができる。
【0208】
＜各種表示画面例＞
次に、生体データ計測システム１８１において表示される各種表示画面について説明する。
次に、生体データ計測システム１８１において表示される各種表示画面について説明する。
【0209】
【0210】
図２４及び図２５は、本実施形態に係る所定回数、セグメント幅及び総加算回数の指定画面を例示する図である。図２４では、スイッチ２１１はセグメント加算を行わないようにオフに設定されている。この場合、セグメント実行間隔とセグメント幅というセグメント加算平均処理に関わるパラメータは編集できないように入力ボックス２１２及び２１３は無効化されている。なお、図２４の入力ボックス２１２及び２１３におけるドットハッチングは入力を受け付けない状態にされていることを意味する。以降ドットハッチングされているボックスは同様に扱う。入力ボックス２１４は、総加算回数が入力されるものである。
【0211】
図２５は、本実施形態に係る所定回数を更新する間隔、セグメント幅及び総加算回数の指定画面を例示する図である。図２５では、スイッチ２１５は、加算回数の自動更新を行うようにオンに設定されている。自動更新が行われるため、入力ボックス２１６には、加算回数の自動更新間隔（更新間隔回数）として「５００回」が入力されている。入力ボックス２１７は、自動更新間隔と同じ値が表示され、この値がセグメント幅として指定される。入力ボックス２１８は、総加算回数が入力されるものである。
【0212】
計測ＷＳ１８３は、加算回数が所定回数に達するたびに、所定回数に加算実行間隔回数の「５００回」を加算して所定回数を更新する。図２５では、総加算回数として４０００回が指定されているため、所定回数は、加算回数が所定回数に達するたびに、自動更新により１０００回、１５００回、・・・、３５００回、４０００回と変化する。
【0213】
生体データ計測システム１８１では、加算回数が所定回数に達するたびにセグメント加算平均処理を行って、セグメント加算データとセグメント加算に係る情報を対応付けてデータ保存サーバ１８５に保存する。
【0214】
【0215】
【0216】
図２６及び図２７は、セグメント加算データの表示結果の一例を説明する図である。図２６は、セグメント加算データを波形データとして評価する図である。図２７はセグメント加算データに対して周波数解析を適用し、周波数スペクトルを用いて評価する図である。
【0217】
図２６は、２つのセグメント加算データを表示させた際の波形データ６３のうち、特に生体信号を強く検出したセンサの波形データを抽出し、縦に並べて重畳表示した例である。実線グラフは加算回数１～５００回のセグメント加算データ、破線グラフは加算回数５０１～１０００回のセグメント加算データをそれぞれ示している。破線グラフの５０１～１０００回のセグメント加算データにおいては、図中矢印で示した箇所で生体信号の強度低下が起きており、何らかの異常が発生していると見なされる。
【0218】
図２７は、１～５００回、５０１～１０００回、１００１～１５００回のセグメント加算データに対して周波数解析が行われた結果を周波数スペクトルとして表示している。１～５００回ではノイズは大きく存在していないが、５０１～１０００回の周波数スペクトル、１００１～１５００回の周波数スペクトルにおいては、特定の周波数にピークが存在しており、何らかのノイズがそれぞれのセグメントに混入していることが分かる。これにより、ユーザは、５０１～１０００回のセグメント、１００１～１５００回のセグメントに異常が発生していると判断できる。解析ＷＳ１８４は、このセグメント加算データに対応した周波数スペクトルも重畳させて表示してもよい。
【0219】
表示部１９７は、取得された加算平均データのリストをディスプレイ５０６に表示させる。図２８は、加算平均データリストの一例を示している。図２８における計測Ａ及び計測Ｂはそれぞれ生体としての被検者を示す情報である。トリガーα及びβは、生体の少なくとも一つの部位にそれぞれ与える刺激に対応するトリガー信号を示し、少なくとも一つのトリガー信号の一例である。異常フラグは、対応したセグメント加算データが異常であるとユーザが指定したことを示している。
【0220】
表示されたセグメント加算データリストの中からユーザが少なくとも一つのセグメント加算データを選択すると、加算平均処理部１９４は、選択された各セグメントに対応した加算回数に基づいて、選択されたセグメント加算データより加算平均処理を行い、加算平均データを生成し、データ保存サーバ１８５に保存する。推定部１９５は、手続きにより生成された加算平均データを用いて推定処理を行うことができる。
【0221】
＜生体データ処理装置１８６の作用効果＞
以上説明したように、生体データ処理装置１８６は、加算平均処理部１８８により、複数のトリガー信号に応じて生体における複数の部位で計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理し、表示させる。生体データ処理装置１８６は、処理結果であるセグメント加算データと、このセグメント加算データにおける加算回数とを対応付けて保存部１９２に保存させる。
以上説明したように、生体データ処理装置１８６は、加算平均処理部１８８により、複数のトリガー信号に応じて生体における複数の部位で計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理し、表示させる。生体データ処理装置１８６は、処理結果であるセグメント加算データと、このセグメント加算データにおける加算回数とを対応付けて保存部１９２に保存させる。
【0222】
その後、生体データ処理装置１８６は、加算平均処理部１８８により、表示されたセグメント加算データに基づいてユーザが異常と判断したセグメントを除外したうえで加算平均を行い、加算平均データを得る。
【0223】
このように、生体データ処理装置１８６は、セグメント加算データを表示、評価及び加算することで、異常なデータを取り除いて加算平均データを得ることができる。
【0224】
また生体データ処理装置１８６は、セグメント加算データの表示結果に基づき、計測装置１８２による計測の中断又は延長を制御する計測制御部１８９を有する。これにより、生体データ処理装置１８６は、計測中の一期間であるセグメントにおいて異常が発生した場合に、このセグメントの計測データを補うように計測を延長させるか、計測そのものに異常が発生したとして計測を中断させることができる。
【0225】
換言すると、生体データ処理装置１８６は、計測の途中にある区間の計測が正常に行われているか否か、計測を延長もしくは中断する必要があるか否かの判断をユーザに途中確認させ、確認結果に応じて計測装置１８２による計測を延長又は中断させることができる。これにより、生体データ処理装置１８６は、異常が発生している区間のデータを取り除いたうえで、当初予定していた加算回数を保証することができる。
【0226】
［第５実施形態］
第５実施形態に係る生体データ計測システム１８１について説明する。本実施形態では、セグメント加算平均データを実行する所定回数とセグメント幅を独立して指定することにより、異常の検出速度を向上させる。生体データ計測システム１８１の全体構成は、第４実施形態と同様であるため、ここでは説明を割愛する。
第５実施形態に係る生体データ計測システム１８１について説明する。本実施形態では、セグメント加算平均データを実行する所定回数とセグメント幅を独立して指定することにより、異常の検出速度を向上させる。生体データ計測システム１８１の全体構成は、第４実施形態と同様であるため、ここでは説明を割愛する。
【0227】
第５実施形態に係る生体データ計測システム１８１は、計測ＷＳ１８３ａと、解析ＷＳ１８４ａと、データ保存サーバ１８５ａと、を有する。計測ＷＳ１８３ａ、解析ＷＳ１８４ａ及びデータ保存サーバ１８５ａは、生体データ処理装置１８６ａを構成する。
【0228】
解析ＷＳ１８４ａは、データ保存サーバ１８５ａを参照してすべてのセグメント加算データを取得し、表示部でセグメント加算データを表現する波形データ、もしくはセグメント加算データに基づいて推定部により推定された体内の電流分布を表示する。解析ＷＳ１８４ａは、表示されたデータに基づいてその中から不正なセグメント加算データを検出し、検出されたセグメント加算データを除外したうえですべてのセグメント加算データを加算平均し、加算平均データを得る。解析ＷＳ１８４ａは、得られた加算平均データに基づいて、活動電流の強さを推定する。解析ＷＳ１８４ａは、解析ＷＳ１８４ａが有するディスプレイに推定結果を表示したり、データ保存サーバ１８５に送信して保存させたり、外部サーバ等の外部装置に送信したりすることができる。
【0229】
第５実施形態に係る生体データ計測システム１８１は、計測ＷＳ１８３ａ、解析ＷＳ１８４ａ及びデータ保存サーバ１８５ａ以外の装置を通信可能に有してもよいし、計測装置１８２以外の他の生体データの計測装置を通信可能に有してもよい。
【0230】
【0231】
（計測ＷＳ１８３ａの機能構成例）
図２９は、計測ＷＳ１８３ａの機能構成の一例を説明するブロック図である。計測ＷＳ１８３ａは、通信部２２１と、加算平均処理部２２２と、計測制御部２２３と、表示部２２４とを有する。
図２９は、計測ＷＳ１８３ａの機能構成の一例を説明するブロック図である。計測ＷＳ１８３ａは、通信部２２１と、加算平均処理部２２２と、計測制御部２２３と、表示部２２４とを有する。
【0232】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図２９は、計測ＷＳ１８３ａが有する主な構成を示すが、計測ＷＳ１８３ａはこれら以外の構成を有してもよい。
【0233】
通信部２２１は、計測装置１８２、解析ＷＳ１８４ａ及びデータ保存サーバ１８５ａのそれぞれとの間でデータ及び信号の送受信を行う。
【0234】
加算平均処理部２２２は、ユーザがキーボード５１１（図３参照）等を用いて入力した所定回数、セグメント幅、及び総加算回数の情報を取得する。総加算回数は、加算平均処理部２２２が磁場データの加算を行う総回数である。なお、加算平均処理部２２２は、ＨＤ５０４等に予め保存された所定回数、セグメント幅及び総加算回数の情報を、ＨＤ５０４等を参照して取得してもよい。
【0235】
加算平均処理部２２２は、計測装置２から通信部２２１を介して、少なくとも１つのトリガー信号を受信する。また加算平均処理部２２２は、少なくとも１つのトリガー信号に対応し、所定回数からセグメント幅分だけ遡って計測装置１８２が計測した磁場データを、通信部２２１を介して受信して加算平均処理する。
【0236】
加算平均処理は、計測装置１８２が計測した磁場データを所定回数からセグメント幅分だけ遡って加算処理した値を加算回数で除算して平均値を算出する処理である。加算平均処理部２２２は、セグメント加算データをセグメント加算データにおける加算回数とセグメントの始点及び終点の加算回数を対応付け、通信部２２１を介してデータ保存サーバ１８５ａに送信する。
【0237】
計測制御部２２３は、解析ＷＳ１８４ａにおけるセグメント加算データの評価結果に基づく指示を、通信部２２１を介して受信し、この推定結果に基づき、計測装置１８２に計測を中断させ、又は延長させることができる。計測制御部２２３は、計測の中断又延長の指示を、指示があった任意のタイミングで割込みデータとして受信可能である。
【0238】
（データ保存サーバ１８５ａの機能構成例）
図３０は、データ保存サーバ１８５ａの機能構成の一例を説明するブロック図である。図３０に示すように、データ保存サーバ１８５ａは、通信部２２５と、保存部２２６とを有する。
図３０は、データ保存サーバ１８５ａの機能構成の一例を説明するブロック図である。図３０に示すように、データ保存サーバ１８５ａは、通信部２２５と、保存部２２６とを有する。
【0239】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図３０は、データ保存サーバ１８５ａが有する主な構成を示すが、データ保存サーバ１８５ａはこれら以外の構成を有してもよい。
【0240】
通信部２２５は、計測ＷＳ１８３ａ及び解析ＷＳ１８４ａのそれぞれとの間でデータ及び信号の送受信を行う。
【0241】
保存部２２６は、通信部２２５を介して計測ＷＳ１８３ａから受信したトリガー回数５２６とエポックデータ５２７とを対応付けて保存する。また保存部２２６は、セグメント加算データのリストも保存する。
【0242】
（解析ＷＳ１８４ａの機能構成例）
次に図３１は、解析ＷＳ１８４ａの機能構成の一例を説明するブロック図である。解析ＷＳ１８４ａは、通信部２２７と、加算平均処理部２２９と、推定部２３０と、表示部２２８と、第５指定受付部２３１と、計測制御部２３２とを有する。
次に図３１は、解析ＷＳ１８４ａの機能構成の一例を説明するブロック図である。解析ＷＳ１８４ａは、通信部２２７と、加算平均処理部２２９と、推定部２３０と、表示部２２８と、第５指定受付部２３１と、計測制御部２３２とを有する。
【0243】
これら各部は、図３に示されている各構成要素の何れかがＲＯＭ５０２からＲＡＭ５０３上に展開されたプログラムに従ったＣＰＵ５０１からの命令によって動作することにより実現される機能又は機能する手段である。なお図３１は、解析ＷＳ１８４ａが有する主な構成を示すが、解析ＷＳ１８４ａはこれら以外の構成を有してもよい。
【0244】
通信部２２７は、計測ＷＳ１８３ａ及びデータ保存サーバ１８５ａのそれぞれとの間でデータ及び信号の送受信を行う。
【0245】
加算平均処理部２２９は、通信部２２７を介して保存部２２６が保存するセグメント加算データのリストを取得し、取得したセグメント加算データのリストを表示部２２８によりディスプレイ５０６等に表示させる。第５指定受付部２３１は、ユーザがセグメント加算データのリストを視認し、選択した結果から加算回数の情報を受け付けることができる。選択された少なくとも一つのセグメント加算データは、それぞれのデータの加算回数に基づいて加算平均処理が行われる。すなわち、提示されたセグメント加算データのリストから異常なセグメント加算データを除外し、加算平均データを取得することができる。
【0246】
推定部２３０は、生体データに基づく処理を行う生体データ処理部の一例である。推定部２３０は、加算平均処理部２２９が生成した磁場データに基づき、生体の活動電流の強さを推定する処理を行う。推定部２３０は、取得した加算平均データを用いて、生体の活動電流の強さを推定する。この推定アルゴリズムには、上述した「Array-Gain Constraint Minimum-Norm Spatial Filter With Recursively Updated Gram Matrix」（例えば、「https://ieeexplore.ieee.org/document/5415622」参照）等を用いることができる。
【0247】
表示部２２８は、推定部２３０による活動電流の強さの推定結果を表示する。例えば表示部２２８は、ディスプレイ５０６に推定結果を表示し、ユーザに視認させることができる。また、表示部２２８は、取得した加算平均データやセグメント加算データを示す波形データ、周波数スペクトルあるいは推定結果等の解析結果を表示することもできる。この一連のデータ取得から表示処理は測定途中であっても行われる。
【0248】
解析ＷＳ１８４ａは、表示部２２８にセグメント加算データを波形データ及び周波数スペクトルとして表示する際には、複数の所定回数に紐づいたセグメント加算データを並列表示させてもよい。解析ＷＳ１８４ａは、表示された波形データを評価し、セグメントにおいて計測が正常に行われた否かを判断し、記録する。
【0249】
計測が正常に行われた否かの判断は、ユーザが目視で行って解析ＷＳ１８４ａに結果を入力してもよいし、解析ＷＳ１８４ａが何らかのアルゴリズムに従って自動的に判断してもよい。例えばユーザは周波数スペクトルを目視し、生体信号以外の周波数の箇所にピークが立っている場合は異常と判断したり、波形データから得られる生体信号の振幅が他のセグメント加算データの振幅と比較して明らかに低下していれば異常と判断したりできる。「明らかに低下」とは、例えば該セグメント加算データに含まれるノイズレベル以上に低下していることをいう。
【0250】
セグメント加算データが異常なデータであると判断された場合には、解析ＷＳ１８４ａは、加算平均データを生成する際にこのセグメント加算データを除外するため、セグメント幅分だけ加算回数を増加させる。解析ＷＳ１８４ａは、ユーザによりキーボード５１１等を用いて指示された計測の延長を受け付ける。計測制御部２３２は、指示に応答して計測装置１８２に計測を延長させることができる。また解析ＷＳ１８４ａは、受け付けた延長の指示を、通信部２２７を介して計測ＷＳ１８３ａに送信し、計測ＷＳ１８３ａが計測制御部２２３によって、計測装置１８２による計測を延長させてもよい。
【0251】
推定部２３０は、測定途中であっても、すでに保存されているセグメント加算データについて活動電流の強さを推定し、ユーザは、表示部２２８が表示した推定部２３０による推定結果を視認できる。
【0252】
また、ユーザは磁場データの数量が不足していると判断した場合には、キーボード５１１等を用いて計測の延長を指示する。計測制御部２３２は、指示に応答して計測装置１８２に計測を延長させることができる。また解析ＷＳ１８４ａは、受け付けた延長の指示を、通信部２２７を介して計測ＷＳ１８３ａに送信し、計測ＷＳ１８３ａが計測制御部２２３によって、計測装置１８２による計測を延長させてもよい。
【0253】
図３２は、生体データ処理装置１８６ａの動作の一例を示すフローチャートである。なお、図３２におけるステップＳ１２１及びステップＳ１２２は、図２３におけるステップＳ１０１、ステップＳ１０２と同じである。図３２におけるステップＳ１２４及びステップＳ１２５は、図２３におけるステップＳ１０３及びステップＳ１０４と同じである。図３２におけるステップＳ１２７からステップＳ１３０は、図２３におけるステップＳ１０６からステップＳ１０９と同じである。以下、図２３との相違点を中心に説明する。
【0254】
ステップＳ１２３において、生体データ処理装置１８６ａは、計測ＷＳ１８３ａにより、通信部２２１を介してエポックデータをデータ保存サーバ１８５ａに送信する。
【0255】
ステップＳ１２６において、生体データ処理装置１８６ａは、ステップＳ１２５にて生成されたセグメント加算データを表示部２２４に送り、表示させる。
【0256】
ステップＳ１３１において、生体データ処理装置１８６ａは、解析ＷＳ１８４ａの第５指定受付部２３１により、通信部２２７を介して保存部２２６が保存する全エポックデータ（以下Ｒａｗデータと呼称）を取得し、このＲａｗデータに含まれるエポックデータのリストを表示部２２８によりディスプレイ５０６等に表示させる。ユーザは、表示されたエポックデータのリストを視認し、キーボード５１１等を用いてエポックデータを選択することができる。このとき、解析ＷＳ１８４ａは、異常セグメントに含まれるエポックデータは除外する。
【0257】
ステップＳ１３２において、解析ＷＳ１８４ａは、選択された少なくとも一つのエポックデータを加算平均処理する。この処理により、解析ＷＳ１８４ａは、少なくとも一つのトリガーに対応した加算平均データを得る。
【0258】
このようにして、解析ＷＳ１８４ａは、活動電流の推定処理に必要な加算平均データを得ることができる。
【0259】
＜各種表示画面例＞
生体データ計測システム１８１により表示される各種表示画面について説明する。
生体データ計測システム１８１により表示される各種表示画面について説明する。
【0260】
（所定回数及び総加算回数の指定画面例）
所定回数とセグメント幅及び総加算回数の設定は、例えばディスプレイに表示される図３３の画面を用いて行われる。
所定回数とセグメント幅及び総加算回数の設定は、例えばディスプレイに表示される図３３の画面を用いて行われる。
【0261】
図３３は、本実施形態に係る所定回数、セグメント幅及び総加算回数の指定画面を例示する図である。図３３では、スイッチ２３３は、加算回数の自動更新を行うようにオン設定されている。自動更新が行われるため、入力ボックス２３４には、加算回数の自動更新間隔（更新間隔回数）として「１００回」が入力されており、入力ボックス２３５にはセグメント幅として「５００回」が入力されている。入力ボックス２３６は、総加算回数を入力するものである。
【0262】
計測ＷＳ１８３ａは、加算回数が所定回数に達するたびに、所定回数に更新間隔回数の「１００回」を加算して所定回数を更新する。図３３では、総加算回数として４０００回が指定されているため、所定回数は、加算回数が所定回数に達するたびに、自動更新により５００回、６００回、７００回、・・・、３７００回、３８００回、３９００回と変化する。
【0263】
入力ボックス２３５はセグメント幅を指定するものである。セグメント幅は入力ボックス２３５を通じて任意のタイミングにより更新されてよい。この時、セグメント加算実行間隔も随時変更されてもよい。また、セグメント幅よりもセグメント加算実行間隔が小さい場合には、最初の所定回数はセグメント幅と同一の値が設定されてもよい。
【0264】
＜生体データ処理装置１８６ａの作用効果＞
以上説明したように、生体データ処理装置１８６ａは、計測ＷＳ１８３ａの加算平均処理部２２２により、複数のトリガー信号に応じて生体における複数の部位で計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理し、表示する。所定回数の間隔はセグメント幅と独立して設定される。
以上説明したように、生体データ処理装置１８６ａは、計測ＷＳ１８３ａの加算平均処理部２２２により、複数のトリガー信号に応じて生体における複数の部位で計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理し、表示する。所定回数の間隔はセグメント幅と独立して設定される。
【0265】
生体データ処理装置１８６ａは、所定回数の間隔をセグメント幅より小さく設定することが可能であり、計測中の最新のデータを第４実施形態よりも早く確認することができる。これにより、生体データ処理装置１８６ａは、計測中に何らかの異常が発生した場合、より早く異常を検出でき、計測の中断や、対策を打ったうえでの計測の延長などの対応を行うことができる。
【0266】
生体データ処理装置１８６ａは、異常が発生した場合には、より早く異常を検出し、対策を取ることで計測の正常性を保証する。また、異常の発生している区間のエポックデータを詳細に除外することができる。
【0267】
［第６実施形態］
第４及び第５実施形態においては、セグメント加算データの評価は計測ＷＳ上で行われていた。計測ＷＳの最も重要な役割は磁場データの収録と記録である。計測ＷＳは、セグメントごとの計測が正常に行われているかを確認しつつも、表示処理などに処理のリソースを割くことにより磁場データの記録漏れなどの悪影響が発生することを回避しなければならない。
第４及び第５実施形態においては、セグメント加算データの評価は計測ＷＳ上で行われていた。計測ＷＳの最も重要な役割は磁場データの収録と記録である。計測ＷＳは、セグメントごとの計測が正常に行われているかを確認しつつも、表示処理などに処理のリソースを割くことにより磁場データの記録漏れなどの悪影響が発生することを回避しなければならない。
【0268】
本実施形態に係る生体データ処理装置１８６ｂは、セグメント加算データの表示機能を解析ＷＳ１８４ｂに集約することで、計測ＷＳ１８３ｂ上において行われる処理を最小限に留めながら第４実施形態と同等の機能を実現する。なお、計測ＷＳ１８３ｂの機能構成は計測ＷＳ１８３ａの機能構成と同じであり、解析ＷＳ１８４ｂの機能構成は解析ＷＳ１８４ａの機能構成と同じであり、データ保存サーバ１８５ｂの機能構成はデータ保存サーバ１８５ａの機能構成と同じであ。
【0269】
＜生体データ処理装置１８６ｂの動作例＞
次に図３４を参照して、生体データ処理装置１８６ｂの動作について説明する。
次に図３４を参照して、生体データ処理装置１８６ｂの動作について説明する。
【0270】
まず、ステップＳ１４１において、生体データ処理装置１８６ｂは、計測ＷＳ１８３ｂの加算平均処理部２２２により、ユーザがキーボード５１１等を用いて入力した所定回数、セグメント幅及び総加算回数の情報を取得する。なお、加算平均処理部２２２は、ＨＤ５０４等に予め保存された所定回数とセグメント幅及び総加算回数の情報をＨＤ５０４等から取得してもよい。
【0271】
続いて、ステップＳ１４２において、生体データ処理装置１８６ｂは、加算平均処理部２２２により、少なくとも一つのトリガー信号と、このトリガー信号に対応して計測された磁場データと、を計測装置１８２から通信部２２１を介して受信して蓄積する。
【0272】
続いて、ステップＳ１４３において、生体データ処理装置１８６ｂは、加算平均処理部２２２により、受信したトリガー信号をカウントしてトリガーごとに加算回数を取得し、加算回数が所定回数に達したか否かを判定する。
【0273】
ステップＳ１４３で達したと判定された場合には（ステップＳ１４３、Ｙｅｓ）、生体データ処理装置１８６ｂは、加算平均処理部２２２により、所定回数からセグメント幅分だけ遡って磁場データを加算平均処理する。一方、達していないと判定された場合には（ステップＳ１４３、Ｎｏ）、生体データ処理装置１８６ｂは、ステップＳ１４２の動作を再度行う。
【0274】
続いて、ステップＳ１４４において、生体データ処理装置１８６ｂは、加算平均処理部２２２により、セグメント加算平均処理結果であるセグメント加算データと、セグメント加算データにおける加算回数の情報とを対応付け、通信部２２１を介してデータ保存サーバ１８５ｂに送信する。データ保存サーバ１８５ｂは、受信したセグメント加算データと、セグメント加算に係る情報と、を対応付けて保存する。セグメント加算データは表示部２２４に送られ、表示される。
【0275】
続いて、ステップＳ１４５において、生体データ処理装置１８６ｂは、計測制御部２２３により、加算回数は最終加算回数（総加算回数）に達したか否かを判定する。なお、この判定は、計測制御部２２３でなく、加算平均処理部２２２が行ってもよい。
【0276】
ステップＳ１４５で、達していないと判定された場合には（ステップＳ１４５、Ｎｏ）、生体データ処理装置１８６ｂは、ステップＳ１４２以降の動作を再度行う。一方、ステップＳ１４５で、達したと判定された場合には（ステップＳ１４５、Ｙｅｓ）、生体データ処理装置１８６ｂは、計測制御部２２３により、計測装置１８２に計測を終了させる。
【0277】
このようにして、生体データ処理装置１８６ｂは、計測ＷＳ１８３ｂにより、加算平均処理を実行し、また計測の中断又は延長の指示に応じた計測装置１８２の制御を行うことができる。
【0278】
続いて、ステップＳ１４６において、生体データ処理装置１８６ｂは、解析ＷＳ１８４ｂの第５指定受付部２３１により、通信部２２７を介して保存部２２６が保存するセグメント加算データのリストを取得する。
【0279】
続いて、ステップＳ１４７において、生体データ処理装置１８６ｂは、取得したセグメント加算データのリストを表示部２２８によりディスプレイ５０６等に表示させる。ユーザは、表示されたセグメント加算データのリストを視認し、任意のセグメント加算データを指定し、表示部２２８に何らかの解析処理を適用した結果である磁場データの波形データもしくは周波数スペクトルや推定結果等を表示させる。
【0280】
続いて、ステップＳ１４８において、生体データ処理装置１８６ｂは、表示部２２８に表示されたセグメント加算データを確認し、正常であるか否かを判断したユーザの指示を受け付ける。
【0281】
ステップＳ１４８において、正常ではないと判断された場合（ステップＳ１４８、Ｎｏ）、ステップＳ１４９において、生体データ処理装置１８６ｂは、セグメントを異常セグメントとして記録し、ステップＳ１５０において最終加算回数の増加、すなわち計測の延長または計測の中断を指示する。この指示は通信部２２７を通じて計測ＷＳ１８３ｂに送られ、即座に指示に応じた処理が実行される。
【0282】
生体データ処理装置１８６ｂは、ステップＳ１４６乃至ステップＳ１５０の処理を最終加算回数に到達するまで繰り返す。
【0283】
一方、ステップＳ１４５で達していると判断された後、ステップＳ１５１において、生体データ処理装置１８６ｂは、第５指定受付部２３１により、通信部２２７を介して保存部２２６が保存するセグメント加算データのリストを取得し、取得したセグメント加算データのリストを表示部２２８によりディスプレイ５０６等に表示させる。このとき、生体データ処理装置１８６ｂは、異常セグメントとして記録されているセグメント加算データを除外する。
【0284】
続いて、ステップＳ１５２において、生体データ処理装置１８６ｂは、選択された少なくとも一つのセグメント加算データをセグメント幅に応じて加算平均処理する。この処理により、少なくとも一つのトリガーに対応した加算平均データを得る。
【0285】
このようにして、生体データ処理装置１８６ｂは、解析ＷＳ１８４ｂにより、生体の活動電流の推定処理に必要な加算平均データを得ることができる。
【0286】
＜生体データ処理装置１８６ｂの作用効果＞
以上説明したように、生体データ処理装置１８６ｂは、加算平均処理部２２２により、複数のトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理を実行する。生体データ処理装置１８６ｂは、処理結果であるセグメント加算データと、このセグメント加算データにおける加算回数とを対応付けて保存部２２６に保存させる。
以上説明したように、生体データ処理装置１８６ｂは、加算平均処理部２２２により、複数のトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびにセグメント加算平均処理を実行する。生体データ処理装置１８６ｂは、処理結果であるセグメント加算データと、このセグメント加算データにおける加算回数とを対応付けて保存部２２６に保存させる。
【0287】
その後、セグメントデータは解析ＷＳ１８４ｂ上において評価され、対応するセグメントにおける計測が正常に実施されたか否かを判断される。計測終了後、生体データ処理装置１８６ｂは、ユーザが異常と判断したセグメントを除外したうえで加算平均を行い、加算平均データを得る。
【0288】
このように、生体データ処理装置１８６ｂは、解析ＷＳ１８４ｂ上においてセグメント加算データを表示、評価及び加算を行うことにより、異常なデータを取り除いて加算平均データを得ることができる。これにより、生体データ処理装置１８６ｂは、計測ＷＳ１８３ｂへの負荷を最小限に留めながら、異常が発生している区間のデータを取り除いたうえで、信号品質を保証することができる。
【0289】
以上に示した実施形態は互いを排さない。
【0290】
以上、実施形態について説明してきたが、本発明は、具体的に開示された上記の実施形態に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。
【0291】
また、実施形態は、プログラムも含む。例えば、プログラムは、少なくとも１つのトリガー信号に応じて計測される生体データの加算回数が所定回数に達するたびに加算平均処理し、少なくとも１つの刺激部位ごとでの加算平均処理結果である加算平均データと、前記加算平均データまたはセグメント加算データと、いずれかのデータにおける前記加算回数と、を対応付けて保存部が保存し、指定された前記加算回数に基づき前記保存部を参照して取得される少なくとも１つの刺激部位ごとでの前記加算平均データを用いて、前記生体データに基づく処理を行う処理をコンピュータに実行させる。このようなプログラムにより、上述した生体データ処理装置と同様の効果を得ることができる。
【0292】
また、上記で説明した実施形態の各機能は、一又は複数の処理回路によって実現することが可能である。ここで、本明細書における「処理回路」とは、電子回路により実装されるプロセッサのようにソフトウェアによって各機能を実行するようプログラミングされたプロセッサや、上記で説明した各機能を実行するよう設計されたASIC(Application Specific Integrated Circuit)、DSP（digital signal processor）、FPGA（field programmable gate array）や従来の回路モジュール等のデバイスを含むものとする。
【符号の説明】
【0293】
１、１８１ 生体データ計測システム
１０、１８６、１８６ａ、１８６ｂ 生体データ処理装置
１００ 被検者（生体の一例）
２、１８２ 計測装置
２００ 磁気センサアレイ
２０１ 磁気センサ
２１０ デュワー
３、１８３、１８３ａ、１８３ｂ 計測ＷＳ
３１、１８７ 通信部
３２、１８８ 加算平均処理部
３３、１８９ 計測制御部
４、１８４、１８４ａ、１８４ｂ 解析ＷＳ
４１、１９３ 通信部
４２、１９５ 推定部（生体データ処理部の一例）
４３ 第１指定受付部
４４ 指示受付部
４５、１９７ 表示部
４６ 判断部
４７ 切替部
４８ 第２指定受付部
４９ 第３指定受付部
５、１８５ データ保存サーバ
５１、１９１ 通信部
５２、１９２ 保存部
５２１ 加算回数
５２２ 加算平均データ
５２４ セグメント加算データ
５２５ セグメント加算に係る情報
５２６ トリガー回数
５２７ エポックデータ
６０ 計測画面
６１ 操作画面
６２ 計測データ画面
６３ 波形データ（磁場データの一例、生体データの一例）
７０ 脊髄位置指定画面
８０、９０、１５０ 推定結果表示画面
１９６ 第４指定受付部
２２１ 通信部
２２２ 加算平均処理部
２２３、２３２ 計測制御部
２２４ 表示部
２２５ 通信部
２２６ 保存部
２２７ 通信部
２２８ 表示部
２２９ 加算平均処理部
２３０ 推定部
２３１ 第５指定受付部
１０、１８６、１８６ａ、１８６ｂ 生体データ処理装置
１００ 被検者（生体の一例）
２、１８２ 計測装置
２００ 磁気センサアレイ
２０１ 磁気センサ
２１０ デュワー
３、１８３、１８３ａ、１８３ｂ 計測ＷＳ
３１、１８７ 通信部
３２、１８８ 加算平均処理部
３３、１８９ 計測制御部
４、１８４、１８４ａ、１８４ｂ 解析ＷＳ
４１、１９３ 通信部
４２、１９５ 推定部（生体データ処理部の一例）
４３ 第１指定受付部
４４ 指示受付部
４５、１９７ 表示部
４６ 判断部
４７ 切替部
４８ 第２指定受付部
４９ 第３指定受付部
５、１８５ データ保存サーバ
５１、１９１ 通信部
５２、１９２ 保存部
５２１ 加算回数
５２２ 加算平均データ
５２４ セグメント加算データ
５２５ セグメント加算に係る情報
５２６ トリガー回数
５２７ エポックデータ
６０ 計測画面
６１ 操作画面
６２ 計測データ画面
６３ 波形データ（磁場データの一例、生体データの一例）
７０ 脊髄位置指定画面
８０、９０、１５０ 推定結果表示画面
１９６ 第４指定受付部
２２１ 通信部
２２２ 加算平均処理部
２２３、２３２ 計測制御部
２２４ 表示部
２２５ 通信部
２２６ 保存部
２２７ 通信部
２２８ 表示部
２２９ 加算平均処理部
２３０ 推定部
２３１ 第５指定受付部
【先行技術文献】
【特許文献】
【0294】
【特許文献1】特許第６５５５８３０号公報
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】
【図34】