Fix interpolation in central pixel
mtazzari opened this issue · comments
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.
Working on this.
Add a warning if dR>dxy/5. It is recommended to have at least 10 points in the central pixel