中国医学物理学杂志
Chinese Journal of Medical Physics 중국의학물리학잡지
- 主管单位： 南方医科大学
- 主办单位： 南方医科大学;中国医学物理学会
- 影响因子： 0.48
- 审稿时间： 3-6个月
- 国际刊号： 1005-202X
- 国内刊号： 44-1351/R
- 论文标题 期刊级别 审稿状态
-
剂量网格分辨率大小对非小细胞肺癌立体定向放射治疗的剂量学影响
目的:定量分析剂量网格分辨率的大小对非小细胞肺癌(NSCLC)立体定向放射治疗(SBRT)计划剂量分布的影响,指导临床选用合适的剂量网格分辨率用于肺癌SBRT计划设计.方法:选取10例NSCLC患者,采用容积旋转调强技术,使用0.20 cm的剂量网格分辨率设计SBRT计划,将治疗计划结果再分别用0.40、0.30、0.25、0.15、0.10 cm的剂量网格分辨率计算终剂量.比较6种不同剂量网格分辨率下的计划靶区(PTV):D2％、Dmean、D98％、均匀性指数(HI)、适形度指数(CI)和危及器官:全肺、胸壁、食管、心脏、脊髓、主动脉、气管树等的相关剂量学参数的差异.结果:与0.20 cm剂量网格分辨率组计划相比较,0.40、0.30、0.25、0.15、0.10 cm剂量网格分辨率计算得到的计划的靶区D2％、Dmean、CI均具有统计学意义(P<0.05);HI除0.15 cm剂量网格分辨率组以外,均具有统计学意义(P<0.05).在危及器官受照剂量方面,与0.20 cm剂量网格分辨率组计划相比较,大于0.20 cm(0.40、0.30、0.25 cm)组计划全肺、胸壁、食管、心脏、脊髓、主动脉、气管树等的相关剂量学参数差异均有统计学意义(P<0.05);小于0.20 cm(0.15、0.10 cm)组计划除全肺V10、V20、V12.5和V13.5有影响以外(P<0.05),胸壁、食管、心脏、脊髓、主动脉、气管树等的相关剂量学参数差异较小(P>0.05).结论:剂量网格分辨率的大小会影响剂量计算的准确性,在NSCLC患者SBRT计划设计时,建议使用0.20 cm或更小的剂量网格分辨率.
-
呼吸运动对Sliding Window调强射野剂量学的影响
目的:研究呼吸运动对Sliding Window(SW)调强射野剂量学的影响.方法:使用胶片分别测量多叶光栅在不同的叶片移动速度(叶片速度)和不同射野宽度条件下,窄束野在呼吸运动平台驱动下胶片上的剂量分布.对运动模式、叶片速度、射野宽度等多个因素分组进行比对分析.结果:不同的运动模式对剂量分布影响差异显著不同;随着叶片速度的增大,射野内剂量波动幅度随之显著增大;随着射野宽度的增大,射野内剂量波动幅度随之减小;射野宽度小叶片速度高,对剂量分布的影响为显著.结论:有呼吸运动影响时,应通过提高射野宽度和(或)减小叶片速度等方法来减小呼吸运动对SW调强射野剂量分布的影响.
-
权重适配的布拉格峰展宽方法
目的:提供一种权重适配的布拉格峰展宽(SOBP)方法,得到平滑的展宽布拉格峰.方法:通过重新拟合质子能量-射程的关系(盖格法则),找出适配函数的函数形式,并对权重进行重新适配,通过求敏感参数k,得到平滑的SOBP,后用蒙特卡洛程序FLUKA进行验证.结果:SOBP的形状对参数k比参数P更加敏感,拟合得到4～32 cm的SOBP,中间平坦区偏差不超过±2％,并解决中间区坍塌的问题.结论:蒙特卡洛模拟检验了权重适配的SOBP方法的有效性.
-
宫颈癌放疗中患者体型改变的特征量与摆位误差的关联
目的:评估宫颈癌患者在螺旋断层放疗(TOMO)中耻骨联合上缘层面的横纵轴宽度比值(简称横纵比)与摆位误差的关联.方法:选取行TOMO的宫颈癌患者30例,采用热塑体膜固定,利用图像引导系统获取分次摆位MVCT图像,并与定位kVCT图像进行配准,记录左右、腹背、头脚方向上的摆位误差.记录每分次患者的横纵比,并将其与摆位误差进行Pearson相关性分析.结果:对于横纵比稳定性>0.8的18例患者,其横纵比与头脚和腹背方向上的随机误差、年龄、体质量指数都呈显著负相关,关联系数R分别为-0.621、-0.488、-0.560、-0.552(P=0.006、0.040、0.016、0.017).横纵比稳定性<0.8的12例患者中,其年龄、体质量指数与δ横纵比的相关性系数R分别为-0.796、-0.566(P=0.002、-0.055).结论:在宫颈癌患者TOMO中,横纵比可以作为一个临床观测指标,评估摆位误差的大小.
-
宫颈癌螺旋断层放疗与静态调强技术的剂量学比较
目的:比较宫颈癌盆腔照射螺旋断层放射治疗(Helical Tomotherapy,HT)计划与七野-静态调强放疗(7F-IMRT)计划的靶区和危及器官的剂量学差异.方法:选取2015年3月～2017年10月行HT治疗的25例宫颈癌患者,均采用俯卧位,憋尿状态.分别进行HT和7F-IMRT两种计划设计,比较两种计划模式之间靶区和危及器官之间的剂量学参数及机器输出跳数.结果:(1)计划靶区(PTV):HT组CI和HI均优于7F-IMRT组,具有统计学差异(0.898±0.017 vs 0.834±0.013,P<0.001;0.062±0.012 vs 0.109±0.019,P<0.001);HT组PTV的D1、D2均低于7F-IMRT组(P<0.001),Dmean、D98、D95均较7F-IMRT组高(P=0.048,P<0.001,P=0.015).(2)对于正常组织,HT组小肠、结肠、直肠和膀胱的Dmax、Dmean、V30、V40以及直肠V50均显著低于7F-IMRT组(P<0.001,P<0.001,P<0.001,P=0.001,P=0.006),其余V50略低,无统计学差异(P=0.130);HT组双侧股骨头V25、D5、Dmean、Dmax也显著低于7F-IMRT组(P<0.001);但是,HT组的机器输出量明显高于7F-IMRT组.结论:HT计划在适形度、均匀性及周围正常组织保护上显著优于7F-IMRT计划,在充分保证靶区剂量的同时又不增加患者放疗毒副反应;但其机器输出量高于7F-IMRT计划,对射束的利用率不高.在后续治疗实施过程中,HT系统作为影像引导放射治疗技术,确保每次治疗过程中的精准度,值得临床应用.
-
碳素纤维床对宫颈癌容积弧形旋转调强放射治疗计划的剂量影响
目的:研究碳素纤维床对宫颈癌容积弧形旋转调强放射治疗(VMAT)计划的剂量影响及其修正方法.方法:使用CT电子密度模体校准大孔径CT的CT值,在Eclipse 10.0计划系统中建立相对电子密度-CT值曲线.将医用电子直线加速器上的碳素纤维床在CT下进行扫描,图像传输至Eclipse计划系统并测量碳素纤维床的CT值.以此CT值为基础在计划系统中建立碳素纤维床的模型,测量碳素纤维床在计划系统中和实际情况下对剂量的衰减系数并进行比较,两者测量条件保持一致.选取宫颈癌ⅠA、ⅡB期共5例患者,使用Eclipse计划系统设计无碳素纤维床计划.之后在患者定位图像上分别建立Thin、Medium和Thick这3种厚度的碳素纤维床模型,将无碳素纤维床治疗计划直接移植到3种不同厚度碳素纤维床图像上,治疗中心不变,机器条数不变,进行剂量计算.终比较有和无碳素纤维床之间、3种不同厚度碳素纤维床之间计划靶区(PTV)与危及器官的剂量差异.结果:Thin和Thick厚度碳素纤维床对实际测量和计划系统计算的剂量衰减系数相差均不超过±1％;不论是PTV还是危及器官,无碳素纤维床与3种不同厚度碳素纤维床的剂量参数结果比较均具有统计学意义(P<0.05);不同厚度碳素纤维床对PTV和危及器官剂量评价的准确性有一定影响.结论:加速器碳素纤维床对宫颈癌VMAT计划的剂量分布有一定影响;在进行VMAT治疗时,应准确建立碳素纤维床的模型参与剂量运算;并且根据靶区与碳素纤维床之间的位置关系选择添加相应厚度的碳素纤维床模型.
-
Geant4中不同物理列表对硼中子俘获治疗剂量分布的影响
Geant4是基于C++编写的开源蒙特卡洛模拟软件,提供多种包含中子与物质相互作用的物理列表.本研究采用Geant4提供的几种物理列表,对沿中子束方向的总吸收剂量、硼剂量和非硼剂量深度分布进行计算,并与FLUKA进行比较,模拟中子能量从0.0253 eV到10 MeV.对于整个模拟中子能段,结果显示添加S(α,β)热模型的高精度中子物理列表(Geant4_HP_T)在总吸收剂量、硼剂量和非硼剂量深度分布上均与FLUKA符合很好,初步验证了Genat4能应用于硼中子俘获治疗(BNCT)相关研究.对于低能中子(<1 MeV),S(α,β)热模型对BNCT剂量深度分布的影响较大,QBBC和QGSP_BERT不适用于BNCT剂量分布计算.
-
术前宫颈癌容积旋转调强与固定野动态调强技术剂量对比
目的:比较容积旋转调强(VMAT)和固定野动态调强放疗(DIMRT)技术在术前宫颈癌放疗中的剂量学差别.方法:选取10例接受放疗的术前宫颈癌患者,勾画靶区,用Eclipse 11.0计划系统设计双弧VMAT计划和5野DIMRT计划,比较两者靶区的适形度指数(CI)和均匀性指数(HI)、危及器官的剂量学差别、加速器跳数和照射时间.结果:双弧VMAT计划与5野DIMRT计划的CI(0.78±0.05,0.84±0.03;P>0.05)和HI(0.05±0.00,0.05±0.00;P>0.05)统计学无显著性差异,膀胱的V40、V45、V50和小肠的V20、股骨头、骨盆剂量有显著性差异之外,其它指标未有显著性差异.结论:VMAT和DIMRT技术两者均能很好达到临床剂量学的要求,但在保护正常器官方面,VMAT比DIMRT技术更优或相当,VMAT技术在减少加速器跳数和照射时间方面优势明显.
-
基于旋转准直器的锥形束CT散射矫正方法
针对锥形束CT(CBCT)图像质量受散射影响比较严重的情况,提出一种基于旋转准直器的CBCT散射矫正方法.该方法在射线源和模体之间放置一个圆形的旋转准直器,并通过准直器的旋转使透过准直器的射线不断沿轴向来回扫描,以获取整个容积图像的投影图像信息,然后利用投影图像的遮挡区域估计整幅图像的散射信息并将其从投影图像中去除,后利用改进FDK算法重建图像.结果表明,与CBCT图像相比,散射矫正后的重建图像CBCT值的均方根误差从16.00％下降为1.18％,杯状伪影从14.005％下降为0.660％,峰值信噪比从16.9594提高到31.4500.CBCT图像质量得到明显提高.
-
胸中段食管癌容积调强放疗计划的验证方法
目的:探讨胸中段食管癌在容积调强放疗(VMAT)中剂量与位置验证的方法.方法:随机挑选10例胸中段食管癌病例在Varian Eclipse 10.0计划系统(TPS)中制定VMAT计划,使用IBA Compass 3.0剂量验证系统进行剂量测量,然后与TPS计划数据进行比较,分析靶区(PTV、CTV与GTV)与危及器官受照剂量和体积参数的差异,并得到其γ通过率.应用锥形束CT(CBCT)验证放疗前摆位误差,1次/周,共6周.结果:γ分析在3 mm/3％标准下,靶区与危及器官通过率在95％以上.靶区D95％与Dmean的测量数据和TPS计算数据相差小于2％.危及器官中,双肺的V20与V30相比较于测量数据,TPS计算数据偏低,差异在1.65％以内.脊髓Dmax差异为2.23％,心脏V30、V40差异小于2％.CBCT位置验证中,前后与左右方向误差大于3 mm例数要多于头脚方向.结论:通过Compass 3.0剂量验证系统与CBCT扫描,是保证胸中段食管癌VMAT安全和可靠的必要手段.
-
基于空间频率与时间序列信息的多类运动想象脑电分类
结合共同空间模式(CSP)、离散小波变换(DWT)和长短期记忆网络(LSTM)方法,提出一种基于空间频率与时间序列信息的多类运动想象脑电特征提取方法.首先利用滑动矩形窗获得时间序列脑电信号,并采用DWT从每一段脑电信号提取运动想象脑电相关的子带小波系数,其次将小波系数通过一对多CSP进一步特征提取,得到的特征作为LSTM的输入,然后对LSTM的时间序列输出在时间步上进行平均,后使用Softmax分类器进行分类.实验结果显示,新算法取得92.23％的准确率,相比CSP特征以及结合频率或时间序列信息的CSP特征有较大提升,表明空间、频率、时间序列信息的互补性和有效性.
-
基于自适应学习的心律失常心拍分类方法
心律失常是因心脏疾病引起的心电活动中的异常症状,早期心室收缩(PVC)是由异位心跳引起的常见心律失常形式.通过心电图(ECG)信号检测PVC对于预测可能的心力衰竭具有重要意义.本文提出一种面向PVC心拍分类的心电信号分类算法,重点研究基于自适应学习的PVC异常心拍分类特征提取模型,通过计算心拍关联后验概率,结合领域专家标注信息训练分类器,提高整体分类效果.实验采用MIT-BIH心律失常数据库的ECG数据,研究结果表明所提方法针对非线性流形结构数据,能够有效提升小样本心拍自适应分类器的准确性.
-
轴流式血泵无位置传感驱动控制系统设计
目的:设计一种无位置传感器轴流式血泵驱动控制系统,实现轴流式血泵的速度闭环控制、无线状态监测及动态调节控制.方法:以STM32F103作为主控芯片,设计六臂全桥驱动电路,采用端电压过零检测实现转子位置和速度检测,并设计PI控制器实现闭环控制,使用无线蓝牙通讯实现泵状态监测及控制.结果:驱动控制系统速度响应快,稳态误差小于1％,在8000 rpm转速下,能够达到人体血液循环所需的流量及压差要求.结论:系统实现了轴流式血泵无位置传感的启动及速度闭环控制,通过蓝牙模块实现对泵运行过程的电压、电流及流量进行监控和运行过程中根据病人生理情况动态调节泵转速.
-
青年与老年昼夜心率信号多尺度化的基本尺度熵分析
目的:本研究采用多尺度化的基本尺度熵方法分析青年组与老年组,白天和黑夜的心率变异性信号.方法:从PhysioBank数据库中,分别提取出青年人和老年人的心率变异性信号,利用多尺度化的基本尺度熵方法对数据进行比较分析.结果:通过计算心率变异性信号的熵值,发现在白天清醒状态下的熵值比睡眠状态下的熵值高,但是昼夜间变化相似;老年组的多尺度化的基本尺度熵相对于青年组的熵值有所偏离,但大致走势一样,可见老年组的生理状态偏离了青年组的佳生理状态.结论:通过熵值的分析,揭示出清醒状态下,心脏系统的自适应性和稳定性较睡眠状态要强;青年与老年的熵值具有昼夜节律相关性;多尺度化的基本尺度熵方法可以用于区别青年与老年,白天与黑夜的心率变异性信号.
-
用于康复治疗的磁振热理疗仪设计
开发一种用于康复治疗的磁振热理疗仪,该理疗仪包括磁热振子与控制系统两部分,磁热振子提供磁、振、热3种物理治疗因子,控制系统实现对时间设定与治疗温度实时监控,能够为患者提供磁场强度、治疗温度、振动按摩等参数多种选择.通过对该理疗仪的磁热振子进行实验,结果显示其含芯线圈的温度在通电5 min后可以达到稳定,且可以通过设置工作电源电压对温度进行精确控制.实验结果证明,该磁振热理疗仪具有多磁场强度、温度独立设置、两种按摩方式等特点,能够很好满足临床上不同患者的需求.
-
振动疗法在康复中的应用及其影响因素
振动疗法利用振动因素使人体局部组织或全身发生周期性或非周期性的变化,从而产生治疗效果,其作为一种物理疗法在康复领域中发挥着重要的作用.应用振动疗法时,主要的影响因素是频率参数的选择,认为存在一种佳振动频率(谐振)是提高疾病康复效果的有效方法.应用佳振动频率刺激作用于机体使之能够在短时间内对神经肌肉性能产生影响,并对其进行时-量-效标准化分析以进一步优化振动疗法,这具有进一步深入研究的意义.本文将从振动疗法的定义、起源、作用原理、在康复治疗中的应用、影响因素以及未来振动疗法的优化等方面进行综述,并提出关于谐振的初步探讨.
-
支气管热成形术治疗难治性哮喘研究进展
难治性哮喘是指大剂量吸入激素(ICS)或口服ICS仍难以控制的一类哮喘,因其长期反复发作,严重影响患者的生活质量,并导致大量的医疗资源消耗.支气管热成形术(BT)作为针对难治性哮喘的一种新型治疗方法,主要通过射频消融减少气道平滑肌数量而发挥作用.如今,越来越多的数据表明BT可有效控制难治性哮喘患者的哮喘症状、降低急性发作率、提高患者生活质量.本文主要对BT的临床应用、影像学在BT中的使用以及BT对气道病理学的影响等进行综述,为临床医生更好选择BT治疗患者、评估患者病情提供帮助.
-
极低频电磁暴露对钙离子跨膜迁移的影响及机理分析
目的:研究极低频(ELF)电磁暴露对细胞钙离子跨膜迁移的影响.方法:以人体肝癌细胞为对象,激光扫描共聚焦显微镜为检测手段,研究细胞钙离子跨膜迁移对ELF电磁暴露的响应.结果:在磁场强度为1.78×10-7 T条件下,与对照组对比,各暴露组均对不同频率和电场强度组合的外加电磁暴露环境有响应,其中f=16 Hz且EP=53 V m、f=45 Hz且EP=53 V m以及f=16 Hz且EP=80 V m的电磁暴露可使钙离子跨膜迁移量显著上升,而f=32 Hz且EP=53 V m、f=60 Hz且EP=53 V m、f=16 Hz且EP=26 V m以及f=16 Hz且EP=87 V m的电磁暴露未使钙离子跨膜迁移量显著上升.结论:不同频率和电场强度组合的ELF电磁暴露均会引起钙离子跨膜迁移量产生变化且呈现出明显差异,这可为ELF电磁环境下生物学应用提供实验依据;注重细胞膜离子通道物理、生物特性基础上探讨ELF电磁暴露对钙离子跨膜迁移影响的机理.
-
近似熵在脑电监测麻醉深度中的应用
目的:探究脑电近似熵在麻醉深度监测中的应用.方法:以2017年1月～2018年1月择期行全麻手术的80例患者为研究对象,测定所有患者麻醉前、麻醉诱导5 min、术中60 min、麻醉苏醒时脑电双频指数(BIS)、脑电近似熵,且观察不同丙泊酚血浆浓度、改良镇静/警觉评分(MOAA/S)BIS值、脑电近似熵值的变化情况,Pearson相关性分析BIS、脑电近似熵值与丙泊酚血浆浓度的相关性.结果:麻醉诱导5 min、术中60 min BIS值、脑电近似熵值比麻醉前均显著降低(P<0.05),且术中60 min的值显著低于麻醉诱导5 min时的值(P<0.05);BIS值、脑电近似熵值随丙泊酚血浆浓度上升而下降,均呈显著负相关(P<0.05);脑电近似熵在MOAA/S评分0～1分变化中有显著差异(P<0.05),而BIS值无显著差异(P>0.05).结论:脑电近似熵对丙泊酚药效变化评价效果与BIS相当,但其相对在麻醉镇静深度判断上有优势.
-
3D打印技术在脊柱外科中的应用进展
总结国内外脊柱外科,特别是在复杂脊柱疾病中(尤其是严重脊柱畸形和脊柱肿瘤)使用3D打印技术的文献,重点关注提到的3D打印种类、具体应用方法、治疗效果以及成本效益分析.对在脊柱外科手术中使用3D打印技术的文献进行电子检索,国内文献使用中国知网数据库进行检索,国外文献使用PubMed以及Web of Science数据库进行检索.检索词为:脊柱外科、3D打印技术、增材制造、快速成型、脊柱畸形、脊柱肿瘤;英文词为:Spine、Scoliosis、Spinal Tumor、3D printing、Additive manufacturing、Rapid prototyping.并对纳入的文献进行分类总结,检索到文献116篇(中文72篇,英文44篇),进行全文阅读分析后纳入文献35篇(中文15篇,英文20篇).根据文献内容将3D打印应用分为手术导板、解剖模型以及个体化植入物.纳入文献提到的优点主要有减少手术时间、出血量、术中照射次数以及手术并发症发生率、改善治疗效果、提高患者满意度等.主要局限在于该技术目前缺乏统一标准、额外的打印费用、增加的学习成本、设计3D打印物耗费的时间以及植入物远期可靠性等.3D打印技术在脊柱疾病的治疗中确有价值,但在常见脊柱疾病中使用3D打印技术的意义和价值并不大,只有在较为复杂的病例中使用这项技术才能符合成本效益原则.
-
基于弥散张量白质网络的阿尔茨海默病研究
目的:利用弥散张量成像(DTI)技术构建大脑白质结构网络,通过网络拓扑参数研究阿尔茨海默病(AD)患者大脑微观结构的异常.方法:重建19例正常人和20例AD患者的结构脑网络,采用双样本t检验,从全脑和特定脑区两个水平对AD患者与正常对照组的网络拓扑参数进行差异性分析.结果:在全脑水平上,AD患者的加权特征路径长度(Lp)值上升,全局效率(Eg)、网络强度(Sp)、局部效率(Elocal)值下降,并且进一步研究发现:AD患者左半球Lp、加权簇系数(Cp)值上升,Eg值下降比右半球明显.在特定脑区水平上,选取楔前叶为感兴趣脑区,发现两侧楔前叶脑区的Lp值上升、Eg值下降,而Cp、Elocal值并没有明显差异.结论:网络拓扑参数可以作为评估AD患者微观结构异常的指标,对早期诊断AD具有重要的指向作用.
-
核磁共振超短回波时间序列技术研究进展
核磁共振成像(MRI)技术拥有良好的软组织分辨率且无电离辐射,在临床和科研方面均得到了广泛应用.超短回波时间序列(UTE)在一定程度上弥补了MRI在短T2组织成像的弱点,使MRI的应用更加广泛.由UTE得到的派生序列有脂肪抑制UTE、单绝热反转恢复UTE、双回波差UTE等.本文介绍核磁共振超短时间回波序列(MR-UTE)技术的发展、原理及其应用,并对MR-UTE技术的发展方向进行展望.
-
基于纤维束空间统计学的吸烟成瘾青少年左右脑白质不对称研究
目的:观察吸烟成瘾青少年左右脑白质结构的不对称性变化.方法:采用基于纤维束空间统计(TBSS)的分析方法,对23名青少年吸烟成瘾者的脑白质弥散张量成像(DTI)数据进行处理,观察分析出现左右脑白质结构异常不对称的脑区.结果:通过分析,我们发现青少年成瘾者的脑白质在胼胝体部、胼胝体膝部、胼胝体压部、前放射冠、内囊前肢等脑区处的左右脑白质不对称差异值较大(FWE校正后P<0.05).Pearson相关分析显示在胼胝体部(r=0.431,P=0.040)、胼胝体膝部(r=0.429,P=0.041)、前放射冠(r=0.477,P=0.021)这些脑区的轴向弥散率(AD)与尼古丁依赖量表(FTND)程度分数呈显著正相关,胼胝体压部与FTND呈显著负相关(r=-0.450,P=0.031).胼胝体膝部的分数各向异性与FTND呈现显著正相关(r=0.678,P=0.000).前放射冠的AD与吸烟指数包年(pack-year)呈现显著正相关(r=0.493,P=0.017).结论:青少年吸烟成瘾者的胼胝体,前放射冠等脑区的脑白质表现出了异常的左右不对称.
-
3D打印技术在骨科关节中的应用进展
随着新材料技术的不断发展,3D打印技术在骨科关节领域中的应用越来越广泛.本文综述了3D打印技术在骨科关节中的研究进展,分别介绍了3D打印在解剖学教育、医生和患者及其家属的术前沟通、临床治疗方案设计的应用,并对关节置换和骨关节软骨修复在组织工程领域的进展进行了阐述.后,对3D打印技术在该领域未来的机遇和挑战进行了预测.
| 年 | 期数 |
| 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 |
