一般来说，本发明涉及利用穿透辐射诸如X-射线观察物体的结构特 征。更具体地说，但不是唯一地，本发明涉及对于物体内部边界特征的X- 射线相衬记录，例如成象。

本申请人的国际专利申请WO95/05725(PCT/AU94/00480)和临时专 利申请PN5811/95中公开了适合于利用硬X-射线进行相衬成象的各种配置 和条件。在苏联专利SU1402871和美国专利US-5319694中公开了其它相关 内容。需要在使用硬X-射线进行相衬成象中可以利用相对简单的条件和配 置，这些简单条件和配置至少在某些实施例中与常规的吸收对比辐射成象 方法更为相关。

根据本发明，提供了一种获得物体的边界图象的方法，所说边界表征 折射率的变化，所说方法包括以下步骤：

用具有强横向空间相干性和垂直于所说折射率变化的传播分量的穿透 辐射照射所说边界；和

在一个成象平面接收至少一部分所说辐射，以便构成所说图象，所说 辐射已经被所说边界折射，使得所说边界在所说图象中表现为相应的强度 变化。

本发明还提供一种用于获得一个物体的边界图象的设备，所说边界表 征折射率的变化，所说设备包括：

一个辐射源，其用具有强横向空间相干性和垂直于所说折射率变化的 传播分量的穿透辐射照射所说边界；和

一个探测器，用于接收至少一部分所说辐射，以便构成所说图象，所 说辐射已经被所说边界折射，使得所说边界在所说图象中表现为相应的强 度变化。

本发明还提供一种获得具有急剧折射率变化或由厚度变化限定的内部 边界的相衬记录的方法，该方法包括以下步骤：

利用具有垂直于所说折射率变化方向或者所说厚度变化方向的显著的 传播向量分量，并且对于折射率或厚度的变化具有足够强的横向空间相干 性的穿透辐射辐照所说边界，以在所说边界处沿辐射波前的定域传播方向 产生可以探测的变化；和

在所说辐射横穿所说边界之后探测和记录至少一部分所说辐射，所说 辐射以这样一种方式横穿所说边界，使得可以观测到沿定域传播方向的所 说变化，从而可以将这些变化记录为辐射强度的定域衰减或急剧变化，这 一记录基本上形成所说边界的图象。

本发明还提供一种获得具有急剧折射率变化或由厚度变化限定的内部 边界的相差衬托记录的设备，所说设备包括：

用X-射线辐射辐照所说边界的装置，所说X-射线辐射具有垂直于所 说折射率变化方向或者所说厚度变化方向的显著的传播向量分量，并且对 于折射率或厚度的变化具有足够强的横向空间相干性，以在所说边界处沿 辐射波前的定域传播方向产生可以探测的变化；

用于在所说辐射横穿所说边界之后探测和记录至少一部分所说辐射的 装置，所说辐射以这样一种方式横穿所说边界，使得可以观测到沿定域传 播方向的所说变化，从而可以将这些变化记录为辐射强度的定域衰减或急 剧变化，这一记录基本上形成所说边界的图象。

本发明还提供一种获得一个物体的边界图象的方法，所说边界表征折 射率变化，所说方法包括以下步骤：

用具有强横向空间相干性和垂直于所说折射率变化的传播分量的穿透 辐射照射所说边界；和

在一个成象平面接收至少一部分所说辐射以构成所说图象，所说辐射 已经被所说边界菲涅耳衍射，从而所说边界在所说图象中表现为相应的强 度变化。

本发明还提供一种用于获得一个物体的边界图象的设备，所说边界表 征折射率的变化，所说设备包括：

一个辐射源，其用具有强横向空间相干性和垂直于所说折射率变化的 传播分量的穿透辐射照射所说边界；和

一个探测器，用于接收至少一部分所说辐射，以便构成所说图象，所 说辐射已经被所说边界菲涅耳衍射，从而所说边界在所说图象中表现为 相应的强度变化。

本发明还提供一种确定一幅图象相位的方法，该方法包括处理所说图 象的相衬图象数据的步骤。

在定域传播方向的强度变化在包含所说记录的图象中最好是可以观察 出的。所说记录，从而所说图象可以是照相构成的或者是电子构成的。因 此，术语“图象”指的是，例如，对一组强度数据，例如一个强度值表或 者强度值的其它存储记录的可观测效果；这个术语不局限于视觉可视的含 义。记录媒体可以包括一个两维象素化探测器，例如一个电子探测器，诸 如电荷耦合器件(CCD)阵列。

