CN104318521A - 多线性子空间学习的医疗图像去噪方法 - Google Patents
多线性子空间学习的医疗图像去噪方法 Download PDFInfo
- Publication number
- CN104318521A CN104318521A CN201410515240.5A CN201410515240A CN104318521A CN 104318521 A CN104318521 A CN 104318521A CN 201410515240 A CN201410515240 A CN 201410515240A CN 104318521 A CN104318521 A CN 104318521A
- Authority
- CN
- China
- Prior art keywords
- denoising
- image
- space
- algorithm
- fritter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明为一种多线性子空间学习的医疗图像去噪方法,将时域和频域信息相结合,并采用多线性子空间学习方法,能够对图像中提取的小块直接去噪。步骤一、空间转换模块将超声图像通过对数变换,变换到对数域后进行后续计算;步骤二、分块模块将超声图像分成相互重叠的相等的小块;步骤三、聚类模块将分好的小块通过聚类方法分成N类,用于后续去噪计算处理;步骤四、去噪模块将每一类的小块直接排成体数据,应用多线性子空间学习方法对每个小块进行去噪;步骤五、聚集模块将去噪后的每一个小块重新返回到图像上的原始位置;步骤六、空间逆转换模块将第五步获得的图像,通过指数变换,变换回最原始的图像空间,所得到的图像即为最终的去噪图像。
Description
技术领域
本发明为一种多线性子空间学习的医疗图像去噪方法,可用于多种医学图像的噪声去除。
背景技术
对器官的精确显示在现代医疗中起到非常重要的作用。在很多医学成像技术中,由于超声成像的非侵入性,硬件的便携性,实时成像以及性价比高等特性已被广泛的应用于医疗辅助诊断中。但相比于其他成像设备,如CT和核磁共振,超声的噪声非常多。通常认为超声的噪声是乘性斑点噪声,其降低了图像质量,模糊了细节信息以及物体边缘,限制了图像的对比度。因此,如何减少斑点噪声并保持图像信息是目前研究的热点。
目前现有的图像去噪方法可根据噪声产生的时间大体分为两类:混合方法和后处理方法。混合方法需要在同一区域产生多幅图像,并把它们结合成一幅图像。此方法并不是本发明所覆盖的范围。本发明的内容是关于后处理的方法研究,也就在超声图像生成后进行图像去噪。到目前为止,已经有很多图像去噪方法被提出来,可分成三类:频域去噪方法,时域去噪方法和视频域去噪方法。目前研究表明,主成分分析法可有效去除图像噪声。首先将图像分成很多相等的小块。由于主成分分析法是线性子空间学习方法,只能对矢量进行处理。因此,首先要将这些小块展开成矢量后,再用主成分分析法进行处理。这一过程,会破坏每个小块中本来的结构特征,另外由于矢量的维数也许大于样本的个数,会造成过度拟合的现象。
发明内容
为了克服以上问题,本发明将时域和频域信息相结合,并采用多线性子空间学习方法,能够对图像中提取的小块直接去噪。
多线性子空间学习的医疗图像去噪方法,包括以下步骤:
步骤一、空间转换模块将超声图像通过对数变换,变换到对数域后进行后续计算;
步骤二、分块模块将超声图像分成相互重叠的相等的小块;
步骤三、聚类模块将分好的小块通过聚类方法分成N类,用于后续去噪计算处理;
步骤四、去噪模块将每一类的小块直接排成体数据,应用多线性子空间学习方法对每个小块进行去噪;
步骤五、聚集模块将去噪后的每一个小块重新返回到图像上的原始位置;
步骤六、空间逆转换模块将第五步获得的图像,通过指数变换,变换回最原始的图像空间,所得到的图像即为最终的去噪图像。
本发明的有益效果:本发明将时域和频域信息相结合,并采用多线性子空间学习方法,对图像中提取的小块直接去噪。由于图像被分成小块,并对相似的小块进行去噪,此过程可用GPU进行加速,以提高算法的计算速度。
附图说明
图1为图像去噪方法流程图;
图2为基于广义多维独立成分分析法的空间转换;
图3为三维体数据分别在第一个模式,第二个模式和第三个模式展开得到的三个矩阵。
具体实施方式
下面结合具体实施例和附图详细说明本发明。
如图1所示,本发明的基于多线性子空间学习方法的快速医疗图像去噪方法包括以下几个步骤:
步骤S101,读取待去噪的图像数据。
步骤S102,将读取的图像数据,进行对数变换,变换到对数空间,将图像空间中乘性噪声转换成加性噪声,使得后续去噪处理更加有效。式1为原始带有乘性噪声的图像数据gor(x,y):
将式1两边同时取对数,得到式2
步骤S103,将对数空间的图像数据,分成大小相等的小块。小块的大小可取为5×5,10×10等,也可根据图像大小及噪声大小自行决定。每个小块之间可以相互重叠,也可以不重叠但相交。
步骤S104,将图像数据中的小块通过聚类方法,或者最近距离方法分成N类。以K-means的聚类方法为例。若每个小块为xq,q=1,2,...,m,m为图像数据中所有小块的个数,通过找到聚类中心C=[c1,c2,...,cn,...,cN],n=1,2,...,N的均值μn使得式3最小
步骤S105,被分为同一类的小块组成一个三阶张量附图2所示为通过多线性子空间学习方法对此三阶张量进行去噪。附图3示出了将三维体数据分别在第一个模式,第二个模式和第三个模式展开得到的三个矩阵。通过求解能量方程式4,广义多维独立成分分析可被用来得到三阶张量每个模式上的变换矩阵Ui(z),z=1,2,3。
以往的研究中,为了去除高频噪声,用户通过经验,保留对应于特征值最大的前J个特征向量组成的Ui(z),用于小块去噪。但由于保留个数J是预先给定的,没有自适应特点,本发明采取一种低秩最小二乘最小均方误差(LMMSE)法,自适应的对每组相似小块进行高频噪声的去除。也就是在广义多维独立成分分析的过程中,对独立成分不降维。而去噪降维的过程由低秩最小二乘最小均方误差法来实现。通过广义多维独立成分分析法得到的核张量为式5
其中是Ui(z)的转置。重建后的小块即为去噪小块。重建方程如式6所示
其中, 噪声方差 Median(HH)是高-高带小波系数;×1,×2,×3分别为将体数据γi在第一个模式,第二个模式和第三个模式展开得到的矩阵上进行乘法运算。
在此步骤运行过程中,为了提高计算效率,采用GPU进行加速。通过对不同类别组小块的并行去噪处理,减少计算时间。
步骤S106,去噪后的小块返回到其原始位置。
步骤S107,将去噪后的的图像数据,进行指数变换,变换到原图像空间。
步骤S108,得到最终去噪后的结果。
虽然参考优选实施例对本发明进行描述,但以上所述实例并不构成本发明保护范围的限定,任何在本发明的精神及原则内的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围内。
Claims (5)
1.多线性子空间学习的医疗图像去噪方法,其特征在于,包括以下步骤:
步骤一、空间转换模块将读取的超声图像进行对数变换,变换到对数空间,从而将图像空间中的乘性噪声转换成对数空间中的加性噪声;
步骤二、分块模块将超声图像分成相互重叠的相等的小块;
步骤三、聚类模块将分好的小块通过聚类方法分成N类,用于后续去噪计算处理;
步骤四、去噪模块将每一类的小块直接排成体数据,应用多线性子空间学习方法对每个小块进行去噪;
步骤五、聚集模块将去噪后的每一个小块重新返回到图像上的原始位置;
步骤六、空间逆转换模块将第五步获得的图像,通过指数变换,变换回最原始的图像空间,所得到的图像即为最终的去噪图像。
2.如权利要求1所述的多线性子空间学习的医疗图像去噪方法,其特征在于,其中步骤二中的小块的大小根据实际情况不同大小取不同,为正方型或长方形,且小块之间为重叠或相交。
3.如权利要求1或2所述的多线性子空间学习的医疗图像去噪方法,其特征在于,其中步骤三中采用以下聚类方法中的一种:K-means算法,K-medoids算法,Clarans算法,BIRCH算法、CURE算法、CHAMELEON算法,DBSCAN算法、OPTICS算法、DENCLUE算法,图论聚类法,STING算法、CLIQUE算法、WAVE-CLUSTER算法。
4.如权利要求1或2所述的多线性子空间学习的医疗图像去噪方法,其特征在于,其中步骤四中通过低秩最小二乘最小均方误差法实现去噪降维。
5.如权利要求4所述的多线性子空间学习的医疗图像去噪方法,其特征在于,其中对步骤四采用GPU进行加速,对不同类别组小块进行并行去噪处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410515240.5A CN104318521A (zh) | 2014-09-29 | 2014-09-29 | 多线性子空间学习的医疗图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410515240.5A CN104318521A (zh) | 2014-09-29 | 2014-09-29 | 多线性子空间学习的医疗图像去噪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104318521A true CN104318521A (zh) | 2015-01-28 |
Family
ID=52373747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410515240.5A Pending CN104318521A (zh) | 2014-09-29 | 2014-09-29 | 多线性子空间学习的医疗图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104318521A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897820A (zh) * | 2018-06-21 | 2018-11-27 | 电子科技大学 | 一种denclue算法的并行化方法 |
CN113313717A (zh) * | 2021-05-18 | 2021-08-27 | 哈尔滨理工大学 | 一种基于深度学习的分节段椎骨ct图像分割方法 |
CN113436093A (zh) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | 一种基于能量域变换的遥感影像去噪方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102664010A (zh) * | 2012-05-04 | 2012-09-12 | 山东大学 | 一种基于多因素频率位移不变特征的鲁棒说话人辨别方法 |
CN103218617A (zh) * | 2013-05-13 | 2013-07-24 | 山东大学 | 一种多线性大间距的特征提取方法 |
-
2014
- 2014-09-29 CN CN201410515240.5A patent/CN104318521A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102664010A (zh) * | 2012-05-04 | 2012-09-12 | 山东大学 | 一种基于多因素频率位移不变特征的鲁棒说话人辨别方法 |
CN103218617A (zh) * | 2013-05-13 | 2013-07-24 | 山东大学 | 一种多线性大间距的特征提取方法 |
Non-Patent Citations (4)
Title |
---|
ANIL K.JAIN: ""Data clustering:50 years beyond K-means"", 《PATTERN RECOGNITION LETTERS》 * |
DANNI AI 等: ""Multiresolution generalized N dimension PCA for ultrasound image denoising"", 《BIOMEDICAL ENGINEERING ONLINE》 * |
TAO LIN 等: ""Survey of hyperspectral image denoising methods based on tensor decompositions"", 《EURASIP JOURNAL ON ADVANCES IN SIGNAL PROCESSING》 * |
XIAN GUO 等: ""Hyperspectral image noise reduction based on rank-1 tensor decomposition"", 《ISPRS JOURNAL OF PHOTOGRAMMETRY AND REMOTE SENSING》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108897820A (zh) * | 2018-06-21 | 2018-11-27 | 电子科技大学 | 一种denclue算法的并行化方法 |
CN108897820B (zh) * | 2018-06-21 | 2022-03-15 | 电子科技大学 | 一种denclue算法的并行化方法 |
CN113313717A (zh) * | 2021-05-18 | 2021-08-27 | 哈尔滨理工大学 | 一种基于深度学习的分节段椎骨ct图像分割方法 |
CN113436093A (zh) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | 一种基于能量域变换的遥感影像去噪方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Choi et al. | Speckle noise reduction technique for SAR images using statistical characteristics of speckle noise and discrete wavelet transform | |
Dey et al. | Parameter optimization for local polynomial approximation based intersection confidence interval filter using genetic algorithm: an application for brain MRI image de-noising | |
Fan et al. | Interference mitigation for synthetic aperture radar based on deep residual network | |
Cozzolino et al. | Nonlocal CNN SAR image despeckling | |
CN111091603B (zh) | 一种超声成像方法、装置、可读存储介质及终端设备 | |
Aggarwal et al. | Line detection in images through regularized Hough transform | |
US7272265B2 (en) | Methods for performing DAF data filtering and padding | |
Dalsasso et al. | SAR image despeckling by deep neural networks: From a pre-trained model to an end-to-end training strategy | |
Zhang et al. | Two-direction nonlocal model for image denoising | |
Yu et al. | PCANet based nonlocal means method for speckle noise removal in ultrasound images | |
CN114677300A (zh) | 一种基于双阶段学习框架的高光谱图像深度降噪的方法及*** | |
CN103955904A (zh) | 一种基于离散分数阶傅里叶变换相位信息的信号重建方法 | |
Huang et al. | Speckle noise reduction in sonar image based on adaptive redundant dictionary | |
Zhou et al. | Deep multi-scale recurrent network for synthetic aperture radar images despeckling | |
Kong et al. | A new 4-D nonlocal transform-domain filter for 3-D magnetic resonance images denoising | |
Liu et al. | Contourlet-CNN for SAR image despeckling | |
Bhardwaj et al. | A Novel Method for Despeckling of Ultrasound Images Using Cellular Automata-Based Despeckling Filter | |
CN104318521A (zh) | 多线性子空间学习的医疗图像去噪方法 | |
Isar et al. | Hyperanalytic wavelet-based robust edge detection | |
CN104978716A (zh) | 一种基于线性最小均方误差估计的sar图像降噪方法 | |
Singh et al. | A homomorphic non-subsampled contourlet transform based ultrasound image despeckling by novel thresholding function and self-organizing map | |
CN104217430A (zh) | 基于l1正则化的图像显著性检测方法 | |
Li et al. | MCSNet: A radio frequency interference suppression network for spaceborne SAR images via multi-dimensional feature transform | |
Zhou et al. | High-resolution hierarchical adversarial learning for OCT speckle noise reduction | |
Kaur et al. | Super resolution of magnetic resonance images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150128 |
|
WD01 | Invention patent application deemed withdrawn after publication |