1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
|
% Lecture de l'image
X = imread('si-13refait.jpg','jpg');
% Réduction de l'image à la zone d'intérêt
BW = X(:,:,1)~=255 & X(:,:,2)~=255 & X(:,:,3)~=255;
idx(1) = find(any(BW,2),1,'first');
idx(2) = find(any(BW,2),1,'last');
X = X(idx(1):idx(2),:,:);
idx(1) = find(any(BW,1),1,'first');
idx(2) = find(any(BW,1),1,'last');
X = X(:,idx(1):idx(2),:);
% Binarisation de l'image
BW = X(:,:,1)==255 & X(:,:,2)==255 & X(:,:,3)==255;
% Séparation en région d'intérêt ROI
L = bwlabel(BW,4);
% Surface des ROI
reg = regionprops(L,'area');
% Sélection des ROI suffisamment grande
idx = find([reg.Area]>500);
L(:) = L(:).*ismember(L(:),idx(:));
% Affichage
figure
imagesc(L);
axis image |
Partager