基于加权最小二乘法的保边缘平滑滤波器（WLS），matlab程序自己编写，更直观

最新推荐文章于 2022-11-29 11:51:48 发布

残月飞雪

最新推荐文章于 2022-11-29 11:51:48 发布

阅读量6.9k

点赞数

CC 4.0 BY-SA版权

分类专栏：图像处理文章标签： WLS

本文链接：https://blog.youkuaiyun.com/majinlei121/article/details/50210939

图像处理专栏收录该内容

79 篇文章

订阅专栏

本文介绍了一种基于Matlab的图像处理算法实现过程，该算法通过计算图像梯度并应用非线性权重来达到平滑和锐化的效果。文章详细展示了如何使用差分、稀疏矩阵和其他数学操作来构建自定义图像处理算子。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

<span style="font-size:18px;">clear all;

g = double(imread('../images/barbara.png'))/255.0;

lambda=1.0;
alpha=1.2;

[m, n]=size(g);

smallNum = 0.0001;
A_x = diff(log(g+eps), 1, 2); 
A_x = 1./(abs(A_x).^alpha + smallNum);
A_x = padarray(A_x, [0 1], 'post');
A_x = A_x(:);
A_x=sparse(1:1:m*n,1:1:m*n,A_x(1:1:m*n)',m*n,m*n);

A_y = diff(log(g+eps), 1, 1);
A_y = 1./(abs(A_y).^alpha + smallNum);
A_y = padarray(A_y, [1 0], 'post');
A_y = A_y(:);
A_y=sparse(1:1:m*n,1:1:m*n,A_y(1:1:m*n)',m*n,m*n);

D_x=sparse([1:1:m*n,1:1:m*n-m],[1:1:m*n,m+1:1:m*n],[-ones(1,m*n),ones(1,m*n-m)],m*n,m*n);
D_y=sparse([1:1:m*n,2:1:m*n],[1:1:m*n,1:1:m*n-1],[-ones(1,m*n),ones(1,m*n-1)],m*n,m*n);
Lg=D_x'*A_x*D_x+D_y'*A_y*D_y;

II=sparse(1:1:m*n,1:1:m*n,ones(1,m*n),m*n,m*n);
A=II+lambda*Lg;

u=A\g(:);
u=reshape(u,m,n);


figure;imshow([g,u]);</span>