clear;
close all;
load Data_Exo_2/SG1.mat

B = log(DataMod(:));
A = [ -Data(:) ones(size(Data(:))) ];

C = [ A B ];
[U, S , V] = svd(C);

param = [ V(1,3) V(2,3) ] / -V(3,3);

load Data_Exo_2/ImSG1.mat

original = ( -log(ImMod) + param(2) ) / param(1);

imshow(I)
title("image de base")
figure
imshow(original)
title("image calculée")

RMSE = mean(mean(sqrt((original - I).^2)))