本发明涉及医疗器械领域中的一种时间校正装置，更具体地涉及一种用于pet系统的时间校正装置。
背景技术：
正电子发射计算机断层成像(positronemissiontomography，以下简称pet)系统是一种非侵入性的功能代谢分子影像设备，它采用的原理如下：放射性核素衰变会产生正电子，正电子运动很短的一段距离(通常为几个毫米)后就会与组织中的电子结合而发生湮灭，其质量转化成一对向相反方向射出的能量均为0.511mev的伽马光子，将标记于分子探针上的正电子核素作为示踪剂，通过探测获取正电子湮灭所产生的伽马光子的时间、能量和位置信息，运用符合技术得到湮灭事件响应线(也简称为lor，lineofresponse)，然后采用图像重建算法成像以反映各部位对示踪剂的摄取程度，进而辅助医生对疾病做出相关诊断。
符合时间分辨率是pet系统的一个重要指标，它表征了pet系统中湮灭事件发生后一对伽马光子到达的不确定度，又称为响应波动。从符合时间分辨率的定义中我们可以知道，在采用晶体阵列的探测器设计中，其响应线会相应的增多。每一条响应线对应着一个时间谱。诸多响应线对应的时间谱组合成了基本单元模块的时间谱。由于晶体阵列中各个晶体条之间存在均一性误差，即晶体阵列中各个晶体条的性能差异(包括光输出、衰减时间常数、光子渡过时间等)，以及光电转换器件区域非线性关系，使得所有时间谱的中心值都不是准确的处在理想中的零时刻，因此，基本单元模块的时间谱的展宽会被拓宽。
现有技术中大多数pet采用的时间校正方法为直接测量校正法，比较经典的直接测量校正法为旋转棒源时间校正法。旋转棒源时间校正法利用透射扫描时使用的棒源，在环绕探测器环进行衰减校正的圆周运动的同时，进行符合数据的采集以及迭代处理等操作来进行时间信息的校正(参见williamsj,luod,k.lm，etal.crystal-basedcoincidencetimingcalibrationmethod,unitedstatepatent7030382,apr.,2006)，该方法利用了经典pet中的衰减校正棒源，对于经典临床pet仪器而言不需要加入额外的装置便可以完成校正。现今市面上逐渐流行pet/ct的引进，pet/ct直接以ct成像反映的成像组织的密度信息为参考，从而生成pet衰减校正需要的衰减因子，使得大多新仪器在进行衰减校正时不再依赖旋转棒源，进一步使得该方法在pet/ct状态下的使用也面临了一定的挑战。
还有一些直接测量校正法通过在探测器环中心放置实心线源、棒源或者注射示踪剂的圆柱假体或圆环假体进行数据采集后，再使用迭代计算的方法进行时间校正(参见williansj.automatedcoincidencetimingcalibrationforapetscanner,unitedstatespatent5272344,dec.,1993)(xiaolili,burr,k.c,gin-chungwang,huinidu,etc.timingcalibrationfortime-of-flightpetusingpositron-emittingisotopesandannihilationtargets,ieeenucl.sci.symp.conf.rec.,pp.1-5,2013)。当使用线源、柱源或者注射示踪剂的直径较小的圆柱假体时，其在进行时间校正时只能涉及到较少响应线，即在圆环上的部分探测器上进行了时间信息校正；而当使用直径较大的射源圆柱假体时，在获取足够多响应线上的时间分布信息的同时，由于其体积较大，大多响应线上的时间分布会变得较为扩散，其中心偏移时刻的确认存在较大的误差。
另外，从系统的复杂程度角度而言，旋转棒源时间校正法中棒源的旋转需要使用额外的机械传动装置，比如用马达带动控制调整棒源的位置，实现精准调控，在机械部件的设计和实现上有着一定难度，增加了整个pet系统的复杂度。旋转棒源时间校正法大多使用迭代算法，基础数据量大，算法复杂度高，算法复杂度为o(n2)，o(n2)为平方阶，在此指算法中的语句执行次数，n为每一个基本单元，也可理解为晶体格的个数，该迭代算法求得最终结果的运算时间较长。旋转棒源一般采用68-ge，该放射源的半衰期是270.95天，一定时间后用户需要重新更换ge源，直接增大了pet的运营成本。
而当使用线源、柱源或注射示踪剂的直径较小的圆柱假体作为校正基础时，虽然使用的校正装置比较简单，但在校正时只能涉及到较少响应线，即圆环上部分响应线上进行了时间信息校正，并不能覆盖全视场fov(fieldofview)内的所有响应线；而当使用直径较大的射源圆柱假体时，在获取足够多响应线上的时间分布信息的同时，由于其体积较大，大多响应线上的时间分布会变得较为扩散，加大了时间分布的半高宽，不能够精确校正视场内所有响应线的时间偏移。
还有一些pet系统使用了圆环假体放在探测器环中心，但是同旋转棒源方法类似，使用了复杂的迭代算法得到最终的校正结果，采集时间虽然短，但是运算复杂度高，为o(n2)，运算时间较长。实际应用中为了提高运算速度，需要额外配备高性能计算服务器，增大了pet系统的成本。
除了上述的方法，还有一种采用时间探针的方法进行时间校正，其将射源与晶体固定在一起，置于fov正中心，以此为基准来获取其与每一个模块的时间偏差值，但是其校正单元主要为基本模块，并未精确到最小基本单元“晶体格”。而且其固定所使用的放射源为68-ge或22-na，68-ge的半衰期为270天，而22-na的2.6年，一定的时间后若想保持pet系统的性能就必须更换放射源。
综上所述可知，现有技术中对pet系统进行时间校正时存在校正效果不彻底、算法复杂化、设备成本高等问题。因此，寻找一种效果优良、运算简单快捷，算法复杂度低且低成本的时间校正方法是当今pet系统所急需的。
技术实现要素：
本发明的目的是提供一种用于pet系统的时间校正装置，从而解决现有技术中pet系统进行时间校正时校正效果不彻底、算法复杂化以及设备成本高的问题。
为了解决上述技术问题，本发明的技术方案是提供一种用于pet系统的时间校正装置，该时间校正装置包括：探测器环，所述探测器环包括若干依序排列的探测器；圆环假体，所述圆环假体位于所述探测器环内，所述圆环假体的圆心与所述探测器环的轴向及径向中心重合；探测模块，所述探测模块位于所述圆环假体内，并且所述探测模块的中心位于所述圆环假体的中心位置；数据获取模块，所述数据获取模块包括相互连接的数据采集模块和能量筛选模块，所述数据采集模块包括所述探测器和所述探测模块，所述能量筛选模块与所述数据采集模块连接并接收所述数据采集模块发送的单事件时间信息，数据符合模块，所述数据符合模块与所述能量筛选模块连接并接收所述能量筛选模块发送的单事件时间信息；时间偏移值计算模块，所述时间偏移值计算模块与所述数据符合模块连接，所述时间偏移值计算模块通过所述单事件时间信息得出所述探测器的偏移值，以及数据校正应用模块，所述数据校正应用模块用于将得到的所述偏移值应用于整个系统，从而对所述单事件时间信息进行校正。
圆环假体的轴向长度不超过所述探测器环的轴向长度，所述圆环假体的外径不超过所述探测器环的内径。
圆环假体的内径大小介于所述探测器环直径的一半和所述探测模块的外径之间。
圆环假体的厚度均匀。
圆环假体内具有射源，所述射源的活度满足条件，即放置射源后pet系统的总体计数率至少为pet系统空扫时计数率的2倍。
圆环假体内射源的活度为30～500uci。
探测模块的时间分辨率高于所述探测器环上所述探测器的时间分辨率。
探测模块为时间分辨率小于1ns的高时间性能探测模块。
探测模块包括硅酸钇镥闪烁晶体、光电转换器件以及电子学读出部分，所述硅酸钇镥闪烁晶体晶体与所述光电转换器件耦合，所述光电转换器件与所述电子学读出部分连接。
光电转换器件为光电倍增管、硅光电倍增管、多像素光子计数器或者盖革模式雪崩二极管中的一种。
数据采集模块通过所述探测器环上的探测器和所述探测模块获取列表式的单事件时间信息，所述能量筛选模块对所述数据采集模块采集到的所述单事件时间信息进行能量窗的解析筛选，滤除掉散射事件。
本发明采用圆环假体而非使用旋转棒源，在很大程度上减小了pet系统的设计，控制了pet系统设计的复杂度，使用圆环假体进行时间校正操作时还有着操作简便的优点，因此可以更好的适应当今流行的pet或者是pet/ct仪器。本发明将探测模块置于视场中心作为探测器环上探测器的时间参考，通过同一参照的方法，降低了算法的复杂度。本发明在校正时只需静态采集一次数据，算法高速简洁有效，相对于现有方法中大量使用的迭代优化校正算法，极大的降低了复杂度，减少了运算时间。
附图说明
图1是根据本发明一个优选实施例的用于pet系统的时间校正装置的布置示意图；
图2是根据本发明一个优选实施例的用于pet系统的时间校正装置的探测模块结构示意图；
图3是根据本发明一个优选实施例的用于pet系统的时间校正装置的流程示意图；
图4是根据本发明一个优选实施例的用于pet系统的时间校正装置的系统示意图；
图5是根据本发明一个优选实施例的用于pet系统的时间校正装置的时间分布示意图；
图6是根据本发明一个优选实施例的用于pet系统的时间校正装置的随机符合对比示意图，其中虚线表示未进行时间校正，实线表示应用了时间校正；
图7是根据本发明一个优选实施例的用于pet系统的时间校正装置的necr曲线对比示意图，其中虚线表示未进行时间校正，实线表示应用了时间校正。
具体实施方式
以下结合具体实施例，对本发明做进一步说明。应理解，以下实施例仅用于说明本发明而非用于限制本发明的范围。
图1为本发明提供的用于pet系统的时间校正装置的布置示意图，该时间校正装置包括探测器环10、圆环假体20和探测模块30，其中，探测器环10由若干探测器11依序排列形成圆环，圆环假体20位于探测器环10内，圆环假体20的圆心与探测器环10的轴向及径向中心重合，探测模块30位于圆环假体20内并且探测模块的中心位于探测器环10和圆环假体20的中心位置。圆环假体20的轴向长度不超过探测器环10的轴向长度，圆环假体20的外径不超过探测器环10的内径，圆环假体20的内径不小于探测模块30的外径。更优选地，圆环假体20的内径介于探测器环10直径的一半和探测模块30的直径之间。圆环假体20的厚度均匀且不易变形。为避免背景和随机的噪声干扰，圆环假体20在进行校正时需要注入射源，射源的活度优选地满足以下条件，即放置射源后pet系统的总体计数率至少为pet系统空扫时计数率的2倍左右或者更高，更优选地为10倍。圆环假体20内射源为pet系统常用的液体状正子源，例如基于18f的18f-fdg、18f-fdopa、18f-flt等，基于15o的15o-h2o、基于11c的11c-乙酸钠、11c-胆碱等放射性正子源。由于仪器灵敏度越高所需的活度越低，本发明中圆环假体20内射源的活度优选为30～500uci，更优选地为300～500uci，其成本低、灵活性好。
在图1所示的实施例中，探测器环10的形状采用的是圆环形，但本发明的探测器环并不局限于圆环形状，对于不同的环状甚至是非环状的探测器同样适用，比如，探测器环10可采用平板状、四边形或者八边形等探测器排布形式。本发明的pet系统上探测模块30采用的是高时间性能探测模块，应当理解的是，若探测模块30的时间性能优于探测器环10上的探测器11的固有时间性能，也即探测模块30的时间分辨率高于探测器环10上探测器11的时间分辨率，在本发明中均可认定为高时间性能探测模块。同时，由于临床上各种pet仪器的时间分辨率存在一定的差异，在本发明中高时间性能探测模块的时间性能小于1ns。另外，相比圆环假体20，探测模块30自身较小，可以为其他形状而不局限于圆形。
图2是根据本发明一个优选实施例的用于pet系统的时间校正装置的探测模块30的结构示意图，其中，探测模块30包括硅酸钇镥闪烁晶体31(简称lyso晶体)、光电转换器件32以及电子学读出部分33，lyso晶体31与光电转换器件32耦合，光电转换器件32连接电子学读出部分33，电子学读出部分33与数据采集部分(图中未示)的连接属于本领域的公知常识，在此不再赘述。在本发明的优选实施例中，光电转换器件32采用光电倍增管(photomultipliertube，pmt)，比如图2的实施例中采用的是滨松r9800型。
图3为根据本发明一个优选实施例的用于pet系统的时间校正装置的流程示意图，本发明的时间校正装置在使用时，包括以下步骤：
步骤s1:在探测器环10中放置并固定含有射源的圆环假体20和探测模块30；
步骤s2:初始化pet系统中的时间修正参数；
步骤s3:进行系统的符合事件获取；
步骤s4:获取各晶体格的时间相对偏移值；
步骤s5:更新并应用系统时间修正参数；其中，
步骤s3中的符合事件获取具体为获取各个晶体格相对于中心的探测模块30的符合事件，从而至少获取包含事件到达的准确时间信息和晶体格位置信息。应当理解的是，探测器中的晶体格是由探测器设计中所采用的晶体阵列规格决定的，在此不做具体限定。步骤s3与现有技术的不同之处在于，本发明的步骤s3中不仅采集了探测器环10本身的符合数据，而且采集了探测器环10叠加位于中心的探测模块30的符合数据。
步骤s4中由探测器环10中的各晶体格相对于探测模块30得到的符合事件获取符合时间分布，并从中得出时间相对偏移值tn，其中n为晶体格编号，n为正整数。传统技术中使用迭代算法进行校正参数的方法相比，本发明的步骤s4使用的是直接测量计算的方法获得每一个晶体格上的时间偏移值，进而得到校正参数，算法高速简洁有效，极大的降低了复杂度，减少了运算时间。
步骤s5中将得到的各晶体格时间相对偏移值tn应用于系统中，具体方法为在单事件获取时对原始数据中的时间信息进行修正，即在原始的时间信息中扣除或添加上得到的时间相对偏移值tn，而具体采取扣除或者添加取决于时间相对偏移值的计算方式。时间相对偏移值的计算方式有两种，第一种是t中心-t圆环，此时校正计算为添加时间相对偏移值tn；第二种是t圆环-t中心，此时校正计算为扣除时间相对偏移值tn。
图4为根据本发明一个优选实施例的用于pet系统的时间校正装置的系统示意图，由图4可知，本发明提供的时间校正装置的系统包括数据获取模块100、数据符合模块200、时间偏移值计算模块300以及数据校正应用模块400，其中，数据获取模块100包括相互连接的数据采集模块110和能量筛选模块120，能量筛选模块120与数据符合模块200连接并向数据符合模块200传递采集到的单事件时间信息。
数据获取模块100用于获取探测模块和探测器环上各探测器的单事件时间信息，并对该单事件时间信息应用能量窗进行部分预处理操作，其中，数据采集模块110用于获取探测模块和探测器环上各探测器的列表式的单事件时间信息，传输并存储至上位机服务器留待处理。能量筛选模块120用于对数据采集模块110采集到的原始的单事件时间信息进行一定能量窗的解析筛选，滤除掉散射事件。例如，在图1的实施例中，数据采集模块110为探测器环10上的各探测器和探测模块30的总和，数据采集模块110获取探测器环10上的探测器11和高性能探测模块30生成的列表式的单事件时间信息，数据采集模块110通过与能量筛选模块120连接，将上述单事件时间信息传输并存储至上位机服务器留待处理，上位机服务器中的能量筛选模块120对数据采集模块110采集到的原始的单事件时间信息进行一定能量窗的解析筛选，滤除掉散射事件。应当理解的是，对于能量窗的选择，本领域技术人员可以根据具体信息采集的需要进行确定，在此不再赘述。
需要注意的是，数据获取模块100中的数据采集模块110和能量筛选模块120的前后顺序在本发明的时间校正装置中可以对调，即可以将能量筛选模块120集成于探测器环10上的各探测器11和探测模块30中，从而可以先通过能量筛选模块120将散射事件滤除，然后再通过数据采集模块110上传到上位机服务器。当射源活动比较高时，该模式下数据采集效率更高。
数据符合模块200用于将数据获取模块100获取到的单事件时间信息根据设定时间窗及其他条件进行数据符合，得到每一个探测器上每一个晶体格相对于探测模块的符合数据。比如，能量筛选模块120设定几个不同的筛选范围，例如350-650kev，250-750kev等，数据符合模块200据此进行数据符合，得到探测器环10上每一个探测器的每一个晶体格对应于高时间性能探测模块30的符合数据。
时间偏移值计算模块300用于对每一个晶体格相对于探测模块的符合数据进行时间分布获取，并得出其偏移值tn，其中n为晶体格编号，n为正整数。
数据校正应用模块400用于将得到的各个晶体格的偏移值tn应用于整个系统，对获取到的每一个单事件中的时间信息进行时间校正操作。具体为对数据获取模块100获取的列表式的单事件时间信息中的时间信息去除所获得的偏移值tn进行时间信息修正。偏移值tn的计算方式以及校正方式参见上文步骤s5，在此不再赘述。
图5是根据本发明一个优选实施例的用于pet系统的时间校正装置的时间分布示意图，其中，虚线e表示时间校正前探测器环10上所有探测器组成的系统的时间分布，其时间分辨率为1.7ns，中心值为121.8ps；实线f表示时间校正后所有探测器组成系统的时间分布，其时间分辨率为1.26ns，中心值为-3.6ps。相对于现有技术中的400-500ps的符合时间分辨率，本发明的时间分辨率更有助于tof-pet技术(time-of-flightpet，也称为飞行时间pet技术)的实现。
实现系统时间校正前后的随机符合以及necr曲线分布分别如图6，7所示，其中，图6是根据本发明一个优选实施例的用于pet系统的时间校正装置的随机符合对比示意图，其中虚线表示未进行时间校正，实线表示应用了时间校正；图7是根据本发明一个优选实施例的用于pet系统的时间校正装置的necr曲线对比示意图，其中虚线表示未进行时间校正，实线表示应用了时间校正。随机符合事件会随着时间窗口的减小而减小，时间窗口的大小又是根据时间分辨率来进行设定的，一般为时间分辨率的两倍。随机符合事件少反映在最终成像上就是图像质量(信噪比)有提升。necr曲线通过necr假体按照nema-nu4标准中的方法测试得出，从图6中可以看出时间校正应用后随机符合分布较校正前有降低，且在图7中可以看出校正前后在真事件没怎么变化的情况下necr有了一个较大的提升，提升幅度可从图6中看出，在整个活度范围下计数率减小了近一倍。
本发明采用圆环假体而非使用旋转棒源在很大程度上减小了pet系统的设计，控制了pet系统设计的复杂度，使用圆环假体进行时间校正操作时还有着操作简便的优点，因此可以更好的适应当今流行的pet或者是pet/ct仪器。
本发明引入了一个探测模块，尤其是高时间性能探测模块置于视场中心作为探测器环上探测器的时间参考，通过同一参照的方法，降低了算法的复杂度。本发明在校正时只需静态采集一次数据，并且可以得知此时的时间复杂度为o(n)，n为探测器环上总的晶体格的数目。算法高速简洁有效，相对于现有方法中大量使用的迭代优化校正算法，极大的降低了复杂度，减小了运算时间。
本发明还可以针对探测器环上的最小探测单元，即小晶体格做晶体时间校正。综合其计算方法，可以提高运算校正效率。
以上所述的，仅为本发明的较佳实施例，并非用以限定本发明的范围，本发明的上述实施例还可以做出各种变化，例如，本发明中的探测模块30的结构可以为晶体/光电转换器件/电子学读出，也可以为晶体/光导/光电和转换器件/电子学读出，或者是光电转换器件/电子学读出czt探测器；再再比如，本发明的高时间性能探测模块所使用的晶体可以是labr3、lso、lyso、luyap、baf2、gso、lfs或者lui3等多种种类，其形状可以为圆柱形，长条形，锥形方形等多种形状；再比如，本发明中高时间性能探测模块所使用的光电转换器件可以为光电倍增管、硅光电倍增管(siliconphotomultiplier,sipm)、多像素光子计数器(multi-pixelphotoncounter,mppc)、盖革模式雪崩二极管(geiger-modeavalanchephotodiode，g-apd)等，只要渡越时间分布tts小于500ps都可以作为本发明中的光电转换器件。总之，凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰，皆落入本发明专利的权利要求保护范围。