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 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

## Discussion

Mon 22 Jul 2024 02:50:09 PM UTC, comment #2:

Joseph Putko suggested a nice reference for such interpolation based warpings: https://pixinsight.com/doc/docs/InterpolationAlgorithms/InterpolationAlgorithms.html

It includes all the equations to simplify the implementation any time we decide to look into it. In particular, I really like its Figure 8 that shows the Moire pattern in noise for the different algorithms (by clicking on the name of the algorithm).

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)!

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>

## Attached Files

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

No files currently attached

## Dependencies

Depends on the following items: None found

Items that depend on this one: None found

Carbon-Copy List
• -email is unavailable- added by makhlaghi
• -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.

## History

Follows 1 latest change.

Date Changed by Updated Field Previous Value => Replaced by
2024-07-22 makhlaghi Carbon-Copy- Added -email is unavailable-