发明名称 一种快速超分辨率X射线荧光CT成像及重构系统及方法
摘要 本发明公开了一种快速超分辨率X射线荧光CT成像及重构系统及方法,其采用源模块、成像模块、处理模块以及移动控制平台组成的成像系统探测样本内部所含高原子序数纳米粒子的浓度及分布;成像模块中后准直器阵列与X射线平板探测器相结合,每个微准直器与样本中特定区域内的纳米粒子相对应,一次成像即可获得样本内部纳米粒子的浓度及分布情况;并基于平板探测器与微准直器在位置以及分辨率方面的关系,进一步微调微准直器阵列的位置后再次成像,多次重复,通过方程求解以实现类似图像插值的结果,提高图像重构质量。
申请公布号 CN103767726B 申请公布日期 2016.01.06
申请号 CN201410064528.5 申请日期 2014.02.25
申请人 重庆大学 发明人 冯鹏;魏彪;何鹏;陈绵毅
分类号 A61B6/03(2006.01)I;G01N23/223(2006.01)I 主分类号 A61B6/03(2006.01)I
代理机构 重庆华科专利事务所 50123 代理人 康海燕
主权项 一种快速超分辨率X射线荧光CT成像及重构系统,其特征在于:所述系统包括源模块、成像模块、处理模块以及移动控制平台;源模块与成像模块之间为放置待测样本的位置;所述源模块由X射线源、滤波器和前准直器组成,X射线源产生的X射线经滤波器滤波、前准直器准直后入射到待测样本,与待测样本中的纳米粒子发生作用,激发后者释放特征荧光X射线;所述成像模块由后准直器与X射线平板探测器构成,荧光X射线经后准直器准直之后被平板探测器探测后形成数据流,传输至高性能计算机进行处理和超分辨率重构;所述后准直器由多个微准直器组成阵列,并与平板探测器级联在一起;所述平板探测器、后准直器均与样本垂直,位于样本正上方,进入每个微准直器的特征荧光X射线恰好与样本内部纳米粒子的分布相对应;所述处理模块采用高性能计算机,对所述数据流进行处理和超分辨率重构;所述移动控制平台是一套通过计算机、运动控制机构实现对样品和后准直器、平板探测器精密移动进行控制的平台,其与样品台以及后准直器、平板探测器相连;所述系统进行快速X射线荧光CT成像及重构的方法步骤如下:1)将样本放置于源模块与成像模块之间,将成像模块与处理模块连接,将控制平台与样品以及成像模块的后准直器阵列、平板探测器相连;其中源模块内部部件保持中心线重合,样本与成像模块保持中心线重合,且样本位于成像模块的正下方;2)启动源模块,使X射线源所发出的X射线经滤波器滤波、前准直器准直后形成扇形束,该扇形束面与样本的中心线垂直,并直接穿过整个样品,与样本内部的纳米粒子相互作用,释放出特征荧光X射线;3)沿样本正上方传播的特征荧光X射线在穿过样本后,被后准直器阵列所准直,并与被X射线平板探测器所探测,探测器每个探测单元的输出与进入该单元的X射线数目成正比;4)探测器依据入射X射线能量的不同,按所属能量区域分别进行记录,设探测器所能探测的X射线能谱范围为E<sub>1</sub>~E<sub>2</sub>(keV),有P×P个探测单元,能谱细分为K段,每段能谱间隔为<img file="FDA0000782447000000021.GIF" wi="465" he="118" />相应会产生K幅P×P大小的图像数据I<sub>k</sub>(m,n):<img file="FDA0000782447000000022.GIF" wi="1174" he="92" />设所探测的特征荧光X射线能量为E<sub>F</sub>(keV),位于第Q个能量段,即(Q‑1)ΔE≤E<sub>F</sub>≤E<sub>F</sub>+QΔE,平板探测器以第Q幅图像数据I<sub>Q</sub>(m,n)作为最终输出图像;5)超分辨率重构:基于平板探测器与微准直器在位置以及分辨率方面的关系,利用移动控制平台微调微准直器阵列的位置,后再次成像,多次重复,然后通过方程组求解,具体地:设微准直器的尺寸为L<sub>s</sub>(mm)×L<sub>s</sub>(mm),平板探测器探测单元尺寸为R(mm)×R(mm),单元个数为P×P,L<sub>s</sub><R,则微型准直器阵列沿先行后列方向移动T<sup>2</sup>‑1次,移动的步长恰好为像元的边长,每移动一次后成一次像,设每次成像后输出图像为:<img file="FDA0000782447000000031.GIF" wi="1268" he="84" />其中T=floor(R/L<sub>s</sub>)重构高分辨率图像为:<img file="FDA0000782447000000032.GIF" wi="891" he="86" />通过下式可进行重构:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>I</mi><mrow><mi>Re</mi><mi>c</mi></mrow></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>=</mo><msup><mi>I</mi><mn>1</mn></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>I</mi><mrow><mi>Re</mi><mi>c</mi></mrow></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>=</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>I</mi><mrow><mi>Re</mi><mi>c</mi></mrow></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>=</mo><msup><mi>I</mi><mi>T</mi></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000782447000000033.GIF" wi="958" he="786" /></maths>对上述方程组进行求解,所得I<sup>Rec</sup>(m,n),1≤m,n≤P即为高分辨率重构图像。
地址 400030 重庆市沙坪坝区沙正街174号