基于双边滤波的图像处理，为什么会输出两个黑框？以下代码

Touuqu

基于双边滤波的图像处理，为什么会输出两个黑框？以下代码

clear all;
4 N. g. {: S% N0 w) c( w' R  T! Eclose all;
6 l0 d+ @& u6 Q6 N; s" q2 X2 eclc;

img=imread('t_2.jpg');
img=mat2gray(img);
7 y# o* U# U# u[m n]=size(img);
%imshow(img);

r=10; %模板半径
imgn=zeros(m+2*r+1,n+2*r+1);
img=imgn(r+1:m+r,r+1:n+r);
imgn(1:r,r+1:n+r)=img(1:r,1:n); %扩展上边界
imgn(1:m+r,n+r+1:n+2*r+1)=imgn(1:m+r,n:n+r); %扩展右边界
9 f( j" q+ F7 a, V8 q, _7 bimgn(m+r+1:m+2*r+1,r+1:n+2*r+1)=imgn(m:m+r,r+1:n+2*r+1); %扩展下边界
9 ]8 J, I' Y6 @. d. y2 Timgn(1:m+2*r+1,1:r)=imgn(1:m+2*r+1,r+1:2*r); %扩展左边界

sigma_d=2;
sigma_r=0.1;
[x,y] = meshgrid(-r:r,-r:r);
; G9 z$ X9 W6 s- f; {w1=exp(-(x.^2+y.^2)/(2*sigma_d^2)); %以距离作为自变量高斯滤波器

h=waitbar(0,'wait...');
) R! N* x/ {/ w! i. i' j% gfor i=r+1:m+r
for j=r+1:n+r
# r. \) H  V% i0 W# k0 T% N: X' E* Y
w2=exp(-(imgn(i-r:i+r,j-r:j+r)-imgn(i,j)).^2/(2*sigma_r^2)); %以周围和当前像素灰度差值作为自变量的高斯滤波器
w=w1.*w2;

    s=imgn(i-r:i+r,j-r:j+r).*w;    imgn(i,j)=sum(sum(s))/sum(sum(w));endwaitbar(i/m);

end
close(h)

figure;
imshow(mat2gray(imgn(r+1:m+r,r+1:n+r)));

7 P, I: P; p& a' b& z' F

relchhiclty

找相关的资料看看