CN101889295A - 生成多尺度对比度增强的图像的方法 - Google Patents
生成多尺度对比度增强的图像的方法 Download PDFInfo
- Publication number
- CN101889295A CN101889295A CN2008801200385A CN200880120038A CN101889295A CN 101889295 A CN101889295 A CN 101889295A CN 2008801200385 A CN2008801200385 A CN 2008801200385A CN 200880120038 A CN200880120038 A CN 200880120038A CN 101889295 A CN101889295 A CN 101889295A
- Authority
- CN
- China
- Prior art keywords
- image
- yardstick
- enhancing
- conversion difference
- value
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 38
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims description 84
- 230000002708 enhancing effect Effects 0.000 claims description 46
- 238000010586 diagram Methods 0.000 claims description 33
- 238000000354 decomposition reaction Methods 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 8
- 230000015556 catabolic process Effects 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 abstract 2
- 230000004048 modification Effects 0.000 abstract 2
- 238000005303 weighing Methods 0.000 abstract 1
- 238000005070 sampling Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 230000001419 dependent effect Effects 0.000 description 6
- 230000007704 transition Effects 0.000 description 6
- 238000003672 processing method Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 206010047571 Visual impairment Diseases 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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/10024—Color image
-
- 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
一个或者多个尺度的图像的至少一个近似图像被创建。通过按像素将尺度s的近似图像的值与所述近似图像的经转换的形式的值相减来创建转换差图像。非线性修正被应用于所述(一个或者多个)转换差图像的值并且特定尺度的至少一个增强的中心差图像通过合并该尺度的或者具有权重Wi,j的更小尺度的被修正的转换差图像而被计算。通过将重构算法应用于所述增强的中心差图像来计算增强的图像。所述转换差图像的值的非线性修正由方向图的值来控制,所述方向图为每个像素包括所关心的局部方位。另外地或者可替代地,通过用由所述方向图控制的权重对所述增强的转换差进行各向异性的加权来计算至少一个增强的中心差图像。
Description
技术领域
本发明涉及用于增强由数字信号表示的图像的图像质量的方法。
背景技术
通常,由数字信号表示的图像(诸如医学图像)在显示或者硬拷贝记录期间或者之前须经图像处理。
灰度值像素到适合用于再现或者显示的值的转换可以包括多尺度图像处理方法(也被称为多分辨率图像处理方法),图像的对比度借助于该方法而被增强。
根据这种多尺度图像处理方法,通过应用以下步骤来处理由像素值的阵列表示的图像。首先,原始图像被分解为多个尺度的细节图像的序列以及偶尔的残留图像。接着,通过将至少一个转换函数应用于这些像素值来修正所述细节图像的像素值。最后,通过将重构算法应用于所述残留图像以及所述被修正的细节图像来计算被处理的图像。
存在对所述转换函数的性能的限制。如果所述转换函数过度地非线性,则图像中的灰度值转变可能被扭曲到外观变得不自然的程度。所述扭曲在显著的灰度级转变的附近更加明显,其可能导致阶跃边缘处的过冲(overshoot)以及面向强阶跃边缘的低变化的区域中的均匀性的损失。产生假象的风险对于CT图像变得更加显著,因为它们具有更急剧的灰度级转变,例如在软组织与造影剂的交界处。人们必须小心地在CT图像上使用多尺度技术。
产生对比度增强的图像同时保持边缘转变的形状的多尺度对比度增强算法已经在2006年12月11日提交的共同待定的欧洲专利申请06125766.3中被描述。
在该方法的一个实施例中,以一个或者多个尺度创建所述图像的至少一个近似图像的转换差图像。,非线性地修正转换差图像。随后,特定尺度的至少一个增强的中心差图像通过合并该尺度的或者更小尺度的被修正的转换差图像来计算。从所述图像获得的空间局部化现象可以被用于创建增强的中心差图像。最后,通过将重构算法应用于所述增强的中心差图像来计算增强的图像。
一般而言,各向同性的滤波器被用在分解和重构过程中并且全方向的增强被应用于细节图像中的系数。各向同性的设计概念在其中图像统计是平稳的(stationary)那些情况下是合理的,意思是图像中的每个小块(patch)都通过与图像的每个其他小块相同的随机过程而被生成。然而,如果观看图像的如下区域:即其中边缘或者线或者均匀区域可能是明显的,则如下是清楚的:基本过程不是平稳的并且从小块到小块变化。
本发明的目的是克服现有技术的限制。
这种限制是没有能力促进对图像特征选择性的检测和增强,例如染色体图像,其中对于指定的DNA分析人们想要以特定尺度并且在某一方向和位置上增强染色体带。
发明内容
通过具有权利要求1所提出的特定的特征的方法来实现上述方面。用于本发明的优选实施例的特定的特征在从属权利要求中被提出。
在本发明的上下文中,特定的术语定义如下:
多尺度分解机制:
图像的多尺度(或者多分辨率)分解是以多个尺度计算所述图像的细节图像的过程。多尺度分解机制一般涉及用于计算所述细节图像的滤波器组。众所周知的技术举例来说有:拉普拉斯金字塔(laplacianpyramid)、伯特金字塔(Burt Pyramid)、拉普拉斯堆(Laplacian stack)、小波分解、QMF滤波器组...。
近似图像:
近似图像是以相同或者更大的尺度,或者以相同或者更低的分辨率表示原始灰度值图像的灰度值图像。特定尺度的近似图像与其中该尺度的所有细节已被省略的原始灰度值图像等同(Mallat S.G.,″多分辨率信号分解理论:小波表示″,IEEE模式识别及机器智能学报,1989年7月,卷11,编号7(″A Theory for Multisolution Signal Decomposition:The Wavelet Representation″,IEEE Trans.On Pattern Analysis andMachine Intelligence,Vol.11,no.7,July 1989))。
细节图像:
细节图像被定义为表示某一尺度的近似图像与更小尺度的近似图像之间的差的像素图。
转换算子:
转换算子是生成所述细节像素值的按像素的修正的算子,作为创建所述灰度值图像的对比度增强的形式的中间步骤。举例来说,这样的算子已经在欧洲专利EP 527 525中被描述。由转换函数定义所述修正并且可以例如将所述修正实现为查找表或者实现为乘法放大。
转换差图像:
尺度s的转换差图像是尺度s的近似图像的每个像素中的基本对比度(elementary contrast)的测量。可以通过求该尺度s的近似图像与经转换的形式之间的差来计算它们。对基本对比度的其他计算是可能的,例如像素与相邻像素的比可以被用于在处理步骤之前是指数变换并且在所述处理步骤之后是对数变换的情况中。
中心差图像:
通过将合并算子(例如求和)应用于所述转换差图像来计算中心差图像。所述合并算子可以是所述转换差图像中对应的像素值的线性或者非线性函数。
方向图:
特定尺度的方向图是对典型地与显著的图像结构相关联的突出的或者局部主导的方位的图像内的每个像素的表示。
方向图的可能的表示是向量图、极坐标表示(每个像素中的优选方位的角度和幅度)、笛卡尔坐标表示(每个像素中的优选方位的水平和垂直分量)。(在更简单的表示中,只有角度被定义并且方位向量的幅度被除去)。
计算方向图的许多实现已经被描述,从诸如第n阶计量坐标(then-th order gauge coordinate)、等照度线(在图像中连接相等强度的点的线)的方位的复杂机制到诸如水平和垂直梯度图像的合并的更直接的实现变化。
所述方向图也可以从细节图像而被计算出。例子是通过边缘小波生成的细节信息,其中所述细节信息以不同尺度表示边缘的局部化。
基于被增强并且被合并到中心差中的转换差图像的多尺度图像处理产生如下的可能性:调节对存在于所述数字图像中的方位的增强。
各向异性的多尺度图像处理基于如上定义的方向图,其在特定尺度上为每个像素定义所关心的优选的局部方位。
可以从所述近似图像以及从所述细节图像来计算出所述方向。
进一步,实现所述各向异性的多尺度图像增强的两个方法被描述:转换差图像的可控制的增强以及增强的转换差图像的各向异性的加权。两种实现可以独立于彼此而被应用或者一起被应用以从所述转换差图像来产生增强的中心差。
本发明一般地被实现为计算机程序产品,当其在计算机上被运行时,所述计算机程序产品适用于执行权利要求中的任意一项权利要求所述的方法并且被存储在计算机可读介质上。
本发明的方法可以被应用于增强医学图像(诸如乳房摄像图像、通过计算层析成像得到的图像等)的图像质量。
本发明的另外的优点和实施例从以下的描述及附图将变得显而易见。
附图说明
图1示出多分辨率图像处理方案,其中从相同尺度的近似图像来计算出中心差图像,
图2是图1的高阶(advanced)增强功能块的细节图,
图3和5示出根据本发明的多分辨率图像处理方法的不同实现,
图4示出图3所示的多分辨率图像处理方法的图像增强步骤,
图6示出图5所示的多分辨率图像处理方法的图像增强步骤,
图7是关于用在上述图像中的符号的图列,
图8示出通过使转换差的增强依赖于它们的方向的对增强的中心差的计算,
图9示出通过使其权重方向非独立的对增强的中心差的计算。
具体实施方式
本发明的对比度增强算法适用于所有多尺度细节表示方法,通过应用逆变换可以根据其来计算原始图像。
其对于可以被计算为转换差图像的加权和的可反转多尺度细节表示是适用的。
所述转换差图像的加权因子和转换偏置可以以如下方式从所述多尺度分解中被减除:使所述转换差图像的最终的加权和与所述细节像素值一致。
对于这些多尺度细节表示,可以通过在所述加权和被计算之前将转换算子应用于所述转换差来增强对比度。
为计算转换差图像的加权和,相同尺度(或者分辨率水平)的近似图像或者更小尺度(或者更精细的分辨率水平)的近似图像可以被使用。
现有技术状况的多尺度对比度增强算法将图像分解为多尺度表示,所述多尺度表示包括以多个尺度表示细节的细节图像以及残留图像。
一些重要的多尺度分解是小波分解、拉普拉斯-高斯(或者LoG分解)、高斯差分(或者DoG)分解以及伯特金字塔。
通过应用随后是子采样步骤的高通及低通滤波器级联来计算所述小波分解。
高通滤波器从特定尺度的近似图像中提取细节信息。
在所述伯特金字塔分解中,通过减去尺度k+1的近似图像的过采样(upsampled)形式从尺度k的近似图像提取出细节信息。
在如EP 527 525中所公开的一个现有技术状况的方法中,通过随后是多尺度重构的细节图像中的像素值的转换来创建图像的对比度增强的形式。
多尺度分解的所有上述实现具有共同的特性。通过合并滑动邻域中的像素值可以从近似图像计算出细节图像中的每个像素值。
在上述情况中,合并函数是加权求和。
对于小波分解,尺度k的细节图像中的像素值被计算为:
dk+1=↓(hd*gk)
gk+1=↓(ld*gk)
其中hd是高通滤波器,ld是低通滤波器,*是卷积算子并且↓是亚采样(subsampling)算子(即省去每第二行和列(every second row andcolumn))。
对于小波重构,尺度k的增强的近似图像被计算为:
hk=lr*(↑hk+1)+hr*(↑f(dk+1))
其中hr是高通滤波器,lr是低通滤波器并且↑是过采样算子(即在任两行和列(any two rows and columns)之间插入具有0值的像素)。
对于伯特分解,尺度k的细节图像中的像素值被计算为:
dk=gk-4g*(↑gk+1)
或者
dk=gk-4g*(↑(↓(g*gk)))
或者
dk=(1-4g*(↑(↓g)))*gk
其中g是高斯低通滤波器并且l是恒等算子。
对于伯特重构,尺度k的增强的近似图像被计算为:
hk=4g*(↑hk+1)+f(dk)
其中f(x)是转换算子。
作为加权和的多尺度细节像素值
假如在伯特多尺度分解中具有系数wk,l的5×5高斯滤波器被使用,其中k=-2,...2并且l=-2,...2,所述亚采样算子去除每第二行和列并且所述过采样算子在任两行和列之间插入具有值0的像素。
在近似图像gk+l中处于位置i,j的像素被计算为:
在过采样的图像uk中处于位置i,j的像素被计算为:
在过采样的、经平滑的图像guk中处于位置i,j的像素被计算为:
最后,在细节图像dk中处于位置i,j的像素被计算为:
一般而言,在细节图像dk中处于位置i,j的像素可以被计算为相同或者更小尺度k,k-1,k-2,...的近似图像中的像素的加权和:
其中l∈{0,...,k}并且r=亚采样因子(1-k)
因为
所以在细节图像dk中处于位置i,j的像素可以被计算为:
项gl(ri,rj)-gl(ri+m,rj+n)被称为转换差。
其表示了近似图像中的中心像素和相邻像素之间的像素值的差。它是局部对比度的测量。
所述转换差的加权和被称为中心差ck(i,j)。
同样地,对于小波分解,细节图像中的像素值可以被计算为近似图像中的像素的加权和。通过将具有滤波器系数hm,n的高通滤波器h应用于近似图像gk来计算细节图像dk+1中的像素值。该滤波操作之后是去除每第二行和列的亚采样算子。
dk+1=↓(h*gk)
在细节图像dk+1中处于位置i,j的像素被计算为:
因为对于高通滤波器,滤波器系数的和等于0:
或者
所以在细节图像dk+1中处于位置i,j的像素可以被计算为:
项gk(2i,2j)-gk(2i+m,2j+n)表示转换差。
以相似的方式,可以被证明的是:其他多尺度分解方法中的细节图像也可以被表示为转换差图像的合并。
转换操作
在像EP 527 525中所公开的一个现有技术状况的方法中,通过将转换算子f(x)应用于细节图像dk来获得对比度增强,或者,等同地为:
这样的转换算子的例子是S形函数(sigmoid function)。这样的转换算子的另一个例子是像在EP 525 527中所公开的一个对比度增强函数。所述转换算子的形状取决于对增强的特定需求,其意图是与高值细节像素相比更多地放大低值细节像素。
转换步骤可能在被重构的、对比度增强的图像中引起边缘转变的形状的变形(deformation)。原因是所述转换函数的非线性。
一般而言,下列各项适用于非线性函数:
f(x+y)≠f(x)+f(y)
现有技术状况的算法首先将细节图像dk中的像素值计算为加权和并且此后应用所述转换步骤。
通过将细节图像dk中的像素值重写为转换差的加权和,有可能在求和之前而不是之后应用所述转换步骤。
现在通过将所述转换步骤应用于所述转换差来得到对比度增强:
以此方式,所述边缘转变的形状在对比度被增强的、重构的图像中更好地被保持。
如果对于每个尺度k,从满分辨率图像g0计算出该尺度的细节图像,并且增强被应用于所述中心差,则所述边缘转变的形状在重构之后最佳地被保持。
本发明的不同实现在图1、3和5中被示出,其中对应的增强步骤在图2、4和6中被示出。
图1示出多尺度图像处理方案,其中从相同尺度的近似图像计算出所述中心差图像。该图示出配备有控制路径的用于多分辨率图像处理的扩展方案。从不同尺度的近似图像,所述中心差图像被计算出并且在受控的增强块E中被增强。块E的输出经受重构过程。
图2示出利用从近似图像获得的方向图的受控的增强块E。
功能块Ak
这些滤波块可以被用于增强近似图像或者从所述近似图像计算出某些特性,典型地为水平和垂直梯度图像。这些特性可以被用于计算所述方向图。
在特定的实施例中,所述近似图像在方向图从其被计算出之前被增强。举例来说,这可以通过对所述近似图像进行额外的平滑以减少噪声及不相关的小的图像结构对方向图计算的影响来执行。
同样地,其他特性可以从所述近似图像计算出,例如水平和垂直梯度图像。在下一个步骤中,即在方向图生成器中(将在下文中被描述),两个图像可以被合并以计算所述方向图。
特定尺度的方向图对于图像内的每个像素是显著的图像结构的突出的或者局部主导的方位的表示。
方向图的可能表示已经在上文中被描述。
各向异性增强
各向异性增强通过(1)所述转换差图像的受控的增强,由此所述控制取决于所述方向图的内容,或者(2)通过增强的转换差图像的各向异性的加权(系数wj)来执行。两种实现可以独立于彼此被应用或者可以一起被应用于从所述转换差图像来创建出增强的中心差。
转换差图像的可控制的增强
每个增强的中心差被计算为增强的转换差图像的合并。
转换差是近似图像中中心像素与在其局部邻域中的像素的像素值的差。
通过使用这两个像素相对于彼此的空间方向,可以将所述增强用于根据所关心的预先定义的方向更多或者更少地修正转换差。
例子(参见图8--3D图像中的不同的层表示不同的值):
假定噪声边缘存在于所述图像内的所关心的区域中。
为了合并边缘的增强而同时减小噪声的影响,可以增强具有与边缘方向垂直的方向的转换差并且衰减平行于所述边缘的转换差。
后一种差图像表示不想要的噪声图像结构,而第一种表示在所述边缘上的局部对比度。具有在这两者之间的方向的转换差图像将经历中间的增强。
从近似图像gk-m开始,通过合并图像gk-m的水平和垂直梯度来计算方向图M。为了减小噪声对梯度计算的影响,使用例如中值滤波器对所述近似图像gk-m进行滤波,产生被滤波的图像Ak(gk-m)。使用水平和垂直方位上的核(kernel){1,-1}来计算所述图像内的垂直和(en)水平一阶梯度。
如果所述方向图M表示例如显著的图像结构的局部方位(或者垂直于其),所述图被计算为垂直和水平梯度的比的反正切函数:
该方向图M被用作所述LUT算子的控制输入。
LUT算子具有两个输入,转换差以及方向图,并且生成增强的转换差作为输出。这种可控制的LUT算子LUT(dk,M(i,j))的例子是转换算子f(dk)与附加的依赖于方向的放大函数Amp(θ)的合并,其中θ是方向图M(i,j)的函数。
所述转换算子f(dk)的例子是S形函数以及像在EP 525 527中所公开的一个对比度增强函数。所述转换算子的形状取决于所述增强的特定的要求,其意图是与高值转换差相比更多地放大低值转换差。
取决于方向的放大函数典型地返回沿优选方位的最大放大因子、用于垂直于优选方位的最小放大因子以及用于中间的方向的逐渐变化的放大因子。
这种依赖于方向的放大函数的例子是与常数放大因子b相乘的余弦函数。作为该函数的输入,可以求由所述方向图M(i,j)规定的所关心的局部方位和处于位置(m,n)的转换差相对于所关心的中心像素的方向之间的差:
Amp(θ)=bcos(θ)
其中
θ=M(i,j)-angle(dk(i+m,j+n)并且angle(dk(i+m,j+n)=tan-1(n/m)
可控制的LUT算子的更高阶的实现是多维度的查找表。
所述LUT算子的结果是增强的转换差以及所述转换差沿由M(i,j)所规定的方位的附加的增强(在图8中由增强的转换差em,n的灰色背景的暗度来指示)。
接着使用各向同性的权重wm,n合并增强的转换差来计算处于位置i,j的增强的中心差。
权重wm,n被计算使得未增强的转换差的加权和刚好匹配将通过伯特分解的直接应用所生成的细节像素值。
增强的转换差图像的各向异性加权
为从处于每个位置(x,y)的转换差图像计算出中心差,求所述转换差的加权和。在2006年12月11日提交的欧洲专利申请06 125 766.3中,权重是各向同性的并且没有定义优先的方位。
通过根据它们相对于中心权重的空间方向来改变权重,可以在总的和中改变沿预先定义的方向的转换差的重要性。
例子(参见图9-3D图像中的不同的层表示不同的值):
假定噪声边缘存在于所述图像内的所关心的区域中。
如果在减小噪声的同时希望得到该边缘的增强,则可以将相对较高的权重应用于具有与所述边缘方向垂直的方向的增强的转换差并且将相对较低的权重应用于平行于所述边缘的增强的转换差。后一种转换差表示不想要的噪声图像结构,而第一种表示在所述边缘上的局部对比度。具有在这两者之间的方向的增强的转换差将以中间的权重被考虑在内。
根据近似图像gk-m对方向图M的计算在第一实施例中详细地被解释。
代替将所述方向图M用作所述LUT算子的控制输入,所述方向图M也可以被用作控制参数来产生各向异性的权重wm,n。
在该实施例中,LUT算子仅取决于未增强的转换差,例如像上文所规定的转换算子f(dk)。
对于图像内的每个位置i,j各向同性的权重wm,n相对于它们到如所述方向图M(i,j)规定的所关心的局部方位的方向而被修正。
作为修正算子,放大函数的机制如在第一实施例中所解释的那样被使用。
所述各向同性的权重wm,n被乘以依赖于方向的放大函数,其典型地返回沿优选方位的最大乘法因子、用于垂直于所述优选方位的方向的最小乘法因子以及用于所述中间的方向的逐渐变化的乘法因子(在图9中用权重wm,n的灰色背景的暗度来指示)。
这种依赖于方向的放大函数的例子是余弦函数。作为该函数的输入,可以求由所述方向图M(i,j)规定的所关心的局部方位和处于位置(m,n)的权重相对于所述中心位置(0,0)的方向之间的差:
Amp(θ)=cos(θ)
其中
θ=M(i,j)-angle(wk(m,n)并且angle(wk(m,n)=tan-1(n/m)
此外,所述各向异性的权重被归一化使得所述权重的和等于1。
在最后一个步骤中,通过使用所述各向异性的权重合并增强的转换差来计算中心差。这为沿由所述方向图M(i,j)规定的所关心的局部方位被定向的增强的转换差产生相对较高的权重。
图8和9所示的增强方法的合并也是可能的。
Claims (11)
1.一种增强由数字信号表示的图像的对比度的方法,其中
-通过将多尺度分解算法应用于所述图像来创建一个或者多个尺度的至少一个近似图像,在表示所述图像的灰度值的尺度的近似图像中该尺度的所有细节已被省略并且由对应的细节图像表示,
-通过按像素将尺度s的近似图像的值与所述近似图像的经转换的形式的值相减来创建转换差图像,
-非线性修正被应用于所述(一个或者多个)转换差图像的值,
-特定尺度的至少一个增强的中心差图像通过合并该尺度的或者具有权重wi,j的更小尺度的所述被修正的转换差图像而被计算,其中根据所述分解算法来选择所述权重wi,j以便如果通过非线性修正的所述增强没有被执行,则所述图像的多尺度分解的等同物将被获得以便所述至少一个中心差图像等于所述特定尺度的细节图像,
-通过将反转所述图像分解算法的图像重构算法应用于所述图像分解来计算增强的图像,所述图像分解在所述特定尺度的所述细节图像已被所述至少一个增强的中心差图像替代之后产生,
其特征在于
所述转换差图像的值的所述非线性修正由方向图的值控制,所述方向图在特定尺度上为每个像素包括所关心的局部方位。
2.一种增强由数字信号表示的图像的对比度的方法,其中
-创建一个或者多个尺度的至少一个近似图像,在表示所述图像的灰度值的尺度的近似图像中该尺度的所有细节已被省略,
-通过按像素将尺度s的近似图像的值与所述近似图像的经转换的形式的值相减来创建转换差图像,
-非线性修正被应用于所述(一个或者多个)转换差图像的值,
-特定尺度的至少一个增强的中心差图像通过合并该尺度的或者具有权重wi,j的更小尺度的所述被修正的转换差图像而被计算,
-通过将重构算法应用于所述增强的中心差图像来计算增强的图像,
其特征在于,通过用由所述方向图控制的权重对被修正的转换差图像的值进行各向异性的加权来计算所述增强的中心差图像中的至少一个。
3.根据权利要求1或者2所述的方法,其特征在于,其中所述转换差图像的值的所述非线性修正由方向图的值控制,所述方向图在特定尺度上为每个像素包括所关心的局部方位并且所述被修正的转换差图像中的至少一个经受由所述方向图控制的各向异性加权。
4.根据权利要求1或者2所述的方法,其特征在于,其中特定尺度的转换差图像从相同尺度的近似图像被计算出。
5.根据权利要求1或者2所述的方法,其特征在于,其中所有所述转换差图像从所述原始图像被计算出。
6.根据权利要求1或者2所述的方法,其特征在于,其中尺度k的转换差图像从尺度m的近似图像被计算出,其中m表示尺度1和尺度k-1之间的尺度。
7.根据权利要求1或者2所述的方法,其特征在于,其中所述中心差图像与所述多尺度细节图像一致。
8.根据前述权利要求中任意一项权利要求所述的方法,其特征在于,其中所述图像是乳房摄像图像。
9.根据前述权利要求中任意一项权利要求所述的方法,其特征在于,其中所述图像是CT图像。
10.一种计算机程序产品,当其在计算机上被运行时,所述计算机程序产品执行前述权利要求中任意一项权利要求所述的方法。
11.一种计算机可读介质,所述计算机可读介质包括计算机可执行的程序代码,所述计算机可执行的程序代码适用于执行前述权利要求中任意一项权利要求所述的步骤。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP07118032A EP2048616A1 (en) | 2007-10-08 | 2007-10-08 | Method of generating a multiscale contrast enhanced image |
EP07118032.7 | 2007-10-08 | ||
PCT/EP2008/063235 WO2009047208A1 (en) | 2007-10-08 | 2008-10-02 | Method of generating a multiscale contrast enhanced image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101889295A true CN101889295A (zh) | 2010-11-17 |
CN101889295B CN101889295B (zh) | 2013-09-18 |
Family
ID=38943830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008801200385A Expired - Fee Related CN101889295B (zh) | 2007-10-08 | 2008-10-02 | 生成多尺度对比度增强的图像的方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8442340B2 (zh) |
EP (2) | EP2048616A1 (zh) |
CN (1) | CN101889295B (zh) |
WO (1) | WO2009047208A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104956396A (zh) * | 2012-12-24 | 2015-09-30 | Lsi公司 | 使用基于来自其它图像的信息的函数的泛函产生目标图像 |
CN105096247A (zh) * | 2014-05-09 | 2015-11-25 | 北京大学 | 图像插值方法和装置 |
CN105723416A (zh) * | 2013-11-18 | 2016-06-29 | 爱克发医疗保健公司 | 图像降噪方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2328124B1 (en) | 2009-11-25 | 2019-05-15 | Agfa Nv | Method of enhancing the contrast of spatially-localized phenomena in an image |
KR102244918B1 (ko) | 2014-07-11 | 2021-04-27 | 삼성전자주식회사 | 시인성을 향상시키고 전력 소모를 줄일 수 있는 디스플레이 컨트롤러와 이를 포함하는 디스플레이 시스템 |
CN113039581A (zh) * | 2018-09-14 | 2021-06-25 | 恩维医疗公司有限公司 | 三维对象的多尺度图像重建 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5196924A (en) | 1991-07-22 | 1993-03-23 | International Business Machines, Corporation | Look-up table based gamma and inverse gamma correction for high-resolution frame buffers |
EP0527525B1 (en) | 1991-08-14 | 1996-10-02 | Agfa-Gevaert N.V. | Method and apparatus for contrast enhancement |
US5526446A (en) * | 1991-09-24 | 1996-06-11 | Massachusetts Institute Of Technology | Noise reduction system |
US6788826B1 (en) * | 1998-11-10 | 2004-09-07 | Agfa-Gevaert | Method for correcting artefacts in a digital image |
JP4197392B2 (ja) * | 2000-10-17 | 2008-12-17 | 富士フイルム株式会社 | ノイズ抑制処理装置並びに記録媒体 |
EP1223553A3 (en) * | 2000-10-17 | 2003-09-24 | Fuji Photo Film Co., Ltd. | Apparatus for suppressing noise by adapting filter characteristics to input image signal based on characteristics of input image signal |
US20060072845A1 (en) * | 2002-12-18 | 2006-04-06 | Koninklijke Philips Electronics N.V. | Method of processing an input image by means of multi-resolution |
CN100420269C (zh) * | 2005-12-09 | 2008-09-17 | 逐点半导体(上海)有限公司 | 一种图像增强处理系统和处理方法 |
US8081832B2 (en) | 2006-12-11 | 2011-12-20 | Agfa Healthcare N.V. | Method of generating multiscale contrast enhanced image |
US8050513B2 (en) | 2006-12-11 | 2011-11-01 | Agfa Healthcare N.V. | Method of generating multiscale contrast enhanced image |
-
2007
- 2007-10-08 EP EP07118032A patent/EP2048616A1/en not_active Withdrawn
-
2008
- 2008-10-02 US US12/681,791 patent/US8442340B2/en active Active
- 2008-10-02 WO PCT/EP2008/063235 patent/WO2009047208A1/en active Application Filing
- 2008-10-02 CN CN2008801200385A patent/CN101889295B/zh not_active Expired - Fee Related
- 2008-10-02 EP EP08805016.6A patent/EP2198402B1/en active Active
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104956396A (zh) * | 2012-12-24 | 2015-09-30 | Lsi公司 | 使用基于来自其它图像的信息的函数的泛函产生目标图像 |
CN105723416A (zh) * | 2013-11-18 | 2016-06-29 | 爱克发医疗保健公司 | 图像降噪方法 |
CN105723416B (zh) * | 2013-11-18 | 2019-01-08 | 爱克发医疗保健公司 | 图像降噪方法 |
CN105096247A (zh) * | 2014-05-09 | 2015-11-25 | 北京大学 | 图像插值方法和装置 |
CN105096247B (zh) * | 2014-05-09 | 2018-02-06 | 北京大学 | 图像插值方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
EP2048616A1 (en) | 2009-04-15 |
US8442340B2 (en) | 2013-05-14 |
CN101889295B (zh) | 2013-09-18 |
EP2198402B1 (en) | 2019-06-12 |
EP2198402A1 (en) | 2010-06-23 |
US20100220911A1 (en) | 2010-09-02 |
WO2009047208A1 (en) | 2009-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101821775B (zh) | 用于生成多尺度的对比度增强的图像的方法 | |
CN101889295B (zh) | 生成多尺度对比度增强的图像的方法 | |
Kollem et al. | Improved partial differential equation-based total variation approach to non-subsampled contourlet transform for medical image denoising | |
CN101821773B (zh) | 增强图像对比度的方法 | |
CN101939763B (zh) | 生成多尺度对比度增强的图像的方法 | |
Voronin et al. | 3-D block-rooting scheme with application to medical image enhancement | |
Zou et al. | 3D filtering by block matching and convolutional neural network for image denoising | |
US8050513B2 (en) | Method of generating multiscale contrast enhanced image | |
US8081832B2 (en) | Method of generating multiscale contrast enhanced image | |
Pham | Kriging-weighted laplacian kernels for grayscale image sharpening | |
EP1933272B1 (en) | Method of generating a contrast enhanced image using multiscale analysis | |
US20220392029A1 (en) | Method and Apparatus for Contrast Enhancement | |
Siddiqi et al. | [Retracted] An Ensembled Spatial Enhancement Method for Image Enhancement in Healthcare | |
Zafeiridis et al. | A New Sharpening Technique for Medical Images using Wavelets and Image Fusion. | |
EP1933273B1 (en) | Generating a contrast enhanced image using multiscale analysis | |
CN101231746B (zh) | 生成多尺度对比增强图像的方法 | |
Tan et al. | A Novel Approach for Image Enhancement via Nonsubsampled Contourlet Transform | |
Azad et al. | Digitizing and sampling of image process | |
Wang et al. | A novel method based on discrete multiple wavelet transform to multispectral image fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190506 Address after: Belgian Mo Patentee after: Agfa Co.,Ltd. Address before: Belgian Mo Patentee before: AGFA HEALTHCARE |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130918 |
|
CF01 | Termination of patent right due to non-payment of annual fee |