资源描述:
《ct图像中轮廓特征的提取new》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、研究论文医疗设备信息CT图像中轮廓特征的提取陈欠根,陈大舵(中南大学机电工程学院,湖南长沙410083)[摘要]本文在分析CT图片的特点的基础上,采用经典的Roberts算子对其进行边缘检测,并对其进行了阈值分割、细化、最后采用轮廓跟踪的方法,在VC6.0的基础上设计了具有选择性的提取出CT图片的外部或内部软组织轮廓的方法。[关键词]边缘检测;轮廓提取;轮廓跟踪[中图分类号]TP391.41[文献标识码]A[文章编号]1007-7510(2003)01-0010-03Thecontourextr
2、actionofCTimageCHENQian-gen,CHENDa-duo(ElectricEngineeringCollegeofSouth-CentralUniversity,WuhanHubei410083,China)Abstract:Thepapermainlyintroducesawayofcontourextraction.FirstRobertsoperatorisusedtoedgede2tection,thenathresholdisusedtodividethepixel
3、valueintotwopartsthatpredigestthecalculation,finallycontourisselectedbyatracearithmetic.ThismethodhasbeenrealizedonVC6.0andabstractedsuccessfullythedifferentcontoursofCTimages.Keywords:edgedetection;contourextraction;contourtrace[2]0引言息。而将多种成像模式的图像结合
4、起来,利用随着科学技术的发展,医学图像种类越来越各自的优点,可以在一幅图像上同时表达来自多多,而且在临床诊断方面也发挥着越来越重要的幅图像上的信息(如融合MRI与PET、CT和作用。目前,医疗图像的种类有计算机断层成像MRI等医学图像进行融合)达到充分利用多种医(CT)、数字减影血管造影(DSA)、单光子发射断疗图像的信息使之互补的目的。医学图像的融合,层成像(SPECT)、磁共振成像(MRI)、数字荧光可有效地丰富图像信息,提高诊断效率。造影(DF)、正电子发射断层成像(PET)等多种成对于图
5、像的融合关键问题在于将多幅图像中像模式的医学图像。这些图像使人体内部的解剖的对应信息进行匹配、对齐,所以对于图像中特征结构、生理功能等多方面的状况通过影像反映出的提取是图像融合要解决的主要问题。本文以一来,从而更加直观地提供了人体解剖、生理及病理幅头部CT为例介绍从CT图片中进行轮廓提取等信息。但是,由于各种成像模式的相关信息内容的过程。对于轮廓的提取可分为图片的预处理、阈各有侧重,如CT、MRI图像主要反映解剖形态结值化处理(细化处理)、轮廓的跟踪提取几步。[1]构;核影像对人体功能状态非常敏
6、感等。更进1图片预处理一步,即使CT、MRI图像反映的解剖信息也是各采用BMP格式510×510的256灰度CT图有侧重,CT对骨骼的显示较为清晰;MRI则对软片(图1)作为研究对象。对于一个物体往往通过组织显示较有优势。总之,单独用某种模式的图像其轮廓就可达到识别的目的,因此本文采用边缘进行诊断往往不能提供医生所需的足够的诊断信作为识别对象。对于图像中反映出来的边缘特征,在灰度图片中实际上就是周围像素灰度值阶跃性[收稿日期]2002-08-23[本文作者]陈欠根副教授。的变化。在CT图像中,骨
7、骼作为图片中的高亮区陈大舵硕士研究生。·10·18卷1期¹2003.1©1994-2008ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.http://www.cnki.net医疗设备信息研究论文和软组织结构之间的灰度值存在明显的跳变,而跟踪算法的跟踪准则为:从图片起始点,定义搜索对于骨骼内部、软组织内部纹理特征灰度值变化方向为左上方(BMP图像的读图顺序从左下角开相对平缓,所以可采用边缘检测算子对外部轮廓始),开始
8、搜索,当检测到相对于起始点的第一个以及内部软组织进行处理。本文选用的是经典的边缘点(灰度值为0)开始按定义的搜索方向进行[3]Roberts算子,对CT图像进行边缘检测,其边下一边缘点的搜索。如果左上方是黑点则定义缘检测算子为为边缘点,如果不是黑点则1222o按所取方向顺时针旋转45g(i,j)={[I(i,j)-I(i+1,j+1)+I(i+1,j)-I(i],j+1)[}直到找到边缘点为止,然后o其中:I(x,y)为像素点在i,j处的灰度值;g再把边缘点在当前搜索的基础上逆时针旋转90,(i