用matlab实现LBP算子

以下是使用MATLAB实现LBP算子的示例代码：

%% LBP算子实现
% 图像读取
img = imread('lena.jpg');
% 灰度化
gray_img = rgb2gray(img);
% 获取图像宽、高
[height, width] = size(gray_img);
% 定义LBP算子半径和邻域大小
radius = 1;
neighbors = 8;
% 初始化LBP图像
lbp_img = zeros(height, width);
% 计算LBP值
for i = radius+1 : height-radius
for j = radius+1 : width-radius
% 获取当前像素的邻域像素值
neighbor_values = zeros(1, neighbors);
for k = 1 : neighbors
theta = 2 * pi * (k-1) / neighbors;
x = round(i + radius * cos(theta));
y = round(j + radius * sin(theta));
neighbor_values(k) = gray_img(x, y);
end
% 计算LBP值
center_value = gray_img(i, j);
lbp_value = 0;
for k = 1 : neighbors
if neighbor_values(k) >= center_value
lbp_value = lbp_value + 2^(k-1);
end
end
% 将LBP值写入LBP图像
lbp_img(i, j) = lbp_value;
end
end
% 显示结果
subplot(1, 2, 1);
imshow(gray_img);
title('原始图像');
subplot(1, 2, 2);
imshow(lbp_img, []);
title('LBP图像');