中国医学物理学杂志
Chinese Journal of Medical Physics 중국의학물리학잡지
- 主管单位： 南方医科大学
- 主办单位： 南方医科大学;中国医学物理学会
- 影响因子： 0.48
- 审稿时间： 3-6个月
- 国际刊号： 1005-202X
- 国内刊号： 44-1351/R
- 论文标题 期刊级别 审稿状态
-
基于宽带检测放疗X-光光声效应仿体实验
目的:搭建一个宽带检测的单振元光声信号采集系统采集仿体的光声信号,探讨不同介质X-光光声信号差异和人体切片体模光声信号的去噪方法.方法:带宽为0.10～1.00 MHz的单振元宽带换能器垂直于X-光照射方向在样本一侧采集超声信号,超声信号通过34 dB带宽为0.05～5.00 MHz的前置放大器放大后由数字采集卡(采样率=40.00 MHz)采集并存储到计算机中.首先通过光声信号在整体包络、信号强度、振荡趋势、上下包络所围成的面积等方面差异对4种仿体进行分析.另外对人体切片模体的X-光声信号在有限长单位冲激响应低通滤波的基础上进行小波平移不变量去噪与只对模体光声信号进行小波平移不变量去噪,比较这两种去噪方式的效果.结果:不同仿体光声信号有差异,这种差异与物质的成分和比例相关.对人体切片体模光声信号通过Db3小波基进行5层分解的小波平移不变量去噪,比有限长单位冲激响应低通滤波能更好地提高信号信噪比.结论:利用6 MV医用直线加速器的短脉冲X-光光声信号对介质的差异性及小波去噪有助于提高其对放疗的研究与应用价值.
-
宫颈癌术后5野调强放疗4个变量组合的佳治疗模式的剂量学
目的:探讨宫颈癌术后5野固定野调强方式(5F-IMRT)应采用的佳治疗模式.方法:首先用改变单一条件的方法,分别比较治疗体位、剂量计算算法、加速器的多叶光栅(MLC)调强方式以及高能X射线能量这4类条件,共计8个变量,包括仰卧位与俯卧位、各向异性算法(AAA)与笔形束卷积算法(PBC)、动态调强即滑窗(SW)与静态调强即静态分段(MSS)、6MVX线与15 MVX线,对靶区、危及器官和正常组织的剂量学差异影响.其次,基于上述比较结果,再比较这4类条件中相对更优的4个变量串联而成的B组和与相对较差的4个变量串联而成的A组对靶区、危及器官和正常组织的剂量学影响.其中剂量学评价指标包括靶区、危及器官和正常组织的剂量受量、靶区剂量适形度指数(CI)、靶区均匀性指数(HI)、机器跳数(MU)、治疗时间(T)等.数据统计采用SPSS 19.0软件,并根据不同的统计对象,选择合适的统计方法(独立样本t检验与配对t检验)进行统计分析,P≤0.05为差异具有统计学意义.结果:在单一条件统计结果中,俯卧位较仰卧位而言,靶区平均剂量(Dmean)、小肠大剂量(Dmax)及膀胱V50(膀胱受到50 Gy剂量的体积占膀胱总体积的百分比)均较低:AAA较PBC而言,计划靶区的CI偏大而HI偏小,膀胱V50、直肠Dmax均偏小;SW较MSS而言,靶区的CI偏大而HI偏小,膀胱V50、直肠Dmax均偏小,但MU和T偏大;高能X射线15MV较6MV而言,靶区Dmean、CI偏大,而HI偏小,危及器官中的脊髓Dmax(考虑摆位误差,实际脊髓外放3mm形成外放脊髓)、膀胱V50、股骨头Dmax均偏小.在多条件变量串联统计结果中,采用俯卧位、AAA、SW以及15 MVX线这4个变量串联的B组与与之对应的仰卧位、PBC、MSS、6 MVX线这4个变量串联的A组相比,B组中靶区Dmean平均降低了1.4％,CI平均增加了5.8％,HI平均降低了53.4％,危及器官中小肠Dmax平均降低了4.4％,直肠Dmax平均降低了2.2％、膀胱V50平均降低了109.1％、股骨头Dmax平均降低了1.5％.结论:宫颈癌术后5F-IMRT的放疗体位应采用俯卧治疗体位、剂量算法应采用AAA、MLC运动方式应采用动态SW、X射线应采用15MV高能X射线这4个变量组成的治疗模式,临床剂量学研究表明该治疗模式为本实验对宫颈癌术后5F-IMRT的佳治疗模式.
-
基于四维CT和锥形束CT确定非小细胞肺癌放疗靶区外放边界
目的:探讨非小细胞肺癌(NSCLC)放疗靶区外放边界的确定.方法:首先随机抽取14例做体部立体定向放疗(SBRT)的NSCLC患者.所有患者序贯行平扫3DCT和四维CT(4DCT)扫描,采用4DCT定位和锥形束CT(CBCT)治疗前扫描评估放疗靶区随呼吸运动的位移和摆位误差,然后结合呼吸运动位移和摆位误差数据计算出NSCLC放疗靶区外放边界.结果:肿瘤中心点左右、前后、上下方向位移分别为(0.18±0.11)、(0.21±0.14)、(0.86±0.85)cm.肿瘤边界3个方向位移分另为(0.48±0.30)、(0.59±0.44)、(0.50±0.33)cm.摆位系统误差3个方向分别为(0.03±0.24)、(0.02±0.43)、(0.02±0.26) cm.利用肿瘤中心点、边界和摆位位移得出NSCLC肿瘤外放边界3个方向大为0.72、0.74、1.30 cm.结论:放射治疗肿瘤组织(RTOG)推荐肺部肿瘤外放边界为左右、前后0.50 cm,上下1.00 cm,与之相比,利用本研究方法计算各个方向的外放边界比原来RTOG提出的外放标准更加精确,给个体化放疗外放边界的选取提供了很好的依据.
-
准直器角度对宫颈癌术后VMAT计划的影响
目的:探讨准直器角度对宫颈癌术后旋转调强(VMAT)计划的影响.方法:选取一名宫颈癌术后病例进行VMAT计划设计(CCW179°～81°、CW181°～79°),采用不同的准直器角度:(0°,0°)、(0°,10°)、(0°,20°)、(0°,30°)、(0°,45°)、(-10°,10°)、(-20°,20°)、(-30°,30°)、(-45°,45°).在Varian Eclipse系统上设计9个计划,从中筛选出剂量学优和差的两个角度,将此两个角度应用于20例宫颈癌术后患者,比较两组计划的靶区适形度(CI),剂量均匀性指数(HI),危及器官的体积剂量V10、V20、V30、V40,平均剂量(Dmean),以及跳数(MU).用SPSS软件的配对t检验方法统计靶区和危及器官的剂量学参数.结果:以折线图表示9个计划的CI和HI,以及膀胱、直肠等其他危及器官的体积剂量参数,剂量学参数优和差计划的准直器角度分别是(0°,30°)、(0°,0°).20例宫颈癌术后的(0°,30°)和(0°,0°)两组计划当中,(0°,30°)计划靶区的HI和CI均优于(0°,0°)计划(HI:t=6.666、6.231、15.464;CI:t=6.650、6.35、14.981,P＜0.05),膀胱的V30、v40、Dmean分别降低了13.5％、3.4％、3.7％(t=4.649、4.511、3.833,P＜0.05);直肠的V30、V40、Dmean分别降低了13.4％、30.3％、6.6％(t=4.945、4.455、8.864,P＜0.05);小肠的V30 V4、Dmean分别降低了30.2％、48.4％、2.7％(t=3.972、4.077、9.604,P＜0.05),加速器执行一次计划的平均MU增加了12.7％(t=7.028,P＜0.05).左右股骨头以及骶骨未见明显规律.结论:准直器角度对宫颈癌术后VMAT计划有明显的影响,在宫颈癌术后VMAT计划设计过程中,应充分考虑准直器角度的因素,以寻求更优质更有执行效率的治疗计划.
-
床面倾斜对剂量验证通过率的影响
目的:基于二维二极管阵列探测器模拟床面倾斜对放疗剂量验证通过率的影响.方法:选择调强放射治疗(IMRT)和三维适形放疗(3DCRT)计划患者各10例的剂量验证计划,在加速器治疗床用边缘充气袋以左右激光灯示值差模拟床面倾斜,倾斜幅度选择1、2、3、5、7、10 mm.选择Gamma(r)分析法,(3％,3mm)绝对剂量分析标准,记录计划和模拟实验结果,用配对t检验分析结果.结果:在无床面倾斜条件下,剂量验证通过率:IMRT为(96.10±1.65)％,3DCRT为(97.50±1.98)％.IMRT患者在倾斜7、l0 mm的通过率分别为(86.50±5.95)％、(85.40±8.53)％,与无倾斜通过率比值分别为0.90±0.05和0.85±0.08;3DCRT患者在倾斜大(10 mm)的通过率为(91.90±4.90)％,与无倾斜比值为0.94±0.05.与无倾斜对比,除了3DCRT在倾斜2 mm没有统计学意义(P=0.366),其余倾斜幅度差异均有统计学意义(P＜0.05).结论:床面倾斜对3DCRT剂量验证通过率影响较小;而IMRT通过率随床面倾斜幅度增大而减小,幅度较大,但变化幅度存在个体性差异,在倾斜大于5 mm时可能对剂量验证通过率产生较大影响.保持床面无倾斜是日常临床工作中一项简单但重要的质量保证工作.
-
容积旋转调强技术用于直肠癌时佳机架角度设置研究
目的:探讨对直肠癌患者应用容积旋转调强(VMAT)技术时是否需要避开膀胱的问题.方法:随机选取12例直肠癌放疗患者的CT图像进行靶区及危及器官(OAR)勾画,并对同一CT图像分别设计“不避,避开30°～330°,40°～320°,50°～310°,60°～300°,70°～290°,80°～280°”7个计划,每个计划采用两个旋转弧进行照射,评估靶区及OAR的剂量分布.利用方差分析不同避开方法的结果是否存在显著性差异,若有就找出好的避开方法,并分析其原因.结果:7个计划的靶区覆盖度均满足临床处方剂量要求,不避计划的靶区适形度和靶区均匀性均较其他设计方案优(P=0.000,P=0.000),而且不避计划处方剂量的95％所包绕的体积V95多(P=0.000),处方剂量的105％包绕的体积V105、Dmean 、Dmax、D5低(P=0.000、P=0.000、P=0.000、P=0.000),膀胱V40、V50、Dmax低(P=0.044、P=0.043、P=0.000),小肠、骨髓、Body的D.低(P=0.039、P=0.000、P=0.000).避开60°～300°的计划的膀胱V20低(P=0.005).角度变化对其他量的影响不显著.结论:在对直肠癌患者应用VMAT技术时,可以不用避开膀胱,用两个全弧照射.
-
基于新的变步长小均方算法在运动环境中实时提取心率的方法
针对光电容积脉搏波容易受运动伪影干扰的问题,为不同运动状态下可穿戴设备的心率提取提供一种新的变步长小均方(HSFLMS)算法.根据不同的运动状态分别进行去噪处理:静止状态采用波峰波谷计数方法,运动状态则采用HSFLMS算法.HSFLMS算法根据双曲正弦函数模型,构造了一种新的关于步长和误差的非线性关系函数,并分析了参数取值对算法性能的影响.根据实验结果可知,相较于其他算法,HSFLMS算法的结构框架更简洁、更高效,能有效提取不同运动状态下的实时主体心率.
-
一种基于虚拟仪器的交互式心电信号源系统
目的:为满足对心电类仪器的研发检测以及对其后续使用进行质量监控、心电检测算法的验证、心电自动诊断研究领域对特殊心电信号的需求,研究设计一种基于虚拟仪器的交互式心电信号源系统.方法:该系统运用软件技术通过绘图和插件组合编辑方式获得所需心电信号,通过DAQ硬件卡和归一化双端输出电路实现以相当于人体真实水平输出.结果:该系统实现了可人为地构建所需心电波形,调用、修改、显示已有心电信号,并可将其以相当于人体真实水平进行输出.结论:该心电信号源系统具有强大的交互式编辑扩展功能,理论上可产生任意所需心电信号,在心电类仪器的研发、质量检测、计量检定校准、心电检测算法的验证、心电自动诊断研究等领域有非常广阔应用前景.
-
陡脉冲电场选择性灭活血液中病原体效果观察
目的:观察纳秒级陡脉冲在10 ns脉宽参数下对病原体的杀伤效应,以及对红细胞的损伤作用,研究高压脉冲电场应用于血液病原体灭活的可能性.方法:以大肠埃细菌为研究对象,采用脉冲参数40～80 kV/cm、脉宽10 ns、重复频率1 Hz的陡脉冲进行实验处理,实验菌株活力以MTT法进行检测,三波长法检测处理前后的血液样品的游离血红蛋白(FHb)浓度.结果:在电场强度40 kV/cm、10个脉冲作用下活性抑制率达到54.5％,不同介质中脉冲数对细菌活性均有显著影响(P＜0.05);不同脉冲电场强度越大其抑制率有所升高,但均无统计学意义(P＞0.05);无菌条件下,脉冲数对红细胞的影响显著(P＜0.05);有菌条件下,脉冲数及电场强度对溶血仍有影响,但统计结果不显著(电场强度P=0.159,脉冲数p=0.278);FHb在10个脉冲数的作用下,没有显著增加,甚至比对照组更低.结论:纳秒级陡脉冲电场具备较强的杀菌能力,且在一定参数条件下不对血细胞产生明显的损伤作用,具备选择性杀伤的潜在作用.
-
3D打印技术在口腔种植领域的应用进展
本文对3D打印技术在口腔种植导板和个性化种植体两个方向上的应用进行阐述,展示这项技术使种植技术发生的转变,并分析了这项技术在应用中存在的缺点和有待解决的问题.3D打印技术的应用将为口腔种植领域的发展带来深远和重大的影响.
-
中药类光敏剂金丝桃素的光谱学特性
本研究系统地研究了金丝桃素(HY)的紫外-可见吸收光谱、光漂白和荧光光谱特性,为HY的临床应用提供一定的参考.实验测量了不同浓度HY-二甲基亚砜溶液的紫外-可见吸收光谱,结果显示HY在可见光范围内存在多个吸收峰,其中,594nm处的吸收强,这表明594nm为HY光动力治疗的首选光源波长;HY的光漂白特性研究表明,HY的初始浓度越高,光漂白所需的时间越长;HY的荧光光谱研究表明,荧光光谱在可见光范围内也存在多峰结构,其中,643 nm处的荧光峰强且半峰宽窄,因此,在光动力诊断中,采用643 nm荧光发射能大大提高检测的灵敏度.
-
超声造影剂微泡次谐波辅助压力估测技术研究进展
超声造影剂微泡在超声激励下可以激发出次谐波,幅度随着环境压力的增加而线性降低,利用这一关系可以非侵入性地测量体内组织的压力水平.目前国内外的研究包括数值计算、体外模拟和体内测量等几个方面,均得到了比较好的结果,表明次谐波辅助压力估测技术是一种比较好的非侵入性测量体内局部压力的方式.通过参数优化、微泡选择、扫描方式改进、数据提取和处理等方面的研究为其临床实际应用奠定了基础.
-
基于组织等效3D打印的多模椎骨体模设计
目的:设计一种基于影像建模、组织等效和3D打印技术的多模态椎骨骨折影像学仿真体模.方法:采集l例健康人腰椎影像样本的腰椎段CT影像序列,将原始DICOM图像导入Mimics交互式的医学影像控制系统软件进行椎骨、肾脏和软组织的三维重建,在重建模型中构建“一”字缝、“Y”字缝的椎骨峡部骨折模型.使用3D打印和组织等效技术打印等X线吸收密度的椎骨骨折模型,经CT、数字X线摄影(DR)多模成像,对模型CT值进行评估,验证体模仿真程度.结果:所获得的椎骨骨折体模在CT、DR下成像清晰,组织X线吸收系数等效于真实人体组织.结论:基于3D打印技术的多模态椎骨骨折体模可在CT和DR影像下获得,并可应用于有限元生物力学分析.模体制做成本低廉且快速,组织等效参数符合正常人体参数范围,能够有效仿真3种不同类型的椎骨峡部裂骨折等人体腹部病变.
-
基于移动设备的三维图像重建架构及其应用
随着移动互联网的兴起,各种移动设备诸如平板电脑和智能手持设备乃至智能手机已渐渐在医疗行业中崭露头角.目前随着移动设备的快速发展,其硬件性能和软件平台越来越强大,在放射科以及其它临床科室中移动设备也逐渐得到了应用,使真正意义上的移动阅片、移动报告成为可能.三维图像重建是影像诊断和临床手术规划中的重要高级功能,以往需要借助专业工作站,目前在移动设备上也已经可以实现.本文主要综述了在移动设备上实现三维图像重建的几种架构方法,包括基于移动设备硬件进行重建、基于Web技术的服务器的远程重建以及基于移动设备的渐进式重建和体绘制.对上述几种方案的关键技术进行了具体描述,并针对几种架构具体分析了其应用场景,为基于移动设备的开发和应用提供借鉴.
-
脑部弥散张量成像数据处理技术在疾病诊断研究中的应用进展
主要介绍弥散张量成像(Diffusion Tensor Image,DTI)技术的发展,介绍脑部DTI数据处理常用的方法应用于疾病诊断的发展概况.本文总结DTI数据处理常用的3种方法,包括基于体素的分析、基于纤维束的空间统计分析方法和白质纤维追踪方法,其中列出了使用各种方法应用于疾病研究所取得的成果,对研究结果进行了分析并对各种分析方法的优缺点进行总结.后对DTI的发展方向提出展望,DTI在科研和临床应用方面都发挥着越来越重要的作用,尤其是科研方面已取得很多重要成果.在临床方面,DTI可用于多种疾病的诊断和治疗,且在某些疾病的诊断中发挥着越来越重要的作用.
-
PET/CT增强在鼻咽癌头颈部淋巴结转移的应用价值
目的:研究PET/CT同机CT增强一站式扫描对鼻咽癌患者头颈部淋巴结转移的检测价值.方法:对61例确诊鼻咽癌患者行PET/CT同机CT对比剂增强一站式扫描,分析PET图像、CT图像和PET/CT增强图像对淋巴结检测的特点,分成3组数据,对比各组对颈部淋巴结转移诊断的差异.结果:61例患者共发现213个病灶,其中CT组发现阳性病灶137个,阴性病灶76个,灵敏度83.8％,特异性86.4％;PET组发现阳性病灶156个,阴性病灶57个,灵敏度95.5％,特异性86.2％;PET/CT增强组发现阳性病灶153个,阴性病灶60个,灵敏度99.3％,特异性98.3％.结论:PET/CT增强扫描在鼻咽癌头颈部淋巴结转移的诊断中,无论是对淋巴结的检出,还是在良恶性判断方面都有重要的应用价值.
-
四维CT扫描参数和扫描模式对CT值的影响
目的:研究四维CT的不同扫描参数和模式对图像CT值的影响.方法:在扫描条件一致的情况下,测量电子密度模体中相应的组织替代材料的CT值均值及标准差,分析曝光时间和旋转时间参数对电影扫描和轴向扫描模式的影响,并比较轴向扫描、螺旋扫描和电影扫描3种模式之间的差异.结果:电影扫描模式下,CT值相差较大的有CB2 50％(含有50％碳酸钙的骨)和致密骨两处:旋转时间分别为0.5 s/rot和1.0 s/rot时,CT值差值分别为14.7 HU和23.1 HU,均小于3.0％;轴向扫描模式下,致密骨处CT值相差较大:旋转时间分别为0.5 s/rot和4.0 s/rot时,CT值差值为27.1 HU(2.4％).相同扫描条件下对同一组织扫描,电影扫描与轴向扫描模式间差别可忽略不计;螺旋扫描与轴向扫描模式,对低密度组织成像CT值差别较大,如呼气态肺和吸气态肺相对密度误差分别为4.3％和2.3％,致密骨的CT值误差较大为46.0 HU,换算成相对电子密度误差为1.6％.结论:旋转时间与曝光时间参数对电影扫描及轴向扫描的影响基本可忽略不计.相同扫描条件下对同一组织扫描,3种扫描模式间差别可忽略不计.
-
X线隧道位测量膝关节骨性关节炎患者股骨髁间窝宽度和高度的临床意义
目的:膝关节X线隧道位测量膝关节骨性关节炎患者与健康成人髁间窝宽度与高度,分析髁间窝高度、宽度与骨性关节炎的关系.方法:选取132例健康成人与97例膝关节骨性关节炎患者,进行膝关节X线隧道位摄片,并将所得结果按性别进行分组、比较.结果:健康成人组髁间窝宽度为(19.45±2.76) mm,高度为(27.93±2.75)mm;骨性关节炎组宽度为(17.74±2.14)mm,高度为(22.04±2.43) mm,二者均有明显统计学差异(P＜0.05).结论:骨性关节炎患者的髁间窝宽度与高度明显小于健康成人,膝关节骨性关节炎同髁间窝的狭窄有关;X线隧道位是显示股骨髁间窝一种有临床价值的影像学方法.
-
彩色多普勒超声联合磁共振胆胰管成像诊断老年胆总管结石270例临床分析
目的:分析彩色多普勒超声(CDU)联合磁共振胆胰管成像(MRCP)诊断老年胆总管结石的价值并讨论其意义.方法:选择2014年1月～2015年12月拟诊断胆总管结石患者270例,男性116例,女性154例,年龄60～85岁.患者均行CDU和MRCP检查,同时以手术或内镜下逆行胰胆管造影术作为金标准,用统计学方法比较单独用CDU或MRCP,以及CDU联合MRCP对诊断老年胆总管结石的价值,同时计算敏感性、特异性和准确性.结果:老年患者270例经手术或ERCP确诊胆总管结石218例,CDU联合MRCP诊断老年胆总管结石的敏感性、特异性和准确性均较单纯CDU检测或MRCP检测高,差异具有统计学意义(P＜0.05).结论:CDU联合MRCP检查可提高老年胆总管结石诊断准确率,且简便安全,可在临床推广应用.
-
基于医学数字X线影像频谱特征的图像检索
目的:针对数字X线影像的频率特征,提出基于频谱相似性比较的图像检索新方法.方法:通过数字X线影像分析,获得其频率分布和幅值特性,根据不同部位(摄影体位)影像频率特性差异明显的特点,对影像进行特征频率曲线的提取、分类和库构建.对任意输入的图像,进行基于库曲线的相似性比较,从而判断图像的摄影部位(摄影体位).结果:采集影像科临床X线影像共计200幅,进行基于频率特征的曲线相似性判断检索,判定影像所属摄影部位.检索结果经临床医生进行主观判断,平均准确率93.78％,平均判断时间0.290 3 s.结论:该方法基于数字X线影像的频谱特征进行图像检索,首次分析不同部位在频率曲线上的特征性,将二维图像检索简化为一维曲线检索,大大缩短了检索时间,具有较高的准确率和鲁棒性.
-
基于四维CT和TumorLoc软件测定肺下叶呼吸动度
目的:利用四维CT(4D-CT)和TumorLoc软件,研究肺下叶(右膈肌层面)距离脊柱不同位置处的呼吸动度.方法:采用放疗专用Philips Brilliance 24排大孔径CT定位机对10例行真空垫固定的患者进行4D-CT模拟定位扫描,将每个呼吸周期的CT图像平均分为10个呼吸时相.通过TumorLoc软件打开每例患者的10个呼吸时相图像,获得肺下叶内(右膈肌层面)距离脊柱40、50、60、70、80、90 mm处血管中心点在三维方向的位移,分析位移变化及左右距离脊柱相同距离位置处三维方向的相关性.结果:左肺下叶(右膈肌层面),距离脊柱40、50、60、70、80、90mm位置处,呼吸动度在Z方向(头脚)分别为(9.5±2.5)mm、(9.7±2.6)mm、(9.5±2.5)mm、(9.3±2.3)mm、(9.7±2.5)mm、(9.5±2.6)mm;右肺下叶(右膈肌层面),距离脊柱40、50、60、70、80、90 mm位置处,呼吸动度在Z方向(头脚)分别为(10.5±2.7) mm、(11.4±3.1)mm、(11.3±3.2) mm、(11.5±3.0) mm、(11.6±4.0) mm、(11.7±4.3) mm;左右相同距离位置处,X方向(左右)、Y方向(前后)差异无统计学意义(P＞0.05).左右相同距离位置处,Z方向(头脚)在40、50、60 mm处差异有统计学意义(P分别为0.005、0.007、0.005);Z方向在70、80、90 mm处差异无统计学意义(P＞0.05).结论:应用4D-CT通过TumorLoc软件可精确测量肺下叶(右膈肌层面)不同位置处在三维方向的呼吸运动度.
| 年 | 期数 |
| 2019 | 01 02 03 |
| 2018 | 01 02 03 04 05 06 07 08 09 10 11 12 |
| 2017 | 01 02 03 04 05 06 07 08 09 10 11 12 |
| 2016 | 01 02 03 04 05 06 07 08 09 10 11 12 |
| 2015 | 01 02 03 04 05 06 |
| 2014 | 01 02 03 04 05 06 |
| 2013 | 01 02 03 04 05 06 |
| 2012 | 01 02 03 04 05 06 |
| 2011 | 01 02 03 04 05 06 |
| 2010 | 01 02 03 04 05 06 |
| 2009 | 01 02 03 04 05 06 |
| 2008 | 01 02 03 04 05 06 |
| 2007 | 01 02 03 04 05 06 |
| 2006 | 01 02 03 04 05 06 |
| 2005 | 01 02 03 04 05 06 |
| 2004 | 01 02 03 04 05 06 |
| 2003 | 01 02 03 04 |
| 2002 | 01 02 03 04 |
| 2001 | 01 02 03 04 |
| 2000 | 01 02 03 04 |
| 1999 | 02 03 04 |