所说辐照装置最好包括一个直径为20微米或更小的X-射线源，其中直 径指的是辐射源在最大强度的一半时的强度分布的全宽度。该装置最好还 包括用于支撑包含被成象的内部边界的样品的平台或支架。

所说穿透辐射，例如X-射线辐射可以是多色的，最好是在硬X-射线范 围内，即在1keV到1MeV的范围内。

可取的是，选择边界与探测装置的间距以提高图象的分辨率。例如， 已经发现通过增加间距可以获得较为清晰的图象，即具有较好对比度的图 象。举例来说，至少相对于0.4米的间距大约1米的间距可以改善对比度。 其部分原因是随着间距的增加背景噪声降低，但是基本保留了由于定域传 播方向的改变而产生的强度变化效果。

本申请中的术语“横向空间相干性”指的是在垂直于波的传播方向的 不同点处的辐射波的复数振幅的相关性。据说当一个波前上的各个点具有 不随时间变化的传播方向时就会产生横向空间相干性。实际上，通过使用 具有较小有效尺度的辐射源或者通过在距离辐射源较大距离处观测辐射光 束可以获得较强的横向空间相关性。例如，对于20keV的X-射线，通常20 微米直径或者更小直径的辐射源是适合的。在辐射源所产生的总辐射通量 足够的前提下，源的尺度越小，对实现本发明的目的越有益。可能需要通 过仔细选择辐射源的X-射线窗口，例如使得该源具有高度一致的厚度和均 匀性来保持横向空间相干性。

下面参照附图介绍本发明的优选实施例，但是仅仅是以示例的方式予 以说明，在所说附图中：

图1为表示被一平行光束辐照的一个圆形横碱面物体的示意图，为了 说明分成三个部分；

图2为表示被一多色光束辐照的一个圆形横截面物体和所产生的相衬 图象的强度的示意图；

图3为表示本发明的一个实施例中X-射线光学配置的示意图；

图4和图5为根据本发明获得的各种边界的X-射线图象，如下所详述。

现在首先概述本发明的数学原理。

样品厚度和X-射线折射率，n(λ)＝1-δ(λ)-iβ(λ)的变化必定会 引起通过样品的X-射线波前形状的改变。折射率n的实部分量1-δ(λ)与折 射程度有关，虚部分量-iβ(λ)与吸收程度有关。更具体地说，对于一个单 元的物质 $\delta \left(\lambda \right)=\frac{{\gamma }_0{\lambda }^2}{2\pi }{N}_O{f}_R---\left(1\right)$ $\beta \left(\lambda \right)=\frac{\lambda }{4\pi }\mu \left(\lambda \right)---\left(2\right)$

其中μ(λ)为线性吸收系数，r0为电子的经典半径，N0为每单位本积中 的原子数量，fR为零散射角度时原子散射因数的实部。系数δ正比于λ2，β 正比于λ4，λ还正比于从辐射源发射的X-射线光子的1/能量。

波前畸变的幅度与垂直于波前传播方向的相位变化梯度相关。在几何 光学近似中，通过一个物体的射线路径的相位差，φ，正比于沿该射线路径 的折射率的实部，δ，减量的积分。对于图1所示的坐标系，这个相位差一 般可以表示为 $\phi (x,z)-k{\int }_0^z\{n(x,z^{\prime })-1\}d{z}^{\prime }---\left(3\right)$

其中k等于2π/λ。定域散射波矢与定域入射波矢的角度偏差Δα正比 于垂直于定域入射波矢方向的相位差梯度。“定域”一词指的是波前上的 一个点(x，y，z)。对于图1所示的坐标系定域散射波矢的数学表示为：

其中s(x，y，z)为波前在点(x，y，z)的法线，如果(φ/x)2+(φ/y)2＜＜ k2，则上述关系式在近轴近似下是有效的。于是角度偏差Δα可以表示为： $\mathrm{\Delta \alpha }\approx \frac{1}{\kappa }\frac{\partial \phi (x,z)}{\partial x}={\int }_0^z\{\frac{\partial n(x,z^{\prime })}{\partial x}-1\}{\mathrm{dz}}^{\prime }----\left(5\right)$

从而角度偏差Δα依赖于垂直于传播波矢的折射率变化，偏差量依赖于 沿波矢方向产生变化的长度，例如样品的厚度。

为了说明该效应的特性，考虑在折射率为n0＝1的介质中置入折射率为 nM的一个球形体，Ω，的情况，如图1和图2所示。

