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
Contents
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