经过这两天的研究，我想有些东西有必要写一下，同时也帮助需要的人
这是一个查找附近的人的一个算法,在网上找了这篇文章 http://blog.csdn.net/dyllove98/article/details/9795815，他的算最小正方形的四个顶点有点问题。
第一步 ，我们试想一下，如果我们要查找附近五公里的人，是不是最先想到根据自己的经纬度(如果连经纬度都不知道怎么回事，我想你应该知道度娘)，再往数据库里面获取其他人的经纬度，根据两个经纬度获取他们的距离，
如果距离小于5公里的话，那就符合条件。
获取距离用两种情况，第一种就是全部查出来在代码里面进行算，第二种情况就是在sql 里面写个函数，在数据库里面算
代码算距离算法：
public static double GetDistance(Degree Degree1, Degree Degree2)
{
double radLat1 = radians(Degree1.Y);
double radLat2 = radians(Degree2.Y);
double a = radLat1 - radLat2;
double b = radians(Degree1.X) - radians(Degree2.X);
double s = 2 * Math.Asin(Math.Sqrt(Math.Pow(Math.Sin(a / 2), 2) +
Math.Cos(radLat1) * Math.Cos(radLat2) * Math.Pow(Math.Sin(b / 2), 2)));
s = s * EARTH_RADIUS;
s = Math.Round(s * 10000) / 10000;
return s;
}
radians是个角度转弧度的算法
private static double degrees(double d)
{
return d * (180 / Math.PI);
}
如果不知道里面的算法是个怎样的原理，你可以百度相关的三角函数公式
sql 距离算法
SET ANSI_NULLS ON
GO
SET QUOTED_IDENTIFIER ON
GO
ALTER function [dbo].[fnGetDistance](@LatBegin float, @LngBegin float, @LatEnd float, @LngEnd float) returns float
as
begin
--距离(千米)
declare @Distance float;
declare @EARTH_RADIUS float;
set @EARTH_RADIUS = 6378.137;
declare @RadLatBegin float,
@RadLatEnd float,
@RadLatDiff float,
@RadLngDiff float;
set @RadLatBegin = @LatBegin * pi()/180.0;
set @RadLatEnd = @LatEnd * pi()/180.;
set @RadLatDiff = @RadLatBegin - @RadLatEnd;
set @RadLngDiff = @LngBegin * pi()/180.0 - @LngEnd * pi()/180.0;
set @Distance = 2 * asin(sqrt(power(sin(@RadLatDiff/2), 2) + cos(@RadLatBegin)
* cos(@RadLatEnd) * power(sin(@RadLngDiff/2), 2)));
set @Distance = @Distance * @EARTH_RADIUS;
--set @Distance = Round(@Distance * 10000) / 10000
return @Distance
end
如果你跟我上面的想法一致的话，我只能对你说处理小数据量还可以。
例如：我们现在要做一个停车管理系统，需车位的可以查附近的供车位的，供车位的可以查附近需车位的，如果数据库里面有一万条供车位的经纬度，需车位访问服务器，服务器要进行十万次距离换算，如果同时有十个人进行访问的话，那服务器被你轻而易举的百万数据量访问了。
那有什么办法可以解决呢？
先看这张图，我们知道原点的坐标也就是我们的经纬度，也在知道半径(附近几公里)，就可以得到圆的最小正方形，从而根据公式可以得到 正方形的四个顶点 就晓得了 最小经纬度和最大经纬度
把我们得到的经纬度往数据库里面一查 ，就能轻然一举的把符合范围内的人查
出来。
根据经纬度和距离，得到四个顶点的距离算法
private static void GetlatLon(double GLON, double GLAT, double distance, double angle, out double newLon,out double newLat)
{
double Ea = 6378137; // 赤道半径
double Eb = 6356725; // 极半径
double dx = distance * Math.Sin(angle * Math.PI / 180.0);
double dy = distance * Math.Cos(angle * Math.PI / 180.0);
//double ec = 6356725 + 21412 * (90.0 - GLAT) / 90.0;
// 21412 是赤道半径与极半径的差
double ec = Eb + (Ea - Eb) * (90.0 - GLAT) / 90.0;
double ed = ec * Math.Cos(GLAT * Math.PI / 180);
newLon = (dx / ed + GLON * Math.PI / 180.0) * 180.0 / Math.PI;
newLat = (dy / ec + GLAT * Math.PI / 180.0) * 180.0 / Math.PI;
}
public static Degree[] GetRectRange(double centorLogitude,double centorlatitude, double distance)
{
double temp = 0.0;
double maxLatitude;
double minLatitude;
double maxLongitude;
double minLongitude;
GetlatLon(centorLogitude, centorlatitude, distance, 0, out temp, out maxLatitude);
GetlatLon(centorLogitude, centorlatitude, distance, 180, out temp, out minLatitude);
GetlatLon(centorLogitude, centorlatitude, distance, 90, out maxLongitude, out temp);
GetlatLon(centorLogitude, centorlatitude, distance, 270, out minLongitude, out temp);
maxLatitude = Math.Round(maxLatitude,6);
minLatitude = Math.Round(minLatitude,6);
maxLongitude = Math.Round(maxLongitude,6);
minLongitude = Math.Round(minLongitude,6);
return new Degree[] {
new Degree(minLongitude,maxLatitude),//left-top
new Degree(minLongitude,minLatitude),//left-bottom
new Degree(maxLongitude,maxLatitude),//right-top
new Degree(maxLongitude,minLatitude) //right-bottom
};
}
这些东西经过测试的 相差不到 1,2米
还是那就话如有不懂的地方 请联系我 QQ:209229923，或在下发为我留言