python - 如何生成粉红噪声图像?
问题描述
我正在尝试用 Python 复制P. Bourke的“景观频率合成”。我以为这将是一个简单的
import numpy as np
from scipy.fft import fft2, ifft2
whitenoise = np.random.uniform(0,1,(256,256,3))
fouriertransformed = np.fft.fftshift(fft2(whitenoise))
pinktransformed = np.reciprocal(fouriertransformed)
pinknoise = ifft2(np.fft.ifftshift(pinktransformed)).real
但这似乎要复杂得多。我怎样才能做到这一点,如何检查结果图像中的功率是否真的下降了 1/f**2?
解决方案
这里的问题是,通过计算,您可以计算振幅pinktransformed = np.reciprocal(fouriertransformed)
的倒数,但您真正想要的是按比例缩放这些振幅,因此您必须将该行替换为1/f**2
pinktransformed = fouriertransformed / f**2
其中f
是一个数组,其中包含与傅立叶变换的每个 bin 对应的频率。
推荐阅读
- javascript - Nodejs Mongoose:尝试在 UTC 时区保存日期,但在本地时区保存
- drupal-7 - 与 drupal 网络表单的 Plupload 集成
- css - 按钮不会与列表项目符号保持一致
- python - Pyinstaller 找不到熊猫钩子
- sql - Sql (Oracle) insert --column 不允许在这里
- sql - SQL:在同一个表中获取重复项
- java - Java XPath:找不到元素
- intellij-idea - 具有共享内容根的 Kotlin/IntelliJ 模块
- sql - SQL - 创建具有重复值的新列
- android - 安卓工作室 3.1.2