# Deblurred and estimate parameters of K blurred lines with high accuracy

• Generate lines with parameters , blurred in an image of size • Launch the minimization algorithm with parameters • From the solution , estimate the lines by Prony

## DATA GENERATION

W=65;               % image width (must be odd W=2M+1)
H=65;               % image height
K=3;                % number of lines
spread=1;           % std of the gaussian blur
noiselevel=0;       % level of noise
randomgen=0;        % boolean if lines are generated randomly of manually
rng(0);             % seed of random numbers generator
plotData=1;         % display images v, x and y
plotComp=0;         % display others comparaisons theo vs. emp

t_k=[pi/6 -pi/5 pi/16]; % array containing angles of lines
a_k=[255 255 255];      % array containing amplitude of lines
p_k=[10 0 -15];         % array containing offset of lines

data_generation;    % Generate the blurred image xstar
% of these K lines with additional
% noise y=xstar+randn(H,W)*noiselevel ## MINIMIZATION ALGORITHM

Nbiter=2000;            % number of iterations
rho=1.9;                % weigth coefficient
tau=150000000;          % primal step (sigma is then computed from tau)
c=sum(a_k./cos(t_k));   % theoritical atomic norm of each line
plotError=1;            % compute and display error
displayIter=0;          % display the counter iteration
algo=2;                 % choose the minimization algorithm to use

minimization_Chambolle; % Find by minimization an approximation wr
% of wstar which is the horizontal Fourier
% transform of the deblurred image sstar

Elapsed time is 258.346062 seconds.   ## IMAGE RECONSTRUCTION AND LINES ESTIMATION

Reconstruct the denoised image xstar from wstar:

resynthesis; Estimate lines parameters and display it

cstar=sum(a_k./cos(t_k));  % theoritical atomic norm of each line
coef=cstar/c;              % factor between experimental and true amplitude
prony;

thetamean = -0.62837     0.17872      0.5198
alphamean = 647.0155      642.9302      633.8648
eta = -0.179142     -14.9009      10.2207
(thetamean-t_k)/t_k = 8.7911e-05   -0.089811  -0.0072628
(alphamean*coef-a_k)/a_k = 0.012517   0.0061234  -0.0080631
eta-p_k = -0.17914    0.099138     0.22073 