CN104376529A - 一种基于glcm的灰度图像彩色化系统和方法 - Google Patents

一种基于glcm的灰度图像彩色化系统和方法 Download PDF

Info

Publication number
CN104376529A
CN104376529A CN201410685330.9A CN201410685330A CN104376529A CN 104376529 A CN104376529 A CN 104376529A CN 201410685330 A CN201410685330 A CN 201410685330A CN 104376529 A CN104376529 A CN 104376529A
Authority
CN
China
Prior art keywords
block
image
pixel
pixels
color
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
Application number
CN201410685330.9A
Other languages
English (en)
Other versions
CN104376529B (zh
Inventor
李超
王涛
盛浩
朱耿良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Hanfei Network Technology Co ltd
Original Assignee
SHENZHEN BEIHANG NEW INDUSTRIAL TECHNOLOGY RESEARCH INSTITUTE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SHENZHEN BEIHANG NEW INDUSTRIAL TECHNOLOGY RESEARCH INSTITUTE filed Critical SHENZHEN BEIHANG NEW INDUSTRIAL TECHNOLOGY RESEARCH INSTITUTE
Priority to CN201410685330.9A priority Critical patent/CN104376529B/zh
Publication of CN104376529A publication Critical patent/CN104376529A/zh
Application granted granted Critical
Publication of CN104376529B publication Critical patent/CN104376529B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/08Projecting images onto non-planar surfaces, e.g. geodetic screens
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • G06T7/49Analysis of texture based on structural texture description, e.g. using primitives or placement rules
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开一种一种基于GLCM的灰度图像彩色化系统和方法,包括图像分块模块、灰度共生矩阵特征提取模块、相似度匹配建立对应关系模块、颜色标记和初步修正模块、最优化着色模块;能够解决以前方法的不足并适应工程实际需求,该方法能加快图像彩色化速度,并能改进图像的彩色化效果。

Description

一种基于GLCM的灰度图像彩色化系统和方法
技术领域
本发明属于图像处理领域,涉及一种灰度图像彩色化方法,特别涉及一种基于GLCM(灰度共生矩阵,Gray-level co-occurrencematrix)的灰度图像彩色化系统和方法。
背景技术
图像是人类认识客观世界最常用的信息载体,图像信息除了包括景物形状、大小等信息更包括颜色信息。没有颜色信息的图像是一种信息表达不完全的形式。在现实中,夜视图像、铅笔手绘图像等都缺少丰富的色彩信息。特别是夜视领域中,微光和热成像技术是当前主流技术,这两种图像都是灰度图像,提供的细节信息有限,而人眼对颜色的分辨率远超对灰度层级的分辨率,如果利用彩色化处理将灰度图像中的细节信息凸现出来,可使人眼对图像的细节有更丰富的认识。老旧照片或铅笔涂鸦绘画中,利用彩色化能生成图像信息,使图像更生动形象。同一场景下的图像实现不同季节下的景色变化等。
实现图像彩色化,需要基于参考图像,在参考图像的映射下实现颜色渲染。基于参考图像的图像彩色化方法一般分为两大类:人工交互法和自动着色法。人工交互方法需要人为绘画涂抹灰度图像,使其变为涂鸦或带有颜色痕迹的图像,然后进行颜色的渲染传播,扩大至整幅图像。Anat Levin等人提出了基于优化扩展的彩色化方法,这种方法主要是基于一种图像分割方法,设立优化条件,寻求最优的图像分割,然后进行着色。这个方法是半自动的方法,所谓半自动方法就是说需要我们人为的在一副灰度图像上画出自己想要添加的颜色线条,根据我们确定的这个颜色线条,结合分割方法,进行着色。这类方法的缺点主要在于需要大量的人力手工标注图像,而且需要很好的观察与细致的工作才能实现好的效果,费时费力。
自动着色法是利用参考图像将色彩迁移到目标图像中。Welsh方法在过去Rein hard色彩传递方法的基础上进行改进,对目标图像的每个像素从源图像上匹配与之对应的像素,利用像素亮度值与邻域统计值作为判断准则,将参考图像的颜色信息传递到目标图像,同时保留目标图像的亮度信息。这类方法虽然减少人为干预,变为自动方法,但色彩渲染效果并不好,只能实现大致的色调迁移,而不能鲜明的区分出各种物体。
发明内容
本发明提供一种基于GLCM的灰度图像彩色化系统和方法,能够解决以前方法的不足并适应工程实际需求,该方法能加快图像彩色化速度,并能改进图像的彩色化效果。
本发明通过以下技术手段实现:
一种基于GLCM的灰度图像彩色化系统,包括图像分块模块、灰度共生矩阵特征提取模块、相似度匹配建立对应关系模块、颜色标记和初步修正模块、最优化着色模块;
所述的图像分块模块,采用像素块网格划分法,根据基本像素块进行图像划分;
所述的灰度共生矩阵特征提取模块,针对每个像素块,分别计算像素块的灰度均值、方差和0°方向、45°方向、90°方向、135°方向的灰度共生矩阵,提取灰度共生矩阵的对比度、相关性、能量、同质性作为特征描述子,组成18维特征向量描述图像块的纹理特征,综合所述各个特征形成表征像素块的特征描述子;
所述的相似度匹配建立对应关系模块,通过计算特征向量的欧氏距离表现像素块的近似程度;欧氏距离最小的灰度目标像素块,即为与之对应的最优匹配块;
所述的颜色标记和初步修正模块,选取像素块的中心窗口位置作为颜色迁移点,然后计算彩色图像块中心窗口内色彩的颜色均值作为迁移色,然后同样选取灰度图像块中心窗口区域作为迁移区域,进行色彩赋值,生成标记涂鸦图像;
所述的最优化着色模块,把颜色痕迹图像的颜色空间从RGB空间转换到YUV空间,构造包含所有像素颜色系数的稀疏矩阵,依据任意像素与其邻域像素颜色加权差的平方和最小原则构建方程组,此方程组用迭代法求解。
其中,所述的基本像素块数目准则是10*10,对于不能整分的图像,采取向下取整的方法。
其中,所述的两个像素块之间相似度匹配的计算除了欧氏距离的叠加外,还要进行灰度均值、方差和灰度共生矩阵特征的权重分配。
其中,所述的像素块的灰度共生矩阵特征、方差和灰度均值的对应权重分别取值0.2,0.1和0.1,对于灰度共生矩阵中各个特征描述子的权重取值分别取0.25,0.25,0.25,0.25。
其中,像素块的中心窗口和灰度目标图像中心窗口的选取规则为:假设{ri}和{tj}是参考图像与目标图像对应的最优匹配块,设彩色参考像素块的尺寸大小为M*N,则选取的中心窗口长度范围是[M/4,3M/4],宽度范围[N/4,3N/4]。
其中,所述的灰度共生矩阵匹配产生误配时,根据相邻像素具有相似的灰度值则具有相似颜色值的准则,进行结果修正,所述的修正过程是将计算灰度目标图像块的各中心窗口的灰度值,进行相邻九宫格区域的比对,如果发现灰度值相似,而迁移色不一致,则统计具有相似灰度值的相邻块的迁移色进行修正。
其中,所述的图像将以相邻的四方格或九宫格像素作为单个像素的方式进行着色,此单个像素记为四方格或九宫格像素块的超像素。计算四方格或九宫格像素块的亮度均值作为超像素的亮度值,此时原图像会缩减至超像素组成的新图像。新图像以超像素为彩色化的基本单元,将超像素代入最优化着色方法求解。求解得到的超像素的色度值作为四方格或九宫格所有像素的色度值。
一种基于GLCM的灰度图像彩色化系统的方法,实现步骤如下:步骤1:图像分块:
即采用像素块网格划分法,根据基本像素块对图像进行像素块划分,然后利用像素块进行特征提取和匹配;
步骤2:灰度共生矩阵特征提取:
即针对每个像素块,分别计算像素块的灰度均值、方差和0°方向、45°方向、90°方向、135°方向的灰度共生矩阵,提取灰度共生矩阵的对比度、相关性、能量、同质性作为特征描述子,一共组成18维特征向量描述图像块的纹理特征,综合所述各个特征形成表征像素块的特征描述子;
步骤3:相似度匹配建立对应关系:
即通过计算特征向量的欧氏距离表现像素块的近似程度,计算出欧氏距离最小的灰度目标像素块,就是与之对应的最优匹配块;
步骤4:颜色标记和初步修正:
即选取像素块的中心窗口位置作为颜色迁移点,然后计算彩色图像块中心窗口内色彩的颜色均值作为迁移色,然后同样选取灰度图像块中心窗口区域作为迁移区域,进行色彩赋值,生成标记涂鸦图像;步骤5:最优化着色:
依据若相邻像素具有相同的亮度值则应具有相同的颜色值的原理。即先把颜色痕迹图像的颜色空间从RGB空间转换到YUV空间,构造包含所有像素颜色系数的稀疏矩阵,依据任意像素与其邻域像素颜色加权差的平方和最小原则构建方程组,此方程组用迭代法求解。
其中,两个像素块之间相似度匹配的计算除了欧氏距离的叠加,还要进行灰度均值,方差和灰度共生矩阵特征的权重分配。
其中,像素块的灰度共生矩阵特征、方差和灰度均值的对应权重分别取值0.2,0.1和0.1;对于灰度共生矩阵中各个特征描述子的权重取值分别取0.25,0.25,0.25,0.25。
其中,所述的像素块的中心窗口和灰度目标图像中心窗口的选取规则为:假设{ri}和{tj}是参考图像与目标图像对应的最优匹配块,设彩色参考像素块的尺寸大小为M*N,则选取的中心窗口长度范围是[M/4,3M/4],宽度范围[N/4,3N/4],灰度目标图像中心窗口的选取规则与之相同。
其中,所述的图像将以相邻的四方格或九宫格像素作为单个像素的方式进行着色,此单个像素记为四方格或九宫格像素块的超像素。计算四方格或九宫格像素块的亮度均值作为超像素的亮度值,此时原图像会缩减至超像素组成的新图像。新图像以超像素为彩色化的基本单元,将超像素代入最优化着色方法求解。求解得到的超像素的色度值作为四方格或九宫格所有像素的色度值。
本发明与现有技术相比所具有的优点是:
(1)加快了图像彩色化速度,能够为灰度图像视频彩色化的实时性打下基础;
(2)图像彩色化效果较传统方法更好,结果图像色彩信息更丰富自然;
(3)实现简单,处理过程明确,使用图像颜色痕迹的自动迁移代替人工手动涂鸦标记,节省时间和劳动力,更容易实现计算机的自动处理。
附图说明
图1为本发明的系统架构图;
图2为本发明的方法流程图;
图3颜色标记结果修正示意图。
具体实施方式
以下将结合附图对本发明的实现过程进行详细描述。
一种基于GLCM的灰度图像彩色化系统,包括图像分块模块、灰度共生矩阵特征提取模块、相似度匹配建立对应关系模块、颜色标记和初步修正模块、最优化着色模块。
基于以上系统的一种新的图像自动彩色化方法,如图1所示,该方法利用具有相似内容的图像作为参考图像,先对参考图像和目标图像分别进行分块,其次计算各块图像灰度共生矩阵、灰度均值、方差等作为特征描述,然后计算各图像块特征向量的欧氏距离,得出最优的像素块匹配关系,接着根据对应关系进行颜色痕迹的迁移,最后实现最优化颜色传播。
其步骤如下:
步骤1:初始图像分块
针对目标图像和参考图像采用像素块网格划分法,基本像素块数目准则是10*10,然后根据图像尺寸大小再次细分,图像越大,分块数目越多。例如800*600的图像,按基本像素块划分为10行10列,每块尺寸大小80*60。对于2400*1200的图像,可划分为20行20列,每块大小120*60。但随着分块数目的增多,相应的匹配过程耗费时间也增多。因此,实验过程中,选取10行10列的分块准则。对于不能整分的图像,为简化运算,我们采取向下取整的方法。如768*512图像,先裁剪成760*510的尺寸,然后进行分块。因为图像的颜色分布呈区域连续性状态,因此进行小范围裁剪并不影响最终的彩色化结果。假如图像尺寸481*321,裁剪取整后为480*320,分块成10行10列,则每块大小48*32,之后利用像素块进行特征提取和匹配。
步骤2:GLCM特征提取
针对每个像素块,分别计算像素块的灰度均值、方差和0°,45°,90°,135°的4个方向的灰度共生矩阵,提取灰度共生矩阵的对比度、相关性、能量、同质性作为特征描述子,这样一共组成18维(1+1+4*4)特征向量描述图像块的纹理特征。
以下为各个特征的计算过程:
(1)灰度均值和方差
像素块的灰度均值是表征像素块的基本特征,图像中不同景物一般也会有不同的灰度值,因此灰度均值能显著区分图像景物的类别。均值计算如公式1所示,
Mean ( i ) = 1 n Σ ( x , y ) ∈ i I ( x , y )          公式
公式中i代表像素块编号,n代表块内像素的总数目,I(x,y)是像素灰度值。
像素块的方差表征了块内的像素变化特征,也能区分不同的图像区域。方差公式如公式2所示,
Std ( i ) = 1 n Σ ( x , y ) ∈ i ( I ( x , y ) - mean ) 2          公式
公式2中mean是像素块的灰度均值,其他符号的代表意义与公式1一致。
(2)对比度
对比度(contrast)反映了图像的清晰度和纹理沟纹深浅的程度。纹理沟纹越深,其对比度越大,视觉效果越清晰;反之,对比度小,则沟纹浅,效果模糊。计算方法如公式3所示。
CON = Σ i , j | i - j | 2 p ( i , j )         公式
公式中,p(i,j)是灰度共生矩阵中(i,j)位置的元素,i,j分别遍历整个矩阵。
(3)相关性
相关性(Correlation)度量空间灰度共生矩阵元素在行或列方向上的相似程度,因此,相关值大小反映了图像中局部灰度相关性。当矩阵元素值均匀相等时,相关值就大;相反,如果矩阵像元值相差很大则相关值小。如果图像中有水平方向纹理,则水平方向矩阵的COR大于其余方向矩阵的COR值。计算如公式4所示:
COR = Σ i , j ( i - μ i ) ( j - μ j ) p ( i , j ) σ i σ j        公式
公式4中,μi和μj分别代表灰度共生矩阵px(i)和py(j)的均值,σi和σj分别为px(i)和py(j)的标准差,px(i)和py(j)分别为灰度共生矩阵行和列方向的元素。
(4)能量
角二阶矩(ASM,angular second moment)是灰度共生矩阵元素值的平方和,也称能量,反映了图像灰度分布均匀程度和纹理粗细度。如果共生矩阵的所有值均相等,则ASM值小;相反,如果其中一些值大而其它值小,则ASM值大。当共生矩阵中元素集中分布时,此时ASM值大。ASM值大表明一种较均匀和规则变化的纹理模式。
ASM = Σ i , j p ( i , j ) 2                公式
上式中,p(i,j)是(i,j)处灰度共生矩阵的元素值。
(5)逆差距
逆差距(IDM,inverse different moment)又称同质性(Homogeneity),反映图像纹理的同质性,度量图像纹理局部变化的多少。其值大则说明图像纹理的不同区域间缺少变化,局部非常均匀。如果灰度共生矩阵对角元素有较大值,IDM就会取较大的值。因此连续灰度的图像会有较大IDM值。
IDM = Σ i , j p ( i , j ) 1 + | i - j |            公式
p(i,j)是(i,j)处灰度共生矩阵的元素值,i和j分别表示矩阵的不同元素位置。
(6)生成特征向量
综合以上各个特征就可以形成表征像素块的特征描述子,4个不同角度灰度共生矩阵的特征如公式7中fj=1,2,3,4(i)所示,最终形成特征描述子为F(i)。
F(i)=[Mean,Std,f1,f2,f3,f4]
                                       公式
fj=1,2,3,4(i)=[CON,COR,ASM,IDM]
公式中,i代表像素块编号。
步骤3:相似度匹配
本方法基于像素块的特征描述子的匹配,因此对于特征向量的计算至关重要。我们通过计算特征向量的欧氏距离表现像素块的近似程度,欧氏距离越小,表明特征描述差距越小,像素块越相似,反之越大,越不相似。两个n维向量an与bn的欧氏距离公式如公式8表示。
d 12 = Σ k = 1 n ( χ 1 k - χ 2 k ) 2            公式
对于两个像素块之间相似度匹配的计算不仅仅是欧氏距离的叠加,还要考虑灰度均值,方差和灰度共生矩阵特征的权重分配,才能更好的描述像素块纹理特征。如公式9所示,dStd和dMean分别表示像素块的灰度共生矩阵特征、方差和灰度均值的欧氏距离。ωj,ωStd和ωmean分别表示对应欧氏距离的权值,考虑到灰度共生矩阵对于纹理匹配的重要性,因此本文试验中对应权重分别取值0.2,0.1和0.1。
D ( i ) = Σ j = 1 4 ω j d f j ( i ) + ω std d Std + ω mean d Mean     公式
如公式10中,对于灰度共生矩阵中各个特征描述子的权重ηi取值分别取0.25,0.25,0.25,0.25。
d f j = 1,2,3,4 ( i ) = η 1 d CON + η 2 d COR + η 3 d ASM + η 4 d IDM     公式
假定{ri}代表彩色参考图像R的第i块像素块,{tj}代表灰度目标图像T的第j块像素块。如公式11所示,D(ri,tj)表示{ri}和{tj}之间加权特征向量的欧氏距离,则与{ri}欧氏距离最小的灰度目标像素块,就是与{ri}对应的最优匹配块。
a = ar g t j min { D ( r i , t j ) }                    公式
步骤4:图像颜色痕迹迁移
实现彩色参考图像和灰度目标图像的像素块匹配之后,就需要将彩色图像的颜色迁移到灰度图像中。在本方法中,我们选取像素块的中心窗口位置作为颜色迁移点,假设{ri}和{tj}是参考图像与目标图像对应的最优匹配块,设彩色参考像素块的尺寸大小为M*N,则选取的中心窗口长度范围是[M/4,3M/4],宽度范围[N/4,3N/4],灰度目标图像中心窗口的选取规则与之相同。然后计算彩色图像块中心窗口内色彩的颜色均值作为迁移色,然后同样选取灰度图像块中心窗口区域作为迁移区域,进行色彩赋值。
步骤5:颜色痕迹结果修正
由于灰度共生矩阵匹配可能产生误配的原因,因此接下来要根据相邻像素具有相似的灰度值则具有相似颜色值的准则,进行结果修正。主要过程是将计算灰度目标图像块的各中心窗口的灰度值,进行相邻九宫格区域的比对,如果发现灰度值相似,而迁移色不一致,则统计具有相似灰度值的相邻块的迁移色进行修正。如图3所示的具有不同灰度值的9块相邻像素块,图中编号5、6、9的位置与其他位置迁移色不同,因此就需要进行色彩修正。统计整个区域中相似灰度值像素块的色彩分布,类似于建立色彩直方图,选取出现次数最多的颜色作为修正色,然后将色彩不同的像素块颜色改正为修正色。图中5、9区域和除6之外的其他区域灰度值相似,因此颜色值要与其他区域保持一致,而6区域的灰度值不与其他区域灰度值一致,因此保持不变。
步骤6:最优化着色
首先,把颜色标记图像的颜色空间从RGB空间转换到YUV空间。接下来构造包含所有点颜色系数的稀疏矩阵,并根据相邻的像素,如果亮度相近则颜色相似的准则,构造任意像素与其邻域像素颜色系数方程组。为加快速度和提高效率,在本方法中采取以相邻的四方格或九宫格像素作为单个像素的方式进行着色。此单个像素记为四方格或九宫格像素块的超像素。计算四方格或九宫格像素块的亮度均值作为超像素的亮度值,此时原图像会缩减至超像素组成的新图像。新图像以超像素为彩色化的基本单元,用超像素依据本文中提到的最优化方法求解。求解到的超像素的色度值作为四方格或九宫格所有像素的色度值。
设任意像素p在YUV空间的色度值是U(p),则其与相邻像素q颜色的加权平均值的之间的差别计算公式如下,
J ( U ) = Σ p ( U ( p ) - Σ q ∈ N ( p ) ω pq U ( q ) ) 2         公式
公式中的ωpq是权函数,当相邻像素p与像素q的亮度值越相近时,权值越大。N(p)表示p的像素相邻域,U表示色度值。权函数ωpq可通过下列公式计算,
ω pq = e - ( Y ( p ) - Y ( q ) ) 2 / 2 σ 2         公式
公式中Y(p),Y(q)是像素的亮度,σ2是p邻域内亮度的方差。按上述公式构造关于图像中所有像素点颜色系数的超大稀疏矩阵,然后求解使J(U)最小的颜色系数U,即图像的颜色。求解方法通过迭代法来完成,最后得到结果图像。
以上所述,仅为本发明部分具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种基于GLCM的灰度图像彩色化系统,其特征在于:包括图像分块模块、灰度共生矩阵特征提取模块、相似度匹配建立对应关系模块、颜色标记和初步修正模块、最优化着色模块;
所述的图像分块模块,采用像素块网格划分法,根据基本像素块进行图像划分;
所述的灰度共生矩阵特征提取模块,针对每个像素块,分别计算像素块的灰度均值、方差和0°方向、45°方向、90°方向、135°方向的灰度共生矩阵,提取灰度共生矩阵的对比度、相关性、能量、同质性作为特征描述子,组成18维特征向量描述图像块的纹理特征,综合所述各个特征形成表征像素块的特征描述子;
所述的相似度匹配建立对应关系模块,通过计算特征向量的欧氏距离表现像素块的近似程度;欧氏距离最小的灰度目标像素块,即为与之对应的最优匹配块;
所述的颜色标记和初步修正模块,选取像素块的中心窗口位置作为颜色迁移点,然后计算彩色图像块中心窗口内色彩的颜色均值作为迁移色,然后同样选取灰度图像块中心窗口区域作为迁移区域,进行色彩赋值,生成标记涂鸦图像;
所述的最优化着色模块,把颜色痕迹图像的颜色空间从RGB空间转换到YUV空间,构造包含所有像素颜色系数的稀疏矩阵,依据任意像素与其邻域像素颜色加权差的平方和最小原则构建方程组,此方程组用迭代法求解。
2.根据权利要求1所述的基于GLCM的灰度图像彩色化系统,其特征在于:所述的基本像素块数目准则是10*10,对于不能整分的图像,采取向下取整的方法。
3.根据权利要求1所述的基于GLCM的灰度图像彩色化系统,其特征在于:所述的两个像素块之间相似度匹配的计算除了欧氏距离的叠加外,还要进行灰度均值、方差和灰度共生矩阵特征的权重分配;所述的像素块的灰度共生矩阵特征、方差和灰度均值的对应权重分别取值0.2,0.1和0.1,对于灰度共生矩阵中各个特征描述子的权重取值分别取0.25,0.25,0.25,0.25。
4.根据权利要求1所述的基于GLCM的灰度图像彩色化系统,其特征在于:
像素块的中心窗口和灰度目标图像中心窗口的选取规则为:假设{ri}和{tj}是参考图像与目标图像对应的最优匹配块,设彩色参考像素块的尺寸大小为M*N,则选取的中心窗口长度范围是[M/4,3M/4],宽度范围[N/4,3N/4]。
5.根据权利要求1所述的基于GLCM的灰度图像彩色化系统,其特征在于:所述的灰度共生矩阵匹配产生误配时,根据相邻像素具有相似的灰度值则具有相似颜色值的准则,进行结果修正,所述的修正过程是将计算灰度目标图像块的各中心窗口的灰度值,进行相邻九宫格区域的比对,如果发现灰度值相似,而迁移色不一致,则统计具有相似灰度值的相邻块的迁移色进行修正。
6.根据权利要求1所述的基于GLCM的灰度图像彩色化系统,其特征在于:所述的图像将以相邻的四方格或九宫格像素作为单个像素的方式进行着色,此单个像素记为四方格或九宫格像素块的超像素。计算四方格或九宫格像素块的亮度均值作为超像素的亮度值,此时原图像会缩减至超像素组成的新图像,新图像以超像素为彩色化的基本单元,将超像素代入最优化着色方法求解,求解得到的超像素的色度值作为四方格或九宫格所有像素的色度值。
7.一种基于GLCM的灰度图像彩色化系统的方法,其特征在于实现步骤如下:
步骤1:图像分块:
即采用像素块网格划分法,根据基本像素块对图像进行像素块划分,然后利用像素块进行特征提取和匹配;
步骤2:灰度共生矩阵特征提取:
即针对每个像素块,分别计算像素块的灰度均值、方差和0°方向、45°方向、90°方向、135°方向的灰度共生矩阵,提取灰度共生矩阵的对比度、相关性、能量、同质性作为特征描述子,一共组成18维特征向量描述图像块的纹理特征,综合所述各个特征形成表征像素块的特征描述子;
步骤3:相似度匹配建立对应关系:
即通过计算特征向量的欧氏距离表现像素块的近似程度,计算出欧氏距离最小的灰度目标像素块,就是与之对应的最优匹配块;
步骤4:颜色标记和初步修正:
即选取像素块的中心窗口位置作为颜色迁移点,然后计算彩色图像块中心窗口内色彩的颜色均值作为迁移色,然后同样选取灰度图像块中心窗口区域作为迁移区域,进行色彩赋值,生成标记涂鸦图像;
步骤5:最优化着色:
依据若相邻像素具有相同的亮度值则应具有相同的颜色值的原理。先把颜色痕迹图像的颜色空间从RGB空间转换到YUV空间,构造包含所有超像素颜色系数的稀疏矩阵,依据任意像素与其邻域像素颜色加权差的平方和最小原则构建方程组,此方程组用迭代法求解。
8.根据权利要求7所述的方法,其特征在于:两个像素块之间相似度匹配的计算除了欧氏距离的叠加,还要进行灰度均值,方差和灰度共生矩阵特征的权重分配。像素块的灰度共生矩阵特征、方差和灰度均值的对应权重分别取值0.2,0.1和0.1;对于灰度共生矩阵中各个特征描述子的权重取值分别取0.25,0.25,0.25,0.25。
9.根据权利要求7所述的方法,其特征在于:所述的像素块的中心窗口和灰度目标图像中心窗口的选取规则为:假设{ri}和{tj}是参考图像与目标图像对应的最优匹配块,设彩色参考像素块的尺寸大小为M*N,则选取的中心窗口长度范围是[M/4,3M/4],宽度范围[N/4,3N/4],灰度目标图像中心窗口的选取规则与之相同。
10.根据权利要求7所述的方法,其特征在于:所述的图像将以相邻的四方格或九宫格像素作为单个像素的方式进行着色,此单个像素记为四方格或九宫格像素块的超像素;计算四方格或九宫格像素块的亮度均值作为超像素的亮度值,以超像素为彩色化的基本单元,代入最优化着色方法求解,求解得到的超像素的色度值作为四方格或九宫格像素块中所有像素的色度值。
CN201410685330.9A 2014-11-25 2014-11-25 一种基于glcm的灰度图像彩色化系统和方法 Active CN104376529B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410685330.9A CN104376529B (zh) 2014-11-25 2014-11-25 一种基于glcm的灰度图像彩色化系统和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410685330.9A CN104376529B (zh) 2014-11-25 2014-11-25 一种基于glcm的灰度图像彩色化系统和方法

Publications (2)

Publication Number Publication Date
CN104376529A true CN104376529A (zh) 2015-02-25
CN104376529B CN104376529B (zh) 2017-08-11

Family

ID=52555420

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410685330.9A Active CN104376529B (zh) 2014-11-25 2014-11-25 一种基于glcm的灰度图像彩色化系统和方法

Country Status (1)

Country Link
CN (1) CN104376529B (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104851074A (zh) * 2015-03-26 2015-08-19 温州大学 基于特征相似性的非局部邻域灰度图像彩色化方法
CN105761202A (zh) * 2016-02-03 2016-07-13 武汉大学 一种彩色图像颜色迁移方法
CN106780583A (zh) * 2016-11-28 2017-05-31 国家海洋局第二海洋研究所 基于矩阵元素插值的任意灰度共生矩阵的计算方法及应用
CN106886791A (zh) * 2016-12-28 2017-06-23 四川木牛流马智能科技有限公司 一种基于条件随机场的二维ct图片中脂肪位置识别方法
CN106934846A (zh) * 2015-12-29 2017-07-07 深圳先进技术研究院 一种布料图像处理方法及系统
CN107018410A (zh) * 2017-05-10 2017-08-04 北京理工大学 一种基于预注意机制和空间依赖性的无参考图像质量评价方法
TWI596572B (zh) * 2016-07-06 2017-08-21 Method of automatically coloring image blocks
CN107146246A (zh) * 2017-05-08 2017-09-08 湘潭大学 一种用于工件加工表面背景纹理抑制方法
CN107705268A (zh) * 2017-10-20 2018-02-16 天津工业大学 一种基于改进的Retinex与Welsh近红外图像增强与彩色化算法
CN107730568A (zh) * 2017-10-31 2018-02-23 山东师范大学 基于权重学习的着色方法和装置
CN107856466A (zh) * 2017-11-10 2018-03-30 扬州市维拉园艺有限公司 一种树皮纹花盆模具的设计方法、模具及花盆
CN109903348A (zh) * 2018-05-14 2019-06-18 西安科技大学 基于分块特征交互图像色彩编辑方法、数字图像处理系统
CN110880003A (zh) * 2019-10-12 2020-03-13 中国第一汽车股份有限公司 一种图像匹配方法、装置、存储介质及汽车
CN112073596A (zh) * 2020-09-18 2020-12-11 青岛大学 特定黑白视频信号的仿真彩色处理方法及系统
CN112950461A (zh) * 2021-03-27 2021-06-11 刘文平 全局与超像素分割相融合的色彩迁移方法
CN114596372A (zh) * 2022-05-07 2022-06-07 武汉天际航信息科技股份有限公司 一种影像颜色迁移方法、影像一致性改善方法和装置
CN115018953A (zh) * 2022-08-04 2022-09-06 广东时谛智能科技有限公司 鞋体颜色和纹理的确定方法、装置、电子设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102779351A (zh) * 2012-06-08 2012-11-14 温州大学 基于局部线性模型优化的交互式灰度图像彩色化方法
CN103489161A (zh) * 2013-09-12 2014-01-01 南京邮电大学 一种灰度图像彩色化方法及装置
US20140071347A1 (en) * 2012-09-07 2014-03-13 Beihang University Method for Editing Propagation of Video and Image Content Based on Local Feature Structure Preservation
CN103839079A (zh) * 2014-03-18 2014-06-04 浙江师范大学 一种基于分类学习的相似图像彩色化算法
CN103985112A (zh) * 2014-03-05 2014-08-13 西安电子科技大学 基于改进多目标粒子群优化聚类的图像分割方法
CN104103082A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种基于区域描述和先验知识的图像显著性检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102779351A (zh) * 2012-06-08 2012-11-14 温州大学 基于局部线性模型优化的交互式灰度图像彩色化方法
US20140071347A1 (en) * 2012-09-07 2014-03-13 Beihang University Method for Editing Propagation of Video and Image Content Based on Local Feature Structure Preservation
CN103489161A (zh) * 2013-09-12 2014-01-01 南京邮电大学 一种灰度图像彩色化方法及装置
CN103985112A (zh) * 2014-03-05 2014-08-13 西安电子科技大学 基于改进多目标粒子群优化聚类的图像分割方法
CN103839079A (zh) * 2014-03-18 2014-06-04 浙江师范大学 一种基于分类学习的相似图像彩色化算法
CN104103082A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种基于区域描述和先验知识的图像显著性检测方法

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104851074A (zh) * 2015-03-26 2015-08-19 温州大学 基于特征相似性的非局部邻域灰度图像彩色化方法
CN104851074B (zh) * 2015-03-26 2017-12-19 温州大学 基于特征相似性的非局部邻域灰度图像彩色化方法
CN106934846B (zh) * 2015-12-29 2020-05-22 深圳先进技术研究院 一种布料图像处理方法及系统
CN106934846A (zh) * 2015-12-29 2017-07-07 深圳先进技术研究院 一种布料图像处理方法及系统
CN105761202A (zh) * 2016-02-03 2016-07-13 武汉大学 一种彩色图像颜色迁移方法
CN105761202B (zh) * 2016-02-03 2018-10-26 武汉大学 一种彩色图像颜色迁移方法
TWI596572B (zh) * 2016-07-06 2017-08-21 Method of automatically coloring image blocks
CN106780583B (zh) * 2016-11-28 2019-11-19 自然资源部第二海洋研究所 基于矩阵元素插值的任意灰度共生矩阵的计算方法及应用
CN106780583A (zh) * 2016-11-28 2017-05-31 国家海洋局第二海洋研究所 基于矩阵元素插值的任意灰度共生矩阵的计算方法及应用
CN106886791A (zh) * 2016-12-28 2017-06-23 四川木牛流马智能科技有限公司 一种基于条件随机场的二维ct图片中脂肪位置识别方法
CN107146246A (zh) * 2017-05-08 2017-09-08 湘潭大学 一种用于工件加工表面背景纹理抑制方法
CN107018410A (zh) * 2017-05-10 2017-08-04 北京理工大学 一种基于预注意机制和空间依赖性的无参考图像质量评价方法
CN107705268A (zh) * 2017-10-20 2018-02-16 天津工业大学 一种基于改进的Retinex与Welsh近红外图像增强与彩色化算法
CN107730568B (zh) * 2017-10-31 2021-01-08 山东师范大学 基于权重学习的着色方法和装置
CN107730568A (zh) * 2017-10-31 2018-02-23 山东师范大学 基于权重学习的着色方法和装置
CN107856466A (zh) * 2017-11-10 2018-03-30 扬州市维拉园艺有限公司 一种树皮纹花盆模具的设计方法、模具及花盆
CN109903348A (zh) * 2018-05-14 2019-06-18 西安科技大学 基于分块特征交互图像色彩编辑方法、数字图像处理系统
CN109903348B (zh) * 2018-05-14 2023-11-17 秦皇岛知聚科技有限公司 基于分块特征交互图像色彩编辑方法、数字图像处理系统
CN110880003A (zh) * 2019-10-12 2020-03-13 中国第一汽车股份有限公司 一种图像匹配方法、装置、存储介质及汽车
CN112073596B (zh) * 2020-09-18 2021-07-20 青岛大学 特定黑白视频信号的仿真彩色处理方法及系统
CN112073596A (zh) * 2020-09-18 2020-12-11 青岛大学 特定黑白视频信号的仿真彩色处理方法及系统
CN112950461A (zh) * 2021-03-27 2021-06-11 刘文平 全局与超像素分割相融合的色彩迁移方法
CN114596372A (zh) * 2022-05-07 2022-06-07 武汉天际航信息科技股份有限公司 一种影像颜色迁移方法、影像一致性改善方法和装置
CN114596372B (zh) * 2022-05-07 2022-07-29 武汉天际航信息科技股份有限公司 一种影像颜色迁移方法、影像一致性改善方法和装置
CN115018953A (zh) * 2022-08-04 2022-09-06 广东时谛智能科技有限公司 鞋体颜色和纹理的确定方法、装置、电子设备及存储介质
CN115018953B (zh) * 2022-08-04 2023-02-10 广东时谛智能科技有限公司 鞋体颜色和纹理的确定方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN104376529B (zh) 2017-08-11

Similar Documents

Publication Publication Date Title
CN104376529A (zh) 一种基于glcm的灰度图像彩色化系统和方法
CN106326937B (zh) 基于卷积神经网络的人群密度分布估计方法
CN107767413A (zh) 一种基于卷积神经网络的图像深度估计方法
CN107862668A (zh) 一种基于gnn的文物图像复原方法
CN106952271A (zh) 一种基于超像素分割和em/mpm处理的图像分割方法
CN103617596A (zh) 一种基于流型转换的图像色彩风格变换方法
CN104182949A (zh) 基于直方图特征点配准的影像匀色与融合方法及系统
CN110675462A (zh) 一种基于卷积神经网络的灰度图像彩色化方法
CN108305266A (zh) 基于条件随机场图结构学习的语义图像分割方法
CN105761202B (zh) 一种彩色图像颜色迁移方法
CN104217440B (zh) 一种从遥感图像中提取建成区的方法
CN104809245A (zh) 一种图像检索方法
CN108510562A (zh) 基于图像分形纹理的数码迷彩图案生成方法
CN103854261A (zh) 色偏图像的校正方法
CN107481183A (zh) 一种基于多维数据关联规则的微光图像彩色化方法
CN110176053A (zh) 一种大规模实景三维整体匀色方法
CN112766102A (zh) 一种基于空谱特征融合的无监督高光谱视频目标跟踪方法
CN112529774A (zh) 一种基于CycleGAN的遥感仿真图像生成方法
CN102542528B (zh) 一种图像转换处理方法及系统
Li et al. Global color consistency correction for large-scale images in 3-D reconstruction
CN102722862B (zh) 采用优化技术的单张图片半自动二维转三维方法和装置
CN112561844B (zh) 一种融合纹理结构的数码迷彩图案自动生成方法
CN104090940B (zh) 一种时序网络与时序数据的多态聚类方法
CN110599416B (zh) 一种基于空间目标图像数据库的非合作目标图像盲复原方法
CN105139331A (zh) 一种面向比特的基于仿生元胞自动机的图像置乱方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20160201

Address after: 518000, room 30, building 3001, Dayun software Town, Longgang District, Shenzhen, Guangdong

Applicant after: Shenzhen Beihang Tianhui Business Incubator Co.,Ltd.

Address before: 518000 room B407, virtual university garden building, South District, Nanshan District hi tech Zone, Guangdong, Shenzhen

Applicant before: SHENZHEN BEIHANG EMERGING INDUSTRIAL TECHNOLOGY Research Institute

GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 518000, room 30, building 3001, Dayun software Town, Longgang District, Shenzhen, Guangdong

Patentee after: Shenzhen Tianhui Business Incubator Co.,Ltd.

Address before: 518000, room 30, building 3001, Dayun software Town, Longgang District, Shenzhen, Guangdong

Patentee before: Shenzhen Beihang Tianhui Business Incubator Co.,Ltd.

CP01 Change in the name or title of a patent holder
TR01 Transfer of patent right

Effective date of registration: 20210419

Address after: 355 200 No. 22 Pengjialan, Zhuxia Village, Tai Lao Shan Town, Fuding City, Ningde City, Fujian Province

Patentee after: Lin Zefeng

Address before: 518000, room 30, building 3001, Dayun software Town, Longgang District, Shenzhen, Guangdong

Patentee before: Shenzhen Tianhui Business Incubator Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210706

Address after: 21f-a, Shenmao business center, 59 Xinwen Road, Meiling community, Lianhua street, Futian District, Shenzhen, Guangdong 518000

Patentee after: Shenzhen Hanfei Network Technology Co.,Ltd.

Address before: No.22, pengjia'ao, Zhuxia village, Taimushan Town, Fuding City, Ningde City, Fujian Province

Patentee before: Lin Zefeng

TR01 Transfer of patent right