CN106056565B - 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 - Google Patents
一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 Download PDFInfo
- Publication number
- CN106056565B CN106056565B CN201610374485.XA CN201610374485A CN106056565B CN 106056565 B CN106056565 B CN 106056565B CN 201610374485 A CN201610374485 A CN 201610374485A CN 106056565 B CN106056565 B CN 106056565B
- Authority
- CN
- China
- Prior art keywords
- image
- mri
- pet
- bright
- dark
- 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.)
- Active
Links
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 33
- 230000000877 morphologic effect Effects 0.000 title claims abstract description 31
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 15
- 238000001914 filtration Methods 0.000 title claims abstract description 8
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 47
- 230000006835 compression Effects 0.000 claims abstract description 29
- 238000007906 compression Methods 0.000 claims abstract description 29
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 230000010339 dilation Effects 0.000 claims description 2
- 230000003628 erosive effect Effects 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract description 15
- 238000000844 transformation Methods 0.000 abstract 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 54
- 230000003595 spectral effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 241000405965 Scomberomorus brasiliensis Species 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000002503 metabolic effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000016776 visual perception Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20028—Bilateral filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明一种基于多尺度形态学双边滤波和对比度压缩的MRI与PET图像融合方法:对PET彩色图像插值运算,使分辨率与MRI图像相同,对R、G、B三通道IHS变换,得I、H、S三通道;对PET图像的I通道和MRI图像多尺度形态学双边滤波分解,得亮特征图像、暗特征图像和基频图像;对MRI图像分解得到的亮特征和暗特征图像对比度压缩、调整灰度范围,使之与I通道分解得到亮特征和暗特征图像灰度范围相同;将对比度压缩后MRI亮特征和暗特征图像与I通道分解得到的亮特征和暗特征图像进行最大值运算,得矫正后亮特征和暗特征图像,矫正后亮特征图像减去矫正后暗特征图像,加上PET图像分解得到的基频图像得到融合亮度图FI;将FI、H和S三通道进行IHS逆变换,得到融合后彩色图像。
Description
(一)技术领域
本发明涉及一种基于多尺度形态学双边滤波分解和对比度压缩的MRI与PET图像融合方法,属于数字图像处理领域,主要涉及数学形态学和图像融合技术。MRI和PET融合得到的图像,既包含MRI图像提供的组织结构解剖信息,又包含PET图像提供的组织结构代谢强度信息,因此更有利于病情的诊断和病灶的定位,具有广阔的运用前景。
(二)背景技术
MRI图像(核磁共振成像)和PET图像(正电子发射型计算机断层成像)融合,是将高分辨率的具有组织解剖信息而不具有功能信息的MRI图像,与低分辨率的具有组织功能信息而不具有组织解剖信息的PET图像融合成一张图像,使该图像同时具有高分辨率的组织解剖信息和功能信息。由于PET是彩色图像,将MRI结构信息注入PET图像时会引起颜色畸变,因此如何在保持PET图像颜色高保真度的同时,将尽可能多的MRI结构信息融合到PET图像中,是MRI与PET图像融合的难点和关键。研究者提出了很多方法进行MRI与PET图像的融合。基于IHS和PCA的融合方法(参见文献:何长涛刘全喜等基于HIS和PCA的多模医学图像融合,能源工程,2010,7:280-285.(He C,Liu Q,Li H,et al.Multimodal medical imagefusion based on IHS and PCA[J].Procedia Engineering,2010,7:280-285.))先采用IHS变换求取PET的亮度图像I,然后用PCA方法将其和MRI分解成两个主成分分量,根据MRI图像和亮度图像I的空域频率融合两主成分分量得到融合后的亮度图像,然后经过IHS逆变换得到融合后的图像。由于该方法采用空域频率全局指标来融合图像,没有对光谱畸变进行处理,因此融合后的图像虽然具有较高的结构信息,但图像整体亮度改变较大,颜色畸变严重。基于小波变换的融合方法(参见文献:哈利巴布等采用小波变换的MRI-PET多模医学图像融合,通讯及其运用,(Haribabu,Maruturi,CH Hima Bindu,and K.Satya Prasad."Multimodal Medical Image Fusion of MRI-PET Using Wavelet Transform."Advancesin Mobile Network,Communication and its Applications (MNCAPPS),2012International Conference on.IEEE,2012.))首先对PET图像进行IHS变换得到其亮度图像,然后分别对亮度图像和MRI图像进行一层离散小波分解,基频采用均值融合的方法,其他高频采用基于空域频率权重的融合方法。由于PET的亮度图像和MRI之间的相关性较低,采用低频均值融合的方法会导致融合图像亮度与PET的亮度图像相差较大,产生较大的颜色畸变。基于混合智能的融合方法(参见文献:卡维莎切拉木图基于混合智能的医学图像融合,运用软计算,2014,20:83-94.(Kavitha C T,Chellamuthu C.Medical imagefusion based on hybrid intelligence[J].Applied Soft Computing,2014,20:83-94.))先通过蚁群算法检测待融合图像的边缘,然后将边缘增强后的图像输入到脉冲耦合神经网络得到对应图像的融合权重图,然后在像素级上选择融合权重较大的输入图像作为融合图像。由于该方法根据MRI像素边沿信息来替换PET图像中对应位置的像素值,因此融合图像中会产生较大的颜色畸变,同时融合图像中不能有效的包含MRI中的区域信息。数学形态学是图像处理的重要工具之一(参见文献:塞拉.图像分析与数学形态学,科学出版社,纽约,1982.(J.Serra,Image analysis and mathematical morphology,Academic Press,New York,1982.)),由MRI图像的成像原理可知,MRI结构信息主要体现为图像中的亮特征和暗特征,在数学形态学的各种运算符中,高帽变换可以用于提取图像中的亮区域和暗区域。不过,由于数学形态学运算可能会严重平滑图像的细节,因此会导致提取的亮特征和暗特征发生畸变。
形态学双边滤波具有一些优良的特性(参见文献:安古洛形态学双边滤波,应用数学学会图像科学杂志,2013,6(3):1790-1822.(Angulo J.Morphological bilateralfiltering[J].SIAM Journal on Imaging Sciences,2013,6(3):1790-1822.))能够在滤除噪声的同时保留图像中大的边缘,因此采用多尺度形态学双边滤波能有效的提取图像中不同尺度的亮特征和暗特征。
(三)发明内容
1、目的:MRI与PET融合得到的图像,既包含MRI图像提供的组织结构解剖信息,又包含PET图像提供的组织结构代谢强度信息,因此更有利于病情的诊断和病灶的定位,具有广阔的运用前景。但是现阶段的融合方法很难做到在保持高光谱保真度的情况下,获得高分辨率的融合图像。
为了解决上述问题并弥补传统方法的不足,本发明提供了一种基于多尺度形态学双边滤波分解和对比度压缩的MRI与PET图像融合方法,通过多尺度形态学双边滤波分解提取MRI和PET亮度图像的亮特征、暗特征和基频图像,根据PET亮特征和暗特征图像压缩MRI的亮特征和暗特征图像的对比度和灰度范围,并与PET图像分解得到的基频图像相加得到融合亮度图像,这样能够在尽可能少的改变PET亮度图像灰度的同时,将MRI图像中尽可能多的视觉结构信息注入到融合图像中,从而获得高分辨率高光谱保真度的融合图像。
2、技术方案:本发明的主要原理为:首先对PET彩色图像进行插值运算,使其分辨率与MRI图像相同,并对其R、G、B三通道进行IHS变换,得到I、H、S三通道;然后分别对PET图像的I通道和MRI图像进行多尺度形态学双边滤波分解,得到对应的亮特征图像、暗特征图像和基频图像;随后分别对MRI图像分解得到的亮特征图像和暗特征图像进行对比度压缩、并调整其灰度范围,使之分别与I通道分解得到的亮特征和暗特征图像的灰度范围相同;其次将对比度压缩后的MRI亮特征和暗特征图像分别与I通道分解得到的亮特征和暗特征图像进行最大值运算,得到矫正后的亮特征图像和暗特征图像,并用矫正后的亮特征图像减去矫正后的暗特征图像,并加上PET分解得到的基频图像得到融合亮度图FI;最后将FI、H和S三通道进行IHS逆变换,得到融合后的彩色图像。
本发明的一种基于多尺度形态学双边滤波分解和对比度压缩的MRI与PET图像融合方法,其步骤为:
步骤一:采用双立方插值的方法对PET图像的R、G、B三通道分别进行插值,使其分辨率和MRI图像的相同,然后对其进行IHS变换、得到I、H、S三通道,其中I表示图像的亮度、H代表图像光谱的波长、S代表光谱的纯度。
步骤二:分别对MRI图像和PET的亮度图像I进行多尺度形态学双边滤波分解。设分解的层数为n,待分解的图像为f0则第i层分解为:
Δi1=max(Wi-1*G(di-1)),i>1的情形;
Δi2=max(Bi-1*G(di-1)),i>1的情形;
di<di-1
Wi=fi-fi-1
Bi=fi-1-fi
其中,
Wi、Bi、fi分别为第i层分解得到的亮特征图像、暗特征图像和基频图像,Wi-1、Bi-1、fi-1为第i-1层分解得到的亮特征图像、暗特征图像和基频图像;Bd,Δ是自适应各向同性扁平结构元素,定义如下
其中d为结构尺度因子,Δ为灰度尺度因子;Δi1、Δi2为第i次分解的灰度尺度因子,di-1、di分别为第i-1和i次分解的结构尺度因子;G(di-1)标准差为di-1的高斯函数,*为卷积运算符,定义如下
其中(x,y)、(m,n)分别为卷积前后图像的像素坐标,M、N为被卷积图像的行列像素数;max()和min()分别为取最大值和最小值函数;和分别为数学形态学开运算、闭运算、膨胀和腐蚀运算符,定义如下
其中(x,y)是图像像素坐标,(u,v)为结构元素中的像素坐标。
对于第一层分解,当i=1时,我们取
W1=f1-f0
B1=f0-f1
其中,M、N为图像f0的像素尺寸。
按照上述方法对输入图像进行n层分解,得到n组基频图像、亮特征图像和暗特征图像,分别为f1、f2……fn,W1、W2……Wn,B1、B2……Bn。设形态学双边滤波分解得到的亮特征、暗特征和基频图像分别为W、B和Base,则其取值如下
Base=fn
令f0为MRI图像,对其进行多尺度形态学双边滤波分解,得到对应的亮特征图像MW、暗特征图像MB;然后令f0为PET的亮度图像I,对其进行多尺度形态学双边滤波分解,得到对应的亮特征图像PW、暗特征图像PB和基频图像PBase。
步骤三:对提取的MRI亮特征图像和暗特征图像分别进行对比度压缩,设待压缩的图像为f,则采用下面的压缩函数对f进行压缩
其中,M、N分别为图像f的行数和列数;式中T和R是输入参数,分别是图像对比度的压缩率和压缩曲线在对数坐标轴下的拐点斜率;参数S、r和p通过联立下列各式求得
其中,
A2=((fmaxfmin)r+Lwa2r)(C-1)
A3=((fmaxfmin)r-Lwa2r)(C-1)
fmax=max(f)
fmin=min(f)
令f分别为MRI图像的亮特征MW和暗特征MB,得到对比度压缩后的亮特征CMW和暗特征CMB。
步骤四:分别用PET分解得到的亮特征图像PW和暗特征图像PB与MRI对比度压缩后的亮特征图像CMW和暗特征图像CMB进行最大值运算,得到新的亮特征图像NCMW和暗特征图像NCMB
其中,(x,y)为图像的像素坐标,如果将灰度图像看作二维标量函数,则Ο1是函数CMW和PW支撑集的交集,则Ο2是函数CMB和PB支撑集的交集。
步骤五:按下式得到融合亮度图像FI
FI=Pbase+NCMW-NCMB
其中Pbase为PET的亮度图像I经过多尺度形态学双边滤波分解后得到的基频图像。
用FI代替PET图像经过IHS变换得到的亮度图像I,并与色调分量H和饱和度分量S一起通过IHS逆变换得到融合后的彩色图像F,从而达到融合高分辨率灰度MRI图像和低分辨率彩色PET图像,得到高分辨率彩色图像的目的。
3、优点及功效:通过多尺度形态学双边滤波分解,能够有效的获得MRI和PET图像的亮特征、暗特征和基频图像,由于MRI基频图像中只包含大尺度的结构,因此只需要将MRI的亮特征和暗特征融合到PET中。通过压缩MRI的亮特征图像和暗特征图像的对比度能够在有效保留MRI视觉结构信息的同时,压缩其动态范围,减少对PET光谱信息的影响。同时MRI亮特征图像和暗特征图像与PET亮特征图像和暗特征图像之间的最大值运算,能够保证得到的亮特征图像和暗特征图像同时包含PET和MRI的结构信息,又避免产生冗余信息,能够进一步减少对PET光谱信息的影响。因此,本发明能够在保证融合图像光谱保真度同时,向融合图像中注入更多的结构信息,获得高光谱保真度高分辨率的融合图像。该发明可广泛应用于MRI和PET图像的融合,具有广阔的市场前景与应用价值。
(四)附图说明
图1为本发明一种基于多尺度形态学双边滤波分解和对比度压缩的MRI与PET图像融合方法的总体原理框图。
图2为多尺度形态学双边滤波分解原理框图。
(五)具体实施方式
为了更好地理解本发明的技术方案,以下结合附图对本发明的实施方式作进一步描述。
本发明的原理框图如图1所示,具体实施步骤如下:
步骤一:对PET图像进行预处理,并将其从RGB颜色空间转换到IHS颜色空间。
首先采用双立方插值的方法对PET图像的R、G、B三通道分别进行插值,使其分辨率和MRI图像的,然后对其进行IHS变换、得到I、H、S三通道,其中I表示图像的亮度、H代表图像光谱的波长、S代表光谱的纯度。
步骤二:分别对MRI图像和PET图像的亮度图像I进行多尺度形态学双边滤波分解,其原理框图如图2所示。
设分解的层数为n,输入的图像的为f0,则多尺度形态学双边滤波分解就是,按照结构尺度d逐渐增大、灰度尺度Δ逐渐减小的次序对输入图像进行形态学双边开运算和闭运算,其中结构尺度和灰度尺度体现在形态学结构元素的大小和形状上。经过形态学双边开运算和闭运算之后,原始图像中的结构尺度和灰度尺度小于给定值的亮特征和暗特征会被消去,因此可以通过原始图像和滤波后的图像之间的差值作为该尺度下的图像亮特征和暗特征。具体方法如下:
d1=3
W1=f1-f0
B1=f0-f1
…
Δi1=max(Wi-1*G(di-1))
Δi2=max(Bi-1*G(di-1))
di=di-1+2
Wi=fi-fi-1
Bi=fi-1-fi
…
Δn1=max(Wn-1*G(dn-1))
Δn2=max(Bn-1*G(dn-1))
dn=dn-1+2
Wn=fn-fn-1
Bn=fn-1-fn
Base=fn
其中,
f0是尺寸为M、N的输入图像;Bd,Δ是自适应各向同性扁平结构元素,其中d为结构元素的最大直径,Δ为最大绝对灰度差;(x,y)是图像像素坐标,(u,v)为结构元素中的像素坐标;G(di-1)是标准差为di-1的高斯函数,*为卷积运算符,定义如下
其中(m,n)为卷积后图像的像素坐标;Wi、Bi、fi分别为第i层分解得到的亮特征图像、暗特征图像和基频图像,W、B、Base分别为经过多尺度双边形态学分解后得到的亮特征图像、暗特征图像和基频图像;和分别是数学形态学开、闭、膨胀和腐蚀运算符;
令f0为MRI图像,对其进行多尺度形态学双边分解,分解的层数n设为6,则得到对应的亮特征图像MW、暗特征图像MB和基频图像;然后令f0为PET亮度图像I,对其进行多尺度双边形态学分解,分解的层数n设为4,则得到对应的亮特征图像PW、暗特征图像PB和基频图像PBase。
步骤三:对提取的MRI亮特征图像和暗特征图像分别进行对比度压缩。
由于通过多尺度形态学双边滤波分解得到亮特征和暗特征图像的灰度范围较大,直接注入到PET图像中会导致严重的颜色畸变,因此需要压缩亮特征图像和暗特征图像的灰度范围。传统的线性灰度变换,并没有考虑人眼视觉感知的特点,灰度经过压缩后,会导致局部细节的丢失。本发明根据人眼对亮度信息不明感而对对比度信息敏感的特性,对图像对比度进行压缩,在压缩特征图像灰度范围的同时,保留图像中尽可能多的视觉结构信息。具体方法如下:
设待压缩的图像为f,则采用下面的压缩函数对f进行压缩
其中,M、N分别为图像f的行数和列数;式中T和R是输入参数,分别是图像对比度的压缩率和压缩曲线在对数坐标轴下的拐点斜率;参数S、r和p通过联立下列各式求得
其中,
A2=((fmaxfmin)r+Lwa2r)(C-1)
A3=((fmaxfmin)r-Lwa2r)(C-1)
fmax=max(f)
fmin=min(f)
令对比度的压缩率T=0.5,压缩曲线在对数坐标轴下的拐点斜率R=1,然后取f分别为MRI图像的亮特征MW和暗特征MB,对其对比度进行压缩,得到对比度压缩后的亮特征CMW和暗特征CMB,
步骤四:分别用PET分解得到的亮特征图像PW和暗特征图像PB与MRI对比度压缩后的亮特征图像CMW和暗特征图像CMB进行最大值运算,得到新的亮特征图像NCMW和暗特征图像NCMW
其中,(x,y)为图像的像素坐标,如果将灰度图像看作二维标量函数,则Ο1是函数NCMW和PW支撑集的交集,则Ο2是函数NCMB和PB支撑集的交集。
步骤五:按下式得到融合后的亮度图像FI
FI=Pbase+NCMW-NCMB
其中Pbase为PET的亮度图像I经过多尺度双边形态学分解后得到的基频图像。
用FI代替PET图像经过IHS变换得到的亮度图像I,并与色调分量H和饱和度分量S一起通过IHS逆变换得到融合后的彩色图像F。从而达到融合高分辨率灰度MRI图像和低分辨率彩色PET图像,得到高分辨率彩色图像的目的。
Claims (1)
1.一种基于多尺度形态学双边滤波分解和对比度压缩的MRI与PET图像融合方法,其特征在以下五个步骤:
步骤一:采用双立方插值的方法对PET图像的R、G、B三通道分别进行插值,使其分辨率和MRI图像的相同,然后对其进行IHS变换、得到I、H、S三通道,其中I表示图像的亮度、H代表图像光谱的波长、S代表光谱的纯度;
步骤二:分别对MRI图像和PET的亮度图像I进行多尺度形态学双边滤波分解;设分解的层数为n,待分解的图像为f0则第i层分解为:
Δi1=max(Wi-1*G(di-1)),i>1的情形;
Δi2=max(Bi-1*G(di-1)),i>1的情形;
di<di-1
Wi=fi-fi-1
Bi=fi-1-fi
其中,Wi、Bi、fi分别为第i层分解得到的亮特征图像、暗特征图像和基频图像,Wi-1、Bi-1、fi-1为第i-1层分解得到的亮特征图像、暗特征图像和基频图像;Bd,Δ是自适应各向同性扁平结构元素,定义如下
其中,d为结构尺度因子,Δ为灰度尺度因子;Δi1、Δi2为第i次分解的灰度尺度因子,di-1、di分别为第i-1和i次分解的结构尺度因子;G(di-1)标准差为di-1的高斯函数,*为卷积运算符,定义如下
其中,(x,y)、(m,n)分别为卷积前后图像的像素坐标,M、N为被卷积图像的行列像素数;max()和min()分别为取最大值和最小值函数;·、和分别为数学形态学开运算、闭运算、膨胀和腐蚀运算符,定义如下
其中,(x,y)是图像像素坐标,(u,v)为结构元素中的像素坐标;
对于第一层分解,当i=1时,取
W1=f1-f0
B1=f0-f1
其中,M、N为图像f0的尺寸;
按照步骤一和步骤二的方法,对输入图像进行n层分解,得到n组基频图像、亮特征图像和暗特征图像,分别为f1、f2……fn,W1、W2……Wn,B1、B2……Bn;设形态学双边滤波分解得到的亮特征、暗特征和基频图像分别为W、B和Base,则其取值如下
Base=fn
令f0为MRI图像,对其进行多尺度形态学双边滤波分解,得到对应的亮特征图像MW、暗特征图像MB;然后令f0为PET的亮度图像I,对其进行多尺度形态学双边滤波分解,得到对应的亮特征图像PW、暗特征图像PB和基频图像PBase;
步骤三:对提取的MRI亮特征图像和暗特征图像分别进行对比度压缩,设待压缩的图像为f,则采用下面的压缩函数对f进行压缩
其中,M、N分别为图像f的行数和列数;式中T和R是输入参数,分别是图像对比度的压缩率和压缩曲线在对数坐标轴下的拐点斜率;参数S、r和p通过联立下列各式求得
其中,
A2=((fmaxfmin)r+Lwa2r)(C-1)
A3=((fmaxfmin)r-Lwa2r)(C-1)
fmax=max(f)
fmin=min(f)
令f分别为MRI图像的亮特征MW和暗特征MB,得到对比度压缩后的亮特征CMW和暗特征CMB;
步骤四:分别用PET分解得到的亮特征图像PW和暗特征图像PB与MRI对比度压缩后的亮特征图像CMW和暗特征图像CMB进行最大值运算,得到新的亮特征图像NCMW和暗特征图像NCMB
其中,(x,y)为图像的像素坐标,将灰度图像看作二维标量函数,则Ο1是函数CMW和PW支撑集的交集,则Ο2是函数CMB和PB支撑集的交集;
步骤五:按下式得到融合亮度图像FI
FI=Pbase+NCMW-NCMB
其中,Pbase为PET的亮度图像I经过多尺度形态学双边滤波分解后得到的基频图像;
用FI代替PET图像经过IHS变换得到的亮度图像I,并与色调分量H和饱和度分量S一起通过IHS逆变换得到融合后的彩色图像F,从而达到融合高分辨率灰度MRI图像和低分辨率彩色PET图像,得到高分辨率彩色图像的目的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374485.XA CN106056565B (zh) | 2016-05-31 | 2016-05-31 | 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374485.XA CN106056565B (zh) | 2016-05-31 | 2016-05-31 | 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106056565A CN106056565A (zh) | 2016-10-26 |
CN106056565B true CN106056565B (zh) | 2018-11-09 |
Family
ID=57172295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610374485.XA Active CN106056565B (zh) | 2016-05-31 | 2016-05-31 | 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106056565B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110675325B (zh) * | 2018-07-03 | 2023-04-28 | 柯鑫 | 一种增强眼底图像的方法和装置 |
CN110335225B (zh) * | 2019-07-10 | 2022-12-16 | 四川长虹电子系统有限公司 | 红外光图像与可见光图像融合的方法 |
CN111445419A (zh) * | 2020-04-09 | 2020-07-24 | 南京图格医疗科技有限公司 | 一种基于数学形态学的医学内窥镜图像增强方法 |
CN111539954B (zh) * | 2020-05-25 | 2024-01-23 | 国网湖南省电力有限公司 | 采用x射线数字影像特征识别电缆缓冲层缺陷的方法、系统及介质 |
CN115699079A (zh) * | 2020-07-07 | 2023-02-03 | 中加健康工程研究院(合肥)有限公司 | 图像处理方法、装置、计算设备和介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184534A (zh) * | 2011-05-25 | 2011-09-14 | 北京航空航天大学 | 一种利用多尺度高帽选择变换的图像融合方法 |
EP2413163A2 (en) * | 2010-04-30 | 2012-02-01 | Industrie University Cooperation Foundation Sogang University | Method for removing noise of PET signal using filtering in PET-MRI fusion device and PET system in PET-MRI fusion device using the same |
CN103679664A (zh) * | 2013-12-30 | 2014-03-26 | 北京航空航天大学 | 一种利用数学形态学交替滤波器的能够保留图像细节的增强方法 |
CN104376536A (zh) * | 2014-11-20 | 2015-02-25 | 北京航空航天大学 | 一种利用数学形态学高帽选择变换构造中心算子的图像增强方法 |
-
2016
- 2016-05-31 CN CN201610374485.XA patent/CN106056565B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2413163A2 (en) * | 2010-04-30 | 2012-02-01 | Industrie University Cooperation Foundation Sogang University | Method for removing noise of PET signal using filtering in PET-MRI fusion device and PET system in PET-MRI fusion device using the same |
CN102184534A (zh) * | 2011-05-25 | 2011-09-14 | 北京航空航天大学 | 一种利用多尺度高帽选择变换的图像融合方法 |
CN103679664A (zh) * | 2013-12-30 | 2014-03-26 | 北京航空航天大学 | 一种利用数学形态学交替滤波器的能够保留图像细节的增强方法 |
CN104376536A (zh) * | 2014-11-20 | 2015-02-25 | 北京航空航天大学 | 一种利用数学形态学高帽选择变换构造中心算子的图像增强方法 |
Non-Patent Citations (2)
Title |
---|
MRI and PET image fusion by combining IHS and retina-inspired models;S Daneshvar et al.;《Information Fusion》;20090520;第114-123页 * |
The multiscale directional bilateral filter and its application to multisensor image fusion;Jianwen Hu et al.;《Information Fusion》;20110126;第196-206页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106056565A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056565B (zh) | 一种基于多尺度形态学双边滤波分解和对比度压缩的mri与pet图像融合方法 | |
CN111583123A (zh) | 一种基于小波变换的融合高低频信息的图像增强算法 | |
JP5628306B2 (ja) | コントラスト改善 | |
CN110232661B (zh) | 基于Retinex和卷积神经网络的低光照彩色图像增强方法 | |
CN103345732B (zh) | 基于Contourlet变换的PCNN图像增强算法及装置 | |
Anand et al. | Mammogram image enhancement by two-stage adaptive histogram equalization | |
CN103034986A (zh) | 一种基于曝光融合的夜视图像增强方法 | |
CN111260580B (zh) | 图像去噪方法、计算机装置及计算机可读存储介质 | |
JP2016046706A (ja) | 画像処理装置、画像処理方法、記憶媒体及びプログラム | |
CN113658057A (zh) | 一种Swin Transformer微光图像增强方法 | |
CN110599418B (zh) | 一种变换域融合的全局色调映射方法 | |
Joshi et al. | Multi-focus image fusion using non-local mean filtering and stationary wavelet transform | |
Liu et al. | Enhancement of low illumination images based on an optimal hyperbolic tangent profile | |
CN117975224B (zh) | 基于指数特征幂次函数的遥感影像融合方法、系统及介质 | |
CN112150356A (zh) | 基于级联框架的单幅压缩图像超分辨率重建方法 | |
Upendra Raju et al. | Dual images in reversible data hiding with adaptive color space variation using wavelet transforms | |
Yu et al. | Decolorization algorithm based on contrast pyramid transform fusion | |
Kumar et al. | A two-level hybrid image fusion technique for color image contrast enhancement | |
Jagtap et al. | Reliable and robust low rank representation based noisy images multi-focus image fusion | |
CN116630198A (zh) | 一种结合自适应伽马校正的多尺度融合水下图像增强方法 | |
Deshmukh et al. | A comparative study of different image fusion techniques for tone-mapped images | |
Sharma et al. | Contrast enhancement using pixel based image fusion in wavelet domain | |
Nikam et al. | Color image enhancement using daubechies wavelet transform and HIS color model | |
Suji et al. | Image fusion-novel approach of contrast enhancement by discrete shearlet transform (DST) | |
Tojo et al. | Image denoising using multi scaling aided double decker convolutional neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |