google map 计算地图面积方法

花了几个小时把js的google计算地图面积的算法改成了c# 的。

  class Program

    {

        static void Main(string[] args)

        {

            // a = new qq.maps.LatLng(39.93, 116.44004334);

            //b = new qq.maps.LatLng(39.93, 116.35421264);

            //c = new qq.maps.LatLng(39.896775, 116.35421264);

            //d = new qq.maps.LatLng(39.89, 116.44004334);

            //29869510.924133233平方米

            List<LatLng> a = new List<LatLng>();

            a.Add(new LatLng() { Latitude = 39.93, Longitude = 116.44004334 });

            a.Add(new LatLng() { Latitude = 39.93, Longitude = 116.35421264 });

            a.Add(new LatLng() { Latitude = 39.896775, Longitude = 116.35421264 });

            a.Add(new LatLng() { Latitude = 39.89, Longitude = 116.44004334 });

            double area = computeArea(a, );

        }

        static double computeArea(List<LatLng> x, double y)

        {

            return Math.Abs(computeSignedArea(x, y));

        }

        static double computeSignedArea(List<LatLng> a, double b)

        {

            double c = ;

            if (b > )

            {

                c = b;

            }

            double g = a.Count - ;

            double e = ;

            LatLng d = a[];

            for (int i = ; i < g; i++)

            {

                e += Zm(d, a[i], a[i + ]);

            }

            return e * c * c;

        }

        static double Zm(LatLng a, LatLng b, LatLng c)

        {

            return Tm(a, b, c) * Um(a, b, c);

        }

        static double Tm(LatLng a, LatLng b, LatLng c)

        {

            List<LatLng> d = new List<LatLng>() { a, b, c, a };

            List<double> e = new List<double>();

            double f = ;

            for (int i = ; i < d.Count - ; i++)

            {

                e.Add(Uf(d[i], d[i + ]));

                f += e[i];

            }

            f /= ;

            double g = Math.Tan(f / );

            for (int j = ; j < d.Count - ; j++)

            {

                g *= Math.Tan((f - e[j]) / );

            }

            return  * Math.Atan(Math.Sqrt(Math.Abs(g)));

        }

        static double Um(LatLng a, LatLng b, LatLng c)

        {

            List<LatLng> d = new List<LatLng>() { a, b, c };

            double[,] e = new double[, ];

            for (int i = ; i < 3; i++)

            {

                LatLng f = d[i];

                double f1 = re(f);

                double f2 = se(f);

                e[i, ] = Math.Cos(f1) * Math.Cos(f2);

                e[i, ] = Math.Cos(f1) * Math.Sin(f2);

                e[i, ] = Math.Sin(f1);

            }

            return  < e[, ] * e[, ] * e[, ] + e[, ] * e[, ] * e[, ] + e[, ] * e[, ] * e[, ] - e[, ] * e[, ] * e[, ] - e[, ] * e[, ] * e[, ] - e[, ] * e[, ] * e[, ] ?  : -;

        }

        static double Uf(LatLng a, LatLng b)

        {

            var c = re(a);

            var d = re(b);

            return  * Math.Asin(Math.Sqrt(Math.Pow(Math.Sin((c - d) / ), ) + Math.Cos(c) * Math.Cos(d) * Math.Pow(Math.Sin((se(a) - se(b)) / ), )));

        }

        static double re(LatLng a)

        {

            return Rd(a.Latitude);

        }

        static double se(LatLng a)

        {

            return Rd(a.Longitude);

        }

        static double Rd(double a)

        {

            return Math.PI /  * a;

        }

    }

    partial struct LatLng

    {

        private double _Longitude;

        public double Longitude

        {

            get

            {

                return _Longitude;

            }

            set

            {

                _Longitude = value;

            }

        }

        private double _Latitude;

        public double Latitude

        {

            get

            {

                return _Latitude;

            }

            set

            {

                _Latitude = value;

            }

        }

    }

 double computeDistanceBetween(LngLat a, LngLat b, double c)

        {

            return Uf(a, b) * c;

        }

        public double ComputeLength(List<LngLat> a, double b = )

        {

            if (a.Count < )

            {

                return ;

            }

            double c = ;

            if (b > )

            {

                c = b;

            }

            double g = a.Count - ;

            double d = ;

            for (var e = ; e < g; ++e)

            {

                d += computeDistanceBetween(a[e], a[e + ], c);

            }

            return d;

        }

google map 计算地图面积方法的相关教程结束。