JP6861414B2

JP6861414B2 - 解析マップ作成装置およびプログラム

Info

Publication number: JP6861414B2
Application number: JP2017070514A
Authority: JP
Inventors: 一郎佐久間; 直輝富井; 宏瀬野; 晴朗本荘; 正俊山▲崎▼
Original assignee: University of Tokyo NUC; Tokai National Higher Education and Research System NUC
Current assignee: University of Tokyo NUC; Tokai National Higher Education and Research System NUC
Priority date: 2017-03-31
Filing date: 2017-03-31
Publication date: 2021-04-21
Anticipated expiration: 2037-03-31
Also published as: US11547340B2; WO2018180796A1; JP2018171194A; US20200037908A1

Description

本発明は、解析マップ作成装置およびプログラムに関し、詳しくは、興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置およびコンピュータをこうした解析マップ作成装置として機能させるプログラムに関する。

この種の技術として、発明者らは、心臓における興奮波の位相マップを作成すると共に位相マップに基づいて各位置（各箇所）における周囲の位相の分散の程度を値０から値１までの数値による位相分散値で示す位相分散マップを作成するものを提案している（非特許文献１参照）。位相分散値は、周囲の位相の分散の程度が大きいほど大きな値となるから、旋回興奮波の旋回中心では大きな値となる。このため、位相分散値が値１に近い閾値（例えば０．８や０．９）を超える位置（箇所）を旋回興奮波の旋回中心として認定することができる。

"Detection Algorithm of Phase Singularity Using Phase Variance Analysis for Epicardial Optical Mapping Data", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 63, NO. 9, SEPTEMBER 2016

このように心臓における興奮波の位相マップに対して位相分散値を演算することにより旋回興奮波の旋回中心を認定することができるが、旋回興奮波の旋回中心は時間の経過に伴って出現したり消滅したり移動したりする。このため、旋回興奮波の旋回中心の時間経過に伴う出現や消滅、移動を解析することが重要な課題となる。

本発明の解析マップ作成装置は、旋回興奮波に関する解析を行なうための解析マップを作成することを主目的とする。また、本発明のプログラムは、コンピュータを旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置として機能させることを主目的とする。

本発明の解析マップ作成装置およびプログラムは、上述の主目的を達成するために以下の手段を採った。

本発明の解析マップ作成装置は、
興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置であって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算する位相分散値演算部と、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成する解析マップ作成部と、
を備えることを要旨とする。

この本発明の解析マップ作成装置では、生体組織の各位置における興奮波の位相に基づいて各位置における周囲の位相の分散の程度を示す位相分散値を演算し、各位置における位相分散値の少なくとも一部の時系列に基づいて解析マップを作成する。位相分散値は、上述したように、周囲の位相の分散の程度が大きいほど大きな値となるから、閾値（例えば０．８や０．９）を超える位置（箇所）を旋回興奮波の旋回中心として認定することができる。閾値を超えないが閾値近傍の位相分散値を有する位置（箇所）は、その時刻より前に旋回中心であったか、或いは、その時刻より後に旋回中心になり得るか、の可能性を見いだすことができる。このため、生体組織の各位置における位相分散値の少なくとも一部の時系列に基づく解析マップは、旋回興奮波に関する解析を行なうための解析マップとなる。これらの結果、旋回興奮波に関する解析を行なうための解析マップを作成することができる。

こうした本発明の解析マップ作成装置において、前記解析マップ作成部は、前記各位置における位相分散値の時間累積に基づいて前記解析マップを作成するものとしてもよい。ここで、位相分散値の時間累積に基づく解析マップとしては、位相分散値を時間累積したマップや、位相分散値を時間平均したマップなどを挙げることができる。こうした解析マップは、旋回興奮波の旋回中心や旋回中心となり得る可能性の高い位置（箇所）の認定などに用いることができる。

