特許 7420419 電気インピーダンストモグラフィを用いた血流力学変数算出方法及び装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-01-15

(45)【発行日】2024-01-23

(54)【発明の名称】電気インピーダンストモグラフィを用いた血流力学変数算出方法及び装置

(51)【国際特許分類】

   A61B 5/0536 20210101AFI20240116BHJP

   A61B 5/021 20060101ALI20240116BHJP

   A61B 5/0295 20060101ALI20240116BHJP

【ＦＩ】

A61B5/0536

A61B5/021

A61B5/0295

【請求項の数】 24

(21)【出願番号】P 2022574273

(86)(22)【出願日】2021-04-09

(65)【公表番号】

(43)【公表日】2023-12-05

(86)【国際出願番号】 KR2021004458

(87)【国際公開番号】W WO2021261721

(87)【国際公開日】2021-12-30

【審査請求日】2022-12-01

(31)【優先権主張番号】10-2020-0077738

(32)【優先日】2020-06-25

(33)【優先権主張国・地域又は機関】KR

(73)【特許権者】

【識別番号】519179372

【氏名又は名称】株式会社バイラップ

【氏名又は名称原語表記】ＢＩＬＡＢ ＣＯ．， ＬＴＤ．

(74)【代理人】

【識別番号】100159547

【弁理士】

【氏名又は名称】鶴谷 裕二

(74)【代理人】

【識別番号】100223365

【弁理士】

【氏名又は名称】大町 真義

(72)【発明者】

【氏名】チャン・グギョン

(72)【発明者】

【氏名】ウイ・フン

(72)【発明者】

【氏名】チャン・ジュン

【審査官】右▲高▼ 孝幸

(56)【参考文献】

【文献】特開２０１９－０２５３３８（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１４／００９４６６４（ＵＳ，Ａ１）

【文献】Martin Proenca，Noninvasive pulmonary artery pressure monitoring by EIT: a model-based feasibility study，Medical & biological engineering & computing，2016年09月17日，vol. 55, no. 6, pp. 949-963

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ５／０５

Ａ６１Ｂ ５／０２

(57)【特許請求の範囲】

【請求項1】

ＥＩＴ（Ｅｌｅｃｔｒｉｃａｌ Ｉｍｐｅｄａｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）に基づいて被検者に対する血流力学（ｈｅｍｏｄｙｎａｍｉｃｓ）変数を算出する方法であって、
被検者の胸部分を対象として前記被検者の心周期（ｃａｒｄｉａｃ ｃｙｃｌｅ）内の離散的時点でのＥＩＴ画像を取得するステップと、
前記ＥＩＴ画像における、心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別するステップであって、前記第１画素の値は前記離散的時点にそれぞれ関連し、前記第２画素の値は前記離散的時点にそれぞれ関連している、識別するステップと、
前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップと、を含む、血流力学変数算出方法。

【請求項2】

前記被検者の胸部分を対象として前記被検者の心周期（ｃａｒｄｉａｃ ｃｙｃｌｅ）内の離散的時点でのＥＩＴ画像を取得するステップは、
前記被検者の胸に複数の電極を取り付けて前記被検者の胸部分に対するインピーダンス（ｉｍｐｅｄａｎｃｅ）データを取得するステップと、前記インピーダンスデータから前記ＥＩＴ画像を復元するステップと、を含む、請求項１に記載の血流力学変数算出方法。

【請求項3】

前記ＥＩＴ画像における、心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別するステップであって、前記第１画素の値は前記離散的時点にそれぞれ関連し、前記第２画素の値は前記離散的時点にそれぞれ関連している、識別するステップは、前記ＥＩＴ画像における前記心臓に該当する領域を第１関心領域（ＲＯＩ：ｒｅｇｉｏｎ ｏｆ ｉｎｔｅｒｅｓｔ）として指定し、前記第１関心領域から前記第１画素を選択し、前記ＥＩＴ画像における前記肺に該当する領域を第２関心領域として指定し、前記第２関心領域から前記第２画素を選択するステップを含む、請求項１に記載の血流力学変数算出方法。

【請求項4】

前記第１時点は心室駆出時間（ｖｅｎｔｒｉｃｕｌａｒ ｅｊｅｃｔｉｏｎ ｔｉｍｅ）であり、前記第２時点は肺灌流時間（ｐｕｌｍｏｎａｒｙ ｐｅｒｆｕｓｉｏｎ ｔｉｍｅ）である、請求項１に記載の血流力学変数算出方法。

【請求項5】

前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップを含む、請求項１に記載の血流力学変数算出方法。

【請求項6】

前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値及び前記最大値に基づいて前記肺動脈圧力を算出するステップを含む、請求項５に記載の血流力学変数算出方法。

【請求項7】

前記第２画素の値は、二つ以上のピークパターンを示し、
前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最大値に関連したピークパターンを除く残りのピークパターンにおける前記第２画素のピーク値、及び前記ピーク値に関連した時点に基づいて前記肺動脈圧力を算出するステップを含む、請求項５に記載の血流力学変数算出方法。

【請求項8】

前記第１画素の値は、二つ以上のピークパターンを示し、
前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最小値に関連したピークパターンを除く残りのピークパターンにおける前記第１画素のピーク値、及び前記ピーク値に関連した時点に基づいて前記肺動脈圧力を算出するステップを含む、請求項５に記載の血流力学変数算出方法。

【請求項9】

前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記被検者の心電図（ＥＣＧ：ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｐｈｙ）、前記被検者の血圧（ｂｌｏｏｄ ｐｒｅｓｓｕｒｅ）、前記被検者のフォトプレチスモグラフィ（ＰＰＧ：ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）及び前記被検者の心理心電図（ＳＣＧ：ｓｅｉｓｍｏｃａｒｄｉｏｇｒａｐｈｙ）に基づいて前記肺動脈圧力を算出するステップを含む、請求項５に記載の血流力学変数算出方法。

【請求項10】

前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記算出された肺動脈圧力を用いて肺血管抵抗（ＰＶＲ：ｐｕｌｍｏｎａｒｙ ｖａｓｃｕｌａｒ ｒｅｓｉｓｔａｎｃｅ）を算出するステップをさらに含む、請求項５に記載の血流力学変数算出方法。

【請求項11】

前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちに最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記最小値及び前記第１時点に基づいて心筋収縮度（ｃｏｎｔｒａｃｔｉｌｉｔｙ）を算出するステップを含む、請求項１に記載の血流力学変数算出方法。

【請求項12】

