program cswt_tr_nsigma ! Uses use s2_types_mod use cswt_tr_mod use cswt_error_mod ! Calls call cswt_error call cswt_tr_const_multiply call cswt_tr_const_subtract call cswt_tr_free call cswt_tr_io_fits_write_wcoeff call cswt_tr_mask_apply call cswt_tr_mask_nonzero call cswt_tr_wcoeff_sigma call parse_options ! Variables character (len=S2_STRING_LEN) :: filename_in character (len=S2_STRING_LEN) :: filename_out character (len=S2_STRING_LEN) :: filename_mask type (cswt_tr) :: tr type (cswt_tr) :: tr_mask integer :: n_dilation integer :: fail real (kind=s2_sp), allocatable, dimension (:) :: sigma real (kind=s2_sp), allocatable, dimension (:) :: mean integer, allocatable, dimension (:,:) :: ncoeff_eff integer, allocatable, dimension (:) :: n_effective_samples integer :: ncoeff_tot ! Subroutines and functions subroutine parse_options () end program cswt_tr_nsigmaConvert wavelet coefficients to nsigma values. The coefficient map for each dilation is considered; the mean for each dilation is subtracted, before dividing by the sigma value computed for that dilation. Note that the mean and sigma are only computed over the non-masked pixels.
Usage: cswt_analysis
Author: J. D. McEwen (mcewen@mrao.cam.ac.uk)
Version: 0.1 - December 2004
character (len=S2_STRING_LEN) :: filename_in
character (len=S2_STRING_LEN) :: filename_out
character (len=S2_STRING_LEN) :: filename_mask
type (cswt_tr) :: tr
type (cswt_tr) :: tr_mask
integer :: n_dilation
integer :: fail
real (kind=s2_sp), allocatable, dimension (:) :: sigma
real (kind=s2_sp), allocatable, dimension (:) :: mean
integer, allocatable, dimension (:,:) :: ncoeff_eff
integer, allocatable, dimension (:) :: n_effective_samples
integer :: ncoeff_tot
subroutine parse_options () ! Calls: getArgument end subroutine parse_optionsParse the options passed when program called. ! @author J. D. McEwen (mcewen@mrao.cam.ac.uk)
Version: 0.1 - November 2004