0

我的图像上有一个滑动窗口。如果该窗口内的平均强度 > 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
4

1 回答 1

0

尝试使用 [col1 row1 32 32] 而不是 [xy 32 32]。我不明白你为什么要以你的方式计算 x 和 y 。您最终会切换 x 和 y，它们位于“单元格”中，而不是像素中。

于 2014-04-24T14:12:24.277 回答