ＥＩＴに基づいて被検者に対する血流力学変数を算出する方法であって、
被検者の胸部分を対象として前記被検者の複数の心周期内の離散的時点でのＥＩＴ画像を取得するステップと、
前記複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び前記複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別するステップであって、前記第１時点は、前記第２時点よりも前記心周期に該当する時間だけ時間的に先立っている、識別するステップと、
前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップと、を含む、血流力学変数算出方法。

【請求項13】

前記第１時点及び前記第２時点は、拡張末期（ｅｎｄ ｄｉａｓｔｏｌｅ）時点であり、
前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて拡張末期心室容積（ＥＤＶＶ：ｅｎｄ－ｄｉａｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するステップを含む、請求項１２に記載の血流力学変数算出方法。

【請求項14】

前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記算出された拡張末期心室容積（ＥＤＶＶ：ｅｎｄ－ｄｉａｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量及び心周期間の一回拍出量（ＳＶ：ｓｔｒｏｋｅ ｖｏｌｕｍｅ）の変化量を用いて心筋収縮度を算出するステップをさらに含む、請求項１３に記載の血流力学変数算出方法。

【請求項15】

前記第１時点及び前記第２時点は、収縮末期（ｅｎｄ ｓｙｓｔｏｌｅ）時点であり、
前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて収縮末期心室容積（ＥＳＶＶ：ｅｎｄ－ｓｙｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するステップを含む、請求項１２に記載の血流力学変数算出方法。

【請求項16】

前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記算出された収縮末期心室容積の変化量及び心周期間の一回拍出量の変化量を用いて駆出率（ＥＦ：ｅｊｅｃｔｉｏｎ ｆｒａｃｔｉｏｎ）を算出するステップをさらに含む、請求項１５に記載の血流力学変数算出方法。

【請求項17】

ＥＩＴに基づいて被検者に対する血流力学変数を算出するための装置であって、
ＥＩＴ画像を記憶する記憶部であって、前記ＥＩＴ画像は、被検者の胸部分を対象として前記被検者の心周期内の離散的時点で取得した画像である、記憶部と、
前記ＥＩＴ画像における心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別し、前記第１画素の値は前記離散的時点にそれぞれ関連し、第２画素の値は前記離散的時点にそれぞれ関連しており、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するように構成された制御部と、を含む、血流力学変数算出装置。

【請求項18】

前記第１時点は心室駆出時間であり、前記第２時点は肺灌流時間である、請求項１７に記載の血流力学変数算出装置。

【請求項19】

前記制御部は、前記第１時点及び前記第２時点に基づいて肺動脈圧力を算出するようにさらに構成される、請求項１７に記載の血流力学変数算出装置。

【請求項20】

前記制御部は、前記第１時点、前記第２時点、前記最小値及び前記最大値に基づいて前記肺動脈圧力を算出するようにさらに構成される、請求項１９に記載の血流力学変数算出装置。

【請求項21】

前記制御部は、前記最小値及び前記第１時点に基づいて心筋収縮度を算出するようにさらに構成される、請求項１７に記載の血流力学変数算出装置。

【請求項22】

ＥＩＴに基づいて被検者に対する血流力学変数を算出するための装置であって、
ＥＩＴ画像を記憶する記憶部であって、前記ＥＩＴ画像は、被検者の胸部分を対象として前記被検者の複数の心周期内の離散的時点で取得された画像である、記憶部と、
前記複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び前記複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別し、前記第１時点は、前記第２時点よりも前記心周期に該当する時間だけ時間的に先立っており、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するように構成された制御部と、を含む、血流力学変数算出装置。

【請求項23】

前記第１時点及び前記第２時点は、拡張末期時点であり、
前記制御部は、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて拡張末期心室容積の変化量を算出するようにさらに構成される、請求項２２に記載の血流力学変数算出装置。

【請求項24】

前記第１時点及び前記第２時点は、収縮末期時点であり、
前記制御部は、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて収縮末期心室容積の変化量を算出するようにさらに構成される、請求項２２に記載の血流力学変数算出装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、医用工学分野に係り、さらに具体的には、被検者に対する血流力学変数を非侵襲的にモニタリングする技術に関する。

【背景技術】

【0002】

肺動脈（ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ）は、心臓の右心室と肺とを繋ぐ動脈であって、人体表面からの接近が難しい。したがって、肺動脈圧力（ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ：ＰＡＰ）を測定するためには、中心静脈を経て心臓の右心房（ｒｉｇｈｔ ａｔｒｉｕｍ）と右心室（ｒｉｇｈｔ ｖｅｎｔｒｉｃｌｅ）を通って肺動脈にカテーテル（ｃａｔｈｅｔｅｒ）を挿入する侵襲的（ｉｎｖａｓｉｖｅ）方法を使用する。この場合は、Ｓｗａｎ－Ｇａｎｚカテーテルを一番多く使用するが、右心房と右心室を通り過ぎるＳｗａｎ－Ｇａｎｚカテーテル手術には熟練した臨床医が必要であり、成功的にカテーテルを挿入しても、心臓に負担を与えて様々な副作用が発生することが知られている。

【0003】

非侵襲的方法でＰＡＰを測定すれば、このような問題を解決することができるので、これについての多くの研究が行われた。特に、心超音波（ｅｃｈｏｃａｒｄｉｏｇｒａｐｈｙ）を用いた多様な方法が提案された。しかし、心超音波は、ユーザが手動で超音波プローブを皮膚に接触させて測定するので、ユーザの熟練度に応じて結果が変動し、同じユーザが繰り返し測定するときの変動も大きい。そして、連続測定が難しいため、モニタリング用途への使用は非常に不便である。

【0004】

かかる問題点を解決するために、電気インピーダンストモグラフィ（ｅｌｅｃｔｒｉｃａｌ ｉｍｐｅｄａｎｃｅ ｔｏｍｏｇｒａｐｈｙ：ＥＩＴ）技術を用いて、ＰＡＰを含む血流力学モニタリング方法についての多くの研究が行われた。人体の胸からＥＩＴデータを測定すると、このデータには、肺の呼吸に対する信号と心臓の血流に対する信号とが混合されている。通常、呼吸信号の大きさは血流信号のｌ０倍以上なので、血流信号を抽出することが非常に難しい。また、血流は呼吸に比べて速く変わるので、ＥＩＴデータの収集速度を上げなければならない。

【0005】

