/******************************************************************************* NAME SPHDZ PURPOSE: This function assigns values to the semimajor axis, semiminor axis, and radius of sphere. If the spheroid code is negative, the first two values in the parameter array (parm) are used to define the values as follows: --If parm[0] is a non-zero value and parm[1] is greater than one, the semimajor axis and radius are set to parm[0] and the semiminor axis is set to parm[1]. --If parm[0] is nonzero and parm[1] is greater than zero but less than or equal to one, the semimajor axis and radius are set to parm[0] and the semiminor axis is computed from the eccentricity squared value parm[1]. This algorithm is given below. --If parm[0] is nonzero and parm[1] is equal to zero, the semimajor axis, radius, and semiminor axis are set to parm[0]. --If parm[0] equals zero and parm[1] is greater than zero, the default Clarke 1866 is used to assign values to the semimajor axis, radius and semiminor axis. --If parm[0] and parm[1] equals zero, the semimajor axis and radius are set to 6370997.0 (This value is represented as the last value in the spheroid code array) and the semiminor axis is set to zero. if a spheroid code is zero or greater, the semimajor and semiminor axis are defined by the spheroid code, listed below in the array assignment, and the radius is set to 6370997.0. If the spheroid code is greater than SPHDCT the default spheroid, Clarke 1866, is used to define the semimajor and semiminor axis and radius is set to 6370997.0. The algorithm to define the semiminor axis using the eccentricity squared value is as follows: semiminor = sqrt(1.0 - ES) * semimajor where ES = eccentricity squared PROGRAMMER DATE ---------- ---- T. Mittan MARCH, 1993 ALGORITHM REFERENCES 1. Snyder, John P., "Map Projections--A Working Manual", U.S. Geological Survey Professional Paper 1395 (Supersedes USGS Bulletin 1532), United State Government Printing Office, Washington D.C., 1987. 2. Snyder, John P. and Voxland, Philip M., "An Album of Map Projections", U.S. Geological Survey Professional Paper 1453 , United State Government Printing Office, Washington D.C., 1989. *******************************************************************************/ /* Assign the radius value to the radius argument if a spheroid is assigned */ #define RADVAL 19 #include "cproj.h" #include "proj.h" /* Semi-Major axis of supported Spheroids */ static double major[SPHDCT] = { 6378206.4, /* 0: Clarke 1866 (default) */ 6378249.145, /* 1: Clarke 1880 */ 6377397.155, /* 2: Bessel */ 6378157.5, /* 3: International 1967 */ 6378388.0, /* 4: International 1909 */ 6378135.0, /* 5: WGS 72 */ 6377276.3452, /* 6: Everest */ 6378145.0, /* 7: WGS 66 */ 6378137.0, /* 8: GRS 1980 */ 6377563.396, /* 9: Airy */ 6377304.063, /* 10: Modified Everest */ 6377340.189, /* 11: Modified Airy */ 6378137.0, /* 12: WGS 84 */ 6378155.0, /* 13: Southeast Asia */ 6378160.0, /* 14: Australian National */ 6378245.0, /* 15: Krassovsky */ 6378270.0, /* 16: Hough */ 6378166.0, /* 17: Mercury 1960 */ 6378150.0, /* 18: Modified Mercury 1968 */ 6370997.0, /* 19: Sphere of Radius 6370997 meters*/ 6377483.865, /* 20: Bessel 1841(Namibia) */ 6377298.556, /* 21: Everest (Sabah & Sarawak) */ 6377301.243, /* 22: Everest (India 1956) */ 6377295.664, /* 23: Everest (Malaysia 1969) */ 6377304.063, /* 24: Everest (Malay & Singapr 1948)*/ 6377309.613, /* 25: Everest (Pakistan) */ 6378388.0, /* 26: Hayford */ 6378200.0, /* 27: Helmert 1906 */ 6378160.000, /* 28: Indonesian 1974 */ 6378160.0, /* 29: South American 1969 */ 6378165.0}; /* 30: WGS 60 */ /* Semi-Minor axis of supported Spheroids */ static double minor[SPHDCT] = { 6356583.8, /* 0: Clarke 1866 (default) */ 6356514.86955, /* 1: Clarke 1880 */ 6356078.96284, /* 2: Bessel */ 6356772.2, /* 3: International 1967 */ 6356911.94613, /* 4: International 1909 */ 6356750.519915, /* 5: WGS 72 */ 6356075.4133, /* 6: Everest */ 6356759.769356, /* 7: WGS 66 */ 6356752.31414, /* 8: GRS 1980 */ 6356256.91, /* 9: Airy */ 6356103.039, /* 10: Modified Everest */ 6356034.448, /* 11: Modified Airy */ 6356752.314245, /* 12: WGS 84 */ 6356773.3205, /* 13: Southeast Asia */ 6356774.719, /* 14: Australian National */ 6356863.0188, /* 15: Krassovsky */ 6356794.343479, /* 16: Hough */ 6356784.283666, /* 17: Mercury 1960 */ 6356768.337303, /* 18: Modified Mercury 1968 */ 6370997.0, /* 19: Sphere of Radius 6370997 meters*/ 6356165.382966, /* 20: Bessel 1841(Namibia) */ 6356097.571445, /* 21: Everest (Sabah & Sarawak) */ 6356100.228368, /* 22: Everest (India 1956) */ 6356094.667915, /* 23: Everest (Malaysia 1969) */ 6356103.038993, /* 24: Everest (Malay & Singapr 1948)*/ 6356108.570542, /* 25: Everest (Pakistan) */ 6356911.946128, /* 26: Hayford */ 6356818.169, /* 27: Helmert 1906 */ 6356774.504086, /* 28: Indonesian 1974 */ 6356774.719, /* 29: South American 1969 */ 6356783.287}; /* 30: WGS 60 */ /* Finds the correct ellipsoid axis ---------------------------------*/ void sphdz(isph,parm,r_major,r_minor,radius) long isph; /* spheroid code number */ double *parm; /* projection parameters */ double *r_major; /* major axis */ double *r_minor; /* minor axis */ double *radius; /* radius */ { double t_major; /* temporary major axis */ double t_minor; /* temporary minor axis */ long jsph; /* spheroid code number */ /* if the spheroid code is a negative number, get the semi-major and semi-minor axis from the projection array --------------------------------------------------------------------------*/ if (isph < 0) { t_major = fabs(parm[0]); t_minor = fabs(parm[1]); if (t_major > 0.0) { /* The semimajor axis and the semiminor axis are in the array, assign them directly --------------------------------------------------------------------*/ if (t_minor > 1.0) { *r_major = t_major; *r_minor = t_minor; *radius = t_major; } /* The semimajor axis and the eccentricity squared values are in the array, therefore, the semiminor axis is computed from the eccentricity squared value parm[1] ----------------------------------------------------------------------*/ else if (t_minor > 0.0) { *r_major = t_major; *radius = t_major; *r_minor = (sqrt(1.0 - t_minor)) * t_major; } /* The semiminor axis is zero or less, assign the semimajor axis to the semiminor axis. -----------------------------------------------------------------*/ else { *r_major = t_major; *radius = t_major; *r_minor = t_major; } } /* The sphroid code is to be used to assign the axis -------------------------------------------------*/ else if (t_minor > 0.0) /* t_major = 0 */ /* Assign Clarke 1866 semi-major and semi-minor axis ---------------------------------------------------*/ { *r_major = major[0]; *radius = major[0]; *r_minor = minor[0]; } else /* Assign Spheroid radius to semi-major and semi-minor axis ---------------------------------------------------------*/ { *r_major = major[RADVAL]; *radius = major[RADVAL]; *r_minor = major[RADVAL]; } } /* Use the spheroid code to assign the semi-major and semi-minor axis -----------------------------------------------------------------*/ else /* isph >= 0 */ { jsph = isph; /* The spheroid code is out of range, assign Clarke 1866 ------------------------------------------------------*/ if (jsph > (SPHDCT - 1)) { p_error("Invalid spheroid selection","INFORMATIONAL"); p_error("Reset to 0","INFORMATIONAL"); jsph = 0; } /* Assign the radius argument to the standard radius value -------------------------------------------------------*/ *r_major = major[jsph]; *r_minor = minor[jsph]; *radius = major[RADVAL]; } return; }