speckle_ps_diff function

Power spectrum calculation