Image Features From Phase Congruency

Peter Kovesi
Department of Computer Science
The University of Western Australia


MATLAB Code for the Algorithms Presented in the Paper :

phasecong.m MATLAB code for calculating Phase Congruency and Phase Symmetry/Asymmetry.

nonmaxsup.m MATLAB code for non-maxima suppression.

hysthresh.m MATLAB code for hysteresis thresholding.


Example MATLAB session:


 >> [im, map] = tiffread('picci.tif');    % Read the image and its colour map
 >> colormap(map);                        % Set colour map
 >> image(im);                            % Display the image
 >> [pc orient] = phasecong(image);       % Calculate phase congruency.
 >> imagesc(pc);                          % Display the phase congruency image
 >> nonmax = nonmaxsup(pc, orient, 1.5);  % Perform non-maximal suppression
 >> features =  hysthresh(nonmax, 0.3, 0.15); % Hysteresis threshold
 >> imagesc(features);                    % Look at the result