判断指定的经纬度坐标点是否落在指定的多边形区域内

描述: 在矢量多边形区域中,一个坐标点的位置是否在区域内 算法: C#代码,适用于任意多边形(凹凸多边形),但是没有考虑实际误差范围的情况(应用在实际问题解决中,接近区域一定范围是可以忽略的,这个误差范围考虑后算法的复杂度会加倍,所以没有做这方面的考虑)

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace p
{
    class Program
    {
        static void Main（string[] args）
        {
            Point p = new Point（112.3839540000 37.9901770000）;

            Point[] ps = new Point[] { new Point（112.4091887500 38.0200324167） new Point（112.4668669727 37.7316897817） new Point（112.6701140430 37.7686085000） new Point（112.6234221484 38.0070487978） };

            Console.WriteLine（“x:“ + p.X + “y:“ + p.Y）;

            Console.WriteLine（“result:“ + （isPtInPoly（p.X p.Y ps） ? “yes“ : “no“））;

            Console.ReadKey（）;
        }

        /// 

        ///  判断指定的经纬度坐标点是否落在指定的多边形区域内
        /// 

        /// 指定点的经度
        /// 指定点的纬度
        /// 指定多边形区域各个节点坐标
        /// True 落在范围内 False 不在范围内
        public static bool isPtInPoly（double ALon double ALat Point[] APoints）
        {
            int iSum iCount iIndex;
            double dLon1 = 0 dLon2 = 0 dLat1 = 0 dLat2 = 0 dLon;

            if （APoints.Length < 3）
            {
                return false;
            }

            iSum = 0;
            iCount = APoints.Length;

            for （iIndex = 0; iIndex < iCount - 1; iIndex++）
            {
                if （iIndex == iCount - 1）
                {
                    dLon1 = APoints[iIndex].X;
                    dLat1 = APoints[iIndex].Y;
                    dLon2 = APoints[0].X;
                    dLat2 = APoints[0].Y;
                }
                else
                {
                    dLon1 = APoints[iIndex].X;
                    dLat1 = APoints[iIndex].Y;
                    dLon2 = APoints[iIndex + 1].X;
                    dLat2 = APoints[iIndex + 1].Y;
                }

                if （（（ALat >= dLat1） && （ALat < dLat2）） || （（ALat >= dLat2） && （ALat < dLat1）））
                {
                    if （Math.Abs（dLat1 - dLat2） > 0）
                    {
                        dLon = dLon1 - （（dLon1 - dLon2） * （dLat1 - ALat）） / （dLat1 - dLat2）;
                        if （dLon < ALon）
                            iSum++;
                    }
                }
            }

            if （（iSum % 2） != 0）
                return true;
            return false;
        }
    }

    /// 

    /// 坐标
    /// 

    public class Point
    {
        private Double x;
        private Double y;

        public Point（）
        {
            this.x = 0;
            this.y = 0;
        }
        public Point（double _x double _y）
        {
            this.x = _x;
            this.y = _y;
        }

        public Double X
        {
            get { return x; }
            set { x = value; }
        }

        public Double Y
        {
            get { return y; }
            

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3164  2012-06-16 14:56  Program.cs