基于单张目标画像的人脸画像合成方法

摘要

本发明公开了一种基于单张目标画像的人脸画像合成方法，在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，将所述第二图像块子集投影到投影矩阵得到第一投影系数集合；然后将测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到匹配的图像块，将查找到的图像块合成得到初始合成图像；最后进行迭代合成最终合成图像。本发明省去费时费力的收集大量照片与画像对的过程，方便合成各种风格的画像。

权利要求

1.一种基于单张目标画像的人脸画像合成方法，其特征在于，该合成方法为：在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，将所述第二图像块子集投影到投影矩阵得到第一投影系数集合；然后将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，将查找到的图像块合成得到初始合成图像；将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像。

2.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，具体为：将CUFS数据库中所有照片与画像对中的图像由彩色图像变成灰度图像，再从照片与画像对集合中选择一对照片与画像中的画像作为目标画像，再将剩下的照片与画像对中的所有照片作为测试照片样本集 ，并从测试照片样本集 中选取一张测试照片 。

3.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，具体为：

3a) 对目标画像建立 层高斯金字塔，高斯金字塔中的每一层都对应着某个尺度下的画像对：

3a1) 构造尺度值为 的高斯函数：

其中， 表示在 尺度值下的高斯函数，，分别表示图像中像素点对应的横、纵坐标值；

3a2) 将图像与3a1)所得的高斯函数进行卷积，得到高斯金字塔第一层的图像；

3a3) 设 ；

3a4) 将尺度值 乘以一个比例系数 作为新的尺度因子来卷积原图像，再将卷积后的图像做比例因子为 的降采样，得到高斯金字塔第 层的图像；

3a5) 设 ；

3a6) 重复步骤3a4)-3a5) 次，即可得到 层高斯金字塔；

3b) 将高斯金字塔中奇数层对应的画像设为字典训练样本集 ，将高斯金字塔中偶数层对应的画像设为合成训练样本集 ;

3c) 将字典训练样本集 和合成训练样本集 中所有的画像划分为相同大小及相同重叠程度的块，从而分别得到字典训练画像块集合 即第一图像块子集和合成训练画像块集合 即第二图像块子集。

4.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，具体为：

4a) 对字典训练画像块集合 即第一图像块子集利用有效稀疏编码方法得到训练画像块特征字典 ：

4a1) 初始化训练画像块特征字典 为高斯矩阵，并将高斯矩阵每一个列向量归一化为单位向量；其中，表示训练画像块特征字典 中字典元素的维数，表示训练画像块特征字典 中字典元素的个数；

4a2) 给定字典训练画像块集合 以及初始化训练画像块特征字典 ，按照下式得到字典训练画像块集合 的初始稀疏表示集合 ：其中，是正则化参数，设值为0.15， 表示2范数， 表示1范数，当给定 时，可以通过特征符号搜索算法进行求解；

4a3) 给定字典训练画像块集合 以及由步骤4a2)得到的初始稀疏表示集合 ，按照下式得到估计的训练画像块特征字典 ：,

当给定 时，通过拉格朗日对偶方法求解出 ；

4a4) 交替迭代步骤4a2)和步骤4a3)，直到以下目标函数收敛，得到最终的训练画像块特征字典 ：,

。

5.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将所述第二图像块子集投影到投影矩阵得到第一投影系数集合，具体为：将合成训练样本集中所有的画像划分为相同大小及相同重叠程度的块，从而得到合成训练画像块集合即第二图像块子集；

5a) 给定训练画像块 和训练画像块特征字典 ；

5b) 构造以下最优化问题来求解训练画像块 对应的稀疏表示 ：

，从而得到合成训练画像块集合 对应的稀疏表示集合 ，其

中，是正则化参数，设值为0.15，上述最优化问题可利用最小角回归方法进行求解；

5c) 对合成训练画像块集合 即第二图像块子集中的每个训练画像块 重复步骤5a)和5b)，得到合成训练画像块集合 对应的稀疏表示集合即第一投影系数集合：；

其中，表示第 个画像块 的稀疏表示，由于每个稀疏表示 都包含稀疏系数值 以及稀疏系数值在稀疏编码中被求解的顺序 ，因此稀疏表示集合 包含有稀疏系数值集合以及对应的稀疏系数值在稀疏编码中被求解的顺序集合 ，其表示如下：；

。

6.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，具体为：

6a) 将测试照片 划分为相同大小及相同重叠程度的测试照片块集合： ，表示测试照片 中的第个照片块；

6b) 将每个测试照片块 用训练画像块特征字典 进行稀疏编码，得到 对应的稀疏表示 即第二投影系数集合；

