技术领域
[0001] 本
发明涉及
磁感应断层成像,特别是涉及用于提高磁感应断层成像系统的灵敏度的方法和系统。
背景技术
[0002] 磁感应断层成像(MIT)是一种应用于工业和医学成像的非侵入且无
接触的成像技术。与其它电学成像技术对比,MIT不要求
传感器与关注对象直接接触来进行成像。
[0003] MIT被用于重构关注对象内部的无源电学属性的空间分布,该无源电学属性为例如电导率σ、电容率ε和磁导率μ。在MIT中,通常介于几kHz直至若干MHz的正弦
电流被应用到发射线圈,从而生成时变
磁场,该磁场通常称为主磁场。由于关注对象(例如
生物组织)导电,该主磁场在关注对象中产生“
涡流电流”。这些涡流电流生成次磁场。这些磁场的组合在接
收线圈中感应形成电
信号,例如
电压。使用若干发射线圈并重复这些测量,得到多组测量数据并将其用于使关注对象的电磁属性随时间的变化
可视化。
[0004] MIT对所有这三种无源电磁属性(电导率、电容率和磁导率)是敏感的。结果,例如,关注对象中的电导率贡献可以被重构。具体而言,MIT适合用于针对生物组织重构图像,因为这种组织的磁导率值μR≈1。
[0005] 由涡流电流感应形成的次磁场带有关于待测量对象的信息。然而,由次磁场感应形成的电压∆V非常小,并且电压∆V与在测量线圈上测量的电压V的比率,例如|∆V/V|在-7一些线圈上可以小至10 。这在当前MIT系统中引入了至少两个问题:第一,系统需要高度精确的ADC,这增加了
硬件成本,且其次,系统对噪声非常敏感且因而对其探测性能产生限制。
[0006] Hermann等人在Physiol. Meas, Vol. 26, pp. S307-S318, 2005中的“A new type of gradiometer for the receiving circuit of magnetic induction tomography (MIT)”这一
现有技术文件公开了一种通过使用用于接收
电路的梯度计来减去一对差分线圈中的信号的方法,这提高了MIT系统的灵敏度。
[0007] 然而,梯度计线圈对线圈装置中的几何对称性非常敏感。当线圈对偏离理想对称形状时,例如,由于机械和/或
温度不稳定造成的
变形的原因,线圈不能完美地相互补偿。
发明内容
[0008] 实现一种具有提高的系统灵敏度的图像重构系统,这将是有利的。
[0009] 为了更好地解决这些事宜中的一个或多个,在本发明的第一方面中,提供了一种用于重构关注对象的图像的系统,且该系统包含:- 至少一个发射线圈,其配置成生成将应用到关注对象的主磁场;以及
- 至少一个测量线圈装置,其配置成测量由次磁场感应形成的
电信号,该次磁场是由关注对象响应于主磁场而生成的;
其中所述至少一个测量线圈装置包含置于基本上同一平面内的多个测量线圈。
[0010] 通过使用置于一个平面内的多个独立测量线圈从而替代传统单个测量线圈,可以独立地计算由关注对象的电导率分布变化造成的次磁场变化,使得当其在例如MIT系统中使用时提高灵敏度。
[0011] 有利的是,该系统还包含处理器,其配置成基于由次磁场感应形成的测量的电信号重构所述关注对象的图像,该处理器具有用于控制所述多个测量线圈中的每一个以测量其上的第一和第二电信号的控制单元。
[0012] 在一个
实施例中,该第一和第二电信号为感应电压,以及该处理器还包含第一选择单元,其配置成从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有所述线圈两端的该第一电压和第二电压之间的差分电压与该第一电压之间的比率的最大绝对值。
[0013] 在另一实施例中,该第一和第二电信号为感应电压,且该处理器还包含第二选择单元,其配置成从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有所述线圈两端的该第一电压和第二电压之间的差分电压的最大绝对值。
[0014] 有利的是,该处理器还包含第一计算器,其配置成基于对应于所选择的测量线圈的差分电压计算该关注对象的电导率分布变化。
[0015] 通过选择对次磁场变化最敏感的测量线圈以及通过在图像重构中使用对应于所选择的测量线圈的差分电压,本发明提高了MIT系统的灵敏度。
[0016] 另外有利的是,该处理器还包含第二计算器,其配置成基于从所述多个第一和第二电压导出的多个加权差分电压计算该关注对象的电导率分布变化。
[0017] 通过对由所述多个测量线圈生成的差分电压进行加权,更多的独立测量数据可以用于图像重构,导致MIT系统的灵敏度增大。
[0018] 依照本发明另一方面,提供了一种重构关注对象的图像的方法,所述方法包含下述步骤:(a)通过至少一个发射线圈生成将应用到关注对象的主磁场;以及
(b)通过至少一个测量线圈装置测量由次磁场感应形成的电信号,所述至少一个测量线圈装置包含置于基本上同一平面内的多个测量线圈,该次磁场是由关注对象响应于主磁场而生成的。
[0019] 下文给出本发明的详细解释以及其它方面。
附图说明
[0020] 根据参照附图的下述详细描述,本发明的上述和其它目的及特征将变得更见清楚,在附图中:图1为传统测量线圈装置的示意图;
图2为依照本发明的测量线圈装置的实施例的示意图;
图3为依照本发明的系统的实施例的示意图;
图4为依照本发明的方法的示意性
流程图。
[0021] 在所有的附图中,相同的附图标记用于表示相似的部分。
具体实施方式
[0022] 在MIT系统中,对象电导率分布和在线圈上的测量的电压之间的关系由麦克斯韦方程建立。测量线圈两端的电压为
电场(E)在该线圈上的积分:(1)。
[0023] 在麦克斯韦方程中,电导率和电场之间的关系按下述导出:(2)。
[0024] 在该方程中, 为在发射线圈上应用的电流
密度,以及且 (3)
其中σ、μ和ε分别为电导率、电容率和磁导率,以及 和 分别为自由空间的磁导率和电容率并且 。空间中任意点处的电场无法直接测量,反而是测量电压V,该电压为电场沿着该线圈的积分。
[0025] 图1为传统MIT系统中使用的包含单个测量线圈100的测量线圈装置的示意图。由于电场在该空间中为3D向量场,以及如果线圈划分为四个段,这四个段上的测量结果是不同的。由于相同的原因,这四个段上的对应于电导率分布变化的电压差∆V将是不同的。
线圈的四个段其中之一上的电压差将对关注对象的电导率分布变化最为敏感。
[0026] 然而,测量每个线圈段两端的电压是困难的。这意味着该线圈只能测量通过由该线圈围成的面积的由次磁场变化感应形成的电压,并且无法识别哪个线圈段对由关注对象的电导率变化造成的次磁场变化最敏感。
[0027] 基于对线圈两端的测量的电压和电导率分布变化之间的关系的理解和认知,本发明提供了一种系统,该系统包含一种通过使用置于基本上同一平面内的多个独立测量线圈从而替代传统单个测量线圈而对电导率分布变化最敏感的新颖测量线圈装置以用于测量电压。
[0028] 图2为依照本发明的测量线圈装置200的实施例的示意图。
[0029] 测量线圈装置包含置于基本上同一平面内的多个独立线圈。在此实施例中,测量线圈装置包含四个独立线圈201、202、203和204。
[0030] 线圈可具有不同形状,例如,它们可以是扇形或方形。相应地,由线圈围成的面积可以基本上相同或不同从而调适到不同的应用。
[0031] 有利的是,线圈印刷在印刷
电路板(PCB)上且
采样通道用于从每个测量线圈读取测量数据。
[0032] 通过使用置于一个平面内的多个独立测量线圈从而替代传统单个测量线圈,可以独立地计算由关注对象的电导率分布变化造成的次磁场变化。
[0033] 图3为依照本发明的系统300的实施例的示意图。
[0034] 系统300包含配置成生成主磁场的发射线圈312、314。主磁场在关注对象301中感应形成涡流电流。关注对象301可以是人的头部或者是一
块导电材料。例如,发射线圈312、314被供应交变电流从而生成主磁场。
[0035] 系统300还包含配置成测量由次磁场感应形成的电信号的至少一个测量线圈装置315、317,该次磁场是由关注对象响应于主磁场而生成的。具体而言,次磁场由关注对象中的涡流电流生成,该涡流电流是由主磁场感应形成的。每个测量线圈装置包含置于一个平面内的多个测量线圈,如图2所示。发射线圈312、314和测量线圈装置315、317可以布置在托架303上。
[0036] 系统300还包含配置成基于测量的电信号(例如线圈两端感应形成的电压)重构图像的处理器320。将参考图2和图3在下文中解释处理器320的详细实施方式。
[0037] 如图3所示,该处理器包含控制单元322,其配置成控制所述多个测量线圈中的每一个,从而在关注对象的电导率分布变化之前和之后测量所述测量线圈上的第一和第二电信号。
[0038] 这意味着每个线圈测量第一电压和第二电信号,例如在测量线圈上感应形成的电压。该第一电压和第二电压之间的差分电压是由关注对象的电导率分布变化引起的。
[0039] 对于如图2所示测量线圈具有四个测量线圈的情形,存在四组测量结果: 和, =1,2,3,4, 和 表示在电导率分布变化之前和之后测量的电压。通常 和为由主磁场和次磁场感应形成的电压的相应总和。
[0040] 应注意,取决于用于收集测量数据的硬件的配置,测量控制可以是顺序的,即每个线圈一个接一个地测量,或者是并行的,即所有线圈同时测量。
[0041] 存在不同的方式来选择/识别其上的感应电压对电导率分布变化最敏感的测量线圈。对于每次测量,测量的电压对于电导率分布变化的灵敏度可以定义为:(4)
其中 表示对应于电导率分布变化的电压变化。
[0042] 由于由次磁场变化生成的电压 与测量的电压 或 相比非常小,灵敏度-7可以小至10 。当每个线圈围成的形状和面积基本相同时,对于不同线圈, 之间的差异可以是非常有限的,且在这种情形中, 可以用于指示灵敏度。
[0043] 在一个实施例中,该处理器还包含第一选择单元324,其配置成从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有所述线圈两端的第一和第二电压之间的差分电压与该第一电压之间的比率的最大绝对值。这意味着具有 的最大绝对值的测量线圈被选择。
[0044] 在另一实施例中,该处理器还包含第二选择单元325,其配置成从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有所述线圈两端的该第一电压和第二电压之间的差分电压的最大绝对值。这意味着具有 的最大绝对值的测量线圈被选择。
[0045] 在又一实施例中,该处理器还包含第一计算器326,其配置成基于对应于所选择的测量线圈的差分电压 ,计算关注对象的电导率分布变化。
[0046] 关注对象的电导率分布变化的计算可以遵从已知的图像重构理论,例如在现有技术文件“Image reconstruction approaches for Philips magnetic induction tomograph”, M. Vauhkonen, M. Hamsch和C.H. Igney, ICEBI 2007, IFMBE Proceedings17, pp. 468-471, 2007中描述的电导率计算和图像重构的方法。
[0047] 这组参数可以依照下述方程来计算,例如所提到的现有技术中的方程(8):(5)
其中 为加权矩阵, 为正则参数且 为正则矩阵, 为复雅可比矩阵的
虚部,为 和 之间的差分电压且对应于电导率分布变化。
[0048] 在另一实施例中,该处理器还包含第二计算器328,其配置成基于从所述多个第一和第二电压导出的多个加权差分电压,计算关注对象的电导率分布变化。在计算电导率分布变化时, 用于取代方程(5)中的 ,其中M为测量线圈装置中测量线圈的数目,例如,当使用图2所示的测量线圈装置时M=4。以此方式,由所述多个测量线圈获得的全部测量结果对电导率分布变化的计算有贡献。加权参数 可以利用该方法来计算。
[0049] (6)雅可比矩阵也用 加权从而得到用于重构的新的雅可比矩阵。
[0050] 图4为依照本发明的方法的示意性流程图。
[0051] 依照本发明,该方法包含步骤410:使用至少一个或多个发射线圈312、314生成主磁场,该主磁场在关注对象301中感应形成涡流电流。
[0052] 该方法还包含步骤420:通过使用至少一个测量线圈装置315、317,测量由涡流电流生成的次磁场感应形成的信号用于图像重构,所述至少一个测量线圈装置包含置于平面200内的多个测量线圈201、202、203、204。
[0053] 有利的是,该方法还包含步骤430:控制所述多个测量线圈中的每一个,从而在关注对象的电导率分布变化之前和之后测量第一和第二电压 、 以用于图像重构。
[0054] 在一个实施例中,该方法还包含步骤440:从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有该第一电压和第二电压之间的差分电压与该第一电压之间的比率的最大绝对值。
[0055] 在另一实施例中,该方法还包含步骤440':从所述多个测量线圈选择一个测量线圈,所选择的测量线圈具有该第一电压和第二电压之间的差分电压的最大绝对值。可以执行步骤440'以替代步骤440。
[0056] 该方法还包含步骤450:使用方程(5),基于对应于所选择的测量线圈的差分电压,计算关注对象的电导率分布变化。
[0057] 在另一实施例中,该方法还包含步骤450':基于从所述多个第一和第二电压导出的多个加权差分电压,计算关注对象的电导率分布变化。这种情况下, 用于取代方程(5)中的 用以计算电导率分布变化,即执行步骤450'从而替代步骤450,且该方法直接从步骤430转到步骤450'。
[0058] 应注意,测量线圈的选择和/或电导率分布变化的计算可以有利地由
计算机程序,和/或与硬件和
软件组合来实施。
[0059] 还应注意,上述实施例说明而非限制本发明且本领域技术人员将能够设计若干可替换实施例而不背离所附
权利要求的范围。在权利要求中,置于括号之间的任何附图标记不应解读为限制权利要求。动词“包含”及其变型的使用不排除在权利要求中或在
说明书中列出的元件或步骤之外的元件或步骤的存在。在元件之前使用不定冠词“一”或“一个”不排除存在多个这种元件。本发明可以借助包含若干不同元件的硬件以及适当编程的计算机来实施。在罗列若干装置的系统权利要求中,这些装置的若干个可以由同一项的硬件或软件来实施。措词第一、第二和第三等的使用不表示任何顺序。这些措词被理解为名称。