USAGE: cepstrum [-s size] [-o overlap] [-r rolloff] [-c constant] [-G] infile outfile
USAGE: cepstrum [-h or -help]

-s ... size ........ percent size of original image to use for the
.................... dimension of regions; 0<integer<100; default=25
-o ... overlap ..... percent overlap of regions; 0<=integer<100; default=50
-r ... rolloff ..... cosine-bell window rolloff in percent for the regions;
.................... 0<integer<100; default=40
-c ... constant .... constant value used in log operation; integer>0;
.................... default=10000
-G ................. create grayscale cepstrum

PURPOSE: To create the cepstrum to depict the type, amount and orientation of camera type blurring in an image.

DESCRIPTION: CEPSTRUM creates the cepstrum of an image to detect its blurring type, amount and orientation. The cepstrum is generate by first computing a running average of the magnitude of the Fourier Transform from a complete set of subsections of the image each of which is windowed using typically about a 40% cosine-bell rolloff to black. The windowing is done to help assure that each subsection conforms to the Fourier Transform assumption that the image is periodic. The negative log is applied to this running average, which is then stretched to span the full dynamic range. Finally, this result is treated as the real component and combined with a black image for the imaginary component and inverse transformed back to the spatial domain. The name cepstrum is derived from the word spectrum, but with the first four letters reversed. It is in some sense a spectrum formed backwards. Any alpha channel on the image will be removed automatically before processing.

For lens defocus, the cepstrum should produce a white ring in a generally dark background. The diameter of the white ring is equivalent to twice the diameter of the defocus amount.

For motion blur, the cepstrum should produce an array of dots on a line oriented at the angle of the motion blur. The distance between the two dots closest and on each side of the center of the cepstrum image is equivalent to twice the amount of motion blur.

ARGUMENTS:

-s size ... SIZE of sub-regions used to generate the cepstrum expressed as percent of original image dimensions. Values are integers such that 0<size<100. The default=25

-o overlap ... OVERLAP is the percent overlap of the sub-regions used to generate the cepstrum. Values are integers such that 0<=overlap<100. Zero indicates that the sub-regions do not overlap. The default=50. The number of regions processed will be approximately ((100/overlap)*(100/size)-1)^2. With size=25 and overlap=50, this amounts to approximate 49 regions.

-r rolloff ... ROLLOFF is the cosine-bell shaped windowing fuction rolloff in percent of the image dimensions used to taper the sub-regions from full image intensities to black near its edges. Values are integers such that 0<rolloff<100. The default=40.

-c constant ... CONSTANT is the multiplier value used in the log enhancement. Smaller values attenuate noise more than larger values. Values are integers>0. The default=10000

-G ... Creates GRAYSCALE cepstrum by converting the input image to grayscale before processing.

REQUIREMENTS: IM version 6.5.4-7 or higher, but compiled with HDRI enabled in any quantum level of Q8, Q16 or Q32. Also requires the FFTW delegate library.

See Fourier Transform with ImageMagick, for more details.

CAVEAT: No guarantee that this script will work on all platforms, nor that trapping of inconsistent parameters is complete and foolproof. Use At Your Own Risk.