血流力学モニタリングの難題の一つは、非侵襲的にＰＡＰを測定することである。非侵襲的にＰＡＰを測定する幾つかの既存の技術が知られているが、これらの技術はいずれも、ＰＡＰ測定において誤差を発生させると評価されている。血流力学分野で知られている知識によると、ＰＡＰは、心臓の右心室が収縮して血液を肺動脈へ送るのにかかる時間である右心室駆出時間（ｒｉｇｈｔ ｖｅｎｔｒｉｃｕｌａｒ ｅｊｅｃｔｉｏｎ ｔｉｍｅ：ＲＶＥＴ）、肺に血液が供給される肺灌流時間（ｐｕｌｍｏｎａｒｙ ｐｅｒｆｕｓｉｏｎ ｔｉｍｅ：ＰＰＴ）、心室の一回拍出量（ｓｔｒｏｋｅ ｖｏｌｕｍｅ：ＳＶ）、一回拍出量の変化量（ＤＳＶ）などによって決定されるが、既存の技術は、これを考慮していない。

【発明の概要】

【発明が解決しようとする課題】

【0006】

本発明の課題は、電気インピーダンストモグラフィを用いて血流力学変数を非侵襲的にモニタリングする、改善された方法及び装置を提供することである。

【0007】

本発明が解決しようとする課題は上述した課題に制限されず、上述していない別の課題は以降の記載から当業者に明確に理解できるであろう。

【課題を解決するための手段】

【0008】

一態様において、ＥＩＴ（Ｅｌｅｃｔｒｉｃａｌ Ｉｍｐｅｄａｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）に基づいて被検者に対する血流力学（ｈｅｍｏｄｙｎａｍｉｃｓ）変数を算出する方法が提供される。本方法は、被検者の胸部分を対象として前記被検者の心周期（ｃａｒｄｉａｃ ｃｙｃｌｅ）内の離散的時点でのＥＩＴ画像を取得するステップと、前記ＥＩＴ画像における、心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別するステップであって、前記第１画素の値は前記離散的時点にそれぞれ関連し、前記第２画素の値は前記離散的時点にそれぞれ関連している、識別するステップと、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップと、を含むことができる。

【0009】

一実施例において、前記被検者の胸部分を対象として前記被検者の心周期（ｃａｒｄｉａｃ ｃｙｃｌｅ）内の離散的時点でのＥＩＴ画像を取得するステップは、前記被検者の胸に複数の電極を取り付けて前記被検者の胸部分に対するインピーダンス（ｉｍｐｅｄａｎｃｅ）データを取得するステップと、前記インピーダンスデータから前記ＥＩＴ画像を復元するステップと、を含む。

【0010】

一実施例において、前記ＥＩＴ画像における、心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別するステップであって、前記第１画素の値は前記離散的時点にそれぞれ関連し、前記第２画素の値は前記離散的時点にそれぞれ関連している、識別するステップは、前記ＥＩＴ画像における前記心臓に該当する領域を第１関心領域（ＲＯＩ：ｒｅｇｉｏｎ ｏｆ ｉｎｔｅｒｅｓｔ）として指定し、前記第１関心領域から前記第１画素を選択し、前記ＥＩＴ画像における前記肺に該当する領域を第２関心領域として指定し、前記第２関心領域から前記第２画素を選択するステップと、を含む。

【0011】

一実施例において、前記第１時点は心室駆出時間（ｖｅｎｔｒｉｃｕｌａｒ ｅｊｅｃｔｉｏｎ ｔｉｍｅ）であり、前記第２時点は肺灌流時間（ｐｕｌｍｏｎａｒｙ ｐｅｒｆｕｓｉｏｎ ｔｉｍｅ）である。

【0012】

一実施例において、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップを含む。

【0013】

一実施例において、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値及び前記最大値に基づいて前記肺動脈圧力を算出するステップを含む。

【0014】

一実施例において、前記第２画素の値は、二つ以上のピークパターンを示し、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最大値に関連したピークパターンを除く残りのピークパターンにおける前記第２画素のピーク値、及び前記ピーク値に関連した時点に基づいて前記肺動脈圧力を算出するステップを含む。

【0015】

一実施例において、前記第１画素の値は、二つ以上のピークパターンを示し、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最小値に関連したピークパターンを除く残りのピークパターンにおける前記第１画素のピーク値、及び前記ピーク値に関連した時点に基づいて前記肺動脈圧力を算出するステップを含む。

【0016】

一実施例において、前記第１時点及び前記第２時点に基づいて肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するステップは、前記第１時点、前記第２時点、前記被検者の心電図（ＥＣＧ：ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｐｈｙ）、前記被検者の血圧（ｂｌｏｏｄ ｐｒｅｓｓｕｒｅ）、前記被検者のフォトプレチスモグラフィ（ＰＰＧ：ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）及び前記被検者の心理心電図（ＳＣＧ：ｓｅｉｓｍｏｃａｒｄｉｏｇｒａｐｈｙ）に基づいて前記肺動脈圧力を算出するステップを含む。

【0017】

一実施例において、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記算出された肺動脈圧力を用いて肺血管抵抗（ＰＶＲ：ｐｕｌｍｏｎａｒｙ ｖａｓｃｕｌａｒ ｒｅｓｉｓｔａｎｃｅ）を算出するステップをさらに含む。

【0018】

一実施例において、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するステップは、前記最小値及び前記第１時点に基づいて心筋収縮度（ｃｏｎｔｒａｃｔｉｌｉｔｙ）を算出するステップを含む。

【0019】

他の態様において、ＥＩＴに基づいて被検者に対する血流力学変数を算出する方法が提供される。本方法は、被検者の胸部分を対象として前記被検者の複数の心周期内の離散的時点でのＥＩＴ画像を取得するステップと、前記複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び前記複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別するステップであって、前記第１時点は、前記第２時点よりも前記心周期に該当する時間だけ時間的に先立っている、識別するステップと、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップと、を含むことができる。

【0020】

一実施例において、前記第１時点及び前記第２時点は拡張末期（ｅｎｄ ｄｉａｓｔｏｌｅ）時点であり、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて拡張末期心室容積（ＥＤＶＶ：ｅｎｄ－ｄｉａｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するステップを含む。

【0021】

【0022】

一実施例において、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記算出された拡張末期心室容積の変化量及び心周期間の一回拍出量（ＳＶ：ｓｔｒｏｋｅｖｏｌｕｍｅ）の変化量を用いて心筋収縮度を算出するステップをさらに含む。

【0023】

一実施例において、前記第１時点及び前記第２時点は、収縮末期（ｅｎｄ ｓｙｓｔｏｌｅ）時点であり、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて収縮末期心室容積（ＥＳＶＶ：ｅｎｄ－ｓｙｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するステップを含む。

【0024】

一実施例において、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するステップは、前記算出された収縮末期心室容積の変化量及び心周期間の一回拍出量の変化量を用いて駆出率（ＥＦ：ｅｊｅｃｔｉｏｎ ｆｒａｃｔｉｏｎ）を算出するステップをさらに含む。

