У меня есть скользящее окно на моем изображении. Если средняя интенсивность в этом окне > 210, то в этом окне рассчитываются признаки GLCM. Если функции GLCM соответствуют условиям в операторах if, то вокруг этого скользящего окна должен быть нарисован прямоугольник. Я пробовал это, используя следующий код, однако прямоугольник не находится в правильном положении на изображении. Я не уверен, поместил ли я код, рисующий прямоугольник, в неправильное место или передал неправильные координаты. Может ли кто-нибудь помочь мне с этим, пожалуйста?
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