Ticket #671 (closed defect: fixed)

Opened 10 years ago

Last modified 7 years ago

Need to add TAN-SIP to WCS formatter

Reported by: ajc Owned by: ktl
Priority: critical Milestone:
Component: afw Keywords:
Cc: ktl,jpizagno,ajc Blocked By:
Blocking: Project: LSST
Version Number:
How to repeat:

We need to output the tan-sip parameters from the fits header (assuming that we will not update these values). The appropriate number of polynomials terms are given by A_ORDER, B_ORDER, AP_ORDER, BP_ORDER. AP and BP are for the inverse transform. For example,

A_ORDER =     3 / polynomial order, axis 1, detector to sky
A_0_2   =     -6.341484E-05 / axis 1 distortion coefficient of v^2
A_0_3   =       -1.458834E-08 / axis 1 distortion coefficient of v^3
A_1_1   =        2.204371E-04 / axis 1 distortion coefficient of u.v
A_1_2   =  1.025434E-07 / axis 1 distortion coefficient of u.v^2
A_2_0   =  -1.001876E-04 / axis 1 distortion coefficient of u^2
A_2_1   =  -2.062207E-09 / axis 1 distortion coefficient of u^2.v
A_3_0   =   8.086817E-08 / axis 1 distortion coefficient of u^3
A_DMAX  =         1.625876 / [pixel] maximum correction
B_ORDER =        3 / polynomial order, axis 2, detector to sky
B_0_2   =       2.511898E-04 / axis 2 distortion coefficient of v^2
B_0_3   =       6.642542E-08 / axis 2 distortion coefficient of v^3
B_1_1   =        -2.703651E-05 / axis 2 distortion coefficient of u.v
B_1_2   =         5.469226E-08 / axis 2 distortion coefficient of u.v^2
B_2_0   =         3.103393E-05 / axis 2 distortion coefficient of u^2
B_2_1   =         4.404894E-08 / axis 2 distortion coefficient of u^2.v
B_3_0   =         3.954436E-08 / axis 2 distortion coefficient of u^3
B_DMAX  =             1.271021 / [pixel] maximum correction
AP_ORDER=                    3 / polynomial order, axis 1, sky to detector
AP_0_1  =         7.738011E-07 / axis 1 distortion coefficient of V
AP_0_2  =         6.345607E-05 / axis 1 distortion coefficient of V^2
AP_0_3  =        -3.144725E-08 / axis 1 distortion coefficient of V^3
AP_1_0  =        -2.824433E-06 / axis 1 distortion coefficient of U
AP_1_1  =        -2.203247E-04 / axis 1 distortion coefficient of U.V
AP_1_2  =         1.779848E-08 / axis 1 distortion coefficient of U.V^2
AP_2_0  =         1.000876E-04 / axis 1 distortion coefficient of U^2
AP_2_1  =        -7.421441E-08 / axis 1 distortion coefficient of U^2.V
AP_3_0  =        -5.308015E-08 / axis 1 distortion coefficient of U^3
BP_ORDER=                    3 / polynomial order, axis 2, sky to detector
BP_0_1  =         9.438048E-07 / axis 2 distortion coefficient of V
BP_0_2  =        -2.510898E-04 / axis 2 distortion coefficient of V^2
BP_0_3  =         6.101012E-08 / axis 2 distortion coefficient of V^3
BP_1_0  =         2.048981E-06 / axis 2 distortion coefficient of U
BP_1_1  =         2.735527E-05 / axis 2 distortion coefficient of U.V
BP_1_2  =        -8.480797E-08 / axis 2 distortion coefficient of U.V^2
BP_2_0  =        -3.104517E-05 / axis 2 distortion coefficient of U^2
BP_2_1  =        -1.130328E-08 / axis 2 distortion coefficient of U^2.V
BP_3_0  =        -4.741393E-08 / axis 2 distortion coefficient of U^3

Description

We need to be able to output TAN-SIP distortion parameters when we output a fits image with the updated crpix values.

Change History

comment:1 Changed 10 years ago by ktl

  • Status changed from new to assigned
  • Owner changed from RayPlante to ktl
  • Component changed from unknown to afw
  • Milestone set to DC3a Apps Framework

comment:2 Changed 10 years ago by ktl

  • Status changed from assigned to needinfo
  • Owner changed from ktl to ajc

Code has been written and checked in as [7626]. Please generate a test case for this.

comment:3 Changed 10 years ago by ktl

Code has now been substantially modified to use the sipA/B/Ap/Bp matrices in Wcs, with the exception of the _DMAX headers. Unfortunately, these do not seem to be set anywhere.

comment:4 Changed 10 years ago by ktl

  • Status changed from needinfo to assigned
  • Owner changed from ajc to ktl

comment:5 Changed 10 years ago by ktl

  • Status changed from assigned to closed
  • Resolution set to fixed

Setting of sipA/B/Ap/Bp matrices from a PropertySet has been implemented. Ticket branch merged to trunk in [8912].

comment:6 Changed 7 years ago by robyn

  • Milestone DC3a Apps Framework deleted

Milestone DC3a Apps Framework deleted

Note: See TracTickets for help on using tickets.