feature_extract.m

libsvm支持向量机

function energy=feature_extract(img)

enhanceresult = histeq(img);
X=enhanceresult;

nlevels = [5] ;        % Decomposition level
pfilter = 'pkva' ;              % Pyramidal filter
dfilter = 'pkva' ;              % Directional filter

% Contourlet decomposition
coeffs = pdfbdec( double(X), pfilter, dfilter, nlevels );

% for i=1:32
%     figure,imshow(coeffs{2}{i},[])
% end
energy=zeros(1,32);
for i=1:32
    energy(i)=sum(sum(coeffs{2}{i}.*coeffs{2}{i}));
end

energysum=sum(energy);

energy=energy/energysum;