发明名称 医学超声图像自适应SUSAN扩散去噪方法
摘要 本发明涉及一种医学超声图像自适应SUSAN扩散去噪方法,按以下步骤进行:(1)自适应选取图像全斑点噪声区域,(2)基于SUSAN算法的图像边缘检测,(3)基于SUSAN扩散的医学超声图像去噪。与现有技术相比具有减轻了人工设定阈值参数的工作量,提高了算法的鲁棒性,更加适合推广应用等特点。
申请公布号 CN102663692A 申请公布日期 2012.09.12
申请号 CN201210085740.0 申请日期 2012.03.28
申请人 汕头大学 发明人 沈民奋;方若宇;孙丽莎;张琼;徐宇贵
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 广州三环专利代理有限公司 44202 代理人 温旭
主权项 1.一种医学超声图像自适应SUSAN扩散去噪方法,按以下步骤进行:(1)自适应选取图像全斑点噪声区域:选取一张图像或图像中的一部分区域,通过图像预处理、特征值分解和K均值法对图像进行自适应全斑点噪声区域的选取:1a  图像预处理:将带噪图像与高斯核函数进行卷积,得到图像的结构张量矩阵,其采用的计算方法如下:<img file="198649DEST_PATH_IMAGE001.GIF" wi="324" he="91" />其中,<img file="335232DEST_PATH_IMAGE002.GIF" wi="57" he="29" />为高斯核函数,<img file="941794DEST_PATH_IMAGE003.GIF" wi="70" he="37" />为图像分别在横轴和纵轴上的梯度;1b  特征值分解:对结构张量矩阵进行特征值分解,得到分解后的矩阵形式为:<img file="935158DEST_PATH_IMAGE004.GIF" wi="294" he="86" />其中,<img file="630100DEST_PATH_IMAGE005.GIF" wi="72" he="26" />为结构张量的特征向量,<img file="304795DEST_PATH_IMAGE006.GIF" wi="23" he="27" />代表的是波动最大的方向,即梯度方向,而<img file="828180DEST_PATH_IMAGE007.GIF" wi="29" he="23" />则指向波动最小的方向,<img file="195708DEST_PATH_IMAGE008.GIF" wi="68" he="31" />则是对应的特征向量方向的扩散强度;1c  基于K均值的自适应全斑点噪声选取:选取模板并遍历图像,求出整个模板中的<img file="375016DEST_PATH_IMAGE009.GIF" wi="65" he="28" />之和<img file="587823DEST_PATH_IMAGE010.GIF" wi="35" he="34" />,通过采用K均值算法找出图像的全斑点噪声区域,计算公式如下:<img file="168977DEST_PATH_IMAGE011.GIF" wi="334" he="34" />其中,<img file="441826DEST_PATH_IMAGE012.GIF" wi="237" he="57" />为模板下的总像素点;(2)基于SUSAN算法的图像边缘检测:找出图像中的全斑点噪声区域后,可计算得到<img file="108431DEST_PATH_IMAGE013.GIF" wi="14" he="26" />值,计算公式如下:<img file="121999DEST_PATH_IMAGE014.GIF" wi="472" he="42" />其中,<img file="619976DEST_PATH_IMAGE015.GIF" wi="187" he="31" />分别代表图像分别在四个邻域方向上的灰度值差分,Median代表取所有数据中的中值;然后,采用SUSAN算法可对图像边缘进行精确检测,方法如下:采用SUSAN圆形模板中心与图像中各像素重合,利用如下判别公式对图像中的各点是否属于USAN区域进行判别<img file="329306DEST_PATH_IMAGE016.GIF" wi="207" he="61" />其中,<img file="217628DEST_PATH_IMAGE017.GIF" wi="41" he="33" />为模板中心像素灰度值,<img file="100133DEST_PATH_IMAGE018.GIF" wi="37" he="33" />为模板中其他像素灰度值,继而,可计算出图像中各像素点的USAN大小:<img file="390300DEST_PATH_IMAGE019.GIF" wi="179" he="76" />采用SUSAN算法进行图像边缘检测就是比较<img file="270532DEST_PATH_IMAGE020.GIF" wi="45" he="36" />与几何阈值的大小,如果<img file="646149DEST_PATH_IMAGE020.GIF" wi="45" he="36" />小于<img file="66767DEST_PATH_IMAGE021.GIF" wi="20" he="32" />,则认为此时模板中心像素点为图像边缘点,反之,则认为该中心像素不是图像的边缘点,其边缘响应计算式如下:<img file="477019DEST_PATH_IMAGE022.GIF" wi="269" he="85" />其中,<img file="528152DEST_PATH_IMAGE023.GIF" wi="150" he="64" />为模板中包含的像素点总数;(3)基于SUSAN扩散的医学超声图像去噪:求得以上的边缘响应后,然后通过以下公式进行非线性扩散去噪:<img file="405714DEST_PATH_IMAGE024.GIF" wi="464" he="184" />其中,常数<img file="567705DEST_PATH_IMAGE025.GIF" wi="248" he="30" />为(2)中的边缘响应,<img file="566885DEST_PATH_IMAGE026.GIF" wi="58" he="26" />代表中值绝对偏差运算,<img file="788919DEST_PATH_IMAGE027.GIF" wi="23" he="31" />为原图像,通过<img file="201446DEST_PATH_IMAGE028.GIF" wi="27" he="24" />次非线性迭代进行去噪,最后便可得到去噪后的图像。
地址 515063 广东省汕头市大学路243号