また、本発明の解析マップ作成装置において、前記解析マップ作成部は、前記各位置における位相分散値のうち所定値以上の位相分散値に対する時間累積に基づいて前記解析マップを作成するものとしてもよい。所定値以上の位相分散値に対する時間累積に基づく解析マップとしては、所定値以上の位相分散値だけを時間累積したマップや、所定値以上の位相分散値だけを時間累積した後に更に所定値で２値化したマップ、所定値以上の位相分散値を時間平均したマップなどを挙げることができる。こうした解析マップは、旋回興奮波の旋回中心のその位置（箇所）における滞留の程度の認定や旋回興奮波の旋回中心の軌跡の認定などに用いることができる。

本発明の解析マップ作成装置において、前記解析マップ作成部は、前記各位置における位相分散値を所定値によって２値化した値の時間累積に基づいて前記解析マップを作成するものとしてもよい。所定値によって２値化した値の時間累積に基づく解析マップとしては、所定値によって２値化した値を時間累積したマップや、所定値によって２値化した値を時間累積した後に更に２値化したマップなどを挙げることができる。こうした解析マップは、旋回興奮波の旋回中心のその位置（箇所）における滞留の程度の認定や旋回興奮波の旋回中心の軌跡の認定などに用いることができる。

本発明の解析マップ作成装置において、前記解析マップ作成部は、前記各位置における位相分散値のうち所定値未満の位相分散値に対する時間累積に基づいて前記解析マップを作成するものとしてもよい。所定値未満の位相分散値に対する時間累積に基づく解析マップとしては、所定値未満の位相分散値を時間累積したマップや、所定値未満の位相分散値を時間平均したマップなどを挙げることができる。こうした解析マップは、旋回興奮波の旋回中心になり易い位置（箇所）の認定や、旋回興奮波の旋回中心にはなり得ない位置（箇所）の認定などに用いることができる。

本発明の解析マップ作成装置において、前記位相分散値は、次式（１）により計算されるものとしてもよい。

本発明のプログラムは、
コンピュータを、興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置として機能させるためのプログラムであって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算するステップと、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成するステップと、
を有することを要旨とする。

この本発明のプログラムでは、コンピュータを、生体組織の各位置における興奮波の位相に基づいて各位置における周囲の位相の分散の程度を示す位相分散値を演算し、各位置における位相分散値の少なくとも一部の時系列に基づいて解析マップを作成する装置として機能させることができる。位相分散値は、上述したように、周囲の位相の分散の程度が大きいほど大きな値となるから、閾値（例えば０．８や０．９）を超える位置（箇所）を旋回興奮波の旋回中心として認定することができる。閾値を超えないが閾値近傍の位相分散値を有する位置（箇所）は、その時刻より前に旋回中心であったか、或いは、その時刻より後に旋回中心になり得るか、の可能性を見いだすことができる。このため、生体組織の各位置における位相分散値の少なくとも一部の時系列に基づく解析マップを用いて解析を行なうことにより、旋回興奮波の旋回中心の時間経過に伴う出現や消滅、移動などの解析をより適正に行なうことができる。

本発明の一実施形態としての解析マップ作成装置２０の構成の概略を示す構成図である。本発明の他の実施形態としての解析マップ作成装置１２０の構成の概略を示す構成図である。位相分散マップの作成処理の一例を示すフローチャートである。位相を濃淡で示す位相マップの一例を示す説明図である。図４中の位置Ａの周囲の位相値の平均ベクトルで示す説明図である。図４中の位置Ｂの周囲の位相値の平均ベクトルで示す説明図である。図４の位相マップに対する位相分散マップを示す説明図である。時間累積位相分散マップの作成処理の一例を示すフローチャートである。旋回中心軌跡マップの作成処理の一例を示すフローチャートである。所定時間間隔の時間Ｔ１〜Ｔ４における位相マップと位相分散マップと位相分散値を閾値で２値化したマップの一例を示す説明図である。所定時間間隔の時間Ｔ５〜Ｔ８における位相マップと位相分散マップと位相分散値を閾値で２値化したマップの一例を示す説明図である。図１０および図１１に示した位相分散マップに対する時間累積位相分散マップを示す説明図である。図１０および図１１に示した位相分散値を閾値で２値化したマップに対する旋回中心軌跡マップを示す説明図である。旋回中心移動予測用マップの作成処理の一例を示すフローチャートである。図１０および図１１に示した位相分散マップに対する旋回中心移動予測用マップを示す説明図である。変形例の旋回中心軌跡マップの作成処理の一例を示すフローチャートである。図１０および図１１に示した位相分散マップに対する変形例の旋回中心軌跡マップを示す説明図である。

