中国生物医学工程学报杂志
Chinese Journal of Biomedical Engineering 중국생물의학공정학보
- 主管单位： 中国科学技术协会
- 主办单位： 中国生物医学工程学会
- 影响因子： 0.61
- 审稿时间： 1-3个月
- 国际刊号： 0258-8021
- 国内刊号： 11-2057/R
- 论文标题 期刊级别 审稿状态
-
基于代价函数掩模法的局部脑损伤fMRI数据处理
功能磁共振数据空间预处理中,个体图像到标准空间的归一化一般采用基于强度小化的自动归一化标准方法.在处理有明显脑结构局灶异常的数据时,由于病灶区与模板的不匹配,该算法会导致图像的严重扭曲和变形,造成数据失真和功能激活区定位错误.采用代价函数掩模方法,归一化处理有明显脑结构局灶异常功能数据,可减少传统归一化方法因脑损伤而扭曲、变形对数据造成的影响.本研究对4例顶叶主运动区有明显脑结构异常(2例脑膜瘤,2例脑胶质瘤),手动功能评价的fMRI砌数据进行处理.首先采用代价函数掩模法对病灶区进行掩盖,然后再进行标准模板的归一化处理,归一化后的结果用SPM软件分析脑激活区.实验表明,基于代价函数掩模法的归一化方法可以有效地减少局灶脑结构异常fMRI数据归一化过程中的扭曲,提高脑功能区定位的准确性,数据处理的结果更符合脑病理情况下脑区活动的真实情况.
-
基于3D医学图像的血管三维分割
血管三维分割在血管疾病(如狭窄或畸形)诊断、手术规划和手术引导等许多实际应用中发挥重要作用.但三维分割的实时性仍是一个难题.本研究提出一种基于水平集的快速三维血管分割方法,该方法用内、外邻域曲面来描述被分割目标的边界,并定义水平集函数为简单的整数符号距离函数.通过扫描内外邻域曲面上的点,使之在速度场的作用下向目标边界移动.该方法的不同之处在于利用简单的模型极大地减小了计算量,分割速度快,大容量的MS-CTA[256<'3>体素]图像可在20s内处理完毕.同时,展示了一些三维血管的分割实例.特别需要指出的是,对于人体的大血管,可在不经血管造影的情况下,直接从CT等三维图像中分割出来.
-
基于经验模式分解和独立分量分析的单导少次EP信号提取
EP信号的单导少次提取一直是生物医学信号处理领域倍受关注的问题.本研究利用经验模式分解(EMD),把单导脑电信号(EP+EEG)分解成多个基本模式分量(IMF)之和,进而选取合适的基本模式分量或者它们的组合,构成1导或多导参考信号,再利用独立分量分析(ICA)成功提取出了期望的EP信号,从而克服了ICA需要多导观测信号的要求.仿真实验证明了本方法的有效性.
-
基于PCA-LDA的HIFU治疗中组织损伤无损检测方法
组织损伤的无损检测是HIFU治疗过程中非常重要的辅助过程.提出了一种通过B超图像无损地检测组织损伤的方法.HIFU辐射后,组织损伤必然会引起B超图像的变化,只对HIFU焦域附近的数字减影图像进行处理和分析.利用新鲜离体猪肉组织和猪肝组织实验,利用主分量分析(PCA)和线性判别分析(LDA)相结合的方法,对减影图像进行特征提取,利用近邻方法进行分类识别.对于离体猪肉组织中损伤的识别率达94.5%,对于离体猪肝组织中组织损伤的识别率达87.5%.PCA-LDA方法的识别率要优于平均灰度作为识别特征等利用单一特征的识别方法.
-
基于模拟座舱的驾驶疲劳实验研究
通过模拟实验研究影响驾驶疲劳的典型因素并探讨驾驶疲劳的评测方法.以模拟驾驶器为核心,结合操控能力、生理,心理测试、主观评价和行为特征,构建了模拟驾驶综合实验系统.设计了一套合理的实验方案,并开展了不同道路状况下长时间模拟驾驶对疲劳的影响、驾驶员不同作业时段的复杂反应、多导脑电的复杂性与驾驶进程关系等实验.实验显示,长时间驾驶对操控能力、反应时间、心率、主观疲劳都有显著性影响;高速公路环境下车速方差、主观疲劳比一般公路区别显著;驾驶员的心理疲劳程度受时段的影响较大;脑电的复杂度反映了驾驶疲劳的变化.研究结果为安全驾驶提供了有意义的指导,同时也验证了实验系统和实验方案的合理性,为进一步的研究工作打下了基础.
-
搜索空间逐步缩小的遗传算法用于尿沉渣图像特征选择的研究
基于尿沉渣图像特征选择问题,提出一种新的特征优选方法,首先引入类空间分层分类思想,将多类成分特征集优选问题转化为两类成分的特征集优选问题以减少特征数,从而减少了后续分类器的维数复杂度并提高了优选后的特征集对成分的识别率;针对红细胞和白细胞的特征集优选问题,采用改进的遗传算法进行处理,先根据统计实验结果,锁定待选特征集中形态特征和纹理特征相应的两个可分度大且相互独立的特征,然后使用基因位逐步锁定技术,结合小生境技术和自适应交叉变异算子,提高了遗传算法的搜索性能;后,为了提高特征集的优选效果和稳定性,引入"多票投选"机制.就多个尿沉渣成分样本进行验证实验,结果表明,该算法优选的特征集与通过其他方式获得的特征集相比,识别率较高,而且明显减少了后续分类器的维数复杂度.
-
改进的熵小方法用于MRI偏差场的校正
灰度不均匀性分布是磁共振成像中经常发生的一种病态混扰,它给医生读片诊断和磁共振图像(MRI)计算机自动化处理带来很大困难.针对偏差场平滑、渐变的特性,本研究提出了一种改进的基于图像信息熵小化的自动校正算法.该方法是一种全自动的无需先验知识的校正方法.在实验过程中,采用模糊C-均值聚类算法来提取反映偏差场变化趋势的同类特征点;通过二次曲面拟合特征点曲面的方法,提取初始化曲面参数;应用具有学习能力的粒子群寻优算法提高信息熵搜索的优化速度.后,在对Brain-Web数据库仿真MRI和实际采集到的MRI校正中,获得了更适合人眼观看的校正图片,为临床诊断提供了准确可信的图像数据.
-
ECM支架构建肺成纤维细胞三维培养模型研究
应用组织工程技术建立新型稳定、可靠的肺成纤维细胞(LF)三维培养模型的方法.应用细胞外基质(ECM)支架为成纤维细胞培养载体,建立成纤维细胞三维培养体系;并同时分别加入IL-13单克隆抗体、TGF-β及LPS.应用细胞荧光免疫、细胞计数及扫描电镜观察肺成纤维细胞三维培养的生长特性及相关物质对细胞生长状况的影响.在此三维培养模型中,细胞、细胞外基质以及IL-13单克隆抗体、TGF-β等调控因素共同构成一个有机整体,成纤维细胞所表达的生物学特性接近于它们在体内肺组织中的特点.本实验应用组织工程技术成功地建立了新型肺成纤维细胞三维培养体系,并研究成纤维细胞的生物学特性,是研究间质性肺病(ILD)领域里一个新的和更为科学的方法.
-
动物血管脱细胞方法及细胞外基质材料评价研究
采用独立设计的反复冻融方法,去除兔颈总动脉中的血管细胞,对所获得的细胞外基质进行组织、生化、力学等分析.在分离兔颈总动脉后,首先用低渗缓冲液处理,然后经低温反复冻融,后用非离子型去污剂处理.所得的支架行组织染色和扫描电镜分析,显示基质的微观结构;荧光染色和基因组DNA PCR分析,检测细胞DNA残留;分别行羟脯氨酸定量分析和轴向拉伸强度分析,检测胶原蛋白和血管生物力学性能的变化;支架没有明显的细胞毒性,溶血率低于医用材料的标准;支架随时间缓慢降解.种植兔的骨髓干细胞,研究细胞在支架上的粘附生长能力.应用本方法能彻底去除血管细胞,没有细胞碎片和DNA的残留,并保留了较完整的细胞外基质和力学性能,具有开放的大孔径结构,细胞生物相容性好,是一种理想的组织工程血管支架材料.
-
功能电刺激腕关节运动的自适应控制技术研究
针对功能电刺激(FES)康复技术中肌肉的时变和高度非线性的特点,提出一种实时调整FES系统输出,实现康复患者自适应训练的方法.考虑到患者问和康复各阶段肌肉对电刺激响应能力的差异,利用神经网络良好的非线性逼近能力,采用神经网络辨识器辨识相关运动肌肉的FES模型,并在此基础上,利用神经网络控制器在线调节控制策略及刺激参数,提高FES系统的自适应能力.通过人体腕关节运动控制的FES实验,验证该方法具有较高的轨迹跟踪控制精度和较强的自适应能力.
-
骨髓间充质干细胞与不同基质修饰的纳米晶胶原基骨体外生物相容性研究
研究大鼠骨髓间充质干细胞(MSCs)诱导培养后与不同基质修饰的纳米晶胶原基骨(nanoHydroxyapatite/collagen,nHAC)的生物相容性,为骨组织工程提供一种新型复合支架材料.SD大鼠MSCs经成骨诱导培养、扩增,进行成骨细胞表征后,种植与支架材料体外复合培养.实验分为4组:实验组A,纤维蛋白(FB)和纤维连接蛋白(FN)修饰的纳米晶胶原基骨((FB+FN)-nHAC);实验组B,纤维蛋白修饰的纳米晶胶原基骨(FB-nHAC);实验组C,纤维连接蛋白修饰的纳米晶胶原基骨(FN-nHAC);对照组D,单纯的纳米晶胶原基骨(nHAC).通过检测支架材料的细胞黏附率、不同时间点(3、7、10、14d)支架材料中细胞数、碱性磷酸酶活性以及扫描电镜观察细胞在材料上的生长状况,比较分析不同支架材料与细胞生物相容性差异.大鼠MSCs经诱导培养14d后,碱性磷酸酶细胞化学染色、Ⅰ型胶原免疫荧光染色及矿化沉积茜素红染色均为阳性;细胞与支架材料黏附率A组高为74.4%;支架材料中细胞数量均随培养时间延长而增长,且A组细胞数增加较快.与相同时相点其他各组材料中细胞数差异有显著性(P<0.05);各时相点细胞碱性磷酸酶活性表达A组高,差异亦有显著性(P<0.05).电镜观察发现4组材料上均有细胞生长,但A组的细胞生长状况明显好于其他组.大鼠MSCs经成骨诱导培养,可表达成骨细胞表型,(FB+FN)-nHAC在体外实验中表现出与细胞优良的生物相容性,可作为较理想的新型复合支架应用于骨组织工程.
-
多孔HA/PU复合支架材料生物相容性的评估
进行了三维多孔立体结构的纳米羟基磷灰石/聚氨酯(HMPU)复合支架材料体外细胞培养和体内肌肉埋植实验研究,评估材料的生物相容性.实验选用SD大鼠的骨髓基质干细胞(BMSCs)和健康的SD雌性大鼠,进行细胞相容性、形态学观察和组织学切片分析.HA/PU支架材料的多孔性为细胞的生长提供了良好的微环境,细胞在内部贴壁爬行、增殖并分化,细胞毒性为零级,材料与周围组织有良好的结合,降解的空间有结缔组织纤维长入.实验表明,HA/PU复合支架材料具有良好的细胞亲和性和组织学相容性,可作为一类新型组织工程支架材料.
-
具有生长和衰老机制的生物神经网络的兴奋特性研究
在具有侧向抑制机制的小世界人工生物神经网络中,对神经元的放电特性引入增长和衰老机制,构建具有模仿实际生物神经网络生长和衰老功能的网络模型.为了研究该网络模型表现出来的发放电的特性与实际生物神经网络行为之间的一些关系,本研究数值仿真分析了该网络模型在不同噪声强度和不同耦合强度下的兴奋特性,发现模拟结果与实际生物神经网络在神经细胞处于生长发育或衰老时,对外界刺激所表现出的兴奋效应及其饱和效应有一定程度上的相似性,研究结果对探索实际生物神经网络的兴奋特性和放电节律具有较好的参考价值.
-
改进的遗传模糊聚类算法及其在MR脑组织分割中的应用
为提高MR图像中脑组织分割的精度,针对目前遗传模糊聚类算法存在的问题,提出了改进的遗传模糊聚类算法.首先,通过完全改变遗传算法的编码方式、变异方式和交叉方式,对现有遗传算法进行改进,从而使遗传算法能获得小的适应度函数值;然后,结合模糊聚类算法,终得到改进的遗传模糊聚类算法.将改进的遗传模糊聚类算法应用于MR脑图像的分割,结果表明,改进的遗传模糊聚类算法的分割质量高于现有的遗传模糊聚类算法和快速模糊聚类算法.改进的遗传模糊聚类算法可以做为一种快速、全自动的MR脑图像分割工具.
-
基于相关分析的fMRI时间序列图像处理
脑与认知科学研究中,对功能磁共振图像的分析,不仅要求脑功能激活区的准确定位,而且要得到脑激活区的动态变化.提出了通过对激活区体元时间序列的分解,构造序列参考波形,用相关分析法检测识别每个时间周期的激活类型,得到脑功能激活区的动态变化图像、生理信号变化时间及信号强度曲线.经视觉试验实测数据建模并检验,该方法可有效检测fMRI数据中与试验密切相关的脑激活区及其周期性变化.
-
心房活动电生理信号特征提取
为从12导联体表心电图中有效地提取反映心房活动的电生理信号,构造一种稀疏表示下的独立成分分析模型.首先,对采集到的12导联体表心电信号进行预处理,然后通过小波变换,实现心电信号小波域的稀疏表示形式,对变换后的心电信号进行独立成分分析,并通过频谱分析确定出反映心房活动的电生理信号源.讨论了平均房颤周长的这一测量指标在房颤患者中的具体应用.实验表明,该方法可以有效提取心房活动的电生理信号,将对心房活动电生理的深入研究产生积极影响.
-
基于CSSD和SVM的抑郁症脑电信号分类
从EEG脑电信号中提取与疾病相关的信息以实现对抑郁症的自动诊断.首先采用共空域子空间分解(CsSD)方法,对躁狂型抑郁症患者与健康人两组的16导联脑电信号进行特征提取,然后用支持向量机(SVM)分类器进行训练和分类测试.实验结果表明,相对于用小波变换提取的频率相关参数为分类特征的分类准确率为88%,采用CSSD方法提取特征参数进行分类可以取得更理想的效果为95%,后者的16导联脑电信号在空间模型上表现出较高的模式可分性.该研究成果对精神抑郁症的物理诊断和研究提供了新的视角.
-
冲击过程中脑组织力学行为的数值模拟
通过有限元法,计算不同加速度载荷下人体脑部颅压变化,并将计算结果与头部损伤准则(HIC)标准进行对比.使用超弹+黏弹本构关系模拟脑组织行为,分析颅骨-脑组织系统在冲击载荷下的力学响应,并对比重力、脑脊液压力等对脑组织力学行为模拟的影响.结果发现,考虑重力与否(重力方向与速度和加速度方向垂直)对求解的影响相对较小,考虑重力因素会使每个颅压波动周期内的小值略有增加.在加速度载荷施加的时间范围内,第一个颅内压力波动周期内的颅压的大值较大,并在随后的波动周期内逐渐降低.对于人体脑部组织保护而言,考虑脑脊液压力的计算更靠近安全设计.对于不同的加速度载荷,HIC标准在判断人体脑部组织损伤时具有不一致的效果.在汽车(机车)的高速制动(加速)过程中,减速(加速)的初始阶段对人体脑组织的危害性大.
-
人工加载对悬吊大鼠骨组织形态计量学影响的实验研究
为探讨力负荷对废用性骨质疏松的影响,针对人工加载对悬吊去负荷SD大鼠胫骨松质骨组织计量学的影响进行了研究.选取9月龄清洁级雌性SD大鼠21只,分为对照组、悬吊组、悬吊加载组,采用自行研制的动态活体加载装置每日对悬吊加载组大鼠加载15rain,频率选取5Hz,外部载荷选取3N.4周后处死大鼠,取左侧胫骨,采用显微CT扫描标本,进行胫骨组织计量学测定分析.结果表明悬吊组大鼠松质骨组织形态计量学参数与对照组及悬吊加载组相比显著改变;对照组与悬吊加载组松质骨组织形态计量学参数无显著差异.说明人工加载能显著改善悬吊大鼠松质骨形态计量学参数.
-
聚焦超声结合血卟啉对S180肿瘤细胞的杀伤作用
初步研究聚焦超声激活血卟啉对S180肿瘤细胞的杀伤作用,并探讨其作用的生物学机制.通过台盼蓝拒染法检测频率为2.2MHz的聚焦超声激活血卟啉作用于S180肿瘤细胞的存活率,筛选出佳的超声处理参数,采用自由基清除剂检测超声激活血卟啉产生的起主要作用的活性氧成分,通过丙二醛含量测定法检测细胞内脂质过氧化含量变化.单纯血卟啉组表现出无明显的细胞毒作用;同等条件下,超声结合血卟啉组对细胞的损伤显著高于单纯超声组,且这种协同杀伤作用可以被组氨酸和甘露醇所抑制,而不被SOD抑制;同时,超声激活血卟啉后细胞内的脂质过氧化含量显著增加.超声结合血卟啉对S180肿瘤细胞有协同杀伤作用,可能由于声动力学处理过程中产生的活性氧自由基作用于细胞膜,增加膜脂质过氧化水平,终导致细胞死亡.
-
循环吸气支气管内瞬态气流运动特性的仿真研究
人体支气管内的瞬态气流运动对气溶胶在支气管内的沉积产生重要影响,深入认识支气管内的瞬态气流运动特性,对于分析气溶胶在人体支气管内的沉积模式具有重要作用.应用计算流体动力学方法,对人体轻度运动强度下Weibel模型的3级到6级支气管内的瞬态气流运动特性进行了数值仿真,讨论了支气管对称平面、横截面内的瞬态气流组织形式.研究结果表明:在支气管的对称平面内,支气管分叉处产生气流分离现象;在吸气加速阶段,在3级支气管横截面内,没有产生涡流运动;在吸气加速阶段,4级支气管内,呈现出具有内壁附近区域速度高、外壁区域附近速度低的倾斜的抛物线特点的速度分布形式;5级支气管内,形成由两个涡流到四个涡流又到两个涡流的涡流演变过程.在6级支气管横截面内,在吸气加速阶段,和4级支气管横截面相似,在对称平面上下形成两个对称旋转的涡流运动,两个涡流中心分别处于上下平面的中心位置.
-
表面改性人工机械心脏瓣膜的动物体内实验研究
通过人工瓣膜临床前羊慢性存活实验,探讨新型国产表面改性人工心脏瓣膜的抗凝血性能和综合性能.在19例羊心脏内植入新型TiO2-x/TiN薄膜材料制成的国产双叶人工瓣膜,与植入热解碳制成的St Jude双叶瓣进行对比,分析凝血指标、血流动力学测定、心脏超声检查和长期生存动物形态学及病理组织学检查,对该瓣膜的综合性能做临床前动物体内评价.成活15例,4例死亡,其中对照组8例、实验组7例,长存活9个月且是实验组.两组凝血指标在围术期及术后24h内相似,无显著性差异(P>0.05);对照组术后1～20周血凝血酶时间均较实验组高,有显著性差异(P<0.05),考虑与华法林抗凝治疗有关,而实验组在术后1周以后恢复术前水平.两组血流动力学围术期变化相似,心脏超声检查及术后半年尸检示两组各瓣膜活动好,均未见血形成,无赘生物,各主要脏器均无栓塞等异常改变.该新型国产人工瓣膜与进口机械瓣有相似的血流动力学功能,获得了满意的长期生存效果,且抗凝血性能明显优于进口机械瓣.
-
不同方法处理牛心包的体外力学和抗钙化性能评价
比较传统GA法(A组)、三氯化铝与乙醇联合处理法(B组)、多聚环氧化物处理法(c组)和光氧化法(D组)四种不同方法处理牛心包后力学与抗钙化性能.对新鲜牛心包随机分为四组,用四种方法处理后,测定组织的极限抗拉强度、断裂伸长率、热皱缩温度,比较其力学性能;建立大鼠皮下植入模型,包埋2月后取出标本测定组织钙含量,并做HE染色及von Kossa染色.极限抗拉强度:B组>A组=C组>D组;断裂伸长率:C组=D组>B组>A组;热皱缩温度:A组=B组>C组>D组;钙含量测定:A组>>C组>B组>D组(P<0.05);HE染色和vonKossa染色结果符合以上测定值.三氯化铝与乙醇联合处理组的组织交联程度与力学强度佳,光氧化处理组与多聚环氧化物(简称PC)处理组的组织柔韧性更好.在抗钙化性能方面,光氧化处理组好,三氯化铝与乙醇联合处理组与多聚环氧化物处理组次之,GA处理组差.
-
经颅磁刺激对握力影响机理的初步研究
目的:探讨在不同刺激强度和不同手部握力输出的条件下,磁刺激前后对握力大小的改变.方法:10名健康受试者,其右手分别以3个不同的背景力量握住握力计,用圆形磁刺激线圈刺激左侧运动皮层,施加的刺激强度从阈值开始,每次递增10%,记录磁刺激施加前后握力输出的变化.结果:对于同一背景握力,刺激强度越大,握力大小的峰峰值及负峰值也越大;而在刺激强度不变的情况下,握力波形的峰峰值和负峰值也随握力的增大而增大;输出握力波形的正峰值和恢复时间与刺激强度、背景握力大小无明显关系.结论:实验结果提示握力峰峰值及负峰值的变化可能与参与肌肉收缩的运动单位数量及其兴奋性有关.
-
基于低频超声和真空负压的微创血糖检测研究
探讨基于低频超声和真空负压微创抽取组织液,预测血糖浓度,实现血糖的无痛、微创、连续检测的可行性.对3位青年健康志愿者的前臂内侧皮肤,采用55kHz的低频超声破坏角质层的脂质双层结构,提高皮肤的通透性,然后采用254mmHg的真空负压透皮抽取组织液,通过检测组织液中的葡萄糖浓度,利用单点校正法预测血糖浓度,后评价系统的检测精度.预测血糖浓度基本能够满足血糖检测标准,有望实现血糖的无痛、微创、连续检测.
-
实验条件下温热环境的热舒适期望
为考察人对温热环境的热感受,在相对湿度适中的条件下,综合考虑温度、平均辐射温度及风速对静坐状态人体的热舒适影响,给出实验条件下温热环境的热舒适期望(TCE)方程,进一步分析人体热舒适与人体生理的联系.结果表明:在实验条件下,温度和风速对人体热舒适有显著影响,辐射温度对热舒适的影响相对较小.该TCE方程组合多个环境热物性参数评价人体的热舒适性,与实验测试中的人体生理与热舒适投票有较高复核度.认为TCE方程可用于实验条件下温热环境的热舒适预测,为现场热舒适评估提供指导,同时也可为非实验条件下的现场热舒适评估提供参考.
-
语言文字加工神经机制的fMRI研究进展
应用血氧合水平依赖功能磁共振成像技术(BOLD-fMRI)进行语言文字加工研究发现,大脑左侧梭状回中部皮质对语言文字视觉加工的任务具有稳定和明显的反应,但是在语言文字加工中该区域的激活是由何种因素引起尚存在不一致的看法.一些研究认为该区域的激活主要由于文字的构词规则、语音和语义的加工等因素所致;但是另外一些研究则认为该区域是通过对语言文字的广泛经验和专业知识构建了专门处理文字视觉加工的皮层模块.作者对近的相关研究进展进行了综述和分析,同时对今后可能的研究方向进行了探讨.
-
生物大分子的电纺及其在药物控释中的应用
随着研究者对于纳米技术以及纳米特性研究兴趣的不断增加,近年来静电纺丝受到越来越多的关注.首先综述了生物大分子静电纺丝的研究进展,主要包括合成高分子和天然生物大分子的静电纺丝以及复合纳米纤维的制备,侧重讨论了各种生物材料的电纺条件以及与生物医用有关的电纺纤维膜性质,如纤维膜形态稳定性、机械性质以及降解行为等.在此基础上,综述了静电纺丝在药物控制释放领域中的应用,针对不同药物的溶解性质,分别讨论了其在电纺纤维中的包埋方法以及释放行为.
| 年 | 期数 |
| 2018 | 01 02 03 04 05 06 |
| 2017 | 01 02 03 04 05 06 |
| 2016 | 01 02 03 04 05 06 |
| 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 05 06 |
| 2002 | 01 02 03 04 05 06 |
| 2001 | 01 02 03 04 05 06 |
| 2000 | 01 02 03 04 |
| 1999 | 01 02 |
| 1998 | 01 02 |
| 1994 | 04 |
