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

qUzalq

% 沿行插值
Sd = cumsum(diff,1);  % 行叠加
! \2 y0 F! g. D+ ^% 沿列插值
: l0 w$ _7 I- W6 bSd = cumsum(Sd,2);    % 列叠加
# |: s- B0 T! s( g" n+ t4 J, s! m为什么要进行叠加？
5 z' U; c$ R7 w2 T/ ^还有这一步也不是很明白，怎么就得到每个像素点的权值了
9 r$ |: U& K; p8 F8 T8 Z. h            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);

NingW

路过

mutougeda

分享这个源代码，应该有用
8 V) g: o9 F: E1 l$ @
, ^! A6 q- |  T& F3 @2 [% O, R% ?function DeNimg = Non_Local_Means(Nimg,PSH,WSH,Sigma)
4 g' ^, C# l+ ^" m/ l; g: c% N' q% Non_Local_Means滤波器
9 x" C0 o. M4 d3 e0 n%函数输入:
" C4 n) f# c; K) t9 R%         Nimg:  输入的图像矩阵 + 带噪声的
%         PSH:  扩展窗尺寸大小
%         WSH: 窗尺寸大小
% s! s9 I# F& h8 ]%         Sigma：方差
; b& R/ X, Q& i%函数输出:
3 q: t" ]7 X% v; D) q! J" F%         DeNimg: 重构滤波图像

if ~isa(Nimg,'double')
; {" o" V1 g' g3 Q5 J. B* S    Nimg = double(Nimg)/255;
/ T+ u. p' q4 {! zend

" F* ?% x4 c  `% v% 图像维数
[Height,Width] = size(Nimg);
u = zeros(Height,Width); % 初始化去噪图像矩阵
% n3 @( v) `7 F" }# IM = u; % 初始化权值矩阵
Z = M; % 初始化叠加权值 accumlated weights
% 避免边界效应
$ I- ?! e& @4 T( j2 g# m) t: d) s! c: ~$ mPP = padarray(Nimg,[PSH,PSH],'symmetric','both');
PW = padarray(Nimg,[WSH,WSH],'symmetric','both');
% padarray使用
% A =
' u% ^) T/ G3 N; O) Z' z%      1     3     4
+ N% S+ Y, b& b/ U%      2     3     4
%      3     4     5
% B = padarray(A, 2 * [1 1], 0, 'both')
6 U9 i/ ~1 c* g8 Q! p& _0 s. p" R%      0     0     0     0     0     0     0
/ `9 H' ]0 f" w3 n3 R! l%      0     0     0     0     0     0     0
: M( o/ t. J- r2 ~4 u- r' c) ]1 B%      0     0     1     3     4     0     0
0 h! o1 k! F) [# H%      0     0     2     3     4     0     0
%      0     0     3     4     5     0     0
2 r- A6 M4 G7 }- o) P3 {9 B%      0     0     0     0     0     0     0
/ b' i# b, Y7 D; w3 H5 V! \%      0     0     0     0     0     0     0
% 主循环
for dx = -WSH:WSH
    for dy = -WSH:WSH
        if dx ~= 0 || dy ~= 0
0 B# ^1 i; M8 O$ \3 R# H            Sd = integral_img(PP,dx,dy);  % 插值图像
            % 获取对应像素点的平方差矩阵
5 v8 `* Q4 u& a/ {$ R            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);      
            % 计算每一个像素点的权值
: m' B" n4 Q* U- x            w = exp(-SDist/(2*Sigma^2));
/ z; Q+ {9 m5 l2 r. ^  N            % 得到相应的噪声点
# y5 ~% l" E3 b            v = PW((WSH+1+dx)WSH+dx+Height),(WSH+1+dy)WSH+dy+Width));
- u4 y/ L6 d5 W4 f9 ?1 Q% o# e            % 更新去噪图像矩阵
. w# L) k) \/ z  M, I# N6 W* V3 y            u = u+w.*v;
0 W, }# _/ I+ E; r$ X            % 更新权值去噪图像矩阵
            M = max(M,w);
            % 更新叠加权值 accumlated weights
            Z = Z+w;
        end
    end
& ?. u2 f$ ]. _: _) E4 A% \end
1 K) T6 m; X, J$ E5 r0 k5 e3 D) f& o% 重构图像        
; C3 q1 [6 M9 Tf = 1;
u = u+f*M.*Nimg;
u = u./(Z+f*M);
DeNimg = u; % 重构去噪图像

; d( S4 W) R* J0 b9 F+ [function Sd = integral_img(v,dx,dy)
% 根据平方差，插值图像
1 V/ J5 {0 T$ |7 M. G' T5 f* m% 变换计算：tx = vx+dx; ty = vy+dy
# o# B2 H+ G; a0 i9 Yt = img_Shift(v,dx,dy);
% 平方差图像
' \- J1 J4 u" t% N* v, R& s. m! \diff = (v-t).^2;
% 沿行插值
+ e  W% z4 u9 ?- \" Y5 QSd = cumsum(diff,1);  % 行叠加
% 沿列插值
! J3 ~* V) N) O5 \& ESd = cumsum(Sd,2);    % 列叠加
9 ^1 c, l% W0 d$ @6 D+ `
! J+ C( j0 f* _) m( Wfunction t = img_Shift(v,dx,dy)
% 在xy坐标系下，进行图像变换操作
t = zeros(size(v));
type = (dx>0)*2+(dy>0);
switch type
, |5 V) g3 T5 ~9 S% h0 P7 h9 q" H    case 0 % dx<0,dy<0: 向右下方移动
        t(-dx+1:end,-dy+1:end) = v(1:end+dx,1:end+dy);
1 M' z8 ~5 E, H8 `. k    case 1 % dx<0,dy>0: 向左下方移动
+ k8 `, S" k. o- K3 _        t(-dx+1:end,1:end-dy) = v(1:end+dx,dy+1:end);
, a; M9 U) z0 T3 K+ B    case 2 % dx>0,dy<0: 向右上方移动
4 t" s7 i, u$ U. p; e, a        t(1:end-dx,-dy+1:end) = v(dx+1:end,1:end+dy);
8 s; p5 K2 Q0 A7 \( f' @    case 3 % dx>0,dy>0: 向左上方移动
        t(1:end-dx,1:end-dy) = v(dx+1:end,dy+1:end);
end