matlab实现的稀疏表示的遥感图像融合

matlab实现的稀疏表示的遥感图像融合

close all
% n is the original signal length
n = 2^12;

% k is number of observations to make
k = 2^10;

% number of spikes to put down
% n_spikes = floor（.01*n）;
n_spikes = 160;


% random +/- 1 signal
f = zeros（n1）;
q = randperm（n）;
f（q（1:n_spikes）） = sign（randn（n_spikes1））;
%f（q（1:n_spikes）） = randn（n_spikes1）;

% measurement matrix
disp（‘Creating measurement matrix...‘）;
R = randn（kn）;

% orthonormalize rows
R = orth（R‘）‘;

if n == 8192  
   % in this case we load a precomputed
   % matrix to save some time
   load Rmatrix_2048_8192.mat
end
%
disp（‘Finished creating matrix‘）;

hR = @（x） R*x;
hRt = @（x） R‘*x;

% noisy observations
sigma = 0; 
y = hR（f） + sigma*randn（k1）;

% regularization parameter
tau = 0.005*max（abs（R‘*y））;

first_tau_factor = 0.8*（max（abs（R‘*y））/tau）;
steps = 5

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

     文件       6148  2007-12-04 11:46  GPSR_5.0\.DS_Store

     文件      65126  2007-01-25 13:54  GPSR_5.0\Camera.tif

     文件       7046  2007-09-23 20:49  GPSR_5.0\demo_continuation.m

     文件       5849  2007-07-13 16:40  GPSR_5.0\demo_image_deblur.m

     文件       8489  2007-10-19 16:40  GPSR_5.0\figures_1_2_3.m

     文件       5225  2007-09-12 18:21  GPSR_5.0\figure_4.m

     文件       5276  2007-12-04 11:46  GPSR_5.0\figure_5.m

     文件       2377  2007-07-12 18:41  GPSR_5.0\figure_6.m

     文件      21851  2007-12-04 11:26  GPSR_5.0\GPSR_Basic.m

     文件      21647  2007-12-04 11:08  GPSR_5.0\GPSR_BB.m

     文件      12630  2007-09-06 15:10  GPSR_5.0\IST.m

     文件       8569  2007-06-16 15:41  GPSR_5.0\l1_ls.m

     文件       3178  2007-12-05 19:06  GPSR_5.0\README.txt

     目录          0  2015-12-18 18:51  GPSR_5.0

----------- ---------  ---------- -----  ----

               173411                    14