for k = 1:nglcm
currentGLCM = glcm(:,:,k);
glcmMean(k) = mean2(currentGLCM);
uX(k) = sum(I.*currentGLCM(sub));
uY(k) = sum(J.*currentGLCM(sub));
sX(k) = sum((I-uX(k)).^2.*currentGLCM(sub));
sY(k) = sum((J-uY(k)).^2.*currentGLCM(sub));
out.contrast(k) = sum(abs(I-J).^2.*currentGLCM(sub));
out.dissimilarity(k) = sum(abs(I - J).*currentGLCM(sub));
out.energy(k) = sum(currentGLCM(sub).^2);
out.entropy(k) = -nansum(currentGLCM(sub).*log(currentGLCM(sub)));
out.inverseDifference(k) = sum(currentGLCM(sub)./( 1 + abs(I-J) ));
out.homogeneity(k) = sum(currentGLCM(sub)./( 1 + (I - J).^2));