其中，表示第 个测试照片块 的稀疏表示，每个稀疏表示 都包含稀疏系数值 以及稀疏系数值在稀疏编码中被求解的顺序 。

7.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，具体为：对合成训练画像块集合 通过贪婪搜索获取每个测试照片块 对应的待选择画像块集合 ：

7a) 设 ， ；

7b) 从顺序集合 中选出与求解顺序 中第 个值相同的所有稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断 中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7c)；

7c) 从稀疏系数值集合 中选出与稀疏系数值 中第 个值欧式距离最相近的 稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7d)；

7d) 令 ，返回步骤7b)；

7e) 将稀疏表示集合 中所有 个稀疏表示对应的合成训练画像块集合 中的 个画像块作为稀疏表示 对应的测试照片块 的待选择画像块集合 ：如果，则利用复制操作将待选择画像块从 个增加到 个，如果 ，则利用测试照片块 的高频信息与待选择画像块的高频信息的欧式距离将待选择画像块从 个减少为个，使待选择画像块集合从 变成 。

8.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将查找到的图像块合成得到初始合成图像，具体为：根据求解马尔可夫随机场模型从待选择画像块集合 中得到每个测试照片块 最终对应的训练画像块，并将这些训练画像块进行融合得到与测试照片 对应的初始合成画像 ， 。

9.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，具体为：提取初始合成画像 、测试照片 以及目标画像的像素值特征以及SURF特征，提取初始合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，获取测试照片块 对应的待选择画像块集合 ，根据所述待选择画像块集合 获取新合成画像。

10.根据权利要求1所述的基于单张目标画像的人脸画像合成方法，其特征在于，所述将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像，具体为：新合成画像代替初始合成画像，迭代执行提取初始合成画像 、测试照片 以及目标画像的像素值特征以及SURF特征，提取初始合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，获取测试照片块 对应的待选择画像块集合 ，根据所述待选择画像块集合 获取新合成画像，直到迭代次数达到指定值 时停止，则新合成画像 ， ，为最终合成画像。

说明书

基于单张目标画像的人脸画像合成方法

技术领域

[0001] 本发明属于图像处理技术领域，具体涉及一种基于单张目标画像的人脸画像合成方法。

背景技术

[0002] 人脸画像合成有许多重要的实际应用，比如数字娱乐。目前推特和微博等社交网络变得越来越流行，很多用户将自己的画像上传到网上作为账户的头像以增加娱乐效果。为了体现个性化和增加趣味性，用户往往更希望自己的画像能有多种风格可供选择。此外，待合成照片的约束越少，用户就会更加亲睐，例如，待合成照片可以是各种角度，各种表情，各种背景以及全身照、合照等等。

[0003] 目前的人脸画像合成方法可以分成以下三大类：一，基于子空间学习的人脸画像合成方法；二，基于贝叶斯推断的人脸画像合成方法；三，基于稀疏表示的人脸画像合成方法。

[0004] Liu等人在文献“Liu Q, Tang X, Jin H, et al. A nonlinear approach for face sketch synthesis and recognition[C]//Computer Vision and Pattern Recognition, 2005. CVPR 2005. IEEE Computer Society Conference on. IEEE, 2005, 1: 1005-1010.”中采用子空间学习的方法来模拟照片转化成画像的非线性过程。具体实施方式是：先将训练集中的照片-画像对及测试照片划分成大小相同及重叠区域相同的图像块，对于测试照片中的每一个照片块在训练集中寻找它的 个近邻照片块并计算重建系数 ，然后再将 个照片块对应的 个画像块用重建系数 进行加权组合得到待合成画像块，最后将所有的待合成画像块融合得到合成画像。但是该方法存在以下缺点：各合成画像块都是单独合成，没有考虑画像块之间的邻域关系，导致融合处的画像不自然，存在噪声。

[0005] Wang等人在文献“Wang X, Tang X. Face photo-sketch synthesis and recognition[J]. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2009, 31(11): 1955-1967.”中提出了一种基于贝叶斯推断的人脸画像合成方法。具体操作步骤如下：首先将训练集中的照片-画像对以及测试照片用相同的方式分块，然后对每一个测试照片块在训练集中寻找它的 个候选照片块，根据测试照片块与 个候选照片块之间的关系以及合成画像中相邻位置的画像块之间的关系，建立马尔可夫随机场模型对每个测试照片块寻找一个最佳的训练画像块作为待合成画像块，最后将所有的待合成画像块融合得到最终的合成画像。该方法存在以下缺陷：每个测试照片块最终只选择一个训练画像块作为待合成画像块进行画像合成，导致合成画像存在块效应和细节缺失的问题。

