判断轨迹点是否在区域内

/**
     * 判断经纬度点是否在区域内
     * @param px
     * @param py
     * @param polygonXA
     * @param polygonYA
     * @return
     */
    public static boolean isPointInPolygon(double px, double py,
            ArrayList<Double> polygonXA, ArrayList<Double> polygonYA) {
        boolean isInside = false;
        double ESP = 1e-9;
        int count = 0;
        double linePoint1x;
        double linePoint1y;
        double linePoint2x = 180;
        double linePoint2y;
 
        linePoint1x = px;
        linePoint1y = py;
        linePoint2y = py;
 
        for (int i = 0; i < polygonXA.size() - 1; i++) {
            double cx1 = polygonXA.get(i);
            double cy1 = polygonYA.get(i);
            double cx2 = polygonXA.get(i + 1);
            double cy2 = polygonYA.get(i + 1);
            // 如果目标点在任何一条线上
            if (isPointOnLine(px, py, cx1, cy1, cx2, cy2)) {
                return true;
            }
            // 如果线段的长度无限小(趋于零)那么这两点实际是重合的，不足以构成一条线段
            if (Math.abs(cy2 - cy1) < ESP) {
                continue;
            }
            // 第一个点是否在以目标点为基础衍生的平行纬度线
            if (isPointOnLine(cx1, cy1, linePoint1x, linePoint1y, linePoint2x,
                    linePoint2y)) {
                // 第二个点在第一个的下方,靠近赤道纬度为零(最小纬度)
                if (cy1 > cy2)
                    count++;
            }
            // 第二个点是否在以目标点为基础衍生的平行纬度线
            else if (isPointOnLine(cx2, cy2, linePoint1x, linePoint1y,
                    linePoint2x, linePoint2y)) {
                // 第二个点在第一个的上方,靠近极点(南极或北极)纬度为90(最大纬度)
                if (cy2 > cy1)
                    count++;
            }
            // 由两点组成的线段是否和以目标点为基础衍生的平行纬度线相交
            else if (isIntersect(cx1, cy1, cx2, cy2, linePoint1x, linePoint1y,
                    linePoint2x, linePoint2y)) {
                count++;
            }
        }
        if (count % 2 == 1) {
            isInside = true;
        }
 
        return isInside;
    }
    public static boolean isPointOnLine(double px0, double py0, double px1,
            double py1, double px2, double py2) {
        boolean flag = false;
        double ESP = 1e-9;// 无限小的正数
        if ((Math.abs(Multiply(px0, py0, px1, py1, px2, py2)) < ESP)
                && ((px0 - px1) * (px0 - px2) <= 0)
                && ((py0 - py1) * (py0 - py2) <= 0)) {
            flag = true;
        }
        return flag;
    }
    public static double Multiply(double px0, double py0, double px1, double py1,
            double px2, double py2) {
        return ((px1 - px0) * (py2 - py0) - (px2 - px0) * (py1 - py0));
    }
    public static boolean isIntersect(double px1, double py1, double px2, double py2,
            double px3, double py3, double px4, double py4) {
        boolean flag = false;
        double d = (px2 - px1) * (py4 - py3) - (py2 - py1) * (px4 - px3);
        if (d != 0) {
            double r = ((py1 - py3) * (px4 - px3) - (px1 - px3) * (py4 - py3))
                    / d;
            double s = ((py1 - py3) * (px2 - px1) - (px1 - px3) * (py2 - py1))
                    / d;
            if ((r >= 0) && (r <= 1) && (s >= 0) && (s <= 1)) {
                flag = true;
            }
        }
        return flag;
    }

调用示例：

ArrayList<Double> arrlon = new ArrayList<Double>();
                    ArrayList<Double> arrlat = new ArrayList<Double>();
                    if(!StringUtils.isEmpty(regionRange.getRegionRange())&&regionRange.getRegionRange()!=null) {
                        JSONArray arrayObj = JSONArray.fromObject(regionRange.getRegionRange());
                        for(int j=0;j<arrayObj.size();j++) {
                            JSONObject objArry= arrayObj.getJSONObject(j);
                            String lon = objArry.getString("lon");
                            String lat = objArry.getString("lat");
                            arrlon.add(Double.parseDouble(lon));
                            arrlat.add(Double.parseDouble(lat));
                        }
                    }

//比较区域坐标点信息和轨迹数据判断车辆是否在区域内
                    boolean areaStatus1;
                    boolean areaStatus2;
                    AlarmBasicInfoEntity basic = new AlarmBasicInfoEntity();
                    for(int k=0;k<trackList.size()-1;k++) {
                        areaStatus1 = LocationUtil.isPointInPolygon(trackList.get(k).getLng(), trackList.get(k).getLat(), arrlon, arrlat);                    
                        areaStatus2 = LocationUtil.isPointInPolygon(trackList.get(k+1).getLng(), trackList.get(k+1).getLat(), arrlon, arrlat);                    
                        if(areaStatus1&&!areaStatus2) {    //areaStatus1为真areaStatus2为假则为出区域
                            if(StringUtils.isEmpty(basic.getInAreaTime())){
                                continue;
                            }
                            basic.setOutAreaTime(trackList.get(k+1).getGpsTime());
                            result.add(basic);
                        }else if(areaStatus2&&!areaStatus1) {    //areaStatus1为假areaStatus2为真则为进区域
                            basic = new AlarmBasicInfoEntity();
                            basic.setInAreaTime(trackList.get(k+1).getGpsTime());
                            basic.setLicense(basicInfo.getCarNo());
                            basic.setVehId(areaList.get(i).getVehId());
                            basic.setRegionName(regionRange.getRegionName());
                        }
                        
                    }