求助：关于非局部均值滤波程序的问题

qUzalq

% 沿行插值
5 v: Y# i! @& T" ?Sd = cumsum(diff,1);  % 行叠加
% 沿列插值
6 P/ p( D: ?& I3 s% p! y5 zSd = cumsum(Sd,2);    % 列叠加
) D  w/ j5 L& u6 O为什么要进行叠加？
4 b* W' c  C/ {. {还有这一步也不是很明白，怎么就得到每个像素点的权值了
" [* D; h' Y$ j# }5 D. |            SqDist = Sd(PatchSizeHalf+1:end-PatchSizeHalf,PatchSizeHalf+1:end-PatchSizeHalf)+Sd(1:end-2*PatchSizeHalf,1:end-2*PatchSizeHalf)-Sd(1:end-2*PatchSizeHalf,PatchSizeHalf+1:end-PatchSizeHalf)-Sd(PatchSizeHalf+1:end-PatchSizeHalf,1:end-2*PatchSizeHalf);
) Z6 A% k  z( @6 b+ t$ f

NingW

路过

mutougeda

分享这个源代码，应该有用
( M  f0 K  A8 F2 ^# R8 u
5 O8 Q6 {8 m' e+ u' xfunction DeNimg = Non_Local_Means(Nimg,PSH,WSH,Sigma)
% Non_Local_Means滤波器
$ J$ J+ L3 {' g% w% J; p%函数输入:
%         Nimg:  输入的图像矩阵 + 带噪声的
) G  y; m$ ?% C9 a3 J%         PSH:  扩展窗尺寸大小
2 _( |5 a/ U' }- `* M%         WSH: 窗尺寸大小
7 {0 ^4 `" S9 w%         Sigma：方差
* b' x# @& y0 E. g+ k% Z$ o%函数输出:
" h- L5 W9 ~/ V# m8 }3 B, v%         DeNimg: 重构滤波图像

if ~isa(Nimg,'double')
    Nimg = double(Nimg)/255;
' ]* U' E" [2 Mend

* z4 L$ ?( e0 o, `/ N% 图像维数
[Height,Width] = size(Nimg);
7 B1 [, O: I$ Y4 @7 a; n+ Eu = zeros(Height,Width); % 初始化去噪图像矩阵
2 n; J. o* p& _- ?) vM = u; % 初始化权值矩阵
$ \  v, X5 @' {  NZ = M; % 初始化叠加权值 accumlated weights
% 避免边界效应
PP = padarray(Nimg,[PSH,PSH],'symmetric','both');
$ A$ T" F4 R$ z# @0 g3 Q* ~( oPW = padarray(Nimg,[WSH,WSH],'symmetric','both');
% padarray使用
% A =
# e0 v0 \0 J' Y%      1     3     4
%      2     3     4
. h& Q. o8 k# V. J; N+ T%      3     4     5
% B = padarray(A, 2 * [1 1], 0, 'both')
%      0     0     0     0     0     0     0
%      0     0     0     0     0     0     0
- [( M' s% c8 w" ~  d7 @%      0     0     1     3     4     0     0
%      0     0     2     3     4     0     0
4 B/ z5 k5 C4 Z( P4 Q, l1 s%      0     0     3     4     5     0     0
% _+ t0 X2 L0 h  N%      0     0     0     0     0     0     0
/ j$ c( I, n+ M3 P$ g9 x) l%      0     0     0     0     0     0     0
9 c3 W7 ^" [2 g. G* [- p, F2 B% 主循环
! F9 d  b" E/ S7 o2 {3 Wfor dx = -WSH:WSH
    for dy = -WSH:WSH
0 t/ c: j4 T/ y        if dx ~= 0 || dy ~= 0
            Sd = integral_img(PP,dx,dy);  % 插值图像
6 G+ e* E- V3 i  |            % 获取对应像素点的平方差矩阵
            SDist = Sd(PSH+1:end-PSH,PSH+1:end-PSH)+Sd(1:end-2*PSH,1:end-2*PSH)-Sd(1:end-2*PSH,PSH+1:end-PSH)-Sd(PSH+1:end-PSH,1:end-2*PSH);      
            % 计算每一个像素点的权值
            w = exp(-SDist/(2*Sigma^2));
            % 得到相应的噪声点
" z' i# h- r  J& A            v = PW((WSH+1+dx)WSH+dx+Height),(WSH+1+dy)WSH+dy+Width));
, X9 s+ _: D+ ^7 ?* Y* f            % 更新去噪图像矩阵
( ~2 n! u4 U! @: R            u = u+w.*v;
            % 更新权值去噪图像矩阵
            M = max(M,w);
, k2 s& D4 f4 \            % 更新叠加权值 accumlated weights
            Z = Z+w;
        end
: u) g, d- I3 @4 M" b    end
9 M* E8 Z( S: Z: u, f9 y+ x: i2 nend
% 重构图像        
f = 1;
u = u+f*M.*Nimg;
0 j; q0 t' u  c$ K7 wu = u./(Z+f*M);
9 I2 N0 H* k( y  |DeNimg = u; % 重构去噪图像
8 k" E' [2 t  K- L
function Sd = integral_img(v,dx,dy)
/ A# A  o: p8 o, f% 根据平方差，插值图像
% 变换计算：tx = vx+dx; ty = vy+dy
( v& d+ m2 L' et = img_Shift(v,dx,dy);
& B' b# s; n. `9 B+ d. U% 平方差图像
0 _$ e& N& d- P0 s$ Pdiff = (v-t).^2;
, M5 ~1 |+ F+ V  e, q! u- O' E6 X% 沿行插值
: t' o7 X) j# R: g$ d  H) w0 CSd = cumsum(diff,1);  % 行叠加
# u& B" X% v" N1 I" ^1 w7 l% 沿列插值
Sd = cumsum(Sd,2);    % 列叠加

function t = img_Shift(v,dx,dy)
% 在xy坐标系下，进行图像变换操作
t = zeros(size(v));
type = (dx>0)*2+(dy>0);
switch type
    case 0 % dx<0,dy<0: 向右下方移动
7 Z$ R1 D( C) [2 s* H        t(-dx+1:end,-dy+1:end) = v(1:end+dx,1:end+dy);
. m3 G" P$ a  Z7 ]$ R    case 1 % dx<0,dy>0: 向左下方移动
        t(-dx+1:end,1:end-dy) = v(1:end+dx,dy+1:end);
/ d% Q5 c5 b2 s    case 2 % dx>0,dy<0: 向右上方移动
        t(1:end-dx,-dy+1:end) = v(dx+1:end,1:end+dy);
    case 3 % dx>0,dy>0: 向左上方移动
        t(1:end-dx,1:end-dy) = v(dx+1:end,dy+1:end);
end