当前位置：首页 > 编程日记 > 正文

python 归一化_只需 45 秒，Python 给故宫画一组手绘图！

编程日记 2024-05-16 11:00:00

作者 | 丁彦军

责编 | 伍杏玲

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

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

看着朋友圈、微博好友都在纷纷晒图，作为程序员的我只能坐在办公室里羡慕不已。

不过普通人得去故宫赏雪，程序员可以直接坐在家里动动手指便可以欣赏到故宫的美丽雪景手绘图了。如何操作呢？可通过Python将故宫的建筑物图片，转化为手绘图(素描效果)。效果图如下：

概念与原理

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

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

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

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

以Sobel 梯度计算来解释：

首先计算出

、

，然后计算梯度角

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

图像的数组形式与变换

其中，需要用到的方法：

Image.open( ): 打开图片
np.array( ) : 将图像转化为数组
convert("L"): 将图片转换成二维灰度图片
Image.fromarray( ): 将数组还原成图像uint8格式

代码如下：

from PIL import Imageimport numpy as npim = Image.open(r"C:甥敳獲AdministratorDesktopgugong微信图片_20190216152248.jpg").convert('L')a=np.asarray(im).astype('float')print(a.shape,a.dtype)(1080, 608) 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=gradgrad_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 ：光源对xyz三方向的影响程度。

3、梯度归一化

构造x和y轴梯度的三维归一化单位坐标系；
梯度与光源相互作用，将梯度转化为灰度。

4、图像生成

具体详情代码如下：

from PIL import Imageimport numpy as npimport osimport joinimport timedef 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(0, 255) im = Image.fromarray(b.astype('uint8')) # 重构图像 im.save(end)def main(): xs=10 start_time = time.clock() startss = os.listdir(r"C:甥敳獲AdministratorDesktopgugong") 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张图片

最终效果图对比：

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

作者：丁彦军，一名痴恋于Python的码农
公众号：恋习Python，在这里我们一起用Python做些有意义的事。
参考资料：
http://www.icourse163.org/learn/BIT-1001870002?tid=1001963001#/learn/announce
代码链接：https://pan.baidu.com/s/1E_aZTRQWOzGV-2GV_iH43w
提取码：64z9
声明：本文为作者投稿，版权归其个人所有。

https://www.dkcj.cn/info/6484.html

python 归一化_只需 45 秒，Python 给故宫画一组手绘图！

相关文章：

Windows平台下程序打包流程

一个apk多个ICON执行入口

vbs之CurrentDirectory

详解javascript: void(0);

读书笔记：编写高质量代码--web前端开发修炼之道（二：5章）

# Qt程序版本号使用

python traceback安装_Python错误：Traceback (most recent call last):感觉是软件出了问题

Directx11教程(6) 画一个简单的三角形(2)

【单调栈前缀和异或】7.21序列求和

hibernate 复合主键根据主键删除_hibernate封装Utils工具类

Linux+Qt 下同一数据空间vfork多进程间通信的一种高效便捷方式（信号槽直接调用）

Eclipse 调试器(引用IT168)

Cisco交换机与路由器的密码恢复_路由交换

投影转换_即插即用，办公投影不用愁：毕亚兹Mini DP转HDMIVGA转换器

事件绑定在IE下this是window的问题

nvidia-jetson系列硬件平台上安装Qt

以后在这里安家

如何用git命令行上传本地代码到github

C# 的快捷键汇总（一）

git修改远程仓库地址

ftp主动和被动模式_【扫盲】FTP基础知识详解

单页面与多页面的区别及优缺点

git使用手册

python 乒乓球_python乒乓球

《深入浅出Windows Phone 8应用开发》

CentOS7.4到Elasticsearch一路坑（五）

轻量级git服务器 Gogs git 服务器搭建

akaze特征匹配怎么去掉不合适的点_SIFT特征点

fopen 中按文本读写与按二进制读写实例

狼奔代码生成工具使用心得