【0025】

他の態様において、ＥＩＴに基づいて被検者に対する血流力学変数を算出するための装置が提供される。本装置は、ＥＩＴ画像を記憶する記憶部であって、前記ＥＩＴ画像は、被検者の胸部分を対象として前記被検者の心周期内の離散的時点で取得された画像であるである、記憶部と、前記ＥＩＴ画像における、心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別し、前記第１画素の値は前記離散的時点にそれぞれ関連し、前記第２画素の値は前記離散的時点にそれぞれ関連しており、前記識別された第１画素の値のうちの最小値に関連した第１時点、及び前記識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出するように構成された制御部と、を含むことができる。

【0026】

一実施例において、前記第１時点は心室駆出時間であり、前記第２時点は肺灌流時間である。

【0027】

一実施例において、前記制御部は、前記第１時点及び前記第２時点に基づいて肺動脈圧力を算出するようにさらに構成される。

【0028】

一実施例において、前記制御部は、前記第１時点、前記第２時点、前記最小値及び前記最大値に基づいて前記肺動脈圧力を算出するようにさらに構成される。

【0029】

一実施例において、前記制御部は、前記最小値及び前記第１時点に基づいて心筋収縮度を算出するようにさらに構成される。

【0030】

別の態様において、ＥＩＴに基づいて被検者に対する血流力学変数を算出するための装置が提供される。本装置は、ＥＩＴ画像を記憶する記憶部であって、前記ＥＩＴ画像は、被検者の胸部分を対象として前記被検者の複数の心周期内の離散的時点で取得された画像である、記憶部と、前記複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び前記複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別し、前記第１時点は、前記第２時点よりも前記心周期に該当する時間だけ時間的に先立っており、前記第１時点に該当する画素の値、及び前記第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するように構成された制御部と、を含むことができる。

【0031】

一実施例において、前記第１時点及び前記第２時点は、拡張末期時点であり、前記制御部は、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて拡張末期心室容積の変化量を算出するようにさらに構成される。

【0032】

一実施例において、前記第１時点及び前記第２時点は、収縮末期時点であり、前記制御部は、前記第１時点に該当する画素の値及び前記第２時点に該当する画素の値に基づいて収縮末期心室容積の変化量を算出するようにさらに構成される。

【発明の効果】

【0033】

本発明の実施例によれば、電気インピーダンストモグラフィを用いて血流力学変数を非侵襲的にモニタリングすることができる技術的効果がある。

【図面の簡単な説明】

【0034】

【図1】ＥＩＴ（Ｅｌｅｃｔｒｉｃａｌ Ｉｍｐｅｄａｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）に基づいて被検者に対する血流力学（ｈｅｍｏｄｙｎａｍｉｃｓ）変数を算出するための装置のブロック図の一実施例を示す図である。

【0035】

【図2a】被検者の胸部分を対象として取得されたＥＩＴ画像の一実施例を示す図である。

【0036】

【図2b】被検者の胸部分を対象として取得されたＥＩＴ画像の一実施例を示す面である。

【0037】

【図3】被検者の１心周期の間に第１画素の値及び第２画素の値が離散的時点に応じて変わるパターンを例示するための図である。

【0038】

【図4】複数の心周期の間における、心室駆出時間（ＴＤ_Ｈ）、肺灌流時間（ＴＤ_Ｌ）、平均通過時間（ＭＴＴ）及び肺動脈圧力（ＰＡＰ）の時間による変化を示す図である。

【0039】

【図5】被検者の２心周期の間に第１画素の値及び第２画素の値が離散的時点に応じて変わるパターンを例示するための図である。

【0040】

【図6】ＥＩＴに基づいて被検者に対する血流力学変数を算出する方法を説明するためのフローチャートの第１実施例を示す面である。

【0041】

【図7】ＥＩＴに基づいて被検者に対する血流力学変数を算出する方法を説明するためのフローチャートの第２実施例を示す図である。

【発明を実施するための形態】

【0042】

本発明の利点、特徴、及びそれらを達成する方法は、添付図面と共に詳細に後述されている実施例を参照すると明確になるだろう。しかし、本発明は、以下で開示する実施例に限定されるものではなく、互いに異なる多様な形態で実現される。本実施例は、単に本発明の開示を完全たるものにし、本発明の属する技術分野における通常の知識を有する者に発明の範疇を完全に知らせるために提供されるものである。本発明は、請求項の範疇によってのみ定められる。

【0043】

本明細書で使用される用語は、単に特定の実施例を説明するために使用されたもので、本発明を限定するものではない。例えば、単数で表現された構成要素は、文脈上明白に単数のみを意味しない限り、複数の構成要素を含む概念と理解されるべきである。また、本発明の明細書において、「含む」又は「有する」などの用語は、明細書上に記載された特徴、数字、ステップ、動作、構成要素、部分品又はこれらの組み合わせが存在することを指定しようとするものであり、このような用語の使用によって一つ又はそれ以上の他の特徴、数字、ステップ、動作、構成要素、部分品又はこれらの組み合わせの存在又は付加可能性が予め排除されるものではない。また、本明細書に記載された実施例において、「モジュール」或いは「部」は、少なくとも一つの機能や動作を行う機能的部分を意味することができる。

【0044】

なお、別に定義しない限り、技術的或いは科学的用語を含めて、ここで使用される全ての用語は、本発明の属する技術分野における通常の知識を有する者によって一般的に理解されるのと同一の意味を有する。一般に使用される辞典に定義されているような用語は、関連技術の文脈上において有する意味と一致する意味であると解釈されるべきであり、本発明の明細書において明白に定義しない限りは、理想的又は過度に形式的な意味で解釈されない。

【0045】

以下、添付図面を参照して本発明の好適な実施例をより詳細に説明する。但し、以下の説明では、本発明の要旨を不要に不明確にするおそれがある場合、広く知られている機能や構成についての具体的な説明は省略する。

【0046】

図１は、ＥＩＴ（Ｅｌｅｃｔｒｉｃａｌ Ｉｍｐｅｄａｎｃｅ Ｔｏｍｏｇｒａｈｙ）に基づいて被検者に対する血流力学（ｈｅｍｏｄｙｎａｍｉｃｓ）変数を算出するための装置のブロック図の一実施例を示す図である。

【0047】

