局部直方图均衡化matlab代码

以下是在Matlab中实现局部直方图均衡化的示例代码：

% 读取图像
img = imread('your_image.jpg');

% 设置每个小块的大小
block_size = 32;

% 获取图像的大小
[img_height, img_width] = size(img);

% 计算图像中小块的数量
num_blocks_x = ceil(img_width / block_size);
num_blocks_y = ceil(img_height / block_size);
num_blocks = num_blocks_x * num_blocks_y;

% 初始化输出图像
out_img = zeros(size(img));

% 对于每个小块
for i = 1:num_blocks_y
    for j = 1:num_blocks_x
        % 计算小块的边界
        x_min = (j - 1) * block_size + 1;
        x_max = min(j * block_size, img_width);
        y_min = (i - 1) * block_size + 1;
        y_max = min(i * block_size, img_height);
        
        % 提取小块
        block = img(y_min:y_max, x_min:x_max);
        
        % 计算小块的直方图均衡化
        block_eq = histeq(block);
        
        % 将均衡化后的小块写回输出图像
        out_img(y_min:y_max, x_min:x_max) = block_eq;
    end
end

% 显示原始图像和输出图像
figure;
subplot(1,2,1); imshow(img); title('原始图像');
subplot(1,2,2); imshow(out_img); title('局部直方图均衡化后的图像');

这段代码将图像分成小块，对每个小块进行直方图均衡化，然后将均衡化后的小块拼接成输出图像。你可以根据需要调整块大小。