本文用 Python 实现 PS 的一种滤镜效果,称为万花筒。也是对图像做各种扭曲变换,最后图像呈现的效果就像从万花筒中看到的一样:

图像的效果可以参考之前的博客:

http://blog.csdn.net/matrix_space/article/details/46789783

import matplotlib.pyplot as plt
from skimage import io
from skimage import img_as_float
import numpy as np
import numpy.matlib
import math file_name='D:/Visual Effects/PS Algorithm/4.jpg';
img=io.imread(file_name) img = img_as_float(img) row, col, channel = img.shape # set the parameters
radius = 100.0
angle = math.pi/3
angle2 = math.pi/4
sides = 10.0
# set the center of the circle, proportion of the image size
centerX = 0.5
centerY = 0.5 iWidth=col
iHeight=row center_x=iWidth*centerX
center_y=iHeight*centerY xx = np.arange (col)
yy = np.arange (row) x_mask = numpy.matlib.repmat (xx, row, 1)
y_mask = numpy.matlib.repmat (yy, col, 1)
y_mask = np.transpose(y_mask) xx_dif = x_mask - center_x
yy_dif = y_mask - center_y r = np.sqrt(xx_dif * xx_dif + yy_dif * yy_dif)
theta = np.arctan2(yy_dif, xx_dif+0.0001) - angle - angle2
temp_theta=theta/math.pi*sides*0.5
temp_r = np.mod(temp_theta, 1.0)
mask_1 = temp_r < 0.5
theta = temp_r * 2 * mask_1 + (1-temp_r) * 2 * (1 - mask_1) radius_c=radius/np.cos(theta) temp_r = np.mod (r/radius_c, 1.0)
mask_1 = temp_r < 0.5
r = radius_c * (temp_r * 2 * mask_1 + (1-temp_r) * 2 * (1 - mask_1)) theta = theta + angle x1_mask = r * np.cos(theta) + center_x
y1_mask = r * np.sin(theta) + center_y mask = x1_mask < 0
x1_mask = x1_mask * (1 - mask)
mask = x1_mask > (col - 1)
x1_mask = x1_mask * (1 - mask) + (x1_mask * 0 + col -2) * mask mask = y1_mask < 0
y1_mask = y1_mask * (1 - mask)
mask = y1_mask > (row -1)
y1_mask = y1_mask * (1 - mask) + (y1_mask * 0 + row -2) * mask img_out = img * 1.0 int_x = np.floor (x1_mask)
int_x = int_x.astype(int)
int_y = np.floor (y1_mask)
int_y = int_y.astype(int) p_mask = x1_mask - int_x
q_mask = y1_mask - int_y img_out = img * 1.0 for ii in range(row):
for jj in range (col):
new_xx = int_x [ii, jj]
new_yy = int_y [ii, jj]
# p = p_mask[ii, jj]
# q = q_mask[ii, jj] img_out[ii, jj, :] = img[new_yy, new_xx, :] plt.figure (1)
plt.imshow (img)
plt.axis('off') plt.figure (2)
plt.imshow (img_out)
plt.axis('off') plt.show()

最新文章

  1. [PHP源码阅读]array_push和array_unshift函数
  2. Android中Service通信(二)——绑定Service进行通信
  3. Codeforces617 E . XOR and Favorite Number(莫队算法)
  4. java考核完的心得
  5. Android Studio快捷键每日一练(6)
  6. JavaScript中call,apply,bind方法的总结。
  7. Harris角点检测算法优化
  8. Win7家庭版开启Administrator管理员帐户的方法
  9. 转载---SQL Server XML基础学习&lt;1&gt;之--FOR XML PATH
  10. HTML5性能优化[转]
  11. redis五种基本类型CRUD操作
  12. 关于HTML5中的sessionStorage和localStorage
  13. 简述DDOS攻击的工作原理
  14. Linux五种IO模型(同步 阻塞概念)
  15. Qt ------ QTabWidget
  16. 未能加载文件或程序集“file:///D:/Program Files (x86)/ArcGIS/DeveloperKit10.0/DotNet/ESRI.ArcGIS.3DAnalyst.dll”或它的某一个依赖项。试图加载格式不正确的程序。 行 129,位置 5。
  17. 查询oracle安装过补丁没有
  18. Nodejs Mocha测试学习
  19. SharePoint中遇到Timeout
  20. linux防火墙过滤规则

热门文章

  1. 使用TensorRT加速yolo3
  2. js中给正则传参、传递变量
  3. assert.deepStrictEqual()
  4. LeetCode 121. Best Time to Buy and Sell Stock (stock problem)
  5. javascript倒计时代码及倒计时弹窗
  6. Codeforces Round #354 (Div. 2)-C. Vasya and String,区间dp问题,好几次cf都有这种题,看来的好好学学;
  7. parse XML &amp; js
  8. 2018/2/15 ES Beats的学习笔记
  9. msp430入门编程03
  10. linux 用户管理、权限管理