判断一个点是否落在多边形内
2009-11-30 21:23
323 查看
一、背景:
如何判断一个指定的经纬度点是否落在一个多边形区域内?
二、实现代码(delphi)
Code
Type
TMyPoint = packed record
X : double;
Y : double;
end;
{*------------------------------------------------------------------------------
判断指定的经纬度坐标点是否落在指定的多边形区域内
@param ALon 指定点的经度
@param ALat 指定点的纬度
@param APoints 指定多边形区域各个节点坐标
@return True 落在范围内 False 不在范围内
------------------------------------------------------------------------------*}
function IsPtInPoly(ALon, ALat: double; APoints: array of TMyPoint): Boolean;
var
iSum, iCount, iIndex: Integer;
dLon1, dLon2, dLat1, dLat2, dLon: double;
begin
Result := False;
if (Length(APoints) < 3) then
begin
Result := False;
Exit;
end;
iSum := 0;
iCount := Length(APoints);
for iIndex :=0 to iCount - 1 do
begin
if (iIndex = iCount - 1) then
begin
dLon1 := APoints[iIndex].X;
dLat1 := APoints[iIndex].Y;
dLon2 := APoints[0].X;
dLat2 := APoints[0].Y;
end
else
begin
dLon1 := APoints[iIndex].X;
dLat1 := APoints[iIndex].Y;
dLon2 := APoints[iIndex + 1].X;
dLat2 := APoints[iIndex + 1].Y;
end;
if ((ALat >= dLat1) and (ALat < dLat2)) or ((ALat>=dLat2) and (ALat < dLat1)) then
begin
if (abs(dLat1 - dLat2) > 0) then
begin
dLon := dLon1 - ((dLon1 -dLon2) * (dLat1 -ALat)) / (dLat1 - dLat2);
if (dLon < ALon) then
Inc(iSum);
end;
end;
end;
if (iSum mod 2 <> 0) then
Result := True;
end;
如何判断一个指定的经纬度点是否落在一个多边形区域内?
二、实现代码(delphi)
Code
Type
TMyPoint = packed record
X : double;
Y : double;
end;
{*------------------------------------------------------------------------------
判断指定的经纬度坐标点是否落在指定的多边形区域内
@param ALon 指定点的经度
@param ALat 指定点的纬度
@param APoints 指定多边形区域各个节点坐标
@return True 落在范围内 False 不在范围内
------------------------------------------------------------------------------*}
function IsPtInPoly(ALon, ALat: double; APoints: array of TMyPoint): Boolean;
var
iSum, iCount, iIndex: Integer;
dLon1, dLon2, dLat1, dLat2, dLon: double;
begin
Result := False;
if (Length(APoints) < 3) then
begin
Result := False;
Exit;
end;
iSum := 0;
iCount := Length(APoints);
for iIndex :=0 to iCount - 1 do
begin
if (iIndex = iCount - 1) then
begin
dLon1 := APoints[iIndex].X;
dLat1 := APoints[iIndex].Y;
dLon2 := APoints[0].X;
dLat2 := APoints[0].Y;
end
else
begin
dLon1 := APoints[iIndex].X;
dLat1 := APoints[iIndex].Y;
dLon2 := APoints[iIndex + 1].X;
dLat2 := APoints[iIndex + 1].Y;
end;
if ((ALat >= dLat1) and (ALat < dLat2)) or ((ALat>=dLat2) and (ALat < dLat1)) then
begin
if (abs(dLat1 - dLat2) > 0) then
begin
dLon := dLon1 - ((dLon1 -dLon2) * (dLat1 -ALat)) / (dLat1 - dLat2);
if (dLon < ALon) then
Inc(iSum);
end;
end;
end;
if (iSum mod 2 <> 0) then
Result := True;
end;
相关文章推荐
- Java 判断一个点是否在多边形区域内【转】
- 如何判断一个指定的经纬度点是否落在一个多边形内
- 判断一个点是否在多边形区域内
- 判断一个点是否在多边形内部 [2] 射线法实现
- 如何判断一个点是否在一个多边形内部(有向面积)
- 判断一个点是否在多边形内部
- openlayers2 开发如何判断一个marker所在的点是否在一个多边形内部
- 判断一个点是否在多边形内部
- 如何判断一个点是否在一个多边形内部(转)
- java实现判断一个经纬度坐标是否在一个多边形内(经自己亲测)
- [转]判断一个点是否落在多边形内
- POJ 3335 判断一个多边形是否存在核
- Mysql中判断一个点是否落在多边形内。
- 判断一个点是否在多边形内(java)
- 利用百度api判断已知经纬度的一个点是否在多边形内
- 判断一个点是否在多边形区域内
- 判断一个点是否位于一个凸多边形内的方法
- geotools判断一个点是否在多边形上
- MySQL中判断一个点是否落在多边形内
- javascript 判断一个点是否在多边形内(向量的使用)