#
#  Original FITS-WCS header:
#
#NAXIS   =                    2 / number of data axes
#NAXIS1  =                  300 / length of data axis 1
#NAXIS2  =                  300 / length of data axis 2
#CRPIX1  =                205.5 / Reference pixel on axis 1
#CRPIX2  =                205.5 / Reference pixel on axis 2
#CRVAL1  =              10.6847 / Value at ref. pixel on axis 1
#CRVAL2  =               41.269 / Value at ref. pixel on axis 2
#CTYPE1  = 'RA---TAN'           / Type of co-ordinate on axis 1
#CTYPE2  = 'DEC--TAN'           / Type of co-ordinate on axis 2
#CDELT1  =                -0.01 / Pixel size on axis 1
#CDELT2  =                 0.01 / Pixel size on axis 2
#PC1_1   =    0.866025405812433 / Transformation matrix element
#PC1_2   =                 -0.5 / Transformation matrix element
#PC2_1   =                  0.5 / Transformation matrix element
#PC2_2   =    0.866025405812433 / Transformation matrix element
#MJD-OBS =     51544.4992571308 / Modified Julian Date of observation
#DATE-OBS= '2000-01-01T11:58:55.816' / Date of observation
#RADESYS = 'FK5     '           / Reference frame for RA/DEC values
#EQUINOX =               2000.0 / [yr] Epoch of reference equinox
#
#
# Transform Set
#  Holds Pixel to World transform spec
  mappings: !<trans:TransformSet>
    node: &_ThA5vHlK4jh3IGS !<trans:TransNode>
      target: *_S0M0RnlbfLIhr2c
    node: !<trans:TransNode>
      parent: *_ThA5vHlK4jh3IGS
      target: *_SXfLKmrjxFP5Vkg
      mapping: *_TSC6b1zgw82amDv

# Pixel Coordinate System
  pixelCoordSys: &_S0M0RnlbfLIhr2c !<coords:PixelCoordSystem>
    pixelSpace !<coords:PixelSpace>
      axis: !<coords:BinnedAxis>
        name: 'Pixel axis 1'
        length: 300
      axis: !<coords:BinnedAxis>
        name: 'Pixel axis 2'
        length: 300

#  Spatial Coordinate System
#   - with Standard Celestial Coordinate space (default)
  spaceCoordSys: &_SXfLKmrjxFP5Vkg !<coords:SpaceSys>
    frame: !<coords:SpaceFrame">
      refPosition: !<coords:StdRefLocation>
        position: 'BARYCENTER'
      spaceRefFrame: 'FK5'
      equinox: 'J2000'

#  The pixel->sky mapping: compose the usual FITS-WCS steps in series.
  pix2world: &_TSC6b1zgw82amDv !<trans:ComposeMap>
    cptMap: !<trans:MapRef>
      ref: *_TnaGlYqhbMMiw8w
    cptMap: !<trans:MapRef>
      ref: *_T0VHFpYY4ZVX5DO
    cptMap: !<trans:MapRef>
      ref: *_T6RF4HGJb8OBrFo
    cptMap: !<trans:MapRef>
      ref: *_TEtpOT0sLvn35tr

#  Shift the pixel origin to the reference point
  shift2D: &_TnaGlYqhbMMiw8w !<trans:ConcatenateMap>
    cptMap: !<trans:MapRef>
      ref: *_TWFdjcHxLauGm4x
    cptMap: !<trans:MapRef>
      ref: *_TEEKYDdiRG5s3LY

  shiftX: &_TWFdjcHxLauGm4x !<trans:Shift>
    offset: -205.5
  shiftY: &_TEEKYDdiRG5s3LY !<trans:Shift>
    offset: -205.5

#  Rotate and scale the shifted pixel axes to get Intermediate World
#  Coordinates.
  scaleRotate: &_T0VHFpYY4ZVX5DO !<trans:Matrix>
    M: 2
    N: 2
    element: !<trans:MatrixElement>
      m: 0
      n: 0
      value: -0.00015114994737346996
    element: !<trans:MatrixElement>
      m: 0
      n: 1
      value: 8.7266462599716482e-05
    element: !<trans:MatrixElement>
      m: 1
      n: 0
      value: 8.7266462599716482e-05
    element: !<trans:MatrixElement>
      m: 1
      n: 1
      value: 0.00015114994737346996

#  De-projection from intermediate world coordinates (the projection
#  plane) to the native spherical coordinates. Note, the document doesn't
#  make it clear whether the forward direction of a "SkyProjection" goes
#  from cartesian to spherical, or the other way round. I've assumed the
#  former here.
  skyProjection: &_T6RF4HGJb8OBrFo !<trans:SkyProjection>
    algorithm: 'TAN'

#  Rotate native spherical coordinates to (RA,Dec). AST uses a 3x3 matrix
#  to represent the required rotations rather than three Euler angles.
#  The following angles have been derived manually from the 3x3 Matrix
#  generated by AST.  We'll probably need to iterate to agree on all the
#  various conventions assumed by this conversion.
  eulerRot: &_TEtpOT0sLvn35tr !<trans:EulerRotation>
    eulerAngle: !<trans:EulerAngle>
      rotationAxis: 'x'
      angle: 0.0
    eulerAngle: !<trans:EulerAngle>
      rotationAxis: 'y'
      angle: -48.731
    eulerAngle: !<trans:EulerAngle>
      rotationAxis: 'z'
      angle: -10.6847