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_nsigma
Convert 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_options
Parse the options passed when program called.
! @author J. D. McEwen (mcewen@mrao.cam.ac.uk)
Version: 0.1 - November 2004