Script to obtain the radial intensity from a 3D image, e.g. a power spectrum.
Usage: just make sure the 3D power spectrum is the front image. The image size in x,y and z should be the same, although this is not checked. Correct the sx as the original pixelsize from the image acquired with the microscope. Choose a step size for the radial plot e.g 0.1 or 0.2. There is no shell width.
The edge variable is for cutting the stack in z for maximum finding, e.g find max from zsize/2 -edge to zsize/2 + edge.