次に、本発明を実施するための形態を実施例を用いて説明する。図１は、本発明の一実施形態としての解析マップ作成装置２０の構成の概略を示す構成図であり、図２は、図１の解析マップ作成装置２０としても機能する他の実施形態としての解析マップ作成装置１２０の構成の概略を示す構成図である。

図１に示す実施形態の解析マップ作成装置２０は、汎用のコンピュータを興奮波が組織内を伝播することにより興奮に伴って活動する生体組織（例えば心臓など）における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置として機能させるためのプログラムを、汎用のコンピュータに組み込むことによって構成されており、生体組織の各位置（箇所）の興奮波の位相を入力して各位置（箇所）の位相分散値を演算する位相分散演算部３０と、位相分散演算部３０により演算された各位置（箇所）の位相分散値の時系列に基づいて解析マップを作成する解析マップ作成部４０と、を備える。解析マップ作成部４０により作成された解析マップは、汎用のコンピュータに接続された（一体の）ディスプレイに表示出力される。

図２に示す実施形態の解析マップ作成装置１２０は、図１に示す実施形態の解析マップ作成装置２０と同様に、汎用のコンピュータを興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置として機能させるためのプログラムを、汎用のコンピュータに組み込むことによって構成されており、生体組織の各位置（箇所）の興奮波のデータを入力して各位置（箇所）の興奮波の位相による位相マップを作成する位相マップ作成部１２２と、位相マップ作成部１２２により作成された位相マップにおける各位置（箇所）の位相分散値を演算して位相分散マップを作成する位相分散マップ作成部１３０と、位相分散マップ作成部１３０により作成された位相分散マップにおける各位置（箇所）の位相分散値の時系列に基づいて解析マップを作成する解析マップ作成部１４０と、を備える。解析マップ作成部１４０により作成された解析マップは、汎用のコンピュータに接続された（一体の）ディスプレイに表示出力される。

図１に示す実施形態の解析マップ作成装置２０と図２に示す実施形態の解析マップ作成装置１２０とは、基本的には、位相分散マップを作成するか否かの相違があるだけで、基本的な処理については同様である。以下、説明の容易のために、図２に示す実施形態の解析マップ作成装置１２０を用いて説明し、必要に応じて図１に示す実施形態の解析マップ作成装置２０について説明する。

位相マップ作成部１２２は、生体組織の各位置（箇所）の興奮波のデータを入力し、同一時刻の各位置における興奮波の位相を示す位相マップを所定時間毎に作成する。生体組織の各位置（箇所）の興奮波のデータとしては、例えば、心臓の各位置に取り付けられたセンサからの電気信号や、心臓の拍動の画像データなどを用いることができる。なお、位相は、興奮波の周期に対して−π〜πを当てはめている。

位相分散マップ作成部１３０は、例えば図３に示す位相分散マップの作成処理により、位相マップ作成部１２２により作成された所定時間毎の位相マップにおける各位置（箇所）の位相を用いて位相分散値を演算して所定時間毎の位相分散マップを作成する。

