Map Projection (XYCONV)

xyconv, nxycnv, utm, ikconv, nikcnv, utmik | cvinit, cviken, cvenik, cvdinit, cvdiken, cvdenik

XYCONV, NXYCNV, UTM, IKCONV, NIKCNV, UTMIK
(UTM/JTM Map Projection)
NAME
    xyconv - Lat.Lon. to Japanese Transverse Mercator (Old Tokyo Datum) X,Y
    nxycnv - Lat.Lon. to Japanese Transverse Mercator (New JGD2000 Datum) X,Y
    utm    - Lat.Lon. to Universal Transverse Mercator X,Y
    ikconv - Japanese Transverse Mercator (Old Tokyo Datum) X,Y to Lat.Lon.
    nikcnv - Japanese Transverse Mercator (New JGD2000 Datum) X,Y to Lat.Lon.
    utmik  - Universal Transverse Mercator X,Y to Lat.Lon.
                    [i: Ido = Latitude ;  k: Keido = Longitude]
 (Note) For UTM, X=0. on the equator, Y=500.(km) on the central meridian.

SYNOPSIS

      call xyconv(gi, gk, fi, fk, x, y)
      call nxycnv(gi, gk, fi, fk, x, y)
      call utm(nc, fi, fk, x, y)
  gi, gk  [float] Lat. and Lon. (in minutes) of origin
  nc      [int]   Zone number of UTM coordinaets (1-60) for Old Tokyo Datum, 
                  or plus 200 (201-260) for New JGD2000 Datum
  fi, fk  [float] Lat. and Lon. (in minutes) of the point to be calculated
  x,  y   [float] Calculated results. Northword(X) and Eastword(Y) coordinates
                  (in km) values are given.

      call ikconv(gi, gk, x, y, fi, fk)
      call nikcnv(gi, gk, x, y, fi, fk)
      call utmik(nc, x, y, fi, fk)
  gi, gk  [float] Lat. and Lon. (in minutes) of origin
  nc      [int]   Zone number of UTM coordinaets (1-60) for Old Tokyo Datum, 
                  or plus 200 (201-260) for New JGD2000 Datum
  x,  y   [float] Northword(X) and Eastword(Y) coordinates (in km)
                  of the point to be calculated.
  fi, fk  [float] Calculated results. Lat. and Lon. (in minutes) are given.


CVINIT, CVIKEN, CVENIK, CVDINIT, CVDIKEN, CVDENIK
(Map Projection)
NAME
    cvinit / cvdinit - Specify projection method
    cviken / cvdiken - Convert Lat.Lon. into plane rectangular coordinates
    cvenik / cvdenik - Convert plane rectangular coordinates into Lat.Lon.

    Two models of Earth's ellipsoid are implemented:
     1. Bessel's ellipsoid : a = 6377.397155 km , flattening = 1/299.152813
     2. GRS80 ellipsoid :    a = 6378.137 km ,    flattening = 1/298.257222

      Following projection numbers are available for Bessel's ellipsoid, while
    the numbers for GRS80 ellipsoid are plus 200 for each projection.
         0    Japanese Transverse Mercator [JTM]
      1-60    Universal Transverse Mercator [UTM] (with zone number)
        61    Universal Polar Stereographic Projection [UPS] North Pole
        62    Universal Polar Stereographic Projection [UPS] South Pole
        65    Same as UTM but non-standard Central Meridian
        70    Mercator projection
        71    Lambert Conformal Conic projection with 1 Standard Pararrel
        72    Lambert Conformal Conic projection with 2 Standard Pararrels
       100    Lambert Azimuthal Equal Area projection from a sphere
               with the same surface area as the ellipsoid
       109    Lambert Azimuthal Equal Area projection from a sphere
               with the same eqatorial radius as the ellipsoid
       199    Latitude/Longitude in minutes are regarded as distance in km

  (Note) Although X=Y=0. at the origin in general,
       for UTM, X=0. on the equator, and Y=500.(km) on the central meridian,
       for UPS, X=Y=2000.(km) at the pole.

SYNOPSIS

      call cvinit(nc, gi, gk, sp1, sp2)
      call cvdinit(nc, dgi, dgk, dp1, dp2)
  nc        [int]     Projection number
  gi,  gk   [float]   Lat. and Lon. (in minutes) of origin (X=Y=0.)
                       (not refered in cases of UTM/UPS)
  dgi, dgk  [double]   [ditto]
  sp1, sp2  [float]   Lat. (in minutes) of two standard pararells
                       (not refered in cases other than projection 72,272)
  dp1, dp2  [double]   [ditto]

      call cviken(fi, fk, ye, xn)
      call cvdiken(dfi, dfk, dye, dxn)
  fi,  fk   [float]   Lat. and Lon. (in minutes) of the point to be calculated
  dfi, dfk  [double]   [ditto]
  ye,  xn   [float]   Calculated results. Eastword(Y) and Northword(X)
                      coordinates (in km) values are given.
  dye, dxn  [double]   [ditto]

      call cvenik(ye, xn, fi, fk)
      call cvdenik(dye, dxn, dfi, dfk)
  ye,  xn   [float]   Eastword(Y) and Northword(X) coordinates (in km)
                      of the point to be calculated.
  dye, dxn  [double]   [ditto]
  fi,  fk   [float]   Calculated results. Lat. and Lon. (in minutes) are given.
  dfi, dfk  [double]   [ditto]

 (Note)  Be aware that the sequence of coordinates values (X, Y) arguments
       in  cviken / cvenik / cvdiken / cvdenik is different from those
       in  utm / utmik / xyconv / ikconv / nxycnv / nikcnv .