X-射线通过样品相对于通过真空的光程差导致产生相位差φ(x)，进而 沿垂直于传播定域方向的方向(图1)产生相位梯度φ/x。以恒定距离平 行于z轴通过物体Ω的光线1与基准光线0之间的相位差为： $\phi (x,y)=\frac{2\pi }{\lambda }{\int }_{-z(x,y)}^0\delta \left(\lambda \right)\mathrm{dl}=\frac{2\pi }{\lambda }\delta \left(\lambda \right)z(x,y)---\left(6\right)$ 其中z(x，y)为光线1与物体Ω的相交长度，表示为 $z(x,y)=2\sqrt{R^2-x^2-y^2}\left(7\right)$

R为物体Ω的半径，δ为折射率系数的减少量。对于一个位于x-z平面 中的圆形截面物体，φ/x和入射光线与相应的折射光线之间的角度偏差 Δα在一个给定点x处的数学关系式为： $\mathrm{\Delta \alpha }=\frac{1}{k}s(x,y,z)$ $=$$\frac{\lambda }{2\pi }\frac{\partial \phi }{\partial x}--\left(8\right)$ $=-2\delta \left(\lambda \right)\left[\frac{x}{\sqrt{R^2-x^2}}\right]$

在公式(8)中，δ(λ)是缓慢变化的，很容易看出该相位梯度收敛在x＝±R 点，在这些点光线可以偏离光轴非常大的角度。在这些极限值处，散射光 束的角度偏差可能非常大，并且导致沿相应的向前方向的强度I的显著损 失，其位置与波长无关，如图2所示多色光B的情形。在轻元素和硬X- 射线的情况下，折射率系数的减少量，δ，通常为10-5到10-6量级，但是无 论怎样，当x接近±R点时，即在样品的边界处或者在内部边界特征点处， 角度偏差Δα都会非常大。

在有关源的大小尺度、物体与源之间距离和物体与图象之间距离、以 及源的光谱分布等不同条件下获得的对比度的特性需要考虑。另一个需要 考虑的影响对比度的因素是物体造成的波前改变的程度。

在平面波的情况下，为了帮助理解这些因素对于形成图象的对比度的 作用，我们可以在一级近似下使用由Cowley(J.M.Cowley，“衍射物理”， 第二版，60页，North Holland出版，1981)推导出的有关相位体的菲涅 耳衍射对比度的公式。根据该公式，可以得出，对于一个产生相差φ(x)的 一维相位体，在波长为λ的平面波照射下，在距离该物体R2处的强度分布 为： $I\left(x\right)=1+\frac{R_2\lambda }{2\pi }\phi ″\left(x\right)\left(9\right)$

假定量值(R2λ/2π)φ(x)很小，则上式在一级近似下是有效的。从这个 非常简单的公式，可以得出一些明确的结论，亦即：

i)对比度直接随R2变化，

ii)图象的构形不依赖于波长λ。只有对比度受到影响。对于一个多色 光源，只需用光谱加权的和值代替式中的λ。

为了取得对于上式在X-射线情况下的有效性范围的认识，我们假定一 个物体具有这样的特征，它使得该物体在10微米的横向距离上产生1弧度 的相位变化。于是φ″～1010m-2，在λ～1，R2～1m的情况下，我们看到 (R2λ/2π)φ″(x)≤1。从而该公式即使对于小的相位体或者相位的合理的迅速 变化也应是有效的。但是，对于斜率的非常陡峭的边沿或者变化，诸如在 计算模拟测试物体(例如，纤维)中经常使用的，φ″可能变得太大(甚至 无穷大)，使得该公式失效。但是即使在这些情况下，图象的基本构形(对 于一个清晰边缘的台阶物体为黑/白线)也是可以再现的，然而无法再现由 于这种不连续性产生的衍射的典型的次要边缘线。另一方面，或许实际上 更为重要的是，对于较小的φ″(x)，即具有较少急剧横向变化的较大特征 物，对比度较低，而且可能限制了实际的可视性。

最近已经有人(P.Cloetens，R.Barrett，J.Baruchel，J.P.Gurgay和 M.Schlenker，J.Phys.D.：Appld.Phys.，1996，29，133- 46；J.P.Guigay，Optik，1977，49，121-5)对于采用平面波形成的这类图象以 菲涅耳衍射为基础进行了更加精确的数学论述。他们的研究得出与上述一 级近似相同的公式。但是更为准确的数学处理得出这样的结论，对于空间 频率为u的光束，当2λR2u2＝1时，至少在相衬射线成象预期条件的正常范 围内，形成最大的对比度。空间频率u与成象物体的结构有关，其中u＝1/A， 其中A为成象物体的富里叶分量的空间周期。

