next up previous contents index
Next: RANDOM_MAP Up: GREG2 Language Internal Help Previous: POLYGON

   
PROJECTION

        [GREG2\]PROJECTION [A0 D0 [Angle]] [/TYPE Ptype]

    Define a projection of the (celestial) sphere  from  point  of  (A0,D0),
    (which  are  Longitude and Latitude respectively) of the specified type.
    Angle is the angle between the Y axis and the North pole.  The  previous
    values are kept if no argument is specified.  All angles are in degrees,
    except if the SYSTEM is  EQUATORIAL  in  which  case  A0  is  the  right
    ascension  and must be specified in hours. Formats like  -dd:mm:ss.s  or
      hh:mm.mmm  in sexagesimal notation up to the point field are  allowed.
    After the point, decimal values are assumed.

    When a projection is active, the User  coordinates  are  assumed  to  be
    projected  coordinates  of  the  sphere,  and hence in the case of small
    field of view where distortion are  negligible,  correspond  to  angular
    offsets  MEASURED  IN  RADIANS.  The  field of view of the projection is
    defined by command LIMITS. For convenience (?), LIMITS may be  specified
    in  a  so-called  ANGLE_UNIT  which  can  be  set to Seconds, Minutes or
    Degrees. Note that only the numbers you give are assumed to be  in  this
    "ANGLE_UNIT" :  they  are immediately converted to the natural projected
    unit (equivalent to Radian in case of small field) by  division  by  the
    appropriate  scaling  factor.   That  is  all what this improperly named
    ANGLE_UNIT means.

    The TYPE can be
      - NONE            Disables the  projection  system.  User  coordinates
        then  loose  their interpretation in terms of projected coordinates.
        The ANGLE_UNIT is then totally ignored.
      - GNOMONIC        Radial projection on the tangent plane. Being R  and
        P the (angular) polar coordinates from the projection point (tangent
        point), the projected coordinates are given by X = Tan(R).Sin(P) and
        Y = Tan(R).Cos(P) .
      - ORTHOGRAPHIC   View  from  infinity.  X  =  Sin(R).Sin(P)  and  Y  =
        Sin(R).Cos(P)
      - AZIMUTHAL      Spherical offsets from the projection  center.   X  =
        R.Sin(P) and Y = R.Cos(P).
      - STEREOGRAPHIC  Uses Tan(R/2) instead of Tan(R),  and  is  thus  less
        distorted  than  the  Gnomonic projection. This is an inversion from
        the opposite pole.
      - LAMBERT        Equal  area   projection.   Projected   distance   is
        2*Sin(R)/Sqrt(2*(1+Cos(R)).
      - AITOFF         Equal area projection. Angle and D0 are ignored
      - RADIO          The standard  radio  astronomy  single  dish  mapping
        "projection", in which X = (A-A0).COS(D0) and Y = D-D0. The Angle is
        obviously ignored.


Gildas manager
1999-03-12