図１に示すように、装置１００は、入力インターフェース１１０、制御部１２０、記憶部１３０、及びディスプレイ部１４０を含むことができる。記憶部１３０は、被検者の胸部分を対象として取得されたＥＩＴ画像を記憶することができる。一実施例において、記憶されたＥＩＴ画像は、被検者の胸部分を対象として被検者の１心周期（Ｃａｒｄｉａｃ Ｃｙｃｌｅ）内の離散的時点（ｄｉｓｃｒｅｔｅ ｐｏｉｎｔｓ ｉｎ ｔｉｍｅ）で取得された画像であり得る。一実施例において、記憶されたＥＩＴ画像は、被検者の胸部分を対象として被検者の複数の心周期における離散的時点で取得された画像であり得る。記憶されたＥＩＴ画像は、被検者の胸に複数の電極を取り付けて被検者の胸部分に対するインピーダンス（ｉｍｐｅｄａｎｃｅ）データを取得し、取得されたインピーダンスデータを復元することにより得られる。被検者の胸部分を対象として取得されたＥＩＴ画像の実施例を示す図２ａ及び図２ｂを参照すると、記憶されたＥＩＴ画像は、肺（ｌｕｎｇ）を示す画像領域２１０、及び画像（ｈｅａｒｔ）を示す画像領域２２０を含むことができる。被検者の心周期は、被検者に対して取得された心電図（ＥＣＧ：ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｐｈｙ）波形におけるＲ波（Ｒ ｗａｖｅｓ）間の周期を測定することにより決定される。ＥＩＴ画像が取得される離散的時点間の周期は、システム解像度によって決定される。例えば、秒あたり１００枚のＥＩＴ画像が取得されるシステムの場合、システム解像度は１０ｍｓであって、１０ｍｓごとに１枚のＥＩＴ画像が取得されるので、ＥＩＴ画像が取得される離散的時点間の周期は１０ｍｓになり得る。

【0048】

記憶部１３０は、開示された技術の様々な実施例による画像処理を行うことにより得られた中間結果の画像データ、開示された技術の様々な実施例による画像処理を行うことにより得られた結果画像データ、開示された技術の様々な実施例による画像処理及び／又は演算処理を行うのに必要な変数値、開示された技術の様々な実施例によってＥＩＴ画像から選択された画素データ、及び開示された技術の様々な実施例による演算処理を行うことにより得られた演算処理結果値を記憶するために使用できる。様々な実施例において、記憶部１３０は、ＥＩＴ画像をＤＩＣＯＭ（Ｄｉｇｉｔａｌ Ｉｍａｇｉｎｇ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ ｉｎ Ｍｅｄｉｃｉｎｅ）形式又は一般的なイメージファイル形式（ＢＭＰ、ＪＰＥＧ、ＴＩＦＦなど）で記憶することができる。記憶部１３０は、制御部１２０の実現に必要なソフトウェア／ファームウェアなどをさらに記憶することができる。記憶部１３０は、フラッシュメモリタイプ（ｆｌａｓｈ ｍｅｍｏｒｙ ｔｙｐｅ）、ハードディスクタイプ（ｈａｒｄ ｄｉｓｋ ｔｙｐｅ）、マルチメディアカード（ＭｕｌｔｉＭｅｄｉａ Ｃａｒｄ：ＭＭＣ）、カードタイプのメモリ（例えば、ＳＤ（Ｓｅｃｕｒｅ Ｄｉｇｉｔａｌ）カード又はＸＤ（ｅＸｔｒｅａｍ Ｄｉｇｉｔａｌ）カードなど）、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＳＲＡＭ（Ｓｔａｔｉｃ Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ Ｅｒａｓａｂｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＰＲＯＭ（Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）、磁気メモリ、磁気ディスク、及び光ディスクのうちのいずれか一つの記憶媒体で実現できるが、当業者であれば、記憶部１３０の実現形態がこれに限定されないことが分かる。

【0049】

入力インターフェース１１０は、開示された技術の様々な実施例による画像処理及び／又は演算処理を行うためにユーザ命令を入力するためのハードウェア及びソフトウェアモジュールで構成できる。入力インターフェース１１０は、必要な各種命令を制御部１２０に入力するか、或いはＥＩＴ画像取得装置によって取得されたＥＩＴ画像データなどの様々な画像データを記憶部１３０に入力するか、或いはディスプレイされた画像の一部又は全部を指示してこれによる様々な画像処理を行うために有利に使用できる。また、入力インターフェース１１０は、ＥＩＴ画像における特定の領域を関心領域（ＲＯＩ：ｒｅｇｉｏｎ ｏｆ ｉｎｔｅｒｅｓｔ）として指定し、指定された関心領域内の任意の箇所を選択するために有利に使用できる。一実施例において、入力インターフェース１１０は、コンピュータのキーボード（ｋｅｙｂｏａｒｄ）、キーパッド（ｋｅｙｐａｄ）、タッチパッド（ｔｏｕｃｈｐａｄ）、マウス（ｍｏｕｓｅ）などを含むことができるが、入力インターフェースの種類は、これに制限されない。例えば、入力インターフェース１１０は、前述した入力装置を用いて制御可能なグラフィックユーザインターフェース（Ｇｒａｐｈｉｃ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）を含むこともできる。ディスプレイ部１４０は、開示された技術の様々な実施例によって様々な画像及び／又は多様なデータを視覚的に表示するためのものであって、ＬＣＤディスプレイ、ＬＥＤディスプレイ、ＡＭＯＬＥＤディスプレイ、ＣＲＴディスプレイなどの様々なディスプレイ装置を含むことができる。

【0050】

制御部１２０は、複数の離散的時点でそれぞれ取得されたＥＩＴ画像における、心臓に該当する領域での第１画素の値、肺に該当する領域での第２画素の値を識別するように構成できる。一実施例において、制御部１２０は、ＥＩＴ画像における心臓を示す領域２２０を第１関心領域２２５として指定し、第１関心領域２２５から第２画素２２７を選択するように構成できる。一実施例において、制御部１２０は、ＥＩＴ画像における肺を示す領域２１０を第２関心領域２１５として指定し、第２関心領域２１５から第２画素２１７を選択するように構成できる。一実施例において、エッジ検出アルゴリズム（ｅｄｇｅ ｄｅｔｅｃｔｉｏｎ ａｌｇｏｒｉｔｈｍ）、画像分割アルゴリズム（ｉｍａｇｅ ｓｅｇｍｅｎｔａｔｉｏｎ ａｌｇｏｒｉｔｈｍ）などのアルゴリズムを用いて制御部１２０が第１関心領域２２５及び第２関心領域２１５を指定することができるが、第１関心領域２２５及び第２関心領域２１５を指定する方式がこれに制限されないことを認識しなければならない。例えば、ユーザが入力インターフェース１１０を用いて第１関心領域２２５及び第２関心領域２１５を手動で指定することも可能である。ここで、第１画素の値及び第２画素の値はいずれも、ＥＩＴ画像が取得された離散的時点に関連している。このような関連関係によって、例えば、第１画素２２７の１番目の値は、１５ｍｓで取得されたＥＩＴ画像から選択された第１画素２２７の値であり、第１画素２２７の二番目の値は、２５ｍｓで取得されたＥＩＴ画像から選択された第１画素２２７の値であるなど、第１画素の値が取得された時点を知ることができる。同様に、同じ関連関係によって、第２画素２１７の値が取得された時点も知ることができる。

