【发布时间】：2010-09-28 03:56:02
【问题描述】：
给定一个表示目标的经度/纬度坐标数据库和一系列纬度/经度表示表示“区域”的多边形的边界坐标。使用简单的光线投射算法将纬度/经度视为笛卡尔坐标以确定目标是否在区域内是否安全？
【问题讨论】：
标签： algorithm geolocation
【发布时间】：2010-09-28 03:56:02
【问题描述】：
如果边缘很短，它应该可以很好地工作。但它不会是“正确的”。
想象一个三角形，在赤道上有两个点 a 和 b，还有一个点在北极附近。显然，（几乎）赤道上方经度在 a 和 b 之间的所有点都在三角形中。但是如果我们尝试将 lats/longs 视为笛卡尔坐标，我们会得到完全不同的东西......
编辑：您需要小心跨经度 = 最大/最小的区域，更要小心那些也包含极点的区域
【讨论】：
-
“短”有多短？几米？公里？
-
@Paul Alexander - 你需要多准确？地球的周长约为 40000 公里，所以如果你说的是 1000 公里，那么地球的曲率将低于 1 度。我会说这相当准确，但这取决于您的应用程序。
你可以试试这个
<html xmlns="http://www.w3.org/1999/xhtml">
<head id="Head1" runat="server">
<title></title>
<script type="text/javascript" src="https://maps.googleapis.com/maps/api/js?key=AIzaSyBfaExFhABnh7I3-vTHJjA4TWqtYnMccKE&sensor=false"></script>
<script type="text/javascript" src="http://www.the-di-lab.com/polygon/jquery-1.4.2.min.js"></script>
<script type="text/javascript">
var map;
var boundaryPolygon;
function initialize() {
google.maps.Polygon.prototype.Contains = function (point) {
// ray casting alogrithm http://rosettacode.org/wiki/Ray-casting_algorithm
var crossings = 0,
path = this.getPath();
// for each edge
for (var i = 0; i < path.getLength(); i++) {
var a = path.getAt(i),
j = i + 1;
if (j >= path.getLength()) {
j = 0;
}
var b = path.getAt(j);
if (rayCrossesSegment(point, a, b)) {
crossings++;
}
}
// odd number of crossings?
return (crossings % 2 == 1);
function rayCrossesSegment(point, a, b) {
var px = point.lng(),
py = point.lat(),
ax = a.lng(),
ay = a.lat(),
bx = b.lng(),
by = b.lat();
if (ay > by) {
ax = b.lng();
ay = b.lat();
bx = a.lng();
by = a.lat();
}
if (py == ay || py == by) py += 0.00000001;
if ((py > by || py < ay) || (px > Math.max(ax, bx))) return false;
if (px < Math.min(ax, bx)) return true;
var red = (ax != bx) ? ((by - ay) / (bx - ax)) : Infinity;
var blue = (ax != px) ? ((py - ay) / (px - ax)) : Infinity;
return (blue >= red);
}
};
var mapProp = {
center: new google.maps.LatLng(37.684, -122.591),
zoom: 10,
mapTypeId: google.maps.MapTypeId.ROADMAP
};
map = new google.maps.Map(document.getElementById("googleMap"), mapProp);
google.maps.event.addListener(map, 'click', function (event) {
if (boundaryPolygon!=null && boundaryPolygon.Contains(event.latLng)) {
alert(event.latLng + " is inside the polygon");
} else {
alert(event.latLng + " is outside the polygon");
}
});
}
function btnCheckPointExistence_onclick() {
var latitude = document.getElementById('txtLattitude').value;
var longitude = document.getElementById('txtLongitude').value;
//alert(latitude); alert(longitude)
var myPoint = new google.maps.LatLng(latitude, longitude);
if (boundaryPolygon == null) {
alert("No Polygon");
}
else {
if (boundaryPolygon.Contains(myPoint)) {
alert(myPoint + "is inside the polygon.");
} else {
alert(myPoint + "is outside the polygon.");
}
}
}
function drawPolygon() {
initialize();
var boundary = '77.702866 28.987153, 77.699776 28.978594 ,77.735996 28.974164 ,77.719946 28.99346 ,77.713423 28.994361 ,77.711706 28.990382 ';
var boundarydata = new Array();
var latlongs = boundary.split(",");
for (var i = 0; i < latlongs.length; i++) {
latlong = latlongs[i].trim().split(" ");
boundarydata[i] = new google.maps.LatLng(latlong[1], latlong[0]);
}
boundaryPolygon = new google.maps.Polygon({
path: boundarydata,
strokeColor: "#0000FF",
strokeOpacity: 0.8,
strokeWeight: 2,
fillColor: 'Red',
fillOpacity: 0.4
});
google.maps.event.addListener(boundaryPolygon, 'click', function (event) {
if (boundaryPolygon.Contains(event.latLng)) {
alert(event.latLng + "is inside the polygon.");
} else {
alert(event.latLng + "is outside the polygon.");
}
});
map.setZoom(14);
map.setCenter(boundarydata[0]);
boundaryPolygon.setMap(map);
}
</script>
</head>
<body onload="initialize();drawPolygon();">
<form id="form1" runat="server">
<input type="text" value="28.982005022927137" name="txtLattitude" id="txtLattitude" placeholder="Enter Latitude">
<input type="text" value="77.72157669067383" name="txtLongitude" id="txtLongitude" placeholder="Enter Langitude">
<input type="button" value="Submit" name="submit" id="submit-text" value="submit" onclick="btnCheckPointExistence_onclick();">
<h3>Check If Point Exist in a Polygon</h3>
<span id="spnMsg" style="font-family: Arial; text-align: center; font-size: 14px; color: red;">this is message</span>
<br />
<br />
<div id="googleMap" style="width: auto; height: 500px;">
</div>
</form>
</body>
</html>
您可以查看Demo here
【讨论】：
不，这样做不安全。您会在两极周围以及 0 经度区域遇到异常情况。
看看 NASA 世界风代码。尽管您必须小心，但它们在地球表面上具有多边形的表示形式-某些表示形式将代表这些异常，而另一些则不会。
【讨论】：