[0006] 高新波等人申请的专利技术“基于稀疏表示的画像-照片生成方法”（申请号：201010289330.9申请日：2010-09-24申请公布号：CN 101958000 A）中公开了一种基于稀疏表示的人脸画像合成方法。具体实施方式是：首先采用已有的画像合成方法生成画像的初始估计，然后再利用稀疏表示的方法合成出合成画像的细节信息，最后将初始估计和细节信息进行线性叠加得到最后的合成画像。该方法存在以下缺陷：合成画像的好坏依赖于初始画像合成方法，此外细节信息的线性叠加使得合成画像不自然。

[0007] 上述三类画像合成方法需要许多照片-画像对作为训练集，然而收集照片-画像对费时费力。此外，当需要合成测试照片对应的某种类型的画像时，上述三类画像合成方法则需要收集许多对应类型的照片-画像对作为训练集，然而在大多数情况下，很难收集到许多指定类型的照片-画像对。此外，上述三类画像合成方法只适用于测试照片与训练集照片背景相同的情况，当测试照片与训练集照片背景不同时，合成画像的背景会存在噪声。当测试照片中存在非人脸部件，例如眼镜、发卡等，而训练集照片中都不存在这些非人脸部件时，合成画像会丢失测试照片中的非人脸部件。

发明内容

[0008] 有鉴于此，本发明的主要目的在于提供一种基于单张目标画像的人脸画像合成方法。

[0009] 为达到上述目的，本发明的技术方案是这样实现的：

[0010] 本发明实施例提供一种基于单张目标画像的人脸画像合成方法，该合成方法为：在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，将所述第二图像块子集投影到投影矩阵得到第一投影系数集合；然后将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，将查找到的图像块合成得到初始合成图像；将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像。

[0011] 上述方案中，所述在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，具体为：将CUFS数据库中所有照片与画像对中的图像由彩色图像变成灰度图像，再从照片与画像对集合中选择一对照片与画像中的画像作为目标画像，再将剩下的照片与画像对中的所有照片作为测试照片样本集 ，并从测试照片样本集 中选取一张测试照片 。

[0012] 上述方案中，所述将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，具体为：

[0013] 3a) 对目标画像建立 层高斯金字塔，高斯金字塔中的每一层都对应着某个尺度下的画像对：

[0014] 3a1) 构造尺度值为 的高斯函数：

[0015]

[0016] 其中， 表示在 尺度值下的高斯函数，，分别表示图像中像素点对应的横、纵坐标值；

[0017] 3a2) 将图像与3a1)所得的高斯函数进行卷积，得到高斯金字塔第一层的图像；

[0018] 3a3) 设 ；

[0019] 3a4) 将尺度值 乘以一个比例系数 作为新的尺度因子来卷积原图像，再将卷积后的图像做比例因子为 的降采样，得到高斯金字塔第 层的图像；

[0020] 3a5) 设 ；

[0021] 3a6) 重复步骤3a4)-3a5) 次，即可得到 层高斯金字塔；

[0022] 3b) 将高斯金字塔中奇数层对应的画像设为字典训练样本集 ，将高斯金字塔中偶数层对应的画像设为合成训练样本集 ;

[0023] 3c) 将字典训练样本集 和合成训练样本集 中所有的画像划分为相同大小及相同重叠程度的块，从而分别得到字典训练画像块集合 即第一图像块子集和合成训练画像块集合 即第二图像块子集。

[0024] 上述方案中，所述将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，具体为：

[0025] 4a) 对字典训练画像块集合 即第一图像块子集利用有效稀疏编码方法得到训练画像块特征字典 ：

[0026] 4a1) 初始化训练画像块特征字典 为高斯矩阵，并将高斯矩阵每一个列向量归一化为单位向量；其中，表示训练画像块特征字典 中字典元素的维数，表示训练画像块特征字典 中字典元素的个数；

[0027] 4a2) 给定字典训练画像块集合 以及初始化训练画像块特征字典 ，按照下式得到字典训练画像块集合 的初始稀疏表示集合 ：

[0028]

[0029] 其中，是正则化参数，设值为0.15， 表示2范数， 表示1范数，当给定 时，可以通过特征符号搜索算法进行求解；

[0030] 4a3) 给定字典训练画像块集合 以及由步骤4a2)得到的初始稀疏表示集合 ，按照下式得到估计的训练画像块特征字典 ：

[0031] ,

[0032]

[0033] 当给定 时，通过拉格朗日对偶方法求解出 ；

[0034] 4a4) 交替迭代步骤4a2)和步骤4a3)，直到以下目标函数收敛，得到最终的训练画像块特征字典 ：

[0035] ,

[0036] 。

