LBP

clear all; close all; clc; %I = im2double(rgb2gray(imread('index_1.bmp')));I = im2double(imread('cameraman.tif')); % imshow(I); [m, n] = size(I); J = zeros(m,n); % I1 = zeros(m,n);for i = 2 : m-1   for j = 2 : n-1      I1(i-1, j-1) = I(i-1, j-1) > I(i, j);       I1(i-1, j) = I(i-1, j) > I(i, j);       I1(i-1, j+1) = I(i-1, j+1) > I(i, j);       I1(i, j-1) = I(i, j-1) > I(i, j);       I1(i, j+1) = I(i, j+1) > I(i, j);       I1(i+1, j-1) = I(i+1, j-1) > I(i, j);       I1(i+1, j) = I(i+1, j) > I(i, j);       I1(i+1, j+1) = I(i+1, j+1) > I(i, j);      J(i, j) = I1(i-1, j-1) * 2^7 + I1(i-1, j) * 2^6 + I1(i-1, j+1) * 2^5 +  I1(i, j+1) * 2^4 + I1(i+1, j+1) * 2^3 + I1(i+1, j) * 2^2 + I1(i+1, j-1) * 2^1 + I1(i, j-1) * 2^0;    endendsubplot(121), imshow(I), title('image'); subplot(122), imshow(uint8(J)), title('image of LBP');