位相分散マップの作成処理では、まず、位相マップを入力し（ステップＳ１００）、入力した位相マップにおける各位置（箇所）における位相分散値を次式（１）により演算する（ステップＳ１１０）。式（１）は、対象の位置の周囲のＮ個の位相値の平均ベクトルのノルムを値１から減ずるものであり、Ｖは位相分散値、Ｎは周囲の位置数、θ_kは位置ｋの位相である。図４は位相を濃淡で示す位相マップの一例を示す説明図であり、図５は図４中の位置Ａの周囲の位相値の平均ベクトルで示す説明図であり、図６は図４中の位置Ｂの周囲の位相値の平均ベクトルを示す説明図である。図５および図６の太矢印は周囲の位相値の平均ベクトルを示す。位置Ａでは、図４および図５から解るように、周囲の位相値は−πからπまで略均等に分散している。このため、平均ベクトルのノルムは値０に近いものとなり、位相分散値Ｖは値１に近いものとなる。一方、位置Ｂでは、図４および図６から解るように、周囲の位相値は偏っている。このため、平均ベクトルのノルムは値１に近いものとなり、位相分散値Ｖは値０に近いものとなる。図７に図４の位相マップに対する位相分散マップを示す。こうした位相分散マップにより、旋回興奮波の旋回中心を認定することができる。なお、実施例では、位相分散値Ｖは、対象の位置に対して一定範囲の窓（図４および図７中位置Ａ，Ｂ近傍の白色の正方形）を設定し、窓内の位置の位相により演算した。

位相マップにおける各位置の位相分散値を演算すると、各位置の位相分散値を位相分散マップとして記憶し（ステップＳ１２０）、次の位相マップ（所定時間後の位相マップ）があるか否かを判定し（ステップＳ１３０）、次の位相マップがある場合にはステップＳ１００の位相マップの入力に戻り、次の位相マップがない場合には、位相分散マップの作成処理を終了する。

解析マップ作成部１４０は、位相分散マップ作成部１３０により作成された所定時間毎の位相分散マップにおける各位置（箇所）の位相分散値の時系列に基づいて解析マップを作成する。例えば、解析マップとして、所定時間毎の位相分散マップにおける各位置（箇所）の位相分散値の時間累積のマップを作成する場合には図８に例示する時間累積位相分散マップの作成処理を実行し、旋回興奮波の旋回中心の軌跡を示すマップを作成する場合には図９に例示する旋回中心軌跡マップの作成処理を実行する。

図８の時間累積位相分散マップの作成処理では、まず、位相分散マップを入力し（ステップＳ２００）、位相分算マップの各位置（箇所）の位相分散値を累積し（ステップＳ２１０）、次の位相分散マップ（所定時間後の位相分散マップ）があるか否かを判定する（ステップＳ２２０）。次の位相分散マップがある場合にはステップＳ２００の位相分散マップの入力に戻り、次の位相分散マップがない場合には、各位置（箇所）の位相分散値の累積値を時間累積位相分散マップとして記憶し（ステップＳ２３０）、時間累積位相分散マップの作成処理を終了する。

図９の旋回中心軌跡マップの作成処理では、まず、位相分散マップを入力し（ステップＳ３００）、各位置（箇所）の位相分散値を閾値（例えば、０．８や０．９）より大きいか否かで２値化する（ステップＳ３１０）。即ち、旋回興奮波の旋回中心を認定するのである。そして、各位置（箇所）の２値化した値を累積し（ステップＳ３２０）、次の位相分散ップ（所定時間後の位相分散マップ）があるか否かを判定する（ステップＳ３３０）。次の位相分散マップがある場合にはステップＳ３００の位相分散マップの入力に戻り、次の位相分散マップがない場合には、各位置（箇所）の２値化した値の累積値を旋回中心軌跡マップとして記憶し（ステップＳ３４０）、旋回中心軌跡マップの作成処理を終了する。

