专利名称：X射线相衬成像的方法和系统的制作方法
技术领域：
本发明涉及X射线无损成像领域，特别涉及一种X射线相衬成像的方法和系统。
背景技术：
X射线相衬成像(XPCI)就是利用X射线透过样品后携带的位相信息对样品内部结构成清晰象，这是近几年发展起来的一种新的X射线成像方法。传统的X射线成像系统都是基于样品不同部位对X射线吸收不同进行的，即所谓的吸收反衬成像。从1895年伦琴发现X射线开始，该方法在众多领域特别是临床医学领域获得了巨大的成功。然而，受原理的限制它很难对密度相差不大的软组织如血管、肿瘤等成清晰象。XPCI技术原则上可分辨密度变化0.0003～0.002g/cm3的不同组织的边界，空间分辨率可达1μm。而目前被认为是软组织检查最好方法的传统CT可分辨的密度变化为0.01g/cm3，空间分辨率仅为亚毫米量级。
较传统的X射线成像方法而言，XPCI方法的优势是明显的。它在生物软组织内部结构研究中有其独特的优越性，其中包括1、血管系统。不用注射造影剂可清晰地观察到动、静脉血管，如心血管系统、毛细血管微循环等；2、呼吸系统，包括气管、支气管和肺的精细结构；3、神经系统，包括传出神经末梢的精细结构；4、包括肝、胃、肠等在内的其他软组织。因此，该方法在生物、医学研究中有重要应用，并很有可能发展成为一种新的临床诊断方法。此外，XPCI在材料科学特别是低Z元素材料研究中有重要应用，在化学反应动力学研究中的应用也已取得突破性的进展。
XPCI对光源性能有较高的要求，具备较好的空间相干性是实现高分辨率成像的必要条件。迄今为止，绝大部分研究都是基于同步辐射光源进行的。由于同步辐射装置是一个多学科的大型研究平台，可供某一项研究的机时有限。这一点在国内显得尤为突出，目前国内仅有北京同步辐射装置能提供XPCI研究的硬X射线。由于它是第一代的兼用机，设计时并未优化到同步辐射，其X射线输出特性无法与国际上普遍采用的第三代同步辐射装置相比，而且每年能提供给同步辐射用户的机时也十分有限。因此，发展一套XPCI实验室成像系统是很有必要的。
XPCI对X射线源的空间相干性要求较高，而对时间相干性则要求不高。S.W.Wilkins等人首次采用微聚焦X射线管的多色辐射实现了相衬成像，这说明建立一套XPCI实验室成像系统是可行的。但他们的微聚焦管X射线管产生的X射线光斑尺寸在10～20微米之间，这个直接影响最终的成像分辨率，同时光源和样品之间的距离必须比较远(30厘米以上)，才能满足成像的相干性要求，带来的问题是X射线光子通亮急剧下降，曝光时间延长(他们使用X射线胶片)，因此干扰因素如震动、温度湿度变化等使得成像分辨率变差，从而严重影响了该类相衬成像系统的应用。
发明内容
本发明的目的是解决上述课题，提供一套实验室使用的新型X射线位相衬度(相衬)成像系统和方法，用于生物软组织、低Z材料样品内部结构的无损成像研究。
本发明的设计原理为
如图1所示，在常规的依次包括产生硬X射线光源的X射线管1、放置样品2的扫描台3和探测器4的X射线成像系统基础上，为实现XPCI，空间相干性是必要条件，也是光源选择时首先要考虑的问题。理想的平行光或点光源是完全空间相干的，对于实际光源只要它的平行度足够好或是点源尺寸足够小，它就具备了一定的空间相干性。同步辐射输出的X射线束呈一个很窄的锥状，其张角与电子的能量成反比，是一个非常小的值。在微弧度量级的张角内辐射出高强度的X射线光子，具有良好的空间相干性和相当高的相干通量，是相衬成像的理想光源。而普通实验室X射线光源的性能则相距甚远，因此选择合适的光源是保证系统成功的关键。
X射线管是最简单也是最通用的实验室X射线光源。普通的X射线管电子束在靶面上的投影呈线状，尺寸一般为1×10mm2左右，这也就是X射线光源的实际尺寸。尽管通过选择不同的引出角度可得到不同的有效光斑，如得到尺度为1mm的方形或圆形光斑。但很难进一步减小光源点尺寸，因为这就意味着减少输出光子数目，对成像来说还减小了视场。因此，普通X射线管并不是XPCI实验室成像系统的理想光源。
目前的技术，X射线管的光源点尺寸已达亚微米，又称为纳聚焦管。微聚焦管通常只采用一个电磁透镜来聚焦电子束至靶，纳聚焦管则在准直单元附近增加了一个电磁聚焦镜，其作用是在电子束到达物镜前进行预聚焦。另一不同之处就是在靠近靶子的位置增加一个光阑过滤杂散电子。采用纳聚焦X射线管可实现更高的成像分辨率，对形成位相衬度而言还可以获得更高的有效通量或相干通量。
但还需合适的成像方式。除X射线全息术以外，实现XPCI的途径主要有1、X射线干涉。利用晶体对X射线的劳厄衍射和透射分束实现两束X光的干涉来探测样品的位相变化，其局限性a)视场小，很难用于临床；b)稳定性要求很高，达亚纳弧度量级，易产生赝象；c)仅适于测量位相变化量小且缓变的样品。2、X射线衍射增强。利用晶体布喇格衍射摇摆曲线随角度变化很灵敏来分辨很小的位相变化梯度，相当于一个带通滤波器。该方法得到的衬度与工作点的选择有关，极易产生赝象，得到的图象较难解释。3、同轴轮廓成像。成像光路与同轴全息一致，类同通常的医学透射成像，具有光路简单、视场大以及图象直观无赝象等特点。在上述各种XPCI成像方法中，全息、干涉和衍射增强方法除要求X射线源具有好的空间相干性外，对光源的时间相干性或单色性也有较高的要求，只有同轴轮廓法允许多色X射线成像。因此，要建立一套XPCI实验室成像系统，同轴轮廓成像方式是唯一的选择。
同轴轮廓成像的衬度形成机理可由菲涅尔衍射来描述，对纯位相物体和点光源的情形，在一阶近似下像面上的光强分布可描述为I(Mx,My;d1+d2,k)=(I0/M2){1-2πrek2d2MΔx,y2∫ρ(x,y,z)dz}]]>其中M＝(d1+d2)/d1为放大倍数，d1，d2分别为样品到光源点和探测平面的距离，k＝2π/λ，λ为波长，re为经典电子半径，ρ为样品内电子密度，I0为物面上的光强值。由此可见，像面上的光强分布正比于电子密度的二阶微分，且与光子能量(k)无关，其中对多色源而言1/k2可用谱分布权重因子替代。这就意味着，无需重构即可直接获得样品内部结构的位相衬度像，而且对光源的单色性要求不高。由于只是在样品内电子密度有突变的位置形成衬度，该方法得到的实际上是一个密度变化的轮廓图，故称为同轴轮廓成像。像面上的强度变化是由X射线通过样品时的位相变化导致的，因而对电子密度变化的探测灵敏度很高。另一方面，只是在密度变化位置附近的X射线相干叠加对衬度的形成有贡献，物平面上只需局域相干照明即可，因此对空间相干性的要求也不高。空间相干长度可表达为Is＝d1λ/σ，其中σ光源尺寸，λ为波长。
由上可见，样品面上入射X射线的空间相干长度ls与d1成正比，增大样品到光源点的距离可改善空间相干性，而d1的增大又会迅速减小放大倍数M，从而提高对探测器空间分辨率的要求，因此，需根据具体实验条件选择一个合适的d1值；而选用的纳聚焦管产生的X射线光斑尺寸可以很小，从而成像分辨率可以很高，但是受到其它因素如探测器等的制约，加之样品与探测器之间的距离d2，与最后的像的分辨率和反衬度相关，故为使得分辨率和反衬度最佳，确定这个最佳距离及选择合适的探测器也是至关重要的。
因此，为实现上述目的，本发明的技术方案为一种X射线相衬成像的方法，其利用依次包括产生硬X射线光源的X射线管，放置样品的扫描台和探测器的成像系统，其特征在于其选择光斑尺寸为0.5-5微米的纳聚焦管来产生X射线光源，所述X射线通过同轴轮廓成像方式将样品成像于该探测器上。
其中，所述的同轴成像方式具体包括以下步骤1)调整该纳聚焦管产生的光源点与扫描台上的样品之间的距离d1，其根据下式得出d1＝lsσ/λ，其中，ls为空间相干长度，其值不少于1微米，σ为光源尺寸，λ为X射线波长；2)调整该样品与探测器之间的距离d2，其可根据下式估计2λd2Mf2≈1，其中，λ是X射线波长，M是放大倍数，M＝(d1+d2)/d1，f是样品需要分辨的空间频率，f＝2π/δ，提样品需要分辨的细节；3)选择探测器进行成像。
所述的距离d1至少为8毫米。
而该距离d2至少为5毫米，一般选取的范围为5毫米至2米之间。
其中，该选择的探测器的分辨率最好符合下式要求分辨率≤δ×M。
另外，本发明的一种X射线相衬成像的系统，依次包括产生硬X射线光源的X射线管，放置样品的扫描台和探测器，其特征在于该X射线管为光斑尺寸为0.5-5微米的纳聚焦管。
其中，该纳聚焦管选择的电压为40kV-200kV，电流强度不低于100μA，出射光锥张角为10-45度。
该探测器的分辨率符合下式要求分辨率≤δ×M，其中δ为样品需要分辨的细节尺寸，M为需放大的倍数。
该探测器可为X射线工业底片或数字探测器。
而所说的数字探测器是指X射线CCD(电荷耦合组件)或成像板。
本发明的积极进步效果在于本发明使用纳聚焦管作为X射线光源，在保证足够的空间相干长度的前提下，其调整样品到光源点的距离最小可为d1＝8mm，样品尺寸一定，入射通量与d1的平方成反比。与已有的S.W.Wilkins的系统相比(d1＝300mm)，仅考虑这一因素时有效通量要增大1400倍。综上述，本发明的方法和系统在空间分辨率和有效通量上已有相当大的提高，可以满足一般实验室研究的需要，并且具备了利用数字探测器开展实时相衬成像的潜力。
本发明的方法和系统的优点在于(1)采用光源点尺寸可达500nm的纳聚焦X射线管，大大提高了X射线相衬成像的有效通量，从而当选用X射线工业底片时可以降低曝光时间，最短可至30秒。当然，具体视底片类型、样品的吸收、d1和d2的距离而定，可从30秒到2小时不等；并且可提高信噪比，为该类成像系统的实用化打下了结实的基础；(2)空间分辨率可达500nm亚微米量级，是生物软组织内部结构高分辨率无损成像的有效手段；(3)采用X射线CCD探测器时，还可以实现实时X射线相衬成像；(4)X射线管输出锥形光束，系统本身具备放大作用，从而降低对探测器分辨率的要求。通过样品台调节样品与光源点的距离，以适应不同种类、不同尺寸的样品成像，适合于大尺寸样品成像。
从而可使普通实验室进行生物软组织、低Z材料样品内部结构的无损成像研究。
图1为本发明X射线相衬成像系统的配置示意图。
图2为应用本发明的X射线相衬成像系统和方法摄制的气垫薄膜中单个气泡的X射线相衬像。
图3为应用本发明摄制的蚊子的X射线相衬像，其中右图为蚊子腹部管道的显微图。
图4(a)为应用本发明一实施例摄制的新鲜蜻蜓的位相衬度图；(b)应用传统的吸收成像原理摄制的新鲜蜻蜓的吸收衬度图。
图5(a)为应用本发明一实施例摄制的新鲜蜻蜓样品眼部的X射线相衬成像的放大图；(b)应用传统的吸收成像原理摄制的新鲜蜻蜓样品眼部的X射线吸收成像的放大图。
图6为应用本发明摄制的小檗碱沉淀的位相衬度成像。
图7为应用本发明摄制的秋水仙碱沉淀的位相衬度成像。
具体实施例方式
为更好地说明本发明系统的配置及本发明方法的实验操作步骤，以及其功能，下面举几个实施例来说明之。其中样品包括气泡、蚊子、蜻蜓等，所有样品都是新鲜、未经任何处理的；选用纳聚焦管(德国Pheonix公司phoenix|x-ray package|analyser 160NF)的光源点尺寸为500nm，并按常规根据样品选择其电压、电流强度及出射光锥张角，下列实施例的出射光锥张角约为30°；采用高分辨率X光胶片(富士IX80或50系列，也可选用其它规格)作为探测器，并按常规视底片类型、样品的吸收、d1和d2的距离而定；扫描台为常规的可以X、Y、Z轴三维平动，还可以转动和扭动扫描台。
操作步骤打开纳聚焦管X射线光源，根据不同的样品选择合适的电压和电流，经过一段时间后X射线趋于稳定，可以用于实验。选择合适的样品到光源的距离d1，通常情况下大于8毫米、尤其是1厘米以上就可以得到足够的空间相干性；接着调整样品与底片的距离d2，对于吸收非常小的样品，即该样品的吸收厚度ts远远大于位相厚度tp(其中ts＝μ-1，μ是样品的吸收率，tp＝2π/φ，φ是单位长度的位相差)时，可直接根据上述公式估算，而对于吸收不是非常小的样品，则根据上述公式估算一个距离后，再作些调整，从而找到一个最佳距离使得成像质量最佳。
实施例1图2给出了气垫薄膜中单个气泡的X射线相衬像，实验中工作电压70kV，电流200μA，曝光时间100s，d1＝50mm，d2＝450mm。在该实验条件下，有机薄膜对X射线的吸收可以忽略，可以看成纯位相物体。实验样品是一种常用的合成缓冲包装材料，在两层塑料薄膜之间采用特殊的方法封入空气，使薄膜之间连续均匀地形成气泡。根据缓冲要求不同，也可制成三层的气垫薄膜，它的缓冲效果比两层更佳。图中给出的是薄膜交叠处的边界，可以看出，该样品应该是缓冲效果较好的三层气垫薄膜。边界四周没有形成旁通，说明它的质量是合格的。
实施例2图3给出了蚊子的X射线相衬像，实验中工作电压70kV，电流100μA，曝光时间120s，d1为130mm，d2为400mm。从蚊子腹部管道的显微图可清晰地分辨出主管道、次级管道、毛细管及其分叉，分辨率可达几个微米。
实施例3为了和吸收衬度像比较，选取蜻蜓样品分别记录了它的位相衬度像和吸收衬度像。图4和图5给出了蜻蜓样品的实验结果。图4(a)是我们利用本发明纳聚焦管硬X射线光源拍摄的位相衬度像(d1＝150mm，d2＝380mm)，除了吸收衬度外，还包含了明显的位相衬度像，图4(b)则主要是吸收衬度像(d1＝150mm，d2＝5mm)。可以非常明显的看出，相比图4(b)，图4(a)显示了蜻蜓的更为明显的细微结构。根据计算可以得到在该实验条件下，当蜻蜓的空间周期为9μm-1时，可以得到比较好的反衬度，这与实验结果是比较一致的，见图4(a)和图5(a)。图5(a)是利用本发明纳聚焦管硬X射线光源拍摄的蜻蜓单个眼睛的位相衬度像(d1＝150mm，d2＝380mm)，包含了明显的位相衬度像。作为比较，图5(b)为相应的吸收衬度像(d1＝150mm，d2＝5mm)。实验中X射线管的工作电压70kV、电流150μA，曝光时间分别为15分钟和2分钟。由X射线衍射的光学传递函数(OTF)随空间频率变化曲线可知，探测器到样品的距离d2值足够小时，记录的实际上只有吸收衬度。另外从现有技术可知，X射线管辐射连续谱强度最大的位置一般在最短波长的1.5倍处，70kV工作电压对应的波长应为0.0266nm，样品到光源点的距离d1＝150mm，样品面上入射X射线空间相干长度为8μm，这说明本实施例实验中X射线光束已具有较高的空间相干性，从图5(a)的位相衬度像中可观察到蜻蜓眼睛内部结构更多的细节，其中的网状结构可能就是蜻蜓的复眼。
实施例4
图6显示小檗碱沉淀的位相衬度成像。选用的靶电压90kV，流强150μA，光源到样品的距离d1为130mm，样品到底片的距离d2为250mm，曝光时间1800秒。从图中可以看到小檗碱沉淀的花瓣状结构。
实施例5图7显示秋水仙碱沉淀的位相衬度成像。其中选用的靶电压45kV，流强150μA，光源到样品的距离d1为130mm，样品到底片的距离d2为200mm，曝光时间1500秒。从图中可以看到针状结构清晰可辨，这样的结果用现有的X射线成像方法和系统尚无法观察到。
显然，如果上述实施例中的探测器如换用X射线CCD，则可以进行实时成像。而一些成像时按常规根据不同的样品来调整的不同参数，如纳聚焦管采用的其它光源尺寸、电压、电流强度、出射光锥张角等在此不一一举例说明。
权利要求
1.一种X射线相衬成像的方法，其利用依次包括产生硬X射线光源的X射线管，放置样品的扫描台和探测器的成像系统，其特征在于其选择光斑尺寸为0.5-5微米的纳聚焦管来产生X射线光源，所述X射线通过同轴轮廓成像方式将样品成像于该探测器上。
2.根据权利要求1所述的方法，其特征在于所述的同轴成像方式具体包括以下步骤1)根据样品，调整该纳聚焦管产生的光源点与扫描台上的样品之间的距离d1，其根据下式得出d1＝lsσ/λ，其中，ls为空间相干长度，其不少于1微米，σ为光源尺寸，λ为X射线波长；2)调整该样品与探测器之间的距离d2，其可根据下式估计2λd2Mf2≈1，其中，λ是X射线波长，M是放大倍数，M＝(d1+d2)/d1，f是样品需要分辨的空间频率，f＝2π/δ，δ是样品需要分辨的细节。3)选择探测器进行成像。
3.根据权利要求2所述的方法，其特征在于该距离d1至少为8毫米。
4.根据权利要求2所述的方法，其特征在于该距离d2至少为5毫米。
5.根据权利要求2所述的方法，其特征在于该选择的探测器的分辨率符合下式要求分辨率≤δ×M。
6.一种X射线相衬成像的系统，依次包括产生硬X射线光源的X射线管，放置样品的扫描台和探测器，其特征在于该X射线管为光斑尺寸为0.5-5微米的纳聚焦管。
7.根据权利要求6所述的X射线相衬成像的系统，其特征在于该纳聚焦管选择的电压为40kV-200kV，电流强度不低于100μA，出射光锥张角为10-45度。
8.根据权利要求6所述的X射线相衬成像的系统，其特征在于该探测器的分辨率符合下式要求分辨率≤δ×M，其中δ为样品需要分辨的细节尺寸，M为需放大的倍数。
9.根据权利要求8所述的X射线相衬成像的系统，其特征在于该探测器为X射线工业底片或数字探测器。
10.根据权利要求9所述的X射线相衬成像的系统，其特征在于该数字探测器为X射线CCD或成像板。
全文摘要
本发明提供一种X射线相衬成像的系统，依次包括产生硬X射线光源的X射线管，放置样品的扫描台和探测器，其特征在于该X射线管为光斑尺寸为0.5－5微米的纳聚焦管。本发明还提供一种X射线相衬成像的方法，其利用上述成像系统通过同轴轮廓成像方式将样品成像于该探测器上，其主要包括以下步骤1)根据样品，调整该纳聚焦管产生的光源点与扫描台上的样品之间的距离d
文档编号G01N23/04GK1725000SQ20041005301
公开日2006年1月25日 申请日期2004年7月21日 优先权日2004年7月21日
发明者肖体乔, 陈敏, 徐洪杰, 杜国浩 申请人:中国科学院上海应用物理研究所