JP2023545859A - 電気インピーダンスイメージング方法、システム、記憶媒体および電子機器 - Google Patents

電気インピーダンスイメージング方法、システム、記憶媒体および電子機器 Download PDF

Info

Publication number
JP2023545859A
JP2023545859A JP2023524710A JP2023524710A JP2023545859A JP 2023545859 A JP2023545859 A JP 2023545859A JP 2023524710 A JP2023524710 A JP 2023524710A JP 2023524710 A JP2023524710 A JP 2023524710A JP 2023545859 A JP2023545859 A JP 2023545859A
Authority
JP
Japan
Prior art keywords
image
matrix
electrical impedance
covariance matrix
measurement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2023524710A
Other languages
English (en)
Other versions
JP7440985B2 (ja
Inventor
可 張
▲キン▼ 張
明涛 管
誼冰 王
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Original Assignee
Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Huarui Boshi Medical Imaging Technology Co Ltd filed Critical Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Publication of JP2023545859A publication Critical patent/JP2023545859A/ja
Application granted granted Critical
Publication of JP7440985B2 publication Critical patent/JP7440985B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0536Impedance imaging, e.g. by tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physiology (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Pulmonology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

電気インピーダンスイメージング方法、システム、記憶媒体および電子機器である。当該方法は、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得する段階(S110)と、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階(S120)と、複数の測定時点に対応する瞬時微分画像に基づいて、画像行列を構築し、当該画像行列中の各列は、瞬時微分画像に対応するベクトルである段階(S130)と、画像行列に基づいて、当該画像行列に対応する共分散行列を決定する段階(S140)と、共分散行列に従って、当該共分散行列の重みベクトルを獲得する段階(S150)と、共分散行列および重みベクトルに基づいて、測定された人体領域の電気インピーダンス状態画像を獲得する段階(S160)とを含む。【選択図】図1

Description

本出願は、2020年10月23日に出願された「電気インピーダンスイメージング方法、システム、記憶媒体および電子機器」という名称の中国特許出願CN202011145265.2の優先権を主張し、その内容全体は、参照により本出願に組み込まれる。
本開示は、電気インピーダンスイメージングの技術分野に属し、特に電気インピーダンスイメージング方法、システム、記憶媒体および電子機器に関する。
現在の電気インピーダンスイメージング技術は、瞬間的な換気または血液灌流によって引き起こされる電気インピーダンス変化画像しか生成できない。しかし、血液灌流画像のような瞬間的な電気インピーダンス変化画像は、表示されると急速に変化し、観察者が画像に対する理解に役立たない。また、瞬間的な電気インピーダンス変化画像は、異なる人間の生理学的サイクル(例えば、心周期)で大きな違いを示すことが多く、従って、観察者が測定された人体領域の機能(例えば、換気または血液灌流)の一定期間内の全体的な状態を理解することを困難にする。
本開示は、瞬間的な電気インピーダンス変化画像が測定された人体の全体的な状態を反映することが困難であるという技術的問題に基づいて、電気インピーダンスイメージング方法、システム、記憶媒体および電子機器を提案する。
第1の態様において、本開示は、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得する段階と、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階と、複数の測定時点での瞬時微分画像に基づいて、画像行列を構築し、当該画像行列中の各列は、一つの瞬時微分画像に対応するベクトルである段階と、画像行列に基づいて、当該画像行列に対応する共分散行列を決定する段階と、共分散行列に従って、当該共分散行列の重みベクトルを獲得する段階と、共分散行列および重みベクトル基づいて、測定された人体領域の電気インピーダンス状態画像を獲得する段階とを含む、電気インピーダンスイメージング方法を提供する。
選択可能に、画像行列は、
Figure 2023545859000002

であり、ここで、Iは、画像行列であり、a(t)は、最初の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、2番目の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、N番目の測定時点での瞬時微分画像に対応するベクトルである。
選択可能に、前記画像行列に基づいて、画像行列に対応する共分散行列を決定する段階は、画像行列に基づいて、第1の予め設定計算式を用いて、画像行列に対応する共分散行列を計算する段階をさらに含み、第1の予め設定計算式は、
Figure 2023545859000003

であり、ここで、Cは、共分散行列であり、Iは、画像行列であり、Tは、行列の転置行列であり、Iは、M*Nサイズの時間平均行列であり、当該時間平均行列中の各列の要素は、
Figure 2023545859000004

であり、ここで、Mは、各フレームの瞬時微分画像中のピクセルの数であり、Nは、測定時点での数であり、a(t)は、i番目の測定時点に対応する瞬時微分画像のベクトルである。
選択可能に、前記共分散行列に従って、当該共分散行列の重みベクトルを獲得する段階は、共分散行列に従って、第2の予め設定計算式を用いて、当該共分散行列の重みベクトルを計算する段階をさらに含み、第2の予め設定計算式は、
Figure 2023545859000005

であり、ここで、wは、重みベクトルであり、wは、M*1サイズの列ベクトルであり、当該列ベクトル中の要素は、1または-1の値を取り、Tは、行列の転置行列であり、Cは、共分散行列である。
選択可能に、前記共分散行列および重みベクトル基づいて、測定された人体領域の電気インピーダンス状態画像を獲得する段階は、共分散行列および重みベクトル基づいて、第3の予め設定計算式を用いて、測定された人体領域の電気インピーダンス状態画像を計算する段階をさらに含み、第3の予め設定計算式は、
Figure 2023545859000006

であり、ここで、aは、電気インピーダンス状態画像であり、Cは、共分散行列であり、wは、重みベクトルである。
選択可能に、前記各測定時点での電気インピーダンス測定信号について、電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階は、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号から予め設定周波数範囲の信号を抽出する段階と、および各抽出された予め設定周波数範囲内の信号について、抽出された予め設定周波数範囲内の信号に基づいて、画像再構成アルゴリズムを使用して、対応する瞬時微分画像を構築する段階とをさらに含む。
選択可能に、画像再構成アルゴリズムは、線形最小二乗法を含む。
第2の態様において、本開示は、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得するように構成される取得モジュールと、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築するように構成される画像構築モジュールと、複数の測定時点での瞬時微分画像に基づいて、画像行列を構築し、当該画像行列中の各列は、一つの瞬時微分画像に対応するベクトルであるように構成される行列構築モジュールと、画像行列に基づいて、当該画像行列に対応する共分散行列を決定するように構成される共分散行列計算モジュールと、共分散行列に従って、当該共分散行列の重みベクトルを獲得するように構成される重みベクトル計算モジュールと、および共分散行列および重みベクトル基づいて、測定された人体領域の電気インピーダンス状態画像を獲得するように構成される状態画像構築モジュールとを含む、電気インピーダンスイメージングシステムを提供する。
第3の態様において、本開示は、記憶媒体を提供し、当該記憶媒体には、プログラムコードが格納され、当該プログラムコードがプロセッサによって実行される場合、上記の実施例のいずれか1項に記載の電気インピーダンスイメージング方法を実現する。
第4の態様において、本開示は、電子機器を提供し、当該電子機器は、メモリおよびプロセッサを含み、メモリには、プロセッサで動作できるプログラムコードが格納され、当該プログラムコードがプロセッサによって実行される場合、上記の実施例のいずれか1項に記載の電気インピーダンスイメージング方法を実現する。
本開示によって提供される電気インピーダンスイメージング方法、システム、記憶媒体および電子機器において、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を介して、測定された人体機能の一定期間内の全体的な状態を反映できる電気インピーダンス状態画像を再構築する。本開示によって提供される電気インピーダンスイメージング方法を使用して、測定された人体機能の一定期間内の全体的な状態を反映する電気インピーダンス状態画像を獲得することができ、従って、観察者が画像を理解し、且つ測定された全体的な人体機能を把握するのに便利であり、測定された人体機能の後続の定性的および定量的分析を容易にすることが分かる。
添付の図面と併せて以下の例示的な実施例の詳細な説明を読むことにより、本開示の範囲をより良く理解することができる。ここで、含まれる添付の図面は、次のとおりである。
本開示の実施例1で提案される電気インピーダンスイメージング方法のフロー概略図を示す。 人間の胸部の3次元血液灌流の電気インピーダンス状態画像の効果図を示す。
本開示の目的、技術的解決策および利点をより明確にするために、以下、添付の図面および実施例と併せて、本開示の実施方法を詳細に説明し、このようにして、本開示が技術的手段を適用して技術的問題を解決し、技術的効果を達成する実現プロセスを完全に理解し、且つそれに応じて実施することができる。
以下の説明では、本開示を完全に理解するために多くの具体的な詳細を説明するが、本開示は、ここに説明されるものとは異なる他の方法で実施することができ、従って、本開示の保護範囲は、以下に開示される具体的な実施例によって限定されない。
(実施例1)
本開示の実施例は、電気インピーダンスイメージング方法を提供し、図1は、実施例1で提案される電気インピーダンスイメージング方法のフロー概略図を示す。図1に示されるように、当該電気インピーダンスイメージング方法は、次のような段階S110~段階S160を含むことができる。
段階S110では、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得する。
ここで、電気インピーダンス測定では、まず測定された人体領域の周囲に電極アレイを固定する必要があり、電極アレイは、複数の電極を含み、次に電極アレイによって測定される人体領域を励起し、且つこれによって生成される応答を測定する。例えば、電極に電流励起を順番に印加し、他の電極でこれによって生成される電圧信号を測定することにより、電気インピーダンス測定信号を獲得する。
複数の測定時点とは、連続した期間内に、測定される人体領域に対して電気インピーダンス測定を複数回実行することにより、複数の測定時点での電気インピーダンス測定信号を獲得することを指す。
段階S120では、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する。
ここで、電気インピーダンス測定信号を獲得した後、電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築することができる。五つの測定時点がある場合、五つの測定時点に対応する瞬時微分画像をそれぞれ獲得することができる。
瞬時微分画像は、瞬時微分画像を再構築する測定時点で参照瞬間(例えば、呼気の終わりに対応する瞬間)に対する測定された人体領域の電気インピーダンス変化を反映することである。
選択可能な実施形態において、段階S120における各測定時点での電気インピーダンス測定信号について、電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階は、次のような段階S121~段階S122を含む。
段階S121では、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号から予め設定周波数範囲の信号を抽出する。
段階S122では、各抽出された予め設定周波数範囲内の信号について、抽出された予め設定周波数範囲内の信号に基づいて、画像再構成アルゴリズムを使用して、対応する瞬時微分画像を構築する。
ここで、段階S121において、信号の時間-周波数特性に従って、電気インピーダンス測定信号から予め設定周波数範囲内の信号を抽出し、当該予め設定周波数の範囲が換気周波数の範囲内の信号であり得る場合、抽出された信号は、換気関連信号であるか、または血液灌流周波数の範囲内の信号であり得る場合、抽出された信号は、血液灌流関連信号である。具体的には、フィルターを使用して、電気インピーダンス測定信号から予め設定周波数範囲内の信号を抽出することができる。
予め設定周波数範囲内の信号を抽出した後、抽出された予め設定周波数範囲内の信号を使用して、画像再構成アルゴリズムによって、瞬時微分画像を構築する。画像再構成アルゴリズムは、線形最小二乗法等の微分再構築アルゴリズムであり得る。
本実施形態では、線形最小二乗法を画像再構成アルゴリズムとして使用して、瞬時微分画像を再構築するが、当業者は、本開示では他の画像再構成アルゴリズムを使用することもできることを理解されたい。
胸部血液灌流の瞬時微分画像再構築を例にすると、その再構築プロセスは、次のとおりであり得る。電気インピーダンス測定信号から予め設定周波数範囲内の信号、即ち、血液灌流関連信号を抽出し、次に当該血液灌流関連信号に基づいて、画像再構成アルゴリズムを用いて画像再構成を実行して、血液灌流画像を獲得する。
段階S130では、複数の測定時点に対応する瞬時微分画像に基づいて、画像行列を構築し、当該画像行列中の各列は、一つの瞬時微分画像に対応するベクトルである。
ここで、複数の測定時点に対応する瞬時微分画像を獲得した後、複数の測定時点に対応する瞬時微分画像を使用して、画像行列を構築する。当該画像行列は、
Figure 2023545859000007

であり、ここで、Iは、画像行列であり、a(t)は、最初の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、2番目の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、N番目の測定時点での瞬時微分画像に対応するベクトルである。
段階S140では、画像行列に基づいて、当該画像行列に対応する共分散行列を決定する。
ここで、画像行列を構築した後、即ち、当該画像行列の共分散行列を計算することができる。
選択可能な実施形態において、前記画像行列に基づいて、画像行列に対応する共分散行列を決定する段階は、次のようなプロセスをさらに含む。
画像行列に基づいて、第1の予め設定計算式を用いて、画像行列に対応する共分散行列を計算し、当該第1の予め設定計算式は、
Figure 2023545859000008

であり、ここで、Cは、共分散行列であり、Iは、画像行列であり、Tは、行列の転置行列であり、Iは、M*Nサイズの時間平均行列であり、当該時間平均行列中の各列の要素は、
Figure 2023545859000009

であり、ここで、Mは、各フレームの瞬時微分画像中のピクセルの数であり、Nは、測定時点での数であり、a(t)は、i番目の測定時点に対応する瞬時微分画像のベクトルである。
ここで、各瞬時微分画像は、列ベクトルa(t)として表されることができ、ここで、tは、i番目の測定時点であり、i=1,2,…,N、Nは、測定時点での数である。ベクトルa(t)中の各要素は、画像中の一つのピクセル値を表す。当該画像行列は、
Figure 2023545859000010

であり、時間平均行列は、M*Nサイズの行列であり、当該行列中の各列の要素は、aである。第1の予め設定計算式によって、画像行列の共分散行列を計算することができる。
段階S150では、共分散行列に従って、当該共分散行列の重みベクトルを獲得する。
選択可能な実施形態において、前記共分散行列に従って当該共分散行列の重みベクトルを獲得する段階は、次のようなプロセスをさらに含む。
共分散行列に従って、第2の予め設定計算式を用いて、当該共分散行列の重みベクトルを計算し、第2の予め設定計算式は、
Figure 2023545859000011

であり、ここで、wは、重みベクトルであり、wは、M*1サイズの列ベクトルであり、当該列ベクトル中の要素は、1または-1の値を取り、Tは、行列の転置行列であり、Cは、共分散行列である。
ここで、重みベクトルwの解は、実際には0-1二次計画問題を解くことである。ここで、Mは、画像のピクセルの総数である。
段階S160では、共分散行列および重みベクトルに基づいて、測定された人体領域の電気インピーダンス状態画像を獲得する。
選択可能な実施形態において、前記共分散行列および重みベクトルに基づいて、測定された人体領域の電気インピーダンス状態画像を獲得する段階は、次のようなプロセスをさらに含む。
共分散行列および重みベクトル基づいて、第3の予め設定計算式を用いて、測定された人体領域的電気インピーダンス状態画像を計算し、第3の予め設定計算式は、
Figure 2023545859000012

であり、ここで、aは、電気インピーダンス状態画像であり、Cは、共分散行列であり、wは、重みベクトルである。
ここで、電気インピーダンス状態画像は、測定された人体領域の一定期間内の全体的な状態を反映できるため、観察者が画像を理解し、且つ測定された全体的な人体機能を把握するのに便利であり、測定された人体機能の後続の定性的および定量的分析を容易にする。
図2は、人間の胸部の3次元血液灌流の電気インピーダンス状態画像の効果図を示す。図2に示されるように、当該電気インピーダンス状態画像は、一定期間内の人間の胸部の3次元血液灌流の電気インピーダンス状態画像を反映し、当該画像に基づいて、観察者が画像を理解し、且つ測定された人体血液灌流の全体的な状況を把握するのに便利であり、測定された人体機能の後続の定性的および定量的分析を容易にする。
(実施例2)
本開示の実施例は、取得モジュール、画像構築モジュール、行列構築モジュール、共分散行列計算モジュール、重みベクトル計算モジュールおよび状態画像構築モジュールを含む、電気インピーダンスイメージングシステムをされに提供する。
取得モジュールは、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得するように構成される。
画像構築モジュールは、各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築するように構成される。
行列構築モジュールは、複数の測定時点での瞬時微分画像に基づいて、画像行列を構築し、当該画像行列中の各列は、一つの瞬時微分画像に対応するベクトルであるように構成される。
共分散行列計算モジュールは、画像行列に基づいて、当該画像行列に対応する共分散行列を決定するように構成される。
重みベクトル計算モジュールは、共分散行列に従って、当該共分散行列の重みベクトルを獲得するように構成される。
状態画像構築モジュールは、共分散行列および重みベクトル基づいて、測定された人体領域の電気インピーダンス状態画像を獲得するように構成される。
(実施例3)
本開示の実施例は、記憶媒体をさらに提供し、当該記憶媒体には、プログラムコードが格納され、当該プログラムコードがプロセッサによって実行される場合、上記の実施例のいずれか1項に記載の電気インピーダンスイメージング方法を実現する。
(実施例4)
本開示の実施例は、電子機器をさらに提供し、当該電子機器は、メモリおよびプロセッサを含み、メモリには、プロセッサで動作できるプログラムコードが格納され、当該プログラムコードがプロセッサによって実行される場合、上記の実施例のいずれか1項に記載の電気インピーダンスイメージング方法を実現する。
以上、添付の図面と併せて本開示の技術的解決策を詳細に説明するが、関連技術では、瞬間的な電気インピーダンス変化画像が測定された人体の全体的な状態を反映することは困難である。そのため、本開示は、電気インピーダンスイメージング方法、システム、記憶媒体および電子機器を提供し、測定された人体領域の複数の測定時点での電気インピーダンス測定信号を介して、測定された人体機能の一定期間内の全体的な状態を反映できる電気インピーダンス状態画像を再構築する。本開示の実施例によって提供される電気インピーダンスイメージング方法を使用することによって、測定された人体機能の一定期間内の全体的な状態を反映する電気インピーダンス状態画像を獲得することができ、従って、観察者が画像を理解し、且つ測定された全体的な人体機能を把握するのに便利であり、測定された人体機能の後続の定性的および定量的分析を容易にすることが分かる。
本出願によって提供されるいくつかの実施例において、開示された装置および方法は、他の方法で実現されることができることを理解されたい。例えば、上記で説明される装置の実施例は、例示にすぎず、例えば、ユニットの分割は、論理的な機能の分割のみにすぎず、実際の実現の場合、別の分割方法が存在する可能性があり、例えば、複数のユニットまたはコンポーネントを別のシステムに結合もしくは統合することができるか、または一部の機能を省略するか、もしくは実行しないことができる。
分離されたコンポーネントとして説明されるユニットは、物理的に分離されているか、または分離されていない場合があり、ユニットとして表示されるコンポーネントは、物理的なユニットである場合もそうでない場合もあり、即ち、1か所に配置されることができるか、または複数のネットワークユニットに分布されることができる。実際の必要性に従って、そのうちの部分または全部を選択して、本開示の実施例の目的を実現することができる。
さらに、本開示の各実施例における各機能ユニットは、一つの処理ユニットに統合されてもよく、各ユニットが物理的に別々に存在してもよく、または二つもしくはそれ以上のユニットが一つのユニットに統合される。上記の統合されたユニットは、ハードウェアの形で、またはソフトウェア機能ユニットの形で実現されることができる。
統合されたユニットがソフトウェア機能ユニットの形で実現され、独立した製品として販売または使用される場合、一つのコンピュータ可読記憶媒体に格納されることができる。このような理解に基づいて、本開示の技術的解決策は、本質的にまたは先行技術への貢献の一部であるか、または当該技術的解決策の全部もしくは一部をソフトウェア製品の形態で具現化することができ、当該コンピュータソフトウェア製品が一つの記憶媒体に格納され、いくつかの命令が含まれて、電子機器(パソコン、サーバー、ネットワーク機器等であり得る)に本開示の各実施例における方法の段階のすべてもしくは一部を実行させる。前述の記憶媒体は、Uディスク、モバイルハードディスク、読み取り専用メモリ(ROM、Read-Only Memory)、ランダムアクセスメモリ(RAM、Random Access Memory)、磁気ディスクまたは光ディスク等のプログラムコードを格納できる様々な媒体を含む。
本開示で開示される実施形態は、上記のとおりであるが、前記内容は、本開示の理解を容易にするための例示であり、本開示を限定するものではない。本開示が属する技術分野の当業者は、本開示に開示される精神および範囲から逸脱することなく、実施の形態および詳細に任意に修正および変更を加えることができるが、本開示の保護範囲は、依然として添付の特許請求の範囲によって定義されなければならない。

Claims (13)

  1. 電気インピーダンスイメージング方法であって、
    測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得する段階と、
    各測定時点での電気インピーダンス測定信号について、前記電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階と、
    複数の測定時点での瞬時微分画像に基づいて、画像行列を構築し、ここで、前記画像行列中の各列は、一つの瞬時微分画像に対応するベクトルである段階と、
    前記画像行列に基づいて、前記画像行列に対応する共分散行列を決定する段階と、
    前記共分散行列に従って、当該共分散行列の重みベクトルを獲得する段階と、
    前記共分散行列および前記重みベクトルに基づいて、前記測定された人体領域の電気インピーダンス状態画像を獲得する段階とを含むことを特徴とする、方法。
  2. 前記測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得する段階は、
    測定された人体領域の周囲に電極アレイを固定し、前記電極アレイは、複数の電極を含む段階と、
    前記電極アレイを介して測定された人体領域を励起し、これによって生成される応答を測定する段階とをさらに含むことを特徴とする
    請求項1に記載の方法。
  3. 前記画像行列は、
    Figure 2023545859000013

    であり、
    ここで、Iは、画像行列であり、a(t)は、最初の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、2番目の測定時点での瞬時微分画像に対応するベクトルであり、a(t)は、N番目の測定時点での瞬時微分画像に対応するベクトルであることを特徴とする
    請求項1に記載の方法。
  4. 前記画像行列に基づいて、前記画像行列に対応する共分散行列を決定する段階は、
    前記画像行列に基づいて、第1の予め設定計算式を用いて、前記画像行列に対応する共分散行列を計算する段階とをさらに含み、ここで、前記第1の予め設定計算式は、
    Figure 2023545859000014

    であり、
    ここで、Cは、共分散行列であり、Iは、画像行列であり、Tは、行列の転置行列であり、Iは、M*Nサイズの時間平均行列であり、当該時間平均行列中の各列は、すべて
    Figure 2023545859000015

    であり、ここで、Mは、各フレームの瞬時微分画像中のピクセルの数であり、Nは、測定時点での数であり、a(t)は、i番目の測定時点に対応する瞬時微分画像のベクトルであることを特徴とする
    請求項1に記載の方法。
  5. 前記共分散行列に従って、当該共分散行列の重みベクトルを獲得する段階は、
    前記共分散行列に従って、第2の予め設定計算式を用いて、当該共分散行列の重みベクトルを計算する段階をさらに含み、ここで、前記第2の予め設定計算式は、
    Figure 2023545859000016

    であり、
    ここで、wは、重みベクトルであり、wは、M*1サイズの列ベクトルであり、当該列ベクトル中の要素は、1または-1の値を取り、Tは、行列の転置行列であり、Cは、共分散行列であることを特徴とする
    請求項1に記載の方法。
  6. 前記共分散行列および前記重みベクトルに基づいて、前記測定された人体領域の電気インピーダンス状態画像を獲得する段階は、
    前記共分散行列および前記重みベクトルに基づいて、第3の予め設定計算式を用いて、前記測定された人体領域の電気インピーダンス状態画像を計算する段階をさらに含み、ここで、前記第3の予め設定計算式は、
    Figure 2023545859000017

    であり、
    ここで、aは、電気インピーダンス状態画像であり、Cは、共分散行列であり、wは、重みベクトルであることを特徴とする
    請求項1に記載の方法。
  7. 前記各測定時点での電気インピーダンス測定信号について、前記電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築する段階は、
    各測定時点での電気インピーダンス測定信号について、当該電気インピーダンス測定信号から予め設定周波数範囲の信号を抽出する段階と、
    各抽出された予め設定周波数範囲内の信号について、抽出された予め設定周波数範囲内の信号に基づいて、画像再構成アルゴリズムを使用して、対応する瞬時微分画像を構築する段階とをさらに含むことを特徴とする
    請求項1に記載の方法。
  8. フィルターを使用して、電気インピーダンス測定信号から予め設定周波数範囲内の信号を抽出することを特徴とする
    請求項7に記載の方法。
  9. 抽出された信号は、換気関連信号または血液灌流関連信号であることを特徴とする
    請求項7に記載の方法。
  10. 前記画像再構成アルゴリズムは、線形最小二乗法を含むことを特徴とする
    請求項7に記載の方法。
  11. 電気インピーダンスイメージングシステムであって、
    測定された人体領域の複数の測定時点での電気インピーダンス測定信号を取得するように構成される取得モジュールと、
    各測定時点での電気インピーダンス測定信号について、前記電気インピーダンス測定信号に基づいて、対応する瞬時微分画像を構築するように構成される画像構築モジュールと、
    複数の測定時点での瞬時微分画像に基づいて、画像行列を構築し、ここで、前記画像行列中の各列は、一つの瞬時微分画像に対応するベクトルであるように構成される行列構築モジュールと、
    前記画像行列に基づいて、前記画像行列に対応する共分散行列を決定するように構成される共分散行列計算モジュールと、
    前記共分散行列に従って、当該共分散行列の重みベクトルを獲得するように構成される重みベクトル計算モジュールと、
    前記共分散行列および前記重みベクトルに基づいて、前記測定された人体領域の電気インピーダンス状態画像を獲得するように構成される状態画像構築モジュールとを含むことを特徴とする、電気インピーダンスイメージングシステム。
  12. 記憶媒体であって、
    前記記憶媒体には、プログラムコードが格納され、前記プログラムコードがプロセッサによって実行される場合、請求項1~10のいずれか1項に記載の電気インピーダンスイメージング方法を実現することを特徴とする、前記記憶媒体。
  13. 電子機器であって、
    前記電子機器は、メモリおよびプロセッサを含み、前記メモリには、前記プロセッサ上で動作できるプログラムコードが格納され、前記プログラムコードが前記プロセッサによって実行される場合、請求項1~10のいずれか1項に記載の電気インピーダンスイメージング方法を実現することを特徴とする、前記電子機器。
JP2023524710A 2020-10-23 2021-08-18 電気インピーダンスイメージング方法、システム、記憶媒体および電子機器 Active JP7440985B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN202011145265.2A CN114469046A (zh) 2020-10-23 2020-10-23 电阻抗成像方法、系统、存储介质以及电子设备
CN202011145265.2 2020-10-23
PCT/CN2021/113160 WO2022083258A1 (zh) 2020-10-23 2021-08-18 电阻抗成像方法、系统、存储介质以及电子设备

Publications (2)

Publication Number Publication Date
JP2023545859A true JP2023545859A (ja) 2023-10-31
JP7440985B2 JP7440985B2 (ja) 2024-02-29

Family

ID=81291597

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2023524710A Active JP7440985B2 (ja) 2020-10-23 2021-08-18 電気インピーダンスイメージング方法、システム、記憶媒体および電子機器

Country Status (5)

Country Link
US (1) US11925449B2 (ja)
EP (1) EP4265179A4 (ja)
JP (1) JP7440985B2 (ja)
CN (1) CN114469046A (ja)
WO (1) WO2022083258A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115067921B (zh) * 2022-06-14 2023-06-13 杭州永川科技有限公司 基于脑电阻抗数据的体外循环术神经功能损伤预测系统
CN116269304B (zh) * 2023-03-09 2024-01-12 南京航空航天大学 一种拼接层预处理方法以及三维电阻抗成像方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010529870A (ja) * 2007-06-07 2010-09-02 デ モントフォート ユニバーシティ 電導性対象物の構造を分析する方法
US20120150458A1 (en) * 2010-12-07 2012-06-14 New York University Apparatus, method and computer-accessible medium for noninvasive determination of electrical properties of tissues and materials
JP2016531618A (ja) * 2013-09-27 2016-10-13 ドレーゲルヴェルク アクチェンゲゼルシャフト ウント コンパニー コマンディートゲゼルシャフト アウフ アクチェンDraegerwerk AG & Co.KGaA 電気インピーダンストモグラフィ装置及び電気インピーダンストモグラフィ方法
CN109758149A (zh) * 2019-01-18 2019-05-17 中国人民解放军第四军医大学 一种结合时差成像和频差成像的电阻抗成像方法
CN111067521A (zh) * 2019-12-31 2020-04-28 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置
JP2021016625A (ja) * 2019-07-22 2021-02-15 国立大学法人千葉大学 生体内物質の可視化装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE345086T1 (de) * 2000-06-09 2006-12-15 Boehm Stephan Dr Verfahren und vorrichtung zur darstellung von information die durch elektrische impedanz- tomografie erhalten ist
US10874324B2 (en) * 2015-10-20 2020-12-29 General Electric Company Detection of physiological changes in impedance monitoring
US10357178B2 (en) * 2015-12-16 2019-07-23 General Electric Company System and method for enhanced electrical impedance tomography
CN106510761B (zh) * 2016-12-12 2019-06-21 重庆大学 一种信噪比后滤波与特征空间融合的最小方差超声成像方法
GB2559157A (en) * 2017-01-27 2018-08-01 Ucl Business Plc Apparatus, method and system for alignment of 3D datasets
KR101948850B1 (ko) * 2017-05-15 2019-02-15 세종대학교산학협력단 자가센싱 콘크리트에 기반하여 구축된 구조물의 손상 모니터링을 위한 전기 임피던스 장치 및 방법
CN109086810A (zh) * 2018-07-19 2018-12-25 天津大学 一种基于模糊聚类的电阻抗层析成像优化算法
CN109745046B (zh) * 2019-01-22 2021-07-02 北京航空航天大学 一种适用于运动状态下的电阻抗成像电极与系统
CN109781791A (zh) * 2019-02-22 2019-05-21 广西大学 基于空间谱估计算法的电阻抗成像方法
CN109864712B (zh) * 2019-04-02 2024-07-05 北京华睿博视医学影像技术有限公司 电阻抗成像设备和方法
CN110251130B (zh) * 2019-05-24 2020-08-25 中国科学技术大学 基于可移动变形组件的电阻抗成像方法、装置及系统
CN110910466B (zh) * 2019-11-22 2022-11-18 东南大学 一种新型多频差分电阻抗层析成像重建算法
CN111192337B (zh) * 2019-12-24 2022-11-15 东南大学 基于稀疏贝叶斯学习的序贯多帧电阻抗成像图像重建算法
CN111281385B (zh) * 2020-03-06 2023-01-31 中国人民解放军第四军医大学 一种基于组织空间分布特征和阻抗随频率变化特性的电阻抗成像方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010529870A (ja) * 2007-06-07 2010-09-02 デ モントフォート ユニバーシティ 電導性対象物の構造を分析する方法
US20120150458A1 (en) * 2010-12-07 2012-06-14 New York University Apparatus, method and computer-accessible medium for noninvasive determination of electrical properties of tissues and materials
JP2016531618A (ja) * 2013-09-27 2016-10-13 ドレーゲルヴェルク アクチェンゲゼルシャフト ウント コンパニー コマンディートゲゼルシャフト アウフ アクチェンDraegerwerk AG & Co.KGaA 電気インピーダンストモグラフィ装置及び電気インピーダンストモグラフィ方法
CN109758149A (zh) * 2019-01-18 2019-05-17 中国人民解放军第四军医大学 一种结合时差成像和频差成像的电阻抗成像方法
JP2021016625A (ja) * 2019-07-22 2021-02-15 国立大学法人千葉大学 生体内物質の可視化装置
CN111067521A (zh) * 2019-12-31 2020-04-28 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置

Also Published As

Publication number Publication date
CN114469046A (zh) 2022-05-13
JP7440985B2 (ja) 2024-02-29
EP4265179A4 (en) 2024-02-14
US20230329575A1 (en) 2023-10-19
EP4265179A1 (en) 2023-10-25
WO2022083258A1 (zh) 2022-04-28
US11925449B2 (en) 2024-03-12

Similar Documents

Publication Publication Date Title
Wang et al. Physiological-model-constrained noninvasive reconstruction of volumetric myocardial transmembrane potentials
JP7440985B2 (ja) 電気インピーダンスイメージング方法、システム、記憶媒体および電子機器
US11363988B2 (en) Systems and methods for accelerated MRI scan
JP5751746B2 (ja) 磁気共鳴灌流画像の生成方法及び磁気共鳴装置
KR101809213B1 (ko) T2-이완 시간의 다중 슬라이스 영상화를 위한 방법 및 장치
US11940578B2 (en) Super resolution in positron emission tomography imaging using ultrafast ultrasound imaging
Weber et al. Predicting tissue conductivity influences on body surface potentials—an efficient approach based on principal component analysis
JP4299781B2 (ja) ソース再構成のための信頼区間を表示する方法及びシステム
US10159424B2 (en) Method and apparatus for generating medical image data records
US10557905B2 (en) System and method for producing high-resolution magnetic relaxation parameter maps
CN112703414B (zh) 分解后的多维图像的存储、显示和分析
WO2024007877A1 (zh) 用于呈现呼吸阻抗与呼吸阻抗变化率关系的方法及装置
Rahimi et al. Sensitivity of noninvasive cardiac electrophysiological imaging to variations in personalized anatomical modeling
US20240057887A1 (en) Three-dimensional ventilation image generation method, and controller and apparatus
CN114748052A (zh) 呼吸电阻抗成像激励参数选择方法、装置、设备与介质
JPH0399630A (ja) 等価双極子測定装置
Adler et al. GREIT: towards a consensus EIT algorithm for lung images
Arshad et al. Motion-robust free-running cardiovascular MRI
KR101297412B1 (ko) 가상 fMRI 데이터 생성 장치 및 방법
JP5369309B2 (ja) 画像処理装置、画像処理プログラム及び磁気共鳴装置
US11497412B2 (en) Combined oxygen utilization, strain, and anatomic imaging with magnetic resonance imaging
JPH10286244A (ja) 生体活動電流源推定装置
Zhang et al. RV-coefficient and its significance test in mapping brain functional connectivity
US20190246915A1 (en) System and method for reconstructing a physiological signal of an artery/tissue/vein dynamic system of an organ in a surface space
Degenhard et al. Consecutive projections onto convex sets

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230421

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20230421

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231017

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231106

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240130

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240208

R150 Certificate of patent or registration of utility model

Ref document number: 7440985

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313114