【0051】

図３は、被検者の１心周期の間に第１画素の値及び第２画素の値が離散的時点に応じて変わるパターンを例示するための図である。

【0052】

図３において、縦軸は画素値を示し、横軸は時間、すなわち離散的時点を示す。図３において、下側の画素値の変化は第１画素２２７の値の変化を示し、上側の画素値の変化は第２画素２１７の値の変化を示す。図３では、図示の便宜上、画素値を連続的な波形で表したが、この画素値は、離散的時点における離散的値であると理解されるべきである。図示の実施例では、豚を被検体として取得したＥＩＴ画像における第１画素２２７の値の変化と第２画素２１７の値の変化を、豚の心周期である８００ｍｓの間に示している。図３において、画素の値は、０～１の間の値としてスケーリングされた値である。また、図３では、０ｍｓでの画素値を０と表したが、これは、ＥＣＧ波形でＲ波が発生する時点（心周期の開始時点又は心臓が収縮を始める時点）を０ｍｓと仮定し、このときの第１画素２２７の値と第２画素２１７の値がいずれも基準値０であると仮定した状態で、時間経過に伴って第１画素２２７の値と第２画素２１７の値が基準値に比べてどれほど変化したかを判別するためである。

【0053】

図３の波形を定性的に分析すると、下側の第１画素２２７の変化は、心臓の容積（ｖｏｌｍｅ）、すなわち心臓に残っている血液の量の変化を示すものであるが、第１画素２２７の値がマイナス側に増加することにより心臓が収縮して血液が心臓から肺動脈を経て肺へ流れることを示す。このような過程を心室駆出（ｖｅｎｔｒｉｃｕｌａｒ ｅｊｅｃｔｉｏｎ）と呼び、第１画素２２７の値が最小値となる時点を心室駆出時間（ｖｅｎｔｒｉｃｕｌａｒ ｅｊｅｃｔｉｏｎ ｔｉｍｅ、ＴＤ_Ｈ）又は収縮末期（ｅｎｄ ｓｙｓｔｏｌｅ）時点と呼ぶ。この時点を基準として、それ以後には心臓がそれ以上収縮するのを止め、心臓に再び血液が増える。図示の如く、この時点以後に、第１画素２２７の値は再び０に近づいてプラス値をとる。一方、上側の第２画素２１７の変化は、肺の容積、すなわち肺に残っている血液の量の変化を示すものであるが、第２画素２１７の値がプラス側に増加することにより、血液が心臓から拍出されて肺へ流入することを示す。このような過程を肺灌流（ｐｕｌｍｏｎａｒｙ ｐｅｒｆｕｓｉｏｎ）と呼び、第２画素２１７の値が最大値となる時点、すなわち心臓が収縮して血液が肺動脈を経由して肺へ供給されることにより、肺にある血液の量が最大となる時点を肺灌流時間（ＴＤ_Ｌ）と呼ぶ。図示されている実施例では、第２画素２１７の値が単調関数で増加するが、他の実施例では、第２画素２１７の値が増加してから減少し、再び増加してから減少し、さらに増加してついに最大値に達する方式で複数のピークパターンを示すことができる。第２画素２１７の値が増加してから減少するのは、心臓が収縮して血液が肺動脈を経由して肺に流入する途中で、さらに血液が逆流して再び心臓へ入る現象によって説明できる。図３において、肺灌流時間ＴＤ_Ｌと心室駆出時間ＴＤ_Ｈとの差異を平均通過時間（ＭＴＴ：ｍｅａｎ ｔｒａｎｓｉｔ ｔｉｍｅ）と呼ぶが、この時間は、心臓から血液が出てくることは中断したが、既に出てきた血液が肺に流れ込む時間を意味する。

【0054】

再び図１の制御部１２０の説明に戻り、制御部１２０は、識別された第１画素２２７の値のうちの最小値に関連した第１時点、及び識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学（ｈｅｍｏｄｙｎａｍｉｃ）変数を算出するようにさらに構成できる。ここで、最小値と最大値は、絶対的な基準における最小値と最大値だけでなく、最小値と最大値からある程度の誤差範囲に入ってくる値を全て包括するものと理解及び解釈されるべきである。図３に示すように、第１時点は心室駆出時間（ＴＤ_Ｈ）であり、第２時点は肺灌流時間（ＴＤ_Ｌ）であり得る。一実施例において、制御部は、第１時点及び第２時点に基づいて、肺動脈圧力（ＰＡＰ：ｐｕｌｍｏｎａｒｙ ａｒｔｅｒｙ ｐｒｅｓｓｕｒｅ）を算出するようにさらに構成される。制御部１２０は、第１時点及び第２時点を変数とする関数であって、肺動脈圧力を算出することができる。一実施例において、制御部１２０は、第１時点及び第２時点の線型結合（ｌｉｎｅａｒ ｃｏｍｂｉｎａｔｉｏｎ）で肺動脈圧力を算出することができる。このような実施例において、線型結合のための第１時点及び第２時点に対する係数は、キャリブレーション（ｃａｌｉｂｒａｔｉｏｎ）手順を経て決定される。ＥＩＴ画像を取得するとき、これと同時に侵襲的方法で肺動脈圧力を測定し、最小二乗法（Ｌｅａｓｔ Ｓｑｕａｒｅ Ｍｅｔｈｏｄ）を適用して、線型結合で表現された肺動脈圧力と侵襲的方法で測定された肺動脈圧力との誤差の二乗を最小化する方向に、第１時点及び第２時点に対する係数を算出することができる。一実施例において、制御部１２０は、第２時点から第１時点を差し引いた値（この場合、第２時点に対する係数は１であり、第１時点に対する係数は－１である）、すなわち、平均通過時間（ＭＴＴ）に基づいて肺動脈圧力を算出するように構成できる。複数の心周期の間に心室駆出時間（ＴＤ_Ｈ）、肺灌流時間（ＴＤ_Ｌ）、平均通過時間（ＭＴＴ）及び肺動脈圧力（ＰＡＰ）の時間による変化を示す図４を参照すると、平均通過時間（ＭＴＴ：ＴＤ_Ｌ－ＴＤ_Ｈ）の波形が肺動脈圧力（ＰＡＰ）の波形に近接して追従することを確認することができる。

