一、图片非局部相似部分Get_PG

此函数实现提取相似块的功能：

function   [Px, Px0] =  Get_PG( im,win, ps ,nlsp,step,delta) %返回值为方差与目标块小的相似块和普通相似块

[h, w, ch]  =  size(im);
S         =  win;
maxr         =  h-ps+1;
maxc         =  w-ps+1;
r         =  1:step:maxr;%目标块不是挨像素取的，一次隔三个像素
r         =  [r r(end)+1:maxr];
c         =  1:step:maxc;
c         =  [c c(end)+1:maxc];
X = zeros(ps^2*ch,maxr*maxc,'double');
Px0 = [];
Px = [];
if nlsp ==1
    k    =  0;
    for l = 1:ch
        for i  = 1:ps
            for j  = 1:ps
                k    =  k+1;
                blk     =  im(r-1+i,c-1+j,l);
                Px(k,:) =  blk(:)';
            end
        end
    end
else
    k    =  0;
    for l = 1:ch
        for i  = 1:ps
            for j  = 1:ps
                k    =  k+1;
                blk  =  im(i:end-ps+i,j:end-ps+j,l);%可能是为了提高运算速度，将需要的数据先存储
                X(k,:) =  blk(:)';%第一列的108个数据即为im（6，6，3）第二个，列移动一格（图）
            end
        end
    end
end
% Index image
Index     =   (1:maxr*maxc);
Index    =   reshape(Index, maxr, maxc);
N1    =   length(r);
M1    =   length(c);
blk_arr   =  zeros(nlsp, N1*M1 );%每列记录10个最相似的编号，行是目标块数目
for  i  =  1 :N1
    for  j  =  1 : M1
        row     =   r(i);%第一行
        col     =   c(j);%第一列
        off     =  (col-1)*maxr + row;%是目标块编号，块编号在X中是按列来的，所以此处变一列，即加
        off1    =  (j-1)*N1 + i;
        
        rmin    =   max( row-S, 1 );
        rmax    =   min( row+S, maxr );
        cmin    =   max( col-S, 1 );
        cmax    =   min( col+S, maxc );
        
        idx     =   Index(rmin:rmax, cmin:cmax);%idx为窗口位置（索引）
        idx     =   idx(:);
        neighbor       =   X(:,idx);
        seed       =   X(:,off);
        dis = sum(bsxfun(@minus,neighbor, seed).^2,1);%对应相减，平方，列求和
        [~,ind]   =  sort(dis);
        indc        =  idx( ind( 1 : nlsp ) );
        blk_arr(:,off1)  =  indc; %indc为最相似十个块的索引
        X_nl = X(:,indc); % or X_nl = neighbor(:,ind( 1 : nlsp ));
        % Removes DC component from image patch group
        DC = mean(X_nl,2);
        X_nl = bsxfun(@minus, X_nl, DC);  %x_nl为10个最相似的块
        % Select the smooth patches
        sv=var(X_nl);
        if max(sv) <= delta
            Px0 = [Px0 X_nl];
        else
            Px = [Px X_nl];
        end
    end
end