JPWO2015002210A1 - Eit測定装置、eit測定方法及びプログラム - Google Patents
Eit測定装置、eit測定方法及びプログラム Download PDFInfo
- Publication number
- JPWO2015002210A1 JPWO2015002210A1 JP2015525246A JP2015525246A JPWO2015002210A1 JP WO2015002210 A1 JPWO2015002210 A1 JP WO2015002210A1 JP 2015525246 A JP2015525246 A JP 2015525246A JP 2015525246 A JP2015525246 A JP 2015525246A JP WO2015002210 A1 JPWO2015002210 A1 JP WO2015002210A1
- Authority
- JP
- Japan
- Prior art keywords
- measurement
- strain gauge
- contour
- point
- eit
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0536—Impedance imaging, e.g. by tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1077—Measuring of profiles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6823—Trunk, e.g., chest, back, abdomen, hip
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/683—Means for maintaining contact with the body
- A61B5/6831—Straps, bands or harnesses
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/7475—User input or interface means, e.g. keyboard, pointing device, joystick
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/02—Operational features
- A61B2560/0242—Operational features adapted to measure environmental factors, e.g. temperature, pollution
- A61B2560/0247—Operational features adapted to measure environmental factors, e.g. temperature, pollution for compensation or correction of the measured physiological value
- A61B2560/0252—Operational features adapted to measure environmental factors, e.g. temperature, pollution for compensation or correction of the measured physiological value using ambient temperature
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0261—Strain gauges
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/043—Arrangements of multiple sensors of the same type in a linear array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/14—Coupling media or elements to improve sensor contact with skin or tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
本願は、2013年7月2日に、日本に出願された特願2013−139164号に基づき優先権を主張し、その内容をここに援用する。
EITは、微弱電流を流すだけで断層画像を取得できるので、X線CT(Computed tomography)と比較して、被曝の問題がなく、小型化や長時間測定、リアルタイムの測定が容易であるという利点がある。
そうすると、EIT測定装置のオペレータ(診療を行う者)は、輪郭の形状や大きさが異なる種々の測定対象について、生成された断層画像に基づいて正確な診断を行うことができない。
以下、第1の実施形態に係るEIT測定装置を、図面を参照して説明する。
図1は、第1の実施形態に係るEIT測定装置の全体構成を示す図である。図1において、EIT測定装置1を示している。
測定用ベルト10は、巻きつける長さの調整が可能な構成となっており、胸部などの他、頭部や腕、脚等にも巻かれて使用されてもよい。測定用ベルト10は、後述するように、EIT測定を行うための電極パッド等が同一のフレキシブル基板に設けられており、これら電極パッド等を一体として扱うことが可能な構成となっている(例えば、特許文献、特願2010−205988参照)。
図2に示すように、測定用ベルト10は、帯状のフレキシブル基板14上において8個の電極パッド12A〜12Hが等間隔の距離Pで並べて配置(周期配列)される構成となっている。また、同じくフレキシブル基板14上において、当該8個の電極パッド12A〜12Hと並列して8個の歪みゲージ13A〜13Hが等間隔の距離Pで周期配列されている。
測定用ベルト10は、これら複数の電極パッド12A〜12H、及び、歪みゲージ13A〜13Hを一体として貼着しながら、測定対象者の測定対象部Xに巻きつけられて使用される。各電極パッド12A〜12Hは、測定用ベルト10が巻きつけられた際に、測定対象部Xの周囲を囲む体表面に接触される仕組みとなっている。
このような構成とすることで、EIT測定装置1のオペレータは、測定用ベルト10を測定対象部Xへ巻きつけて取り付ける作業を行うのみで測定を開始することができるので、電極パッドの取り付けの手間を大幅に削減し、EIT測定の作業効率を改善することができる。
なお、測定用ベルト10は、測定時における測定対象者の負担を軽減する等の目的で、フレキシブル基板14が、さらに非導電性のベルト状の布に覆われた態様を成していてもよい。
また、この場合において、測定用ベルト10は、上記布ベルトのうち電極パッド12A〜12Hと接する部分が導電性ゲルや導電性繊維電極で構成されるとともに、測定装置1は、電極パッド12A〜12Hと体表面との間に当該導電性ゲルまたは導電性繊維電極を介在させながらEIT測定を行ってもよい。このようにすることで、測定用ベルト10の構成を、布ベルト部分のみを取り外して洗濯あるいは廃棄できるような構成とすることができ、衛生面での使い勝手を一層向上させることができる。
また、測定用ベルト10は、歪みゲージ13A〜13Hそれぞれの裏面の同じ位置に、同じ特性のひずみゲージを有していてもよい。これら測定箇所毎の一対の歪みゲージを一つのブリッジ回路に接続する2アクティブゲージ法を用いることで、歪みゲージ13A〜13Hの温度補正自動化と高感度化、高精度化を実現することができる。
或いは、測定用ベルト10は、歪みゲージ13A〜13Hそれぞれの裏面の同じ位置に、温度センサを有していてもよい。そして測定装置1は、歪みゲージ13A〜13Hによって取得される曲率データを、当該温度センサによって取得される温度データに基づいて、曲率データの温度補正を行ってもよい。このようにすることで、歪みゲージ13A〜13Hの高感度化、高精度化を実現することができる。
本実施形態に係るEIT測定本体部20は、汎用のパーソナルコンピュータ(パソコン)及び一般的な周辺機器(パソコンモニタ等)で構成される。
また、操作入力部212は、例えばマウス、キーボード、タッチパネル等から構成され、オペレータによる各種操作の入力を受け付ける。画像表示部213は、液晶ディスプレイ等であって、EIT測定時において必要な情報や、取得された断層画像等を表示する。
外部インターフェイス214は、外部装置との通信を行うための通信インターフェイスであり、特に本実施形態においては、専用の通信ケーブルを介して測定用ベルト10と接続され、測定用ベルト10から種々の信号を取得する機能部である。
CPU200、RAM210、HDD211、操作入力部212、画像表示部213、外部インターフェイス214は、システムバス215を介して相互に電気的に接続されている。
EIT測定制御部201は、複数の電極パッド12A〜12Hへの通電、及び、当該電極パッド12A〜12H間に生じる電圧信号の取得をしながら、測定対象部Xの断層画像を取得する。
また、輪郭推定部202は、歪みゲージ13A〜13Hを介して取得される曲率データに基づいて、測定対象部Xの輪郭形状及びその大きさを推定する。
また、図5は、第1の実施形態に係るEIT測定制御部の機能を説明する第2の図である。
以下、図4、図5を参照しながら、EIT測定制御部201の機能について説明する。図4は、測定対象部Xの例として、例えば測定対象者の胸部に測定用ベルト10が巻かれ、体表面に電極パッド12A〜12Hが接触した状態を示している。図4に示すように、電極パッド12A〜12Hの間には電流源Iと電圧計Vが接続される構成となっており、また、EIT測定制御部201は、この電流源Iと電圧計Vを制御する機能を有している。
図5は、EIT測定制御部201が生成した断層画像の例を示している。図5は、EIT測定制御部201によって取得された測定対象者の胸部における断層画像であり、電気インピーダンスが高い領域ほど濃い色合いで示している。図5によれば、空気の存在によって電気インピーダンスが高く測定される肺野が、左右に存在している様子を知ることができる。
なお、EIT測定制御部201が断層画像を生成する手法としては、上述した逆投影法の他に有限要素法(Finite Element Method:FEM)を用いる手法、または、有限要素法と逆投影法を組み合わせる手法が考えられる。逆投影法を用いた場合、EIT測定制御部201は、ある状態を基準とした相対的な変化のみしか画像化することができないが、有限要素法を用いることで、断層面における絶対的な電気抵抗率[Ωm]に基づいた断層画像を形成することができる。
したがって、EIT測定制御部201が取得した断層画像は、断層画像における各座標位置と、測定対象における実際の断層面における各位置との絶対的な位置関係が特定されたものとはなっていないため、オペレータは、この断層画像に基づいて正確な診断を行うことができない。
そこで、本実施形態に係るEIT測定装置1は、まず、測定用ベルト10に設けられた歪みゲージ13A〜13Hと、輪郭推定部202とにより、測定対象部Xの断層面の輪郭の形状を推定する処理を行う。そして、上述したEIT測定制御部201は、輪郭推定部202によって推定された輪郭の形状に基づいた断層画像を生成する。
また、図7〜図15は、各処理の具体的な内容を説明する第1〜第9の図である。
次に、本実施形態に係る輪郭推定部202が、測定対象部Xの断層面の輪郭形状を推定する処理を具体的に説明する。
図6に示すように、輪郭推定部202は、相対位置特定処理S1、形状特定処理S2、サイズ特定処理S3を経て、測定対象部Xの輪郭形状及び大きさを推定する。以下、各処理S1〜S3の内容について、図6及び図7〜図15を参照しながら、詳細に説明する。
まず、輪郭推定部202は、相対位置特定処理S1において、歪みゲージ13A〜13Hごとの相対的な位置関係を特定する。
具体的には、図7に示すように、輪郭推定部202は、第1ステップS10(図6)において、まず所定のxy座標を設定するとともに、測定用ベルト10の長手方向に沿って並べて配置される歪みゲージ13A〜13Hのうち、一つおきに指定される歪みゲージ(歪みゲージ13A、13C、13E、13G)それぞれについてのxy座標上の座標位置を仮想的に示す基準点N13A、N13C、N13E、N13Gを、所定のxy座標上に仮決めする処理を行う。
なお、相対的な位置関係によって関連付けられた1つの基準点(例えばN13A)と、当該基準点の両隣の2つの従属点(例えばN13B1、N13H2)との組み合わせを一組の「束」と定義する。
まず、HDD211は、既知のデータとして、歪みゲージ13Aと歪みゲージ13Bとの間隔が距離Pであることを記憶している。そして輪郭推定部202は、HDD211を参照して、歪みゲージ13Aと歪みゲージ13Bとの間隔を示す距離Pを読み取る。また、輪郭推定部202は、歪みゲージ13Aと歪みゲージ13Bによって検出された曲率データから、各々が配される箇所における曲率半径を取得する。
ここで、歪みゲージ13Aが配される箇所における曲率半径がRa、歪みゲージ13Bが配される箇所における曲率半径がRbであったとする。
なお、輪郭推定部202は、微小距離dPを十分に小さく設定した場合、点Q1と点Q2の間の微小距離dPを円弧の長さとしてではなく、点Q1と点Q2の直線距離と近似して点Q2の座標位置を特定する処理を行ってもよい。
次に、輪郭推定部202は、点Q2から微小距離dPだけ従属点N13B1に近づいた点Q3の座標位置を特定する。この点Q3の座標位置の特定において輪郭推定部202は、点Q2における曲率半径がR2であることから、点Q2と中心O1とを結ぶ直線上において、点Q2を基準とした曲率半径R2によって定まる中心O2を特定し、その中心O2,半径R2で定まる円弧上の点であって、点Q2から従属点N13B1の方向に微小距離dPだけ離れた点Q3の座標位置を特定する(図8参照)。
以上の処理により、基準点N13Aと従属点N13B1との相対的な位置関係を示す、従属点N13B1についての相対座標値が特定される。同様に、輪郭推定部202は、従属点N13H2についての相対座標値を特定する。輪郭推定部202は、このようにして、相対的な位置関係が特定された基準点N13Aと、2つの従属点N13B1、N13H2と、を一組の束として関連付ける。
輪郭推定部202は、以上のような処理により、実際の体表面の形状を精度よく反映する仮定(歪みゲージ間の各点における曲率半径がなだらかに変化するものと仮定)のもと、基準点に対する従属点の相対的な位置関係を特定することができる。
図9に示すように、輪郭推定部202は、基準点N13Aの座標位置(xa、ya)を基準とした相対座標値(dxb、dyb)、(dxh、dyh)から、2つの従属点N13B1、N13H2の座標位置(xb、yb)、(xh、yh)を特定する。そして輪郭推定部202は、基準点N13Aの座標位置(xa、ya)と、2つの従属点N13B1、N13H2の座標位置(xb、yb)、(xh、yh)とを、束T1の情報として関連付けてRAM210に一時的に記録する。従属点N13B1の座標位置(xb、yb)は、この段階において、xb=xa+dxb、yb=xa+dybで算出できる。
このように、束T1においては、基準点N13Aと、基準点N13Aに関連付けられた2つの従属点N13B1、N13H2の相対的な位置関係が特定されている。
同様に、輪郭推定部202は束T2の情報として、基準点N13Cの座標位置の情報と、基準点N13Cに関連付けられた2つの従属点N13D1、N13B2の各座標位置の情報とを対応付けてRAM210に一時的に記録する。
さらに、輪郭推定部202は束T3、束T4の情報として、基準点N13E、N13Gの座標位置の情報と、基準点N13E及びN13Gに関連付けられた2つの従属点N13F1、N13D2及びN13H1、N13F2の各座標位置の情報とを対応付けてRAM210に一時的に記録する。
今、基準点の位置に対応する歪みゲージそれぞれの間に一つのみの電極パッドが存在し、体に貼着している電極パッドの数が合計8個であるとすると、異なる束(T1またはT2)に属する2つの従属点N13B1、N13B2は、同一の歪みゲージ13Bの位置を示す。したがって、この場合、2つの従属点N13B1(第1従属点)、N13B2(第2従属点)は、本来、同一の座標位置で表されるものと考えられる。
また、図9において、束T2に属する従属点N13D1と束T3に属する従属点N13D2は、同一の歪みゲージ13Dの位置を示している。同様に、従属点N13F1とN13F2は、同一の歪みゲージ13Fの位置を示すもの、そして、従属点N13H1とN13H2は、同一の歪みゲージ13Hの位置を示している。
したがって、これらの従属点N13D1とN13D2、N13F1とN13F2、及び、N13H1とN13H2についても同様の考えが成り立つ。
同様に、輪郭推定部202は、束T2をx軸に沿って原点Oに平行移動させるように座標位置を変更する処理を行う。
ここで、輪郭推定部202は、束T1、T2をそれぞれy軸及びx軸に沿って平行移動させた場合に、従属点N13B1とN13B2とが重なる一致点K(図10)の座標位置を特定する。ここで、輪郭推定部202は、従属点N13B1の現時点における座標位置が(xb1、yb1)であって、従属点N13B2の現時点における座標位置が(xb2、yb2)であった場合に、一致点Kの座標位置を(xb1、yb2)と特定する。
そして、輪郭推定部202は、各従属点N13B1と一致点Kとのy座標値についての差分(yb1−yb2)を算出し、束T1に含まれる基準点N13Aのy座標値について、差分(yb1−yb2)を減算する処理を行う。同様に、輪郭推定部202は、各従属点N13B2と一致点Kとのx座標値についての差分(xb1−xb2)を算出し、束T2に含まれる基準点N13Cのx座標値について、差分(xb1−xb2)を加算する処理を行う。
このようにして、輪郭推定部202は、従属点N13B1、N13B2の座標位置を一致させるように、各束T1、T2の座標位置を変更する。
図11に示すように、輪郭推定部202は、束T3をy軸に沿って平行移動させた際に、束T3の従属点N13D2が、束T2の従属点N13D1に最も接近するように、束T3各点の座標位置を変更する処理を行う。なお、この段階では束T2の座標位置は既に固定されているので、従属点N13D1とN13D2の座標位置は完全には一致せず、幾分の誤差を有するものとなる。
この処理の後、輪郭推定部202は、さらに、束T4をx軸に沿って平行移動させた際に、束T4の従属点N13F2が束T3の従属点N13F1に最も接近するように、束T4の各点の座標位置を変更する処理を行う。この段階では、上記と同様に、束T3の座標位置は既に固定されているので、従属点N13F1とN13F2の座標位置は完全に一致しない。また、束T4の従属点N13H1と束T1の従属点N13H2の座標位置も一致しない。
したがって、この段階では、各従属点の対、N13B1とN13B2、N13D1とN13D2、N13F1とN13F2、N13H1とN13H2、における座標位置の誤差量が不均一な状態となっている。よって、輪郭推定部202は、各従属点の対における座標位置の誤差量が均一化されるように、各束T1〜T4の座標位置を再度微調整する。例えば、輪郭推定部202は、各束T1〜T4のそれぞれについて、束T1〜T4に含まれる各点の相対的な位置関係を維持したまま微小移動を行い、当該微小移動ごとに誤差量を算出しつつ、その誤差量が均一化された時点で微小移動を停止する。
図12に示すように、各従属点の対、N13B1とN13B2、N13D1とN13D2、N13F1とN13F2、N13H1とN13H2における誤差量が均一化されたところで、各束T1〜T4の最終的な位置が特定される。
そして、輪郭推定部202は、従属点N13B1とN13B2各々が示す座標位置の中点を算出し、ここで算出された中点の座標位置を、歪みゲージ13Bの位置を示す点N13Bとして特定する(図6、第4ステップS13)。輪郭推定部202は、同様の処理を、他の従属点の対(N13D1とN13D2、N13F1とN13F2、N13H1とN13H2)に対しても行う。
図13に示すように、輪郭推定部202は、相対位置特定処理S1(第1ステップS10〜第3ステップS12)の処理を経て、各歪みゲージ13A〜13Hの相対的な位置関係を示す点N13A〜N13Hを、座標上において特定する。
なお図10〜13を用いて説明した処理は、より一般的に説明すると、各束T1〜T4が示す基準点(N13A〜N13G)の座標位置それぞれを原点O方向に移動させた場合において、当該基準点と従属点との相対的な位置関係が当該移動前の相対的な位置関係を保持した状態の各束における同一の歪みゲージを示す各従属点の座標位置の距離が最も近くなるような、各束を構成する基準点と従属点の移動後の座標位置を算出する。そして、輪郭推定部202はその同一の歪みゲージを示す従属点間の中点を、当該歪みゲージの位置を示す点として特定する。
本実施形態に係る輪郭推定部202は、相対位置特定処理S1を完了して、各歪みゲージ13A〜13Hの相対的な位置関係を特定した後、これらを所定の関数曲線で結び合わせながら、測定対象部Xの輪郭形状を特定する形状特定処理S2を実行する(図6)。
図14Aに示すように、輪郭推定部202は、歪みゲージ13Aの位置を示す点N13Aと、歪みゲージ13Aに隣接する歪みゲージ13Bの位置を示す点N13Bとを結ぶ曲線ABを特定する。ここで輪郭推定部202が特定する曲線ABは、歪みゲージ13Aと歪みゲージ13Bの間における輪郭の形状を推定するものとなる。
まず、輪郭推定部202は、点N13Aと原点Oとを結ぶ長さrAの線分と、点N13Bと原点Oとを結ぶ長さrBの線分と、が原点Oにおいて成す角度θABを特定する(図14A参照)。
次に、輪郭推定部202は、HDD211を参照して、予め設定された微小角度dθを特定する。この微小角度dθは、角度θABを均等にm等分(mは2以上の整数)して微小に分割した角度を示している。そして、輪郭推定部202は、点N13Aと原点Oとを結ぶ線分を微小角度dθだけ傾けた直線上において補完点H1を特定する。次いで、輪郭推定部202は、さらに微小角度dθだけ傾けた直線上において補完点H2を特定する。このように、輪郭推定部202は、微小角度dθごとに特定される原点Oを通る複数の直線の各々について、補完点H1、H2、・・・を特定する(図14A参照)(図6、ステップS20)。
また、輪郭推定部202は、相対位置特定処理S1を経て点N13A、点N13Bの座標位置を特定しているので、その座標位置から点N13Aと原点Oとの距離rA、点N13Bと原点Oとの距離rBを算出する。そして、輪郭推定部202は、f(0)=rA、f(θAB)=rBとなる拘束条件を満たすように2次関数f(θ)を設定する(図14B参照)。
輪郭推定部202は、上記処理を繰り返して、補完点H1、H2、・・・の全てについての位置を特定する(図6、ステップS21)。
本実施形態に係る輪郭推定部202は、形状特定処理S2を完了して、測定対象となる部分の輪郭の形状を推定した後、当該輪郭の形状の大きさが現実の測定対象となる部分の大きさと一致するように、そのサイズを特定するサイズ特定処理S3を実行する(図6)。
ここで別途測定された周囲長とは、EIT測定装置1を用いたEIT測定とは別に、例えばメジャー等を用いて実測された、測定対象部Xの周囲長である。
輪郭推定部202は、形状特定処理S2により形状が推定された輪郭の周囲長を算出し、これを、メジャー等を用いて別途測定された測定対象部Xの周囲長と比較する。そして、推定された輪郭の周囲長と、上記実測された周囲長とが異なる場合には、推定された輪郭形状を維持したまま、全体のサイズを拡大または縮小して、推定された輪郭の周囲長と、上記実測された周囲長とが一致するように調整する。
図16は、第1の実施形態に係る輪郭推定部の処理結果を示す図である。
図16は、輪郭推定部202が上記各処理S1〜S3を経て推定した測定対象者の胸郭の輪郭形状Aと、当該胸郭について別途X線CTによって取得されたCT画像と、を重ね合わせた図を示している。図16に示すように、輪郭推定部202が、測定用ベルト10に備えられた歪みゲージ13A〜13Hに基づいて推定した輪郭形状Aは、CT画像に示される胸郭の輪郭形状と概ね一致している。
ここで、EIT測定制御部201が、輪郭推定部202によって推定された輪郭形状Aに基づいて断層画像を生成する具体的な手法について簡単に説明する。ここで、EIT測定制御部201が上述した有限要素法に基づいて断層画像を生成する場合を考える。
まず、EIT測定制御部201は、予めHDD211に記憶された所定のFEMモデルβ(図17A左)を指定する。このFEMモデルβは、他の被測定者によるCT画像等に基づいて予め作成されたものであって、EIT測定における測定対象部Xと同一の部位(例えば胸部)が指定される。そして、EIT測定制御部201は、輪郭推定部202によって推定された輪郭形状Aと、指定したFEMモデルβの輪郭形状B(図17A右)と、を重ね合わせて比較する(図17B参照)。
次いで、EIT測定制御部201は、算出された伸縮比率(rx、ry)を、最初に指定したFEMモデルβに適用し、同比率でX方向、Y方向にそれぞれ伸縮されたFEMモデルβ’を生成する。
そして、EIT測定制御部201は、電極パッド12A〜12Hを介して取得した電圧信号に対し、FEMモデルβ’についての有限要素法を適用することで、測定対象部Xの輪郭形状及び大きさに即した断層画像を取得する。
したがって、オペレータは、その断層画像と、測定対象とした部分の絶対的な位置関係を特定することができるので、測定対象ごとに形状や大きさが異なっていたとしても、当該断層画像に基づいて正確な診断を行うことができる。
例えば、輪郭推定部202が平行移動させる束T1〜T4を動かす順番は任意であり、図9に示す状態から図12に示す状態へと移行することができれば、その間の処理内容は特に限定されない。
またこの場合において、EIT測定本体部20の機能が実装された測定回路11は、さらに外部機器との無線または有線通信機能を有しているとともに、取得された断層画像を、種々の端末装置(スマートフォン、タブレット型コンピュータ、小型ゲーム機等)に無線あるいは有線で送信する機能を備えていてもよい。
次に、第2の実施形態に係るEIT測定装置1について、図面を参照しながら説明する。なお、第1の実施形態と同一の機能構成については、同一の符号を付してその説明を省略する。
なお、図18において、図面が煩雑となるのを避けるため、電極パッド12A〜12H及び歪みゲージ13A〜13Hの表記を省略しているが、実際には、第1の実施形態と同様に、フレキシブル基板14上に周期配列された構成となっている。
図18に示すように、本実施形態に係る測定用ベルト10は、その長手方向に沿って(電極パッド12A〜12H等と並列して)、間隔αで周期配列されながら測定用ベルト10に貼着された首位長計測用電極パッド301、302、・・・、30fを備えている。
図19に示すように、本実施形態に係るEIT測定本体部20は、さらに、周囲長測定部203を備えている。周囲長測定部203は、周囲長計測用電極パッド301〜30fを介して取得される電圧信号に基づいて、測定対象部Xの周囲長を計測する機能を有している。周囲長測定部203は、第1の実施形態におけるEIT測定制御部201、輪郭推定部202と同様に、EIT測定本体部20全体の動作を司るCPU200が、測定用プログラム実行時においてその機能を発揮する態様であってもよい。
例えば、図19に示すような状態のとき、周囲長計測用電極パッド301と周囲長計測用電極パッドとが密接しているため、当該周囲長計測用電極パッド301、30b間の電気インピーダンスは、両電極の容量結合に基づいて低下する。他の電極パッド対である302と30cの間、303と30eの間、304と30fの間についても同様の結果となる。一方、周囲長計測用電極パッド30a、305は、他の電極パッドとの容量結合が成されないため、周囲長測定部203は、周囲長計測用電極パッド30a、305については高い電気インピーダンス値を検出する。
なお周囲長計測用電極パッドの間隔αは、必ずしも電極パッドの間隔と同じ必要はなく、間隔αを小さくすることで周囲長測定精度を向上させることができる。
よって、第2の実施形態に係るEIT測定装置によれば、EIT測定を行う手順を一層簡素化することができる。
また、上述の説明において、各実施形態に係る測定用ベルト10は、電極パッド12A〜12H及び歪みゲージ13A〜13Hを長手方向に8個ずつ周期配列した構成であるとし、また、輪郭推定部202における各処理S1〜S3の内容は、歪みゲージ13A〜13Hの配列数が8個であることを前提として説明した。
しかし、他の実施形態に係るEIT測定装置1は、例えば、電極パッド及び歪みゲージを16個または32個備える態様としてもよい。例えば測定用ベルト10が伸縮可能な機構を備えていない場合には、測定対象部Xの大きさによって予め測定用ベルト10に複数設けられた電極パッド12A〜12Hや歪みゲージ13A〜13Hのうち測定対象部Xの体表面に接触する電極パッド12A〜12Hや、これに並列して配置される歪みゲージ13A〜13Hの数が変化することが考えられる。
このような態様の場合、例えば、輪郭推定部202は、測定用ベルト10上の周期配列において一つ以上の歪みゲージを隔てた所定間隔ごとに指定される歪みゲージの位置を示す、予め定められた基準点の座標位置を、HDD211から読み取り、その読み取った座標位置を初期座標値と特定する(第1ステップS10)。
そして、当該基準点の座標位置に対する、当該基準点が示す歪みゲージの中間に配される何れかの歪みゲージの位置を示す従属点の座標位置を特定する相対座標値を、歪みゲージを介して取得した曲率データに基づいて算出する(第2ステップS11)。
以下、第1の実施形態で説明した第3ステップS12及び第4ステップS13を行うことで、各歪みゲージの相対的な位置関係を特定することができる。
上述したように、第1の実施形態に係る相対位置特定処理S1では、一つおきに指定される歪みゲージ(歪みゲージ13A、13C、13E、13G(図2))のxy座標上における仮想位置を仮決めする際、EIT測定装置1は、基準点N13A、N13E、及び、基準点N13C、N13Gを、それぞれ、互いに直交するx軸上及びy軸上に配する(図7参照)。すなわち、第1の実施形態に係るEIT装置1は、測定用ベルト10が測定対象部Xに巻かれる際、当該測定用ベルト10上に周期配列される歪みゲージのうちの4つ(歪みゲージ13A、13C、13E、13G)は、常に、互いに直交する軸(x軸、y軸)上に配されることを前提として計算を行っている。
図20に示すように、測定用ベルト10は、帯状のフレキシブル基板14上において、16個の歪みゲージ13A〜13Pが等間隔の距離Pで周期配列されている。なお、図20においては、図面が煩雑となるのを避けるため、測定回路11及び電極パッド12A〜12Hの表記を省略しているが、実際には、第1の実施形態と同様に設けられている。
図21に示すように、測定用ベルト10は、当該測定用ベルト10の端部に配される歪みゲージ13Aを測定対象物Xの胸部の中心(胸骨体)に合わせながら巻き付けられている。この場合、各歪みゲージ13A〜13Pごとの間隔(距離P)の整数倍(本実施形態では16倍)が胸囲L1に一致しないため、測定対象物X(測定対象者)の前後方向の対称軸にx軸を、左右方向の対称軸にy軸を定義した場合に、歪みゲージ13A以外の歪みゲージがx軸上、y軸上からずれた位置に配される。
したがって、第1の実施形態に係るEIT測定装置1のように、歪みゲージのうちの4つが、常に、互いに直交するx軸及びy軸上に配されることを前提として輪郭(胸郭)の推定処理を行うと誤差が生じてしまう。
そこで、本実施形態に係る輪郭推定部202は、このような場合であっても、高精度で測定対象物Xの輪郭を推定可能とすべく以下の処理を実行する。
本実施形態に係るEIT装置1を扱うオペレータは、まず、測定用ベルト10の一端側に配された歪みゲージ(本実施形態では、歪みゲージ13A)を第1基点αとして、測定対象物X(測定対象者)の胸部の中心(胸骨体上)に配する。
次に、オペレータは、測定用ベルト10における第1基点αを測定対象物Xの胸骨体上に固定したまま、当該測定用ベルト10を測定対象物Xに巻き付ける。
なお、この場合、測定用ベルト10の他端側の余りの部分に配される歪みゲージ13Pは、輪郭の形状の推定処理には用いられない。
次に、本実施形態に係る輪郭推定部202が、測定対象部Xの断層面の輪郭形状を推定する処理を具体的に説明する。
第1の実施形態と同様に、輪郭推定部202は、相対位置特定処理S1において、歪みゲージ13A〜13Oごとの相対的な位置関係を特定する。
具体的には、図23に示すように、輪郭推定部202は、第1ステップS10(図6)において、まず所定のxy座標を設定するとともに、第1基点αに配される歪みゲージ13Aがy軸上に配されるものとして、当該歪みゲージ13Aについてのxy座標上の座標位置を示す基準点N13Aを、xy座標上に仮決めする。
なお、本実施形態に係る輪郭推定部202は、上記相対的な位置関係によって関連付けられた1つの基準点(例えば、基準点N13A)と、当該基準点の両側に、直列に連続して配される従属点(例えば、N13B1、N13C1、N13D1・・・等)との組み合わせを一組の「束」と定義する。
また、仮想の歪みゲージ13Qの両側に隣接して配される歪みゲージは、歪みゲージ13Hと、歪みゲージ13Iである(図21、図22参照)。この場合において、輪郭推定部202は、基準点N13Qの位置(xq、yq)を基準とした歪みゲージ13Hの相対的な位置を示す従属点を従属点N13H2、基準点N13Qの位置を基準とした歪みゲージ13Iの相対的な位置を示す従属点を従属点N13I2として、各従属点N13H2、N13I2のxy座標上の位置を特定する。
ただし、輪郭推定部202は、基準点N13Aに対する従属点N13O1の相対的位置を特定する場合においては、基準点N13Aから従属点N13O1までの距離が、オペレータによって入力された距離P1B(<距離P)であることを利用する。具体的には、輪郭推定部202は、図8において、基準点N13Aからの微小距離dPの積算量が距離P1Bとなった地点を従属点N13O1とする。なお、輪郭推定部202は、基準点N13Aに対する従属点N13B1の相対的位置を特定する場合は、基準点N13Aから基準点N13B1までの距離が距離Pであることを利用する。
更に、輪郭推定部202は、基準点N13Qから従属点N13I2までの距離が距離(P−P1A)であることを利用して、従属点N13I2との相対的位置を特定する。具体的には、輪郭推定部202は、図8において、基準点N13Qからの微小距離dPの積算量が距離(P−P1A)となった地点を従属点N13I2とする。
なお、上記の場合において、本実施形態に係る輪郭推定部202は、第2基点βに配される仮想の歪みゲージ13Qにより取得される曲率がゼロ(曲率半径が無限大)であると仮定した上で、基準点N13Qに対する従属点N13H2、N13I2の相対的位置を特定する。
図24に示すように、輪郭推定部202は、基準点N13Aの位置(xa、ya)を基準とした従属点N13B1の相対的な位置(xb、yb)、及び、従属点N13O1の相対的な位置(xo、yo)を特定する。同様に、輪郭推定部202は、基準点N13Qの位置(xq、yq)を基準とした従属点N13H2の相対的な位置(xh、yh)、及び、従属点N13I2の相対的な位置(xi、yi)を特定する。
輪郭推定部202は、上記処理を繰り返して、基準点N13Aの一方側において直列に従属する従属点N13B1〜N13E1の各々の位置、及び、基準点N13Aの他方側において直列に従属する従属点N13L1〜N13O1の各々の位置を特定する。これにより、基準点N13Aを中心とした束T1の位置が仮決めされる(図25参照)。
輪郭推定部202は、上記処理を繰り返して、基準点N13Qの一方側において直列に従属する従属点N13D2〜N13H2の各々の位置、及び、基準点N13Qの他方側において直列に従属する従属点N13I2〜N13M2の各々の位置を特定する。これにより、基準点N13Qを中心とした束T2の位置が仮決めされる(図25参照)。
つまり、異なる束(T1またはT2)に属する2つの従属点N13D1、N13D2は、同一の歪みゲージ13Dの位置を示す。また、異なる束に属する2つの従属点N13E1、N13E2は、同一の歪みゲージ13Eの位置を示す。同様に、異なる束に属する2つの従属点N13L1、N13L2は、同一の歪みゲージ13Lの位置を示す。そして、異なる束に属する2つの従属点N13M1、N13M2は、同一の歪みゲージ13Mの位置を示す。
このとき、輪郭推定部202は、同一の歪みゲージ13Dの位置を示す従属点N13D1、N13D2、同一の歪みゲージ13Eの位置を示す従属点N13E1、N13E2、同一の歪みゲージ13Lの位置を示す従属点N13L1、N13L2、同一の歪みゲージ13Mの位置を示す従属点N13M1、N13M2の各々の座標位置の誤差の合計が最も小さくなる位置に移動する。
図26に示すように、束T1と束T2とがオーバーラップする領域に属する各従属点の対であるN13D1とN13D2、N13E1とN13E2、N13L1とN13L2、N13M1とN13M2における位置の誤差の合計が最小となるところで、各束T1、T2の最終的な位置が特定される。
図27は、上述した第4ステップS13の処理が完了した直後の状態を示している。なお、図27において、従属点N13B〜N13Oは、それぞれ、歪みゲージ13B〜13Oに対応するxy座標上の位置を示している。ここで、従属点N13B、N13C、N13N、N13Oは、それぞれ、束T1に属していた従属点N13B1、N13C1、N13N1、N13O1の位置と同一の位置とされる。また、従属点N13F、N13G、N13H、N13I、N13J、N13Kは、それぞれ、束T2に属していた従属点N13F2、N13G2、N13H2、N13I2、N13J2、N13K2と同一の位置とされる。そして、束T1と束T2とがオーバーラップする領域に属する従属点N13D、N13E、N13L、N13Mは、それぞれ、束T1、T2の各々に属する従属点N13D1、N13D2、従属点N13E1、N13E2、従属点N13L1、N13L2、従属点N13M1、N13M2の中点の位置とされる。
輪郭推定部202は、以上の第1ステップS10〜第4ステップS13の処理を経て、15個の歪みゲージ13A〜13O(及び仮想の歪みゲージ13Q)のxy座標上における位置を特定する。
以上のように、第3の実施形態に係る輪郭推定部202は、相対位置特定処理S1の第1ステップS10において、複数の歪みゲージ13A〜13Pのうち、測定対象物Xの対称軸(胸骨体と背骨突起を結ぶ軸)上の所定位置(第1基点α)に配される歪みゲージ13Aの座標位置を予め定められた初期座標値と特定する。そして、この場合における初期座標値は、上記測定対象物Xの対称軸に対応する軸(y軸)上の所定位置とされる。
これにより、当該歪みゲージ13Aが測定対象物Xの対称軸に対応する軸上に配されることを拘束条件として相対位置特定処理S1を実行する際に、精度よく輪郭の推定を行うことができる。
このようにすることで、輪郭推定部202は、当該仮想の歪みゲージ13Qが測定対象物Xの対称軸に対応する軸上に配されることを拘束条件として相対位置特定処理S1を実行する。これにより、測定用ベルト10を測定対象物Xに巻き付けた場合に、測定対象物Xの対称軸上に歪みゲージが配されていない場合であっても、精度よく輪郭の推定を行うことができる。
図28は、測定対象物Xとして胸囲L2の胸部に測定用ベルト10を巻き付けた場合を示している。ここで、胸囲L2は、測定用ベルト10の一端に配される歪みゲージ13Aから他方の端部に配される歪みゲージ13Pまでの長さよりも長い長さである。
この場合も、測定用ベルト10は、当該測定用ベルト10の端部に配される歪みゲージ13Aを測定対象物Xの胸部の中心(胸骨体)に合わせながら巻き付けられている。また、この場合、各歪みゲージ13A〜13Pごとの間隔(距離P)の整数倍(本実施形態では16倍)が胸囲L2に一致しないため、測定対象物X(測定対象者)の前後方向の対称軸にx軸を、左右方向の対称軸にy軸を定義した場合に、歪みゲージ13A以外の歪みゲージがx軸上、y軸上からずれた位置に配される。
したがって、第1の実施形態に係るEIT測定装置1のように、歪みゲージのうちの4つが、常に、互いに直交するx軸及びy軸上に配されることを前提として輪郭(胸郭)の推定処理を行うと誤差が生じてしまう。
しかし、第3の実施形態の変形例に係るEIT測定装置1は、上記距離P1A、P1Bを自動で取得する態様としてもよい。
EIT測定装置1は、測定対象物Xに巻き付けられた測定用ベルト10が、当該測定対象物Xと密接する範囲を、当該周期配列される電極パッド間に生じる電気インピーダンスの変化に基づいて検出する。ここで、上記電極パッド間のインピーダンスは、当該電極パッドに対し、生体である測定対象物Xが密接しているか否かに応じて変化する。
これにより、EIT測定装置1は、測定対象物Xに密接する電極パッドと歪みゲージ13A〜13Pとの相対的な位置関係に基づいて、一端側の歪みゲージ(歪みゲージ13A)に隣接して配される他端側の歪みゲージ、及び、当該他端側の歪みゲージとの距離P1Bを自動的に特定することができる。
10 測定用ベルト
11 測定回路
12A〜12H 電極パッド
13A〜13P 歪みゲージ
14 フレキシブル基板
19 信号ケーブル
20 EIT測定本体部
200 CPU
201 EIT測定制御部
202 輪郭推定部
203 周囲長測定部
210 RAM
211 HDD
212 操作入力部
213 画像表示部
214 外部インターフェイス
301〜30f 周囲長計測用電極パッド
Claims (10)
- 並べて配置される複数の電極パッドと、当該複数の電極パッドと並列して配置される複数の歪みゲージと、を一体として貼着しながら、生体の測定対象とする部分に巻きつけられて使用される測定用ベルトと、
前記複数の電極パッドへの通電、及び、当該電極パッド間に生じる電圧信号の取得をしながら、前記測定対象とする部分の断層画像を取得するEIT測定制御部と、
前記歪みゲージを介して取得される曲率データに基づいて、前記測定対象とする部分の輪郭形状及びその大きさを推定する輪郭推定部と、
を備え、
前記輪郭推定部は、
前記曲率データに基づいて、前記歪みゲージごとの相対的な位置関係を特定する相対位置特定処理と、当該相対的な位置関係が特定された歪みゲージ間を、所定の関数曲線で結び合わせながら前記輪郭形状を特定する形状特定処理と、を実行するEIT測定装置。 - 前記輪郭推定部は、
前記形状特定処理の後に、さらに、前記推定された輪郭形状の周囲長が、別途測定された前記測定対象とする部分の周囲長と一致するように、前記推定された輪郭形状を拡大または縮小するサイズ特定処理を実行する
請求項1に記載のEIT測定装置。 - 前記輪郭推定部は、前記相対位置特定処理において、
並べて配置される複数の前記歪みゲージのうち、一つ以上の前記歪みゲージを隔てた所定間隔ごとに指定される歪みゲージの位置を示す基準点の座標位置を、予め定められた初期座標値と特定する第1ステップと、
前記基準点の座標位置に対する、当該基準点が示す歪みゲージの中間に配される何れかの歪みゲージの位置を示す従属点の座標位置を示す相対座標値を、前記歪みゲージを介して取得される曲率データに基づいて算出する第2ステップと、
一の前記基準点の座標位置に対する相対座標値から座標位置が特定される第1従属点と、他の基準点の座標位置に対する相対座標値から座標位置が特定される従属点であって前記第1従属点と同一の前記歪みゲージの位置を示す第2従属点と、の座標位置が最も近くなるように、前記一の基準点及び前記他の基準点の座標位置を変更する第3ステップと、
前記第3ステップによって変更した後の前記第1従属点と前記第2従属点との中点を、前記基準点が示す歪みゲージの両隣に配される歪みゲージの位置を示す2つの従属点の座標位置と特定する第4ステップと、
を実行する請求項1または請求項2に記載のEIT測定装置。 - 前記輪郭推定部は、
前記相対位置特定処理の前記第1ステップにおいて、前記複数の歪みゲージのうち、前記測定対象とする部分の対称軸上に配される歪みゲージを示す基準点の座標位置を、予め定められた初期座標値と特定する
請求項3に記載のEIT測定装置。 - 前記輪郭推定部は、
前記測定対象とする部分に巻き付けられた前記測定用ベルト上の、当該測定対象とする部分の対称軸上に配される位置に歪みゲージが配されない場合において、当該測定対象とする部分の対称軸上に配される位置に歪みゲージが配されているものとみなして、前記相対位置特定処理を実行する
請求項4に記載のEIT測定装置。 - 前記輪郭推定部は、前記形状特定処理において、
一の歪みゲージの位置と、当該一の歪みゲージに隣接する他の歪みゲージの位置と、を結ぶ曲線を特定する複数の補完点を設定するとともに、前記複数の補完点の原点からの距離は、当該補完点と、原点と、前記一の歪みゲージの位置と、が成す角度についての所定の関数によって定められる
請求項1から請求項5の何れか一項に記載のEIT測定装置。 - 前記複数の電極パッドと並列して配置されながら、前記測定用ベルトに貼着される周囲長計測用電極パッドと、
前記周囲長計測用電極パッドを介して取得される電圧信号に基づいて、前記測定対象とする部分の周囲長を計測する周囲長測定部と、
を備える請求項1から請求項6の何れか一項に記載のEIT測定装置。 - 並べて配置される複数の電極パッドと、当該複数の電極パッドに沿って配置される複数の歪みゲージと、が一体として貼着された測定用ベルトを、生体の測定対象とする部分に巻きつけ、
EIT測定制御部が、前記複数の電極パッドへの通電、及び、当該電極パッド間に生じる電圧信号を取得しながら、前記測定対象とする部分の断層画像を取得し、
輪郭推定部が、前記歪みゲージを介して取得される曲率データに基づいて、前記歪みゲージごとの相対的な位置関係を特定する相対位置特定処理と、当該相対的な位置関係が特定された歪みゲージ間を、所定の関数曲線で結び合わせながら前記測定対象とする部分の輪郭形状を特定する形状特定処理と、を実行し、前記輪郭形状及びその大きさを推定するEIT測定方法。 - 並べて配置される複数の電極パッドと、当該複数の電極パッドに沿って配置される複数の歪みゲージと、を一体として貼着するとともに、生体の測定対象とする部分に巻きつけられて使用される測定用ベルトと、を備えるEIT測定装置のコンピュータを、
前記複数の電極パッドへの通電、及び、当該電極パッド間に生じる電圧信号を取得しながら、前記測定対象とする部分の断層画像を取得するEIT測定制御手段、
前記歪みゲージを介して取得される曲率データに基づいて、前記測定対象とする部分の輪郭形状及びその大きさを推定する輪郭推定手段、
として機能させ、
前記輪郭推定手段は、
前記曲率データに基づいて、前記歪みゲージごとの相対的な位置関係を特定する相対位置特定処理と、当該相対的な位置関係が特定された歪みゲージ間を、所定の関数曲線で結び合わせながら前記輪郭形状を特定する形状特定処理と、を実行するプログラム。 - 並べて配置される複数の電極パッドと、当該複数の電極パッドと並列して配置される複数の曲率センサと、を一体として貼着しながら、生体の測定対象とする部分に巻きつけられて使用される測定用ベルトと、
前記複数の電極パッドへの通電、及び、当該電極パッド間に生じる電圧信号の取得をしながら、前記測定対象とする部分の断層画像を取得するEIT測定制御部と、
前記曲率センサを介して取得される曲率データに基づいて、前記測定対象とする部分の輪郭形状及びその大きさを推定する輪郭推定部と、
を備え、
前記輪郭推定部は、
前記曲率データに基づいて、前記曲率センサごとの相対的な位置関係を特定する相対位置特定処理と、当該相対的な位置関係が特定された曲率センサ間を、所定の関数曲線で結び合わせながら前記輪郭形状を特定する形状特定処理と、を実行するEIT測定装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013139164 | 2013-07-02 | ||
JP2013139164 | 2013-07-02 | ||
PCT/JP2014/067590 WO2015002210A1 (ja) | 2013-07-02 | 2014-07-01 | Eit測定装置、eit測定方法及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2015002210A1 true JPWO2015002210A1 (ja) | 2017-02-23 |
JP6385929B2 JP6385929B2 (ja) | 2018-09-05 |
Family
ID=52143782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015525246A Active JP6385929B2 (ja) | 2013-07-02 | 2014-07-01 | Eit測定装置、eit測定方法及びプログラム |
Country Status (4)
Country | Link |
---|---|
US (1) | US10285618B2 (ja) |
EP (2) | EP4154809B1 (ja) |
JP (1) | JP6385929B2 (ja) |
WO (1) | WO2015002210A1 (ja) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11457832B2 (en) | 2013-07-02 | 2022-10-04 | School Juridical Person Kitasato Institute | Measurement device, shape estimation device, measurement method, shape estimation method, and non-transitory recording medium recording program |
EP3054262B1 (en) * | 2013-09-30 | 2019-01-16 | School Juridical Person Kitasato Institute | Length measurement device, length measurement method, and program |
WO2017124378A1 (zh) * | 2016-01-21 | 2017-07-27 | 深圳迈瑞生物医疗电子股份有限公司 | 生理传感装置及包括该生理传感装置的生理监护设备 |
JP6754110B2 (ja) * | 2016-02-05 | 2020-09-09 | 学校法人北里研究所 | 運転手監視装置 |
EP3245942A1 (en) * | 2016-05-19 | 2017-11-22 | Nokia Technologies Oy | Location device for identifying location(s) for physiological measurements and methods using the same |
US10368781B2 (en) * | 2016-06-16 | 2019-08-06 | Maureen Shifflett | Device for measuring the circumference of a child's head |
EP3472571B1 (en) * | 2016-06-20 | 2022-03-02 | Koninklijke Philips N.V. | Medical coupling unit and sensor-side connector |
US11793418B2 (en) * | 2016-11-11 | 2023-10-24 | Sentec Ag | Sensor belt and positioning aid for electro-impedance tomography imaging in neonates |
DE102016014252B4 (de) | 2016-11-30 | 2023-02-02 | Drägerwerk AG & Co. KGaA | Vorrichtung und Verfahren zur Bestimmung einer Umfangsform einer Elektrodenanordnung zur Elektro-Impedanz-Tomographie |
EP3549515A4 (en) * | 2016-12-02 | 2020-08-05 | Bilab Co., Ltd. | ELECTRODE BELT DEVICE FOR MEASURING A BIOS SIGNAL |
US11412946B2 (en) * | 2017-11-14 | 2022-08-16 | Timpel Medical B.V. | Electrical impedance tomography device and system having a multi-dimensional electrode arrangement |
DE102018108643A1 (de) * | 2018-04-11 | 2019-11-14 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Positionsbestimmungsvorrichtung zum Bestimmen einer Position eines Gegenstands innerhalb einer tubulären Struktur |
WO2020017013A1 (ja) * | 2018-07-19 | 2020-01-23 | Posh Wellness Laboratory株式会社 | 検出装置、シートベルト、及び運転手監視システム |
GB2580164A (en) * | 2018-12-21 | 2020-07-15 | Imperial College Sci Tech & Medicine | A sensor |
JP7417251B2 (ja) | 2019-12-23 | 2024-01-18 | 国立大学法人千葉大学 | 生体内物質の可視化装置 |
JP7186458B2 (ja) * | 2020-08-06 | 2022-12-09 | 学校法人北里研究所 | 監視装置、監視方法、プログラム |
WO2022137566A1 (ja) * | 2020-12-25 | 2022-06-30 | Posh Wellness Laboratory株式会社 | ベルト型センサ、検出装置、監視システム、および推定方法 |
CN113080925B (zh) * | 2021-03-03 | 2024-03-01 | 桂林电子科技大学 | 一种电阻抗断层成像的三维边界测量结构及方法 |
CN113069100B (zh) * | 2021-03-03 | 2024-04-02 | 桂林电子科技大学 | 一种可旋转电阻抗断层成像的三维边界测量结构及方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003339658A (ja) * | 2002-05-30 | 2003-12-02 | Kao Corp | 体脂肪測定装置 |
US20050059903A1 (en) * | 2003-09-11 | 2005-03-17 | Tanita Corporation | Impedance type thickness measurement device |
JP2005152061A (ja) * | 2003-11-21 | 2005-06-16 | Matsushita Electric Works Ltd | 体脂肪測定装置 |
WO2007070997A1 (en) * | 2005-12-20 | 2007-06-28 | Dixtal Biomedica Industria E Comercio Ltda. | Electrode assembly for electrical impedance tomography |
WO2009118701A1 (en) * | 2008-03-27 | 2009-10-01 | Koninklijke Philips Electronics N.V. | Method and system for measuring an object of interest |
US20100198101A1 (en) * | 2007-09-24 | 2010-08-05 | Xubo Song | Non-invasive location and tracking of tumors and other tissues for radiation therapy |
WO2012045188A1 (en) * | 2010-10-07 | 2012-04-12 | Swisstom Ag | Sensor device for electrical impedance tomography imaging, electrical impedance tomography imaging intrument and electrical impeance tomography method |
JP2012090880A (ja) * | 2010-10-28 | 2012-05-17 | Kitasato Institute | 生体の電気的インピーダンス断層像測定装置 |
JP2012228514A (ja) * | 2011-04-23 | 2012-11-22 | Draeger Medical Gmbh | 電気インピーダンストモグラフィ装置および電気インピーダンストモグラフィのための方法 |
WO2012168836A2 (en) * | 2011-06-10 | 2012-12-13 | Koninklijke Philips Electronics N.V. | Dynamic constraining with optical shape sensing |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6519862B1 (en) * | 2000-10-04 | 2003-02-18 | The United States Of America As Represented By The Secretary Of The Navy | Device for the acquisition of contoured human body surface vibration signals |
KR100598146B1 (ko) | 2004-04-28 | 2006-07-07 | 메디게이트(주) | 체임피던스 측정용 벨트 전극 장치 |
JP4710881B2 (ja) * | 2007-06-29 | 2011-06-29 | パナソニック電工株式会社 | 体脂肪測定装置 |
JP2010205988A (ja) | 2009-03-04 | 2010-09-16 | Panasonic Corp | 窒化物半導体素子及びその製造方法 |
CN101524273B (zh) | 2009-04-03 | 2010-07-28 | 姜枫 | 环状电极终端装置及人体电阻率或体液含量的测量方法 |
JP5469571B2 (ja) | 2010-09-14 | 2014-04-16 | 学校法人北里研究所 | 生体の電気的インピーダンス断層像測定装置 |
JP2013139164A (ja) | 2011-12-28 | 2013-07-18 | Sumitomo Rubber Ind Ltd | 空気入りタイヤ |
-
2014
- 2014-07-01 JP JP2015525246A patent/JP6385929B2/ja active Active
- 2014-07-01 EP EP22206215.0A patent/EP4154809B1/en active Active
- 2014-07-01 WO PCT/JP2014/067590 patent/WO2015002210A1/ja active Application Filing
- 2014-07-01 EP EP14820626.1A patent/EP3017758B1/en active Active
- 2014-07-01 US US14/901,445 patent/US10285618B2/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003339658A (ja) * | 2002-05-30 | 2003-12-02 | Kao Corp | 体脂肪測定装置 |
US20050059903A1 (en) * | 2003-09-11 | 2005-03-17 | Tanita Corporation | Impedance type thickness measurement device |
JP2005152061A (ja) * | 2003-11-21 | 2005-06-16 | Matsushita Electric Works Ltd | 体脂肪測定装置 |
WO2007070997A1 (en) * | 2005-12-20 | 2007-06-28 | Dixtal Biomedica Industria E Comercio Ltda. | Electrode assembly for electrical impedance tomography |
US20100198101A1 (en) * | 2007-09-24 | 2010-08-05 | Xubo Song | Non-invasive location and tracking of tumors and other tissues for radiation therapy |
WO2009118701A1 (en) * | 2008-03-27 | 2009-10-01 | Koninklijke Philips Electronics N.V. | Method and system for measuring an object of interest |
WO2012045188A1 (en) * | 2010-10-07 | 2012-04-12 | Swisstom Ag | Sensor device for electrical impedance tomography imaging, electrical impedance tomography imaging intrument and electrical impeance tomography method |
JP2012090880A (ja) * | 2010-10-28 | 2012-05-17 | Kitasato Institute | 生体の電気的インピーダンス断層像測定装置 |
JP2012228514A (ja) * | 2011-04-23 | 2012-11-22 | Draeger Medical Gmbh | 電気インピーダンストモグラフィ装置および電気インピーダンストモグラフィのための方法 |
WO2012168836A2 (en) * | 2011-06-10 | 2012-12-13 | Koninklijke Philips Electronics N.V. | Dynamic constraining with optical shape sensing |
Non-Patent Citations (1)
Title |
---|
J.M. KHOR ET AL.: "Development of a Sensor Network for Dynamic Boundary Measurement in Neonatal Electrical Impedance To", IFMBE PROCEEDINGS, vol. 第25/II巻, JPN7017003687, 2009, pages 386 - 389, ISSN: 0003782646 * |
Also Published As
Publication number | Publication date |
---|---|
JP6385929B2 (ja) | 2018-09-05 |
EP4154809A1 (en) | 2023-03-29 |
EP3017758B1 (en) | 2022-12-14 |
US20160302690A1 (en) | 2016-10-20 |
WO2015002210A1 (ja) | 2015-01-08 |
US10285618B2 (en) | 2019-05-14 |
EP3017758A1 (en) | 2016-05-11 |
EP3017758A4 (en) | 2017-08-16 |
EP4154809B1 (en) | 2024-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6385929B2 (ja) | Eit測定装置、eit測定方法及びプログラム | |
JP4905197B2 (ja) | 内臓脂肪測定装置 | |
JP6583829B2 (ja) | 画像生成装置 | |
US20070252586A1 (en) | Instrument and Method for Measuring Three-Dimensional Motion | |
US20110262015A1 (en) | Image processing apparatus, image processing method, and storage medium | |
US9345415B2 (en) | Method and apparatus for measuring multiple ECG leads using a device with embedded leads | |
EP3035846B1 (en) | Apparatus and method for estimating shape | |
CN102131459B (zh) | 内脏脂肪测定装置 | |
JP2014064922A (ja) | 電極電位マッピング | |
US11457832B2 (en) | Measurement device, shape estimation device, measurement method, shape estimation method, and non-transitory recording medium recording program | |
CN116802695A (zh) | 生成解剖腔体的解剖模型 | |
CN116824048B (zh) | 一种传感器、雅可比矩阵求解方法、三维成像系统及方法 | |
JP2016093472A (ja) | 呼吸関連信号の計測システム及びその計測方法 | |
JP6602319B2 (ja) | 「高解像度」の電気インピーダンスイメージングのための装置および方法 | |
WO2023135625A1 (ja) | 評価装置、評価方法、および記録媒体 | |
JP4600735B2 (ja) | 画像重ね合せ装置 | |
JP6745495B2 (ja) | 形状推定装置、スキャニング装置、動作検出装置、形状推定方法、スキャニング方法、動作検出方法、プログラム | |
JP6518084B2 (ja) | 生体内信号源位置検出装置及び生体内信号源位置検出方法 | |
JP2016116629A (ja) | 体脂肪測定装置、体脂肪測定方法及びプログラム | |
JP2014016168A (ja) | 荷重分布計測システム、荷重分布計測装置 | |
CN114224313B (zh) | 电阻抗成像方法及计算机可读存储介质 | |
US20210353212A1 (en) | Magnetic resonance imaging apparatus, method, and storage medium | |
JP4617748B2 (ja) | 変換装置 | |
WO2014041580A1 (ja) | 荷重分布計測システム、情報処理装置及び情報処理方法 | |
JP2636757B2 (ja) | 3次元電極位置・頭部モデル推定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170227 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20171121 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180424 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180619 |
|
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: 20180731 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180808 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6385929 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |