title: Power spectral density of 2D field - Python
body: <p>I would like to use Welch's method for calculating the power spectral density of a 2D field. There is an implementation available in Scipy, but according to the <a href="https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.welch.html" rel="noreferrer">docs</a> it will only work for 1D timeseries.</p>
  
  <p>Is anyone aware of an implementation that would work for 2D fields?</p>
  
  <p>So far I have been using the 2D Fourier transform to do the calculation, but I have heard that this method is sensitive to noise and that Welch's method is more robust.</p>