diff --git a/plots/NGC1068_x274020/NGC1068_FOC.png b/plots/NGC1068_x274020/NGC1068_FOC.png index d9d86bb..6ba0955 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC.png and b/plots/NGC1068_x274020/NGC1068_FOC.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_IQU.png new file mode 100644 index 0000000..26099eb Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_P.png b/plots/NGC1068_x274020/NGC1068_FOC_P.png index ac7441b..e550b10 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_P_err.png index 1bfffef..d32c4d8 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_SNRi.png index 008e286..24fd7d8 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_SNRp.png index 6ec23ff..41eb676 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_center_image.png b/plots/NGC1068_x274020/NGC1068_FOC_center_image.png index efca298..d026d69 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_center_image.png and b/plots/NGC1068_x274020/NGC1068_FOC_center_image.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png index 76015d6..d72f600 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_IQU.png new file mode 100644 index 0000000..286280d Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png index 36777ea..7f2bc3b 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png index c6670de..81f7ce0 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png index 5b03843..8ef875e 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png index 28c5217..0bb8a03 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png index cc68add..0ad1526 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_IQU.png new file mode 100644 index 0000000..2985fb2 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png index 82343b1..03bd9f8 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png index 236a8d6..0f50954 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png index 2be854c..5feb260 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png index b3f78cb..d53ae0e 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot.png index ad5271d..0f15e5b 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_IQU.png new file mode 100644 index 0000000..5a10037 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P.png index 0cfc9e0..efa341a 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P_err.png index b5c05c1..64ed298 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRi.png index 721ab03..033f729 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRp.png index c551f0d..78ae955 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM010_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot.png new file mode 100644 index 0000000..ff74553 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_IQU.png new file mode 100644 index 0000000..7e749fe Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P.png new file mode 100644 index 0000000..d888469 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P_err.png new file mode 100644 index 0000000..75bd958 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRi.png new file mode 100644 index 0000000..56e154a Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRp.png new file mode 100644 index 0000000..f65b260 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_FWHM020_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot.png index 1f8ff24..33e1e67 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_IQU.png new file mode 100644 index 0000000..f77091d Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P.png index ea61e83..4025623 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P_err.png index 0254a86..66d8a7b 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRi.png index a44c5cc..b2f1cd0 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRp.png index c68744f..94050e1 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM010_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot.png new file mode 100644 index 0000000..72a3801 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_IQU.png new file mode 100644 index 0000000..e14375c Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P.png new file mode 100644 index 0000000..7ef4fd5 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P_err.png new file mode 100644 index 0000000..218a7d3 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRi.png new file mode 100644 index 0000000..29e47ce Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRp.png new file mode 100644 index 0000000..9d2798e Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_gaussian_after_FWHM020_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_rot.png index ff7e0a9..35729c2 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_rot_IQU.png new file mode 100644 index 0000000..84a75b7 Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_rot_IQU.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_rot_P.png index eece091..e929ad5 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_rot_P_err.png index c035c11..c1f2aab 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRi.png index 4167b85..2baf60d 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRp.png index acf2599..5360236 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_rot_SNRp.png differ diff --git a/src/FOC_reduction.py b/src/FOC_reduction.py index 317a256..4967b53 100755 --- a/src/FOC_reduction.py +++ b/src/FOC_reduction.py @@ -16,20 +16,20 @@ import lib.plots as proj_plots #Functions for plotting data def main(): ##### User inputs ## Input and output locations -# globals()['data_folder'] = "../data/NGC1068_x274020/" -# infiles = ['x274020at.c0f.fits','x274020bt.c0f.fits','x274020ct.c0f.fits', -# 'x274020dt.c0f.fits','x274020et.c0f.fits','x274020ft.c0f.fits', -# 'x274020gt.c0f.fits','x274020ht.c0f.fits','x274020it.c0f.fits'] -# globals()['plots_folder'] = "../plots/NGC1068_x274020/" + globals()['data_folder'] = "../data/NGC1068_x274020/" + infiles = ['x274020at.c0f.fits','x274020bt.c0f.fits','x274020ct.c0f.fits', + 'x274020dt.c0f.fits','x274020et.c0f.fits','x274020ft.c0f.fits', + 'x274020gt.c0f.fits','x274020ht.c0f.fits','x274020it.c0f.fits'] + globals()['plots_folder'] = "../plots/NGC1068_x274020/" # globals()['data_folder'] = "../data/NGC1068_x14w010/" # infiles = ['x14w0101t_c0f.fits','x14w0102t_c0f.fits','x14w0103t_c0f.fits', # 'x14w0104t_c0f.fits','x14w0105p_c0f.fits','x14w0106t_c0f.fits'] # globals()['plots_folder'] = "../plots/NGC1068_x14w010/" - globals()['data_folder'] = "../data/3C405_x136060/" - infiles = ['x1360601t_c0f.fits','x1360602t_c0f.fits','x1360603t_c0f.fits'] - globals()['plots_folder'] = "../plots/3C405_x136060/" +# globals()['data_folder'] = "../data/3C405_x136060/" +# infiles = ['x1360601t_c0f.fits','x1360602t_c0f.fits','x1360603t_c0f.fits'] +# globals()['plots_folder'] = "../plots/3C405_x136060/" # globals()['data_folder'] = "../data/CygnusA_x43w0/" # infiles = ['x43w0101r_c0f.fits', 'x43w0102r_c0f.fits', 'x43w0103r_c0f.fits', @@ -91,22 +91,22 @@ def main(): # Data binning rebin = True if rebin: - pxsize = 0.50 + pxsize = 0.10 px_scale = 'arcsec' #pixel or arcsec rebin_operation = 'sum' #sum or average # Alignement align_center = 'image' #If None will align image to image center display_data = False # Smoothing - smoothing_function = 'combine' #gaussian_after, gaussian or combine - smoothing_FWHM = 1.00 #If None, no smoothing is done + smoothing_function = 'gaussian' #gaussian_after, gaussian or combine + smoothing_FWHM = 0.10 #If None, no smoothing is done smoothing_scale = 'arcsec' #pixel or arcsec # Rotation rotate_stokes = True #rotation to North convention can give erroneous results rotate_data = False #rotation to North convention can give erroneous results # Polarization map output - figname = '3C405_FOC' #target/intrument name - figtype = '_combine_FWHM100_rot' #additionnal informations + figname = 'NGC1068_FOC' #target/intrument name + figtype = '_gaussian_FWHM010_rot' #additionnal informations SNRp_cut = 3 #P measurments with SNR>3 SNRi_cut = 30 #I measurments with SNR>30, which implies an uncertainty in P of 4.7%. step_vec = 1 #plot all vectors in the array. if step_vec = 2, then every other vector will be plotted