多帧图像超分辨率重构代码和分析

代码包括低分辨率图像的获取和矫正和重构，已经解释和测试

clear all;
clc;
tic;
addpath（genpath（ ‘D:\已完成仿真\FRSR\重建初值‘））;%加入路径
addpath（genpath（ ‘D:\已完成仿真\FRSR\产生待重建的低分辨率图像‘））;%加入路径
number=16;
scale=4;
sd=[7.5035e-04;0.0034;0.0023;0.0039;0.0042;1.6262e-04;0.0044;0.0031];
sp=[0.0147;0.0188;0.0166;0.0058;0.0259;0.0088;0.0133;0.0245 ];
noisesd=[sd;sp];
delta=[1.2459 0.2544;1.0818 1.5652;3.2187 0.1709;...
          1.1261 0.8884;3.7900 4.4700;0.8561 1.4863;...
          1.3039 0.0585;4.3434 0.7524;1.6179 0.9584;...
         2.4087 4.1360;3.9513 1.0876;2.4135 1.9259;...
         3.3602 2.5350;4.9883 4.2239;3.2744 0.6202;...
         3.0191 1.4393];
 %phi=rand（number1）*5;
phi=[2.1312;4.0317;1.2249;2.5983;1.4762;0.4953;0.1746;...
          3.7099;0.1464;4.2639;1.4812;2.3967;3.7673;4.3981;...
          1.7816;1.8918];

y=cell（number1）;
f=cell（number1）;
%g=[0.00010.00020.00030.00040.00050.0030.0040.0050.0040.0030.00050.00040.00030.00020.0001];




x=imread（‘lena.tiff‘）;
%x=x（::3）;
[mn]=size（x）;
LR=imread（‘初值X0.bmp‘）;%读入待重建的低分辨率图像
%插值
s=imresize（LRscale‘bicubic‘）;%双三次插值，效果对比
%f0=deconvreg（spsfm*n*0.002）;
f{1}=double（s）;
%f{1}=imresize（f0scale‘bicubic‘）;
%s=imresize（f0scale‘bicubic‘）;
s1=imresize（LRscale‘bilinear‘）;%双线性插值
x=double（x）;

for i=1:number
    imageName = [ ‘LR1_‘ num2str（i） ‘.bmp‘];
   y{i}=double（imread（imageName））;
end%读入16幅同一场景下的LR图像
%Sc=generate_L（f{1}number）;
D_f=cell（number1）;
deD_f=cell（number1）;
sadd=zeros（size（x））;
iterative_number=15;%设定迭代次数
RMS=zeros（1iterative_number）;
PSNR=zeros（1iterative_number）;
sigma_d=12;sigma_r=0.25;
for i=1:iterative_number
    i
Sc=generate_L（f{i}numberdeltaphisdspscale）;%在原有的退化模型下退化
delta=-delta;phi=-phi;
  for k=1:number
      n=k;
     % y{k}=double（y{k}）;
     Sc{k}=255*Sc{k};
%size（Sc{k}）
%size（y{k}）
    D_f{k}=sign（Sc{k}-y{k}）;    
    deD_f{k}=huifu（D_f{k}nscalenoisesddeltaphi）;    
    sadd=sadd+deD_f{k};
  end
 sadd=sadd/number;
%frec=fswt（f{i}）;
frec=composite_cs（f{i}）;
 beta=80;
f{i+1}=f{i}-8*sadd-frec;%公式见文章中的迭代算法

   for m1=1:m
     for n1=1:n
        if f{i+1}（m1n1）<0
             f{i+1}（m1n1）=0;
        elseif f{i+1}（m1n1）>255
                  f{i+1}（m1n1）=255;
        end
      end
   end%将像素值统一到0-255
   
[MN]=size（f{i+1}）;

%f{i+1}（M-1:）=f{i+1}（M-2:）; 
%f{i+1}（M:）=f{i+1}（M-2:）; 
  %f{i} = mybilateral（f{i}sigma_dsigma_r）;
  [PSNR（i）RMS（i）]=indicator（xf{i}）
   figure;imshow（f{i}/255[]）;
   xlabel（[‘第‘num2str（i）‘次迭代重建图像 PSNR=‘num2str（PSNR（i））]）;
   delta=-delta;phi=-phi;
 
end
   


toc
figure
imshow（x/255[]）;xlabel（‘原始HR图像‘）
%双三次插值结果
[PSNR1RMS1]=indicator（xdouble（s））
figure
imshow（s[]）;xlabel（[‘双三次插值 PSNR=‘num2str（PSNR1）]）
%双线性插值结果
[PSNR2RMS2]=indicator（xdouble（s1））
figure
imshow（s1[]）;xlabel（[‘双线性插值 PSNR=‘num2str（PSNR2）]）

figure
imshow（LR[]）;xlabel（‘输入LR图像‘）
%保存评价指标
dict_path=[‘PSNR+RMSE.mat‘];
save（dict_path ‘PSNR‘ ‘RMS‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-12-28 09:49  像素上超分辨重构\
     文件         899  2012-04-09 09:20  像素上超分辨重构\Affine_trans.asv
     文件     5315639  2011-12-26 16:58  像素上超分辨重构\Fast and robust multiframe super resolution.pdf
     文件        3120  2012-07-25 10:43  像素上超分辨重构\chongjian.m
     文件         405  2012-04-10 20:35  像素上超分辨重构\composite_cs.m
     文件        2639  2012-04-09 08:36  像素上超分辨重构\dic_training.m
     文件        1487  2012-02-21 22:10  像素上超分辨重构\generate_L.m
     文件        1249  2012-07-25 10:24  像素上超分辨重构\get_LR.m
     文件         196  2011-08-31 17:53  像素上超分辨重构\goback.m
     文件         483  2012-04-10 15:17  像素上超分辨重构\huifu.m
     文件         310  2012-04-10 20:31  像素上超分辨重构\indicator.m
     文件        2886  2012-02-07 09:28  像素上超分辨重构\shift_fun1.m
     目录           0  2018-03-27 16:24  像素上超分辨重构\产生待重建的低分辨率图像\
     文件      298097  2011-07-26 21:38  像素上超分辨重构\产生待重建的低分辨率图像\1.gif
     文件         730  2012-04-09 09:45  像素上超分辨重构\产生待重建的低分辨率图像\generate_L1.m
     文件         427  2017-09-19 20:09  像素上超分辨重构\产生待重建的低分辨率图像\input_LR.m
     文件      249438  2011-07-26 21:38  像素上超分辨重构\产生待重建的低分辨率图像\lena.tiff
     文件         345  2012-07-25 10:38  像素上超分辨重构\产生待重建的低分辨率图像\para.mat
     文件       17462  2012-07-25 10:38  像素上超分辨重构\产生待重建的低分辨率图像\初值X0.bmp
     目录           0  2018-03-27 16:24  像素上超分辨重构\测试图像\
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\1.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\10.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\11.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\12.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\13.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\14.gif
     文件      298097  2011-07-26 21:38  像素上超分辨重构\测试图像\15.gif
     文件      257214  2011-07-26 21:38  像素上超分辨重构\测试图像\16.tiff
     文件      224348  2011-07-26 21:38  像素上超分辨重构\测试图像\17.tiff
     文件      398064  2011-07-26 21:38  像素上超分辨重构\测试图像\18.tiff
     文件      292638  2011-07-26 21:38  像素上超分辨重构\测试图像\19.tiff
............此处省略373个文件信息