search

JMMC-OpenDev / oimaging

This module contains the OImaging java application developed by the JMMC technical team.
GNU General Public License v3.0
4 stars 3 forks source link

Beam #24

Open FerreolS opened 2 years ago

FerreolS commented 2 years ago

From the oifits files it is possible to estimate the beam that is: minor, major axes and PA of an ellipse describing the resolution of the interferometer.

FerreolS commented 2 years ago

There is 2 ways to estimate the beam

FerreolS commented 2 years ago

I propose to compute the beam by estimating the covariance matrix of the (u,v) coordinate of the baseline. Its eigen values will be the axis and the eigenvector will gives the angle of the ellipse in Fourier domain. Then convert to image domain.

bourgesl commented 4 months ago

To compute the beam, should we use U-V frequencies from VIS2 tables only or also consider VIS and T3 tables ? Some OIFITS files may not have VIS2 tables, then no beam estimation. As there is redundancy between tables, how to normalize properly ?

FerreolS commented 4 months ago

Good question All the UV frequencies has to be taken only once. This IMO drives out all T3 tables has I believe that if 3 telescopes are cophased then there must be corresponding VIS2 or VIS for the 3 baselines. (even if, in gravity, the recorded uv of the T3 differs do to integration difference). For VIS and VIS2, it depends of the instrumental mode. If there is phase reference, we should take the UV of the VIS mainly (no need of V2 in that case), if the VIS contains only differencial quantities, the UV coordinates are probably counted twice in VIS and VIS2 table. As you can have files with mix of these modes, depending of how the tables were built, we probably need to retrieve all the non-redondant UV coordinate for both table.

FerreolS commented 4 months ago

I propose to compute the beam by estimating the covariance matrix of the (u,v) coordinate of the baseline. Its eigen values will be the axis and the eigenvector will gives the angle of the ellipse in Fourier domain. Then convert to image domain.

To keep track of the information the Julia code, the documentation and a demo for that is here

bourgesl commented 4 months ago

Here is OImaging beta release that computes the beam on VIS2 only and displays the ellipse at the image center (overlay): https://www.jmmc.fr/apps/beta/OImaging/

bourgesl commented 4 months ago

What algorithm can be used to discard redundant UV coordinate ? it seems tricky to take precision into account (gridding ?)

FerreolS commented 4 months ago

Nice! Should we add a button to display/hide it and maybe move it to a corner.

FerreolS commented 4 months ago

What algorithm can be used to discard redundant UV coordinate ? it seems tricky to take precision into account (gridding ?)

I don't know. Probably the redondant UV coordinates in VIS and VIS2 have exactly the same value. This has to be checked in practice

FerreolS commented 4 months ago

Here is OImaging beta release that computes the beam on VIS2 only and displays the ellipse at the image center (overlay): https://www.jmmc.fr/apps/beta/OImaging

✔️ I have checked that the plotted ellipse in the exactly the same than in my notebook for RCar

bourgesl commented 3 months ago

I updated OImaging beta: https://www.jmmc.fr/apps/beta/OImaging/

Changes:

I left the beam at the image center as the plot is zoomable

jean-philippeberger commented 3 months ago

👍 Ca serait pas mal d’avoir des retours de validation (je vais tester ça dès que je me libère). Je sais que Ferréol a déjà testé la fonctionalité.

JPB

Jean-Philippe Berger Institut de Planétologie et d’Astrophysique de Grenoble http://ipag.osug.fr/ Univ. Grenoble Alpes, CNRS Tél: +33 (0)4 76 74 86 68
414 Rue de la Piscine 38400 Saint Martin d’Hères
France Email: @.***

On 4 Jun 2024, at 10:11, Laurent Bourgès @.***> wrote:

I updated OImaging beta: https://www.jmmc.fr/apps/beta/OImaging/

Changes:

Added the toogle button 'Beam' on image panels to show or hide the beam info (displayed in its tooltip) Modified the beam estimation to use distinct spatial frequencies from both OI_VIS and OI_VIS2 tables (with small rounding of spatial frequencies) I left the beam at the image center as the plot is zoomable

— Reply to this email directly, view it on GitHub https://github.com/JMMC-OpenDev/oimaging/issues/24#issuecomment-2146885186, or unsubscribe https://github.com/notifications/unsubscribe-auth/AKM4IURUYGFGBZT4AAYMBLTZFVZDDAVCNFSM5FT5EVS2U5DIOJSWCZC7NNSXTN2JONZXKZKDN5WW2ZLOOQ5TEMJUGY4DQNJRHA3A. You are receiving this because you are subscribed to this thread.