[0037] 上述方案中，所述将所述第二图像块子集投影到投影矩阵得到第一投影系数集合，具体为：将合成训练样本集 中所有的画像划分为相同大小及相同重叠程度的块，从而得到合成训练画像块集合 即第二图像块子集；

[0038] 5a) 给定训练画像块 和训练画像块特征字典 ；

[0039] 5b) 构造以下最优化问题来求解训练画像块 对应的稀疏表示 ：

[0040] ，从而得到合成训练画像块集合 对应的稀疏表示集合 ，其中，是正则化参数，设值为0.15，上述最优化问题可利用最小角回归方法进行求解；

[0041] 5c) 对合成训练画像块集合 即第二图像块子集中的每个训练画像块 重复步骤5a)和5b)，得到合成训练画像块集合 对应的稀疏表示集合即第一投影系数集合：

[0042] ；

[0043] 其中，表示第 个画像块 的稀疏表示，由于每个稀疏表示 都包含稀疏系数值以及稀疏系数值在稀疏编码中被求解的顺序 ，因此稀疏表示集合 包含有稀疏系数值集合 以及对应的稀疏系数值在稀疏编码中被求解的顺序集合 ，其表示如下：

[0044] ；

[0045] 。

[0046] 上述方案中，所述将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，具体为：

[0047] 6a)  将测试照片 划分为相同大小及相同重叠程度的测试照片块集合：，表示测试照片 中的第个照片块；

[0048] 6b) 将每个测试照片块 用训练画像块特征字典 进行稀疏编码，得到 对应的稀疏表示 即第二投影系数集合；

[0049] 其中，表示第个测试照片块 的稀疏表示，每个稀疏表示 都包含稀疏系数值 以及稀疏系数值在稀疏编码中被求解的顺序 。

[0050] 上述方案中，所述根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，具体为：对合成训练画像块集合 通过贪婪搜索获取每个测试照片块 对应的待选择画像块集合 ：

[0051] 7a) 设 ， ；

[0052] 7b) 从顺序集合 中选出与求解顺序 中第个值相同的所有稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断 中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7c)；

[0053] 7c) 从稀疏系数值集合 中选出与稀疏系数值 中第个值欧式距离最相近的 稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7d)；

[0054] 7d) 令 ，返回步骤7b)；

[0055] 7e) 将稀疏表示集合 中所有 个稀疏表示对应的合成训练画像块集合 中的个画像块作为稀疏表示 对应的测试照片块 的待选择画像块集合 ：如果，则利用复制操作将待选择画像块从 个增加到 个，如果 ，则利用测试照片块的高频信息与待选择画像块的高频信息的欧式距离将待选择画像块从 个减少为 个，使待选择画像块集合从 变成 。

[0056] 上述方案中，所述将查找到的图像块合成得到初始合成图像，具体为：根据求解马尔可夫随机场模型从待选择画像块集合 中得到每个测试照片块 最终对应的训练画像块，并将这些训练画像块进行融合得到与测试照片 对应的初始合成画像 ， 。

[0057] 上述方案中，所述将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，具体为：提取初始合成画像 、测试照片 以及目标画像的像素值特征以及SURF特征，提取初始合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，获取测试照片块 对应的待选择画像块集合 ，根据所述待选择画像块集合 获取新合成画像。

[0058] 上述方案中，所述将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像，具体为：新合成画像代替初始合成画像，迭代执行提取初始合成画像 、测试照片以及目标画像的像素值特征以及SURF特征，提取初始合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，获取测试照片块 对应的待选择画像块集合 ，根据所述待选择画像块集合 获取新合成画像，直到迭代次数达到指定值 时停止，则新合成画像 ，，为最终合成画像。

[0059] 与现有技术相比，本发明的有益效果：

[0060] 第一，本发明由于只需要单张目标画像作为训练集，从而省去费时费力的收集大量照片与画像对的过程，因此能方便的合成各种风格的画像；

[0061] 第二，本发明由于采用全局搜索来寻找近邻，因此能很好的合成眼镜、发卡等非人脸部件；

[0062] 第三，本发明由于用稀疏表示中的稀疏系数值及稀疏系数值被求解的顺序这两个量来寻找近邻块，因此能有效的匹配拥有不同背景的图像块，从而有效解决测试照片与训练集照片背景不同的问题。

附图说明

[0063] 图1为本发明的实现流程图；

[0064] 图2为三种不同类型的画像；

[0065] 图3为本发明与现有的三类方法在CUHK student数据库、AR数据库和XM2VTS数据库上的合成画像的对比结果图；

[0066] 图4为本发明用三种不同类型的画像作为训练集时在CUHK student数据库、AR数据库和XM2VTS数据库上的合成画像的对比结果图。

