专利名称：电阻抗断层功能图像与解剖图像的合成方法
技术领域：
本发明属于生物医学成像技术领域，尤其是一种电阻抗断层功能图像与解剖图像 的合成方法。
背景技术：
医学成像已经成为现代医疗领域不可或缺的一部分，其应用贯穿于整个临床工 作，不仅广泛用于疾病诊断，而且在外科手术和放射治疗等的计划设计、方案实施以及疗效 评估方面发挥着重要作用。目前，医学图像分为解剖图像和功能图像两个部分。解剖图像 主要描述人体形态信息，包括X射线透射成像、CT、MRI、US以及各类内窥镜(如腹腔镜及喉 镜)获取的序列图像等；功能图像主要描述人体代谢信息，包括PET、SPECT、fMRI等。电阻 抗断层成像(EIT)技术是近三十年新兴的无创功能成像技术，它具有无核素、无辐射、体积 小、价格低及可连续成像等特点，是一种新的功能成像技术。电阻抗断层成像与解剖成像是两种不同的成像模式，前者提供的是功能分布信 息，后者提供的是解剖结构分布信息，不同成像模式提供的图像信息是具有互补性的，为了 综合使用多种成像模式以提供更全面的信息，有必要将两种类型图像的有效信息进行整 合。目前，功能图像与解剖图像进行合成的仅有PET与CT或MRI的图像合成，而没有将电 阻抗断层功能图像与解剖图像合成的方法。
发明内容
本发明的目的在于克服现有技术的不足，提供一种电阻抗断层功能图像与解剖图 像的合成方法，该方法解决了电阻抗断层功能图像信息与临床解剖图像信息的有效合成问 题，为临床诊断和治疗提供更加充分、全面的解剖和功能图像信息。本发明解决现有的技术问题是采取以下技术方案实现的一种电阻抗断层功能图像与解剖图像的合成方法，包括如下步骤(1)获取相同时刻的电阻抗断层功能图像和解剖图像；(2)将电阻抗断层功能图像和解剖图像进行边界配准；(3)提取电阻抗断层功能图像中的电阻抗敏感区域的图像特征；(4)将电阻抗断层功能图像与解剖图像使用图像分割方法进行图像合成得到伪彩 色显示的合成图像。而且，在步骤(1)和(3)之间还包括采用共轭梯度方法重构电阻抗断层功能图像 的步骤。而且，所述重构电阻抗断层功能图像的具体步骤(1)采用canny算子提取解剖图像的边界特征；(2)利用解剖图像的边界特征建立电阻抗断层成像数学模型进行正问题求解；(3)采用共轭梯度法重建具有与解剖图像相同边界特征的电阻抗断层功能图像。而且，所述步骤(4)中的图像分割方法是将电阻抗敏感区域的图像特征用颜色映射到解剖图像中的对应位置，从而得到伪彩色显示的合成图像。本发明的优点和积极效果是1、本方法有效地将电阻抗断层功能图像信息与解剖图像信息合成生成伪彩色显 示图像，充分发挥了电阻抗断层图像及高分辨率的解剖图像各自的特点，将两种图像进行 优势互补、相互融合，既能反映解剖结构又能反映测量区域电导率变化情况，为临床诊断和 治疗提供更加充分、全面的解剖和功能图像信息。2、本方法解决了电阻抗断层功能图像与解剖图像的有效合成问题，其与目前的 PET/CT技术相比具有处理方法简单、易于实现以及可连续多次成像的优势。3、本发明将电阻抗断层功能图像与解剖图像进行合成生成伪彩色显示图像，解决 了电阻抗断层功能图像信息与临床解剖图像信息的有效融合问题，为临床诊断和治疗提供 更加充分、全面的解剖和功能图像信息，具有处理方法简单、性能高且易于实现等特点。
具体实施例方式以下结合附图
对本发明实施例做进一步详述。一种电阻抗断层功能图像与解剖图像的合成方法，是将电阻抗断层功能图像与解 剖图像信息合成在一起从而获得伪彩色显示的合成图像，为临床诊断和治疗提供一种更加 充分、全面的解剖和功能图像信息。电阻抗断层功能图像(fEIT功能图像)和解剖图像(CT 解剖图像)是采用不同的成像系统在同一时刻获取的图像。一种电阻抗断层功能图像与解剖图像的合成方法，包括如下步骤(1)获取相同时刻的电阻抗断层功能图像和解剖图像；(2)采用共轭梯度方法重构电阻抗断层功能图像；重构电阻抗断层功能图像的目的在于提高电阻抗断层功能成像质量，共轭梯度方 法的具体步骤为①采用canny算子提取解剖图像的边界特征；②利用解剖图像的边界特征建立电阻抗断层成像数学模型进行正问题求解；电阻抗成像电阻抗断层成像数学模型可以用如下线性系统方程Ax = b其中，A为Jacobian矩阵(灵敏度系数矩阵)；b为边界测量数据；x为内部场域 的电导率分布。考虑到边界特征和接触阻抗，数学模型描述为
v> = 0，场域内
J>#ds = I,注入区
E+ dn
炉+2^￡ = 11,测量区 on
式中，n为边界上任意一点的外法向量，识为场域内的电势分布，、为电导率，s为 接触面积，z为接触阻抗，I为注入电流值，U为边界测量电压值。该问题是在给定第三类边 界条件下求电势识的泊松方程的定解。解剖图像边界划分得到的不同区域具有不同的电导 率，所以求得的解包含了解剖信息，而且该解与边界电流之间的线性映射关系就是第一式 c^ Jacobian
③采用共轭梯度法重建具有与解剖图像相同边界特征的电阻抗断层功能图像；电阻抗断层功能图像重建就是已知Jacobian矩阵A和测量数据b，求解X。由于 A不是对称正定矩阵，因此，需要对原方程进行转换。AtAx = ATb由于A并非列满秩矩阵，因此，ATA也不是对称正定矩阵，可以通过增加正则化矩阵 的办法来解决这个问题即可得到(AtA+ A I) x = ATb此时，系数矩阵为对称正定矩阵，可以运用共轭梯度算法进行求解。共轭梯度法是依次沿着非零共轭方向进行搜索达到其极小点而进行求解的，当迭 代误差小于设定阈值，那么迭代结束。迭代过程经过第i次迭代得到的解Xi连续逼近真解，迭代误差巧=b_AXi，搜索方向设定为 Pi，迭代次数与迭代误差的更新为1 = APi，从而有Xi = Xi_!+a jpj,、是标量，使得HcOTiHa)最小，其中r(a)=ri_1-aqi_1。所以有
权利要求
一种电阻抗断层功能图像与解剖图像的合成方法，其特征在于包括如下步骤(1)获取相同时刻的电阻抗断层功能图像和解剖图像；(2)将电阻抗断层功能图像和解剖图像进行边界配准；(3)提取电阻抗断层功能图像中的电阻抗敏感区域的图像特征；(4)将电阻抗断层功能图像与解剖图像使用图像分割方法进行图像合成得到伪彩色显示的合成图像。
2.根据权利要求1所述的电阻抗断层功能图像与解剖图像的合成方法，其特征在于 在步骤(1)和(3)之间还包括采用共轭梯度方法重构电阻抗断层功能图像的步骤。
3.根据权利要求2所述的电阻抗断层功能图像与解剖图像的合成方法，其特征在于 所述重构电阻抗断层功能图像的具体步骤(1)采用carmy算子提取解剖图像的边界特征；(2)利用解剖图像的边界特征建立电阻抗断层成像数学模型进行正问题求解；(3)采用共轭梯度法重建具有与解剖图像相同边界特征的电阻抗断层功能图像。
4.根据权利要求1所述的电阻抗断层功能图像与解剖图像的合成方法，其特征在于 所述步骤(4)中的图像分割方法是将电阻抗变化区域的特征参数值用颜色映射到解剖图 像中的对应位置，从而得到伪彩色显示的合成图像。
全文摘要
本发明涉及一种电阻抗断层功能图像与解剖图像的合成方法，包括如下步骤(1)获取相同时刻的电阻抗断层功能图像和解剖图像；(2)将电阻抗断层功能图像和解剖图像进行边界配准；(3)提取电阻抗断层功能图像中的电阻抗敏感区域的图像特征；(4)将电阻抗断层功能图像与解剖图像使用图像分割方法进行图像合成得到伪彩色显示的合成图像。本发明将电阻抗断层功能图像与解剖图像进行合成生成伪彩色显示图像，解决了电阻抗断层功能图像信息与临床解剖图像信息的有效合成问题，为临床诊断和治疗提供更加充分、全面的解剖和功能图像信息，具有处理方法简单、性能高且易于实现等特点。
文档编号G06T7/00GK101984462SQ20101027585
公开日2011年3月9日 申请日期2010年9月8日 优先权日2010年9月8日
发明者王化祥, 陈晓艳 申请人:天津科技大学