圖像分析初步上機實驗報告姓名：班級：專業(yè)：指導(dǎo)老師：2023.11靜態(tài)圖象的分析題目要求靜態(tài)圖象的分析：根據(jù)1組靜態(tài)圖象，利用一組2幅或3幅圖象分析圖象中的目標的根本形狀，空間關(guān)系。題目分析如下圖的運動器材，結(jié)構(gòu)較為復(fù)雜，各個局部的位置關(guān)系各不相同，三幅圖像突出表現(xiàn)的部位不同，通過處理三幅圖像，分別得到不同部位的位置關(guān)系。通過對原始圖像進行膨脹腐蝕等運算，得到目標區(qū)域的連通域，通過目標區(qū)域的中軸線關(guān)系，得到目標區(qū)域的位置關(guān)系。算法設(shè)計程序代碼對p06-1.JPG圖像的處理程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%圖片1處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I=imread('E:\material\learning\Image\tuxiang\HOMEWORK\p06-1.jpg');I1=rgb2gray(I);%彩色轉(zhuǎn)灰度圖figure(1);imshow(I1);SE=strel('square',15);X=imclose(I1,SE);Y=imopen(X,SE);J=im2bw(Y,0.5);%灰度轉(zhuǎn)二值圖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);%彩色轉(zhuǎn)灰度圖figure(1);imshow(I);SE=strel('square',11);X=imclose(I1,SE);Y=imopen(X,SE);J=im2bw(Y,0.5);%灰度轉(zhuǎn)二值圖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);%彩色轉(zhuǎn)灰度圖figure(1);imshow(I);SE=strel('square',15);X=imclose(I1,SE);Y=imopen(X,SE);J=im2bw(Y,0.5);%灰度轉(zhuǎn)二值圖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圖像的處理結(jié)果如上圖所示，座椅與底座的夾角為19.1705度。對p06-2.JPG圖像的處理結(jié)果如上圖所示，靠背與座椅的夾角為121.6206度，座椅與把手的夾角為129.7893度。對p06-3.JPG圖像的處理結(jié)果如上圖所示，把手區(qū)域兩連通域的夾角為55.2599度?？偨Y(jié)處理靜態(tài)圖像主要應(yīng)用到的算法是，膨脹腐蝕算法得到包含目標區(qū)域的連通域，通過對各個連通域的面積進行比擬，選取自己想得到的目標區(qū)域，各部位的位置關(guān)系是通過曲線擬合得到連通域的對稱軸，求對稱軸的夾角，即得到各部件的位置關(guān)系。因目標區(qū)域顏色都偏深，圖像二值化以后，目標區(qū)域都是黑色，所以通過imcomplement函數(shù)進行圖像反相；通過regionprops函數(shù)獲得各個連通域面積；通過polyfit函數(shù)擬合對稱軸曲線。動態(tài)圖象的分析一、題目要求提取足球場上黃色球衣和蘭色球衣兩個同學(xué)的空間運動軌跡。二、題目分析首先對視頻進行圖片分割，然后通過膨脹腐蝕等算法進行圖片處理，找到目標區(qū)域的質(zhì)心，將質(zhì)心連接起來就是該目標任務(wù)的運動軌跡。三、算法設(shè)計程序代碼視頻別離程序%%讀取視頻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圖片處理程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%函數(shù)名稱: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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%函數(shù)名稱: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);