网库网(www.wangkuwang.com)精品网站源码,织梦建站模版,游戏源代码分享平台

模板户源码

当前位置:首页 -> CMS教程 -> 其它教程 正文

只需45秒,Python给故宫画一组手绘图!

时间:2019-03-18 07:55:28 [整站源码]作者:zhaopulei


丁彦军 CDA数据分析师

作者 丁彦军

来源 | 恋习Python(ID:sldata2017


13日早晨,当北京市民拉开窗帘时发现,窗外雪花纷纷扬扬在空中飘落,而且越下越大,树上、草地、屋顶、道路上,都落满雪花。京城银装素裹,这是今冬以来北京迎来的第二场降雪。

一下雪,北京就变成了北平,故宫就变成了紫禁城。八万张门票在雪花飘下来之前,便早已预订一空。

(图片来源:故宫官网 版权归故宫官网所有)



看着朋友圈、微博好友都在纷纷晒图,小编只能羡慕不已。


不过,突然想到,可以通过Python将故宫的建筑物图片,转化为手绘图(素描效果)。效果图如下:



一、概念与原理


我们都知道手绘图效果的特征主要有:


  • 黑白灰色;边界线条较重;相同或相近色彩趋于白色;略有光源效果


核心原理:利用像素之间的梯度值和虚拟深度值对图像进行重构,根据灰度变化来模拟人类视觉的模拟程度


把图像看成二维离散函数,灰度梯度其实就是这个二维离散函数的求导,用差分代替微分,求取图像的灰度梯度。常用的一些灰度梯度模板有:Roberts 梯度、Sobel 梯度、Prewitt 梯度、Laplacian 梯度。


以Sobel 梯度计算来解释:


首先计算出 ,然后计算梯度角 梯度方向及图像灰度增大的方向,其中梯度方向的梯度夹角大于平坦区域的梯度夹角。如下图所示,灰度值增加的方向梯度夹角大,此时梯度夹角大的方向为梯度方向。对应在图像中寻找某一点的梯度方向即通过计算该点与其8邻域点的梯度角,梯度角最大即为梯度方向。



二、图像的数组形式与变换



其中,需要用到的方法:


  • Image.open( ): 打开图片

  • np.array( ) : 将图像转化为数组

  • convert("L"): 将图片转换成二维灰度图片

  • Image.fromarray( ): 将数组还原成图像uint8格式


代码如下:


from PIL import Image
import numpy as np

im = Image.open(r"C:UsersAdministratorDesktopgugong微信图片_20190216152248.jpg").convert(L)
a=np.asarray(im).astype(float)
print(a.shape,a.dtype)
(1080608) float64
#(1080, 608)分别表示高度,宽度


三、图像的手绘效果处理


实现思路步骤:


1、梯度的重构


numpy的梯度函数的介绍


np.gradient(a) : 计算数组a中元素的梯度,f为多维时,返回每个维度的梯度 

离散梯度: xy坐标轴连续三个x轴坐标对应的y轴值:a, b, c 其中b的梯度是(c-a)/2 


而c的梯度是: (c-b)/1


当为二维数组时,np.gradient(a) 得出两个数组,第一个数组对应最外层维度的梯度,第二个数组对应第二层维度的梯度。 


代码如下:


grad=np.gradient(a)
grad_x,grad_y=grad
grad_x = grad_x * depth / 100.#对grad_x值进行归一化
grad_y = grad_y * depth / 100.#对grad_y值进行归一化


2、构造guan光源效果


设计一个位于图像斜上方的虚拟光源
光源相对于图像的视角为Elevation,方位角为Azimuth
建立光源对各点梯度值的影响函数
运算出各点的新像素值



其中:

np.cos(evc.el) : 单位光线在地平面上的投射长度

dx,dy,dz :光源对x,y,z三方向的影响程度


3、梯度归一化


  • 构造x和y轴梯度的三维归一化单位坐标系;

  • 梯度与光源相互作用,将梯度转化为灰度。


4、图像生成


具体详情代码如下:


from PIL import Image
import numpy as np
import os
import join
import time

def image(sta,end,depths=10):
    a = np.asarray(Image.open(sta).convert(L)).astype(float)
    depth = depths  # 深度的取值范围(0-100),标准取10
    grad = np.gradient(a)  # 取图像灰度的梯度值
    grad_x, grad_y = grad  # 分别取横纵图像梯度值
    grad_x = grad_x * depth / 100.#对grad_x值进行归一化
    grad_y = grad_y * depth / 100.#对grad_y值进行归一化
    A = np.sqrt(grad_x ** 2 + grad_y ** 2 + 1.)
    uni_x = grad_x / A
    uni_y = grad_y / A
    uni_z = 1. / A
    vec_el = np.pi / 2.2  # 光源的俯视角度,弧度值
    vec_az = np.pi / 4.  # 光源的方位角度,弧度值
    dx = np.cos(vec_el) * np.cos(vec_az)  # 光源对x 轴的影响
    dy = np.cos(vec_el) * np.sin(vec_az)  # 光源对y 轴的影响
    dz = np.sin(vec_el)  # 光源对z 轴的影响
    b = 255 * (dx * uni_x + dy * uni_y + dz * uni_z)  # 光源归一化
    b = b.clip(0255)
    im = Image.fromarray(b.astype(uint8))  # 重构图像
    im.save(end)

def main():
    xs=10
    start_time = time.clock()
    startss = os.listdir(r"C:UsersAdministratorDesktopgugong")
    time.sleep(2)
    for starts in startss:
        start = .join(starts)
        sta = C:/Users/Administrator/Desktop/gugong/ + start
        end = C:/Users/Administrator/Desktop/gugong/ + HD_ + start
        image(sta=sta,end=end,depths=xs)

    end_time = time.clock()
    print(程序运行了  ---- + str(end_time - start_time) +    秒)
    time.sleep(3)

main()
程序运行了  ----43.01828205879955   秒  #一共35张图片


最终效果图对比:



最后,你自己动手试试吧?通过此代码为自己画一张手绘图,也可以为自己的家乡或母校画。


参考资料:

http://www.icourse163.org/learn/BIT-1001870002?tid=1001963001#/learn/announce

故宫雪景图片素材链接:

https://pan.baidu.com/s/1E_aZTRQWOzGV-2GV_iH43w

提取码:64z9



CDA 课程咨询丨赵老师

联系电话:13381275813
扫描二维码





更多精彩文章

免费!10本必读的机器学习书籍(附下载)

Python 告诉你绝不知道的1983-2018 春晚

第九届CDA考试状元访谈 | 看看状元们都说了啥

猫眼5万条评论告诉你《飞驰人生》值不值得看?!

盘点 | 过去一年,AI 圈都发生了哪些事儿(附年度报告下载)

Python程序员!变身!这有一张最高效成长路线图(附资源)


文章转载自公众号

恋习Python 恋习Python
文章已于修改

    发送中

    本文标签:AutoTags插件服务端需要您提供购买者的账号和密码才能继续访问  折翼天使  莎莎源码  吾爱源码  其他源码 

    转载请注明来源:PHP手机端发卡多种支付商业版源码

    本文永久链接地址:https://www.suibianlu.com/11942.html

    郑重声明:
    本站所有内容均由互联网收集整理、网友上传,并且以计算机技术研究交流为目的,仅供大家参考、学习,不存在任何商业目的与商业用途。
    若您需要商业运营或用于其他商业活动,请您购买正版授权并合法使用。 我们不承担任何技术及版权问题,且不对任何资源负法律责任。
    如无法链接失效或侵犯版权,请给我们来信:admin@suibianlu.com

    栏目导航
    最新文章
    热门文章
    Top