I have a sliding window on my image. If the average intensity within that window is > 210 then GLCM features are calculated on that window. If the GLCM features meet the conditions in the if statements then a rectangle should be drawn around this sliding window. I have tried this using the following code however the rectangle is not in the correct position on the image. I am unsure if I have put the code which draws the rectangle in the wrong place or if I have passed in the wrong coordinates. Can anyone help me with this please?
N2=8;
info2 = repmat(struct, ceil(size(Z, 1) / N2), ceil(size(Z, 2) / N2));
for row1 = 1:N2:size(Z, 1)%loop through each pixel in the 8x8 window
for col1 = 1:N2:size(Z, 2)
x = (row1 - 1) / N2 + 1;
y = (col1 - 1) / N2 + 1;
imgWindow2 = Z(row1:min(end,row1+N2-1), col1:min(end,col1+N2-1));
average2 = mean(imgWindow2(:));
window2(x,y).average=average2;
if average2>210
offsets0 = [0 1];
glcms = graycomatrix(imgWindow2,'Offset',offsets0);
stats = graycoprops(glcms,'all');
correlation=[stats.Correlation];
contrast=[stats.Contrast];
homogeneity=[stats.Homogeneity];
energy=[stats.Energy];
%if these conditions are met then this window contains an ROI
if (homogeneity > 0.9)
if (contrast<0.2)
if (energy>0.6)
%%show the ROI on the original image%%
Z1 = insertShape(Z, 'rectangle', [x y 32 32]);
figure(2);
imshow(Z1);
end
end
end
end % end if>210
end
end % window