Perlin噪声-526互联

title: Perlin噪声
date: 2023-04-18 11:10:14
tags:
    - 噪声
    - perlin
categories:
    - 充电学习

perlin基本信息

Perlin噪声（Perlin noise，又称为柏林噪声）指由Ken Perlin发明的自然噪声生成算法，具有在函数上的连续性，并可在多次调用时给出一致的数值。在电子游戏领域中可以透过使用Perlin噪声生成具连续性的地形；或是在艺术领域中使用Perlin噪声生成图样。

由于一些历史原因，Simplex噪声和分形噪声（texture synthesis）都曾在学术论文中被单独称作Perlin噪声。

经典Perlin噪声

perline噪声是基于晶格的方法。它属于梯度噪声，其原理就是将坐标系划分成一块一块的晶格，之后在晶格的顶点出生成一个随机梯度，通过与晶格顶点到晶格内点的向量进行点乘加权计算后得到噪声。

### 实现perlin噪声的过程需要插值 * 对于一维：插值使用的是一个在0处为1，在1处为0，在0.5处为0.5的连续单调递减函数。例如对，设$c_{0}$,$c_{1}$为左右两个整数点的数值，t为该点距离左边点的距离，使用$(1-t)$作为插值函数，则该点的值为$c_{1}(1-t)+c_{0}t$。 ![](https://img2023.cnblogs.com/blog/1218402/202304/1218402-20230424181648148-382393529.jpg)

但是$(1-t)$是线性插值，人工痕迹比较严重，并且在整数点上不连续。Perlin建议使用$3t^{2}-2t$作为插值函数。后来建议使用$6t^{5}-15t+10t^{{3}$作为插值函数。事实上，只有在区间[0,1]内的连续函数$f$，有$f(0)=1,f(1)=0$且$f}{}(0)=f^{}(1)=0$的函数皆可作为插值函数。

对于二维：对于点$(x,y)$，令$i=\lfloor x \rfloor,j=\lfloor y \rfloor$，它所在的晶格的四个顶点分别为$(i,j)、(i+1,j)、(i+1,j+1)、(i,j+1)$。令$u=x-i,v=y-j$，这四个顶点对点$(x,y)$的贡献可以使用它们的梯度$(g_{00},g_{10},g_{11},g_{01})$和$(x,y)$点与这四个顶点的方向$((u,v),(u-1,v),(u-1,v-1),(u,v-1))$进行点积获得。但是在二维的情况下，插值更为复杂。首先需要对 $(i,j)$ 和$(i+1,j)$ 两点在$x$方向插值，得到点$(x,j)$的值；之后对$(i,j+1)$ 和 $(i+1,j+1)$两点在$x$方向插值，得到点$(x,j+1)$的值；最后对$(x,j)$和$(x,j+1)$ 在$y$方向插值，得到$(x, y)$的值。

## code [完整代码](https://github.com/DreamOneYou/Perline) ```python def interpolant(t): return t*t*t*(t*(t*6 - 15) + 10)

def generate_perlin_noise_2d(
shape, res, tileable=(False, False), interpolant=interpolant
):
delta = (res[0] / shape[0], res[1] / shape[1])
d = (shape[0] // res[0], shape[1] // res[1])
grid = np.mgrid[0:res[0]:delta[0], 0:res[1]:delta[1]]
.transpose(1, 2, 0) % 1
# Gradients
angles = 2np.pinp.random.rand(res[0]+1, res[1]+1)
gradients = np.dstack((np.cos(angles), np.sin(angles)))
if tileable[0]:
gradients[-1,:] = gradients[0,:]
if tileable[1]:
gradients[:,-1] = gradients[:,0]
gradients = gradients.repeat(d[0], 0).repeat(d[1], 1)
g00 = gradients[ :-d[0], :-d[1]]
g10 = gradients[d[0]: , :-d[1]]
g01 = gradients[ :-d[0],d[1]: ]
g11 = gradients[d[0]: ,d[1]: ]
# Ramps
n00 = np.sum(np.dstack((grid[:,:,0] , grid[:,:,1] )) * g00, 2)
n10 = np.sum(np.dstack((grid[:,:,0]-1, grid[:,:,1] )) * g10, 2)
n01 = np.sum(np.dstack((grid[:,:,0] , grid[:,:,1]-1)) * g01, 2)
n11 = np.sum(np.dstack((grid[:,:,0]-1, grid[:,:,1]-1)) * g11, 2)
# Interpolation
t = interpolant(grid)
c0 = n00(1-t[:,:,0]) + t[:,:,0]n10
c1 = n01(1-t[:,:,0]) + t[:,:,0]n11
return np.sqrt(2)((1-t[:,:,1])c0 + t[:,:,1]*c1)


## 结果
* Perlin 2D噪声 ![](/2023/Perlin噪声/generate_fractal_noise_2d.png)
* 3D perline噪声

![](https://img2023.cnblogs.com/blog/1218402/202304/1218402-20230424181928096-1537958653.gif)

* 给人物图像添加perlin噪声

<table><tr>
<td><img src="https://img2023.cnblogs.com/blog/1218402/202304/1218402-20230424181758778-370438299.jpg" border=0 width=400px height=400px></td>
<td><img src="https://img2023.cnblogs.com/blog/1218402/202304/1218402-20230424181808567-521996403.png" border=0 width=400px height=400px></td>
</tr></table>