search

mtazzari / galario

Gpu Accelerated Library for Analysing Radio Interferometer Observations
https://mtazzari.github.io/galario/
GNU Lesser General Public License v3.0
31 stars 15 forks source link

Fix interpolation in central pixel #129

Closed mtazzari closed 6 years ago

mtazzari commented 6 years ago

Interpolating to R=0 in the central pixel leads to erroneous high values if the ints(R) array is very steep for R tending to 0. We should instead use flux conservation inside the central pixel, as I did in pyvfit: compute the total brightness inside the cell. It requires a simple integral (e.g. trapezoidal) on the points falling inside the pixel and dividing by the area.

mtazzari commented 6 years ago

Working on this. Add a warning if dR>dxy/5. It is recommended to have at least 10 points in the central pixel