这些论述都是以用理想的平面波照射为前提的。光束的发散性会使图 象模糊，这个模糊量正比于R2(在这方面其影响与常规的射线照相相同)。 于是上述作者(Cloetens等人)在兼顾对比度和分辨率的情况下，给出总 的最佳R2值为：

R2≤2λ/α2                              (10)

其中α＝s/R1为光源相对物体的角度，并且涉及(几乎是)平面波情况。 应当指出Cloetens等人特别要求需要使用高度单色的X-光源，并且在本申 请所述的优选实施例中仅仅考虑平面波情况。

正如所指出的，尽管我们从原理上阐述了与常规的射线照相更为接近 的球面波情况，但是上面的论述主要涉及平面波情况。为了帮助理解空间 波的情况，我们现在考虑两者之间的关系，通过对于用距离物体R1的一个 点源来成象一个物体(球面波情况)的菲涅耳-克契霍夫表示式的简单分 析可以建立这种关系。这表明对于球面波情况存在一个简单的关系式，其 中包含平面波项，但是修正了物体到成象平面的距离R，即 $\frac{1}{R}=\frac{1}{R_1}+\frac{1}{R_2}---\left(11\right)$

并且将图象放大了(R1+R2)/R1倍。根据简单的几何理论，按照射线光学， 在球面波情况下由于辐射源尺度造成的对比度或者分辨率的损失将不成为 问题，因为图象和源尺度都被放大了，后者放大了R2/R1倍，对于较大的R2， 两者渐近于相同的因数。在球面波情况下影响对比度的因素是(在与射线 照相相关的能量和空间分辨率范围)2λR2(1+R2/R1)u2应当较大(但是通常 小于1)。这个表示式的结果由于R2较大，或者λ较大，或者空间频率u 较大而较大。作为一个示例性说明，对于实际的射线照相可以给出下列值 λ＝0.2；u≤2×105(对应于20微米或更大的空间周期)，从而R2≈2.5m (假设R2/R1＝3)对于最高的空间频率将给出最大对比度。较大的R2值适合 于在较低空间频率下得到最大的对比度。

可以指出，函数φ″倾向于增强在一个图象中相物体的边沿和边界。如 果还存在该物体的一个吸收分量，则至少在一级近似下该分量会直接增大 图象对比度(例如参见Guigay论文所述的方程7，1977)。该技术能够 补充和增强常规的射线图象，并且得到新的信息。还可以指出，对于包含 不同的相衬(包含φ的拉普拉斯函数)的图象中的对比度的适当的处理需要 借助例如强度方程的传递解对图象进行数值处理(参见 T.E.Gureyev，A.Roborts和K.Nugent，Journal of Optical Society of America，Vol.A12，1932页和1942页，1995年，此文以引用方式结合 在本申请中)，以修正相位φ(x)。

我们现在讨论应用从这些推定中得出的概念的实际设置。在第一实施 例中(图3)，使用了一个具有强空间相干性的辐射源S和一个X-射线成 象探测器D，例如底片、可光作用的荧光板(例如富士成象板)、或者一 个两维电子探测器。垂直于传播方向的折射率急剧变化、或者沿传播方向 的厚度变化的区域能够引起通过这些区域的波前的定域传播方向的明显改 变。因此，从点源S发射的球面波前W1在经过物体0时畸变为W2。通过 在距样品足够距离处记录波前的强度，可以探测到由于样品中折射率的急 剧变化或者厚度的变化产生的强度变化，并且将它们的位置记录在图象 中。这相当于相衬成象的一种形式。成象探测器的位置是经过选择的，从 而使探测器的空间分辨率足以分辨由于波前的急剧畸变产生的强度差别， 并且如上所述，使对比度在实际应用中达到最佳。

通常，折射率或厚度的急剧跃变在图象中成象为相应点处强度的急剧 损失或迅速变化。在图象中给定点处强度损失或迅速变化的这种特性基本 上与波长无关，所以即使在使用多色光时也会在图象中产生非常清晰的对 比度变化。

这种配置的特征在于，对于一个环形源分布，图象的空间分辨率对于 各个方向是相同的，并且基本上由源的大小确定。它的优点还在于能够显 著地放大图象，从而可以使用象富士成象板这样的记录媒体，这种成象板 具有许多有益的特性，比如较宽的动态范围和较高的灵敏度，但是不具有 较高的空间分辨率。

