astronomy-laboratory / pyoof

pyoof is a Python package which performs out-of-focus (OOF) holography on astronomical beam maps for single-dish radio telescopes.
BSD 3-Clause "New" or "Revised" License
0 stars 1 forks source link

Create reader Artemis data #1

Open tcassanelli opened 1 year ago

tcassanelli commented 1 year ago

A new and smarter function to read in .fits files or even HDF5 (.h5) should be implemented. Notice that the number of observed frequencies and number of beam maps should also be expanded here.

https://github.com/AstroLab-UChile/pyoof/blob/94d1e324837ededf2b1886ed1ebdfcebd2fa7474/pyoof/aux_functions.py#L22-L97

tcassanelli commented 1 year ago

It may be important to talk to Rodrigo and try to have a common way to store the data. Also to check here would be whether we need another extra step of data reduction prior storing the files? i.e., will there be a manual basket weaving? I'd rather to keep this extra step out of pyoof. Breaking this down in #5 .

VicentePenaLet commented 1 year ago

I was able to perform holography with the SEPIA data that Rodrigo sent. I had issues at first trying to get the beam maps for the different $d_z$ to have the same shape. After some trial and error I managed to get it working properly.

Looking at the computed beam in comparison with the real beam we get a general outline of true shape but the detail is not captured well. I'm trying different parameters to see if the results change a bit. below is an example of the real beam, computed beam, and the phase map.

sepia

fitbeam_n5_0 000658 fitphase_n5_0 000658

tcassanelli commented 1 year ago

I'm surprised that fit do not look that bad, I thought that weirdly shaped beam will affect more the output. Interesting ... Since the beams look sort of offset from the center, would it be a good idea perhaps to include the illumination position as a parameter?

tcassanelli commented 1 year ago

Can you also show the residual of the beam maps? by doing so we can check whether the pulse is off-center or note.