コード例 #1
0
        /// <summary>
        /// Initializes the transform using the parameters from the specified coordinate system information
        /// </summary>
        /// <param name="projInfo">A ProjectionInfo class contains all the standard and custom parameters needed to initialize this transform</param>
        protected override void OnInit(ProjectionInfo projInfo)
        {
            if (projInfo.Parameters.ContainsKey("W"))
            {
                _w = projInfo.ParamD("W");
                if (_w <= 0) throw new ProjectionException(27);
            }
            else
            {
                _w = .5;
            }
            if(projInfo.Parameters.ContainsKey("M"))
            {
                _m = projInfo.ParamD("M");
                if(_m <= 0)throw new ProjectionException(27);
            }
            else
            {
                _m = 1;
            }
            _rm = 1/_m;
            _m /= _w;


        }
コード例 #2
0
        /// <summary>
        /// Initializes the transform using the parameters from the specified coordinate system information
        /// </summary>
        /// <param name="projInfo">A ProjectionInfo class contains all the standard and custom parameters needed to initialize this transform</param>
        protected override void OnInit(ProjectionInfo projInfo)
        {
            if ((_h = projInfo.ParamD("h")) <= 0) throw new ProjectionException(-30);
	        if (Phi0 == 0) throw new ProjectionException(-46);
            _radiusG = 1 + (_radiusG1 = _h / A);
            _c = _radiusG * _radiusG - 1.0;
	        if (IsElliptical) {
                _radiusP = Math.Sqrt(OneEs);
                _radiusP2 = OneEs;
                _radiusPInv2 = ROneEs;

	        } else {
                _radiusP = _radiusP2 = _radiusPInv2 = 1.0;

	        }
        }
コード例 #3
0
 /// <summary>
 /// Initializes the transform using the parameters from the specified coordinate system information
 /// </summary>
 /// <param name="projInfo">A ProjectionInfo class contains all the standard and custom parameters needed to initialize this transform</param>
 protected override void OnInit(ProjectionInfo projInfo)
 {
     double t=0;
     if (projInfo.StandardParallel1 != null) t = projInfo.StandardParallel1.Value * Math.PI / 180;
     t = projInfo.ParamD("lat_ts")*Math.PI/180;
     if((K0 = Math.Cos(t)) < 0) throw new ProjectionException(-24);
     if (!IsElliptical) return;
     t = Math.Sin(t);
     K0 /= Math.Sqrt(1- Es * t * t);
     E = Math.Sqrt(Es);
     _apa = Proj.Authset(Es);
     _qp = Proj.Qsfn(1, E, OneEs);
 }
コード例 #4
0
        /// <summary>
        /// Initializes the transform using the parameters from the specified coordinate system information
        /// </summary>
        /// <param name="projInfo">A ProjectionInfo class contains all the standard and custom parameters needed to initialize this transform</param>
        protected override void OnInit(ProjectionInfo projInfo)
        {
            double con;
            double f;
            double d;
            double toRadians = projInfo.GeographicInfo.Unit.Radians;
            _rot = projInfo.ParamI("no_rot") == 0;
            bool azi = projInfo.ParamD("alpha") != 0.0;
            if (azi)
            {
                _lamc = projInfo.ParamD("lonc")*toRadians;
                _alpha = projInfo.ParamD("alpha") *toRadians;
                if (Math.Abs(_alpha) < Tol ||
                    Math.Abs(Math.Abs(Phi0) - HalfPi) <= Tol ||
                    Math.Abs(Math.Abs(_alpha) - HalfPi) <= Tol)
                    throw new ProjectionException(32);
            }
            else
            {
                _lam1 = projInfo.GetLam1();
                _phi1 = projInfo.GetPhi1();
                _lam2 = projInfo.GetLam2();
                _phi2 = projInfo.GetPhi2();
                if(Math.Abs(_phi1 - _phi2) <= Tol ||
                    (con = Math.Abs(_phi1)) <= Tol ||
                    Math.Abs(con - HalfPi) <= Tol ||
                    Math.Abs(Math.Abs(Phi0) - HalfPi) <= Tol ||
                    Math.Abs(Math.Abs(_phi2) - HalfPi) <= Tol)
                {
                    throw new ProjectionException(33);
                }
            }
            _ellips = Es > 0;
            double com = _ellips ? Math.Sqrt(OneEs) : 1;
            if (Math.Abs(Phi0) > EPS10)
            {
                double sinph0 = Math.Sin(Phi0);
                double cosph0 = Math.Cos(Phi0);
                if (_ellips)
                {
                    con = 1 - Es*sinph0*sinph0;
                    _bl = cosph0*cosph0;
                    _bl = Math.Sqrt(1 + Es*_bl*_bl/OneEs);
                    _al = _bl*K0*com/con;
                    d = _bl*com/(cosph0*Math.Sqrt(con));
                }
                else
                {
                    _bl = 1;
                    _al = K0;
                    d = 1/cosph0;
                }


                if ((f = d*d - 1) <= 0)
                {
                    f = 0;
                }
                else
                {
                    f = Math.Sqrt(f);
                    if (Phi0 < 0) f = -f;
                }
                _el = f += d;
                if (_ellips)
                {
                    _el *= Math.Pow(Proj.Tsfn(Phi0, sinph0, E), _bl);
                }
                else
                {
                    _el *= TSFN0(Phi0);
                }
            }
            else
            {
                _bl = 1/com;
                _al = K0;
                _el = d = f = 1;
            }
            if(azi)
            {
                _gamma = Math.Asin(Math.Sin(_alpha)/d);
                Lam0 = _lamc - Math.Asin((.5*(f - 1/f))*Math.Tan(_gamma))/_bl;
            }
            else
            {
                double h;
                double l;
                if(_ellips)
                {
                    h = Math.Pow(Proj.Tsfn(_phi1, Math.Sin(_phi1), E), _bl);
                    l = Math.Pow(Proj.Tsfn(_phi2, Math.Sin(_phi2), E), _bl);
                }
                else
                {
                    h = TSFN0(_phi1);
                    l = TSFN0(_phi2);
                }
                f = _el/h;
                double p = (l - h)/(l + h);
                double j = _el*_el;
                j = (j - l*h)/(j + l*h);
                if ((con = _lam1 - _lam2) < -Math.PI)
                {
                    _lam2 -= Math.PI*2;
                }
                else if (con > Math.PI)
                {
                    _lam2 += Math.PI*2;
                }
                Lam0 = Proj.Adjlon(.5*(_lam1 + _lam2) - Math.Atan(j*Math.Tan(.5*_bl*(_lam1 - _lam2))/p)/_bl);
                _gamma = Math.Atan(2*Math.Sin(_bl*Proj.Adjlon(_lam1 - Lam0))/(f - 1/f));
                _alpha = Math.Asin(d*Math.Sin(_gamma));
            }
            _singam = Math.Sin(_gamma);
            _cosgam = Math.Cos(_gamma);
            if(projInfo.ParamI("rot_conv") != 0)
            {
                f = _gamma;
            }
            else
            {
                f = _alpha;
            }
            _sinrot = Math.Sin(f);
            _cosrot = Math.Cos(f);
            if(projInfo.ParamI("no_uoff") != 0)
            {
                _u0 = 0;
            }
            else
            {
                _u0 = Math.Abs(_al*Math.Atan(Math.Sqrt(d*d - 1)/_cosrot)/_bl);
               
            }
            if (Phi0 < 0) _u0 = -_u0;
        }