圖像特征提取matlab程序

圖像特征提取matlab程序圖像特征提取matlab程序圖像特征提取matlab程序xxx公司圖像特征提取matlab程序文件編號(hào)：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計(jì)，管理制度%直接幀間差分，計(jì)算閾值并進(jìn)行二值化處理（效果不好）clc;clear;Im1=double(imread(''));%讀取背景圖片Im2=double(imread('');%讀取當(dāng)前圖片[XYZ]=size(Im2);%當(dāng)前圖片的各維度值DIma=zeros(X,Y);fori=1:Xforj=1:YDIma(i,j)=Im1(i,j)-Im2(i,j); %計(jì)算過(guò)幀間差分值endendfigure,imshow(uint8(DIma))%顯示差分圖像title('DIma')med=median(DIma);%計(jì)算二值化閾值：差值圖像中值mad=abs(mean(DIma)-med);%中值絕對(duì)差T=mean(med+3**mad)%初始閾值Th=5*T;%調(diào)整閾值BW=DIma<=Th;%根據(jù)閾值對(duì)圖像進(jìn)行二值化處理figure,imshow(BW)%se=strel('disk',2);%膨脹處理%BW=imopen(BW,se);%figure,imshow(BW)%title('BW')[XXYY]=find(BW==0);%尋找有效像素點(diǎn)的最大邊框handle=rectangle('Position',[min(YY),min(XX),max(YY)-min(YY),max(XX)-min(XX)]);set(handle,'EdgeColor',[000]);hei=max(XX)-min(XX);%邊框高度mark=min(YY)+1;whilemark<max(YY)-1%從邊框左邊開(kāi)始到右邊物質(zhì)循環(huán)，尋找各個(gè)人體邊緣left=0;right=0;forj=mark:max(YY)-1ynum=0;fori=min(XX)+1:max(XX)-1ifBW(i,j)==0;ynum=ynum+1;endendifynum>*hei%如果該列上像素點(diǎn)個(gè)數(shù)大于邊框高度的某個(gè)閾值left=j%認(rèn)為找到了左邊界mark=j;breakendendifleft~=0%如果找到了左邊界forj=left+1:max(YY)-1%則從左邊界到邊框右邊界尋找人體右邊界ynum=0;fori=min(XX)+1:max(XX)-1ifBW(i,j)==0;ynum=ynum+1;endendifynum<*hei%如果該列上像素點(diǎn)個(gè)數(shù)小于邊框高度的某個(gè)閾值right=j%認(rèn)為找到了左邊界mark=j;breakendendifleft~=right%如果左右邊界值不相等，即同時(shí)找到了左右邊界，繼續(xù)以左右邊界為界尋找上下邊界[top,bottom]=Find_edge(BW,min(XX)+1,max(XX),left,right,;handle=rectangle('Position',[left,top,right-left,bottom-top]);