判斷點是否在面內

	public class Point {
	/**
	* 是否有 橫斷<br/> 參數爲四個點的座標
	*
	* @param px1
	* @param py1
	* @param px2
	* @param py2
	* @param px3
	* @param py3
	* @param px4
	* @param py4
	* @return
	*/
	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;
	}
	/**
	* 目標點是否在目標邊上邊上<br/>
	*
	* @param px0
	* 目標點的經度座標
	* @param py0
	* 目標點的緯度座標
	* @param px1
	* 目標線的起點(終點)經度座標
	* @param py1
	* 目標線的起點(終點)緯度座標
	* @param px2
	* 目標線的終點(起點)經度座標
	* @param py2
	* 目標線的終點(起點)緯度座標
	* @return
	*/
	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));
	}
	/**
	* 判斷目標點是否在多邊形內(由多個點組成)<br/>
	*
	* @param px
	* 目標點的經度座標
	* @param py
	* 目標點的緯度座標
	* @param polygonXA
	* 多邊形的經度座標集合
	* @param polygonYA
	* 多邊形的緯度座標集合
	* @return
	*/
	public static boolean isPointInPolygon(double px, double py,
	List<Double> polygonXA, List<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;
	}
	}