GNU 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: |
Mohammad Akhlaghi <makhlaghi> |
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> |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
Carbon-Copy List
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
Powered by Savane 3.13-f8d8.
Corresponding source code
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)!