一种内窥式X射线发光断层成像装置，包括：激发源单元(1)，用于从外部激发源激发受检对象内部的探针；信号采集单元(2)，用于采集受检对象的X射线图像和所述探针受激发后发出的近红外光信号图像；计算单元(4)，用于对采集的X射线图像和红外光信号图像进行预处理，对预处理后的X射线图像进行稀疏重建，以获取受检对象的结构信息和信号采集单元中的内窥探头在受检对象内部的位置信息，以及利用获取的结构信息和位置信息，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部的靶向目标的位置信息和分布信息。
激发源单元(1)，用于从外部激发源激发受检对象内部的探针；
信号采集单元(2)，用于采集受检对象的X射线图像和所述探针受激发后发出的近红外光信号图像；
计算单元(4)，用于对采集的X射线图像和红外光信号图像进行预处理，对预处理后的X射线图像进行稀疏重建，以获取受检对象的结构信息和信号采集单元中的内窥探头在受检对象内部的位置信息，以及利用获取的结构信息和位置信息，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部的靶向目标的位置信息和分布信息。
2.根据权利要求1所述的成像装置，其中，所述信号采集单元包括近红外光信号采集模块(22)，所述近红外光信号采集模块(22)包括气球式成像导管(221)、内窥探头(222)、信号传输光纤束(223)和光电倍增管阵列探测器(224)。
3.根据权利要求2所述的成像装置，其中所述气球式成像导管(221)包括医用导管，所述医用导管的一端连接于信号传输光纤束(223)中的空气输送管(2232)，另一端固定于信号传输光纤束(223)的外表面，并将内窥探头(222)包裹其中，用于采集近红外光信号图像时充入空气以撑起空腔类器官的腔壁。
操纵单元(3)，包括C型旋转臂(31)，用于多角度地采集受检对象的X射线投影图像；以及成像体固定床(32)，用于固定受检对象，并进行调整以使受检对象位于X射线成像的视场内。
5.根据权利要求2所述的成像装置，其中，所述内窥探头包括N对鼓形透镜(2221)和棒透镜(2222)，每对鼓形透镜(2221)和棒透镜(2222)分别与信号传输光纤束(223)中的信号传输光纤(2231)相连，用于采集所述探针受激后发出的近红外光信号图像。
6.根据权利要求2所述的成像装置，其中，所述信号传输光纤束(223)包括N根信号传输光纤(2231)，用于将采集的近红外光信号图像传输到光电倍增管阵列探测器(224)；和空气输送管(223)，用于向气球式成像导管(221)中充入空气。
7.根据权利要求1-6之一所述的成像装置，其中，所述探针是荧光纳米颗粒探针。
8.根据权利要求1所述的成像装置，其中，所述计算单元显示光学三维重建的结果，包括对获得的受检对象内部靶向目标的重建结果和所述结构信息进行图像融合，以便对融合图像进行三维显示。
9.根据权利要求1所述的成像装置，其中，计算单元包括图像预处理模块，用于对X射线图像进行预处理，包括噪声去除、坏点补偿以及图像校正。
10.根据权利要求1或9所述的成像装置，其中，计算单元包括图像预处理模块，对近红外光信号图像进行预处理，包括噪声去除、感兴趣区域提取以及坏点补偿。
11.根据权利要求1所述的成像装置，其中，计算单元包括图像重建模块，用于：构建前向光传输混合模型；
根据前向光传输混合模型，建立内窥式X射线发光断层成像的总系统方程，获得近红外光信号图像的计算值向量；
根据获得的近红外光信号图像的计算值向量与采集的近红外光信号图像的测量值向量之间的误差，结合探针的稀疏分布和非负约束，建立稀疏正则化目标函数；
对建立的目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而获得受检对象内部的靶向目标的位置信息和分布信息。
12.根据权利要求11所述的成像装置，其中，构建前向光传输混合模型包括：基于获得的内窥探头在受检对象内部的位置信息，根据受检对象生物组织的光学特性参数，将生物组织划分为低散射特性组织、高散射特性组织、空腔；分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程；构造边界耦合条件来耦合这三个光传输方程，同时结合X射线激发探针的过程，构建内窥式X射线发光断层成像的前向光传输混合模型。
13.根据权利要求11或12所述的成像装置，其中，建立总系统方程包括：利用下式，基于有限元方法对构建的前向光传输混合模型进行离散化，组装每个离散点上的子系统方程：Akρ＝Pk，
其中，Ak是第k个X射线照射角度对应的系统矩阵，取决于生物体内三大类组织的分布、光学特性参数和X射线照射到体表的位置；ρ是荧光纳米颗粒探针的浓度；Pk是与第k个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值；
通过联立M个X射线照射角度对应的系统方程，建立内窥式X射线发光断层成像的总系统方程：Aρ＝P，
其中，A是总系统矩阵，定义为A＝[A1，A2，…，AM]T；
P是与所有M个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值向量，定义为P＝[P1，P2，…，PM]T；
14.根据权利要求13所述的成像装置，其中，利用下式建立稀疏正则化目标函数：其中，Pm是与所有M个X射线照射角度对应的、所述采集的近红外光信号图像的测量值向量；
对所建立的稀疏正则化目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而间接获取受检对象内部的靶向目标的位置信息和分布信息。
15.根据权利要求12所述的成像装置，其中，按照如下规则对生物组织的类型进行划分：将荧光纳米颗粒探针所在组织、体表所属组织以及约化散射系数与吸收系数比值小于
10的组织划分为低散射特性组织，将约化散射系数与吸收系数比值大于等于10的组织划分为高散射特性组织，将近红外光信号图像采集时内窥探头所处的组织划分为空腔。
16.根据权利要求12所述的成像装置，其中，按照如下规则分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程：采用三阶简化球谐波方程描述近红外光在低散射特性组织中的传输过程，采用扩散方程描述近红外光在高散射特性组织中的传输过程，采用自由空间光传输方程描述近红外光从空腔与散射特性组织之间的边界到内窥探头的传输过程。
17.根据权利要求12所述的成像装置，其中，应用下式来耦合低散射特性组织和高散射特性组织的光传输方程：Ф(r)＝α(r)(β1(r)φ1(r)+β2(r)φ2(r))，其中，r是低散射特性组织与高散射特性组织相交的边界上的任意一点；
φ1(r)和φ2(r)是采用三阶简化球谐波方程计算得到的点r处的光通量分量；
β1(r)和β2(r)是三阶简化球谐波方程的边界折射率不匹配系数。
18.根据权利要求17所述的成像装置，其中应用下式耦合散射特性组织与空腔的光传输方程：式中，r是散射特性组织与空腔相交的边界上的任意一点；
σ是散射特性组织的标识因子，若低散射特性组织与空腔相交，则σ的值取0；若高散射特性组织与空腔相交，则σ的值取1。
(a)采集受检对象内部的探针受激后发射的近红外光信号图像，并采集受检对象的X射线图像；
(b)对采集的近红外光信号图像和X射线图像进行图像预处理；
(c)对预处理后的X射线图像进行稀疏重建，以获得受检对象的结构信息和内窥探头在受检对象内部的位置信息；
(d)利用获取的结构信息和位置信息，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部的靶向目标的位置信息和分布信息。
20.根据权利要求19所述的方法，还包括显示光学三维重建的结果，包括对获得的受检对象内部靶向目标的重建结果和所述结构信息进行图像融合，以便对融合图像进行三维显示。
21.根据权利要求19所述的方法，其中，对X射线图像进行的图像预处理包括噪声去除、坏点补偿以及图像校正。
22.根据权利要求19所述的方法，其中，对近红外光信号图像进行的图像预处理包括噪声去除、感兴趣区域提取以及坏点补偿。
23.根据权利要求19所述的方法，其中，所述步骤(d)包括(d1)构建前向光传输混合模型；
(d2)根据前向光传输混合模型，建立内窥式X射线发光断层成像的总系统方程，获得近红外光信号图像的计算值向量；
(d3)根据获得的近红外光信号图像的计算值向量与采集的近红外光信号图像的测量值向量之间的误差，结合探针的稀疏分布和非负约束，建立稀疏正则化目标函数；
(d4)对建立的目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而获得受检对象内部的靶向目标的位置信息和分布信息。
24.根据权利要求23所述的方法，其中，步骤(d1)包括：基于获得的内窥探头在受检对象内部的位置信息，根据受检对象生物组织的光学特性参数，将生物组织划分为低散射特性组织、高散射特性组织、空腔；分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程；构造边界耦合条件来耦合这三个光传输方程，同时结合X射线激发探针的过程，构建内窥式X射线发光断层成像的前向光传输混合模型。
25.根据权利要求23或24所述的方法，其中，步骤(d2)包括利用有限元方法对构建的前向光传输混合模型进行离散化，组装每个离散点上的子系统方程，建立第k个X射线照射角度对应的系统方程：Akρ＝Pk，
其中，Ak是第k个X射线照射角度对应的系统矩阵，取决于生物体内三大类组织的分布、光学特性参数和X射线照射到体表的位置；ρ是荧光纳米颗粒探针的浓度；Pk是与第k个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值；
通过联立M个X射线照射角度对应的系统方程，建立内窥式X射线发光断层成像的总系统方程：Aρ＝P，
其中，A是总系统矩阵，定义为A＝[A1，A2，…，AM]T；
P是与所有M个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值向量，定义为P＝[P1，P2，…，PM]T；
26.根据权利要求25所述的方法，其中，所述步骤(d3)包括：建立下列稀疏正则化目标函数：其中，Pm是与所有M个X射线照射角度对应的、所述采集的近红外光信号图像的测量值向量；
所述步骤(d4)包括对所建立的目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而间接获取受检对象内部的靶向目标的位置信息和分布信息。
27.根据权利要求24所述的方法，其中，按照如下规则对生物组织的类型进行划分：将荧光纳米颗粒探针所在组织、体表所属组织以及约化散射系数与吸收系数比值小于10的组织划分为低散射特性组织，将约化散射系数与吸收系数比值大于等于10的组织划分为高散射特性组织，将近红外光信号图像采集时内窥探头所处的组织划分为空腔。
28.根据权利要求24所述的方法，其中，按照如下规则分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程：采用三阶简化球谐波方程描述近红外光在低散射特性组织中的传输过程，采用扩散方程描述近红外光在高散射特性组织中的传输过程，采用自由空间光传输方程描述近红外光从空腔与散射特性组织之间的边界到内窥探头的传输过程。
29.根据权利要求28所述的方法，还包括应用下式来耦合低散射特性组织和高散射特性组织的光传输方程：Ф(r)＝α(r)(β1(r)φ1(r)+β2(r)φ2(r))，其中，r是低散射特性组织与高散射特性组织相交的边界上的任意一点；
φ1(r)和φ2(r)是采用三阶简化球谐波方程计算得到的点r处的光通量分量；
β1(r)和β2(r)是三阶简化球谐波方程的边界折射率不匹配系数。
30.根据权利要求29所述的方法，其中，应用下式耦合散射特性组织与空腔的光传输方程：式中，r是散射特性组织与空腔相交的边界上的任意一点；
σ是散射特性组织的标识因子，若低散射特性组织与空腔相交，则σ的值取0；若高散射特性组织与空腔相交，则σ的值取1。
内窥式X射线发光断层成像装置及方法
技术领域
[0001] 本发明涉及医学影像技术领域，更进一步，涉及一种内窥式X射线发光断层成像装置以及一种内窥式X射线发光断层成像方法。
背景技术
[0002] 由于能够三维地重建生物体内荧光探针的空间位置和浓度分布，扩散荧光断层成像(Fluorescence diffuse optical tomography，FDOT)已经成为一种很有应用潜力的分子影像技术。在其潜在的应用中，FDOT已经广泛地应用于小动物或表层器官成像，如脑部和乳腺等。这种情况下，荧光信号可以穿透成像物体，并被外部的探测部件接收。但是，对于深埋于例如大型动物体或人体的受检对象的内部器官，由于荧光在体内的严重衰减，FDOT技术无法提供准确的图像，甚至于无法测量体表的荧光信号。
[0003] 内窥成像技术通过生物体内的空腔类器官，如肠道、食管道、胃道、膀胱等，将探测部件深入体内并靠近感兴趣区域，可以获得内部器官的精确成像。鉴于内窥成像技术的优势，Piao(Daqing Piao，Hao Xie，Weili Zhang，and Jerzy S.Krasinski，“Endoscopic，rapid near-infrared optical tomography，”Optics Letters 2006，31(19)，2876-2878)等提出了一种内窥式扩散荧光断层成像(endoscopic fluorescence diffuse optical tomography，eFDOT)系统，包括eFDOT装置和方法，用于对生物体内部器官的成像。
[0004] 然而，这种eFDOT系统在实际应用过程中存在一些问题，严重制约着进一步的临床应用。
[0005] 首先，由于外部近红外激发光源的照射，生物体内将会产生很强的自体荧光，这将严重影响重建图像的质量。自体荧光的去除也是FDOT技术中存在的棘手问题，一般采用图像后处理算法将其去除；但这并未从根本上消除自体荧光，并且严重依赖于图像后处理算法的性能。
[0006] 其次，在已有的eFDOT装置中，采用直径在厘米级的微型化内窥探头，其中激发光纤和探测光纤等间隔地、交叉分布于内窥探头的顶部侧面。由于内窥探头直径很小，因此激发光纤与探测光纤之间的距离很小。如此小的激发源-探测器距离，将会导致信号间的强干扰，从而使成像结果不再准确。
[0007] 第三，在已有的eFDOT方法中，采用扩散方程(Diffusion Equation，DE)作为前向光传输模型。由于自身的局限性，在低散射区域、光源和折射率不匹配组织边界附近，DE不能准确描述光的传输过程。因此，在eFDOT方法中，DE将不再适用。
[0008] 近年来，一种名为X射线发光断层成像(X-ray luminescence computed tomography，XLCT)的分子影像技术用于小型动物成像。XLCT采用X射线激发小型动物体内的荧光纳米颗粒(Phosphor Nanoparticles，PNP)探针，能够发射近红外光；通过CCD相机采集发射的近红外光，采用近红外光学断层成像技术可以重建小型受检对象内部的PNP探针分布。这种成像技术由于采用X射线作为激发源，能够激发处于生物体组织深部的PNP探针，并且能够去除自体荧光。然而，由于发射的近红外光穿透深度的限制，XLCT技术仍只能用于小型动物体和表层器官的成像。
发明内容
[0009] 因此，本发明实施例提供了一种内窥式X射线发光断层成像装置以及一种内窥式X射线发光断层成像方法，以实现对受检对象内部的靶向目标的精确成像。
[0010] 根据本发明实施例的一方面，公开了一种内窥式X射线发光断层成像装置，包括：激发源单元，用于从外部激发源激发受检对象内部的探针；信号采集单元，用于采集受检对象的X射线图像和所述探针受激发后发出的近红外光信号图像；计算单元，用于对采集的X射线图像和近红外光信号图像进行预处理，对预处理后的X射线图像进行稀疏重建，以获取受检对象的结构信息和信号采集单元中的内窥探头在受检对象内部的位置信息，以及利用获取的结构信息和位置信息，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部的靶向目标的位置信息和分布信息。
[0011] 根据优选实施例，信号采集单元包括近红外光信号采集模块，所述近红外光信号采集模块包括气球式成像导管、内窥探头、信号传输光纤束和光电倍增管阵列探测器。
[0012] 根据优选实施例，气球式成像导管包括医用导管，医用导管的一端连接于信号传输光纤束中的空气输送管，另一端固定于信号传输光纤束的外表面，并将内窥探头包裹其中，用于采集近红外光信号图像时充入空气以撑起空腔类器官的腔壁。
[0013] 根据优选实施例，成像装置还包括：操纵单元，包括C型旋转臂，用于多角度地采集受检对象的X射线投影图像；以及成像体固定床，用于固定受检对象，并进行调整以使其位于X射线成像的视场内。
[0014] 根据优选实施例，内窥探头包括N组鼓形透镜和棒透镜对，每对鼓形透镜和棒透镜均与信号传输光纤束中的信号传输光纤相连，用于采集所述探针受激后发出的近红外光信号图像。
[0015] 根据优选实施例，信号传输光纤束包括N根信号传输光纤，用于将采集的近近红外光信号图像传输到光电倍增管阵列探测器；和空气输送管，用于向气球式成像导管中充入空气。
[0016] 根据优选实施例，探针是荧光纳米颗粒探针。
[0017] 根据优选实施例，计算单元显示光学三维重建的结果，包括对获得的受检对象内部靶向目标的重建结果和所述结构信息进行图像融合，以便对融合图像进行三维显示。
[0018] 根据优选实施例，计算单元包括图像预处理模块，用于对X射线图像进行预处理，包括噪声去除、坏点补偿以及图像校正。
[0019] 根据优选实施例，计算单元包括图像预处理模块，对近红外光信号图像进行预处理，包括噪声去除、感兴趣区域提取以及坏点补偿。
[0020] 根据优选实施例，计算单元包括图像重建模块，用于：构建前向光传输混合模型；根据前向光传输混合模型，建立内窥式X射线发光断层成像的总系统方程，获得近红外光信号图像的计算值向量；根据获得的近红外光信号图像的计算值向量与采集的近红外光信号图像的测量值向量之间的误差，结合探针的稀疏分布和非负约束，建立稀疏正则化目标函数；对建立的目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而获得受检对象内部的靶向目标的位置信息和分布信息。
[0021] 根据优选实施例，构建前向光传输混合模型包括：基于获得的内窥探头在受检对象内部的位置信息，根据受检对象生物组织的光学特性参数，将生物组织划分为低散射特性组织、高散射特性组织、空腔；分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程；构造边界耦合条件来耦合这三个光传输方程，同时结合X射线激发探针的过程，构建内窥式X射线发光断层成像的前向光传输混合模型。
[0022] 根据优选实施例，建立总系统方程包括：利用下式，基于有限元方法对构建的前向光传输混合模型进行离散化，组装每个离散点上的子系统方程：
[0023] Akρ＝Pk
[0024] 其中，Ak是第k个X射线照射角度对应的系统矩阵，取决于生物体内三大类组织的分布、光学特性参数和X射线照射到体表的位置；ρ是荧光纳米颗粒探针的浓度；Pk是与第k个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值；
[0025] 通过联立M个X射线照射角度对应的系统方程，建立内窥式X射线发光断层成像的总系统方程：
[0026] Aρ＝P
[0027] 其中，A是总系统矩阵，定义为A＝[A1，A2，…，AM]T；
[0028] T是矩阵转置操作符；
[0029] P是与所有M个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值向量，定义为P＝[P1，P2，…，PM]T；
[0030] M是选取的X射线照射角度的个数。
[0031] 根据优选实施例，利用下式建立稀疏正则化目标函数：
[0032]
[0033] 其中，Pm是与所有M个X射线照射角度对应的、所述采集的近红外光信号图像的测量值向量；
[0034] λ是稀疏正则化参数；
[0035] ||F||1定义为求解矩阵F的L1范数；以及
[0036] 对所建立的稀疏正则化目标函数进行求解，获得探针在受检对象内部的空间位置和浓度分布，从而间接获取受检对象内部的靶向目标的位置信息和分布信息。
[0037] 根据优选实施例，按照如下规则对生物组织的类型进行划分：将荧光纳米颗粒探针所在组织、体表所属组织以及约化散射系数与吸收系数比值小于10的组织划分为低散射特性组织，将约化散射系数与吸收系数比值大于等于10的组织划分为高散射特性组织，将近红外光信号图像采集时内窥探头所处的组织划分为空腔。
[0038] 根据优选实施例，按照如下规则分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在这三大类组织中的传输过程：采用三阶简化球谐波方程描述近红外光在低散射特性组织中的传输过程，采用扩散方程描述近红外光在高散射特性组织中的传输过程，采用自由空间光传输方程描述近红外光从空腔与散射特性组织之间的边界到内窥探头的传输过程。
[0039] 根据优选实施例，应用下式来耦合低散射特性组织和高散射特性组织的光传输方程：
[0040] Φ(r)＝α(r)(β1(r)φ1(r)+β2(r)φ2(r))
[0041] 其中，r是低散射特性组织与高散射特性组织相交的边界上的任意一点；
[0042] Φ(r)是采用扩散方程计算得到的点r处的光通量；
[0043] φ1(r)和φ2(r)是采用三阶简化球谐波方程计算得到的点r处的光通量分量；
[0044] α(r)是扩散方程的边界折射率不匹配系数；
[0045] β1(r)和β2(r)是三阶简化球谐波方程的边界折射率不匹配系数。
[0046] 根据优选实施例，应用下式耦合散射特性组织与空腔的光传输方程：
[0047]
[0048] 式中，r是散射特性组织与空腔相交的边界上的任意一点；
[0049] J(r)是点r处、指向空腔内部的光流率；
[0050] σ是散射特性组织的标识因子，若低散射特性组织与空腔相交，则σ的值取0；若高散射特性组织与空腔相交，则σ的值取1。
[0051] 根据本发明的另一方面，提供了一种内窥式X射线发光断层成像方法，包括如下步骤：(a)采集受检对象内部的探针受激后发射的近红外光信号图像，并采集受检对象的X射线图像；(b)对采集的近红外光信号图像和X射线图像进行图像预处理；(c)对预处理后的X射线数据进行稀疏重建，以获得受检对象的结构信息和内窥探头在受检对象内部的位置信息；(d)利用获取的结构信息和位置信息，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部的靶向目标的位置信息和分布信息。
[0052] 有利地，根据本发明的内窥式X射线断层成像装置和方法，其优点在于：
[0053] 第一，由于采用X射线作为外部激发源激发受检对象内部的荧光纳米颗粒探针，克服了现有技术中采用近红外激发光源激发荧光探针带来的强自体荧光干扰的问题，能够获取高质量的近红外光信号图像，从而有效地改善重建图像的质量。
[0054] 第二，由于采用X射线作为外部激发源激发受检对象内部的荧光纳米颗粒探针，同时采用内窥探头作为内部探测器采集荧光纳米颗粒探针受激后发出的近红外光，克服了现有技术中将激发源和探测器同时集成于微型化内窥探头前端而带来的信号干扰和成像不准确问题，能够改善激发模式下的内窥式光学断层成像技术的信号干扰和成像准确性。
[0055] 第三，由于采用光传输混合模型作为内窥式X射线发光断层成像技术的前向模型，并且同时考虑荧光纳米颗粒探针的稀疏分布特性和内窥测量数据相对于生物体整体区域的稀疏特性来建立稀疏正则化目标函数，克服了现有技术中采用扩散方程和L2范数正则化带来的成像精度和分辨率不足的问题，能够有效地改善内窥式光学断层成像的精度和分辨率。
[0056] 第四，由于采用内窥探头作为内部探测器来采集荧光纳米颗粒探针受激后发出的近红外光，克服了现有X射线发光断层成像技术仅局限于小型动物体和表层器官成像的问题，能够实现对例如大型动物体或人体的受检对象内部的深层靶向目标的准确成像。
附图说明
[0057] 将结合在附图的辅助下对本发明优选实施例进行的说明来对教导的优选实施例和其它改进进行说明。在附图中：
[0058] 图1是根据本发明实施例的内窥式X射线发光断层成像装置的组成框图；
[0059] 图2是根据本发明实施例的内窥式X射线发光断层成像装置的示意图；
[0060] 图3是根据本发明实施例的内窥式X射线发光断层成像装置的近红外光信号采集子模块剖面组成示意图；
[0061] 图4是根据本发明实施例的内窥式X射线发光断层成像装置的内窥探头剖面示意图；
[0062] 图5是根据本发明实施例的内窥式X射线发光断层成像方法的流程图；
[0063] 图6是根据本发明实施例的内窥式X射线发光断层成像方法的内部靶向目标的三维重建子流程图。
具体实施方式
[0064] 下面结合附图来详细说明本发明实施例，所描述的实例仅旨在便于对本发明的理解，而不应理解为对本发明的限制。
[0065] 根据本发明实施例，利用X射线光管作为外部激发源激发受检对象内部的探针，同时利用微型内窥探头作为探测器来采集探针受激后发出的近红外光。然后利用与X射线光管、成像体固定床共线的X射线平板探测器获取受检对象的三维结构成像信息。基于内窥式X射线发光断层成像装置及其成像过程，建立内窥式X射线发光断层成像的物理模型和前向光传输模型；在该模型基础上，利用有限元方法，同时针对受检对象内部探针的稀疏分布特性和测量数据的严重不足性，构造稀疏正则化目标函数；采用合适的优化算法求解目标函数，实现对受检对象内部器官病变或深层靶标的精确成像。
[0066] 下面结合图1、图2、图3和图4对根据本发明实施例的内窥式X射线发光断层成像装置进行详细描述。根据本发明实施例的内窥式X射线发光断层成像装置可以包括：激发源单元1、信号采集单元2、操纵单元3和计算单元4。
[0067] 激发源单元1可以包括X射线光管11和X射线光管电源12，作为外部激发源激发受检对象内部的荧光纳米颗粒探针，同时还用于扫描受检对象以获取其解剖结构信息和内窥探头在体内的位置信息。在实施例中以荧光纳米颗粒探针为例进行说明，本领域技术人员可以理解，探针也可以是例如稀土探针等。
[0068] 信号采集单元2可以包括X射线信号采集模块21和近红外光信号采集模块22。X射线信号采集模块21可以包括X射线CMOS平板探测器211，用于采集受检对象的X射线投影图像；近红外光信号采集模块22可以包括气球式成像导管221、内窥探头222、信号传输光纤束223和光电倍增管阵列探测器224，用于内窥式地采集荧光纳米颗粒探针受激后发出的近红外光信号图像。
[0069] 操纵单元3可以包括C型旋转臂模块31和成像体固定床模块32。该C型旋转臂模块31可以包括C型臂311、旋转控制装置312、C型臂支撑柱313和步进电机314，用于多角度地采集受检对象的X射线投影图像；成像体固定床模块32包括成像体固定床321、平动滑轨322、升降台323和步进电机324，用于将受检对象固定于成像体固定床上，并在X射线成像视场内对受检对象的位置进行调整。
[0070] 计算单元4可以包括系统操作控制模块41、图像预处理模块42以及图像重建模块43。
[0071] 例如，上述各部件的结构关系如下所述：
[0072] 将X射线光管11和X射线CMOS平板探测器211固定于C型臂311的两端，相对放置。X射线光管11通过电源线连接于X射线光管电源12，另外通过控制线连接于计算单元4，用于控制X射线光管11发射X射线的开关和信号强弱；X射线CMOS平板探测器211通过控制线和数据线与计算单元4相连，用于控制X射线CMOS平板探测器211对X射线信号的采集以及采集的信号到计算单元4的数据传输。C型臂311与旋转控制装置312相连，旋转控制装置312固定于C型臂支撑柱313上。此外，旋转控制装置312通过控制线与步进电机314相连，步进电机314通过控制线与计算单元4相连；通过计算单元4的系统操作控制模块41控制步进电机314，使C型臂311在旋转控制装置312上例如等弧度间隔滑动。
[0073] 成像体固定床321固定于平动滑轨322上，而平动滑轨322通过控制线与步进电机324相连。此外，平动滑轨322固定于升降台323上，而升降台323亦通过控制线与步进电机
324相连；步进电机324通过控制线与计算单元4相连。
[0074] 操作中，将受检对象固定于成像体固定床321。之后，通过计算单元4的系统操作控制模块41控制步进电机324，调节升降台323，使成像体固定床321沿着垂直方向升降；同时调节平动滑轨322，使成像体固定床321沿水平方向滑动，最终使成像体固定床321与X射线光管11、X射线CMOS平板探测器211共线，并保证C型臂311等弧度滑动时，成像体固定床321始终处于X射线光管11辐射和X射线CMOS平板探测器211成像范围之内。
[0075] 在近红外光信号采集模块22中，气球式成像导管221由医用导管组成，一端连接于信号传输光纤束223中的空气输送管2232，另一端固定于信号传输光纤束223的外表面，并将内窥探头222包裹其中，用于在采集近红外光信号图像时被充入空气以撑起空腔类器官的腔壁。内窥探头222可以包括N组鼓形透镜2221和棒透镜2222对(例如，一般不少于60个，本例中取60组)，每对鼓形透镜2221和棒透镜2222均与信号传输光纤束223中的信号传输光纤2231相连，用于采集荧光纳米颗粒探针受激后发出的近红外光信号图像。
[0076] 近红外光信号经鼓形透镜2221到达棒透镜2222，再经棒透镜2222作用后耦合到信号传输光纤2231上。信号传输光纤束223采用双层环状阶梯模式，外层由等间隔均匀分布的N根信号传输光纤2231(例如，一般不少于60个，本例中取60根)构成，内层由一根空气输送管2232构成，内层略长于外层2厘米，并且两层之间和最外层由抗干扰材料进行包裹；其中，信号传输光纤2231的一端连接于内窥探头222中的鼓形透镜2221和棒透镜2222，另一端连接于光电倍增管阵列探测器224，用于向光电倍增管阵列探测器224传送内窥探头222采集的近红外光信号；空气输送管2232的一端连接于气球式成像导管221，另一端连接于空气充气机2233，并在长于外层信号传输光纤2231的部分开个圆形小孔，用于向气球式成像导管221中充入空气。光电倍增管阵列探测器224由方阵排列的一组(一般不少于60个，本例中取
60)光电倍增管组成，并通过控制线与光电倍增管阵列探测器电源箱2241相连，同时通过控制线和数据线与计算单元4相连。
[0077] 计算单元4中的系统操作控制模块41用于控制操纵单元3、激发源单元1和信号采集单元2的整体操作，以完成对C型臂311的旋转控制装置312、成像体固定床321、X射线光管11、X射线CMOS平板探测器211和近红外光信号采集模块22的控制，实现对X射线投影图像和近红外光信号图像的采集。
[0078] 图像预处理模块42用于对采集的X射线投影图像和近红外光信号图像进行预处理，例如对采集的X射线投影图像进行坏点补偿、去噪、图像校正等处理，对采集的近红外光信号图像进行坏点补偿、去噪和感兴趣区域提取等处理。
[0079] 图像重建模块43用于对预处理后的X射线投影图像进行稀疏重建，获取受检对象的解剖结构信息和内窥探头的位置信息，并且对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部深层靶向目标的位置信息和分布信息。
[0080] 下面参照图5来描述根据本发明实施例的内窥式X射线断层成像方法。该方法可以包括：
[0081] 步骤S1，多角度数据采集
[0082] (1)利用计算单元4中的系统操作控制模块41调节成像体固定床321，使之位于适合X射线发光断层成像的位置：
[0083] 具体地，首先将受检对象固定于成像体固定床321上，通过计算单元4的系统操作控制模块41控制步进电机324，来调节升降台323，使成像体固定床321沿着高度方向升降至C型臂311的中心高度处；然后，调节平动滑轨322，使成像体固定床321沿水平方向滑动，保证成像体固定床321与X射线光管11、X射线CMOS平板探测器211共线放置；最后，微调平动滑轨322和升降台323，保证C型臂311等弧度滑动时，成像体固定床321始终处于X射线光管11辐射和X射线CMOS平板探测器211最佳成像范围之内。
[0084] (2)利用计算单元4中的系统操作控制模块41采集多个角度的X射线投影图像和近红外光信号图像：
[0085] 首先，通过计算单元4的系统操作控制模块41控制步进电机314，使C型臂311在旋转控制装置312上等间隔小角度滑动(一般不大于3°，本例中取3°)。本领域技术人员可以理解，也可以按照可变间隔来采集图像。
[0086] 然后，在每个角度上，通过计算单元4的系统操作控制模块41开启X射线光管11，照射受检对象；最后，利用信号采集模块2中的近红外光信号采集子模块22采集受检对象内部的荧光纳米颗粒探针受激后发射的近红外光信号图像，同时利用X射线信号采集子模块21采集X射线投影图像。
[0087] 步骤S2数据预处理
[0088] (S21)利用计算单元4中的图像预处理模块42对采集的多角度近红外光信号图像进行预处理，例如包括但不限于噪声去除、感兴趣区域提取以及坏点补偿等；
[0089] (S22)利用计算单元4中的图像预处理模块42对采集的多角度X射线投影数据进行预处理，例如包括但不限于噪声去除、坏点补偿以及图像校正等。
[0090] 步骤S3解剖结构信息的获取
[0091] 利用计算单元4中的图像重建模块43对预处理后的X射线投影数据进行稀疏重建，获得受检对象的三维体数据和内窥探头在体内的位置信息；然后，利用例如3DMED软件中的人机交互式分割方法对获得的三维体数据进行器官分割，获取其解剖结构信息。本领域技术人员可以理解，也可以通过其他方法来通过三维体数据得到受检对象的解剖结构信息。
[0092] 步骤S4内部靶向目标的三维重建
[0093] 利用图像重建模块43，结合预处理后的X射线投影图像，对预处理后的近红外光信号图像进行光学三维重建，获取受检对象内部靶向目标的精确位置和分布信息。
[0094] 参照图6，本步骤可以具体实现如下：
[0095] (S41)前向光传输混合模型的构建
[0096] 基于获取的受检对象的解剖结构信息、各组织的光学特性参数以及内窥探头222在体内的位置信息，采用折射率边界耦合条件，构建基于混合三阶简化球谐波方程-扩散方程-自由空间光传输方程的内窥式X射线发光断层成像的前向光传输混合模型。
[0097] 具体地，基于步骤S3获取的受检对象的解剖结构信息以及内窥探头在受检对象内部的位置信息，结合组织的光学特性参数，将生物组织划分为低散射特性组织、高散射特性组织、空腔三大类。划分过程可以按如下两步规则进行：
[0098] 首先，根据组织的特定位置进行粗划分。将荧光纳米颗粒探针所在组织和体表所属组织划分为低散射特性组织，将近红外光信号图像采集时内窥探头所处的组织划分为空腔；
[0099] 第二步，根据组织的光学参数特性进行细划分。应用式(1)和下面的准则将生物组织划分为不同散射特性组织：
[0100]
[0101] 式中， 是不同散射特性组织的划分标准因子，μ′s是生物组织的约化散射系数，μa是生物组织的吸收系数。可以按照如下准则进行划分：如果 则将生物组织划分为高散射特性组织；如果 则将生物组织划分为低散射特性组织。
[0102] 然后，分别采用三阶简化球谐波方程、扩散方程和自由空间光传输方程来描述近红外光在上述步骤中划分的三大类组织中的传输过程。
[0103] 该描述过程可以按如下规则进行：采用三阶简化球谐波方程描述近红外光在低散射特性组织中的传输过程，采用扩散方程描述近红外光在高散射特性组织中的传输过程，采用自由空间光传输方程描述近红外光从空腔与散射特性组织之间的边界到内窥探头的传输过程。具体的描述方法是本领域公知的，此处不再赘述。
[0104] 然后，构造边界耦合条件以耦合(couple)这三个光传输方程，同时结合X射线激发荧光纳米颗粒探针过程，构造内窥式X射线发光断层成像的前向光传输混合模型。例如，X射线激发荧光纳米颗粒探针的过程可以用下面的公式(7)来描述，即，可以根据公式(7)建立起荧光纳米颗粒探针浓度、X射线强度和发射的近红外光强度的关系，然后在光传输方程中结合公式(7)。
[0105] 该构造过程可以按如下步骤进行：
[0106] 首先，建立三阶简化球谐波方程与扩散方程之间的耦合关系。由于三阶简化球谐波方程和扩散方程描述的均是散射特性生物组织中的光传输过程，亦即光在耦合边界两侧的生物组织中均以漫射光传输，因此两个光传输方程在耦合边界处获得的出射光流率相同。因此，三阶简化球谐波方程在耦合边界上获得的出射光流率 为：
[0107]
[0108] 式中，r是低散射特性组织与高散射特性组织相交的耦合边界上的任意一点，φ1(r)和φ2(r)是采用三阶简化球谐波方程计算得到的点r处的光通量分量，β1(r)和β2(r)是三阶简化球谐波方程的边界折射率不匹配系数。
[0109] 扩散方程在耦合边界上获得的出射光流率 为：
[0110]
[0111] 式中，Φ(r)是采用扩散方程计算得到的点r处的光通量，α(r)是扩散方程的边界折射率不匹配系数。
[0112] 联立上述两个方程，并令 建立三阶简化球谐波方程与扩散方程之间的物理量在耦合边界上的转换关系，以实现低散射特性组织和高散射特性组织光传输方程的耦合：
[0113] Φ(r)＝α(r)(β1(r)φ1(r)+β2(r)φ2(r)) (4)
[0114] 此外，可以建立三阶简化球谐波方程或扩散方程与自由空间光传输方程之间的耦合关系。
[0115] 在内窥式X射线发光断层成像技术中，近红外光从荧光纳米颗粒探针发出后，先在散射特性生物组织中传输，然后到达散射特性组织与空腔的边界，最后穿过空腔到达内窥探测器。因此，在散射特性组织与空腔的边界两侧，光由漫射光变换为非漫射光。根据自由空间光传输理论，漫射光将在两者边界上形成朗伯源(Lambertian source)，而形成的朗伯源的光流率即为从散射特性组织出射的光流率。因此，可通过下式构建散射特性组织与空腔的光传输方程之间的耦合关系：
[0116]
[0117] 式中，r是散射特性组织与空腔相交的边界上的任意一点；J(r)是点r处、指向空腔内部的光流率；σ是散射特性组织的标识因子，可以按如下规则取值：若低散射特性组织与空腔相交，则其值取0；若高散射特性组织与空腔相交，则其值取1。
[0118] 最后，基于上述建立的边界耦合条件，同时结合X射线激发荧光纳米颗粒探针过程，建立内窥式X射线发光断层成像的前向光传输混合模型：
[0119]
[0120] 式中，ΩLD和 是低散射特性组织及其边界；ΩHD和 是高散射特性组织及其边界；BLH是低散射特性组织与高散射特性组织相交的耦合边界；BV是散射特性组织与空腔相交的耦合边界；R是内窥探头的采集点集合；μa(r)，μa1(r)，μa2(r)，μa3(r)和D(r)是生物组织的光学特性参数；A1，A2，B1，B2，C1，C2，D1和D2分别是生物组织边界折射率相关参数；v是生物组织边界的外法向；μ′a是空腔的吸收系数；P(rd)是内窥探测器上接收的光功率；θr′和分别是点r′，rd处的表面法向量与两者之间的方向向量之间的夹角；ξ(r′，rd)是点r′与点rd之间的可见度因子；S(r)是荧光纳米颗粒探针受激后发出的近红外光形成的近红外光光源的能量密度，其由照射到荧光纳米颗粒探针的X射线强度I(r)、荧光纳米颗粒探针的浓度ρ(r)和荧光纳米颗粒探针向近红外光的转化效率ε共同确定：
[0121] S(r)＝εI(r)ρ(r) (7)
[0122] (S42)系统方程的建立
[0123] 首先，利用有限元方法对步骤(S41)建立的前向光传输混合模型进行离散化，组装每个离散点上的子系统方程，建立第k个X射线照射角度对应的系统方程：
[0124] Akρ＝Pk (8)
[0125] 式中，Ak是第k个X射线照射角度对应的系统矩阵，依赖于生物体内三大类组织的分布、光学特性参数和X射线照射到体表的位置；ρ是荧光纳米颗粒探针的浓度；Pk是第k个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值。
[0126] 然后，通过联立M个X射线照射角度对应的系统方程，建立内窥式X射线发光断层成像的总系统方程：
[0127] Aρ＝P (9)
[0128] 式中，A是总系统矩阵，定义为A＝[A1，A2，…，AM]T；T是矩阵转置操作符；P是与所有M个X射线照射角度对应的、内窥探头上近红外光信号图像的计算值向量，定义为P＝[P1，P2，…，PM]T；M是选取的X射线照射角度的个数。一般选取1～5个X射线照射角度下采集的近红外光信号图像进行光学三维重建。
[0129] (S43)目标函数的建立
[0130] 根据步骤S42获得的近红外光信号图像的计算值向量与针对原始采集的近红外光信号图像得到的测量值向量之间的误差，结合荧光纳米颗粒探针的稀疏分布和非负约束，建立下列稀疏正则化目标函数：
[0131]
[0132] 式中，Pm是与M个X射线照射角度对应的，由内窥探头采集的原始近红外光信号图像的测量值向量；
[0133] λ是稀疏正则化参数；
[0134] ||F||1定义为求解矩阵F的L1范数。
[0135] (S44)目标函数的优化求解
[0136] 选用合适的稀疏正则化问题求解算法对步骤(S43)建立的目标函数进行求解，获得荧光纳米颗粒探针在受检对象内部的空间位置和浓度分布，从而间接获取靶向目标的位置信息和分布信息。本领域技术人员可以理解，这里的靶向目标的位置信息和分布信息通常对应于受检对象的病灶或异常组织的位置信息和分布信息。本领域技术人员还可以理解，分布信息通常提供了靶向目标的强度等信息。
[0137] (S45)三维重建结果显示。
[0138] 对步骤(S44)获得的受检对象内部的靶向目标的重建结果和步骤S3获取的解剖结构信息进行图像融合，获得重建的靶向目标在受检对象中的三维图像，并显示重建的靶向目标的三维图像。
[0139] 本领域技术人员将受益于前述说明书和相关附图中呈现的教导而想到本文所述的本发明的多种改型和其它实施例。因此，应该理解本发明不限于所公开的具体实施例，并且应该理解，本发明旨在将变型和其它实施例包括在由所附权利要求限定的范围内。尽管本文采用了具体术语，但是只是从描述的角度来使用它们，而并非限制。
