图像分析初步实验报告2015年.docxVIP

图像分析初步上机实验报告姓名：班 级：专业：指导老师：2015.11静态图象的分析题目要求静态图象的分析：根据1组静态图象，利用一组2幅或3幅图象分析图象中的目标的基本形状，空间关系。题目分析如图所示的运动器材，结构较为复杂，各个部分的位置关系各不相同，三幅图像突出表现的部位不同，通过处理三幅图像，分别得到不同部位的位置关系。通过对原始图像进行膨胀腐蚀等运算，得到目标区域的连通域，通过目标区域的中轴线关系，得到目标区域的位置关系。算法设计程序代码对p06-1.JPG图像的处理程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%图片1处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I = imread(E:\material\learning\Image\tuxiang\HOMEWORK\p06-1.jpg);I1 = rgb2gray(I); % 彩色转灰度图figure(1);imshow(I1);SE=strel(square,15);X = imclose(I1,SE);Y = imopen(X,SE);J = im2bw(Y,0.5); %灰度转二值图P = imcomplement(J); %图像反相figure(2);imshow(P);%%%%%%%%%%%%%%%%%%% 保留目标区域连通域%%%%%%%%%%%%%%%%%%%%%%%%%%%%L=bwlabeln(P,8); % 连通域分析S=regionprops(L,Area); % 获得各连通域面积aera3=cat(1,S.Area);P1=max(aera3);A=find([S.Area]==P1);aera3(A)=0;P2=max(aera3);B = find([S.Area]==P2);aera3(B) = 0;P3 = max(aera3);C = find([S.Area]==P3);bw1=ismember(L,find([S.Area]==P3)); % 保留面积次大的连通域aera3(C) = 0;P4 = max(aera3);D = find([S.Area]==P4);bw2=ismember(L,find([S.Area]==P4)); % 保留面积次大的连通域bw4 = bw1+bw2;SE1 = strel(line,150,0);bw3 = imopen(bw4,SE1);%%%%%%%%%%%%%%%%%%%%进一步优化保留目标区域连通域%%%%%%%%%%%%%%%%%%%%%%%%L1=bwlabeln(bw3,4); % 连通域分析S1=regionprops(L1,Area); % 获得各连通域面积aera2=cat(1,S1.Area);Q1=max(aera2);A1 = find([S1.Area]==Q1);aera2(A1)=0;Q2=max(aera2);bw6 = ismember(L1,find([S1.Area]==Q1)); bw7 = ismember(L1,find([S1.Area]==Q2)); bw8 = bw6 + bw7;figure(3);imshow(bw8);%%%%%%%%%%%%%%%%%%%%%%%%%%曲线拟合%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%hold onbw5=imclearborder(bw8);bw5=imopen(bw5,strel(disk,2));bw5=bwareaopen(bw5,150);cc=bwconncomp(bw5);t=1:size(bw5,2);k = [0 0];for i=1:cc.NumObjects [y,x]=ind2sub(size(bw5),cc.PixelIdxList{i}); p=polyfit(x,y,1); plot(t,p(1)*t+p(2),r) k(i)= p(1);endangle1 = atan(k(1))*180/pi; %连通域1的倾斜角angle2 = atan(k(2))*180/pi; %连通域2的倾斜角angle = abs(angle1) + abs(angle2); %两个连通域的夹角对p06-2.JPG图像的处理程序%%%%%%%%%%%%%%%%%%%%%%%%图片1处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I = imread(E