原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
%直接帧间差分，计算阈值并进行二值化处理（效果不好）
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,botto
文档评论（0）