该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
function X1 = LBP(Resimage)
% Returns a rotation invariant LBP (uniform patterns) histogram (10 bins)
% of picture X.
% the size of picture X must be at least 3x3 pixels
w1 = (1/sqrt(2))^2;%根号1/2平方
w2 = (1-1/sqrt(2))*(1/sqrt(2));
[sy sx] = size(Resimage);
Xi = zeros(sy+2,sx+2);%xi为变量产生全0矩阵扩大图像,便于处理
Xi(2:sy+1,2:sx+1) = Resimage;%赋值为0
Xi2 = zeros(sy+2,sx+2);
Xi4= zeros(sy+2,sx+2);
Xi6 = zeros(sy+2,sx+2);
Xi8 = zeros(sy+2,sx+2);
p1 = zeros(sy+2,sx+2);
p2 = zeros(sy+2,sx+2);
p3 = zeros(sy+2,sx+2);
p4 = zeros(sy+2,sx+2);
p5 = zeros(sy+2,sx+2);
p6 = zeros(sy+2,sx+2);
p7 = zeros(sy+2,sx+2);
p8 = zeros(sy+2,sx+2);
p9 = zeros(sy+2,sx+2);
p1(3:sy+2,3:sx+2) = Resimage;
p2(3:sy+2,2:sx+1) = w2*double(Resimage);
p3(3:sy+2,1:sx) = Resimage;
p4(2:sy+1,3:sx+2) = w2*double(Resimage);
p5(2:sy+1,2:sx+1) = (1-1/sqrt(2))^2*double(Resimage);
p6(2:sy+1,1:sx) = w2*double(Resimage);
p7(1:sy,3:sx+2) = Resimage;
p8(1:sy,2:sx+1) = w2*double(Resimage);
p9(1:sy,1:sx) = Resimage;
Xi1 = w1*p1+ p2+p4 + p5 + 0.000001; %从右到下Xi1 to the right and down from X
Xi2(3:sy+2,2:sx+1) = Resimage;
Xi3 = w1*p3 + p2 + p6 + p5 + 0.000001;
Xi4(2:sy+1,1:sx) = Resimage;
Xi5 = w1*p9 + p8 + p6 + p5 + 0.000001;
Xi6(1:sy,2:sx+1) = Resimage;
Xi7 = w1*p7 + p8 + p4 + p5 + 0.000001;
Xi8(2:sy+1,3:sx+2) = Resimage;
Xi = (Xi4>=Xi)+2*(Xi5>=Xi)+4*(Xi6>=Xi)+8*(Xi7>=Xi)+16*(Xi8>=Xi)+32*(Xi1>=Xi)+64*(Xi2>=Xi)+128*(Xi3>=Xi);
X1 = Xi(3:sy,3:sx);%取3*3间隔
end