目地是进行视频中人体动作的识别,这个代码好像是计算人体边界轮廓特征参数,与模版进行比对。 我知道里面有面积 周长 宽 高 ,可是那个旋转计算交点个数那个参数是什么啊 看了很多论文都没有看到相关的 。
代码如下:
im=im2double(imread('bk.jpg'));
figure;imshow(im,[]);title('原图1')
if length(size(im))==3
im=rgb2gray(im);
end
figure;imshow(im,[]);title('原图1灰度图')
im2=im2double(imread('box10.jpg'));
figure;imshow(im2,[]);title('原图2')
if length(size(im2))==3
im2=rgb2gray(im2);
end
figure;imshow(im2,[]);title('原图2灰度图')
dif=abs(im2-im);
figure;imshow(dif,[]);title('两帧图之间的差异图')
[idx,idy]=find(dif>0.1);
imk=zeros(size(im));
imk(min(idx):max(idx),min(idy):max(idy))=im2(min(idx):max(idx),min(idy):max(idy));
figure;imshow(imk,[]);title('通过帧间差阈值提取的目标图')
tem=im2(min(idx):max(idx),min(idy):max(idy)); %截图目标图
thr=graythresh(tem);
bw=im2bw(tem,thr);
bw=1-bw; %对目标图进行二值化
figure;imshow(bw);title('二值图')
Area=sum(bw(:)); %面积
[m,n]=size(bw);
mk=zeros(m,n);
for ss=2:m-1
for tt=2:n-1
if bw(ss,tt)+bw(ss,tt+1)==1 %两个相邻点之和等于1,则该点为边界点,保存下来
mk(ss,tt)=1;
end
if bw(ss,tt)+bw(ss+1,tt)==1 %两个相邻点之和等于1,则该点为边界点,保存下来
mk(ss,tt)=1;
end
end %提取二值图像人物轮廓
end
figure;imshow(mk,[])
ZC=sum(mk(:)); %周长
[idx,idy]=find(mk==1);
hei=max(idx)-min(idx)+1; %高
wid=max(idy)-min(idy)+1; %宽
flag=1;
for ss=1:10:180 %围绕中心间隔10度旋转一次,对旋转之后的图像,计算中心的一条线上的交点个数
tem=imrotate(mk,ss,'nearest','loose'); %最邻近线性插值
[m,n]=size(tem);
Par(flag)=sum(tem(round(m/2),:));
flag=flag+1;
end
P1=[Area,ZC,hei,wid];
Parameter=[Par,P1] %将参数拼接起来
代码如下:
im=im2double(imread('bk.jpg'));
figure;imshow(im,[]);title('原图1')
if length(size(im))==3
im=rgb2gray(im);
end
figure;imshow(im,[]);title('原图1灰度图')
im2=im2double(imread('box10.jpg'));
figure;imshow(im2,[]);title('原图2')
if length(size(im2))==3
im2=rgb2gray(im2);
end
figure;imshow(im2,[]);title('原图2灰度图')
dif=abs(im2-im);
figure;imshow(dif,[]);title('两帧图之间的差异图')
[idx,idy]=find(dif>0.1);
imk=zeros(size(im));
imk(min(idx):max(idx),min(idy):max(idy))=im2(min(idx):max(idx),min(idy):max(idy));
figure;imshow(imk,[]);title('通过帧间差阈值提取的目标图')
tem=im2(min(idx):max(idx),min(idy):max(idy)); %截图目标图
thr=graythresh(tem);
bw=im2bw(tem,thr);
bw=1-bw; %对目标图进行二值化
figure;imshow(bw);title('二值图')
Area=sum(bw(:)); %面积
[m,n]=size(bw);
mk=zeros(m,n);
for ss=2:m-1
for tt=2:n-1
if bw(ss,tt)+bw(ss,tt+1)==1 %两个相邻点之和等于1,则该点为边界点,保存下来
mk(ss,tt)=1;
end
if bw(ss,tt)+bw(ss+1,tt)==1 %两个相邻点之和等于1,则该点为边界点,保存下来
mk(ss,tt)=1;
end
end %提取二值图像人物轮廓
end
figure;imshow(mk,[])
ZC=sum(mk(:)); %周长
[idx,idy]=find(mk==1);
hei=max(idx)-min(idx)+1; %高
wid=max(idy)-min(idy)+1; %宽
flag=1;
for ss=1:10:180 %围绕中心间隔10度旋转一次,对旋转之后的图像,计算中心的一条线上的交点个数
tem=imrotate(mk,ss,'nearest','loose'); %最邻近线性插值
[m,n]=size(tem);
Par(flag)=sum(tem(round(m/2),:));
flag=flag+1;
end
P1=[Area,ZC,hei,wid];
Parameter=[Par,P1] %将参数拼接起来
