针对软物质高密度分辨的X射线相位衬度增强方法与流程
本发明属于硬x射线成像与检测领域，具体涉及弱吸收软物质样品的x射线高密度分辨的相位衬度增强方法。
背景技术：
近年来，随着x射线探测器和新型光源技术的发展，特别是20世纪90年代第三代同步辐射（synchrotronradiation，sr）光源的投入使用，把x射线成像带到一个全新的高度。同步辐射装置由电子直线加速器、增强器、储存环、光束线站等部分组成，电子经直线加速器加速后进入增强器进一步加速，当电子达到预定能量后注入到储存环，在储存环内高速运动。在储存环内有弯转磁铁及各种插入件，可以产生磁场使得电子偏转，这时电子束沿其运动的切线方向向外发射出同步辐射光。同步辐射产生的x射线波段宽广，是探测物质微观结构的有力工具，同时同步辐射光源具有常规光源不可比拟的优良性能，如高亮度，高准直性，高偏振性，高相干性和宽频谱范围等优点，已成为当今众多学科基础研究和高技术开发应用研究的新型高质量光源。
基于上海同步辐射光源（电子储存环能量3.5gev，周长为432米）的x射线成像与生物医学应用光束线站的x射线同轴成像系统，具有装置简单、成像直观且分辨率高的特点，如图1所示，根据研究样品的特性，同步辐射装置wiggler产生的高亮度同步辐射光（能量范围8-72.5kev），经过双晶单色器dcm选通后，成为一束高准直单色光入射到样品上，采用高分辨探测器cmos（c11440–22c,hamamatsuphotonics,japan）采集高质量的x射线与样品相互作用后产生的吸收和相移信号，通过样品台（samplestage）的180°的旋转扫描进行投影成像，利用滤波反投影ct重构技术，获取被研究样品的三维微观原位结构信息。
当x射线穿过样品时，它们对x射线的作用不仅因其吸收因子的存在，改变x射线的强度，还由于其相位因子的存在，可以改变x射线的相位。x射线相位相衬成像（x-rayphase-contrastimaging:xpci）原理是利用x射线透过样品后满足菲涅尔衍射传播规律，在像空间的光强分布如图2（a）所示：根据样品距探测器的距离（sdd）不同，在像空间范围内即可以采集到吸收像（sdd=0），也可以采集到包含吸收衬度和相位衬度的相衬像（sdd>0），图2（b）和（c）是sr-basedil-pci方法对生物软组织成像（兔肺）特点比较，可以明显看出差异性和相衬优势。在硬x射线区域内（10-100kev），对于主要由较轻元素组成的软物质、软组织、高分子聚合材料等弱吸收样品而言，其相位因子较其吸收因子要高3-5个数量级，既其相位衬度远远大于其吸收衬度，ilpci虽然是利用样品的相位信息成像，il-pci得到的是样品的边缘增强图像，成像衬度差、密度分辨率不高，严重限制了该技术用于软物质样品研究。如何进一步从il-pci采集到的图像中重构出增强型相位衬度信息并用于定量研究，受到越来越多研究团队的密切关注。
技术实现要素：
本发明目的是提供了一种针对软物质探测的x射线相位衬度增强方法，以提高x射线成像衬度及密度分辨率，满足基于由低原子序数元素（低z元素）构成的软物质定量研究的需求。
为了解决上述技术问题，本发明提供了如下的技术方案：
针对软物质高密度分辨的x射线相位衬度增强方法，包括如下步骤：
（1）采用硬x射线同轴相衬成像法获取低z软物质样品的相衬数据；
（2）按式i对步骤（1）获取的相衬数据进行相衬增强算法运算，在不损失空间分辨的基础上得到高密度分辨的软物质信息的空间分布
其中，
（3）对步骤（2）得到的高相衬信息的空间分布，采用滤波反投影算法进行图像三维特征结构重构。
优选地，在10-100kev的硬x射线区域内，所述低z元素软物质样品的相位因子较吸收因子高3-5个数量级。
本发明方法有效突破传统x射线同轴成像在弱吸收软物质样品的微结构分析方面的局限性，开创性地将x射线同轴成像法与本发明的高密度分辨的相衬增强算法结合，大幅提高软物质样品的特征结构和形态的探测能力，同时具有高信噪比，降低辐照剂量的特点，从而将x射线同轴相衬成像的应用领域拓展到高分子材料、生物医学软组织、中药材等弱吸收物体样品的微米级无损、高分辨三维成像，实现低z软物质样品的高衬度定量研究。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制，在附图中：
图1是同步辐射硬x射线同轴成像系统结构示意图；
图2是硬x射线同轴相位衬度成像的原理示意图；
图3是用于增强算法的成像系统mtf测量；
图4是针对标准高分子聚合物混合样品（软物质）的增强成像实验结果；
图5是针对大鼠肝包虫模型样品（软组织）的增强成像实验结果；
图6是基于相衬增强算法的软组织结构与形态的特征分割和提取结果，包括肝包虫体、血供系统等。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
针对低z元素软物质探测的x射线相位衬度增强方法，包括如下步骤：
（1）采用同步辐射硬x射线同轴相衬成像法获取软物质样品的相衬实验数据；
（2）按式i对步骤（1）获取的相衬数据进行衬度增强算法运算，在不损失空间分辨的基础上得到高密度分辨的样品空间信息分布
其中，
（3）对步骤（2）得到的高衬度信息的空间分布，采用滤波反投影算法进行图像三维重构。
数据采集利用imageproplus6.0软件进行处理。
同步辐射x射线同轴相衬成像采用上海同步辐射光源（电子储存环能量3.5gev，周长为432米）的生物医学成像线站bl13w1的x射线同轴相位衬度成像系统。具体操作：根据低z元素样品的透过率，确定同轴相衬成像的实验参数，包括光子能量、物像距离、探测器像素大小等。设定曝光时间和投影间隔，在[0，π]范围内旋转样品，采集样品同轴相衬成像投影数据；每旋转30°，将样品移出探测器视场，采集2张背景图，数据采集全部结束时，关闭光闸，采集5张暗电流噪声图，以校正投影图的背景。
根据相衬增强算法获得的样品相位信息的空间分布，采用滤波反投影重构算法，进行切片重构、三维堆积、渲染，最终得到软物质样品的高密度分辨三维特征形态。
基于上述方法开展针对软物质探测的x射线高密度分辨的相衬增强方法的测试实验。
一、测量该成像系统的mtf曲线
采用jima分辨率标靶，首先采集图像数据，然后计算出不同空间频率下的mtf函数，并绘制出mtf曲线。
二、实验样品以高分子聚合物混合样品（软物质）为例
高分子聚合物混合样品由pmma棒材一根（直径为1mm，物质密度为1.19g/cm3）、pa6管材三根（外径为1.34mm，内径为1mm，管壁为0.17mm，物质密度为1.13g/cm3）、ps套管一根（内径为3.5mm，物质密度为1.04g/cm3）组成，将pmma棒材和pa6管材放入ps套管中构成，以上各样品材料均由英国goodfellow公司生产。
如图4所示：（a）为硬x射线相衬成像的三维切片重构图（未经过增强算法的重构切片），由于边缘增强效应仅能够分辨实验样品的轮廓信息，对于弱吸收混合样品基本无法体现其不同组分的密度分布和差异，甚至无法与背景分离，导致无法进一步开展定量分析。（b）为经过增强算法的三维切片重构图，不仅有效各组分样品的轮廓信息而且很好地体现了其不同的密度分布，并且与背景也有明显区别，为样品的定量分割和分析奠定了良好地基础。（c）是（a）图的直方图分析，可以看出只有一个灰度分布峰，背景和不同组分材料难以区别。（d）是（b）图的直方图分析，可以看出其灰度分布明显不同，根据其密度不同，从右至左各峰依次为：pmma，pa6和ps，因此，该增强方法可以有效提高软物质密度的分辨能力。（e）是x射线相衬成像的样品三维立体图（未经过增强算法），可以看出仅重构了样品的轮廓三维信息，无法区别其内在物质分布，与真实样品不符。（f）是经过增强算法处理的样品三维立体图，可以明显看出，样品形貌分布和物质密度分布均与真实样品内外表现一致，完全符合样品的真实表现形式。
三、实验样品以大鼠肝包虫模型（软组织）为例
动物致病模型的建立，是用人为的方法，使动物在一定的致病因素（物理的、化学的、生物的）作用下，造成动物软组织、器官或全身一定损害，出现某些类似人类疾病的功能、代谢、形态结构方面的变化或各种疾病，为研究疾病诊断、治疗策略和致病机制研究奠定技术基础。
本发明大鼠肝包虫模型的建立参考“肝囊型包虫病小鼠动物模型的优化”，叶建蔚等，疾病预防控制通报，2012年，第27卷第5期，第1-7页。
图5是采用软物质探测的x射线相位衬度增强方法得到的大鼠肝泡状包虫病理模型样品三维立体重构图（高衬度渲染），可明显观察到的不同三维微观病理形态学表现：棘球蚴区（如箭头1所示）；钙化区（如箭头2所示）；炎症带（如箭头3所示）；原头蚴（如箭头4所示）；以及囊结构和隔膜形态分布（单囊和多囊泡结构：如箭头5所示）。
图6是基于该增强方法获得的高密度分辨的软组织三维定量特征分割和提取结果（高衬度渲染），箭头所指为肝泡状棘球蚴病灶区的三维特征分布，其它部分为钙化区分布；棘球蚴（肝包虫幼体）在血管内外的三维特征分布，证实肝包虫具有通过血供系统侵犯血管及浸润转移的行为方式。说明该增强算法可有效区分并提取不同病灶区域的三维形态与分布特征，并可以开展定量分析。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
- 还没有人留言评论。精彩留言会获得点赞！