hfs_blanking

Back to the main HfS page

hfs_blanking is a procedure to flag the fit values in the output out and FITS files from hfs_cube, according to any line parameter value, error, or |value|/error.

The procedure reads the values and errors of the line fit parameters from the out files, and you can choose which line parameter to use for blanking the fit. Several passes can be done, using different line fit parameters or blanking criteria.
The output files are the blanked files <parfile>_comp#_blank.out and <parfile>_<fit_par>_comp#_blank.fits, without the fit values flagged out by hfs_blanking.

Example of a hfs_blanking run:

$ hfs_blanking i22test.par ________________________________________________ HfS_blanking. Blanking fit values from hfs_cube Robert Estalella, 2017/05_________________________________________________ Par file: i22test.par dim= [512, 512, 4, 9] Storage allocated successfully for arrays Reading out file i22test_comp1.out Reading out file i22test_comp2.out
Number of pixels Total: 49 Fitted: 49 Comp 1: 48 Comp 2: 49
HfS Blanking Menu ____________________________________ 0. Quit without saving 1. Save blanked files and exit 2. Delta_V 3. V_lsr 4. A*tau_m 5. tau_m Choose option or parameter to use for blanking (0-5): 2
Delta_V Value Error |Val|/Err Comp Min Max Min Max Min Max 1 2.50E-02 2.94E+00 5.94E-02 6.18E+00 2.30E-02 8.33E-01 2 2.63E-02 2.15E+00 7.88E-02 6.69E+00 4.99E-02 1.81E+00
Blanking Criterion Menu_______________________________ 0. Back 1. Blanking for Error > cutoff 2. Blanking for Error < cutoff 3. Blanking for |Value|/Error < cutoff 4. Blanking for Value > cutoff 5. Blanking for Value < cutoff Choose option (0-5): 1 Blanking (1). Enter cutoff: 1.0
Blanking pass 1 Blanked Left Pixels: 4 45 Comp 1: 14 34 Comp 2: 12 37
Delta_V Value Error |Val|/Err Comp Min Max Min Max Min Max 1 2.50E-02 3.98E-01 5.94E-02 9.04E-01 5.70E-02 5.38E-01 2 2.63E-02 1.26E+00 7.88E-02 9.25E-01 4.99E-02 1.81E+00
Blanking Criterion Menu_______________________________ 0. Back 1. Blanking for Error > cutoff 2. Blanking for Error < cutoff 3. Blanking for |Value|/Error < cutoff 4. Blanking for Value > cutoff 5. Blanking for Value < cutoff Choose option (0-5): 0 HfS Blanking Menu ____________________________________ 0. Quit without saving 1. Save blanked files and exit 2. Delta_V 3. V_lsr 4. A*tau_m 5. tau_m Choose option or parameter to use for blanking (0-5): 1
Writing blanked out files i22test_comp1_blank.out i22test_comp2_blank.out
Writing blanked fits files maps/i22test_dvline_comp1_blank.fits maps/i22test_vlsr_comp1_blank.fits maps/i22test_ataum_comp1_blank.fits maps/i22test_taum_comp1_blank.fits maps/i22test_dvline_comp2_blank.fits maps/i22test_vlsr_comp2_blank.fits maps/i22test_ataum_comp2_blank.fits maps/i22test_taum_comp2_blank.fits

Input

<parfile>.par, the same parameter file used as input to hfs_cube. From the information in this file, hfs_blanking reads the corresponding <parfile>_comp#.out file for each velocity component (created by hfs_cube), with the parameters of the fitted spectra.

Output

<parfile>_comp#_blank.out: an ASCII file for each velocity component, without the fit values flagged out by hfs_blanking.
maps/<parfile>_dvline_comp#_blank.fits, maps/<parfile>_vlsr_comp#_blank.fits, maps/<parfile>_ataum_comp#_blank.fits, maps/<parfile>_taum_comp#_blank.fits: FITS files for each component in folder maps, with non-flagged pixels, and NaN for pixels flagged out by hfs_blanking. Each FITS file has two planes, the first plane with the values of the parameter, and the second plane with the uncertainties.