[0067] 图5为本发明用三种不同类型的画像作为训练集时在网上下载的合照上的合成画像的对比结果图。

具体实施方式

[0068] 下面结合附图和具体实施方式对本发明进行详细说明。

[0069] 本发明实施例提供一种基于单张目标画像的人脸画像合成方法，图1所示，该合成方法为：在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，将所述第二图像块子集投影到投影矩阵得到第一投影系数集合；然后将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，将查找到的图像块合成得到初始合成图像；将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像。

[0070] 所述在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，具体为：将CUFS数据库中所有照片与画像对中的图像由彩色图像变成灰度图像，再从照片与画像对集合中选择一对照片与画像中的画像作为目标画像，再将剩下的照片与画像对中的所有照片作为测试照片样本集 ，并从测试照片样本集 中选取一张测试照片 。

[0071] 所述在数据库样本集内随机选择一对照片与画像对中的画像作为目标画像，其他剩余照片与画像对中的照片作为测试照片样本集，之后，该方法还包括：对目标画像的多尺度特征进行主成分分析来获取目标画像的多尺度特征的主成分。

[0072] 所述将对所述目标画像进行分块后得到的图像块进行划分获得第一图像块子集和第二图像块子集，具体为：

[0073] 3a) 对目标画像建立 层高斯金字塔，高斯金字塔中的每一层都对应着某个尺度下的画像对：

[0074] 3a1) 构造尺度值为 的高斯函数：

[0075]

[0076] 其中， 表示在 尺度值下的高斯函数，，分别表示图像中像素点对应的横、纵坐标值；

[0077] 3a2) 将图像与3a1)所得的高斯函数进行卷积，得到高斯金字塔第一层的图像；

[0078] 3a3) 设 ；

[0079] 3a4) 将尺度值 乘以一个比例系数 作为新的尺度因子来卷积原图像，再将卷积后的图像做比例因子为 的降采样，得到高斯金字塔第 层的图像；

[0080] 3a5) 设 ；

[0081] 3a6) 重复步骤3a4)-3a5) 次，即可得到 层高斯金字塔；

[0082] 3b) 将高斯金字塔中奇数层对应的画像设为字典训练样本集 ，将高斯金字塔中偶数层对应的画像设为合成训练样本集 ；

[0083] 3c) 将字典训练样本集 和合成训练样本集 中所有的画像划分为相同大小及相同重叠程度的块，从而分别得到字典训练画像块集合 即第一图像块子集和合成训练画像块集合 即第二图像块子集。

[0084] 所述将所述获得的第一图像块子集进行矩阵分解得到投影矩阵，具体为：

[0085] 4a) 对字典训练画像块集合 即第一图像块子集利用有效稀疏编码方法得到训练画像块特征字典 ：

[0086] 4a1) 初始化训练画像块特征字典 为高斯矩阵，并将高斯矩阵每一个列向量归一化为单位向量；其中，表示训练画像块特征字典 中字典元素的维数，表示训练画像块特征字典 中字典元素的个数；

[0087] 4a2) 给定字典训练画像块集合 以及初始化训练画像块特征字典 ，按照下式得到字典训练画像块集合 的初始稀疏表示集合 ：

[0088]

[0089] 其中，是正则化参数，设值为0.15， 表示2范数， 表示1范数，当给定 时，可以通过特征符号搜索算法进行求解；

[0090] 4a3) 给定字典训练画像块集合 以及由步骤4a2)得到的初始稀疏表示集合 ，按照下式得到估计的训练画像块特征字典 ：

[0091] ,

[0092]

[0093] 当给定 时，通过拉格朗日对偶方法求解出 ；

[0094] 4a4) 交替迭代步骤4a2)和步骤4a3)，直到以下目标函数收敛，得到最终的训练画像块特征字典 ：

[0095] ,

[0096] 。

[0097] 所述的有效稀疏编码方法，见文献“Lee H, Battle A, Raina R, et al. Efficient sparse coding algorithms[C]//Advances in neural information processing systems. 2006: 801-808.”，文献中包含了特征符号搜索算法和拉格朗日对偶方法。

[0098] 获取训练画像块特征字典 的方法不局限于上述的有效稀疏编码方法，还可用K-SVD算法，见文献“Aharon M, Elad M, Bruckstein A. k-svd: An algorithm for designing overcomplete dictionaries for sparse representation[J]. Signal Processing, IEEE Transactions on, 2006, 54(11): 4311-4322.”。

[0099] 所述将所述第二图像块子集投影到投影矩阵得到第一投影系数集合，具体为：将合成训练样本集 中所有的画像划分为相同大小及相同重叠程度的块，从而得到合成训练画像块集合 即第二图像块子集；

[0100] 5a) 给定训练画像块 和训练画像块特征字典 ；

[0101] 5b) 构造以下最优化问题来求解训练画像块 对应的稀疏表示 ：

[0102] ，从而得到合成训练画像块集合 对应的稀疏表示集合 ，其中，是正则化参数，设值为0.15，上述最优化问题可利用最小角回归方法进行求解。

[0103] 最小角回归方法见文献“Efron B, Hastie T, Johnstone I, et al. Least angle regression[J]. The Annals of statistics, 2004, 32(2): 407-499.”；

[0104] 5c) 对合成训练画像块集合 即第二图像块子集中的每个训练画像块 重复步骤5a)和5b)，得到合成训练画像块集合 对应的稀疏表示集合即第一投影系数集合：

[0105] 。

[0106] 其中，表示第 个画像块 的稀疏表示，由于每个稀疏表示 都包含稀疏系数值以及稀疏系数值在稀疏编码中被求解的顺序 ，因此稀疏表示集合 包含有稀疏系数值集合 以及对应的稀疏系数值在稀疏编码中被求解的顺序集合 ，其表示如下：

[0107] ；

[0108] 。

[0109] 所述将从测试照片样本集中选择的测试照片分块后得到的图像块投影到投影矩阵得到第二投影系数集合，具体为：

[0110] 6a)  将测试照片 划分为相同大小及相同重叠程度的测试照片块集合：，表示测试照片 中的第个照片块；

[0111] 6b) 将每个测试照片块 用训练画像块特征字典 进行稀疏编码，得到 对应的稀疏表示 即第二投影系数集合；

[0112] 其中，表示第个测试照片块 的稀疏表示，每个稀疏表示 都包含稀疏系数值 以及稀疏系数值在稀疏编码中被求解的顺序 。

[0113] 所述根据相似性尺度距离所述第二投影系数集合在第一投影系数集合中查找到距离最近的前K个系数，根据获得的前K个系数在第二图像块子集查找到匹配的图像块，具体为：对合成训练画像块集合 通过贪婪搜索获取每个测试照片块 对应的待选择画像块集合 ：

[0114] 7a) 设 ， ；

[0115] 7b) 从顺序集合 中选出与求解顺序 中第个值相同的所有稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断 中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7c)；

[0116] 7c) 从稀疏系数值集合 中选出与稀疏系数值 中第个值欧式距离最相近的 稀疏表示 ，并将这些选出的稀疏表示记为集合 ，并用 值替代稀疏表示集合 ，再判断中稀疏表示的个数 ：如果 ，则转到步骤7e)，如果 ，则保留替代前的稀疏表示集合 ，并转到步骤7e)，如果 ，则转到步骤7d)；

[0117] 7d) 令 ，返回步骤7b)；

[0118] 7e) 将稀疏表示集合 中所有 个稀疏表示对应的合成训练画像块集合 中的个画像块作为稀疏表示 对应的测试照片块 的待选择画像块集合 ：如果，则利用复制操作将待选择画像块从 个增加到 个，如果 ，则利用测试照片块的高频信息与待选择画像块的高频信息的欧式距离将待选择画像块从 个减少为 个，使待选择画像块集合从 变成 。

[0119] 所述将查找到的图像块合成得到初始合成图像，具体为：根据求解马尔可夫随机场模型从待选择画像块集合 中得到每个测试照片块 最终对应的训练画像块，并将这些训练画像块进行融合得到与测试照片 对应的初始合成画像 ， ；

[0120] 所述将对所述目标画像，测试照片以及初始合成图像进行分块得到第三图像块子集，将所述测试照片分块后得到的图像块在第三图像块子集内根据图像块的多尺度特征查找匹配的图像块，将所述查找到的图像块构成新合成图像，具体为：提取初始合成画像 、测试照片 以及目标画像的像素值特征以及SURF特征，提取初始合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，获取测试照片块 对应的待选择画像块集合 ，根据所述待选择画像块集合 获取新合成画像。

[0121] 所述提取合成画像 、测试照片 以及目标画像的像素值特征以及SURF特征，具体为

[0122] 9a)  将合成画像 划分为相同大小及相同重叠程度的画像块集合：，表示合成画像 中的第个画像块；

[0123] 9b)  将测试照片 划分为相同大小及相同重叠程度的测试照片块集合：，表示测试照片 中的第个照片块；

[0124] 9c)  将目标画像划分为相同大小及相同重叠程度的训练画像块集合：，表示目标画像中的第个画像块；

[0125] 9d) 对画像块集合 中的每个块提取像素值特征和SURF特征；

[0126] 9e) 对测试照片块集合 中的每个块提取像素值特征和SURF特征；

[0127] 9f) 对训练画像块集合 中的每个块提取像素值特征和SURF特征；

[0128] 所述提取合成画像 、测试照片 以及目标画像的多尺度特征的投影系数，具体为：

[0129] 10a) 对合成画像 建立高斯金字塔 和拉普拉斯金字塔 ，再对高斯金字塔 的每层画像求解水平一阶导数 和二阶导数 以及竖直一阶导数 和二阶导数 。因此，合成画像 拥有6个多尺度特征：合成画像 的高斯金字塔 ，合成画像 的拉普拉斯金字塔 ，合成画像 的水平一阶导数金字塔 ，合成画像 的水平二阶导数金字塔 ，合成画像 的竖直一阶导数金字塔 ，合成画像 的竖直二阶导数金字塔 ；

[0130] 10b) 对目标画像建立高斯金字塔 和拉普拉斯金字塔 ，再对高斯金字塔 的每层画像求解水平一阶导数 和二阶导数 以及竖直一阶导数 和二阶导数 。因此，目标画像拥有6个多尺度特征：目标画像的高斯金字塔 ，目标画像的拉普拉斯金字塔 ，目标画像的水平一阶导数金字塔 ，目标画像的水平二阶导数金字塔 ，目标画像的竖直一阶导数金字塔 ，目标画像的竖直二阶导数金字塔 ；

[0131] 10c) 对测试照片 建立高斯金字塔 和拉普拉斯金字塔 ，再对高斯金字塔 的每层照片求解水平一阶导数 和二阶导数 以及竖直一阶导数 和二阶导数 。因此，测试照片 拥有6个多尺度特征：测试照片 的高斯金字塔 ，测试照片 的拉普拉斯金字塔 ，测试照片 的水平一阶导数金字塔 ，测试照片 的水平二阶导数金字塔 ，测试照片 的竖直一阶导数金字塔 ，测试照片 的竖直二阶导数金字塔 ；

[0132] 10d) 将步骤10a)-10c)中的各多尺度特征投影到步骤2中的主成分以获取各多尺度特征的投影系数；

[0133] 所述获取测试照片块 对应的待选择画像块集合 ，具体为：

[0134] 11a) 将每个测试照片块 从测试照片块集合 中利用像素值特征搜索求得 个最近邻；则画像块集合 中 个对应的画像块 即为搜索
得到的待选择画像块；

[0135] 11b) 将每个测试照片块 从训练画像块集合 中利用多尺度特征的投影系数搜索求得 个最近邻；

[0136] 11c) 将每个测试照片块 对应的画像块 从训练画像块集合中利用多尺度特征的投影系数搜索求得 个最近邻；

[0137] 11d) 对于每个测试照片块 ，当测试照片块 的高频分量大于某个指定值时，将测试照片块 在步骤11a)-11c)中获得的 个最近邻利用测试照片块 的各特征和 个最近邻的各特征的欧式距离进行最优化重新排序：

[0138] 11d1) 利用提取的图像块多尺度特征的投影系数以及图像块像素值特征和SURF特征，先求解每个测试照片块 与 个最近邻对应特征的欧式距离：

[0139] ， ；

[0140] 其中， 表示测试照片块 与第 个近邻画像块在第个特征上的欧式距离，表示图像块所提取的特征个数；

[0141] 11d2) 对步骤11d1)得到的 ， ，通过求解以下最优化问题得到测试照片块 与第 个最近邻所有对应特征的欧式距离 的权重 。

[0142]

[0143] 其中，是对角矩阵，对角元素为 ；

[0144] 11d3) 利用步骤11d1)得到的欧式距离 与步骤11d2)得到的权重 ，按照得到每个测试照片块 与 个最近邻新的距离从而进行重新排序。

[0145] 11e) 对于每个测试照片块 ，当测试照片块 的高频分量小于等于某个指定值时，直接将测试照片块 在测试照片块集合 中利用像素值特征搜索求得个最近邻；

[0146] 所述获取新合成画像，具体为：

[0147] 从步骤11a)-11e) 得到的测试照片块 对应的 个最近邻中选择前 个最近邻，通过求解马尔可夫随机场模型得到每个测试照片块 最终对应的训练画像块，然后经过融合得到测试照片 对应的新合成画像 ， ；

[0148] 所述将新合成图像替代初始合成图像，重复迭代上述步骤合成最终合成图像，具体为：新合成画像代替初始合成画像，迭代执行步骤9a）11e），直到迭代次数达到指定值~时停止，则新合成画像 ， ，为最终合成画像。

[0149] 实验数据：

[0150] 本发明的效果可以通过以下仿真实验做进一步的说明。

[0151] 1.仿真条件

