add data rotation (instead of stokes rotation) and add sentinels

This commit is contained in:
Thibault Barnouin
2021-06-10 20:02:53 +02:00
parent f7b6d9dfaf
commit 1dc40f96c4
41 changed files with 178 additions and 43 deletions

View File

@@ -134,7 +134,7 @@ def polarization_map(Stokes, SNRp_cut=3., SNRi_cut=30., step_vec=1,
pang_err = Stokes[np.argmax([Stokes[i].header['datatype']=='Pol_ang_err' for i in range(len(Stokes))])]
pivot_wav = Stokes[0].header['photplam']
convert_flux = Stokes[0].header['photflam']
convert_flux = 1.#Stokes[0].header['photflam']
wcs = WCS(Stokes[0]).deepcopy()
#Compute SNR and apply cuts