在这种配置中除了包含辐射源和探测器以外，在样品和探测器之间还 可以插入一个高分辨率角度分析仪。该高分辨率角度分析仪可以是例如一 个适当弯曲的晶体，所说晶体具有劳埃几何形状，其曲率经过选择以适合 辐射源的某些特征波长。对这种方法改进的目的在于分辨比在前述配置中 可以观测到的样品折射率和厚度变化更弱的变化。

可以指出，非常显著地放大图象是可能的，从而即使使用空间分辨率 非常低的探测器诸如富士成象板也可以使图象具有非常高的空间分辨率。 还可以指出，由于这种成象方法基本上与X-射线能量无关，所以可以让光 源在高管电压下工作，从而对样品产生较少的吸收剂量，这在临床应用中 十分重要。

利用上述技术记录的相衬图象的一些实例表示在图4和图5中。图4 表示了一个10微米厚的塑料薄膜的边缘的图象，这种塑料片与 Davis，Gao，Gureyev，Stevenson和Wilkins所用的是相同的(参见 Phys.Rev.Letters，1995，Vol.74，3173页)，并且相当于一个纯粹的相位 体。图5表示了在一个聚合物基质中的一个气泡和玻璃纤维的图象，该样 品类似于Davis，Gao，Gureyev，Stevenson和Wilkins所报告使用的样品(参 见Nature Vol.373，595-598页，1995)，并且相当于一个近乎纯粹的相 位体。在每种情况下都能够看到具有比在常规的吸收对比度图象中预期的 更清晰的对比度。特别是，在图4中薄膜的边缘表现为黑/白对比度特征是 清晰可见的，气泡和纤维的边缘也是如此。所使用的辐射源为常规的10微 米直径微焦距源(Kevex Model PXS)，其铜阳极工作电压为40kV。对于 图4所示的配置，源到样品和样品到底片的距离都是700毫米，而图5中 相应的距离分别为120毫米和1000毫米。应当指出，由于源具有强空间相 干性，在这些实施例中的对比度几乎是完全可见的。这种对比度主要是强 度损失对比度，从这个意义上说类似于常规的吸收对比度，但是不同之处 在于它表示由于在物体边界处发生的折射散射(或菲涅耳衍射)造成的强 度损失，如方程(8)所示。一个常规的直径为0.1毫米的微焦距源具有 大约0.1毫米长的投影，如在照片上的标尺刻度所示，因此使得对比度大 大模糊。

为了比较在本申请中所述的相衬成象与标准的吸收成象，下列的表给 出了在不同的源能量E条件下，实现65％吸收所需的碳样品的吸收厚度ta 和实现φ的2π的相位变化所需的样品的相位厚度tp。

表1     EkeV     λ()     ta(μm)     tp(μm)     50     0.25     435000     133     12     1     5000     30     1.2     10     4     3     0.25     50     1.3     1.2

该表的结果表明可以如何将相衬成象应用于利用高能光源对非常小的 物体成象。

可取的是，在样品与探测器之间的光路上可以设置具有X-射线可透过 的窗口的真空管或者类似器件以减少空气散射的影响，并且确保它们的光 学质量对于X-射线束的相干性不具有不利的影响。

本发明的方法尤其适合于对各种物质中的诸如裂缝、空隙和托层特征 进行成象，因为这些特征可使X-射线折射率产生最大偏差，并且空间变化 非常急剧。为了给出可以观测到的对比度，光源最好具有非常小的有效尺 度，例如说小于20微米的量级，探测器最好是一个高分辨率成象探测器， 比如X-射线底片或两维电子探测器，例如CCD阵列。该方法还有助于在医 疗射线成象中显著提高重要特征的对比度。

本申请概述了使用穿透辐射诸如硬X-射线进行差分相衬成象的一些简 化条件和配置，它们特别适合于医疗和工业应用。这些新方法与用于吸收 对比度射线照相的传统方法更为接近，并且应当比我们先前在上述的 WO95/05725和PN5811/95中所述的方法更容易实施，特别是在大面积照射 的情况下。在给定源功率的情况下，由于它们可以使用非常宽光谱的源， 它们还比早先的单色方法具有明显较短的曝光时间。

在整个说明书中，除了需要采用其它的含义，单词“包括”(单数形 式)，或变型诸如“包括”(复数形式)或“包括”(分词形式)应当理 解为意味着包含所述的整体或整体的集合，而不是意味着排除其它的整体 或整体的集合。