; REGRD
;---------------------------------------------------------------
;! Regrids an image from one co-ordinate frame to another
;# TASK IMAGING COORDINATES
;-----------------------------------------------------------------------
;; Copyright (C) 1995, 2002, 2009, 2013
;; Associated Universities, Inc. Washington DC, USA.
;;
;; This program is free software; you can redistribute it and/or
;; modify it under the terms of the GNU General Public License as
;; published by the Free Software Foundation; either version 2 of
;; the License, or (at your option) any later version.
;;
;; This program is distributed in the hope that it will be useful,
;; but WITHOUT ANY WARRANTY; without even the implied warranty of
;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
;; GNU General Public License for more details.
;;
;; You should have received a copy of the GNU General Public
;; License along with this program; if not, write to the Free
;; Software Foundation, Inc., 675 Massachusetts Ave, Cambridge,
;; MA 02139, USA.
;;
;; Correspondence concerning AIPS should be addressed as follows:
;; Internet email: aipsmail@nrao.edu.
;; Postal address: AIPS Project Office
;; National Radio Astronomy Observatory
;; 520 Edgemont Road
;; Charlottesville, VA 22903-2475 USA
;-----------------------------------------------------------------------
REGRD LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC
REGRD: Transforms image coordinate system and geometry
INNAME Input name
INCLASS Input class
INSEQ 0.0 9999.0 Input sequence, 0 -> high
INDISK 0.0 9.0 Input disk, 0 -> any
OUTNAME Output name
OUTCLASS Output class
OUTSEQ -1.0 9999.0 Output sequence
0 -> highest unique
OUTDISK 0.0 9.0 Output disk
0 -> highest with room
IMSIZE 0.0 16384.0 Output image size
APARM Input map parameters, (if not
specified in the map header)
1) epoch of mean coordinates
2) epoch prefix
1: Julian (eg J2000.0)
2: Besselian (eg B1950.0)
3: Besselian without
E-terms (eg b1950.0)
BPARM Output map parameters
1) coordinate system
1: equatorial
2: galactic
3: ecliptic
2) epoch of mean coordinates
3) epoch prefix
1: Julian (eg J2000.0)
2: Besselian (eg B1950.0)
3: Besselian without
E-terms (eg b1950.0)
4) geometry
1: SIN 7: AIT
2: TAN 8: STG
3: ARC 9: CAR
4: NCP 10: MOL
5: GLS 11: PAR
6: MER
5) blanking control
0: "magic blanking"
1: zeros
CPARM Output axis specification
1-5): first axis (see HELP)
1: hour (or degree)
2: minute (or arcmin)
3: second (or arcsec)
4: reference pixel
5: coord increment (arcsec)
6-10): second axis similarly
See HELP for important information concerning
transforming coordinates, and usage of CPARM.
----------------------------------------------------------------
REGRD
Type: Task
Use: REGRD will regrid an image from one coordinate frame
and geometry to another. This includes precession as
described below.
COORDINATE SYSTEMS.
Coordinate transformations between the IAU1976 and
Bessel-Newcomb systems are done with full precision
assuming zero proper motion, parallax, and recessional
velocity at J2000.0
Specifying a Julian epoch 'J' to REGRD implies that the
coordinates are referenced to the new IAU1976/FK5 system.
Specifying a Besselian epoch to REGRD implies that the
coordinates are referenced to the old Bessel-Newcomb/FK4
system. An epoch prefix of 'B' indicates the convention
that the coordinates include the effect of the E-terms,
whereas 'b' indicates that they have already been removed.
FK4 catalogue cooordinates were not corrected for the
elliptic terms of aberration (E-terms) except for
positions within 10 degrees of the pole. Most earlier
catalogues did not correct for them.
The default behaviour here is to assume that the E-terms
are included in all Besselian coordinates (including near
the pole). This can be defeated if it is known that the
input coordinates have already been corrected, or if it is
required that the output coordinates not contain them.
See the EXPLAIN section for a brief description of the
algorithm.
Adverbs:
INNAME......Input image name, standard defaults.
INCLASS.....Input image class, standard defaults.
INSEQ.......Input image sequence number, 0 -> highest.
INDISK......Input disk drive number, 0 -> any.
OUTNAME.....Output image name, standard defaults.
OUTCLASS....Output image class, standard defaults.
OUTSEQ......Output image sequence number, 0 -> highest unique
OUTDISK.....Output disk drive number, 0 -> highest with
space.
IMSIZE......Output image size (pixels), maximum 16384.
APARM.......Input map parameters to be used if not defined in
the image header.
1) epoch of mean equatorial or ecliptic
coordinates, e.g. 1950, 2000.
2) epoch prefix,
1 = Julian (as J2000.0)
2 = Besselian (as B1950.0)
3 = Besselian without E-terms (eg b1950.0)
If not defined in the header or APARMs, the epoch
defaults to 2000.0. The prefix defaults to B if
the epoch is 1950.0, otherwise to J,
BPARM.......Coordinate frame and geometry of the output map
1) coordinate frame
1: equatorial (mean of epoch)
2: galactic
3: ecliptic (mean of epoch)
anything else defaults to the input map value.
2) epoch of mean equatorial or ecliptic
coordinates, e.g. 1950, 2000.
If negative or zero, defaults to the input
value if defined, otherwise 2000.0
3) epoch prefix,
1 = Julian (as J2000.0)
2 = Besselian (as B1950.0)
3 = Besselian without E-terms (eg b1950.0)
Anything else defaults to the input value if
defined, otherwise if the epoch is 1950.0 it
defaults to B, otherwise to J.
4) geometry: projective ones
1: SIN, sine (orthographic) - NEW default
2: TAN, tangent (gnomonic)
3: ARC, arc (zenithal equidistant)
4: NCP, north celestial pole tangent
5: STG, stereographic
----- all sky types ----------------
6: GLS, global sinusoid (Sanson-Flamsteed)
7: MER, Mercator
8: AIT, Hammer-Aitov
9: CAR, Plate Carree ("cartesion")
10: MOL, Molweide's
11: PAR, Parabolic (Craster)
These should have ref latitude = 0 or they will be
"oblique" which you probably do not want.
5) output blanking control
0: "magic" blanking
1: zero
CPARM.......Output axis specification: all zero causes the input
reference pixel to be used (as transformed).
1-5) apply to the first axis
1-3) specify the coordinate reference pixel.
IF THE VALUE SPECIFIED IS OUTSIDE THE RANGE
-24HR TO +24HR, OR -360 TO +360 DEGREES, THE
(TRANSFORMED) COORDINATES OF THE CENTRE OF
THE INPUT MAP WILL BE USED.
1: hour for equatorial, degree for the others
2: minute for equatorial, arcmin for the others
3: second for equatorial, arcsec for the others
4: coordinate reference pixel. If zero, the
centre of the map is assumed.
5: coordinate increment (arcsec per pixel,
should be negative). If zero, the input
coordinate increment will be scaled
according to the input and output image
sizes.
6-10) apply to the second axis as for the first
except the range is -90 to +90 degrees
----------------------------------------------------------------
REGRD: Transforms image coordinate system and geometry
Author: Mark Calabretta
Related tasks: GEOM, HGEOM, LGEOM, PGEOM, COMB
Algorithm
~~~~~~~~~
For each pixel in the output image:
1) compute its sky coordinates
2) transform to sky coordinates on the input map -
a) remove E-terms (equatorial FK4 coordinates only)
b) perform the spherical coordinate rotation specified by
three Euler angles
c) apply E-terms (equatorial FK4 coordinates only)
3) compute the pixel coordinates on the input map
4) interpolate the pixel value from the nearest nine pixels -
a) quadratic interpolation in X for each of the three rows
b) quadratic interpolation in Y of the result
Parameters used for the transformation (step 2) are recorded
in the history file.
E-terms
~~~~~~~
The E-terms are recomputed for every pixel (steps 2a and 2c,
computation of the E-terms is not done iteratively).
1) In RA: (E1*COS(RA) + E2*SIN(RA))/COS(DEC)
2) In DEC: (E2*COS(RA) - E1*SIN(RA))*SIN(DEC) + E3*COS(DEC)
Euler angles
~~~~~~~~~~~~
For a spherical coordinate rotation from SYSTEM1 to SYSTEM2:
1) PHI0: Longitude of the ascending node in SYSTEM1.
Of the two points of intersection of the equators of SYSTEM1
and SYSTEM2, the ascending node is the one where the equator
of SYSTEM2 crosses from south to north as viewed in SYSTEM1.
2) THETA: The angle between the poles of the two systems.
Positive for a positive rotation about the ascending node.
3) PHI: Longitude of the ascending node in SYSTEM2.
Blanking
~~~~~~~~
Blank pixels are fully accounted for in the sense that one
blank pixel in the input map produces only one blank pixel in
the output map. The basic criteria is that the output pixel
will be blank if and only if the pixel (P0) on the input map
nearest the position computed at step 3 above is blank.
If P0 is not blank and any of the eight pixels surrounding
it are, then the quadratic interpolation reduces to a linear or,
if necessary, a constant interpolation or extrapolation. In the
worst case where all of the neighbouring pixels are blank, the
"interpolated" value would be the value at P0.
----------------------------------------------------------------