【0055】

一実施例において、制御部１２０は、第１時点、第２時点、第１画素２２７の値のうちの最小値及び第２画素２１７の値のうちの最大値に基づいて肺動脈圧力を算出するようにさらに構成できる。前述した実施例と同様に、制御部１２０は、第１時点、第２時点、前記最小値及び前記最大値の線型結合で肺動脈圧力を算出することができる。このような実施例においても、最小二乗法を適用して、線型結合で表現された肺動脈圧力と侵襲的方法で測定された肺動脈圧力との誤差の二乗を最小化する方向に、第１時点、第２時点、前記最小値及び前記最大値に対する係数を算出することができる。制御部１２０は、第２画素２１７の値が二つ以上のピークパターンを示す場合、第１時点、第２時点、前記最小値、前記最大値、前記２つ以上のピークパターンのうちの最大値に関連したピークパターンを除く残りのピークパターンにおける第２画素２１７のピーク値、及びこれらのピーク値に関連した時点に基づいて肺動脈圧力を算出するようにさらに構成できる。この場合も、制御部１２０は、前述した実施例と同様に、最小二乗法を適用して、線型結合で表現された肺動脈圧力と侵襲的方法で測定された肺動脈圧力との誤差の二乗を最小化する方向に、第１時点、第２時点、前記最小値、前記最大値、第２画素２１７のピーク値、及びこれらのピーク値に関連した時点に対する係数を算出することができる。制御部１２０は、第１画素２２７の値が二つ以上のピークパターンを示す場合、第１時点、第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最小値に関連したピークパターンを除く残りのピークパターンにおける第１画素２２７のピーク値、及びこれらのピーク値に関連した時点に基づいて肺動脈圧力を算出するようにさらに構成できる。この場合も、制御部１２０は、前述した実施例と同様に、最小二乗法を適用して、線型結合で表現された肺動脈圧力と侵襲的方法で測定された肺動脈圧力との誤差の二乗を最小化する方向に、第１時点、第２時点、前記最小値、前記最大値、第１画素２１７のピーク値、及びこれらのピーク値に関連した時点に対する係数を算出することができる。制御部１２０は、第１時点、第２時点、被検者の心電図、被検者の血圧（ｂｌｏｏｄ ｐｒｅｓｓｕｒｅ）、被検者のフォトプレチスモグラフィ（ＰＰＧ：ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）及び被検者の心理心電図（ＳＣＧ：ｓｅｉｓｍｏｃａｒｄｉｏｇｒａｐｈｙ）に基づいて肺動脈圧力を算出するようにさらに構成できる。肺動脈圧力を算出するにあたり、複数の種類のデータを複合的に反映するほどより正確な肺動脈圧力算出値を得ることが可能となる。 制御部１２０は、算出された肺動脈圧力を用いて肺血管抵抗（ＰＶＲ：ｐｕｌｍｏｎａｒｙ ｖａｓｃｕｌａｒ ｒｅｓｉｓｔａｎｃｅ）を算出するようにさらに構成できる。制御部１２０は、前記最小値及び第１時点、すなわち心室駆出時間（ＴＤ_Ｈ）に基づいて心臓がどれほど収縮することができるかを示す尺度である心筋収縮度（ｃｏｎｔｒａｃｔｉｌｉｔｙ）を算出するようにさらに構成できる。一実施例において、制御部１２０は、前記最小値を第１時点の値で割った値に基づいて心筋収縮度を算出するように構成される。

【0056】

制御部１２０は、複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別するようにさらに構成できる。ここで、第１時点は、第２時点よりも心周期に該当する時間だけ時間的に先立っている。制御部１２０は、第１時点に該当する画素の値及び第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出するようにさらに構成できる。被検者の二心周期の間に第１画素の値及び第２画素の値が離散的時点に応じて変わるパターンを例示する図５に示されているように、一実施例において第１時点及び第２時点は、記号（Ｔ１＿１及びＴ２＿１）で表した拡張末期（ｅｎｄ ｄｉａｓｔｏｌｅ）時点である。周知の如く、拡張末期時点は、心臓に血液が一番多く充填されたときの時点であって、ＥＣＧ波形においてＲ波が発生する時点から約１００ｍｓ以前の時点又はＰ波が終わる時点である。このような実施例において、制御部１２０は、第１時点に該当する画素の値及び第２時点に該当する画素の値に基づいて拡張末期心室容積（ＥＤＶＶ：ｅｎｄ－ｄｉａｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するようにさらに構成される。一実施例において、制御部１２０は、算出された拡張末期心室容積の変化量及び心周期間の一回拍出量（ＳＶ：ｓｔｒｏｋｅ ｖｏｌｕｍｅ）の変化量を用いて心筋収縮度を算出するようにさらに構成される。また、図５に示されているように、一実施例において、第１時点及び第２時点は、記号（Ｔ１＿２及びＴ２＿２）で表した収縮末期（ｅｎｄ ｓｙｓｔｏｌｅ）時点である。前述したように、収縮末期時点は、心室駆出時間（ＴＤ_Ｈ）に該当する時点であって、心臓から血液が一番多く抜け出したときの時間、すなわち心室の収縮が止まる時間を意味する。このような実施例において、制御部１２０は、第１時点に該当する画素の値及び第２時点に該当する画素の値に基づいて収縮末期心室容積（ＥＳＶＶ：ｅｎｄ－ｓｙｓｔｏｌｉｃ ｖｅｎｔｒｉｃｕｌａｒ ｖｏｌｕｍｅ）の変化量を算出するようにさらに構成される。一実施例において、制御部１２０は、算出された収縮末期心室容積の変化量及び心周期間の一回拍出量の変化量を用いて駆出率（ＥＦ：ｅｊｅｃｔｉｏｎ ｆｒａｃｔｉｏｎ）を算出するようにさらに構成される。

【0057】

制御部１２０は、ハードウェア的な面で、特定用途向け集積回路（Ａｐｐｌｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔｓ：ＡＳＩＣｓ）、デジタルシグナルプロセッサ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒｓ：ＤＳＰ）、デジタルシグナルプロセッサデバイス（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｄｅｖｉｃｅｓ：ＤＳＰＤｓ）、プログラマブルロジックデバイス（Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅｓ：ＰＬＤ）、フィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙｓ：ＦＰＧＡｓ）、プロセッサ（ｐｒｏｃｅｓｓｏｒｓ）、制御器（ｃｏｎｔｒｏｌｌｅｒｓ）、マイクロコントローラ（ｍｉｃｒｏ－ｃｏｎｔｒｏｌｌｅｒｓ）及びマイクロプロセッサ（ｍｉｃｒｏｐｒｏｃｅｓｓｏｒｓ）のうちの少なくとも一つを用いて実現できる。また、制御部１２０は、前述したハードウェアプラットフォーム（ｐｌａｔｆｏｒｍ）上で実行可能なファームウェア（ｆｉｒｍｗａｒｅ）／ソフトウェアモジュールで実現できる。この場合、ファームウェア／ソフトウェアモジュールは、適切なプログラム（ｐｒｏｇｒａｍ）言語で書かれた一つ又はそれ以上のソフトウェアアプリケーション（ｓｏｆｔｗａｒｅ ａｐｐｌｉｃａｔｉｏｎｓ）によって実現できる。

