凹坑划痕检测17042100.rar_凹坑_凹坑 检测_凹坑划痕_凹坑检测_凹坑；划痕

%%%%%%%%%%%%%%%%%%%%改进的Canny小面积去除%%%%%%%%%%%%%% clc;clear all; I=imread('17041300-04.bmp'); I=rgb2gray(I); figure,imshow(I); B5=edge(I,'canny',0.15); figure,imshow(B5); %%%%%%%%%%%%%%%%%%%%%%%%%%%形态学处理%%%%%%%%%%%%% se=strel('disk',3); %先膨胀后腐蚀称为闭合； 先腐蚀后膨胀称为开启 %se = strel('square',5)% 创建6*6的正方形 %se=strel('periodicline',2,[-1 -1])% 创建平坦线型结构元素 %se= strel('diamond',2)% 创建菱形2是结构元素原点到菱形最远的距离 K1=imclose(B5,se); figure,imshow(K1); figure,imshow(K1); %%%%%%%%%%%%%%%%标记所有连通域%%%%%%%%%%%%%% img_reg1 = regionprops(K1, 'area', 'boundingbox','Eccentricity'); rects1= cat(1, img_reg1.BoundingBox); for i = 1:size(rects1, 1) rectangle('position', rects1(i, :), 'EdgeColor', 'r'); end %%%%%%%%%%%%%%%%去掉极大和极小的连通域%%%%%%%%%%%%%% Ibw = bwareaopen(im2uint8(K1),250);%将小于XX像素的单元去掉 L= bwlabel(Ibw ); img_reg2 = regionprops(L, 'area'); areas2 = [img_reg2.Area]; idx = find(areas2<3000);%%%% BW1 = ismember(L,idx); figure,imshow(BW1); %%%%%%%%%%%%%%%%%%%%%填充图像的空洞%%%%%%%%%%%%%%%%%% I3 = imfill(BW1,'holes');%%%%%%%% %%%%%%%%%%%%%%提取缺陷和划痕%%%%%%%%%%%%%%% L1= bwlabel(I3); img_reg3 = regionprops(L1, 'area'); areas3 = [img_reg3.Area]; idx = find(areas3<5000);%%%% BW2 = ismember(L1,idx); img_reg6 = regionprops(BW2 , 'Centroid'); Centroids6 = [img_reg6.Centroid]; Centroids61=reshape(Centroids6,15,2); figure,imshow(BW2); %%%%%%%%%%%%%% 求解形状特征参数 %%%%%%%%%%%%%%% I40 = edge(BW2); se=strel('disk',3); %先膨胀后腐蚀称为闭合； 先腐蚀后膨胀称为开启 I41=imclose(I40,se); figure,imshow(I41); L = bwlabel(BW2); L1 = bwlabel(I41); L = bwlabel(BW2); Ar = zeros(1,max(L(:))); perimeter = zeros(1,max(L(:))); metric = zeros(1,max(L(:))); X = zeros(1,max(L(:))); Pwl = zeros(1,max(L(:))); Pr = zeros(1,max(L(:))); idx1=[]; idx2=[]; for i=1:max(L(:)) Ar(i)=sum(BW2(L==i));%面积 perimeter(i)=sum(I41(L1==i));%周长 metric(i) = 4*pi*Ar(i)/perimeter(i)^2;%圆形度 X(i)=Ar(i)/perimeter(i);%线度 if metric(i)<=0.1 idx1=[idx1;i]; else idx2=[idx2;i]; end [y,x]=find(L==i); x0=min(x(:)); x1=max(x(:)); y0=min(y(:)); y1=max(y(:)); hold on if x1-x0>=y1-y0 Pwl(i)=(x1-x0)/(y1-y0); else Pwl(i)=(y1-y0)/(x1-x0);%细长度 end Pr(i)=Ar(i)/((y1-y0)*(x1-x0));%占空比 end BW3 = ismember(L,idx1); BW4 = ismember(L,idx2); figure,imshow(BW3); figure,imshow(BW4); shapefeatures=[Ar;perimeter;metric;Pwl;Pr;X]; %%%%%%%%%%%%%%%%标记缺陷位置%%%%%%%%%%%%%%%%%%% figure,imshow(I); img_reg4 = regionprops(BW3, 'area', 'boundingbox','Centroid'); Centroids4 = [img_reg4.Centroid]; rects4 = cat(1, img_reg4.BoundingBox); for i = 1:size(rects4, 1) rectangle('position', rects4(i, :), 'EdgeColor', 'r'); end figure,imshow(I); img_reg5 = regionprops(BW4, 'area', 'boundingbox','Centroid'); Centroids5 = [img_reg5.Centroid]; rects5 = cat(1, img_reg5.BoundingBox); % show all the largest connected region for i = 1:size(rects5, 1) rectangle('position', rects5(i, :), 'EdgeColor', 'r'); end