gnuastro-devel
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[task #15676] Calculating SIP to TPV coversions


From: Sachin Kumar Singh
Subject: [task #15676] Calculating SIP to TPV coversions
Date: Sun, 7 Jun 2020 18:25:44 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:76.0) Gecko/20100101 Firefox/76.0

Update of task #15676 (project gnuastro):

                  Status:             In Progress => Done                   
        Percent Complete:                      0% => 100%                   
             Open/Closed:                    Open => Closed                 

    _______________________________________________________

Follow-up Comment #1:

The task is complete and now SIP->PV:-). The output of the sample file used to
create the SIP fits file, which was used to make the final PV file are shown
below:


>From test Gnuastro

XTENSION= 'IMAGE   '           / IMAGE extension                              
 
BITPIX  =                  -32 / number of bits per data pixel                
 
NAXIS   =                    2 / number of data axes                          
 
NAXIS1  =                 3352 / length of data axis 1                        
 
NAXIS2  =                 2532 / length of data axis 2                        
 
PCOUNT  =                    0 / required keyword; must = 0                   
 
GCOUNT  =                    1 / required keyword; must = 1                   
 
                                                                              
 
                      / World Coordinate System (WCS)                         
 
WCSAXES =                    2 / Number of coordinate axes                    
 
CRPIX1  =        1763.09429932 / Pixel coordinate of reference point          
 
CRPIX2  =         612.79893494 / Pixel coordinate of reference point          
 
PC1_1   =   -0.015946182574045 / Coordinate transformation matrix element     
 
PC1_2   =    0.999917332233850 / Coordinate transformation matrix element     
 
PC2_1   =   -0.999829574614539 / Coordinate transformation matrix element     
 
PC2_2   =   -0.015870824365013 / Coordinate transformation matrix element     
 
CDELT1  =  0.00078807879448870 / [deg] Coordinate increment at reference point
 
CDELT2  =  0.00078813807299668 / [deg] Coordinate increment at reference point
 
CUNIT1  = 'deg'                / Units of coordinate increment and value      
 
CUNIT2  = 'deg'                / Units of coordinate increment and value      
 
CTYPE1  = 'RA---TPV'           / TAN (gnomonic) projection + distortions      
 
CTYPE2  = 'DEC--TPV'           / TAN (gnomonic) projection + distortions      
 
CRVAL1  =       113.3064285324 / [deg] Coordinate value at reference point    
 
CRVAL2  =        31.8788976712 / [deg] Coordinate value at reference point    
 
LONPOLE =                180.0 / [deg] Native longitude of celestial pole     
 
LATPOLE =        31.8788976712 / [deg] Native latitude of celestial pole      
 
MJDREFI =                  0.0 / [d] MJD of fiducial time, integer part       
 
MJDREFF =                  0.0 / [d] MJD of fiducial time, fractional part    
 
RADESYS = 'FK5'                / Equatorial coordinate system                 
 
EQUINOX =               2000.0 / [yr] Equinox of equatorial coordinates       
 
PV1_1   =   9.999868977236E-01 / TPV coefficient: x                           
 
PV1_2   =  -8.842015406258E-06 / TPV coefficient: y                           
 
PV1_4   =   1.822783853090E-04 / TPV coefficient: xx                          
 
PV1_5   =   1.365776728868E-04 / TPV coefficient: xy                          
 
PV1_6   =   2.257022520460E-05 / TPV coefficient: yy                          
 
PV1_7   =   1.920222380358E-04 / TPV coefficient: xxx                         
 
PV1_8   =  -3.830056501464E-05 / TPV coefficient: xxy                         
 
PV1_9   =  -1.347022772508E-04 / TPV coefficient: xyy                         
 
PV1_10  =  -2.471225520868E-05 / TPV coefficient: yyy                         
 
PV2_0   =   4.672062523823E-04 / TPV coefficient: 1                           
 
PV2_1   =   1.000878428319E+00 / TPV coefficient: y                           
 
PV2_2   =  -4.973468945505E-04 / TPV coefficient: x                           
 
PV2_4   =  -1.411294363996E-04 / TPV coefficient: yy                          
 
PV2_5   =   7.439750263035E-04 / TPV coefficient: yx                          
 
PV2_6   =  -1.224454154415E-03 / TPV coefficient: xx                          
 
PV2_7   =  -7.014289760876E-04 / TPV coefficient: yyy                         
 
PV2_8   =   3.351859178423E-04 / TPV coefficient: yyx                         
 
PV2_9   =  -4.470842932691E-04 / TPV coefficient: yxx                         
 
PV2_10  =   9.945076014401E-04 / TPV coefficient: xxx                         
 
                                                                              
 
                      / Versions and date                                     
 
DATE    = '2020-06-07T21:30:48' / file creation date (YYYY-MM-DDThh:mm:ss UT) 
 
CFITSIO = '3.47    '           / CFITSIO version.                             
 
WCSLIB  = '7.1     '           / WCSLIB version.                              
 
GSL     = '2.6     '           / GNU Scientific Library version.              
 
GNUASTRO= '0.12    '           / GNU Astronomy Utilities version.             
 
COMMIT  = 'e00f02a-dirty'      / Git's commit description in running dir.    