所定時間間隔の時間Ｔ１〜Ｔ８における位相マップと位相分散マップと位相分散値を閾値で２値化したマップとを図１０および図１１に示す。図１０および図１１に示した位相分散マップに対する時間累積位相分散マップを図１２に示し、図１０および図１１に示した位相分散値を閾値で２値化したマップに対する旋回中心軌跡マップを図１３に示す。図１２では、位相分散値の時間累積値（白さの程度）は旋回興奮波の旋回中心の発生しやすさの程度を示すと考えられる。したがって、この図１２のマップは、旋回興奮波の旋回中心や旋回中心となり得る可能性の高い位置（箇所）の認定などに用いることができる。図１３では、２値化した値の時間累積値（白さの程度）は旋回興奮波の旋回中心の滞留の程度を示すと考えられ、２値化した値の時間累積値の分布（白い位置の分布）は旋回興奮波の旋回中心の移動を示すと考えられる。したがって、この図１３のマップは、旋回興奮波の旋回中心のその位置（箇所）における滞留の程度の認定や旋回興奮波の旋回中心の軌跡の認定などに用いることができる。

また、解析マップ作成部１４０では、所定時間毎の位相分散マップにおける各位置（箇所）の位相分散値が閾値（例えば０．８や０．９）未満の値の時間累積のマップを旋回中心移動予測用マップとして作成するものとしてもよい。この場合、図１４に例示する旋回中心移動予測用マップの作成処理を実行すればよい。旋回中心移動予測用マップの作成処理では、まず、位相分散マップを入力し（ステップＳ４００）、位相分算マップの各位置（箇所）の位相分散値が閾値（例えば０．８や０．９）以上の位相分散値を値０とする（ステップＳ４１０）。そして、各位置（箇所）の位相分散値を累積し（ステップＳ４２０）、次の位相分散ップ（所定時間後の位相分散マップ）があるか否かを判定する（ステップＳ４３０）。次の位相分散マップがある場合にはステップＳ４００の位相分散マップの入力に戻り、次の位相分散マップがない場合には、各位置（箇所）の位相分散値の累積値を旋回中心移動予測用マップとして記憶し（ステップＳ４４０）、旋回中心移動予測用マップの作成処理を終了する。図１０および図１１に示した位相分散マップに対する旋回中心移動予測用マップを図１５に示す。図１５では、位相分散値の時間累積値（白さの程度）は旋回興奮波の旋回中心が移動する可能性の程度を示すと考えられる。したがって、この図１５のマップは、旋回興奮波の旋回中心になり易い位置（箇所）の認定や、旋回興奮波の旋回中心にはなり得ない位置（箇所）の認定などに用いることができる。

以上説明した実施形態の解析マップ作成装置１２０では、生体組織の各位置（箇所）の興奮波の所定時間毎の位相マップにおける各位置（箇所）の位相分散値を演算して位相分散マップを作成し、位相分散マップにおける各位置（箇所）の位相分散値の時系列に基づく解析マップを作成する。こうして作成された解析マップは、旋回興奮波の旋回中心の発生や消滅、移動などの解析や、旋回興奮波の旋回中心の発生の可能性の程度の解析などに用いることができる。これらの結果、旋回興奮波の旋回中心の時間経過に伴う出現や消滅、移動などの解析に資することができる。

実施形態の解析マップ作成装置１２０では、図８の時間累積位相分散マップの作成処理により各位置（箇所）の位相分散値の時間累積のマップを作成するものとしたが、各位置（箇所）の位相分散値の時間平均のマップを作成するものとしてもよい。この場合のマップは、図１２と同様になる。また、図１４の旋回中心移動予測用マップの作成処理により各位置（箇所）の位相分散値が閾値未満の値の時間累積のマップを作成するものとしたが、各位置（箇所）の位相分散値が閾値未満の値の時間平均のマップを作成するものとしてもよい。この場合のマップは、図１５と同様になる。

