局部均值分解matlab代码

以下是局部均值分解的MATLAB代码示例：

function [low_freq, high_freq] = local_mean_decomposition(img, block_size, sigma)

% 均值分解函数
% img: 输入图像
% block_size: 每个块的大小
% sigma: 高频成分的方差

% 将图像分块
[h, w] = size(img);
block_h = floor(h / block_size);
block_w = floor(w / block_size);
blocks = mat2cell(img, block_size * ones(1, block_h), block_size * ones(1, block_w));

% 对每个块进行均值分解
low_freq = zeros(h, w);
high_freq = zeros(h, w);
for i = 1:block_h
    for j = 1:block_w
        block = blocks{i, j};
        mean_val = mean(block(:));
        low_freq((i-1)*block_size+1:i*block_size, (j-1)*block_size+1:j*block_size) = mean_val;
        high_freq((i-1)*block_size+1:i*block_size, (j-1)*block_size+1:j*block_size) = block - mean_val;
    end
end

% 对高频成分进行平滑处理
h = fspecial('gaussian', 5, sigma);
high_freq = imfilter(high_freq, h, 'symmetric');

end

该函数将输入图像分成多个块，然后对每个块进行均值分解，将每个块的均值视为低频成分，将块减去均值得到高频成分。最后，对高频成分进行高斯平滑处理，以减少噪声。函数输出低频和高频成分。