! @(#)wa_lucy.prg	17.1.1.1 (ES0-DMD) 01/25/02 17:21:33
! @(#) @(#)wa_lucy.prg 17.1.1.1 02/01/25 ESO @(#)
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++
!.IDENTIFICATION: wa_lucy.prg
!.PURPOSE:    deconvolve  an image with a regularization in the wavelet space
!.USE:        LUCY/WAVE Image_In Psf Imag_Out [residual]  [Nb_Scales]
!                [N_Sigma, Noise] [Epsilon] [Max_Iter]
!.AUTHOR:     Jean-Luc Starck
!.VERSION:    02/01/25
! KEYWORDS:   Wavelet transform, Multiresolution analysis, restoration
! ----------------------------------------------------------------------

DEFINE/PARAM P1 ? IMAGE "Enter Image:  "
DEFINE/PARAM P2 ? IMAGE "Enter Psf:  "
DEFINE/PARAM P3 ? IMAGE "Enter Result:  " 
DEFINE/PARAM P4 "residual" IMAGE "Enter Residual:  " 
DEFINE/PARAM P5 4 N "Enter the number of scales:  " 1,20
DEFINE/PARAM P6 4.,0 N "Enter N_sigma, Noise  " 
DEFINE/PARAM P7 0.001 N "Convergence parameter:"
DEFINE/PARAM P8 100 N "Enter the maximun iteration number :  " 1,10000

WRITE/KEYW IN_A/C/1/60 {P1}
WRITE/KEYW IN_B/C/1/60 {P2}
WRITE/KEYW OUT_A/C/1/60  {P3}
WRITE/KEYW OUT_B/C/1/60  {P4}

WRITE/KEYW INPUTI/I/1/1  {P5}
WRITE/KEYW INPUTI/I/2/1  {P8}

WRITE/KEYW INPUTR/R/1/2  {P6}
WRITE/KEYW INPUTR/R/3/1  {P7}

RUN CON_EXE:wa_lucy.exe