图像特征提取matlab程序

%直接帧间差分，计算阈值并进行二值化处理（效果不好）

clc;

clear;

Im1 = double(imread('lena.TIF')); %读取背景图片

Im2 = double(imread('lena.TIF'); %读取当前图片

[X Y Z] = size(Im2); %当前图片的各维度值

DIma = zeros(X,Y);

for i = 1:X

for j = 1:Y

DIma(i,j) =Im1(i,j) - Im2(i,j); %计算过帧间差分值

end

end

figure,imshow(uint8(DIma)) %显示差分图像

title('DIma')

med = median(DIma); %计算二值化阈值：差值图像中值

mad = abs(mean(DIma) - med); %中值绝对差

T = mean(med + 3*1.4826*mad) %初始阈值

Th =5*T; %调整阈值

BW = DIma <= Th; %根据阈值对图像进行二值化处理

figure,imshow(BW)

%se = strel('disk',2); %膨胀处理

%BW = imopen(BW,se);

%figure,imshow(BW)

%title('BW')

[XX YY] = find(BW==0); %寻找有效像素点的最大边框

handle = rectangle('Position',[min(YY),min(XX) ,max(YY)-min(YY),max(XX)-min(XX)]); set(handle,'EdgeColor',[0 0 0]);

hei = max(XX)-min(XX); %边框高度

mark = min(YY)+1;

while mark < max(YY)-1 %从边框左边开始到右边物质循环，寻找各个人体边缘

left = 0;right = 0;

for j = mark:max(YY)-1

ynum = 0;

for i = min(XX)+1 : max(XX)-1

if BW(i,j) == 0;

ynum = ynum + 1;

end

end

if ynum > 0.1*hei %如果该列上像素点个数大于边框高度的某个阈值left = j %认为找到了左边界

mark = j;

break

end

end

if left~= 0 %如果找到了左边界

for j = left+1 : max(YY)-1 %则从左边界到边框右边界寻找人体右边界ynum = 0;

for i = min(XX)+1 : max(XX)-1

if BW(i,j) == 0;

ynum = ynum + 1;

end

end

if ynum < 0.1*hei %如果该列上像素点个数小于边框高度的某个阈值

right= j %认为找到了左边界

mark = j;

break

end

end

if left~=right %如果左右边界值不相等，即同时找到了左右边界，继续以左右边界为界寻找上下边界

[top,bottom] = Find_edge(BW,min(XX)+1,max(XX),left,right,0.1);

handle = rectangle('Position',[left,top,right-left,bottom-top]);

set(handle,'EdgeColor',[0 0 0]); %显示图像中人体有效边界end

end

mark = mark + 1;

end % while