[bug #66588] Arithmetic's trim operator gives 2x2 image when all image is blank

[Mohammad Akhlaghi]

URL:
  <https://savannah.gnu.org/bugs/?66588>

                 Summary: Arithmetic's trim operator gives 2x2 image when all
image is blank
                   Group: GNU Astronomy Utilities
               Submitter: makhlaghi
               Submitted: Thu 19 Dec 2024 09:10:02 PM CET
                Category: Arithmetic
                Severity: 3 - Normal
              Item Group: Output not reasonable
                  Status: In Progress
                 Privacy: Public
             Assigned to: makhlaghi
             Open/Closed: Open
         Discussion Lock: Any


    _______________________________________________________

Follow-up Comments:


-------------------------------------------------------
Date: Thu 19 Dec 2024 09:10:02 PM CET By: Mohammad Akhlaghi <makhlaghi>
As of Gnuastro 0.23, when the input to the 'trim' operator of Arithmetic is
all blank, it will produce an output that has 2x2=4 pixels (all still blank!)!
But this is not correct: it should be 0x0 (there are no pixels with a
non-blank value!).

To reproduce the problem, run the following commands:


$ astarithmetic 10 10 2 makenew nan + --output input.fits
$ astarithmetic input.fits trim -o trimmed.fits
$ astfits trimmed.fits --keyvalue=NAXIS1,NAXIS2 --quiet
2      2









    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?66588>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/

signature.asc
Description: PGP signature