CN109862827B - 用于皮下微生物通道的定位的电阻抗测量和eit成像 - Google Patents
用于皮下微生物通道的定位的电阻抗测量和eit成像 Download PDFInfo
- Publication number
- CN109862827B CN109862827B CN201780047721.XA CN201780047721A CN109862827B CN 109862827 B CN109862827 B CN 109862827B CN 201780047721 A CN201780047721 A CN 201780047721A CN 109862827 B CN109862827 B CN 109862827B
- Authority
- CN
- China
- Prior art keywords
- impedance
- electrical
- electrical impedance
- human limb
- acupoint
- 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.)
- Active
Links
- 238000007920 subcutaneous administration Methods 0.000 title description 13
- 238000002847 impedance measurement Methods 0.000 title description 11
- 230000000813 microbial effect Effects 0.000 title description 11
- 238000003384 imaging method Methods 0.000 title description 4
- 230000004807 localization Effects 0.000 title description 3
- 230000037361 pathway Effects 0.000 title description 3
- 230000008859 change Effects 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012546 transfer Methods 0.000 claims abstract description 20
- 238000002604 ultrasonography Methods 0.000 claims description 58
- 230000006870 function Effects 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 5
- 239000003990 capacitor Substances 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 4
- 238000001816 cooling Methods 0.000 claims description 2
- 238000010438 heat treatment Methods 0.000 claims description 2
- 210000001519 tissue Anatomy 0.000 description 32
- 230000004069 differentiation Effects 0.000 description 24
- 210000004027 cell Anatomy 0.000 description 20
- 230000003834 intracellular effect Effects 0.000 description 18
- 239000006185 dispersion Substances 0.000 description 17
- 238000005259 measurement Methods 0.000 description 16
- 239000000523 sample Substances 0.000 description 15
- 230000002159 abnormal effect Effects 0.000 description 10
- 238000010606 normalization Methods 0.000 description 10
- 239000012528 membrane Substances 0.000 description 9
- 210000000170 cell membrane Anatomy 0.000 description 7
- 230000001419 dependent effect Effects 0.000 description 6
- 238000001467 acupuncture Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 239000012620 biological material Substances 0.000 description 3
- 210000003850 cellular structure Anatomy 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000002593 electrical impedance tomography Methods 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 229920000642 polymer Polymers 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 241001334141 Rugopharynx alpha Species 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 239000006285 cell suspension Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 210000003722 extracellular fluid Anatomy 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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/0531—Measuring skin impedance
- A61B5/0532—Measuring skin impedance specially adapted for acupuncture or moxibustion
-
- 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/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4887—Locating particular structures in or on the body
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
-
- 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/046—Arrangements of multiple sensors of the same type in a matrix array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4854—Diagnosis based on concepts of traditional oriental medicine
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Chemical & Material Sciences (AREA)
- Artificial Intelligence (AREA)
- Dermatology (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Alternative & Traditional Medicine (AREA)
- Pain & Pain Management (AREA)
- Electrochemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
一种用于对对象的一个或更多个电阻抗特性的变化进行定位的方法,该方法包括以下步骤:(i)获得对象在不同位置处的电阻抗数据;(ii)使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;以及(iii)识别由多个电阻抗特性中的一个或更多个电阻抗特性的变化所识别的位置。
Description
发明领域
本发明的实施例涉及用于皮下微生物通道(micro bio-channel)的定位的电阻抗测量。
发明背景
电阻抗断层摄影(EIT)是一种已知的成像技术,其特别是在医学和用于检测潜在的异常形态(例如乳腺癌)的其他应用中使用。一般地,多个电极附着到待成像的对象上。将输入电压施加在一小组“输入”电极的两端,且在“输出”电极处测量输出电流,或者将输入电流施加在一小组“输入”电极之间,且在“输出”电极处或在输出电极对之间测量输出电压。例如,当在一小组“输入”电极之间施加非常小的交流电流时,测量在输出电极之间或在“输出”电极对之间的电位差。然后在不同的一小组“输入”电极之间施加电流,并测量在输出电极之间或在“输出”电极对之间的电位差。然后,可以使用适当的图像重建技术来构建基于电阻抗的变化的电阻抗图像。
发明的简要描述
根据本发明的一些实施例,提供了一种用于对对象的一个或更多个电阻抗特性的变化进行定位的方法,该方法包括以下步骤:
(i)获得对象在不同位置处的电阻抗数据;
(ii)使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;以及
(iii)确定由多个电阻抗特性中的一个或更多个电阻抗特性的变化所识别的位置。
与所测量的电阻抗数据相关的电阻抗特性可以从所测量的电阻抗数据导出,并且这些电阻抗特性可以用于定位人体形态,诸如例如皮下的微生物通道和/或低液压阻力的皮下微体积(sub-cutaneous micro-volume)和/或经络和/或穴位。
根据本发明的一些实施例,提供了一种用于对对象的一个或更多个参数阻抗值的变化进行定位的方法,该方法包括以下步骤:
(i)获得对象在不同位置处在一定频率范围内的电阻抗数据;
(ii)使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;
(iii)构造性地组合(constructively combining)所确定的多个电阻抗特性中的所选择的电阻抗特性,以提供对象的随位置变化的至少一个参数阻抗值;以及
(iii)确定由对象的至少一个参数阻抗值的空间变化识别的位置。
电气模型可以设想串联连接的第一和第二阻抗,它们与第三阻抗并联连接。电气模型可以设想电容器和串联连接的电阻器,它们与另一电阻器并联连接。电气模型可以是在任何分辨率下可用的分形模型。
电阻抗特性可选自包括下列项的组:
在频率下限处的阻抗,
在频率上限处的阻抗,
弛豫频率fr,在该频率处存在阻抗的变化,
在该弛豫频率处的阻抗,以及
在该弛豫频率处的阻抗梯度。
电气模型可以设想电容和串联连接的电阻,它们与并联电阻并联连接以形成具有弛豫频率的模型电路,其中用于成像的参数阻抗值是电容、弛豫频率、串联电阻和并联电阻中的两个或更多个的组合。
电气模型可以设想“膜”电容和串联连接的细胞内电阻,它们与细胞外电阻并联连接,其中参数阻抗值包括下列项之一:膜阻抗、膜电导率、细胞内阻抗乘积、细胞内阻抗差值、细胞内阻抗归一化差值、细胞内阻抗微分、细胞内阻抗归一化微分、细胞内电导率乘积、细胞内电导率差值、细胞内电导率归一化差值、细胞内电导率微分、细胞内电导率归一化微分、细胞内时间常数、细胞内频率常数、细胞外阻抗乘积、细胞外阻抗差值、细胞外阻抗归一化差值、细胞外阻抗微分、细胞外阻抗归一化微分、细胞外电导率乘积、细胞外电导率差值、细胞外电导率归一化差值、细胞外电导率微分、细胞外电导率归一化微分、细胞外时间常数、细胞外频率常数、细胞外-内阻抗乘积、细胞外-内阻抗差值、细胞外-内阻抗归一化差值、细胞外-内微分、细胞外-内归一化微分、细胞外-内电导率乘积、细胞外-内电导率差值、细胞外-内电导率归一化差值、细胞外-内电导率微分、细胞外-内电导率归一化微分、由色散梯度α修改的任何一个前述参数。
电气模型可以设想第一阻抗和串联连接的第二阻抗,它们与第三阻抗并联连接以形成具有弛豫频率的模型电路,其中参数阻抗值是第一阻抗、弛豫频率、第二阻抗和第三阻抗中的两个或更多个的组合。
电气模型可以设想内含物边界阻抗和串联连接的内含物内阻抗(intra-inclusion impedance),它们与内含物间阻抗(inter-inclusion impedance)并联连接,其中参数阻抗值包括下列项之一:内含物边界阻抗、内含物边界电导率、内含物内阻抗乘积、内含物内阻抗差值、内含物内阻抗归一化差值、内含物内阻抗微分、内含物内阻抗归一化微分、内含物内电导率乘积、内含物内电导率差值、内含物内电导率归一化差值、内含物内电导率微分、内含物内电导率归一化微分、内含物内时间常数、内含物内频率常数、内含物间阻抗乘积、内含物间阻抗差值、内含物间阻抗归一化差值、内含物间阻抗微分、内含物间阻抗归一化微分、内含物间电导率乘积、内含物间电导率差值、内含物间电导率归一化差值、内含物间电导率微分、内含物间电导率归一化微分、内含物间时间常数、内含物间频率常数、内含物间-内阻抗乘积、内含物间-内阻抗差值、内含物间-内阻抗归一化差值、内含物间-内阻抗微分、内含物间-内归一化微分、内含物间-内电导率乘积、内含物间-内电导率差值、内含物间-内电导率归一化差值、内含物间-内电导率微分、内含物间-内电导率归一化微分、由色散梯度α修改的任何一个前述参数。
频率范围可以在0和20MHz之间。频率范围可以在0和100MHz之间。
根据本发明的一些实施例,提供了一种用于对对象的参数阻抗值的变化进行定位的方法,该方法包括以下步骤:
(i)获得对象在不同位置处的电阻抗数据;
(ii)分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;
(iii)构造性地组合来自所述多个电阻抗特性中的所选择的电阻抗特性,以提供对象的随位置变化的参数阻抗值;
(iv)确定由对象的一个或更多个参数阻抗值的空间变化识别的位置。
对象的电阻抗数据对于生物材料可以在0和20MHz-100MHz之间的频率带宽的情况下被收集。
步骤(iii)可以包括根据阻抗加强算法(impedance emphasising algorithm)来组合预定电阻抗特性。
步骤(i)可以包括获得对象在多个频率处的电阻抗数据。对于生物材料,传递函数可以在0-100MHz的频率范围内由Cole-Cole公式[Cole,1920;Cole,1924]给出。
该方法可用于分析具有细胞结构或细胞状结构的导电对象,并且步骤(ii)可包括使用等效电阻抗电路来模拟该结构,例如Cole-Cole模型[Cole,1920;Cole,1924]。
在限制情况下,等效电阻抗电路可以包括细胞膜电容(C)、细胞内电阻(Ri)和细胞外电阻(Re)。其中等效电阻抗电路包括与细胞内电阻(Ri)串联的细胞膜电容(C),细胞膜电容(C)和细胞内电阻(Ri)与细胞外电阻(Re)或等效电路并联。
电阻抗特性可以选自由Ri(细胞/组(Group)内电阻)、Re(细胞/组外电阻)、C(细胞/组电容)、fr(细胞/组弛豫频率)和α(细胞/组弛豫因子)。
步骤(iii)可以包括通过乘法来组合fr(弛豫频率)和C(细胞/组电容),这可以提供参数阻抗值。
超声换能器可以通过向身体组织施加第一超声信号、接收身体组织的超声响应信号特征以及提供表示超声响应信号的第二输出信号来执行超声检测。
用于获得电阻抗数据的超声换能器和电极阵列可以安装在设备的可移动元件上。
间隔构件可介于可移动元件和受试者之间。
间隔构件可包括一个或更多个孔。间隔构件的孔可以被配置成在使用时与电极阵列的电极对准。间隔构件可以是超声透明的。间隔构件可以是不导电的。
可旋转元件可以具有用于超声换能器的窗口(例如,在可旋转元件中的孔)。电极阵列的电极可以被支撑在窗口上或者可以不存在于窗口上。
根据本发明的一些实施例,提供了一种计算机程序,其提供用于使处理器执行上面定义的方法中的任一个的指令。
根据本发明的一些实施例,提供了一种包括用于执行上面定义的方法中的任一个的装置的系统或设备。
附图的简要描述
为了更好地理解本发明,现在将仅通过示例参考附图,其中:
图1是电阻抗测量设备的图解说明;
图2A和2B示出了在单个或多个色散的情况下的作为频率的函数的所测量的电阻抗的曲线图;
图3示出了在“微观标度”下具有细胞或细胞状结构的对象的示例电阻抗电路模型;
图4示出了在“宏观标度”下具有细胞或细胞状结构的对象的一般电阻抗电路模型;
图5是类似于图1所示的电阻抗测量设备但包括超声检测器的电阻抗测量设备的图解说明;
图6A、6B、6C、6D示出了电极阵列的不同示例;
图7是用于定位皮下微生物通道、经络或穴位的电阻抗测量设备的示例;以及
图8是用于定位皮下微生物通道、经络或穴位的电阻抗测量设备的示例。
发明的实施例的详细描述
图1图解地示出了用于测量对于负载12的阻抗数据的电阻抗测量或电阻抗断层摄影(EIT)设备10。负载12包括附着有多个电极的导电对象。术语“导电的”意味着对象能够传导电流,但是它不一定需要很好地传导电流。对象可以是被称为受试者的活体动物,例如,被称为人类受试者的人类。
设备10还包括信号控制器112,信号控制器112包括信号源14和信号检测器16以及计算机18。在一个实施例中,信号源提供电流作为输入信号,且信号检测器检测电压作为输出信号。在另一个实施例中,信号源提供电压作为输入信号,且信号检测器检测电流作为输出信号。
计算机18通常至少包括处理器和存储器。存储器存储当被加载到处理器中时控制计算机的计算机程序。
使用源14经由电极来将输入信号施加到对象,并且使用检测器16来测量存在于相同或其它电极处的合成输出信号。对于不同频率的输入信号,重复这个过程。例如,电信号可以由信号源14在0Hz(直流)和20MHz或100MHz之间的多个频率下被施加,以使能够获得对象的频率相关电阻抗数据。
用于阻抗测量的电极的间隔确定对象被分析时的分辨率或标度。可以在预期的感兴趣标度(例如,微米或毫米范围)下获得电阻抗测量。作为感兴趣标度的示例,对于生物对象,我们可能对单细胞或组细胞水平或组织或组织水平感兴趣。随后,将使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性。所使用的电气模型可以取决于阻抗测量的分辨率/标度。
设备10还包括用于将由控制器112使用的电极移动到相对于负载12的不同位置处的系统13。该移动可以是物理地将电极平移到不同的物理位置的物理移动,或者可以包括重新配置由电极阵列所使用的电极,使得存在电极阵列的有效移动。
在一个位置处确定对象的多个电阻抗特性,然后改变位置,并再次确定对象的多个电阻抗特性。因此,对象的多个电阻抗特性可以在多个不同位置中的每一个位置处被自动确定。设备10然后根据多个电阻抗特性中的一个或更多个电阻抗特性的空间变化来识别位置。例如,皮下微生物通道、经络或穴位通过多个电阻抗特性中的一个或更多个电阻抗特性的预定空间变化来识别。
参考图2A和2B,可以将使用上述方法获得的电阻抗数据绘制为频率的函数。该曲线22表示对象的阻抗变化与频率的关系或传递函数。计算机18可操作来执行适当的算法以分析获得的阻抗传递函数或频率相关阻抗特性,并从而确定对象在不同位置处的多个电阻抗特性。
电阻抗特性通常包括下列项中的一个或更多个:
a)在极限ω->0(下限)处的阻抗
b)在极限ω->∞(上限)处的阻抗
c)(i)存在阻抗的变化时的弛豫频率
(ii)在该变化频率处的阻抗
(iii)阻抗的变化的梯度,特别是在弛豫频率处;
例如,如果在所使用的频率范围内存在包括生物材料的α、β和γ色散的N个色散[Cole K S,Permeability and impermeability of cell membranes for ions.ColdSpring Harbor Symp.Quant.Biol.8 pp110-22,1940],其中N>1,则色散频率ω1,ω2,...ωN-1,ωN被识别出,并且对于特定色散m的电阻抗特性通常包括下列项中的一个或更多个:
a)对于m=1,在下(全局)限ω->0处的阻抗
对于m>1,在下(局部)限ω->ωm-a处的阻抗,其中a<(ωm-ωm-1)
并且可能是1/2(ωm-ωm-1)
b)对于m=N,在上(全局)限ω->∞处的阻抗
对于m<N,上(局部)限ω->ωm+b处的阻抗,其中b<(ωm+1-ωm),
且可能是b~1/2(ωm+1-ωm)
c)(i)存在阻抗的变化时的弛豫频率ωm(frm)
(ii)在该变化频率处的阻抗
(iii)变化的梯度
这些阻抗特性中的一个或更多个的变化量可用于分析由于细胞内/外或细胞状物内/外相关变化引起的对象的结构。
在一些实施例中,使用等效电阻抗电路来对在分析下的对象建模。可以使用图3所示的等效电阻抗电路20来对对象建模。可以使用等效电阻抗电路20来建模的对象可以包括人体或动物组织。
在所示实施例中,等效电阻抗电路20包括与细胞外部分23并联的细胞部分21。细胞部分21具有串联的电容C和电阻Ri。电阻C与细胞膜/边界相关联,且电阻Ri与细胞的内部相关联。细胞外部分23具有电阻Re。电阻Re与在细胞外部的结构相关联。电阻Re与串联连接的电容C和电阻Ri并联连接。
对于该电路的单个色散阻抗传递函数的非限制性示例是:
在极限ω->0中,Z->Ri
在极限ω->∞中,Z->Ri//Re,即RiRe/(Ri+Re)
在频率fr处和阻抗Zr存在具有梯度α的变化(色散)。
对于在生物组织中的多个色散的传递模型可由Cole-Cole方程(Cole K S 1940,Cole K S 1941,McAdams E T等人,1995)建模如下:
Z=R∝+(R0-R∝)/(1+(jf/fr))(1-α)
通常,如果三元件等效电路用于简单建模细胞悬液(Fricke和Morse,1925)或组织,那么这个方程可以被改写为下面的方程:
Z=Re·Ri/(Re+Ri)+(Re-Re·Ri/(Re+Ri))/(1+(jf/fr))(1-α)
其中,R∝是使Re和Ri并联的结果。
有在频率fri处和阻抗Zri具有梯度αi的变化(色散)。
如上面所指示的,计算机18可操作来执行适当的算法以分析所测量的阻抗数据,并提取在每个位置处的在分析下的对象的多个电阻抗特性。例如,基于所测量的阻抗数据,该算法可操作来绘制作为频率的函数的阻抗数据点,并使用该模型来产生最佳拟合线22,以形成图2所示的传递函数。从该传递函数中,计算机18能够确定对象的多个单独阻抗特性。这些阻抗特性可包括:
a)在极限ω->0处的阻抗,其给出Re
b)在极限ω->∞处的阻抗,其给出Ri Re/(Ri+Re)
c)(i)存在阻抗的变化时的弛豫频率fr
(ii)在该变化频率处的传递函数的阻抗Zr
(iii)给出弛豫因子的变化的梯度α。
阻抗特性可用于使用模型来确定另外的阻抗特性。
例如,如果Re和RiRe/(Ri+Re)都是已知的,那么可以确定Ri。
在变化(色散)频率fr处的传递函数的阻抗Zr是电容器支配传递特性的情况,因为随着频率的每一次小的增加,它明显更好地进行减小阻抗。在变化(色散)频率fr处的阻抗Zr可以被建模为1/(j.2πfr.C)。因此,C可以被确定为1/(j.2πfr.Zr)。
单独阻抗特性(Re,Ri,fr,Zr,α,C)随位置的变化可用于定位形态,如皮下微生物通道、经络线或穴位。
然而,单独阻抗特性的变化量可能不足以实现按位置对结构的准确微分分析。例如,细胞膜电容(C)或弛豫频率(fr)的变化量可能不足以例如基于那些单独的阻抗特性而容易地可检测。
在本发明的实施例中,所选择的预定阻抗特性被“构造性地”组合以提供对象的参数阻抗值。以这种方式提供参数阻抗值的阻抗特性的构造性组合加强单独电阻抗特性的变化。这使对象的结构能够更准确地被分析。
举一个简单的例子,如果存在电阻抗特性之一的10%的增加(例如细胞膜电容(C))从初始值C1到1.1C1)和另一电阻抗特性的10%的增加(例如弛豫频率(fr)从初始值fr1到1.1fr1),则这些单独的10%增加可能不足以容易地可检测,例如在这些单独电阻抗特性的测量中可辨别。然而,通过乘法来组合这些单独电阻抗特性以提供参数阻抗值将导致21%(1.21fr1C1)的更大增加,这是更容易检测的。
阻抗特性可以具有与特定形态的正相关、中性相关或负相关。正相关意味着它在形态存在时增加,尽管可能不明显。负相关意味着它在形态存在时减小,尽管可能不明显。中性相关意味着它在形态存在时不改变。具有正相关的阻抗特性可以通过取反数来转换成具有负相关的阻抗特性(反之亦然)。
用于检测特定形态的阻抗特性的构造性组合意味着通过乘法(或加权相加)来组合在相同意义上对该形态关联的阻抗特性以创建参数阻抗值,以及通过除法(或加权相减)来组合在相反意义上对该形态关联的阻抗特性。
任何所确定的阻抗特性可以以任何期望的方式构造性地组合以提供参数阻抗值,该参数阻抗值比任何组成阻抗特性具有对形态变化的更大敏感度。
在图3中描述的极限水平处的阻抗特性的组合的非限制性示例:
组合参数测量
(组合的内/外/膜阻抗/电导率)
a)膜阻抗/电导率和相关量:
膜阻抗:Zm=1/2π*fr*C
膜电导率:σm=2π*fr*C
b)组合的细胞内阻抗/电导率:
乘积或商:Ri*Zm或Ri/Zm或Ri/C或RiC
或:σi*σm
差值/归一化差值:
a*Ri-b*Zm
或:c*σi-d*σm
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Ri-b*Zm)/Zm
或:(a*Ri-b*Zm)/Ri
可选地:(c*σi-d*σm)/σm
或:(c*σi-d*σm)/σi
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
细胞内时间常数:Ri*C
或:细胞内频率常数1/Ri*C
c)组合的细胞外阻抗/电导率:
乘积或商:Re*Zm或Re/Zm或Re/C或ReC
或:σx*σm
差值/归一化差值:
a*Re-b*Zm
或:c*σx-d*σm
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Re-b*Zm)/Zm
或:(a*Re-b*Zm)/Re
可选地:(c*σx-d*σm)/σm
或:(c*σx-d*σm)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
细胞外时间常数:Re*C
或:细胞外频率常数1/Re*C
d)组合的细胞外至内阻抗/电导率:
乘积:Re*Ri
或:σx*σi
商:Re/Ri或Ri/Re
差值/归一化差值:
a*Re-b*Ri
或:c*σx-d*σi
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Re-b*Ri)/Ri
或:(a*Re-b*Ri)/Re
可选地:(c*σx-d*σi)/σi
或:(c*σx-d*σi)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
e)三重组合:
ReRiC 或 ReRi/C或
(Ri/Re)C 或 (Ri/Re)/C或
(Ri*Re)Zm 或 (Ri*Re)/Zm或
具有异常色散特性(α)的组合集成细胞参数测量
a)“异常”膜阻抗/电导率和相关量:
“异常”膜阻抗:
α*Zm
或:α/Zm
“异常”膜电导率:
α*σm
或:α/σm
b)组合的“异常”细胞内阻抗/电导率:
乘积: α*Ri*Zm
或:σi*σm
差值/归一化差值:
α*(a*Ri-b*Zm)
或:α*(c*σi-d*σm)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Ri-b*Zm)/Zm
或:α*(a*Ri-b*Zm)/Ri
可选地:α*(c*σi-d*σm)/σm
或:α*(c*σi-d*σm)/σi
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
细胞内时间常数:α*(Ri*C)
或:细胞内频率常数α*(1/Ri*C)
c)组合的细胞外阻抗/电导率:
乘积:α*Re*Zm
或者:α*σx*σm
差值/归一化差值:
α*(a*Re-b*Zm)
或:α*(C*σx-d*σm)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Re-b*Zm)/Zm
或:α*(a*Re-b*Zm)/Re
可选地:α*(C*σx-d*σm)/σm
或:α*(c*σx-d*σm)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
细胞外时间常数:α*Re*C
或:细胞外频率常数α*(1/Re*C)
d)组合的细胞外至内阻抗/电导率:
乘积:α*Re*Ri
或:α*σx*σi
差值/归一化差值:
α*(a*Re-b*Ri)
或:α*(c*σx-d*σi)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Re-b*Ri)/Ri
或:α*(a*Re-b*Ri)/Re
可选地:α*(C*σx-d*σi)/σi
或:α*(C*σx-d*σi)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
e)三重组合:
α*ReRiC 或 α*ReRi/C 或
α*(Ri/Re)C 或 α*(Ri/Re)/C 或
α*(Ri*Re)Zm 或 α*(Ri*Re)/Zm 或
具有色散频率Fr的组合集成细胞参数测量。用Fr直接代替上面的在a)至e)中的α。
具有色散频率Fr和色散特性(α)的组合集成细胞参数测量。用α*Fr直接代替上面的在a)至e)中的α。
合适的阻抗加强算法可由计算机18实现,以选择用于组合的最佳电阻抗特性及它们的组合方式,以最大化因而得到的参数阻抗值的变化。
对象的参数阻抗值在一个位置处被获得后,它在其他位置处被获得。该方法然后通过对象的至少一个参数阻抗值的变化来识别位置。例如,皮下微生物通道、经络或穴位通过对象的一个或更多个参数阻抗值的预定空间变化来识别。
图4示出了在分析下的对象的更一般的模型。在所示实施例中,等效电阻抗电路30包括与内含物间部分33并联的内含物部分31。内含物部分31具有串联的阻抗Z1和阻抗Z2。阻抗Z1可以与内含物边界相关联(可以代表一组细胞的膜相关成分),以及阻抗Z2可以与内含物的内部相关联(可以代表一组细胞的细胞内相关成分)。内含物间部分33具有阻抗Z3。阻抗Z3与在内含物外部的结构相关联(可以代表一组细胞的细胞外成分)。阻抗Z3与串联连接的阻抗Z1和Z2并联连接。
对于该电路的阻抗传递函数为:
在图4中描述的水平处的阻抗特性的组合的非限制性示例:
组合参数测量
a)内含物边界阻抗/电导率和相关量:
内含物边界阻抗:Zm=1/2π*fr*Z2
内含物边界电导率:σm=2π*Fr*Z2
b)组合的内含物内阻抗/电导率:
乘积:Z1*Zm
或:σ1*σm
差值/归一化差值:
a*Z1-b*Zm
或:c*σ1-d*σm
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Z1-b*Zm)/Zm
或:(a*Z1-b*Zm)/Z1
可选地:(c*σ1-d*σm)/σm
或:(c*σ1-d*σm)/σ1
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
内含物内时间常数:Z1*Z2
或:内含物内频率常数 1/Z1*Z2
c)组合的内含物间阻抗/电导率:
乘积:Z3*Zm
或:σ3*σm
差值/归一化差值:
a*Z3-b*Zm
或:c*σ3-d*σm
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
d)组合的内含物间到内阻抗/电导率:
乘积:Re*Ri
或:σx*σi
差值/归一化差异:
a*Re-b*Ri
或:c*σx-d*σi
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Re-b*Ri/Ri
或:(a*Re-b*Ri)/Re
可选地:(c*σx-d*σi)/σi
或:(c*σx-d*σi)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
(a*Z3-b*Zm)/Zm
或:(a*Z3-b*Zm)/Z3
可选地:(c*σ3-d*σm)/σm
或:(c*σ3-d*σm)/σ3
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
内含物间时间常数:Z3*Z2
或:内含物间频率常数1/Z3*Z2
具有异常色散特性(α)的组合集成参数测量
a)“异常”内含物边界阻抗/电导率和相关量:
“异常”内含物边界阻抗:
α*Zm
或:α/Zm
“异常”内含物边界电导率:
α*σm
或:α/σm
b)组合的“异常”内含物内阻抗/电导率:
乘积:α*Z1*Zm
或:σ1*σm
差值/归一化差值:
α*(a*Z1-b*Zm)
或:α*(c*σ1-d*σm)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Z1-b*Zm)/Zm
或:α*(a*Z1-b*Zm)/Z1
可选地:α*(c*σ1-d*σm)/σm
或:α*(c*σ1-d*σm)/σ1
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
内含物内时间常数:α*(Z1*Z2)
或:内含物内频率常数α*(1/Z1*Z2)
c)组合的内含物间阻抗/电导率:
乘积:α*Z3*Zm
或:σ3*σm
差值/归一化差值:
α*(a*Z3-b*Zm)
或:α*(c*σ3-d*σm)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Z3-b*Zm)/Zm
或:α*(a*Z3-b*Zm)/Z3
可选地:α*(c*σ3-d*σm)/σm
或:α*(c*σ3-d*σm)/σ3
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
内含物间时间常数:α*Z3*Z2
或:内含物间频率常数α*(1/Z3*Z2)
d)组合的内含物间到内阻抗/电导率:
乘积:α*Re*Ri
或:α*σx*σi
差值/归一化差值:
α*(a*Re-b*Ri)
或:α*(c*σx-d*σi)
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量;
微分/归一化微分:
α*(a*Re-b*Ri)/Ri
或:α*(a*Re-b*Ri)/Re
可选地:α*(c*σx-d*σi)/σi
或:α*(c*σx-d*σi)/σx
其中,系数a、b、c和d是常数(-∞-+∞),以用于匹配待使用的量。
这个模型是如在US6856824中先前所述的分形模型。阻抗Z1、Z2、Z3中的每一个可以使用电路30或者在极限水平处表示,其中Z1等效于Ri,Z2等效于C,以及Z3等效于Re。术语“分形”用来表达下面的事实:无论在我们看结构的哪个维度级别处,模型都是一样的。
参考图5,示出了用于检测包括电极阵列101和超声探头102的身体组织的信号特征的设备10。
电极阵列101包括设置在电极板105的表面104上的多个电极103。在使用中,身体组织(未示出)相邻于电极板105,相邻于电极板105的表面104,或者与表面104接触,或者与表面104间隔开。电极103能够在身体组织上的电阻抗测量期间向身体组织施加第一电信号。电极103电耦合到第一控制器111,第一控制器111用于将第一电信号发送到电极103以施加到身体组织,并且用于从电极103接收第一输出信号,该第一输出信号取决于在电极103处接收的电响应信号(身体组织的特征)。
超声探头102包括设置在超声探头102的表面106上的多个超声换能器107。超声换能器107能够在身体组织上的超声检查期间向身体组织施加第一超声信号。超声换能器107电耦合到第二控制器112,其用于向超声换能器107提供使超声换能器107将第一超声信号施加到身体组织的、通常以电脉冲的形式的第二输入信号,并且用于从超声换能器107接收第二输出信号,该第二输出信号取决于在超声换能器107处接收的超声响应信号(身体组织的特征)。
超声探头102的设置有超声换能器107的表面106与电极板105相邻,并且在电极板105的与电极板105的设置有电极103的表面104相对的一侧上。因此,如果电极板105水平地放置,其中电极板105的其上设置电极103的表面104朝上,那么超声探头102在电极板105下方,其中超声探头102的其上设置了超声换能器107的表面106也朝上。因此,超声换能器107布置在实质上平行于电极板105的平面中。这使电信号和超声信号能够在实质上彼此平行的方向上施加到身体组织。
电极板105的移动使在比电极103被部署于其上的电极板105的区域更大的身体组织的区域上的不同位置处能够施加第一电信号且能够检测电响应信号。系统13可以引起电极板105和超声探头102相对于受试者的平移,如前面参考图1所述的。
在这个示例中,超声探头102也相对于受试者旋转。超声探头102和电极板105机械地耦合,由此,超声探头102围绕实质上垂直于电极板105的轴108可旋转。对于给定尺寸的身体组织的区域,可以部署更少的电极103,这可以降低电连接的复杂性。电极板105的旋转还使用比电极103上的物理间隔更接近地间隔开的电极103的增量位置以精细分辨率实现电测量。
电极103耦合到第一端口109,以及超声换能器107耦合到第二端口110。第一端口109是双向的,用于向电极103传送信号和从电极103传送信号。第二端口110也是双向的,用于向超声换能器107传送信号和从超声换能器107传送信号。为了清楚起见,在图1中未示出在第一端口109和电极103之间以及在第二端口110和超声换能器107之间的连接。这些连接可以例如定位在电极板105的与表面104相对的表面上,或者可以在电极板105的内部。
存在耦合到第一端口109的第一控制器111。第一控制器111生成经由第一端口109传送到电极103中的一个或更多个电极的第一输入信号,其中响应于第一输入信号,第一电信号被传输到身体组织。第一电信号穿过身体组织并在电极103中的另一个处被接收。在本说明书和所附权利要求中,这些所接收的信号被称为电响应信号。根据电响应信号,第一输出信号经由第一端口109被传送到第一控制器111。
存在耦合到第二端口110的第二控制器112。第二控制器112生成经由第二端口110传送到超声换能器107的第二输入信号。第二输入信号可以是例如电信号或光信号。超声换能器107将第二输入信号转换成传输到身体组织的第一超声信号。第一超声信号在身体组织中被反射。在本说明书和所附权利要求中,这些反射被称为超声响应信号。超声响应信号由超声换能器107检测,超声换能器107将超声响应信号转换成经由第二端口110传送到第二控制器112的第二输出信号。
第一控制器111和第二控制器112耦合到计算机18。计算机18基于第一输出信号来生成电阻抗数据,并且基于第二输出信号来生成超声数据。超声数据和电阻抗数据是身体组织的特征。
可选地,计算机18可以组合超声数据和电阻抗数据,且显示器114(可选的)可以显示表示组合的超声数据和电阻抗数据的图像。通过这种手段,电阻抗数据和超声数据可以被组合以提供增强的图像,其可以帮助检测和定位身体组织的特征,例如皮下微生物通道、经络线和穴位的特征。仅从电阻抗数据或超声数据可能不明显的身体组织的特征在电阻抗数据和超声数据的组合后可能变得明显。图像可以是二维或三维的。
尽管图5所示的电极板105是平坦的,但这不是本发明的本质特征,并且电极板105或者至少表面104可以是非平坦的。例如,表面104可以以类似于身体组织的形状形成轮廓。这使身体组织的形状的变形能够被减少或避免。超声探头102的表面106和超声换能器107的布置可以被形成轮廓以与相邻电极板105的形状互补。
通过采用与身体组织的形状互补的形状,可以减小在电极板105和身体组织之间以及在超声换能器107和身体组织之间的信号路径的长度,导致在检测响应信号时设备的提高的灵敏度。
在图4所示的实施例中,电极板105是圆形的。这不是本发明的本质特征,且可以使用其它形状。
经由第一端口109和第二端口110传送的信号可以是电流或电压,或者可以是光信号。此外,它们可以是模拟或数字信号。在光信号被使用的情况下,在光信号和电信号之间的转换可以由超声换能器107、由电极103以及由第一控制器111和第二控制器112执行。数模转换和模数转换可以由超声换能器107、由电极103以及由第一控制器111和第二控制器112执行。超声换能器107、电极103以及第一控制器111和第二控制器112可以包括信号处理,例如放大和滤波。第一控制器111可以与电极板105构成整体,以及第二控制器112可以与超声探头102构成整体,在这种情况下,第一端口109和第二端口110中的一个或两个可以分别在电极板105或超声探头102的内部。可选地,第一控制器111可以通过电缆与电极板105间隔开,和/或第二控制器112可以通过电缆与超声探头102间隔开。
第一控制器111和第二控制器112可以耦合,并且实际上可以是公共控制器。这使第一和第二信号的生成能够同步。例如,可以控制第一和第二信号的相对定时和/或幅度。
设备10可以包括用于将在评估下的对象与电极载板105间隔开的间隔构件。当电极载板105旋转时,间隔构件和对象不旋转。以这种方式,对象被保护免受来自旋转载板105的旋转力,并且对患者的不适可以被减少或消除。可旋转载板105可具有用于所附接的超声探头106的超声窗口。窗口可以是穿过板105的孔。间隔膜800可以由超声透明膜(诸如例如聚合物,例如高电阻抗聚合物)制成。在该示例中,每个电极可以具有穿过间隔构件的孔。
图6A、6B、6C、6D公开了电极103的电极阵列101的示例。这些电极阵列101是可旋转的并且具有旋转对称性。间隔构件用于将在评估下的对象与旋转电极阵列101间隔开。载板105是不导电的,并且可以由例如塑料材料制成。电极103部署在电极载板105的整个平坦表面上,并且优选地凹入电极载板105中,使得它们不与放置在电极载板105上的对象产生物理接触。每个电极由一个点表示。
电极在包括重复的单元的布置中被布置在电极载板105上,重复的单元在电极载板105上重复,并且具有小于或等于90°的旋转对称角。电极被部署成使得能够使用通过连续位置旋转小于90°的旋转位移的电极的图案来进行电阻抗的测量。
在图6A中,电极103在正方形矩阵中等距地布置,使得电极位于布置在连续体中的正方形的角处。电极布置在正方形的棋盘形格局的每个正方形的一个或更多个角处。这种布置实现旋转对称,其为90°的倍数。在这种布置中,除了相邻于该布置的边界的电极之外,每个电极具有布置在正方形中的四个最接近的邻近电极。可选地,可以通过将每个正方形细分成四个较小的正方形来提供更密集的矩阵。
电极103在X列×N行的阵列中,其中X=1、2、3、...N;Y=1、2、…M;例如3×16阵列、3×24阵列;5×16阵列、5×24阵列,在电极之间的间隙可以是任何可能的值,例如0.01mm、0.02mm、0.05mm、0.1mm、0.5mm、1mm等。
在图6B中,电极103在三角形矩阵中等距地布置,使得电极位于在连续体中布置的等边三角形的角处。电极布置在三角形的棋盘形格局的每个三角形的一个或更多个角处。三角形可以是等边三角形。此外,三角形可具有相等的尺寸。这种布置实现旋转对称,其为60°的倍数。在这种布置中,除了相邻于该布置的边界的电极之外,每个电极具有布置在六边形中的六个最接近的邻近电极。可选地,可以通过将每个等边三角形细分成四个较小的等边三角形来提供更密集的三角形矩阵。
电极103在一半占据的X列×N行阵列中,X=1、2、3、...N;Y=1、2、…M;例如3×16阵列、3×24阵列;5×16阵列,5×24阵列,在电极之间的间隙可以是任何可能的值,例如0.01mm、0.02mm、0.05mm、0.1mm、0.5mm、1mm等。
在图6C和6D中,电极103布置在穿过公共点O的五条或更多条径向延伸的线上。在二维中,N条径向线具有360°/N的规则间隔,并且在电极103之间沿径向线的间隙可以是任何可能的值,例如0.01mm、0.02mm、0.05mm、0.1mm、0.5mm、1mm等。
每个电极的直径可以是1-5μm,但不限于此范围。在最接近的电极之间的中心距离可以在例如0.01-1mm之间。
电极103可以被配置成在皮下的高达5-10mm,通常以1-2mm分辨率在皮下的2-4mm进行测量。
图7示出了除电阻抗测量之外还可以使用或不使用超声测量的设备10的示例。电极载板105可以相对于对象(受试者)横向移动200。该移动可以由系统13(未示出)引起,该系统例如可以包括伺服马达。设备10移动电极103的阵列,并在每个位置处进行测量。
设备:
(i)获得对象在不同位置处的电阻抗数据;
(ii)使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;以及
(iii)确定由多个电阻抗特性中的一个或更多个的变化所识别的位置。
在识别了多个电阻抗特性中的一个或更多个电阻抗特性的变化的位置后,设备10移动电极103的载板105,使得在载板105的中间的孔202位于对应于皮下微生物通道、经络或穴位的位置上。这允许针灸师确定在孔202内的穴位。
设备10可以被配置为沿着肢体移动的带或绳。然后横向运动越过肢体。
对象的电阻抗数据可以在多个不同位置中的每一个处在频率范围内被获得。
对象的随位置变化的参数阻抗值可以通过构造性地组合所确定的多个电阻抗特性中的所选择的特性来形成,并且感兴趣的位置可以通过对象的至少一个参数阻抗值的变化来识别。
超声传感器可以位于孔202内。
图8示出了除电阻抗测量之外还可以使用或不使用超声测量的设备10的示例。电极载板105被分成两个部分105A和105B。
部分105A是电极103的阵列101。该阵列可以是固定的或可移动的。例如,它可以相对于对象(受试者)横向移动200。该移动可以由系统13(未示出)引起,该系统例如可以包括伺服马达。
部分105B是电极103的阵列101。该阵列可以是可移动的,并且操作读取头210。在所示的示例中,存在两个读取头210,例如,它可以相对于对象(受试者)横向和/或纵向移动。该移动可以由系统13(未示出)引起,该系统例如可以包括伺服马达。
设备10移动电极103的阵列,并在每个位置处进行测量。
设备10:
(i)获得对象在不同位置处的电阻抗数据;
(ii)使用假定的电气模型的传递函数来分析所获得的电阻抗数据,以确定对象的多个电阻抗特性随位置的变化;以及
(iii)确定由多个电阻抗特性中的一个或更多个电阻抗特性的变化所识别的位置。
在识别了多个电阻抗特性中的一个或更多个电阻抗特性的变化的位置后,设备10移动读取头210以确定该位置。
对象的电阻抗数据可以在多个不同位置中的每一个处在一定频率范围内被获得。
对象的随位置变化的参数阻抗值可以通过构造性地组合所确定的多个电阻抗特性中的所选择的电阻抗特性来形成,并且感兴趣的位置可以通过对象的至少一个参数阻抗值的变化来识别。
在前述示例中,温度控制子系统可用于例如经由电极板105或独立地提供受试者的恒定或差温加热/冷却。因此,在不同位置处的测量可以在受控的固定或可变温度条件下进行。
在前面的示例中,在不同位置处的测量可以在没有针灸的情况下进行,或者可以在针灸期间进行。
通过应用时变参数,例如温度变化和针灸的缺乏/存在,可能确定电阻抗特性的时间变化。这可以允许确定参数,例如在针灸刺激期间沿经络线的速度、传输方向和传输感觉。动态测量从在针灸被给与时的数据的变化来检测通道中的间质流体流动。该装置还可以根据红细胞的1-3MHz的特征频率来辨别在经络通道和血管之间的差异。
数据库可用于将受试者的所测量的电阻抗特性与在不同性别、年龄、正常和异常组下的数据进行比较。
尽管在前面的段落中已经参考各种非限制性示例描述了本发明的实施例,但是应当认识到,可以做出对给定示例的修改而不偏离如所要求保护的本发明的范围。
虽然在前述说明书中努力将注意力引到本发明的被认为具有特别的重要性的那些特征,但是应理解,申请人要求关于在上文中提到和/或在附图中示出的任何可取得专利权的特征或特征的组合的保护。
Claims (11)
1.一种用于通过对对象的一个或更多个电阻抗特性的变化进行定位而将载板中的孔定位在穴位上方的方法,所述方法包括以下步骤:
(i)使用一个或更多个伺服马达跨人类肢体物理地移动由所述载板承载的电极阵列的至少一部分或全部,以获得所述对象在不同位置处的电阻抗数据;
(ii)使用处理器分析所获得的电阻抗数据,以使用设想的电气模型的传递函数来基于所述对象的多个电阻抗特性随位置的变化确定所述穴位的位置;以及
(iii)响应于基于所述人类肢体的所述多个电阻抗特性随位置的变化来确定所述穴位的位置,使用所述一个或更多个伺服马达物理地移动由所述载板承载的所述电极阵列的至少一部分或全部,使得所述载板中的孔被定位在所述人类肢体上的所述穴位的所确定的位置上方,所述所确定的位置通过所述多个电阻抗特性的变化而被识别为穴位,其中所述电极阵列的电极不存在于所述孔上方,并且所述孔允许针灸师识别并使用所述孔内的所述穴位。
2.如权利要求1所述的方法,包括:使用用于执行电阻抗检测的电极阵列以通过向所述人类肢体施加第一电信号、接收所述人类肢体的电响应信号特征以及提供表示所述电响应信号的第一输出信号来获得所述人类肢体的所述电阻抗数据。
3.如权利要求1或2所述的方法,包括在沿所述人类肢体的纵向方向和跨所述人类肢体的横向方向二者上物理地移动由所述载板承载的所述电极阵列的至少一部分或全部,以在不同位置处执行电阻抗检测。
4.如权利要求1至3中的任一项所述的方法,其中,所述电极阵列被分成被配置用于相对于彼此移动的两个部分。
5.如权利要求4所述的方法,包括移动所述电极阵列的所述部分之一以向用户确定所述人类肢体上的所述穴位的位置。
6.如任一前述权利要求所述的方法,还包括提供加热/冷却使得所述人类肢体的在不同位置处所获得的电阻抗数据发生在所控制的固定温度条件下。
7.如权利要求1至6中的任一项所述的方法,包括使用超声换能器以通过向所述人类肢体施加第一超声信号、接收所述人类肢体的超声响应信号特征以及提供表示所述超声响应信号的第二输出信号来执行超声检测。
8.如任一前述权利要求所述的方法,其中,所述电气模型是电路模型,所述电路模型设想第一阻抗和第二阻抗串联连接为与第三阻抗并联连接的串联组合,或者
其中,所述电路模型设想电容器和电阻器串联连接为与另一电阻器并联连接的串联组合,
其中所述电阻抗特性选自包括下列项的组:
在频率下限处的阻抗,
在频率上限处的阻抗,
存在阻抗的变化时的弛豫频率fr,
在所述弛豫频率处的阻抗,以及
在所述弛豫频率处的阻抗梯度。
9.如权利要求8所述的方法,其中,所述电气模型是分形模型,并且能够在任何分辨率下使用。
10.如任一前述权利要求所述的方法,其中,获得所述对象在不同位置处的电阻抗数据包括:
在一定频率范围内获得所述对象在不同位置处的电阻抗数据,
其中,所述频率范围在0和100MHz之间并包括1-3MHz。
11.一种用于将承板中的孔定位在穴位上方的装置,包括:
电极阵列,其包括电极;
载板,其包括孔并且承载所述电极阵列;
一个或更多个伺服马达;
处理器和存储器,其存储用于基于人类肢体的电阻抗特性的变化将所述载板中的孔定位在穴位上方的计算机程序,其中所述处理器、所述存储器和所述计算机程序使所述装置能够:
控制所述一个或更多个伺服马达在沿所述人类肢体的纵向方向和跨所述人类肢体的横向方向上物理地移动由所述载板承载的所述电极阵列的至少一部分或全部,以获得所述人类肢体在不同位置处的电阻抗数据;
使用设想的电气模型的传递函数来分析电阻抗数据,以基于所述设想的电气模型的所述人类肢体的多个电阻抗特性随位置的变化来确定所述穴位的位置;和
响应于基于所述设想的电气模型的所述人类肢体的所述多个电阻抗特性随位置的变化确定所述穴位的位置,控制所述一个或更多个伺服马达物理地移动由所述载板承载的所述电极阵列的至少一部分或全部,使得所述载板中的孔被定位在所述人类肢体上的所述穴位的所确定的位置上方,所述所确定的位置通过所述设想的电气模型的所述多个电阻抗特性的变化而被识别为穴位,其中所述电极阵列的电极不存在于所述孔上方,并且所述孔允许针灸师识别并使用所述孔内的所述穴位。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB1609551.5 | 2016-05-31 | ||
GB1609551.5A GB2553261B (en) | 2016-05-31 | 2016-05-31 | Electrical impedance measurement and EIT image for location of a micro bio-channel under skin |
PCT/GB2017/051514 WO2017207972A1 (en) | 2016-05-31 | 2017-05-26 | Electrical impedance measurement and eit image for location of a micro bio-channel under skin |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109862827A CN109862827A (zh) | 2019-06-07 |
CN109862827B true CN109862827B (zh) | 2024-03-22 |
Family
ID=56410807
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780047721.XA Active CN109862827B (zh) | 2016-05-31 | 2017-05-26 | 用于皮下微生物通道的定位的电阻抗测量和eit成像 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20210219863A1 (zh) |
EP (1) | EP3461254A1 (zh) |
JP (2) | JP2019523428A (zh) |
KR (1) | KR20190020701A (zh) |
CN (1) | CN109862827B (zh) |
GB (1) | GB2553261B (zh) |
WO (1) | WO2017207972A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112220472B (zh) * | 2020-09-25 | 2023-09-15 | 重庆大学 | 基于交叉四电极法的胸腔电阻抗无创检测方法 |
CN114880863B (zh) * | 2022-05-20 | 2024-03-01 | 合肥工业大学 | 分布式可再生能源集群阻抗聚合模型自适应分频降阶方法 |
GB2621542A (en) * | 2022-06-23 | 2024-02-21 | Wang Wei | Electrical impedance mammography |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6856824B1 (en) * | 1998-08-28 | 2005-02-15 | De Montfort University | Apparatus and method for detecting abnormalities in bodily matter |
CN101754717A (zh) * | 2007-06-07 | 2010-06-23 | 德蒙特福特大学 | 用于分析电传导对象的结构的方法 |
WO2012076881A1 (en) * | 2010-12-07 | 2012-06-14 | The University Of Sussex | Electrical impedance detection and ultrasound scanning of body tissue |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3859983A (en) * | 1973-06-07 | 1975-01-14 | Albert A Dohring | Acupuncture locating device and needle guide |
JPS51159490U (zh) * | 1975-06-13 | 1976-12-18 | ||
US4016870A (en) * | 1975-10-14 | 1977-04-12 | Chuck Lock | Electronic acupuncture point finder |
JPH0342918Y2 (zh) * | 1985-10-16 | 1991-09-09 | ||
DE10217765A1 (de) * | 2002-04-20 | 2003-11-20 | Hti Gmbh | Verfahren zum Auffinden und/oder Klassifizieren von Akupunkturpunkten auf der Haut |
GB0228375D0 (en) * | 2002-12-05 | 2003-01-08 | Innovation And Entpr Off Of | Wound mapping |
EP1613212B1 (en) * | 2003-03-25 | 2016-08-24 | Fresenius Medical Care Holdings, Inc. | Device and method for performing electrical impedance tomography |
JP4247969B2 (ja) * | 2003-05-15 | 2009-04-02 | 株式会社 東北テクノアーチ | 経穴位置の評価装置 |
WO2006044868A1 (en) * | 2004-10-20 | 2006-04-27 | Nervonix, Inc. | An active electrode, bio-impedance based, tissue discrimination system and methods and use |
JP4881574B2 (ja) * | 2005-04-21 | 2012-02-22 | パナソニック株式会社 | 経穴位置評価装置 |
CN100508880C (zh) * | 2007-09-17 | 2009-07-08 | 中国人民解放军第四军医大学 | 实时多信息提取的电阻抗扫描检测系统 |
US8198998B1 (en) * | 2008-01-23 | 2012-06-12 | Adaptive Networks, Inc. | Identification of power line locations by power line signatures |
KR100965351B1 (ko) * | 2009-11-23 | 2010-06-22 | 박문서 | 인체내 임피던스 측정을 위한 전극 장치를 이용한 인체내 임피던스 측정과 시술 장치 |
GB2524470B (en) * | 2014-02-16 | 2019-04-17 | Wang Wei | An apparatus and method for "high-resolution" electrical impedance imaging |
-
2016
- 2016-05-31 GB GB1609551.5A patent/GB2553261B/en active Active
-
2017
- 2017-05-26 US US16/305,507 patent/US20210219863A1/en active Pending
- 2017-05-26 JP JP2019516278A patent/JP2019523428A/ja active Pending
- 2017-05-26 CN CN201780047721.XA patent/CN109862827B/zh active Active
- 2017-05-26 EP EP17727346.3A patent/EP3461254A1/en active Pending
- 2017-05-26 WO PCT/GB2017/051514 patent/WO2017207972A1/en unknown
- 2017-05-26 KR KR1020187037916A patent/KR20190020701A/ko not_active Application Discontinuation
-
2023
- 2023-05-17 JP JP2023081865A patent/JP2023111912A/ja active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6856824B1 (en) * | 1998-08-28 | 2005-02-15 | De Montfort University | Apparatus and method for detecting abnormalities in bodily matter |
CN101754717A (zh) * | 2007-06-07 | 2010-06-23 | 德蒙特福特大学 | 用于分析电传导对象的结构的方法 |
WO2012076881A1 (en) * | 2010-12-07 | 2012-06-14 | The University Of Sussex | Electrical impedance detection and ultrasound scanning of body tissue |
Also Published As
Publication number | Publication date |
---|---|
EP3461254A1 (en) | 2019-04-03 |
JP2023111912A (ja) | 2023-08-10 |
GB2553261B (en) | 2022-05-04 |
GB2553261A (en) | 2018-03-07 |
GB201609551D0 (en) | 2016-07-13 |
CN109862827A (zh) | 2019-06-07 |
KR20190020701A (ko) | 2019-03-04 |
JP2019523428A (ja) | 2019-08-22 |
US20210219863A1 (en) | 2021-07-22 |
WO2017207972A1 (en) | 2017-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chitturi et al. | Spatial resolution in electrical impedance tomography: A topical review | |
Bera et al. | Electrical impedance spectroscopy (EIS)-based evaluation of biological tissue phantoms to study multifrequency electrical impedance tomography (Mf-EIT) systems | |
JP2023111912A (ja) | 皮下のマイクロ生体チャンネルの場所の特定のための電気インピーダンス測定及びeitイメージ | |
Malone et al. | Stroke type differentiation using spectrally constrained multifrequency EIT: evaluation of feasibility in a realistic head model | |
EP2073702B2 (en) | Data collection for electrical impedance tomography | |
CN109745046B (zh) | 一种适用于运动状态下的电阻抗成像电极与系统 | |
CN103027673A (zh) | 刺激目标范围标定方法 | |
US9927385B2 (en) | Active control guards and rationometric calibration and reconstruction for use with electrical capacitance volume tomography | |
Meir et al. | Electrical impedance tomographic imaging of a single cell electroporation | |
US11992683B2 (en) | Bioimpedance measurement method and apparatus with electrical stimulation performance | |
CN102138790A (zh) | 一种三维电阻抗断层成像系统性能测试方法与装置 | |
Yang et al. | A 3-D reconstruction solution to current density imaging based on acoustoelectric effect by deconvolution: a simulation study | |
US8928332B2 (en) | Electrical impedance imaging | |
JP2019523428A5 (zh) | ||
Canali et al. | Electrical impedance tomography methods for miniaturised 3D systems | |
Liu et al. | Design of a microscopic electrical impedance tomography system using two current injections | |
Zarafshani et al. | A 3D multi-frequency response electrical mesh phantom for validation of the planar structure EIT system performance | |
Sun et al. | A New Method for Electrical Impedance Tomography with Incomplete Electrode Array | |
US20220287567A1 (en) | System and method for measuring tissue parameters by use of capacitive tactile sensor | |
Liu et al. | Detection of pressure ulcers using electrical impedance tomography | |
Lehti-Polojärvi et al. | Multifrequency electrical impedance tomography in biological applications: A multimodal perspective | |
Xiang et al. | Image Reconstruction for Multi-frequency Electromagnetic Tomography based on Multiple Measurement Vector Model | |
He et al. | A combined regularization algorithm for electrical impedance tomography system using rectangular electrodes array | |
Busch et al. | Realization of a multi-layer EIT-system | |
CN110988043B (zh) | 基于多频加权频差的多介质分离成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40003386 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant |