CN112819922B - 一种基于连续线条的人物肖像画生成方法 - Google Patents

一种基于连续线条的人物肖像画生成方法 Download PDF

Info

Publication number
CN112819922B
CN112819922B CN202110143469.0A CN202110143469A CN112819922B CN 112819922 B CN112819922 B CN 112819922B CN 202110143469 A CN202110143469 A CN 202110143469A CN 112819922 B CN112819922 B CN 112819922B
Authority
CN
China
Prior art keywords
value
sum1
pixel
portrait
point
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
Application number
CN202110143469.0A
Other languages
English (en)
Other versions
CN112819922A (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.)
Zhengzhou University of Light Industry
Original Assignee
Zhengzhou University of Light Industry
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 Zhengzhou University of Light Industry filed Critical Zhengzhou University of Light Industry
Priority to CN202110143469.0A priority Critical patent/CN112819922B/zh
Publication of CN112819922A publication Critical patent/CN112819922A/zh
Application granted granted Critical
Publication of CN112819922B publication Critical patent/CN112819922B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/40—Filling a planar surface by adding surface attributes, e.g. colour or texture
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Image Processing (AREA)

Abstract

本发明提出了一种基于连续线条的人物肖像画生成方法，步骤为：采集彩色的人物肖像图像；采用平均值法得到灰度图像；生成包围肖像区域的圆形蒙版；在圆形蒙版的圆形轮廓上均匀设定多个锚点并计算坐标；得到肖像图像，在肖像图像上选取两个锚点并连线，计算连线经过的像素点的像素值的均值，将最小值对应的锚点作为起始点；将起始点存入数组中，依次计算并顺序存入数组中，以黑色的连线顺序连接数组中的锚点得到初始线肖像画；去除冗余保留连线交点得到黑色线肖像画；检测人物肖像图像的皮肤区域，顺序连接皮肤区域内的锚点并叠加得到彩色线肖像画。本发明自动化程度高，计算速度快，人物肖像还原精度高，为人物肖像生成提供了技术支撑。

Description

