function LocDensEst=spattempLBP(refimg,curimg,backimg,radius,mbSize,density,u1)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%refimg,curimg,backimg:Three consecutive frames
%Rradius:sample radius
%neighbors:The number of sampling points=8
%u1:in order to get the size of LocDensEst
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
neighbors=8;
%the size of the input image.
% refimg=double(refimg);
% curimg=double(curimg);
% backimg=double(backimg);
[ysize xsize] = size(refimg);
%sample
a = 2*pi/neighbors;
spoints=zeros(neighbors,2);
for i = 1:neighbors
spoints(i,1) =round(-radius*sin((i-1)*a));%dx
spoints(i,2) =round(radius*cos((i-1)*a));%dy
end
%compute the size of neighbor block
miny=min(spoints(:,1));%The first column minima
maxy=max(spoints(:,1));%The first column maximum
minx=min(spoints(:,2));%The second column minima
maxx=max(spoints(:,2));%The second column maximum
% Block size, each LBP code is computed within a block of size bsizey*bsizex
bsizey=maxy-miny;%get the size of bolck -1,even
bsizex=maxx-minx;
%Edge supplement 0
Refimg=zeros(ysize+bsizey,xsize+bsizex);
Curimg=zeros(ysize+bsizey,xsize+bsizex);
Backimg=zeros(ysize+bsizey,xsize+bsizex);
for i=1:ysize+bsizey%%%%
for j=1:xsize+bsizex%%%
if (i<=bsizey/2 || j<=bsizex/2 || i>ysize+bsizey/2 || j>xsize+bsizex/2 )
Refimg(i,j)=0;
Curimg(i,j)=0;
Backimg(i,j)=0;
else
Refimg(i,j)=refimg(i-bsizey/2,j-bsizex/2);
Curimg(i,j)=curimg(i-bsizey/2,j-bsizex/2);
Backimg(i,j)=backimg(i-bsizey/2,j-bsizex/2);
end
end
end
%对前后三帧图像邻域中像素值与中间像素值进行比较得到0、1码,
code1=zeros(1,26);
code2=zeros(1,26);
[r,w]=size(u1);%u1:in order to get the size of LocDensEst
LocDensEst=zeros(r,w);
a=1;b=1;
for i=bsizey/2+1:mbSize*density:ysize+bsizey/2
for j=bsizex/2+1:mbSize*density:xsize+bsizex/2
%spiral order
code1(1)=cmp(Refimg(i,j),Curimg(i,j));%spiral order
code1(2)=cmp(Refimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));
code1(3)=cmp(Refimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code1(4)=cmp(Refimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code1(5)=cmp(Refimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code1(6)=cmp(Refimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code1(7)=cmp(Refimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code1(8)=cmp(Refimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code1(9)=cmp(Refimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code1(10)=cmp(Curimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));%third frame
code1(11)=cmp(Curimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code1(12)=cmp(Curimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code1(13)=cmp(Curimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code1(14)=cmp(Curimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code1(15)=cmp(Curimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code1(16)=cmp(Curimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code1(17)=cmp(Curimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code1(18)=cmp(Backimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));%third frame
code1(19)=cmp(Backimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code1(20)=cmp(Backimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code1(21)=cmp(Backimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code1(22)=cmp(Backimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code1(23)=cmp(Backimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code1(24)=cmp(Backimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code1(25)=cmp(Backimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code1(26)=cmp(Backimg(i,j),Curimg(i,j));
%zigzag order
code2(1)=cmp(Refimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));
code2(2)=cmp(Curimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));
code2(3)=cmp(Backimg(i+spoints(3,1),j+spoints(3,2)),Curimg(i,j));
code2(4)=cmp(Backimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code2(5)=cmp(Curimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code2(6)=cmp(Refimg(i+spoints(2,1),j+spoints(2,2)),Curimg(i,j));
code2(7)=cmp(Refimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code2(8)=cmp(Curimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code2(9)=cmp(Backimg(i+spoints(1,1),j+spoints(1,2)),Curimg(i,j));
code2(10)=cmp(Backimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code2(11)=cmp(Curimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code2(12)=cmp(Refimg(i+spoints(8,1),j+spoints(8,2)),Curimg(i,j));
code2(13)=cmp(Refimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code2(14)=cmp(Curimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code2(15)=cmp(Backimg(i+spoints(7,1),j+spoints(7,2)),Curimg(i,j));
code2(16)=cmp(Backimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code2(17)=cmp(Curimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code2(18)=cmp(Refimg(i+spoints(6,1),j+spoints(6,2)),Curimg(i,j));
code2(19)=cmp(Refimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code2(20)=cmp(Curimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code2(21)=cmp(Backimg(i+spoints(5,1),j+spoints(5,2)),Curimg(i,j));
code2(22)=cmp(Backimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code2(23)=cmp(Curimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code2(24)=cmp(Refimg(i+spoints(4,1),j+spoints(4,2)),Curimg(i,j));
code2(25)=cmp(Refimg(i,j),Curimg(i,j));
code2(26)=cmp(Backimg(i,j),Curimg(i,j));
%spectrum analysis
B1=SSTLBPSpectrum(code1);
B2=SSTLBPSpectrum(code2);
LocDensEst(a,b)=(B1+B2)/2;
b=b+1;
code1=zeros(1,26);
code2=zeros(1,26);
end
a=a+1;
b=1;
end
没有合适的资源?快使用搜索试试~ 我知道了~
spattempLBP.zip_neighbors_三值模式_像素时空信息_图像锯齿
![preview](https://csdnimg.cn/release/downloadcmsfe/public/img/white-bg.ca8570fa.png)
共1个文件
m:1个
![preview-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/scale.ab9e0183.png)
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 134 浏览量
2022-07-14
21:56:32
上传
评论
收藏 1KB ZIP 举报
温馨提示
时空局部二值模式,不仅采用的空间信息,也采用了时间信息,采用前后三帧图像的8领域像素与中心像素比较,大于为一,反之为0,然后分别采用锯齿形和螺旋形方式进行排列得到灰度值
资源推荐
资源详情
资源评论
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![doc](https://img-home.csdnimg.cn/images/20210720083327.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![application/x-rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
收起资源包目录
![package](https://csdnimg.cn/release/downloadcmsfe/public/img/package.f3fc750b.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
共 1 条
- 1
资源评论
![avatar-default](https://csdnimg.cn/release/downloadcmsfe/public/img/lazyLogo2.1882d7f4.png)
![avatar](https://profile-avatar.csdnimg.cn/3542a6830a364bf6a0adb34bfd72ce9d_weixin_42662605.jpg!1)
APei
- 粉丝: 65
- 资源: 1万+
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助
![voice](https://csdnimg.cn/release/downloadcmsfe/public/img/voice.245cc511.png)
![center-task](https://csdnimg.cn/release/downloadcmsfe/public/img/center-task.c2eda91a.png)
安全验证
文档复制为VIP权益,开通VIP直接复制
![dialog-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/green-success.6a4acb44.png)