I have the following function which calculates a GLCM and then a given statistic parameter. I would like to pass this function to NLFILTER to do the calculation for a whole image (in small windows, e.g. convolution). I already have then NLFILTER set up to run using the parallel computing toolbox, so I would really like to convert the function I have below:
function [s]=glcm(img,meth)
%GLCM calculates a Gray Level Co-occurence matrix & stats for a given sub
% image.
% Input: Sub Image (subI) and a method (meth)...
% 'Contrast','Correlation','Energy','Homogeneity'
%
subI=uint8(img);
m=graycomatrix(img,'Offset',[0 1],'NumLevels',8,'Symmetric',true);
if meth(1:3)=='con'
s=graycoprops(m,'Contrast');
s=s.Contrast;
elseif meth(1:3)=='cor'
s=graycoprops(m,'Correlation');
s=s.Correlation;
elseif meth(1:3)=='ene'
s=graycoprops(m,'Energy');
s=s.Energy;
elseif meth(1:3)=='hom'
s=graycoprops(m,'Homogeneity');
s=s.Homogeneity;
else
error('No method selected.')
end
I am really stuck on how to convert this to a function handle suitable for use with NLFILTER. Any ideas? Thanks.
When you create an anonymous function, you can pass additional, static arguments in the function definition: