PW L - Power spectrum - Log10

PURPOSE

Generates full, unscrambled 0.5 * logarithm of the power spectrum from complex Fourier transform for 2-D or 3-D pictures. The input 2D image can be real if it fits into the memory. Operations supports multiple radix FFT, thus the input image does not have to have power-of-two dimensions The resulting data are real and can be displayed like a normal image.

SEE ALSO

PW	[Amplitudes of Fourier transform]
FT	[Fourier Transform ||]

USAGE

.OPERATION: PW L

.INPUT FILE: FOU001
[Enter the file name of the Fourier transform.]

.OUTPUT FILE: POW001
[Enter the name of the file where the 0.5*LOG10(power spectrum) is to be kept.]

NOTES

1. The center of 2-D power spectrum is placed in (NSAM/2+1, NROW/2+1);

2. The center of 3-D power spectrum is placed in (NSAM/2+1, NROW/2+1, NSLICE/2+1).

3. The term |F(0)| is replaced by term |F(1)| to facilitate visualization.

SUBROUTINES

PW2SDR, PW2SR, PW3SDR, PW3SR

CALLER

FOUR1 , FOUR1B