基于均值漂移图像分割技术去除CT图像金属伪影
|
摘要
CT系统在诊断时,金属作为一种高密度的物体在射线通过时的衰减系数比人体其他组织高很多,从而引起X射线被这些物体作用后急剧衰减,导致对应的投影数据失真,重建时图像上会出现伪影。这些由于金属置入物带来的伪影被统称为金属伪影。而传统的方法不能完全准确分割,所以不满足临床要求。应用均值漂移分割图像技术分割金属图像,并应用临床数据进行了验证,结果证明均值漂移分割技术用于CT金属伪影去除分割过程效果比较理想。
In the diagnosis of CT system, as a high density of object, metal has a higher attenuation coefficient than other tissues of human body, which greatly attenuates X-ray, resulting in data distortion of corresponding projection, and artifact on reconstruction image. These artifacts caused by metal implants are referred as metal artifacts. Conventional segmentation methods are not accurate enough and can not be used as clinical analysis. In this paper, a state of the art Mean-Shift technique is employed to improve the accuracy of metallic object segmentation. Real patient datasets are studied to evaluate the efficacy of the method. The result shows that the method can effectively reduce CT metal artifact.
In the diagnosis of CT system, as a high density of object, metal has a higher attenuation coefficient than other tissues of human body, which greatly attenuates X-ray, resulting in data distortion of corresponding projection, and artifact on reconstruction image. These artifacts caused by metal implants are referred as metal artifacts. Conventional segmentation methods are not accurate enough and can not be used as clinical analysis. In this paper, a state of the art Mean-Shift technique is employed to improve the accuracy of metallic object segmentation. Real patient datasets are studied to evaluate the efficacy of the method. The result shows that the method can effectively reduce CT metal artifact.
引文
[1]Gu JW,Zhang L,Chen ZQ,et al.Qinghua Daxue Xuebao:Ziran Kexueban.2006;46(6):825-828.谷建伟,张丽,陈志强,等.基于投影的CT图像金属伪影非线性权重校正[J].清华大学学报:自然科学版,2006,46(6):825-828.
[2]Gu JW,Zhang L,Chen ZQ,et al.Hedianzixue yu Tance Jishu.2006;26(6):905-907.谷建伟,张丽,陈志强,等.常用的几种插值方法在CT系统金属伪影校正中的性能比较[J].核电子学与探测技术,2006,26(6):905-907.
[3]Gonzalez RC,Woods RE,Eddins SL.Digital Image Processing Using MATLAB,Pearson Education.2004.
[4]Bal M,Spies L.Metal artifact reduction in CT using tissue-class modeling and adaptive prefiltering.Med Phys.2006;33(8):2852-2859.
[5]Oehler M,Buzug TM.A Sinogram-Based Metal Artifact Suppression Strategy for Transmission Computed Tomography.Geometriebestimmung mit industrieller Computertomographie,PTB-Bericht PTB-F-54,Bartscher(Ed.),Brauschweig.2007:255-262.
[6]Ziegler A,Nielsen T,Grass M.Iterative reconstruction of a region of interest for transmission tomography.Med Phys.2008;35(4):1317-1327.
[7]Lu QW,Chen WF.Medical image automatically opti-mizes segmentation based on difference of mutual information.Sci China Ser E Inform Sci.2006;36(6):657-667.
[8]Li T,Li X,Wang J,et al.Nonlinear sinogram smoothing for low-dose x-ray CT.IEEE Trans Nuc Sci.2004;51:2505-2513.
[9]Aki K,Lee WH.Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes.Geophys Res.1976;81(23):4381-4399.
[10]Oehler M,Buzug TM.Maximum-Likelihood-Ansatzzur Metall artefact reduction beider Computertomographie.Bildverarbeitung für die Medizin.Mai,Springer,Berlin,2006:36-40.
[11]Oehler M,Buzug TM.Sinogram inpainting for metal artifact reduction in CT images.Medical&Biological Engineering&Computing.2008.
[12]Siddiqui KM,Johnson JP,Reiner BI,et al.Discrete cosine transform JPEG compression vs.2D JPEG2000compression:JNDmetrix visual discrimination model image quality analysis.Proc SPIE.2005;5748:202-207.
[13]Huang Y,Wang WX,Li WS.Anisotropic filter based modified canny algorithm//Fourth International Conference on Fuzzy Systems and Knowledge Discovery.2007;1:736-740.
[14]Lampert CH.Anisotropic Gaussian filtering using fixed point arithmetic//Proc IEEE International Conference in Image Processing.2006:1565-1568.
[15]Buades A,Chien A,Morel JM,et al.Topology preserving linear filtering applied to medical imaging.Recent UCLA Computational and Applied Mathematics Reports.2007.
[16]LüQW,Chen WF.Jisuanji Xuebao.2006;29(2):296-301.吕庆文,陈武凡.基于互信息量的图像分割[J].计算机学报,2006,29(2):296-301.
[17]Wang J,Li T,Lu H,et al.Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose x-ray computed tomography.IEEE Trans Med Imaging.2006;25(10):1272-1283.
[18]Whiting BR.Signal statistics in x-ray computed tomography.Pro SPIE Med Imaging.2002;4682:53-60.
[19]Mahnken AH,Raupach R,Wildberger JE,et al.A new algorithm for metal artifact reduction in computed tomography:in vitro and in vivo evaluation after total hip replacement.Invest Radiol.2003;38:769-775.
[20]Zhao DP.New advance of seismic tomography and its application to subduction zones and earthquake fault zones:A review.Island Arc.2001;10(3):68-84.
[21]Guo B.Bingjing:Zhongguo Dizhenju Dizhi Yanjiusuo Shoushi Lunwen.2003.郭飚.流动宽频带地震台阵:地震定位与地层层析成像[D].北京:中国地震局地质研究所硕士论文,2003.
[22]Lin ZC,Shi QY.Zhongguo Tuxiang Tuxing Xuelao.2001;6(2):142-147.林宙辰,石青云.用四次多项式插值消除医用X射线CT中的金属伪影[J].中国图像图形学报,2001,6(2):142-147.
[23]Nikou C,Galatsanos NP,Likas AC.A class-adaptive spatially variant mixture model for image segmentation.IEEE Trans Image Process.2007;16(4):1121-1130.
[24]Zhao S,Robertson DD,Wang G,et al.X-ray CT metal artifact reduction using wavelets:an application for imaging total hip prostheses.IEEE Trans Med Imaging.2000;19(12):1238-1247.
[25]Ma JH,Yang D,Huang J,et al.Jisuanji Gongcheng yu Yingyong.2008;44(24):168-171.马建华,杨迪,黄静,等.基于NL-PF和MIMS的CT金属伪影消除算法[J].计算机工程与应用,2008,44(24):168-171.
[26]Yang YF,Su ZX,Hu JY.Zhongguo Tuxiang Tuxing Xuebao.2005;10(10):1248-1251.杨云峰,苏志勋,胡金燕.一种保持边缘特征的图像插值方法[J].中国图象图形学报,2005,10(10):1248-1251.
[27]Gu JW,Zhang L,Chen ZQ,et al.CT Lilun yu Yingyong Yanjiu.2005;14(3):24.谷建伟,张丽,陈志强,等.工业CT图像的伪影成因和校正方法综述[J].CT理论与应用研究,2005,14(3):24.
[28]Hou Y,Guo BL.Jisuanji Gongcheng yu Yingyong.2007;43(36):29-30.侯叶,郭宝龙.基于等周算法的图像分割[J].计算机工程与应用,2007,43(36):29-30.
[29]Takahashi Y,Mori S,Kozuka T,et al.Preliminary.Study of correction of original metal artifacts due to I-125seeds in postimplant dosimetry for prostate permanent implant.Radiat Med.2006;24(2):133-138.
[2]Gu JW,Zhang L,Chen ZQ,et al.Hedianzixue yu Tance Jishu.2006;26(6):905-907.谷建伟,张丽,陈志强,等.常用的几种插值方法在CT系统金属伪影校正中的性能比较[J].核电子学与探测技术,2006,26(6):905-907.
[3]Gonzalez RC,Woods RE,Eddins SL.Digital Image Processing Using MATLAB,Pearson Education.2004.
[4]Bal M,Spies L.Metal artifact reduction in CT using tissue-class modeling and adaptive prefiltering.Med Phys.2006;33(8):2852-2859.
[5]Oehler M,Buzug TM.A Sinogram-Based Metal Artifact Suppression Strategy for Transmission Computed Tomography.Geometriebestimmung mit industrieller Computertomographie,PTB-Bericht PTB-F-54,Bartscher(Ed.),Brauschweig.2007:255-262.
[6]Ziegler A,Nielsen T,Grass M.Iterative reconstruction of a region of interest for transmission tomography.Med Phys.2008;35(4):1317-1327.
[7]Lu QW,Chen WF.Medical image automatically opti-mizes segmentation based on difference of mutual information.Sci China Ser E Inform Sci.2006;36(6):657-667.
[8]Li T,Li X,Wang J,et al.Nonlinear sinogram smoothing for low-dose x-ray CT.IEEE Trans Nuc Sci.2004;51:2505-2513.
[9]Aki K,Lee WH.Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes.Geophys Res.1976;81(23):4381-4399.
[10]Oehler M,Buzug TM.Maximum-Likelihood-Ansatzzur Metall artefact reduction beider Computertomographie.Bildverarbeitung für die Medizin.Mai,Springer,Berlin,2006:36-40.
[11]Oehler M,Buzug TM.Sinogram inpainting for metal artifact reduction in CT images.Medical&Biological Engineering&Computing.2008.
[12]Siddiqui KM,Johnson JP,Reiner BI,et al.Discrete cosine transform JPEG compression vs.2D JPEG2000compression:JNDmetrix visual discrimination model image quality analysis.Proc SPIE.2005;5748:202-207.
[13]Huang Y,Wang WX,Li WS.Anisotropic filter based modified canny algorithm//Fourth International Conference on Fuzzy Systems and Knowledge Discovery.2007;1:736-740.
[14]Lampert CH.Anisotropic Gaussian filtering using fixed point arithmetic//Proc IEEE International Conference in Image Processing.2006:1565-1568.
[15]Buades A,Chien A,Morel JM,et al.Topology preserving linear filtering applied to medical imaging.Recent UCLA Computational and Applied Mathematics Reports.2007.
[16]LüQW,Chen WF.Jisuanji Xuebao.2006;29(2):296-301.吕庆文,陈武凡.基于互信息量的图像分割[J].计算机学报,2006,29(2):296-301.
[17]Wang J,Li T,Lu H,et al.Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose x-ray computed tomography.IEEE Trans Med Imaging.2006;25(10):1272-1283.
[18]Whiting BR.Signal statistics in x-ray computed tomography.Pro SPIE Med Imaging.2002;4682:53-60.
[19]Mahnken AH,Raupach R,Wildberger JE,et al.A new algorithm for metal artifact reduction in computed tomography:in vitro and in vivo evaluation after total hip replacement.Invest Radiol.2003;38:769-775.
[20]Zhao DP.New advance of seismic tomography and its application to subduction zones and earthquake fault zones:A review.Island Arc.2001;10(3):68-84.
[21]Guo B.Bingjing:Zhongguo Dizhenju Dizhi Yanjiusuo Shoushi Lunwen.2003.郭飚.流动宽频带地震台阵:地震定位与地层层析成像[D].北京:中国地震局地质研究所硕士论文,2003.
[22]Lin ZC,Shi QY.Zhongguo Tuxiang Tuxing Xuelao.2001;6(2):142-147.林宙辰,石青云.用四次多项式插值消除医用X射线CT中的金属伪影[J].中国图像图形学报,2001,6(2):142-147.
[23]Nikou C,Galatsanos NP,Likas AC.A class-adaptive spatially variant mixture model for image segmentation.IEEE Trans Image Process.2007;16(4):1121-1130.
[24]Zhao S,Robertson DD,Wang G,et al.X-ray CT metal artifact reduction using wavelets:an application for imaging total hip prostheses.IEEE Trans Med Imaging.2000;19(12):1238-1247.
[25]Ma JH,Yang D,Huang J,et al.Jisuanji Gongcheng yu Yingyong.2008;44(24):168-171.马建华,杨迪,黄静,等.基于NL-PF和MIMS的CT金属伪影消除算法[J].计算机工程与应用,2008,44(24):168-171.
[26]Yang YF,Su ZX,Hu JY.Zhongguo Tuxiang Tuxing Xuebao.2005;10(10):1248-1251.杨云峰,苏志勋,胡金燕.一种保持边缘特征的图像插值方法[J].中国图象图形学报,2005,10(10):1248-1251.
[27]Gu JW,Zhang L,Chen ZQ,et al.CT Lilun yu Yingyong Yanjiu.2005;14(3):24.谷建伟,张丽,陈志强,等.工业CT图像的伪影成因和校正方法综述[J].CT理论与应用研究,2005,14(3):24.
[28]Hou Y,Guo BL.Jisuanji Gongcheng yu Yingyong.2007;43(36):29-30.侯叶,郭宝龙.基于等周算法的图像分割[J].计算机工程与应用,2007,43(36):29-30.
[29]Takahashi Y,Mori S,Kozuka T,et al.Preliminary.Study of correction of original metal artifacts due to I-125seeds in postimplant dosimetry for prostate permanent implant.Radiat Med.2006;24(2):133-138.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心