Sample file used for testing:

XTENSION= 'IMAGE   '           / IMAGE extension                              
 
BITPIX  =                  -32 / number of bits per data pixel                
 
NAXIS   =                    2 / number of data axes                          
 
NAXIS1  =                 3352 / length of data axis 1                        
 
NAXIS2  =                 2532 / length of data axis 2                        
 
PCOUNT  =                    0 / required keyword; must = 0                   
 
GCOUNT  =                    1 / required keyword; must = 1                   
 
                                                                              
 
                      / Versions and date                                     
 
DATE    = '2020-05-13T21:59:29' / file creation date (YYYY-MM-DDThh:mm:ss UT) 
 
CFITSIO = '3.47    '           / CFITSIO version.                             
 
WCSLIB  = '6.4     '           / WCSLIB version.                              
 
GSL     = '2.6     '           / GNU Scientific Library version.              
 
GNUASTRO= '0.11.57-561c'       / GNU Astronomy Utilities version.             
 
HISTORY   Astrometric solution by SCAMP version 2.6.7 (2020-04-09)            
 
COMMENT   (c) 2010-2016 IAP/CNRS/UPMC                                         
 
COMMENT                                                                       
 
EQUINOX =        2000.00000000 / Mean equinox                                 
 
RADESYS = 'ICRS    '           / Astrometric system                           
 
CTYPE1  = 'RA---TPV'           / WCS projection type for this axis            
 
CTYPE2  = 'DEC--TPV'           / WCS projection type for this axis            
 
CUNIT1  = 'deg     '           / Axis unit                                    
 
CUNIT2  = 'deg     '           / Axis unit                                    
 
CRVAL1  =   1.133064285324E+02 / World coordinate on this axis                
 
CRVAL2  =   3.187889767116E+01 / World coordinate on this axis                
 
CRPIX1  =   1.763094299320E+03 / Reference pixel on this axis                 
 
CRPIX2  =   6.127989349400E+02 / Reference pixel on this axis                 
 
CD1_1   =  -1.256684833965E-05 / Linear projection matrix                     
 
CD1_2   =   7.880136457752E-04 / Linear projection matrix                     
 
CD2_1   =  -7.880037542618E-04 / Linear projection matrix                     
 
CD2_2   =  -1.250840093191E-05 / Linear projection matrix                     
 
PV1_0   =   9.714784609182E-05 / Projection distortion parameter              
 
PV1_1   =   9.999868977236E-01 / Projection distortion parameter              
 
PV1_2   =  -8.842015406255E-06 / Projection distortion parameter              
 
PV1_4   =   1.822783853090E-04 / Projection distortion parameter              
 
PV1_5   =   1.365776728868E-04 / Projection distortion parameter              
 
PV1_6   =   2.257022520459E-05 / Projection distortion parameter              
 
PV1_7   =   1.920222380359E-04 / Projection distortion parameter              
 
PV1_8   =  -3.830056501464E-05 / Projection distortion parameter              
 
PV1_9   =  -1.347022772508E-04 / Projection distortion parameter              
 
PV1_10  =  -2.471225520868E-05 / Projection distortion parameter              
 
PV2_0   =   4.672062523823E-04 / Projection distortion parameter              
 
PV2_1   =   1.000878428319E+00 / Projection distortion parameter              
 
PV2_2   =  -4.973468945505E-04 / Projection distortion parameter              
 
PV2_4   =  -1.411294363997E-04 / Projection distortion parameter              
 
PV2_5   =   7.439750263035E-04 / Projection distortion parameter              
 
PV2_6   =  -1.224454154415E-03 / Projection distortion parameter              
 
PV2_7   =  -7.014289760875E-04 / Projection distortion parameter              
 
PV2_8   =   3.351859178423E-04 / Projection distortion parameter              
 
PV2_9   =  -4.470842932692E-04 / Projection distortion parameter              
 
PV2_10  =   9.945076014402E-04 / Projection distortion parameter              
 
FGROUPNO=                    1 / SCAMP field group label                      
 
ASTIRMS1=   0.000000000000E+00 / Astrom. dispersion RMS (intern., high S/N)   
 
ASTIRMS2=   0.000000000000E+00 / Astrom. dispersion RMS (intern., high S/N)   
 
ASTRRMS1=   2.056759573495E-04 / Astrom. dispersion RMS (ref., high S/N)      
 
ASTRRMS2=   2.493773027890E-04 / Astrom. dispersion RMS (ref., high S/N)      
 
ASTINST =                    1 / SCAMP astrometric instrument label           
 
FLXSCALE=   0.000000000000E+00 / SCAMP relative flux scale                    
 
MAGZEROP=           0.00000000 / SCAMP zero-point                             
 
PHOTIRMS=           0.00000000 / mag dispersion RMS (internal, high S/N)      
 
PHOTINST=                    1 / SCAMP photometric instrument label           
 
PHOTLINK=                    F / True if linked to a photometric field        
 
END



The code for the same is pushed here
<https://gitlab.com/sachinkumarsingh092/gnuastro-test-files>. This task can be
closed now.

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/task/?15676>

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




reply via email to

[Prev in Thread] Current Thread [Next in Thread]