Digital Signal Processing Reference
In-Depth Information
SigmaNoise = 5.
s1 = s + n* SigmaNoise
; Denoising using the undecimated WT on the sphere at 4sigma
Nsig = 4.
mrs_wtfilter, s1, fwt4, nsigma= Nsig, nbrscale=5, SigmaNoise=SigmaNoise
; Denoising using the curvelet transform
mrs_curfilter, s1, fct4, nsigma= Nsig, nbrscale=5, SigmaNoise=SigmaNoise
; Denoising using the combined denoising
mrs_cbfilter, s1, fcb4, nsigma= Nsig, nbrscale=5, SigmaNoise=SigmaNoise
; Display and write the figure to the disk
tvs, s, /log, tit='Synchrotron emission', png='fig_sync.png'
tvs, s1 > 30, /log, tit='Synchrotron emission + noise', $
png='fig_sync_noise5.png'
tvs , fwt4 > 30, /log, title='Undecimated Wavelet Denoising (4sigma)', $
png='fig_sync_wtfilter5.png'
tvs , fct4 > 30, /log, title='Curvelet Denoising (4sigma)',
$
png='fig_sync_curfilter5.png'
tvs , fcb4 > 30, /log, title='Combined Filtering (4sigma)',
png='fig_sync_cbfilter5.png'
tvs , s1- fwt4,
title='Residual undecimated Wavelet Denoising (4sigma)',$
png='fig_sync_resi_wtfilter5.png'
tvs , s1 - fct4,
title='Residual
curvelet Denoising (4sigma)',
$
png='fig_sync_resi_curfilter5.png'
tvs , s1
- fcb4,
title='Residual
combined Filtering (4sigma)',
$
png='fig_sync_resi_cbfilter5.png'
; Print the standard deviation (error) between the true image
; and the denoised images
print, 'Err WT = ', sigma(s-fwt4) , ', Err Cur
= ', sigma(s-fct4),
',
$
Err Comb = ', sigma(s-fcb4)
We find the outcome here to be
==>
Err WT = 1.25,
Err Cur
= 1.07,
Err Combined = 0.86
