main_interpolation.m

% Created by: Ph.D Nelson Eduardo Díaz Díaz
% Post-doctorador Pontíficia Universidad Católica de Valparaíso (PUCV)
% Date: 17 May 2021

% Solution using The 3-D N2-Queens problem
clear;
close all
clc
addpath(genpath('./utils'));
addpath(genpath('./results'));
addpath(genpath('./CodedApertureDesignOptimal'));
addpath(genpath('./dataset'))
addpath(genpath('./Reconstruction'))

%profile on
N = 256;  % spatial size
[a,b]=DDDRSNNP(N);

G = zeros(N,N);
mask = zeros(N,N,N);
temp0 =[];
temp1 =[];
a1 = a(1);
b1 = b(1);
x = 0;
for i = 1:N
    for j=1:N
        G(j,i) = mod(x+(i*a1) + (j*b1),N) +1;
        mask(j,i,G(j,i))=1;
    end
end

[T] = distan(G);
m = max(T(:));
minimum = min(T(:));


for k=1:1
    alldataset = {'wine2','wave','splash'};
    dataset = alldataset{k};
    load(dataset+".mat")
end

Y = sampling(data,mask,N); % compute compressive measurements with designed mask
figure(1)
imagesc(Y), title('Compressive Measurements using Designed Coded Apertures')
colormap('gray')

[J] = JointCodedApertures(Y,mask,N);

% Reconstruction using designed coded apertures
[Xrec] = interpolation(J);
disp('Interpolation of Temporal Rolling Shutter Compressive Measurements Captured using Designed Coded Apertures')
[p,s,r] = metrics(data,Xrec);

implay(Xrec),title('Interpolation with Designed Coded Apertures')

[mask] = random_Mask(N);
Y = sampling(data,mask,N); % compute compressive measurements with random mask
figure(2)
imagesc(Y), title('Compressive Measurements using Random Coded Apertures')
colormap('gray')

[J] = JointCodedApertures(Y,mask,N);
% Reconstruction using random coded apertures
[Xrec] = interpolation(J);
disp('Interpolation of Temporal Rolling Shutter Compressive Measurements Captured using Random Coded Apertures')
[p,s,r] = metrics(data,Xrec);
implay(Xrec),title('Interpolation with Random Coded Apertures')