圖像分析初步實驗報告

圖像分析初步上機實驗報告姓名：班級：專業(yè)：指導老師：2023.11靜態(tài)圖象的分析題目要求靜態(tài)圖象的分析：根據1組靜態(tài)圖象，利用一組2幅或3幅圖象分析圖象中的目標的根本形狀，空間關系。題目分析如下圖的運動器材，結構較為復雜，各個局部的位置關系各不相同，三幅圖像突出表現的部位不同，通過處理三幅圖像，分別得到不同部位的位置關系。通過對原始圖像進行膨脹腐蝕等運算，得到目標區(qū)域的連通域，通過目標區(qū)域的中軸線關系，得到目標區(qū)域的位置關系。算法設計程序代碼對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);%%%%%%%%%%%%%%%%%%%保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%%%%%%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);%%%%%%%%%%%%%%%%%%%%進一步優(yōu)化保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%%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);%%%%%%%%%%%%%%%%%%%%%%%%%%曲線擬合%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%holdonbw5=imclearborder(bw8);bw5=imopen(bw5,strel('disk',2));bw5=bwareaopen(bw5,150);cc=bwconncomp(bw5);t=1:size(bw5,2);k=[00];fori=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:\material\learning\Image\tuxiang\HOMEWORK\p06-2.jpg');I1=rgb2gray(I);%彩色轉灰度圖figure(1);imshow(I);SE=strel('square',11);X=imclose(I1,SE);Y=imopen(X,SE);J=im2bw(Y,0.5);%灰度轉二值圖P=imcomplement(J);%圖像反相SE1=strel('line',25,100);bw2=imopen(P,SE1);figure(2);imshow(bw2);%%%%%%%%%%%%%%%%%%保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%%%%%L=bwlabeln(bw2,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);bw=ismember(L,find([S.Area]==P2));bw1=ismember(L,find([S.Area]==P3));SE1=strel('square',50);bw2=imopen(bw,SE1);%%%%%%%%%%%%%%%%%進一步優(yōu)化保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%%%%L1=bwlabeln(bw2,8);%連通域分析S1=regionprops(L1,'Area');%獲得各連通域面積aera2=cat(1,S1.Area);Q1=max(aera2);A1=find([S1.Area]==Q1);aera2(A1)=0;Q2=max(aera2);B1=find([S1.Area]==Q2);aera2(B1)=0;Q3=max(aera2);bw3=ismember(L1,find([S1.Area]==Q2));bw4=ismember(L1,find([S1.Area]==Q3));bw5=bw3+bw4+bw1;figure(3);imshow(bw5);holdonbw5=imclearborder(bw5);bw5=imopen(bw5,strel('disk',2));bw5=bwareaopen(bw5,150);cc=bwconncomp(bw5);t=1:size(bw5,2);k=[000];fori=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;%靠背的傾斜角angle2=atan(k(2))*180/pi;%座椅的傾斜角angle3=atan(k(3))*180/pi;%把手的傾斜角angle4=180-abs(angle1)-abs(angle2);%兩個連通域的夾角angle5=180-abs(angle2)-abs(angle3);%兩個連通域的夾角對p06-3.JPG圖像的處理程序%%%%%%%%%%%%%%%%%%%%%%%%%%%圖片3處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I=imread('E:\material\learning\Image\tuxiang\HOMEWORK\p06-3.jpg');I1=rgb2gray(I);%彩色轉灰度圖figure(1);imshow(I);SE=strel('square',15);X=imclose(I1,SE);Y=imopen(X,SE);J=im2bw(Y,0.5);%灰度轉二值圖P=imcomplement(J);%圖像反相figure(2);imshow(P);%%%%%%%%%%%%%%%%%%%%%%保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%%%%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);bw1=ismember(L,find([S.Area]==P3));SE1=strel('line',50,100);bw2=imopen(bw1,SE1);%%%%%%%%%%%%%%%%%%%進一步優(yōu)化保存目標區(qū)域連通域%%%%%%%%%%%%%%%%%%%%%%%L1=bwlabeln(bw2,4);%連通域分析S1=regionprops(L1,'Area');%獲得各連通域面積aera2=cat(1,S1.Area);Q1=max(aera2);A1=find([S1.Area]==Q1);aera2(A1)=0;Q2=max(aera2);bw3=ismember(L1,find([S1.Area]==Q1));bw4=ismember(L1,find([S1.Area]==Q2));bw5=bw3+bw4;figure(3);imshow(bw5);%%%%%%%%%%%%%%%%%%%%%%%%%%曲線擬合%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%holdonbw5=imclearborder(bw5);bw5=imopen(bw5,strel('disk',2));bw5=bwareaopen(bw5,150);cc=bwconncomp(bw5);t=1:size(bw5,2);k=[00];fori=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=180-abs(angle1)-abs(angle2);%兩個連通域的夾角運行效果1.對p06-1.JPG圖像的處理結果如上圖所示，座椅與底座的夾角為19.1705度。對p06-2.JPG圖像的處理結果如上圖所示，靠背與座椅的夾角為121.6206度，座椅與把手的夾角為129.7893度。對p06-3.JPG圖像的處理結果如上圖所示，把手區(qū)域兩連通域的夾角為55.2599度?？偨Y處理靜態(tài)圖像主要應用到的算法是，膨脹腐蝕算法得到包含目標區(qū)域的連通域，通過對各個連通域的面積進行比擬，選取自己想得到的目標區(qū)域，各部位的位置關系是通過曲線擬合得到連通域的對稱軸，求對稱軸的夾角，即得到各部件的位置關系。因目標區(qū)域顏色都偏深，圖像二值化以后，目標區(qū)域都是黑色，所以通過imcomplement函數進行圖像反相；通過regionprops函數獲得各個連通域面積；通過polyfit函數擬合對稱軸曲線。動態(tài)圖象的分析一、題目要求提取足球場上黃色球衣和蘭色球衣兩個同學的空間運動軌跡。二、題目分析首先對視頻進行圖片分割，然后通過膨脹腐蝕等算法進行圖片處理，找到目標區(qū)域的質心，將質心連接起來就是該目標任務的運動軌跡。三、算法設計程序代碼視頻別離程序%%讀取視頻video_file='E:\material\learning\Image\tuxiang\HOMEWORK\m06.avi';video=VideoReader(video_file);frame_number=floor(video.Duration*video.FrameRate);%%別離圖片fori=1:frame_numberimage_name=strcat('m06-',num2str(i));image_name=strcat(image_name,'.jpg');I=read(video,i);%讀出圖片imwrite(I,image_name,'jpg');%寫圖片I=[];end圖片處理程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%函數名稱:TrackBlue%功能:基于顏色的目標追蹤%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functionTrackBlueforn=1:26imageName=strcat('m06-',num2str(n*5+330));imageName=strcat(imageName,'.jpg');[Xpos,Ypos]=Tracking(imageName);XposA(n)=Xpos;YposA(n)=Ypos;endbgim=imread('E:\material\learning\Image\tuxiang\HOMEWORK\m06-1.jpg');bgimG=bgim(280:479,150:550);figure(2);imshow(bgimG);%背景圖holdonplot(XposA,YposA,'y-o','LineWidth',2,'MarkerFacecolor','y');%運動路徑標記axison;gridon;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%函數名稱:Tracking%功能:目標追蹤%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function[Xpos,Ypos]=Tracking(imageName)I=imread(imageName);I1=I(280:479,150:550);R=I(:,:,1);%紅色分量G=I(:,:,2);%綠色分量B=I(:,:,3);%藍色分量J1=B*4-G*2-R*2;%突出藍色分量J=J1(280:479,150:550);%圖像截取SE1=strel('disk',2);SE2=strel('square',15);X=imerode(J,SE1);%圖像腐蝕Y=imdilate(X,SE2);%圖像膨脹K=im2bw(Y,0);