BUTTERWORTH_BANDREJECT_CVIP

butterworth_bandreject_cvip() - performs butterworth band reject filtering operation.

SYNTAX
DESCRIPTION
REFERENCE
EXAMPLE
CREDITS

SYNTAX

out_s = butterworth_bandreject_cvip( spectrum, block_size, transform_type,order, fc_low, fc_high )

Input parameters include :

spectrum - An image of size [m n d] which is the spectrum of an image. It can be created by any of the 6 possible spectrums availabe in cviptools matlab toolbox: DCT, FFT, Haar, Walsh,Hadamard, Wavelet.
block_size - The same parameter block_size used in the transforms.
transform_type - The user should specify which spectrum is being given to this function. FFT has the origin at the center while the other transforms have their origin located at upper left corner of the image.transform_type is a string. If given spectrum is FFT, then set transform_type = 'fft' or 'center'.Else, the origin would be located at the upper left corner.
order - The butterworth filters parameter.An integer greater than or equal to 1.As order increases the filter frequency response approches the ideal filter.
fc_low - The low cutoff frequency for band pass filter.
fc_high - The high cutoff frequency for band pass filter.

Output parameter include :

out_s - output image.

DESCRIPTION

This function performs 2D bandreject filtering on the spectrum of an Image. It applies butterowrth filter of order specified by the user to the input spectrum and passes the frequencies that lie outside the specified frequency band.

REFERENCE

1.Scott E Umbaugh DIGITAL IMAGE PROCESSING AND ANALYSIS: Applications with MATLAB and CVIPtools 3rd Edition.

EXAMPLE

% Create an spectrum by available transforms

 spectrum = ones(128,256,3);

% Block size

 block_size1  = [];
 block_size2 = [64 128];

% Transform type

 transform_type1 =' fft';
 transform_type2 ='non-fft';
 transform_type3 ='center';

% Order

 order = 6;

% Low frequency

 fc_low1 = 32;
 fc_low2 = 23;
 fc_low3 = 16;

% High frequency

 fc_high1 = 45;
 fc_high2 = 55;
 fc_high3 = 32;

% Call function

 out_s1 = butterworth_bandreject_cvip( spectrum,block_size1,transform_type1,order,fc_low1,fc_high1);
 out_s2 = butterworth_bandreject_cvip( spectrum,block_size2, transform_type2,order, fc_low2,fc_high2);
 out_s3 = butterworth_bandreject_cvip( spectrum,block_size2 , transform_type3,order,fc_low3,fc_high2);

% Display output

 figure; imshow(out_s1,[]);title(' Output with FFT Transform');
 figure; imshow(out_s2,[]);title(' Output with Non FFT Transform')
 figure; imshow(out_s3,[]);title('with center');

BUTTERWORTH_BANDREJECT_CVIP

Contents

SYNTAX

DESCRIPTION

REFERENCE

EXAMPLE

CREDITS