我正在尝试用 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?