taskGNU Astronomy Utilities - Tasks: task #15801, Point-based interpolation method...


You are not allowed to post comments on this tracker with your current authentication level.

task #15801: Point-based interpolation method for rebinning

Submitter:  Raul Infante-Sainz <infantesainz>
Submitted:  Mon 26 Oct 2020 10:20:36 AM UTC
Should Start On:  Mon 26 Oct 2020 12:00:00 AM UTC Should be Finished on:  Mon 26 Oct 2020 12:00:00 AM UTC
Category:  Warp Priority:  5 - Normal
Item Group:  New feature Status:  None
Privacy:  Public Assigned to:  None
Percent Complete:  0% Open/Closed:  Open
Effort:  0.00

Mon 26 Oct 2020 11:04:20 AM UTC, comment #1: 

This is part of what we discussed privately but can be relevant:

Most existing interpolation methods (for example those in Scipy, that Raul initially referenced) are interpolation-based: assuming each pixel as a point in its center on the coordinate grid. Then finding the new grid values based on the positions of the centers of the pixels of the new grid.

But so far, Warp doesn't use this type of interpolation (even though it is so much easier to implement!) because it is terrible for low signal-to-noise images and generally interpolation always adds so many systematics and adds noise.

What Warp currently does is to put the new pixel over the old one, find the polygon within each input pixel that overlaps with each output pixel. calculate that polygon's area, and multiply the fractional area by the pixel value (to get the fractional value) and finally to put that fractional value into the value of the output pixel. There is thus no interpolation involved and everything is robust and accurate and the signal-to-noise increases, its like a small convolution.

But I am not against using a point-based interpolation method of rebinning if someone else is interested to add it (I just don't have time to implement it myself)!

Mohammad Akhlaghi <makhlaghi>
Group administrator
Mon 26 Oct 2020 10:20:36 AM UTC, original submission:  

I ask to Mohammad about having a way of rebinning images in which we could specify the operator (mean, median, etc.). After a small discussion we ended wit the conclusion that having this kind of rebinning method (point-based/interpolation) in Gnuastro would be good.

Raul Infante-Sainz <infantesainz>
Group Member


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by makhlaghi (Posted a comment)
  • -email is unavailable- added by infantesainz (Submitted the item)

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.


    No changes have been made to this item

    Back to the top

    Powered by Savane 3.12.
    Corresponding source code