【0058】

図６は、ＥＩＴに基づいて被検者へ対する血流力学変数を算出する方法を説明するためのフローチャートの第１実施例を示す図である。

【0059】

図６に示されているように、本方法は、被検者の胸部分を対象として被検者の心周期内の離散的時点でのＥＩＴ画像を取得するステップ（Ｓ６０５）から始まる。一実施例において、ＥＩＴ画像は、被検者の胸に複数の電極を取り付けて被検者の胸部分に対するインピーダンスデータを取得し、取得されたインピーダンスデータを復元することにより得られる。ステップ（Ｓ６１０）では、複数の離散的時点でそれぞれ取得されたＥＩＴ画像における心臓に該当する領域での第１画素の値、及び肺に該当する領域での第２画素の値を識別する。本ステップでは、ＥＩＴ画像における心臓を示す領域２２０を第１関心領域２２５として指定し、第１関心領域２２５から第１画素２２７を選択し、ＥＩＴ画像における肺を示す領域２１０を第２関心領域２１５として指定し、第２関心領域２１５から第２画素２１７を選択することができる。ここで、第１画素の値及び第２画素の値はいずれも、ＥＩＴ画像が取得された離散的時点に関連している。ステップ（Ｓ６１５）では、識別された第１画素の値のうちの最小値に関連した第１時点、及び識別された第２画素の値のうちの最大値に関連した第２時点に基づいて、少なくとも一つの血流力学変数を算出する。一実施例において、第１時点は心室駆出時間であり、第２時点は肺灌流時間である。本ステップでは、第１時点及び前記第２時点に基づいて肺動脈圧力を算出することができる。本ステップでは、第１時点、第２時点、前記最小値及び前記最大値に基づいて肺動脈圧力を算出することができる。本ステップでは、第２画素の値が二つ以上のピークパターンを示す場合、第１時点、第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最大値に関連したピークパターンを除く残りのピークパターンにおける第２画素のピーク値、及びこれらのピーク値に関連した時点に基づいて肺動脈圧力を算出することができる。本ステップでは、第１画素の値が二つ以上のピークパターンを示す場合、第１時点、第２時点、前記最小値、前記最大値、前記二つ以上のピークパターンのうちの前記最小値に関連したピークパターンを除く残りのピークパターンにおける第１画素のピーク値、及びこれらのピーク値に関連した時点に基づいて肺動脈圧力を算出することができる。本ステップでは、第１時点、第２時点、被検者の心電図、被検者の血圧、被検者のフォトプレチスモグラフィ及び被検者の心理心電図（ＳＣＧ）に基づいて肺動脈圧力を算出することができる。本ステップでは、算出された肺動脈圧力を用いて肺血管抵抗を算出することができる。また、本ステップでは、前記最小値及び第１時点に基づいて心筋収縮度を算出することができる。

【0060】

図７は、ＥＩＴに基づいて被検者に対する血流力学変数を算出する方法を説明するためのフローチャートの第２実施例を示す図である。

【0061】

図７に示されているように、本方法は、被検者の胸部分を対象として被検者の複数の心周期内の離散的時点でのＥＩＴ画像を取得するステップ（Ｓ７０５）から始まる。ステップ（Ｓ７１０）では、複数の心周期のうちの第１心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第１時点に該当する画素の値、及び複数の心周期のうちの第２心周期内の離散的時点で取得されたＥＩＴ画像における心臓に該当する領域での画素の値のうち第２時点に該当する画素の値を識別する。ここで、第１時点は、第２時点よりも心周期に該当する時間だけ時間的に先立っている。ステップ（Ｓ７１５）では、第１時点に該当する画素の値、及び第２時点に該当する画素の値に基づいて少なくとも一つの血流力学変数を算出する。一実施例において、第１時点及び第２時点は拡張末期時点であり、このような実施例において、本ステップでは、第１時点に該当する画素の値及び第２時点に該当する画素の値に基づいて拡張末期心室容積の変化量を算出することができる。本ステップでは、算出された拡張末期心室容積の変化量及び心周期間の一回拍出量の変化量を用いて心筋収縮度を算出することができる。一実施例において、第１時点及び第２時点は収縮末期時点であり、このような実施例において、本ステップでは、第１時点に該当する画素の値及び第２時点に該当する画素の値に基づいて収縮末期心室容積の変化量を算出することができる。また、本ステップでは、算出された収縮末期心室容積の変化量及び心周期間の一回拍出量の変化量を用いて駆出率を算出することができる。

【0062】

以上の説明において、ある構成要素が他の構成要素に接続又は結合されるという記載の意味は、当該構成要素が他の構成要素に直接接続又は結合されるという意味だけでなく、それらがその間に介在した一つ又はそれ以上の別の構成要素を介して接続又は結合されることができるという意味を含むと理解されるべきである。その他にも、構成要素間の関係を記述するための用語（例えば、「の間に」など）も類似な意味で解釈されるべきである。

【0063】

本明細書に開示された実施例において、図示されている構成要素の配置は、発明が実現される環境又は要求事項によって変わり得る。例えば、一部の構成要素が省略されるか、或いは幾つかの構成要素が一つに統合されて実施され得る。また、幾つかの構成要素の配置順序及び連結が変更できる。

【0064】

以上では、本発明の様々な実施例について図示及び説明したが、本発明は、上述した特定の実施例に限定されず、上述した実施例は、添付する特許請求の範囲で請求する本発明の要旨から外れることなく、当該発明の属する技術分野における通常の知識を有する者によって様々に変形実施できるのはもとより、それらの変形実施例が本発明の技術的思想や範囲とは別個に理解されてはならない。よって、本発明の技術的範囲は、添付する特許請求の範囲のみによって定められるべきである。

【産業上の利用可能性】

【0065】

本発明は、電気インピーダンストモグラフィを用いて血流力学変数を非侵襲的にモニタリングすることができる。

【図1】

【図2a】

【図2b】

【図3】

【図4】

【図5】

【図6】

【図7】