[0152] 本发明是在中央处理器为Intel(R) Core(TM) i3-2100 3.10GHz、内存4G、WINDOWS 7操作系统上，运用美国Mathworks公司开发的MATLAB R2013b进行仿真。

[0153] 数据库采用CUFS数据库，该数据库包含CUHK student数据库、AR数据库和XM2VTS数据库。

[0154] AR数据库见文献“Martinez A M. The AR face database[J]. CVC Technical Report, 1998, 24.”。

[0155] XM2VTS数据库见文献“Messer K, Matas J, Kittler J, et al. XM2VTSDB: The extended M2VTS database[C]//Second international conference on audio and video-based biometric person authentication. 1999, 964: 965-966.”。

[0156] 2.仿真内容

[0157] 取CUHK student数据库中1个男性画像作为目标画像。测试照片样本集 则为CUFS数据库当中CUHK student数据库剩下的187张照片、AR数据库中123张照片和XM2VTS数据库中295张照片。

[0158] 为了证明本发明能合成多种风格的画像，目标画像除了取CUFS数据库中的画像，还从网上下载了两张不同风格的画像作为目标画像，如图2所示。其中，图2(1)和图2(2)分别为网上下载的两张不同风格的画像，图2(3)为CUHK student数据库中1个男性画像。

[0159] 利用本发明和现有基于局部线性的方法LLE、基于马尔可夫随机场模型的方法MRF以及基于马尔可夫权重场模型的方法MWF，在CUFS数据库上进行照片到画像的合成仿真实验，结果如图3所示。其中：

[0160] 图3(a)为原始照片，

[0161] 图3(b)为现有基于局部线性的方法LLE合成的画像，

[0162] 图3(c)为现有基于马尔可夫随机场模型的方法MRF合成的画像，

[0163] 图3(d)为现有基于马尔可夫权重场模型的方法MWF合成的画像，

[0164] 图3(e)为本发明方法用图2(3)作为目标画像时合成的画像。

[0165] 由图3实验结果可见，图3(b) 图3(d)中的第一张画像都没有发卡，第二张画像都~没有刘海，其余的四张画像都没有眼镜，而且图3(b) 图3(d)中的第三张画像到第六张画像~
的背景都存在噪声。而图3(e)所示的本发明的第一张画像有发卡，第二张画像有刘海，其余的四张画像有眼镜，而且所有合成画像的背景都很干净。

[0166] 在用图2中的三种不同风格的画像依次作为目标画像的情况下，本发明方法在CUFS数据库上进行照片到画像的合成仿真实验，结果如图4所示。其中：

[0167] 图4(a)为原始照片，

[0168] 图4(b)为本发明方法用图2(1)作为目标画像时合成的画像，

[0169] 图4(c)为本发明方法用图2(2)作为目标画像时合成的画像，

[0170] 图4(d)为本发明方法用图2(3)作为目标画像时合成的画像。

[0171] 由图4实验结果可见，图4(b)的风格与图2(1)相似，图4(c)的风格与图2(2)相似，图4(d)的风格与图2(3)相似，此外图4(b) 图4(d)所示的本发明的第一张画像有发卡，第二~张画像有刘海，其余的四张画像有眼镜，而且所有合成画像的背景都很干净。

[0172] 在用图2中的三种不同风格的画像依次作为目标画像的情况下，本发明方法在网上下载的合照上进行照片到画像的合成仿真实验，结果如图5所示。其中：

[0173] 图5(a)为原始照片，

[0174] 图5(b)为本发明方法用图2(1)作为目标画像时合成的画像，

[0175] 图5(c)为本发明方法用图2(2)作为目标画像时合成的画像，

[0176] 图5(d)为本发明方法用图2(3)作为目标画像时合成的画像。

[0177] 由图5实验结果可见，图5(b)的风格与图2(1)相似，图5(c)的风格与图2(2)相似，图5(d)的风格与图2(3)相似，此外图5(b) 图5(d)所示的本发明的合成画像均不受测试图~像的大小以及人脸数量的影响。

[0178] 综上所述：本发明由于在用多尺度金字塔作为特征时借助了迭代的思想，并在以稀疏表示中稀疏系数值和稀疏系数值在稀疏编码中被求解的顺序作为特征的情况下，利用贪婪搜索的方法，使得在以单张目标画像作为训练集的情况下，在测试照片与训练集照片背景不同时，能忽略背景的影响而合成测试照片对应的画像；在测试照片中存在训练集照片中都不存在的非人脸部件时，能对测试照片合成包括非人脸部件在内的对应画像；此外，还能忽略目标画像的风格和测试照片的大小以及测试照片中所包容的人脸数目等因素进行画像合成，验证了本发明的实用性。

[0179] 以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。