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

Touuqu

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

clear all;
" }6 v% d1 V+ Z# a/ ?; r* Rclose all;
. s* s& A4 ]- N# f. \clc;

img=imread('t_2.jpg');
  ^  k3 x3 ~1 x  m; aimg=mat2gray(img);
0 O' L% E+ b1 c9 @: N[m n]=size(img);
" b  s5 f4 k, P( Q%imshow(img);

r=10; %模板半径
) w* j+ V( i! ?3 A2 Dimgn=zeros(m+2*r+1,n+2*r+1);
img=imgn(r+1:m+r,r+1:n+r);
. h) ?, \) u  S' f) t  fimgn(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); %扩展右边界
imgn(m+r+1:m+2*r+1,r+1:n+2*r+1)=imgn(m:m+r,r+1:n+2*r+1); %扩展下边界
( u" ~; q6 P. |' e3 J6 z- Oimgn(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);
w1=exp(-(x.^2+y.^2)/(2*sigma_d^2)); %以距离作为自变量高斯滤波器

h=waitbar(0,'wait...');
for i=r+1:m+r
for j=r+1:n+r
: I1 K2 c$ |" ~& P9 _
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
0 K' H  a0 h. k6 \# D; Gclose(h)

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

0 x% E6 {! V/ m  y  m

4 y+ x* o0 i! D6 Y  |5 A

" Y) M& K; y. m% D3 m  c8 i) [, n6 C

relchhiclty

找相关的资料看看