图像分块BCS稀疏表示与重建

该程序的功能是实现图像的稀疏分解并重建图像，采用压缩感知理论，利用小波分析，在小波域分块观测，可以用很少的观测值就能重建原图像。程序包中已经整合了小波包和测试图像，没有多余文件，程序经过亲测并修改，能在MATLAB中直接运行得到很好的结果，直接运行main_msbcsspl，速度非常快！

%
% function X = BlockDWT（x block_sizes）
%
%   This function performs wavelet-domain block-based sampling of the
%   spatial-domain image x. block_sizes gives the size of blocks in
%   each level of the DWT.
%

%
% MS-BCS-SPL: Multiscale Block Compressed Sensing - Smooth Projected Landweber
% Copyright （C） 2012  James E. Fowler
% http://www.ece.mstate.edu/~fowler
%
% This program is free software; you can redistribute it and/or
% modify it under the terms of the GNU General Public License
% as published by the Free Software Foundation; either version 2
% of the License or （at your option） any later version.
%
% This program is distributed in the hope that it will be useful
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program; if not write to the Free Software
% Foundation Inc. 675 Mass Ave Cambridge MA 02139 USA.
%


function X = BlockDWT（x block_sizes）

num_levels = length（block_sizes）;

[num_rows num_cols] = size（x）;

x = waveletcdf97（x num_levels）;

subband_num_rows = num_rows / 2^num_levels;
subband_num_cols = num_cols / 2^num_levels;

X{1} = x（1:subband_num_rows 1:subband_num_cols）;

for level = 1:num_levels
    
    block_size = block_sizes（level）;
    
    subband1 = x（1:subband_num_rows ...
        （subband_num_cols + 1）:（2 * subband_num_cols））;
    X{level + 1}{1} = im2col（subband1 [block_size block_size] ‘distinct‘）;
    subband2 = x（（subband_num_rows + 1）:（2 * subband_num_rows） ...
        1:subband_num_cols）;
    X{level + 1}{2} = im2col（subband2 [block_size block_size] ‘distinct‘）;
    subband3 = x（（subband_num_rows + 1）:（2 * subband_num_rows） ...
        （subband_num_cols + 1）:（2 * subband_num_cols））;
    X{level + 1}{3} = im2col（subband3 [block_size block_size] ‘distinct‘）;
    
    subband_num_rows = subband_num_rows * 2;
    subband_num_cols = subband_num_cols * 2;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     786486  2007-03-21 15:41  MSBCS\Images\baboon.bmp

     文件      66614  2018-06-13 21:23  MSBCS\Images\Barbara256.bmp

     文件      66614  2006-12-03 11:15  MSBCS\Images\lena256.bmp

     文件     196662  2018-05-25 20:47  MSBCS\Images\Mandrill256.bmp

     文件       2058  2018-05-28 23:29  MSBCS\MSBCS\BlockDWT.m

     文件       6580  2018-05-29 00:48  MSBCS\MSBCS\DDWT_Decoder.m

     文件       3010  2018-07-10 02:58  MSBCS\MSBCS\Decoder.m

     文件       1703  2012-08-17 22:43  MSBCS\MSBCS\Encoder.m

     文件       3984  2018-07-11 00:12  MSBCS\MSBCS\GenerateProjection.m

     文件       1805  2012-08-17 22:43  MSBCS\MSBCS\GetSubrates.m

     文件       2149  2018-05-28 23:15  MSBCS\MSBCS\GPSR_Encoder.m

     文件       3746  2018-06-16 08:59  MSBCS\MSBCS\GPSR_GenerateProjection.m

     文件       2184  2018-05-28 23:12  MSBCS\MSBCS\InverseBlockDWT.m

     文件       3637  2019-02-13 22:06  MSBCS\MSBCS\main_msbcsspl.m

     文件       1329  2012-08-14 02:54  MSBCS\MSBCS\PSNR.m

     文件       2733  2012-08-17 22:43  MSBCS\MSBCS\README

     文件       1312  2012-08-14 02:54  MSBCS\MSBCS\RMS.m

     文件      30728  2006-07-09 01:34  MSBCS\waveletcdf97\demo\palm.jpg

     文件       5604  2006-07-10 00:33  MSBCS\waveletcdf97\demo\waveletcdf97_demo.m

     文件      15973  2006-07-10 01:17  MSBCS\waveletcdf97\doc\waveletcdf97.html

     文件       2429  2006-07-09 21:26  MSBCS\waveletcdf97\doc\waveletcdf97_01.png

     文件      31276  2006-07-09 22:19  MSBCS\waveletcdf97\doc\waveletcdf97_02.jpg

     文件       1766  2006-07-09 22:54  MSBCS\waveletcdf97\doc\waveletcdf97_03.png

     文件       2760  2006-07-09 23:12  MSBCS\waveletcdf97\doc\waveletcdf97_04.png

     文件       3403  2006-07-10 00:33  MSBCS\waveletcdf97\test\test_main.m

     文件       5844  2006-07-10 01:13  MSBCS\waveletcdf97\waveletcdf97.m

     文件       6148  2010-02-09 10:37  MSBCS\WaveletSoftware\.DS_Store

     文件        886  2003-07-17 16:31  MSBCS\WaveletSoftware\afb.m

     文件        870  2018-05-29 00:11  MSBCS\WaveletSoftware\afb2D.m

     文件        988  2018-05-29 00:56  MSBCS\WaveletSoftware\afb2D_A.m

............此处省略85个文件信息