BM3D BM3D完整代码（不是.p文件）

BM3D完整代码（不是.p文件），可以看到全部代码，供学习参考

% 图像分块，并且做变换，为找相似块做准备
% k:块大小，p：块移动步长，lambda2D，delta 收缩阈值
% block 返回的块，transform_block 变换后的块
% block2row_idxblock2col_idx 为保存的 块索引 与 块左上角在图像中坐标的 对应关系
function [blocktransform_blockblock2row_idxblock2col_idx] =im2block（imgkplambda2Ddelta）
    [rowcol]=size（img）;
    %这个阈值该用什么公式呢？？
    thres=lambda2D*delta*sqrt（2*log（row*col））;
    % r_num:行方向 上 应该有 多少个块
    r_num=floor（（row-k）/p）+1;
    c_num=floor（（col-k）/p）+1;
    block=zeros（kkr_num*c_num）;
    block2row_idx=[];
    block2col_idx=[];
    cnt=1;
    for i=0:1:r_num-1
        rs=1+i*p;
        for j=0:1:c_num-1
            cs=1+j*p;
            block（::cnt）=img（rs:rs+k-1cs:cs+k-1）;
            block2row_idx（cnt）=rs;
            block2col_idx（cnt）=cs;
            %该用什么变换呢？？
            tr_b=fft2（block（:

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-02-03 15:01  BM3D\
     文件        1107  2015-02-03 15:01  BM3D\im2block.m
     文件        1084  2015-02-03 14:57  BM3D\inv_transform_3d.m
     文件        6506  2015-02-03 14:56  BM3D\main.m
     文件        1928  2015-02-03 14:58  BM3D\search_similar_block.m
     文件        1265  2015-02-03 14:58  BM3D\testmatch.m
     文件         144  2015-02-03 14:59  BM3D\thres_shrink.m
     文件        2660  2015-02-03 15:01  BM3D\transform_3d.m