hall3
Contents
read image
small_area = 500;
rgb = imread('img4.jpg');
sz = size(rgb)
if length(sz)>2
g = im2double(rgb2gray(rgb));
else
g = im2double(rgb);
end
imshow(rgb);
sz =
960 1280 3
Warning: Image is too big to fit on screen; displaying
at 50%
black = g<0.3;
imshow(black);
Warning: Image is too big to fit on screen; displaying
at 50%
white = (g>0.5);
imshow(white);
Warning: Image is too big to fit on screen; displaying
at 50%
se1 = strel('square',7);
b1 = imdilate(black,se1);
w1 = imdilate(white,se1);
bw1 = w1&b1;
bw2 = bwareaopen(bw1,small_area);
imshow(bw2);
Warning: Image is too big to fit on screen; displaying
at 50%
CC = bwconncomp(bw2);
disp(CC);
L = labelmatrix(CC);
imshow(label2rgb(L));
Connectivity: 8
ImageSize: [960 1280]
NumObjects: 2
PixelIdxList: {[6346x1 double] [60174x1 double]}
Warning: Image is too big to fit on screen; displaying
at 50%
s = regionprops(CC,'Image','FilledImage','Area','EulerNumber');
areas = [s.Area];
[maxarea idx] = max(areas);
imshow(s(idx).Image);
bw3 = s(idx).FilledImage&~s(idx).Image;
bw4 = bwareaopen(bw3,small_area);
CC2 = bwconncomp(bw4);
L2 = labelmatrix(CC2);
imshow(label2rgb(L2));
s2 = regionprops(CC2,'Area','Centroid');
area = cat(1,s2.Area);
fprintf('number of squares %d\n',length(area));
fprintf('max area %d min area %d\n',max(area),min(area));
fprintf('Euler %d\n',-s(idx).EulerNumber);
number of squares 33
max area 12177 min area 530
Euler 32
pts = reshape([s2.Centroid],2,length(s2));
plot(pts(1,:),pts(2,:),'ko');
hold off
axis ij
axis equal
grid
rgbm = label2RGB(L,'jet','k');
mask = min(rgbm,[],3)==0;
gm = g.*mask;
rgbn = cat(3,gm,gm,gm)+im2double(rgbm);
imshow(rgbn);
Warning: Image is too big to fit on screen; displaying
at 50%
