speckle_ps function

Power spectrum calculation