bugGNU Astronomy Utilities - Bugs: bug #62607, Converting to TPV and then SIP...

 
 

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

bug #62607: Converting to TPV and then SIP crashes for some images

Submitter:  Pedram Ashofteh-Ardakani <pedram>
Submitted:  Thu 09 Jun 2022 08:50:18 AM UTC
   
 
Category:  Fits Severity:  3 - Normal
Item Group:  Crash Status:  None
Privacy:  Public Assigned to:  None
Open/Closed:  Open

Tue 14 Jun 2022 09:57:16 AM UTC, comment #1: 

== update ==

Using 'crop' before changing the WCS distortion will cause the same crash:


$ wget "https://fits.gsfc.nasa.gov/registry/tpvwcs/tpv.fits" -O tpv.fits # size 1 Mb
$ astcrop tpv.fits -h0 --center=52.5529049,-28.7293062 --widthinpix --width=10 -otpv-crop.fits
$ astfits tpv-crop.fits -h1 --wcsdistortion SIP -o tpv-crop-sip.fits
gsl: svd.c:285: ERROR: svd of MxN matrix, M<N, is not implemented
Default GSL error handler invoked.
Aborted (core dumped)


Pedram Ashofteh-Ardakani <pedram>
Group Member
Thu 09 Jun 2022 08:50:18 AM UTC, original submission:  

Problem


I've tried converting between TPV and SIP multiple times. Usually it works fine. But there are some cases that it won't work.

Why I did this


Sometimes, opening the TPV image converted from SIP alongside another frame in the ds9 makes the program malfunction. One of the frames just shoots out of the pane! I thought maybe the converted WCS has some bugs in it.

Therefore, I tried converting multiple times to see if gnuastro might have problems with the converted image too. For some images, it worked just fine no matter how many times I converted the image back and forth.

However, for the image attached below, that is not the case. I get the following output:


$ astfits img-sip-tpv.fits -h1 --wcsdistortion SIP -o img-sip-tpv-sip.fits
gsl: svd.c:285: ERROR: svd of MxN matrix, M<N, is not implemented
Default GSL error handler invoked.
Aborted (core dumped)


  • 'img-sip.fits' is the original image
  • 'img-sip-tpv.fits' is the image converted from SIP to TPV
  • 'img-sip-tpv-sip.fits' would be the image converted from SIP, to TPV, and back to SIP: which fails


What to expect


Not a core dump! A slight offset from the original SIP can be expected since it is converted back and forth using multiplication and trigonometric functions.

Steps to reproduce



$ astfits img-sip.fits --wcsdistortion TPV -o img-sip-tpv.fits
$ astfits img-sip-tpv.fits --wcsdistortion SIP -o img-sip-tpv-sip.fits


Pedram Ashofteh-Ardakani <pedram>
Group Member

 

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

Attached Files
file #53287:  img-sip.fits added by pedram (17KiB - application/fits)
file #53288:  img-sip-tpv.fits added by pedram (11KiB - application/fits)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by pedram (Submitted the item)
  • -email is unavailable- added by pedram
  • -email is unavailable- added by pedram
  •  

    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.

     

    Follow 4 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2022-06-09 pedram Attached File- Added img-sip.fits, #53287
        Attached File- Added img-sip-tpv.fits, #53288
        Carbon-Copy- Added makhlaghi
        Carbon-Copy- Added sks_15

    Back to the top

    Powered by Savane 3.13-4448.
    Corresponding source code