一种基于连续线条的人物肖像画生成方法
技术领域
本发明涉及计算机图形学、图像处理的技术领域，尤其涉及一种基于连续线条的人物肖像画生成方法。
背景技术
肖像从不同角度上有不同的理解，美术意义(或摄影)上的肖像，是指通过绘画、摄影、文字等艺术手段，使肖像权人的人物形象在物质载体下再现的一种观赏造型作品。近年来，数字媒体技术正逐渐融入生活，极大地影响着人们的生活和工作方式，利用计算机生成人脸及肖像的图形在很多领域具有广泛的应用，如视频会议、网络游戏、移动数字娱乐中，使用特定算法生成的人脸、肖像等图形代替真实人脸，在数据处理、数据传送和营造气氛等方面具有重要的实用价值和意义，不仅可以提高数据传送的速度，而且能够很好地营造轻松、和谐的气氛。
在艺术设计等领域，采用特定的创新方法、肖像生成算法以及简易装置生成具有特定风格的肖像，并以实物的方式进行再现，通常能够以较高的价格进行销售并获利，具有较强的实际需求。
人脸及肖像生成能简单而形象地对人物进行描绘，涵盖了多个学科和领域，包括计算机视觉、计算机图形学和心理学等，一直以来都是图像处理领域一个热门的研究方向，人脸及肖像等的生成及风格化自诞生以来，就与艺术创作紧密的联系。随着计算机及图像处理技术的发展，人们在风格化研究上做了很多尝试，涌现出大量算法及应用。
申请号为201710868480.7的发明专利提出一种基于特征提取的人物肖像漫画风格化生成方法，提取人脸特征点后将人脸分解为五官、脸型、发型、服饰等组件，根据这些组件的特征与组件素材库中已有的组件素材进行匹配，最后将得到的各匹配组件合成为人物肖像，并对其进行夸张渲染；采用的肖像风格化算法将肖像的结构特征与风格化信息分开处理。申请号为201610895094.2的发明专利公开了一种自动描绘人脸线条画的方法，可以利用机器人自动绘制一幅边缘连续的人脸肖像画，并能够在肖像画上进行个性化签名，是一种自动化的、形式新颖的绘制方法。申请号为201010163881.0的发明专利提供了一种基于人脸图像的卡通人脸图形生成方法和装置，其中的方法具体包括：识别所述人脸图像，得到人脸基准线；基于所述人脸基准线以及人脸图形构件库，匹配组合得到卡通人脸图形，用以使得生成的卡通人脸图形既能满足实时性需求，又能够贴合真实的人脸图像。已有的人脸肖像生成方法，大多基于模板匹配、特征提取和边缘检测的方法，数据库中的模板数量有限，不能很好的表征现实世界中众多的人脸肖像，特征提取的参数有限，所能生成并表征的人物肖像数量有限。
发明内容
针对现有的人脸肖像生成方法所能生成并表征的人物肖像数量有限，不能很好的表征现实世界中众多的人脸肖像的技术问题，本发明提出一种基于连续线条的人物肖像画生成方法，思路简单，操作方便，在一定程度上能够贴合真实的人脸图像，并能输出所有的锚点编号、锚点坐标及锚点的连接顺序，便于后续的机械臂及人工连线；且具有自动化程度高、计算速度快和人物肖像还原精度高等优点。
为了达到上述目的，本发明的技术方案是这样实现的：一种基于连续线条的人物肖像画生成方法，其步骤如下：
步骤一：通过相机采集或下载待处理的彩色的人物肖像图像；
步骤二：采用平均值法对步骤一中的人物肖像图像进行图像灰度化，得到灰度图像；
步骤三：对灰度图像进行调整，并生成与调整后灰度图像大小相同的包围肖像区域的圆形蒙版；
步骤四：在圆形蒙版的圆形轮廓上均匀设定多个锚点，并计算各个锚点的坐标；
步骤五：根据步骤三得到的圆形蒙版和灰度图像得到肖像图像，在肖像图像上的圆形轮廓上选取任意两个锚点并连线，计算连线经过的像素点的像素值的均值，将均值中的最小值对应的锚点作为起始点；
步骤六：将步骤五得到的起始点存入数组中，依次计算后续连接点并顺序存入数组中，在与肖像图像大小相同的圆形空白画布上，以黑色的连线将数组中的起始点和后续连接点顺序连接，得到初始线肖像画；
步骤七：去除步骤六中初始线肖像画上的冗余，保留连线的交点，得到黑色线肖像画；
步骤八：检测步骤一中人物肖像图像的皮肤区域并作为感兴趣区域，循环步骤二-步骤七以肤色的连线顺序连接表征皮肤区域的起始点和后续连接点，得到肤色线肖像画，将肤色线肖像画与步骤七得到的黑色线肖像画叠加得到彩色线肖像画。
所述步骤二中平均值法的实现方法为：提取人物肖像图像的红绿蓝三个通道的平均值作为灰度值：
Gray＝(Rc+Gc+Bc)/3；
其中，Gray表示灰度图像的像素值，Rc、Gc、Bc分别表示人物肖像图像的红、绿、蓝三个通道的像素值。
所述步骤三中对灰度图像进行调整的方法为：
以人物肖像区域为中心，以双线性插值的方式将灰度图像的大小调整为700×700像素；
所述圆形蒙版生成的方法为：
生成一个大小为700×700像素的画布图像，以画布图像的中心为原点，350个像素的宽度为半径在画布图像上画圆，并将圆形区域内的像素值设置为0，圆形区域外的像素值设置为1，得到包围肖像区域的圆形蒙版。
所述步骤四的实现方法为：在圆形蒙版的圆形轮廓的圆周上每隔1.8°均匀的设定200个锚点，以圆形轮廓的顶部最高位置处的点为初始锚点，并按顺时针方向依次计算其他锚点的坐标，依次存储到集合pinsXY(i1,j1)中，其中，i1表示锚点的编号，取值范围为1～200；j1表示锚点的x坐标或者y坐标，取值范围为1～2，当j1＝1时表示x坐标，当j1＝2时表示y坐标。
所述肖像图像是将圆形蒙版中的值取非，然后将圆形蒙版与灰度图像进行点乘得到的；
所述计算连线经过的像素点的像素值的均值的方法为：
从半径为350的圆形轮廓上的200个锚点中任意挑选两个锚点，将这两个锚点连线，统计连线所经过的肖像图像的像素点的数量Numij及像素点的像素值之和Sumij，其中，i的取值为1～199，j的取值为(i+1)～200；如果数量Numij>30，则计算该连线所经过的像素点的像素值的均值Aveij＝Sumij/Numij；否者j＝j+1；统计连线所经过像素点的像素值的均值的集合Ave＝{Aveij}中的最小值Avemn，其中，m∈i，n∈j，且m<n，则选取锚点m作为线肖像画的起始点。
所述统计连线所经过的肖像图像的像素点的数量Numij及像素点的像素值之和Sumij时，计算任意两个锚点A和锚点B的连线所经过肖像图像的像素点的坐标的方法为：锚点A的坐标为(pinsXY(i,1)，pinsXY(i,2))，锚点B的坐标为(pinsXY(j,1)，pinsXY(j,2))，锚点A和锚点B之间连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)；则：
(1)当锚点A的y坐标yA＝pinsXY(i,2)与锚点B的y坐标yB＝pinsXY(j,2)相等，则锚点A和锚点B处于同一行，判断锚点A和锚点B的x坐标值：
若xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA-xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA-xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1-1、y坐标的值为AB_xy(mdsum1,2)＝yA；
若xB>xA，则锚点A和锚点B之间连线经过的像素点的数量为(xB-xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB-xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1-1、y坐标的值为AB_xy(mdsum1,2)＝yA；
(2)锚点A的x坐标xA＝pinsXY(i,1)与锚点B的x坐标xB＝pinsXY(j,1)相等，则锚点A和锚点B处于同一列，判断锚点A和锚点B的y坐标值：
若yA>yB，则锚点A和锚点B之间连线经过的像素点的数量为(yA-yB+1)个，连线经过的像素点的坐标的集合AB_xy中，索引mdsum1的取值为1～(yA-yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1-1、x坐标的值为AB_xy(mdsum1,1)＝xA；
若yB>yA，则锚点A和锚点B之间连线经过的像素点的数量为(yB-yA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yB-yA+1)，索引mdsum2的取值为1～2；连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1-1、x坐标值为AB_xy(mdsum1,1)＝xA；
(3)锚点A和锚点B在x方向的偏移xOffset＝|xA–xB|，锚点A和锚点B在y方向的偏移yOffset＝|yA–yB|，若xOffset≥yOffset，锚点A和锚点B的连线的斜率的绝对值小于等于1，则：
若yA>yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA–xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA–yB)/(xA–xB)|)，其中，round()为四舍五入取整函数；
若yA>yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB–xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA yB)/(xA–xB)|)；
若yA<yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA–xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)；
若yA<yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB–xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)；
(4)锚点A和锚点B在x方向的偏移为xOffset＝|xA–xB|，在y方向的偏移为yOffset＝|yA–yB|，若xOffset<yOffset，锚点A和锚点B的连线的斜率的绝对值大于1，则：
若yA>yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yA–yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1–1；
若yA>yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yA–yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标值为AB_xy(mdsum1,2)＝yB+mdsum1–1；
若yA<yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yB–yA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标值为AB_xy(mdsum1,2)＝yA+mdsum1–1；
若yA<yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yB–yA+1)，mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1–1；
根据连线上各个像素点的x坐标和y坐标的值，通过Pr(x,y)提取像素点的像素值即灰度值，将锚点A和锚点B之间连线经过的像素点的灰度值相加得到像素值之和Sumij。
所述步骤六的实现方法为：
Step1：通过步骤五的处理，得到锚点m为线画肖像的起始点，锚点n是后续连接点，且锚点m和锚点n的连线经过的所有像素点的灰度值的均值最小；将数值m存入数组sequence，将肖像图像的锚点m和锚点n的连线经过的所有像素点的灰度值设定为255；
Step2：以锚点n为出发点，分别计算锚点n到锚点1～锚点(m–1)、锚点(m+1)～锚点(n–1)号及锚点(n+1)～锚点200的连线经过的像素的个数，并排除其中连线所经过像素数少于30的锚点；然后采用步骤五的方法计算得到锚点pi是后续连接点，使得锚点n和锚点pi的连线经过的所有像素的像素值的均值最小；将数值n存入数组sequence中，将肖像图像的锚点n和锚点pi的连线经过的所有像素点的灰度值设定为255；
Step3：同理，以得到的后续连接点pi作为出发点，循环Step2，得到该出发点的后续连接点pi+1，数值pi+1存入数组sequence；将肖像图像的pi和pi+1号锚点的连线所经过的所有像素点的灰度值设定为255；
Step4：以黑色的连线依次在半径为350的圆形空白画布上顺序连接数组sequence中的锚点，得到初始线肖像画。
所述步骤七中黑色线肖像画的获取方法为：提取初始线肖像画的第2～699行、第2～699列构成感兴趣区域，依次对感兴趣区域中的每一个像素点Pixel(m1,n1)的八邻域进行分析，将像素点Pixel(m1,n1)的八个邻域的像素值求和得到邻域值Sneighbourhood，其中，m1的取值范围为2～699，n1的取值范围为2～699；如果像素点Pixel(m1,n1)的像素值为0，并且邻域值Sneighbourhood>4，则将像素点Pixel(m1,n1)的像素值设定为1，得到黑色线肖像画。
所述步骤八中检测人物肖像图像的皮肤区域的方法为：
将人物肖像图像从RGB颜色空间转换为YCbCr颜色空间，其中，Y是指亮度分量，Cb指蓝色色度分量，而Cr指红色色度分量，且归一化的转换矩阵HI为：
则亮度分量Y＝HI(1,1)*Rc+HI(1,2)*Gc+HI(1,3)*Bc+16，蓝色色度分量Cb＝HI(2,1)*Rc+HI(2,2)*Gc+HI(2,3)*Bc+128，红色色度分量Cr＝HI(3,1)*Rc+HI(3,2)*Gc+HI(3,3)*Bc+128，Rc、Gc、Bc分别表示原始的彩色的人物肖像图像的红、绿、蓝三个通道的像素值；
若亮度分量Y、蓝色色度分量Cb或红色色度分量Cr的像素的值大于236，则将像素值修改为236；若亮度分量Y、蓝色色度分量Cb或红色色度分量Cr的像素的值小于15，则将像素值修改为15，得到YCbCr图像；
遍历YCbCr图像中的每个像素点，计算各个像素点的肤色值skin为：
skin(ii,jj)＝(cos(2.54)*(Cb(ii,jj)–109.39)+sin(2.54)*(Cr(ij,jj)-152.03)–1.61)^2/(a^2)+
(-sin(2.54)*(Cb(ij,jj)-109.39)+cos(2.54)*(Cr(ij,jj)-152.03)–2.42)^2/(b^2)；
其中，Cb(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的蓝色色度分量Cb的值；Cr(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的红色色度分量Cr的值；a、b为自定义的椭圆肤色模型的长轴、短轴，且a＝25.39；b＝14.03；
如果skin<＝1，则将第ii行、第jj列的像素点标记为皮肤点，所有皮肤点的集合为皮肤区域。
所述肤色线肖像画的实现方法为：采用步骤二-步骤七中生成黑色线肖像画的方法，再次生成表示皮肤区域的线画肖像的起点及后续连接点，并以人物肖像图像中的皮肤区域的像素点的R通道的均值作为连线的R值，以人物肖像图像中的皮肤区域的像素点的G通道的均值作为连线的G值，以人物肖像图像中的皮肤区域的像素点的B通道的均值作为连线的B值，以三个通道的R值、G值、B值为颜色的连线，依次连接表示皮肤区域的线画肖像的起点及后续连接点，构成肤色线肖像画。
与现有技术相比，本发明的有益效果：对所获取的人物彩色肖像图像进行处理，对灰度图像计算得到起点和后续连接点，以黑色线依次连接起点和后续连接点得到线画肖像；提出了皮肤检测方法，将肖像中的皮肤区域作为感兴趣区域，计算得到皮肤区域的起点和后续连接点，以肤色作为连线的颜色，依次连接皮肤区域的起点和后续连接点，并叠加到黑色线的线画肖像上，最后实现彩色线画肖像。本发明具有下述有益效果：(1)自动化程度高，对人物彩色肖像图像自动生成彩色线画肖像；(2)计算速度快，在安装window 10企业版的ThinkPad T460P笔记本电脑上(Intel(R)Core(TM)i7-6700HQ CPU@2.6GHz,内存为16GB)，运行时间小于5分钟；(3)人物肖像还原精度高，经过多位观察员的人工判读，能够有效的配准原始的人物彩色肖像图像和所生成的线肖像画。本发明为人物肖像生成提供了技术支撑，同时研究思路与方法也可以推广应用到其它字符、标志等图案的生成，具有很高的实用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的流程图。
图2为原始的人物彩色肖像图像的灰度图。
图3为本发明图2的灰度图像。
图4为本发明包围肖像的圆形蒙版。
图5为本发明在图4的轮廓上设置锚点的示意图。
图6为本发明的感兴趣区域的肖像图像。
图7为本发明查找线肖像画起点的流程图。
图8为本发明两个锚点处于同一行的示意图。
图9为本发明两个锚点处于同一列的示意图。
图10为本发明两个锚点的连线的斜率的绝对值小于等于1的示意图。
图11为本发明两个锚点的连线的斜率的绝对值大于1的示意图。
图12为本发明锚点之间的1000次顺序连线构成的线肖像画的示意图。
图13为本发明的不同次数的顺序连线所构成的线肖像画的示意图，其中，(a)为100次，(b)为300次，(c)为500次，(d)为1500次。
图14为本发明对图12去除冗余的肖像画。
图15为本发明中图3中皮肤区域的示意图。
图16为本发明中皮肤区域的线肖像画。
图17为本发明中彩色的线肖像画。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，一种基于连续线条的人物肖像画生成方法，其步骤如下：
步骤一：通过相机采集或下载待处理的彩色的人物肖像图像。
人物肖像图像的来源及格式要求：可应用于本发明的肖像的来源非常广泛，可以通过相机采集目标人员的图像，从而获得人物肖像图像；也可以从网络上下载相关人物肖像的图像，也可以是来自于人物雕塑、油画等形式的图像。常见的jpg、png、bmp、tif等格式的肖像图像均可作为本发明的数据源。肖像在图像中所占的幅面不限，但人物肖像应尽可能清晰的显示在图像中。图像的宽度和高度尽量相同，确保肖像的宽高比在合理的范围内。背景中的干扰物应尽量少，以保证生成线画肖像的效果和质量。待处理的彩色的人物肖像图像的原始图像Ori_img如图2所示，原始图像Ori_img的大小800×800。
步骤二：采用平均值法对步骤一中的人物肖像图像进行图像灰度化，得到灰度图像。
将待处理的人物肖像图像即原始图像Ori_img进行图像灰度化，常用的图像灰度化方法是平均值法，即提取红绿蓝三个通道的平均值作为灰度值：
Gray＝(Rc+Gc+Bc)/3 (1)
其中，Gray表示灰度图像的像素值，Rc、Gc、Bc分别表示原始的彩色的人物肖像图像的红、绿、蓝三个通道的像素值，从而得到灰度图像Gray_img，如图3所示。图像灰度化从彩色图像的红、绿、蓝三个通道转变为灰度这一个通道，减少了计算量。
步骤三：对灰度图像进行调整，并生成与调整后灰度图像大小相同的包围肖像区域的圆形蒙版。
以人物肖像区域为中心进行调整，采用matlab中imresize()函数实现。以双线性插值的方式将灰度图像Gray_img的大小调整为700×700像素。这样全部流程的整体的计算量适中，计算时间较短，且最终在画布上生成的线画肖像的图像质量也能满足人工评判的要求。
设定一个大小为700×700像素的画布图像，以该画布图像的中心为原点，350个像素的宽度为半径画圆，并将圆形区域内的像素值设置为0，即用黑色表示；圆形区域外的像素值设置为1，即用白色表示，如图4所示，该图像即为灰度图像Gray_img的蒙版，命名为圆形蒙版Mask_img。方便后续在圆形画布上进行处理和分析，最终的结果也可以用于指导圆形画框上的实物连接，形成可供销售的线画肖像商品。
步骤四：在圆形蒙版的圆形轮廓上均匀设定多个锚点，并计算各个锚点的坐标。
在圆形蒙版Mask_img的圆形轮廓(圆周)上每隔1.8°均匀的设定200个锚点，并计算各个锚点的坐标。以圆形轮廓的顶部最高位置处，即圆形轮廓的零点方向的点为初始锚点，即1号锚点(Anchor_point1)，并按顺时针方向依次计算后序锚点的坐标，依次存储到集合pinsXY(i1,j1)中，其中，i1表示锚点的编号，取值范围为1～200；j1表示锚点的x坐标或者y坐标，取值范围为1～2，当j1＝1时表示x坐标，当j1＝2时表示y坐标，如图5所示。
步骤五：根据步骤三得到的圆形蒙版和灰度图像得到肖像图像，在肖像图像上的圆形轮廓上选取任意两个锚点并连线，计算连线经过的像素点的像素值的均值，将均值中的最小值对应的锚点作为线肖像画的起始点。
将圆形蒙版Mask_img中的值取非，即将圆形蒙版Mask_img中值为1的改为0，将值为0的改为1。然后将圆形蒙版Mask_img与灰度图像Gray_img进行点乘，即可得到感兴趣区域的肖像图像Portrait_area，如图6所示。
从半径为350的圆形轮廓上的200个锚点中任意挑选两个锚点，将这两个锚点连线，统计该连线所经过的肖像图像Portrait_area的像素点的数量Numij及像素点的像素值之和Sumij，其中，i的取值为1～199，j的取值为(i+1)～200。如果数量Numij>30，则计算该连线所经过的像素点的像素值的均值，即Aveij＝Sumij/Numij。阈值设定为30，防止出现相邻或相近的锚点间的连线的长度太短，从而影响生成先肖像画的效果。统计连线所经过像素点的像素值的均值的集合Ave中的最小值Avemn，其中，m∈i，n∈j，且m<n，则选取m号锚点Anchor_point m作为线肖像画的起始点，线肖像画的起始点的计算方法如图7所示。最小值Avemn对应的是两个锚点连线经过的像素点的灰度值的均值最小，该连线可最优先的体现头发、眉毛、眼睛、鼻子和嘴等灰度值较低的区域。
其中，任意两个锚点A、B的连线所经过肖像图像Portrait_area的像素点的计算方法如下：锚点A的坐标为(pinsXY(i,1)，pinsXY(i,2))，锚点B的坐标为(pinsXY(j,1)，pinsXY(j,2))，即锚点A的x坐标为xA＝pinsXY(i,1)，锚点A的y坐标为yA＝pinsXY(i,2)，锚点B的x坐标为xB＝pinsXY(j,1)，锚点B的y坐标为yB＝pinsXY(j,2)。两个锚点A、B之间连线经过的像素点的坐标的集合为AB_xy。则两个锚点A、B的位置关系包含以下四种情况：处于同一行、处于同一列、两点的连线的斜率的绝对值小于等于1、两点的连线的斜率的绝对值大于1，从而对这两个锚点A、B的连线经过的像素点进行统计和分析。
(1)两个锚点A、B处于同一行
如果锚点A的y坐标yA＝pinsXY(i,2)与锚点B的y坐标yB＝pinsXY(j,2)相等，则两个锚点A、B处于同一行。然后判断锚点A和锚点B的x坐标值：
如果xA>xB，示意图如图8所示，则两个锚点之间连线经过的像素点的数量为(xA-xB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xA-xB+1)，mdsum2的取值为1～2；对应的连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1-1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA。
如果xB>xA，则两个锚点之间连线经过的像素点的数量为(xB-xA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xB-xA+1)，mdsum2的取值为1～2；对应的连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1-1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA。
(2)两个锚点A、B处于同一列
如果锚点A的x坐标xA＝pinsXY(i,1)与锚点B的x坐标xB＝pinsXY(j,1)相等，则两个锚点A、B处于同一列，然后判断锚点A和锚点B的y坐标值：
如果yA>yB，示意图如图9所示，则两个锚点之间连线经过的像素点的数量为(yA-yB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yA-yB+1)，mdsum2的取值为1～2；对应的连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1-1，各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA。
如果yB>yA，则两个锚点之间连线经过的像素点的数量为(yB-yA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yB-yA+1)，mdsum2的取值为1～2；对应的各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1-1，各个像素点的x坐标值为AB_xy(mdsum1,1)＝xA。
(3)两个锚点A、B的连线的斜率的绝对值小于等于1
设定两个锚点A、B在x方向的偏移为xOffset＝|xA–xB|，在y方向的偏移为yOffset＝|yA–yB|。如果xOffset≥yOffset，则表示两个锚点A、B的连线的斜率的绝对值小于等于1，斜率的坡度较为舒缓。
如果yA>yB且xA>xB，示意图如图10所示，则两个锚点A、B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xA–xB+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA–yB)/(xA–xB)|)，其中，round()为四舍五入取整函数。
如果yA>yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xB–xA+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA yB)/(xA–xB)|)，其中round()函数的功能为通过四舍五入取整。
如果yA<yB且xA>xB，则两个锚点A、B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xA–xB+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)，其中round()函数的功能为通过四舍五入取整。
如果yA<yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(xB–xA+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)，其中，round()函数的功能为通过四舍五入取整。
(4)两个锚点A、B的连线的斜率的绝对值大于1
设定两个锚点A、B在x方向的偏移为xOffset＝|xA–xB|，在y方向的偏移为yOffset＝|yA–yB|。如果xOffset<yOffset，则表示两个锚点A、B的连线的斜率的绝对值大于1，斜率的坡度较为陡峭。
如果yA>yB且xA>xB，示意图如图11所示，则两个锚点A、B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yA–yB+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1–1，其中，round()函数的功能为通过四舍五入取整。
如果yA>yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yA–yB+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，各个像素点的y坐标值为AB_xy(mdsum1,2)＝yB+mdsum1–1，其中，round()函数的功能为通过四舍五入取整。
如果yA<yB且xA>xB，则两个锚点A、B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yB–yA+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，各个像素点的y坐标值为AB_xy(mdsum1,2)＝yA+mdsum1–1，其中，round()函数的功能为通过四舍五入取整。
如果yA<yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yB–yA+1)，mdsum2的取值为1～2；对应的各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1–1，其中，round()函数的功能为通过四舍五入取整。
根据连线上各个像素点的x坐标和y坐标的值，通过Pr(x,y)提取像素点的像素值(即灰度值)；将锚点A和锚点B之间连线经过的像素点的灰度值相加即得到像素值之和Sumij。
步骤六：将步骤五得到的起始点存入数组中，依次计算线肖像画的后续连接点并顺序存入数组中，在与肖像图像大小相同的圆形空白画布上，以黑色的连线将数组中的起始点和后续连接点顺序连接，得到初始线肖像画。
Step1：通过上面的处理过程，得到了第m号锚点Anchor_point m为线画肖像的起点Start_point，并得到n号锚点是后续连接点，并且m和n号锚点的连线所经过的所有像素点的灰度值的均值最小。将数值m存入数组sequence。将肖像图像Portrait_area的m和n号锚点的连线所经过的所有像素点的灰度值设定为255，以减小该条连线对后续所要产生的连线的影响。
Step2：以n号锚点为出发点，分别计算该锚点到1号～(m–1)号、(m+1)号～(n–1)号及(n+1)号～200号共198个锚点(排除了m号锚点，因为m号锚点刚和n号锚点连线)的连线所经过的像素的个数，并排除其中连线所经过像素数少于30的锚点。然后采用步骤五的方法，并计算得到pi号锚点是后续连接点，使得n号和pi号锚点的连线所经过的所有像素的像素值的均值最小。将数值n存入数组sequence。将肖像图像Portrait_area的n和pi号锚点的连线所经过的所有像素点的灰度值设定为255，以减小该条连线对后续所要产生的连线的影响。
Step3：同理，以新得到的后续连接点pi作为出发点，重复Step2，得到该出发点的后续连接点pi+1，数值pi+1存入数组sequence。将肖像图像Portrait_area的pi和pi+1号锚点的连线所经过的所有像素点的灰度值设定为255，以减小该条连线对后续所要产生的连线的影响。重复执行996次。
通过以上步骤，得到数组sequence，该数组共存储1001个锚点编号，以黑色的连线，依次在半径为350的圆形空白画布上顺序连接数组sequence中的锚点，如图12所示。数组sequence中存储的锚点的数量并非越多越好，数量太少，不能较好的刻画人物肖像的细节；数量太多，则生成的线画肖像太过浓密，过多的线条会遮盖细节。本发明设定数组sequence存储的锚点的个数在1000左右，做到了非常好的均衡，整体上可以非常好的还原人物肖像，细节上也能很好的体现出来。
设定Step1～Step3循环的次数，分别为100、300、500和1500次，则分别得到不同大小的数组sequence，对应的就会产生图13(a)、(b)、(c)所示的分别为100、300、500和1500次顺序连线所构成的线肖像画。
步骤七：去除步骤六中初始线肖像画上的冗余，保留连线的交点，得到黑色线肖像画。
针对图12，提取该图的第2～699行、第2～699列所构成感兴趣区域，依次对该区域中的每一个像素点Pixel(m1,n1)的八邻域进行分析，将Pixel(m1,n1)的八个邻域的像素值求和得到邻域值Sneighbourhood，其中，m1的取值范围为2～699，n1的取值范围为2～699。如果像素点Pixel(m1,n1)的像素值为0，并且邻域值Sneighbourhood>4，则将Pixel(m1,n1)的像素值设定为1。对比图14和图12，明显发现图14在一定程度上减少了无关线段的干扰，使得肖像画更为精炼准确。
步骤八：检测步骤一中人物肖像图像的皮肤区域，循环步骤二-步骤七以肤色的连线顺序连接皮肤区域内的起始点和后续连接点，得到肤色线肖像画，将肤色线肖像画与步骤七得到的黑色线肖像画叠加得到彩色线肖像画。
检测原始的人物肖像图像中的皮肤区域；并以同样的方法针对皮肤区域生成起始点和199个后续连接点，并依次连线，其中连线的颜色为肤色，即生成表征皮肤区域的肤色连线的线肖像画。
将原始的彩色的人物肖像图像从RGB颜色空间转换为YCbCr颜色空间，其中，Y是指亮度分量，Cb指蓝色色度分量，而Cr指红色色度分量。设定归一化的转换矩阵HI为：
则亮度分量Y＝HI(1,1)*Rc+HI(1,2)*Gc+HI(1,3)*Bc+16，蓝色色度分量Cb＝HI(2,1)*Rc+HI(2,2)*Gc+HI(2,3)*Bc+128，红色色度分量Cr＝HI(3,1)*Rc+HI(3,2)*Gc+HI(3,3)*Bc+128，Rc、Gc、Bc分别表示原始的彩色的人物肖像图像的红、绿、蓝三个通道的像素值。通过以上方法计算得到的Y、Cb、Cr这三个分量中的像素的值，如果大于236，则将像素值修改为236；Y、Cb、Cr这三个分量中的像素的值，如果小于15，则将像素值修改为15，得到YCbCr图像。遍历YCbCr图像中的每个像素点，采用公式(3)计算各个像素点的肤色值skin为：
skin(ii,jj)＝(cos(2.54)*(Cb(ii,jj)–109.39)+sin(2.54)*(Cr(ij,jj)-152.03)–1.61)^2/(a^2)+
(-sin(2.54)*(Cb(ij,jj)-109.39)+cos(2.54)*(Cr(ij,jj)-152.03)–2.42)^2/(b^2)；(3)
其中，Cb(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的蓝色色度分量Cb的值；Cr(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的红色色度分量Cr的值；a、b为自定义的椭圆肤色模型的长轴、短轴，a＝25.39；b＝14.03。如果skin<＝1，则将第ii行、第jj列的像素点标记为皮肤点，可以用于表示肖像中的皮肤，所有皮肤点的集合为皮肤区域，皮肤区域如图15中的白色区域所示。
采用上述公式(2)、(3)的方法可计算并标记出原始彩色的人物肖像图像中的所有可以用于表示皮肤的像素点。将所有可用于表示皮肤的像素点的集合作为感兴趣区域，采用上述的生成黑色的线画肖像的方法，再次生成表示皮肤区域的线画肖像的起点及后续连接点，并以原始的人物彩色肖像图像中的皮肤区域的像素点的R通道的均值作为连线的R值，以原始的人物彩色肖像图像中的皮肤区域的像素点的G通道的均值作为连线的G值，以原始的人物彩色肖像图像中的皮肤区域的像素点的B通道的均值作为连线的B值，则以R值、G值、B值三个值为颜色的连线，依次连接表示皮肤区域的线画肖像的起点及后续连接点，即可构成肤色线肖像画，如图16所示，可以看出以肤色作为连线的颜色，多条连线很好的覆盖了皮肤区域，对还原人物肖像的肤色具有一定的帮助；将该肤色的线肖像画叠加到黑色线肖像画上，最终构成了彩色的线肖像画，如图17所示，可以看出黑色线肖像画对人物的头发、眉毛、眼睛、鼻子、嘴巴和脸部轮廓具有很好的表征，在此基础上叠加的肤色的线肖像画很好的还原了皮肤区域，最终生成的线肖像画具有一定的层次感，且细节还原度高。
本发明首先获取原始的人物彩色肖像图像，然后进行图像灰度化，生成包围肖像的圆形蒙版，在圆形蒙版的圆周上均匀的布设锚点，并计算各个锚点的坐标；计算线肖像画的起点，提出两个锚点的连线所经过的像素的查找方法；依次计算线肖像画的后续连接点，实现锚点的顺序存储，以黑色连线依次连接起点和后续连接点，生成黑色的线肖像画；去除黑色的线肖像画的冗余，保留交点，得到较为精细的黑色线肖像画；检测原始的人物彩色肖像图像中的皮肤区域，并以同样的方法针对皮肤区域生成起点和200个后续连接点，并依次连线，其中设定连线的颜色为肤色，肤色连线叠加到黑色线肖像画上，得到彩色线肖像画。本发明的方法具有自动化程度高、计算速度快和人物肖像还原精度高等有益效果，为人物肖像画生成提供了技术支撑，同时该研究思路与方法也可以推广应用到其它字符、标志等图案的生成，具有很高的实用价值。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于连续线条的人物肖像画生成方法，其特征在于，其步骤如下：
步骤一：通过相机采集或下载待处理的彩色的人物肖像图像；
步骤二：采用平均值法对步骤一中的人物肖像图像进行图像灰度化，得到灰度图像；
步骤三：对灰度图像进行调整，并生成与调整后灰度图像大小相同的包围肖像区域的圆形蒙版；
所述步骤三中对灰度图像进行调整的方法为：
以人物肖像区域为中心，以双线性插值的方式将灰度图像的大小调整为700×700像素；
所述圆形蒙版生成的方法为：
生成一个大小为700×700像素的画布图像，以画布图像的中心为原点，350个像素的宽度为半径在画布图像上画圆，并将圆形区域内的像素值设置为0，圆形区域外的像素值设置为1，得到包围肖像区域的圆形蒙版；
步骤四：在圆形蒙版的圆形轮廓上均匀设定多个锚点，并计算各个锚点的坐标；
步骤五：根据步骤三得到的圆形蒙版和灰度图像得到肖像图像，在肖像图像上的圆形轮廓上选取任意两个锚点并连线，计算连线经过的像素点的像素值的均值，将均值中的最小值对应的锚点作为起始点；
步骤六：将步骤五得到的起始点存入数组中，依次计算后续连接点并顺序存入数组中，在与肖像图像大小相同的圆形空白画布上，以黑色的连线将数组中的起始点和后续连接点顺序连接，得到初始线肖像画；
步骤七：去除步骤六中初始线肖像画上的冗余，保留连线的交点，得到黑色线肖像画；
步骤八：检测步骤一中人物肖像图像的皮肤区域并作为感兴趣区域，循环步骤二-步骤七以肤色的连线顺序连接表征皮肤区域的起始点和后续连接点，得到肤色线肖像画，将肤色线肖像画与步骤七得到的黑色线肖像画叠加得到彩色线肖像画。
2.根据权利要求1所述的基于连续线条的人物肖像画生成方法，其特征在于，所述步骤二中平均值法的实现方法为：提取人物肖像图像的红绿蓝三个通道的平均值作为灰度值：
Gray＝(Rc+Gc+Bc)/3；
其中，Gray表示灰度图像的像素值，Rc、Gc、Bc分别表示人物肖像图像的红、绿、蓝三个通道的像素值。
3.根据权利要求1所述的基于连续线条的人物肖像画生成方法，其特征在于，所述步骤四的实现方法为：在圆形蒙版的圆形轮廓的圆周上每隔1.8°均匀的设定200个锚点，以圆形轮廓的顶部最高位置处的点为初始锚点，并按顺时针方向依次计算其他锚点的坐标，依次存储到集合pinsXY(i1,j1)中，其中，i1表示锚点的编号，取值范围为1～200；j1表示锚点的x坐标或者y坐标，取值范围为1～2，当j1＝1时表示x坐标，当j1＝2时表示y坐标。
4.根据权利要求1或3所述的基于连续线条的人物肖像画生成方法，其特征在于，所述肖像图像是将圆形蒙版中的值取非，然后将圆形蒙版与灰度图像进行点乘得到的；
所述计算连线经过的像素点的像素值的均值的方法为：
从半径为350的圆形轮廓上的200个锚点中任意挑选两个锚点，将这两个锚点连线，统计连线所经过的肖像图像的像素点的数量Numij及像素点的像素值之和Sumij，其中，i的取值为1～199，j的取值为(i+1)～200；如果数量Numij>30，则计算该连线所经过的像素点的像素值的均值Aveij＝Sumij/Numij；否者j＝j+1；统计连线所经过像素点的像素值的均值的集合Ave＝{Aveij}中的最小值Avemn，其中，m∈i，n∈j，且m<n，则选取锚点m作为线肖像画的起始点。
5.根据权利要求4所述的基于连续线条的人物肖像画生成方法，其特征在于，所述统计连线所经过的肖像图像的像素点的数量Numij及像素点的像素值之和Sumij时，计算任意两个锚点A和锚点B的连线所经过肖像图像的像素点的坐标的方法为：锚点A的坐标为(pinsXY(i,1)，pinsXY(i,2))，锚点B的坐标为(pinsXY(j,1)，pinsXY(j,2))，锚点A和锚点B之间连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)；则：
(1)当锚点A的y坐标yA＝pinsXY(i,2)与锚点B的y坐标yB＝pinsXY(j,2)相等，则锚点A和锚点B处于同一行，判断锚点A和锚点B的x坐标值：
若xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA-xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA-xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1-1、y坐标的值为AB_xy(mdsum1,2)＝yA；
若xB>xA，则锚点A和锚点B之间连线经过的像素点的数量为(xB-xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB-xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1-1、y坐标的值为AB_xy(mdsum1,2)＝yA；
(2)锚点A的x坐标xA＝pinsXY(i,1)与锚点B的x坐标xB＝pinsXY(j,1)相等，则锚点A和锚点B处于同一列，判断锚点A和锚点B的y坐标值：
若yA>yB，则锚点A和锚点B之间连线经过的像素点的数量为(yA-yB+1)个，连线经过的像素点的坐标的集合AB_xy中，索引mdsum1的取值为1～(yA-yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1-1、x坐标的值为AB_xy(mdsum1,1)＝xA；
若yB>yA，则锚点A和锚点B之间连线经过的像素点的数量为(yB-yA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yB-yA+1)，索引mdsum2的取值为1～2；连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1-1、x坐标值为AB_xy(mdsum1,1)＝xA；
(3)锚点A和锚点B在x方向的偏移xOffset＝|xA–xB|，锚点A和锚点B在y方向的偏移yOffset＝|yA–yB|，若xOffset≥yOffset，锚点A和锚点B的连线的斜率的绝对值小于等于1，则：
若yA>yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA–xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA–yB)/(xA–xB)|)，其中，round()为四舍五入取整函数；
若yA>yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB–xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yB+mdsum1×|(yA yB)/(xA–xB)|)；
若yA<yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(xA–xB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xA–xB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xB+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)；
若yA<yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(xB–xA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(xB–xA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝xA+mdsum1–1、y坐标的值为AB_xy(mdsum1,2)＝round(yA+mdsum1×|(yA–yB)/(xA–xB)|)；
(4)锚点A和锚点B在x方向的偏移为xOffset＝|xA–xB|，在y方向的偏移为yOffset＝|yA–yB|，若xOffset<yOffset，锚点A和锚点B的连线的斜率的绝对值大于1，则：
若yA>yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yA–yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yB+mdsum1–1；
若yA>yB且xA<xB，则锚点A和锚点B之间连线经过的像素点的数量为(yA–yB+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yA–yB+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标值为AB_xy(mdsum1,2)＝yB+mdsum1–1；
若yA<yB且xA>xB，则锚点A和锚点B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合AB_xy(mdsum1,mdsum2)中，索引mdsum1的取值为1～(yB–yA+1)，索引mdsum2的取值为1～2；连线上各个像素点的x坐标值为AB_xy(mdsum1,1)＝round(xB+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标值为AB_xy(mdsum1,2)＝yA+mdsum1–1；
若yA<yB且xA<xB，则两个锚点A、B之间连线经过的像素点的数量为(yB–yA+1)个，连线经过的像素点的坐标的集合为AB_xy(mdsum1,mdsum2)，其中，mdsum1的取值为1～(yB–yA+1)，mdsum2的取值为1～2；连线上各个像素点的x坐标的值为AB_xy(mdsum1,1)＝round(xA+mdsum1×|(xA–xB)/(yA–yB)|)，连线上各个像素点的y坐标的值为AB_xy(mdsum1,2)＝yA+mdsum1–1；
根据连线上各个像素点的x坐标和y坐标的值，通过Pr(x,y)提取像素点的像素值即灰度值，将锚点A和锚点B之间连线经过的像素点的灰度值相加得到像素值之和Sumij。
6.根据权利要求5所述的基于连续线条的人物肖像画生成方法，其特征在于，所述步骤六的实现方法为：
Step1：通过步骤五的处理，得到锚点m为线画肖像的起始点，锚点n是后续连接点，且锚点m和锚点n的连线经过的所有像素点的灰度值的均值最小；将数值m存入数组sequence，将肖像图像的锚点m和锚点n的连线经过的所有像素点的灰度值设定为255；
Step2：以锚点n为出发点，分别计算锚点n到锚点1～锚点(m–1)、锚点(m+1)～锚点(n–1)号及锚点(n+1)～锚点200的连线经过的像素的个数，并排除其中连线所经过像素数少于30的锚点；然后采用步骤五的方法计算得到锚点pi是后续连接点，使得锚点n和锚点pi的连线经过的所有像素的像素值的均值最小；将数值n存入数组sequence中，将肖像图像的锚点n和锚点pi的连线经过的所有像素点的灰度值设定为255；
Step3：同理，以得到的后续连接点pi作为出发点，循环Step2，得到该出发点的后续连接点pi+1，数值pi+1存入数组sequence；将肖像图像的pi和pi+1号锚点的连线所经过的所有像素点的灰度值设定为255；
Step4：以黑色的连线依次在半径为350的圆形空白画布上顺序连接数组sequence中的锚点，得到初始线肖像画。
7.根据权利要求6所述的基于连续线条的人物肖像画生成方法，其特征在于，所述步骤七中黑色线肖像画的获取方法为：提取初始线肖像画的第2～699行、第2～699列构成感兴趣区域，依次对感兴趣区域中的每一个像素点Pixel(m1,n1)的八邻域进行分析，将像素点Pixel(m1,n1)的八个邻域的像素值求和得到邻域值Sneighbourhood，其中，m1的取值范围为2～699，n1的取值范围为2～699；如果像素点Pixel(m1,n1)的像素值为0，并且邻域值Sneighbourhood>4，则将像素点Pixel(m1,n1)的像素值设定为1，得到黑色线肖像画。
8.根据权利要求1或7所述的基于连续线条的人物肖像画生成方法，其特征在于，所述步骤八中检测人物肖像图像的皮肤区域的方法为：
将人物肖像图像从RGB颜色空间转换为YCbCr颜色空间，其中，Y是指亮度分量，Cb指蓝色色度分量，而Cr指红色色度分量，且归一化的转换矩阵HI为：
则亮度分量Y＝HI(1,1)*Rc+HI(1,2)*Gc+HI(1,3)*Bc+16，蓝色色度分量Cb＝HI(2,1)*Rc+HI(2,2)*Gc+HI(2,3)*Bc+128，红色色度分量Cr＝HI(3,1)*Rc+HI(3,2)*Gc+HI(3,3)*Bc+128，Rc、Gc、Bc分别表示原始的彩色的人物肖像图像的红、绿、蓝三个通道的像素值；
若亮度分量Y、蓝色色度分量Cb或红色色度分量Cr的像素的值大于236，则将像素值修改为236；若亮度分量Y、蓝色色度分量Cb或红色色度分量Cr的像素的值小于15，则将像素值修改为15，得到YCbCr图像；
遍历YCbCr图像中的每个像素点，计算各个像素点的肤色值skin为：
skin(ii,jj)＝(cos(2.54)*(Cb(ii,jj)–109.39)+sin(2.54)*(Cr(ij,jj)-152.03)–1.61)^2/(a^2)+(-sin(2.54)*(Cb(ij,jj)-109.39)+cos(2.54)*(Cr(ij,jj)-152.03)–2.42)^2/(b^2)；
其中，Cb(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的蓝色色度分量Cb的值；Cr(ii,jj)表示YCbCr图像中的第ii行、第jj列的像素点的红色色度分量Cr的值；a、b为自定义的椭圆肤色模型的长轴、短轴，且a＝25.39；b＝14.03；
如果skin<＝1，则将第ii行、第jj列的像素点标记为皮肤点，所有皮肤点的集合为皮肤区域。
9.根据权利要求8所述的基于连续线条的人物肖像画生成方法，其特征在于，所述肤色线肖像画的实现方法为：采用步骤二-步骤七中生成黑色线肖像画的方法，再次生成表示皮肤区域的线画肖像的起点及后续连接点，并以人物肖像图像中的皮肤区域的像素点的R通道的均值作为连线的R值，以人物肖像图像中的皮肤区域的像素点的G通道的均值作为连线的G值，以人物肖像图像中的皮肤区域的像素点的B通道的均值作为连线的B值，以三个通道的R值、G值、B值为颜色的连线，依次连接表示皮肤区域的线画肖像的起点及后续连接点，构成肤色线肖像画。
CN202110143469.0A 2021-02-02 2021-02-02 一种基于连续线条的人物肖像画生成方法 Active CN112819922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110143469.0A CN112819922B (zh) 2021-02-02 2021-02-02 一种基于连续线条的人物肖像画生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110143469.0A CN112819922B (zh) 2021-02-02 2021-02-02 一种基于连续线条的人物肖像画生成方法

Publications (2)

Publication Number Publication Date
CN112819922A CN112819922A (zh) 2021-05-18
CN112819922B true CN112819922B (zh) 2024-03-12

Family

ID=75860625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110143469.0A Active CN112819922B (zh) 2021-02-02 2021-02-02 一种基于连续线条的人物肖像画生成方法

Country Status (1)

Country Link
CN (1) CN112819922B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1818869A1 (en) * 2004-11-25 2007-08-15 NEC Corporation Face image synthesis method and face image synthesis device
CN104200505A (zh) * 2014-08-27 2014-12-10 西安理工大学 一种人脸视频图像的漫画式动画生成方法
CN106127841A (zh) * 2016-06-22 2016-11-16 丁焱 一种基于人脸照片生成个性卡通动态图的方法
CN107730573A (zh) * 2017-09-22 2018-02-23 西安交通大学 一种基于特征提取的人物肖像漫画风格化生成方法
DE102017011132A1 (de) * 2017-12-01 2019-06-06 Mühlbauer Gmbh & Co. Kg Verfahren zur Herstellung eines Personenporträts für ein Ausweisdokument

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1818869A1 (en) * 2004-11-25 2007-08-15 NEC Corporation Face image synthesis method and face image synthesis device
CN104200505A (zh) * 2014-08-27 2014-12-10 西安理工大学 一种人脸视频图像的漫画式动画生成方法
CN106127841A (zh) * 2016-06-22 2016-11-16 丁焱 一种基于人脸照片生成个性卡通动态图的方法
CN107730573A (zh) * 2017-09-22 2018-02-23 西安交通大学 一种基于特征提取的人物肖像漫画风格化生成方法
DE102017011132A1 (de) * 2017-12-01 2019-06-06 Mühlbauer Gmbh & Co. Kg Verfahren zur Herstellung eines Personenporträts für ein Ausweisdokument

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖像漫画绘制机器人技术研究;倪菲;付庄;曹其新;赵言正;;自然杂志;20070815(04);全文 *

Also Published As

Publication number Publication date
CN112819922A (zh) 2021-05-18

Similar Documents

Publication Publication Date Title
Žeger et al. Grayscale image colorization methods: Overview and evaluation
US9449253B2 (en) Learning painting styles for painterly rendering
CN104282002B (zh) 一种数字图像的快速美容方法
US20110292051A1 (en) Automatic Avatar Creation
CN112950775A (zh) 一种基于自监督学习的三维人脸模型重建方法及系统
JP2010507854A (ja) ビデオイメージシーケンスを仮想シミュレーションする方法及び装置
KR20230097157A (ko) 개인화된 3d 헤드 모델 변형을 위한 방법 및 시스템
CN113723385B (zh) 视频处理方法及装置、神经网络的训练方法及装置
JP7462120B2 (ja) ２次元（２ｄ）顔画像から色を抽出するための方法、システム及びコンピュータプログラム
KR102353556B1 (ko) 사용자 얼굴기반 표정 및 포즈 재현 아바타 생성장치
CN112995534B (zh) 视频生成方法、装置、设备及可读存储介质
CN116997933A (zh) 用于构造面部位置图的方法和系统
CN112221145B (zh) 游戏脸模生成方法及装置、存储介质及电子设备
CN117157673A (zh) 用于形成个性化的3d头部和面部模型的方法和系统
CN111243051A (zh) 基于肖像照片的简笔画生成方法、系统及存储介质
CN110992247A (zh) 一种人像照片头发拉直特效的实现方法及系统
Zhao et al. Research on the application of computer image processing technology in painting creation
Li et al. Line drawing guided progressive inpainting of mural damages
Tresset et al. Generative portrait sketching
CN112819922B (zh) 一种基于连续线条的人物肖像画生成方法
Luo et al. Robot artist performs cartoon style facial portrait painting
KR20010084996A (ko) 단일 이미지를 이용한 3차원 아바타 제작 방법 및 이를이용한 자판기
Ye et al. Hybrid scheme of image’s regional colorization using mask r-cnn and Poisson editing
Laishram et al. High-quality face caricature via style translation
Xue et al. User-guided chinese painting completion–a generative adversarial network approach

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant