C-Jonn 2017-05-08 03:31 采纳率: 0%
浏览 1382

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

clear all;
close all;
clc;

img=imread('t_2.jpg');
img=mat2gray(img);
[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); %扩展右边界
imgn(m+r+1:m+2*r+1,r+1:n+2*r+1)=imgn(m:m+r,r+1:n+2*r+1); %扩展下边界
imgn(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

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));

end
waitbar(i/m);

end
close(h)

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

  • 写回答

1条回答 默认 最新

  • 白萝卜。 2017-05-08 05:38
    关注
    评论

报告相同问题?

悬赏问题

  • ¥15 本题的答案是不是有问题
  • ¥15 关于#r语言#的问题:(svydesign)为什么在一个大的数据集中抽取了一个小数据集
  • ¥15 C++使用Gunplot
  • ¥15 这个电路是如何实现路灯控制器的,原理是什么,怎么求解灯亮起后熄灭的时间如图?
  • ¥15 matlab数字图像处理频率域滤波
  • ¥15 在abaqus做了二维正交切削模型,给刀具添加了超声振动条件后输出切削力为什么比普通切削增大这么多
  • ¥15 ELGamal和paillier计算效率谁快?
  • ¥15 蓝桥杯单片机第十三届第一场,整点继电器吸合,5s后断开出现了问题
  • ¥15 file converter 转换格式失败 报错 Error marking filters as finished,如何解决?
  • ¥15 Arcgis相交分析无法绘制一个或多个图形