デュアルフィルタmatlab

ダイレクトコード

function img = myBilateralFilter(Image, kerSize, delta)
% Image      
% kerSize      
% delta    
% img     

%%
% c，r    kerSize          
c = floor(kerSize(1)/2);
r = floor(kerSize(2)/2);

%       
padImage = padarray(Image, [c, r], 'symmetric');
img = zeros(size(Image));

% G distance weight
siz = (kerSize-1)/2;
std = delta;
[x, y] = meshgrid(-siz(2): siz(2), -siz(1): siz(1));
arg    = -(x.*x + y.*y)/(2*std*std);
G      = exp(arg);
G(G

関数呼び出し

delta = 5;
bilateralImg = myBilateralFilter(Image, [5, 5], delta);
bilateralImg = uint8(bilateralImg);