CN109589126A - 基于宽光束小步长扫描方式的x射线发光断层成像方法 - Google Patents
基于宽光束小步长扫描方式的x射线发光断层成像方法 Download PDFInfo
- Publication number
- CN109589126A CN109589126A CN201811214285.3A CN201811214285A CN109589126A CN 109589126 A CN109589126 A CN 109589126A CN 201811214285 A CN201811214285 A CN 201811214285A CN 109589126 A CN109589126 A CN 109589126A
- Authority
- CN
- China
- Prior art keywords
- ray
- imaging
- scanning
- beams
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 67
- 238000004875 x-ray luminescence Methods 0.000 title claims abstract description 8
- 230000005284 excitation Effects 0.000 claims abstract description 10
- 238000010276 construction Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 36
- 239000011159 matrix material Substances 0.000 claims description 17
- 239000013307 optical fiber Substances 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000003325 tomography Methods 0.000 claims description 11
- 239000000835 fiber Substances 0.000 claims description 10
- 230000004907 flux Effects 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 9
- 238000013519 translation Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 239000000523 sample Substances 0.000 abstract 1
- 238000002474 experimental method Methods 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004904 shortening Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000004613 tight binding model Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 241001522296 Erithacus rubecula Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 235000013405 beer Nutrition 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000003068 molecular probe Substances 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种基于宽光束小步长扫描方式的X射线发光断层成像方法,包括下列步骤:①数字X射线成像系统和XLCT成像参数设置。③将注入纳米荧光颗粒的待测物体放到成像腔。④X射线束的设置,实现宽光束小步长的扫描激发方式。⑤对待测物体进行扫描,而后通过光纤探测器接收近红外光子,用光电倍增管记录光子数量。⑥正向模型构建。⑦逆向模型构建。
Description
技术领域
本发明属于生物医学工程及医学影像学领域,涉及一种基于宽光束小步长扫描方式的X射线发光断层成像新方法。
背景技术
光学分子影像作为分子成像技术的重要分支,可以对活体状态下的生物体进行细胞和分子水平上的定性和定量成像,该技术利用特定的分子探针跟踪生物体的病理变化,在体外采用光学分子成像系统采集标记物产生的光信号,这种光信号携带了大量的信息,根据成像模型可以对获得的信号进行三维重建,完成特异性成像,这对于早期疾病检测、治疗、药物研发以及药效评估都有着重大的研究意义。目前,光学分子影像方法主要包括生物断层发光成像(Bioluminescence Tomography,BLT)、荧光分子断层成像(FluorescenceMolecular Tomography,FMT)、扩散光学断层成像(Diffuse Optical Tomography,DOT)和X射线发光断层成像(X-ray Luminescence Computed Tomography,XLCT)等技术。而XLCT作为一种新型的光学分子成像技术,因X射线穿透能力具有更深的成像深度,可以实现深部组织测量且具有较高成像分辨率,因而受到广泛的关注。XLCT对于肿瘤的诊断和治疗有着突破性的意义,它可以利用纳米荧光分子对肿瘤细胞进行标定,通过X射线照射激发纳米荧光分子发射近红外光,然后近红外光穿透生物组织后被探测器接收到,根据接收到的携带位置和浓度信息的光信号进行建模与重建,从而反演出纳米荧光分子在生物体内的三维分布,确定肿瘤的位置和大小,该技术为肿瘤诊治提供一种新型有效的分子影像技术,具有十分重要的意义。
目前国内外在XLCT技术研发方面还在初级阶段,按照X射线激发模式划分,可以分为锥束和窄束XLCT这两类系统。锥束XLCT成像时间短,但是具有较为严重的欠定性,往往会导致重建图像的质量差,而窄束XLCT在成像分辨率方面具备优势,但由于其成像时间长,阻碍了该方法的临床推广应用。针对窄束XLCT成像时间长的问题,国外学者提出一种多光束XLCT成像的方法。该方法通过采用多个光束进行扫描,在牺牲一些分辨率的前提下,成倍缩短了成像时间,该方法的提出推动了窄束XLCT的进一步发展。例如采用四条X射线窄束对仿体同时扫描,可以使每条光束扫描的距离降低到单光束扫描的四分之一,因此扫描时间也可以降低到单光束扫描的四分之一。多光束的XLCT成像方法虽然成倍缩短了成像时间,但与锥束XLCT相比,由于仍采用X射线窄束作为激发光,在有限时间内用于激发纳米荧光分子的X射线光子数有限,限制了激发产生近红外光的光子数,由此导致所采集的近红外信号强度有限,易受噪声干扰的问题。窄束XLCT具有高空间分辨率的特点,但也面临着成像时间长的弊端,制约了该技术在临床应用上的价值,因此缩短成像时间并保证成像质量已经成为该研究面临的一个关键问题。多窄束XLCT能够缩短成像时间,但是随着X线束数的增加,成像模型的信息完备性降低,导致重建图像的分辨率也随之大大下降。
发明内容
本发明提出一种基于宽光束小步长扫描方式的X射线发光断层成像新方法,可以在缩短成像时间的同时保证成像质量,其中,采用宽光束增强了受激产生的信号强度,由此缩短了成像时间,小步长扫描增强了信息完备性,解决宽光束带来的分辨率下降问题。本发明的技术方案如下:
一种基于宽光束小步长扫描方式的X射线发光断层成像方法,包括下列步骤:
①数字X射线成像系统和XLCT成像参数设置:将系统开启,让X微焦点光源和CMOS探测器预热,X射线光源位于物体的正上方且到物体的距离为80cm,利用光纤接收从待测物体表面发出的近红外光。
③将注入纳米荧光颗粒的待测物体放到成像腔,调整位移台进行校准操作,保证待测物体在实验过程中可以全部被扫描到。
④X射线束的设置:根据待测物体尺寸设计合适的准直器实现多光束X射线准直,设采用k条光束,光束初始位置沿成像腔直径方向将成像腔k等分,射线宽度为0.8mm,平移步长为0.1mm,即实现宽光束小步长的扫描激发方式。
⑤对待测物体进行扫描,而后通过光纤探测器接收近红外光子,用光电倍增管记录光子数量,即测量量Φ,该角度扫描完成后,将线性平移台调整到初始位置,依次旋转60°,120°,再次按上述步骤进行扫描,并将采集到的信号保存,用于后续图像重建。
⑥正向模型构建:利用有限元法将待测物体进行有限元剖分建模,剖分单元采用四面体,将信号采集光纤的位置配准到该模型,将光纤探测器所在位置设置点光源,代入扩散方程后求解每个光纤探测器对应的系统矩阵Φi,i=1,2,…,nd,nd,为探测器数量,则获得整体系统矩阵A=[Φn1,Φn2,...,Φnd]Τ,然后根据选择性激发向量对整体系统矩阵A进行修正,选择性激发向量定义为:位置r处的节点noder在X射线束照射范围ΩX内,其值设置为1,否则设置为0,即得到I是每一个投影角度下线性扫描步长的数量,J是投影角度的个数,I×J为总的采集数据的次数,采用k条宽光束进行小步长扫描,扫描I*J次,得到待测物体表面的光通量密度的测量量b。
⑦逆向模型构建:对纳米发光粒子密度的空间分布η的求解作无约束最优化求解b是待测物体表面的光通量密度的测量量,采用ART迭代算法进行迭代处理,以求得η的最小范数-最小二乘解β为松弛因子,取β=0.2。
本发明通过增大光束宽度,增强了X射线的强度,从而提高了X射线剂量的利用率,相比于窄束XLCT更利于成像而且缩短了成像时间;同时,通过小步长的方法提高了重建图像的质量,达到了窄束XLCT重建图像的效果。
附图说明
图1多窄束XLCT成像系统结构图
图2三个光束扫描模式
图3 delta=1.2mm的情况下,(a-e)步长分别为:0.5mm,0.4mm,0.3mm,0.2mm,0.1mm
图4 delta=1.0mm的情况下,(a-e)步长分别为:0.5mm,0.4mm,0.3mm,0.2mm,0.1mm.
图5 delta=0.8mm的情况下,(a-e)步长分别为:0.5mm,0.4mm,0.3mm,0.2mm,0.1mm.
图6 delta=1.2,1.0,0.8mm下三个评价指标的折线图(a)DICE(b)CNR(c)RMSE
图7(a1-c1)为窄光束重建的图像,步长分别为0.3,0.2,0.1,单位是mm;
(a2-c2)为宽光束小步长重建的图像,步长同样分别为0.3,0.2,0.1,单位是mm
具体实施方式
下面从几个方面介绍本发明。
1多窄束XLCT成像系统
实验成像系统为Pixarray 100小动物数字放射成像系统,由美国BIOPTICS公司制造。该系统的探测器为1024×1024的CCD阵列,像素大小为50μm×50μm,14级灰度。横向及纵向的空间分辨率均为每毫米20像素。X射线管的焦斑尺寸为50μm。探测器点扩散函数的半高宽为110μm。实验中,X射线源的工作电压为40kV,工作电流为0.5mA。成像物体模拟小鼠的生理参数设计而成。实验设置X射线源位于物体的正上方且到物体的距离为80cm,半圆弧状的光纤架位于物体的正下方,光纤用铅皮包裹,而后固定在光纤架上,用于接收从待测物体表面发出的光流量。待测物体放于直径为30mm,高为50mm的成像腔内,成像腔固定在平移台上。通过控制器可以实现成像腔的线性平移和旋转,成像系统结构图如图1所示。
2宽光束小步长扫描方式的XLCT成像原理
在XLCT成像过程中,X射线源发出的X射线穿透生物组织,照射到纳米发光材料,激发产生近红外光,X射线的强度在整个过程中会发生变化,根据朗伯比尔定律,可得到在生物组织中X射线的强度变化:
其中,X(r0)是X射线在初始位置r0处的强度,μ(τ)是X射线在位置τ处的衰减系数。
在成像过程中,纳米发光材料经过X射线照射后,受激发出近红外光,激发出的光源强度可表示为:
S(r)=εX(r)η(r). (2)
其中,S(r)是近红外光强度,η(r)是纳米荧光粒子的密度,ε是光子产额。
光子在生物软组织中的传播过程可以用辐射传输方程(RTE)来描述,由于生物软组织对近红外光具有高散射性和低吸收性,RTE一般被近似为扩散方程,因此成像模型可以表达为:
-▽·[D(r)▽Φ(r)]+μa(r)Φ(r)=S(r) r∈Ω (3)
其中,D(r)=(3(μa(r)+(1-g)μs(r)))-1是扩散系数,μa(r)是吸收系数,g是各向异性参数,μs(r)是散射系数,Φ(r)是光通量密度。
采用稳态扩散方程描述光子传播,在边界处可由罗宾边界条件进行补充:
其中,k(r,n,n′)是边界不匹配因子,它取决于Ω内的折射率n和周围环境的折射率n′,v(r)是外向单位法向向量。
有限元方法已被广泛应用于求解扩散方程,特别是对任意几何图形的求解。基于有限元理论,将式方程(3)(4)联立,用形状函数ψi和ψj进行离散化,构建矩阵方程:
M·Φ=F·S (5)
其中,
方程(5)中的矩阵M是正定的,因此可以得到:
Φ=M-1·F·S=M-1·F·ε·X·η=A·η (6)
实际应用时,矩阵A的获得依靠求解格林函数:将探测器所在位置设置为点光源位置,代入扩散方程后可以求解出每个点光源照射后的物体内部的光通量密度分布Φi(i=1,2,…,nd,nd是探测器的个数),则A=[Φn1,Φn2,...,Φnd]Τ。由于X射线是有间距的多条宽光束,在一次照射下只能照射到待测物体的部分区域,然后通过窄步长平移多次才能将待测物体全部照射到,因此在每次采集的数据都需要用选择性矩阵来修正系数矩阵A,我们设定选择性激发向量如下:
当位置r处的节点(noder)在X射线束照射范围ΩX内,其值设置为1,否则设置为0。其中i=1,...,I×J为总的采集数据的次数,I是每一个投影角度下线性扫描步长的数量,J是投影角度的个数。修正后的系数矩阵A为:
其中,nd是探测器的个数,Φj(j=n1,...,nd)为第j个探测光纤对应的格林函数。然后采用3条宽光束进行小步长扫描,扫描I*J次,得到待测物体表面的光通量密度的测量量b。
由于上述的离散过程导致了实际测量数据个数远远小于待求的未知量个数,并且方程的解很容易受噪声干扰,因此很难用直接的矩阵求逆的方法得到纳米发光粒子密度的空间分布η。在此情形下,只能通过对原问题的求解过程转化为无约束最优化求解:
正规化方法主要有截断SVD法、Tikhonov-Mill法,迭代正规法(Landweber法)、代数重建技术(algebraic reconstruction technique,ART)等。由于ART是广泛应用于医学成像领域的一种基于矩阵行的迭代重建技术,因此本发明采用该技术进行原问题求解:
式中,j=1,...,m(m为节点个数),β称为松弛因子,其取值范围为[0,2],当β取0.1~0.3时,重建结果比较稳定但收敛速度慢,因此这里取值为0.2。
3基于宽光束小步长扫描方式的X射线发光断层成像新方法技术流程
本发明的基于宽光束小步长扫描方式的X射线发光断层成像新方法的流程描述如下:
1、数字X射线成像系统设置:将系统开启,让X微焦点光源和CMOS探测器预热,为保证实验过程中成像系统的稳定性,在实验前至少需要让机器预热30分钟。
2、设置XLCT成像参数:X射线光源位于物体的正上方且到物体的距离为80cm,半圆弧状的光纤架位于物体的正下方,光纤用铅皮包裹,而后固定在光纤架上,用于接收从待测物体表面发出的近红外光。
3、将注入纳米荧光颗粒的待测物体放到成像腔,成像腔的直径为30mm,高为50mm,调整位移台进行校准操作,保证待测物体在实验过程中可以全部被扫描到。
4、X射线束的设置:根据待测物体尺寸设计合适的准直器实现多光束X射线准直,这里采用3条光束,光束初始位置沿成像腔直径方向将成像腔3等分,射线宽度为0.8mm,平移步长为0.1mm,即实现宽光束小步长的扫描激发方式。
5、设置数字放射成像系统的曝光参数,在上述成像设置下,对待测物体进行扫描,而后通过光纤探测器接收近红外光子,用光电倍增管记录光子数量,即测量量Φ。该角度扫描完成后,将线性平移台调整到初始位置,依次旋转60°,120°,再次按上述步骤进行扫描,并将采集到的信号保存,用于后续图像重建。
6、正向模型构建:正向模型构建:利用有限元法将待测物体进行有限元剖分建模,剖分单元采用四面体,将信号采集光纤的位置配准到该模型,将光纤探测器所在位置设置点光源,代入扩散方程后求解每个光纤探测器对应的系统矩阵(或格林函数)Φi(i=1,2,…,nd,nd为探测器数量),则获得整体系统矩阵A=[Φn1,Φn2,...,Φnd]Τ,然后根据选择性激发向量对整体系统矩阵A进行修正,选择性激发向量定义为:位置r处的节点(noder)在X射线束照射范围ΩX内,其值设置为1,否则设置为0,即得到(I是每一个投影角度下线性扫描步长的数量,J是投影角度的个数,实验中取J=3,I×J为总的采集数据的次数),实验中采用3条宽光束进行小步长扫描,扫描I*J次,得到待测物体表面的光通量密度的测量量b。
7、逆向模型构建:由于实际测量数据个数远远小于待求的未知量的个数,因此需要对纳米发光粒子密度的空间分布η的求解作无约束最优化求解(b是待测物体表面的光通量密度的测量量),然后我们采用ART迭代算法进行迭代处理,以求得η的最小范数-最小二乘解(β为松弛因子,本发明取β=0.2)。根据需要,可以从η的三维重建结果中取任意截面的图像,用于分析该平面上纳米发光粒子密度的分布。
该方法与目前的窄光束XLCT相比,缩短了成像时间,增大了X射线的强度,提高了X射线剂量的利用率,而且在重建图像的质量上,达到了窄光束XLCT重建图像的效果。本发明采用三个光束进行照射仿体,如图2所示,三条光束均分整个仿体,在数值仿真和模拟实验中,扫描的距离是单光束的三分之一,因此在扫描时间上成倍的缩短,同时在此基础上增加光束宽度又进一步减少了扫描距离,从而进一步缩短了扫描时间。在实际的仿体实验中,形成窄光束的准直器不易制作,而且会额外增加成本,而使用宽光束,降低了准直器的制作难度,而且得到的数据与仿真的数据差异性小,这样更利于成像,而且缩短了实验扫描时间,这对临床研究提供了可能性,具有重要的研究意义和广泛的应用前景。
本发明采用3个光束小步长扫描方法,扫描模式如图2所示,delta为射线宽度,distance为相邻射线束之间的距离,step为步长,投影角度选0°,60°和120°三个方向,每个投影角度下平移的步长总数为I=2R/step/n,R为仿体直径,n为射线束的个数,目标体objects直径是1.6mm。在图1的扫描模式下,我们分别研究了X射线宽度为delta=1.2mm,delta=1.0mm和delta=0.8mm三种情况下,步长大小对成像质量的影响,步长step设置均依次采用0.5mm,0.4mm,0.3mm,0.2mm,0.1mm这5种情况,重建的图像如图3、4、5所示。
分别记录了图3、4、5不同步长下重建图像的三个评价指标:相似度(DICE),对比度噪声比(CNR)和均方根误差(RMSE),其中,DICE是指建区域中强度值大于最大值的50%的体元的面积与实际目标体的面积之比,一般认为DICE越接近1,越好;CNR是描述重建区域与背景之间的差异的,一般认为CNR越大,重建图像的质量越好;RMSE是描述重建目标体与实际目标体之间的差异的,一般认为RMSE越小,重建图像的质量越好。根据三个指标的数据,绘制了折线图,如图6所示:
从DICE的折线图可以看出,随着步长step的减小,三种射线宽度的DICE均呈上升趋势,表明,重建区域与实际目标体的大小越来越接近;从CNR的折线图可以看出,随着步长step的减小,三种射线宽度的CNR均呈上升趋势,表明,重建区域与背景之间的差异越来越明显;从RMSE的折线图可以看出,随着步长step的减小,三种射线宽度的RMSE均呈下降趋势,表明,重建区域浓度值的大小与实际目标体区域浓度值的大小的差异越来越小。三个评价指标的变化均表明了,在宽光束下,随着步长的减小,重建图像的质量越来越好,尤其是在小步长下(如step=0.1mm),重建图像的质量明显得到提高。
基于宽光束小步长的XLCT,在重建图像的质量上明显得到提高,为了检验重建质量的效果,我们做了三组与窄光束成像的对比实验。在窄光束实验中,束宽delta和步长step相等,分别为0.3mm,0.2mm,0.1mm,其他实验条件与宽光束窄步长实验一致,为了进行很好的对比,我们分别取光束宽度delta为0.8mm情况下,步长step分别为0.3mm,0.2mm,0.1mm的重建图像进行对比,如图7所示:
为了对重建图像的质量进行更好的对比,我们对宽光束和窄光束重建图像的三个评价指标进行了对比,如表1所示:
表1宽光束窄步长与窄光束窄步长评价指标的对比
对比(a1)和(a2),可以发现,(a2)的DICE和CNR均比(a1)的稍微大些,且(a2)的RMSE比(a1)小,表明(a2)重建图像的质量比(a1)好;对比(b1)和(b2),可以发现,(b2)与(b2)的DICE相差很小,但(b2)的CNR比(b1)稍微小一些,而(b2)的RMSE比(b1)小得多,从整体上来看,(b1)和(b2)重建图像的质量相差不大;对比(c1)和(c2),可以发现(c2)的DICE和CNR均比(c1)大,且(c2)的RMSE比(c1)小,这表明,(c2)重建图像的质量要比(c1)好一些。通过三组实验的对比,我们可以得到以下结论:小步长可以解决宽光束带来的分辨率下降问题,提高重建图像的质量,而且可以用来实现窄光束的效果。
Claims (1)
1.一种基于宽光束小步长扫描方式的X射线发光断层成像方法,包括下列步骤:
①数字X射线成像系统和XLCT成像参数设置:将系统开启,让X微焦点光源和CMOS探测器预热,X射线光源位于物体的正上方且到物体的距离为80cm,利用光纤接收从待测物体表面发出的近红外光。
③将注入纳米荧光颗粒的待测物体放到成像腔,调整位移台进行校准操作,保证待测物体在实验过程中可以全部被扫描到;
④X射线束的设置:根据待测物体尺寸设计合适的准直器实现多光束X射线准直,设采用k条光束,光束初始位置沿成像腔直径方向将成像腔k等分,射线宽度为0.8mm,平移步长为0.1mm,即实现宽光束小步长的扫描激发方式;
⑤对待测物体进行扫描,而后通过光纤探测器接收近红外光子,用光电倍增管记录光子数量,即测量量Φ,该角度扫描完成后,将线性平移台调整到初始位置,依次旋转60°,120°,再次按上述步骤进行扫描,并将采集到的信号保存,用于后续图像重建;
⑥正向模型构建:利用有限元法将待测物体进行有限元剖分建模,剖分单元采用四面体,将信号采集光纤的位置配准到该模型,将光纤探测器所在位置设置点光源,代入扩散方程后求解每个光纤探测器对应的系统矩阵Φi,i=1,2,…,nd,nd,为探测器数量,则获得整体系统矩阵A=[Φn1,Φn2,...,Φnd]Τ,然后根据选择性激发向量对整体系统矩阵A进行修正,选择性激发向量定义为:位置r处的节点noder在X射线束照射范围ΩX内,其值设置为1,否则设置为0,即得到I是每一个投影角度下线性扫描步长的数量,J是投影角度的个数,I×J为总的采集数据的次数,采用k条宽光束进行小步长扫描,扫描I*J次,得到待测物体表面的光通量密度的测量量b;
⑦逆向模型构建:对纳米发光粒子密度的空间分布η的求解作无约束最优化求解b是待测物体表面的光通量密度的测量量,采用ART迭代算法进行迭代处理,以求得η的最小范数-最小二乘解β为松弛因子,取β=0.2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811214285.3A CN109589126B (zh) | 2018-10-18 | 2018-10-18 | 基于宽光束小步长扫描方式的x射线发光断层成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811214285.3A CN109589126B (zh) | 2018-10-18 | 2018-10-18 | 基于宽光束小步长扫描方式的x射线发光断层成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109589126A true CN109589126A (zh) | 2019-04-09 |
CN109589126B CN109589126B (zh) | 2022-08-19 |
Family
ID=65957865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811214285.3A Active CN109589126B (zh) | 2018-10-18 | 2018-10-18 | 基于宽光束小步长扫描方式的x射线发光断层成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109589126B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111951400A (zh) * | 2020-08-06 | 2020-11-17 | 天津大学 | 基于U-net网络的窄束X射线激发发光断层成像方法 |
CN113057582A (zh) * | 2021-03-05 | 2021-07-02 | 中国人民解放军空军军医大学 | 基于mlem的自适应fista初始图像的定量锥束x射线发光计算机断层扫描方法 |
CN113288188A (zh) * | 2021-05-17 | 2021-08-24 | 天津大学 | 基于分组注意力残差网络的锥束x射线发光断层成像方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5986278A (en) * | 1997-06-25 | 1999-11-16 | Quantum Imaging Corporation | Filmless photon imaging apparatus |
WO2000067062A1 (en) * | 1999-05-05 | 2000-11-09 | Genomic Solutions Inc. | Laser-based optical scanning system |
US20040131983A1 (en) * | 2002-12-19 | 2004-07-08 | Tomoyuki Ohzeki | Photothermographic material and image forming method using same |
US20060064000A1 (en) * | 2004-09-21 | 2006-03-23 | Vizard Douglas L | Apparatus and method for multi-modal imaging |
US20090324048A1 (en) * | 2005-09-08 | 2009-12-31 | Leevy Warren M | Method and apparatus for multi-modal imaging |
US20110251484A1 (en) * | 2010-04-09 | 2011-10-13 | Carpenter Colin M | Molecular imaging using radioluminescent nanoparticles |
CN103239255A (zh) * | 2013-05-20 | 2013-08-14 | 西安电子科技大学 | 一种锥束x射线发光断层成像方法 |
US20150157286A1 (en) * | 2013-12-06 | 2015-06-11 | Rensselaer Polytechnic Institute | Stored luminescence computed tomography |
CN105640582A (zh) * | 2016-03-02 | 2016-06-08 | 中国人民解放军第四军医大学 | 一种深部组织x射线激发多光谱断层成像系统及方法 |
CN106725565A (zh) * | 2016-11-18 | 2017-05-31 | 天津大学 | 一种稀疏投影下的锥束xct成像质量评估方法 |
-
2018
- 2018-10-18 CN CN201811214285.3A patent/CN109589126B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5986278A (en) * | 1997-06-25 | 1999-11-16 | Quantum Imaging Corporation | Filmless photon imaging apparatus |
WO2000067062A1 (en) * | 1999-05-05 | 2000-11-09 | Genomic Solutions Inc. | Laser-based optical scanning system |
US20040131983A1 (en) * | 2002-12-19 | 2004-07-08 | Tomoyuki Ohzeki | Photothermographic material and image forming method using same |
US20060064000A1 (en) * | 2004-09-21 | 2006-03-23 | Vizard Douglas L | Apparatus and method for multi-modal imaging |
US20090324048A1 (en) * | 2005-09-08 | 2009-12-31 | Leevy Warren M | Method and apparatus for multi-modal imaging |
US20110251484A1 (en) * | 2010-04-09 | 2011-10-13 | Carpenter Colin M | Molecular imaging using radioluminescent nanoparticles |
CN103239255A (zh) * | 2013-05-20 | 2013-08-14 | 西安电子科技大学 | 一种锥束x射线发光断层成像方法 |
US20150157286A1 (en) * | 2013-12-06 | 2015-06-11 | Rensselaer Polytechnic Institute | Stored luminescence computed tomography |
CN105640582A (zh) * | 2016-03-02 | 2016-06-08 | 中国人民解放军第四军医大学 | 一种深部组织x射线激发多光谱断层成像系统及方法 |
CN106725565A (zh) * | 2016-11-18 | 2017-05-31 | 天津大学 | 一种稀疏投影下的锥束xct成像质量评估方法 |
Non-Patent Citations (10)
Title |
---|
YUEMING ZHANG等: "Method for improving the spatial resolution of narrow x-ray beam-based x-ray luminescence computed tomography imaging", 《JOURNAL OF BIOMEDICAL OPTICS》 * |
ZHANG WEI等: "Multiple pinhole collimator based X-ray luminescence computed tomography", 《BIOMEDICAL OPTICS EXPRESS》 * |
ZHANG WEI等: "X-ray luminescence computed tomography:X-ray luminescence computed tomography:an emerging molecular imaging modality", 《INSTRUMENTATION》 * |
张月明: "提高X射线发光断层成像系统性能的优化设计与实验研究", 《中国优秀硕士学位论文全文数据库》 * |
张翔: "X光分幅相机磁透镜结构对电子束成像质量的影响", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
徐晶: "X射线荧光光谱分析法在钨矿检测中的应用与研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
钞田田: "窄束X射线激发发光成像系统研制和重建方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
陈冬梅: "X射线发光断层成像研究", 《中国博士学位论文全文数据库信息科技辑》 * |
高鹏: "高精度锥束X射线发光断层成像关键问题研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
黄盛阳等: "X线激发荧光断层成像的研究进展", 《医疗卫生装备》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111951400A (zh) * | 2020-08-06 | 2020-11-17 | 天津大学 | 基于U-net网络的窄束X射线激发发光断层成像方法 |
CN113057582A (zh) * | 2021-03-05 | 2021-07-02 | 中国人民解放军空军军医大学 | 基于mlem的自适应fista初始图像的定量锥束x射线发光计算机断层扫描方法 |
CN113288188A (zh) * | 2021-05-17 | 2021-08-24 | 天津大学 | 基于分组注意力残差网络的锥束x射线发光断层成像方法 |
CN113288188B (zh) * | 2021-05-17 | 2023-10-10 | 天津大学 | 基于分组注意力残差网络的锥束x射线发光断层成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109589126B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pogue et al. | Maps of in vivo oxygen pressure with submillimetre resolution and nanomolar sensitivity enabled by Cherenkov-excited luminescence scanned imaging | |
Chen et al. | Cone beam x‐ray luminescence computed tomography: A feasibility study | |
US10285659B2 (en) | Stored luminescence computed tomography | |
Graves et al. | A submillimeter resolution fluorescence molecular imaging system for small animal imaging | |
CN108577876B (zh) | 一种多边形静止ct及其工作方法 | |
JP7440406B2 (ja) | 蛍光x線測定のための方法および測定装置 | |
US8676302B2 (en) | Systems and methods for multi-spectral bioluminescence tomography | |
CN109589126B (zh) | 基于宽光束小步长扫描方式的x射线发光断层成像方法 | |
CN111915733B (zh) | 基于LeNet网络的三维锥束X射线发光断层成像方法 | |
CN1278153A (zh) | 检测组织形态的方法 | |
CN104013387B (zh) | 一种太赫兹快速断层成像系统及方法 | |
Zhang et al. | Assessing indocyanine green pharmacokinetics in mouse liver with a dynamic diffuse fluorescence tomography system | |
CN109646032B (zh) | 一种基于加权调制的多光束x射线激发发光断层成像方法 | |
CN109709127A (zh) | 低散射x射线荧光ct成像方法 | |
Shi et al. | Reconstruction of x-ray fluorescence computed tomography from sparse-view projections via L1-norm regularized EM algorithm | |
Cong et al. | X-ray fan-beam luminescence tomography | |
CN113288188B (zh) | 基于分组注意力残差网络的锥束x射线发光断层成像方法 | |
CN103750824B (zh) | 一种针对小动物荧光层析成像系统的信息提取方法 | |
JP4567318B2 (ja) | 生体組織部分内の領域の位置測定方法および装置 | |
Jakub et al. | Studies of J-PET detector to monitor range uncertainty in proton therapy | |
Stuker | Hybrid imaging: combining fluorescence molecular tomography with magnetic resonance imaging | |
Feng et al. | X-ray fluorescence microtomography based on polycapillary-focused X-rays from laboratory source | |
CN107374588B (zh) | 一种基于同步聚类的多光源荧光分子断层成像重建方法 | |
CN113057582B (zh) | 基于mlem的自适应fista初始图像的定量锥束x射线发光计算机断层扫描方法 | |
Chen et al. | Cone Beam X‐Ray Luminescence Tomography Imaging Based on KA‐FEM Method for Small Animals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |