一种基于像素值量化和小波分析的全自动照片油画化方法技术领域
本方法属于数字图像/视频信号处理领域，尤其涉及一种基于像素值量化和小波分析的全
自动照片油画化方法。
背景技术
随着生活水平的不断提高，如今人们不仅关注衣食住行等物质层面上的基本需要，而且
越来越注重精神层面上的艺术享受。油画作为一种历史悠久、表现力强的艺术形式，受到了
广大艺术爱好者的青睐。然而，绘制油画需要具备专业绘画技能，不具备专业绘画技能的人
无法自己绘制油画。近年来，照片油画化技术的出现突破了专业绘画技能的限制。该技术以
真实照片为基础，采用图像处理等手段使照片具有类似于油画的显示风格。照片油画化技术
使不具备任何专业绘画技能的人也能制作出极具艺术气息的油画，在很大程度上满足了人们
的艺术创作愿望。
申请号为CN200910148465.0的专利公开了一种通过图像制备油画的方法。该方法将图像
按照不同颜色进行分块，然后采用人工或软件处理方式修正分块线框图，最后将该线框图打
印并进行人工上色。该方法虽然降低了油画的绘画难度，但仍需要制作者具备一定的绘画技
巧，整体的油画制作工序也较为复杂。
申请号为CN201010142541.X的专利公开了一种图像的油画风格化方法。该方法首先根
据能量值将图像划分成面积不等的区域，再根据区域的面积在区域内布置相应尺寸的笔划，
然后根据原始图像的笔划方向场来渲染布置的笔划，最后为渲染结果添加光照效果。该方法
虽然能够避免人工上色过程，但是制作过程中的阈值选择、笔划刷设定都依赖人工来完成，
参数的设置过程较为繁琐，对油画制作者的要求较高。
申请号为CN201110385358.7的专利公开了一种基于图像处理的照片转化为数字油画的
方法。该方法通过滤波、取样、聚类过程将图像中的颜色转化为数字油画允许的颜色，从而
完成油画化处理。但是，聚类中心需要人工选取。而且聚类结果很可能不理想，此时甚至需
要反复选取聚类中心，来获得符合条件的结果。该技术方案的油画制作过程仍然比较复杂。
申请号为CN201410616110.0的专利公开了一种基于图像的数字油画画布的自动生成方
法。该方法虽然能够自动地划分出不同的油画颜色区域，但制作油画时仍然依赖人工上色，
不能自动生成油画。
现有的照片油画化技术存在的不足可以大致总结为以下几点：
1.只能根据照片生成相应的填涂框架，仍需要进行后续的人工上色才能完成油画制作；
2.需要手工设定的参数较多。为了获得较好的效果，油画制作者首先需要对各步骤中的
参数逐个了解才能进行具体操作；
3.油画生成过程步骤繁复，油画制作耗时较长。
发明内容
针对上述技术存在的问题，本发明提出了一种基于像素值量化和小波分析的全自动照片
油画化方法。该方法利用小波分析技术将数字照片的内容划分为平滑区域和纹理区域，然后
分别采用不同的量化步长对平滑区域和纹理区域分别进行像素值量化，最后在小波域中进行
图像纹理融合，从而得到具有油画风格的图像。
为实现上述目的，本发明采用以下技术方案：
本发明提出的基于像素值量化和小波分析的全自动照片油画化方法的总体流程图如图1
所示，具体包括以下步骤：
步骤1，计算均匀度。
步骤1.1，图像灰度化。将原始的彩色照片图像转化为灰度图像。
步骤1.2，获得灰度直方图。统计步骤1.1中得到的灰度图像的灰度直方图。
步骤1.3，计算均匀度A。统计直方图中频数超过总像素数0.7％的灰度级数目，该数值
即定义为均匀度A。均匀度A能够近似表征直方图的分布均匀程度。灰度直方图分布越均匀，
则均匀度A的数值越大。
步骤1.4，均匀度A下限控制。若步骤1.3计算出的均匀度A小于25，则将均匀度A设
定为25；若步骤1.3计算出的均匀度A大于或等于25，则跳过此步骤。
步骤2，获得细节图。
步骤2.1，基于小波分析的边缘检测。对步骤1.1中得到的灰度图像进行小波变换，然后
将小波域中的低频子带系数全部置零，最后进行小波反变换，得到边缘检测结果。
步骤2.2，图像二值化。将步骤2.1中得到的边缘检测结果进行二值化处理，二值化阈值
取0.4倍的A。
步骤2.3，中值滤波。对步骤2.2中得到的二值化图像进行中值滤波去噪，中值滤波窗口
尺寸设置为5×5。
步骤2.4，形态学膨胀。对步骤2.3中得到的滤波结果进行形态学膨胀运算，获得细节图，
膨胀运算半径设置为3。
步骤3，像素值量化。
步骤3.1，区域划分。根据步骤2.4中得到的细节图对原始的彩色照片图像进行区域划分。
细节图中为白色的部分对应原始彩色照片图像的纹理区域；细节图中为黑色的部分对应原始
彩色照片图像的平滑区域。
步骤3.2，纹理区域像素值量化。对原始彩色照片图像中的纹理区域的RGB三个颜色通
道的像素值进行量化，量化步长选为A。
步骤3.3，平滑区域像素值量化。对原始彩色照片图像中的平滑区域的RGB三个颜色通
道的像素值进行量化，量化步长选为3倍的A。
步骤3.4，自动合成。完成步骤3.2和步骤3.3后，纹理区域和平滑区域的像素值都已经
过量化处理，它们自动合成得到完整的原始照片图像量化结果。
步骤4，纹理融合。
步骤4.1，画布纹理提取。对没有任何绘画内容的油画画布图像进行小波变换，然后将低
频子带系数置零，只保留细节子带系数。
步骤4.2，量化图像的小波变换。对步骤3.4中得到的量化后图像进行小波变换，得到各
个子带的小波系数。
步骤4.3，小波域纹理融合。将步骤4.1中得到的小波系数和步骤4.2中得到的小波系数
对应相加，然后进行小波反变换，得到最终的照片油画化结果。
与现有技术相比，本发明具有以下明显的优势和有益的效果。
1.照片油画化的所有步骤均由计算机算法实现，无需额外的人工处理步骤；
2.算法实现的过程是全自动的，在算法执行过程中不需要人工设置任何技术参数，仅需
给定输入的照片图像即可获得结果，操作简单；
3.计算量较小，算法运行速度快。
附图说明
图1为本发明所涉及的基于像素值量化和小波分析的全自动照片油画化方法的总体流程
图；
图2为步骤1的详细流程图；
图3为步骤2的详细流程图；
图4为步骤3的详细流程图；
图5为步骤4的详细流程图；
图6为三张油画化处理前的原始照片图像；
图7为采用本发明的技术方案得到的三张照片油画化处理结果。
具体实施方式
以下结合附图和具体实施方式对本发明做进一步说明。
如图1所示，本发明所涉及的基于像素值量化和小波分析的全自动照片油画化方法包含
4个主要步骤：计算均匀度、获得细节图、像素值量化和纹理融合。接下来对每一个步骤的
具体实施方式进行详细介绍：
1.计算均匀度
该步骤的详细流程如图2所示。首先，将原始的照片图像转化为灰度图像，接着得到灰
度图像的灰度直方图，然后统计直方图中频数超过总像素数0.7％的灰度级数目，该数目即定
义为均匀度A。均匀度A可以近似表征直方图的分布均匀程度。灰度直方图分布越均匀，则
均匀度A的数值越大。均匀度A将作为后续计算步骤的重要参考指标。为保证后续步骤的有
效进行，设置均匀度A下限值为25。若按照上述方法计算出的均匀度A不足25，则均匀度A
按25计算。
2.获得细节图
该步骤的详细流程如图3所示。将步骤1中得到的灰度图像利用小波变换进行边缘检测，
然后对边缘检测结果进行二值化处理，二值化阈值取0.4A，之后对二值化图像进行中值滤波
去噪，最后进行形态学膨胀运算得到最终的细节图。细节图中，白色部分代表纹理区域，黑
色部分代表平滑区域。
3.像素值量化
该步骤的详细流程如图4所示。步骤2中得到的细节图能够将原始照片图像划分为纹理
区域与平滑区域。对于纹理区域，采用大步长进行量化，用来模拟油画局部的粗糙感；对于
平滑区域，采用小步长进行量化，在保持概貌基本不变的情况下增加颜色的分层效果，提升
层次感。本实施例取纹理区域的量化步长为3A，取平滑区域的量化步长为A。
4.纹理融合
该步骤的详细流程如图5所示。此步骤旨在使处理结果表现出类似于油画画布的粗糙质
感。本实施例选取了一幅空白的油画画布图像，然后对画布图像进行小波变换，并且将小波
的低频系数置零。之后对步骤3得到的照片图像量化结果进行小波变换，并把画布的小波系
数与该变换结果进行叠加。最后进行小波反变换，从而得到照片油画化的处理结果。
图6所示的是油画化处理前的原始照片图像。
图7所示的是采用本发明所述的技术方案得到的照片油画化处理结果。
本说明书实施例描述的内容仅仅是对发明构思的实现形式的举例说明，本发明的保护范
围不应当被视为仅限于实施例所述的具体形式。本发明的保护范围也及于本领域技术人员根
据本发明构思所能够想到的等同技术手段。