実施形態の解析マップ作成装置１２０では、旋回中心軌跡マップについては、各位置（箇所）の位相分散値を閾値（例えば０．８や０．９）で２値化した値を累積することにより作成するものとしたが、２値化した値を累積する際に値１を上限として累積するものとしてもよい。この場合、作成される旋回中心軌跡マップは、図１３のような濃淡のないマップとなる。また、各位置（箇所）の位相分散値のうち閾値（例えば０．８や０．９）以上の位相分散値だけを累積することにより作成するものとしてもよい。この場合、図１６に例示する旋回中心軌跡マップの作成処理を実行すればよい。図１６の旋回中心軌跡マップの作成処理では、まず、位相分散マップを入力し（ステップＳ５００）、各位置（箇所）の位相分散値のうち閾値（例えば、０．８や０．９）未満の位相分散値を値０とする（ステップＳ５１０）。そして、各位置（箇所）の位相分散値を累積し（ステップＳ５２０）、次の位相分散ップ（所定時間後の位相分散マップ）があるか否かを判定する（ステップＳ５３０）。次の位相分散マップがある場合にはステップＳ５００の位相分散マップの入力に戻り、次の位相分散マップがない場合には、各位置（箇所）の位相分散値の累積値を旋回中心軌跡マップとして記憶し（ステップＳ５４０）、旋回中心軌跡マップの作成処理を終了する。この処理による図１０および図１１に示した位相分散マップに対する旋回中心軌跡マップを図１７に示す。図１７は図１３とほぼ同様である。各位置（箇所）の位相分散値のうち閾値以上の位相分散値だけの時間累積のマップを作成するのに代えて、各位置（箇所）の位相分散値のうち閾値以上の位相分散値だけの時間平均のマップを作成するものとしてもよい。この場合のマップは、図１７と同様になる。

また、旋回中心軌跡マップの作成手法としては、所定時間毎の位相マップの各位置（箇所）の位相分散値のうち閾値以上となる位置（箇所）のカウント数をその位置と共に示すことによって旋回中心軌道マップを作成するものとしてもよい。この場合、旋回中心軌道マップは、位相分散マップを作成せずに作成することになるため、図２に示す解析マップ作成装置１２０ではなく、図１に示す解析マップ作成装置２０によって作成されることになる。このように、位相分散マップを作成せずに解析マップを作成するものとしてもよい。

以上、本発明を実施するための形態について実施例を用いて説明したが、本発明はこうした実施例に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において、種々なる形態で実施し得ることは勿論である。

本発明は、旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置の製造産業などに利用可能である。

２０，１２０解析マップ作成装置、３０位相分散演算部、４０，１４０解析マップ作成部、１２２位相マップ作成部、１３０位相分散マップ作成部。

Claims

興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置であって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算する位相分散値演算部と、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成する解析マップ作成部と、
を備え、
前記解析マップ作成部は、前記各位置における位相分散値の時間累積に基づいて前記解析マップを作成する、
解析マップ作成装置。
興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置であって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算する位相分散値演算部と、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成する解析マップ作成部と、
を備え、
前記解析マップ作成部は、前記各位置における位相分散値のうち所定値以上の位相分散値に対する時間累積に基づいて前記解析マップを作成する、
解析マップ作成装置。
興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置であって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算する位相分散値演算部と、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成する解析マップ作成部と、
を備え、
前記解析マップ作成部は、前記各位置における位相分散値を所定値によって２値化した値の時間累積に基づいて前記解析マップを作成する、
解析マップ作成装置。
興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置であって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算する位相分散値演算部と、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成する解析マップ作成部と、
を備え、
前記解析マップ作成部は、前記各位置における位相分散値のうち所定値未満の位相分散値に対する時間累積に基づいて前記解析マップを作成する、
解析マップ作成装置。
請求項１ないし４のうちのいずれか１つの請求項に記載の解析マップ作成装置であって、
前記位相分散値は、式（１）により計算される、
解析マップ作成装置。
コンピュータを、興奮波が組織内を伝播することにより興奮に伴って活動する生体組織における旋回興奮波に関する解析を行なうための解析マップを作成する解析マップ作成装置として機能させるためのプログラムであって、
前記生体組織の各位置における興奮波の位相に基づいて前記各位置における周囲の位相の分散の程度を示す位相分散値を演算するステップと、
前記各位置における位相分散値の少なくとも一部の時系列に基づいて前記解析マップを作成するステップと、
を有し、
前記解析マップを作成するステップは、前記各位置における位相分散値の時間累積に基づいて前記解析マップを作成するステップである、
プログラム。