2019独角兽企业重金招聘Python工程师标准>>> 
实现原理，根据球面计算公式，计算两点间的距离。
由于地球是一个不规则的椭球体，当前计算方法，还是有很大的误差在。但能基本实现附近查询。
函数：返回两个坐标之间的公里数
/****mssql函数：计算经纬度****************************/
create FUNCTION [GetDistance]
(
@GPSLng DECIMAL(12,6),
@GPSLat DECIMAL(12,6),
@Lng DECIMAL(12,6),
@Lat DECIMAL(12,6)
)
RETURNS DECIMAL(12,4)
AS
BEGIN
DECLARE @result DECIMAL(12,4)
SELECT @result=6371.004*ACOS(SIN(@GPSLat/180*PI())*SIN(@Lat/180*PI())+COS(@GPSLat/180*PI())*COS(@Lat/180*PI())*COS((@GPSLng-@Lng)/180*PI()))
RETURN @result
END使用示例：
--查找数据库中5公里内的记录
select * from (
SELECT dbo.GetDistance(36.6513750000,116.9046970000,Latitude,Longitude)AS dis FROM cards
) B
where dis<5
order by dis asc;
sql查询语句：
SELECT *,GetDistance(某一点的经度,某一点的纬度,数据库中经度,数据库中纬度)AS dis FROM 表名 dis<5
查找数据库中5公里内的记录
注：
6371.004 地球半径
6371.004*ACOS(SIN(@GPSLat/180*PI())*SIN(@Lat/180*PI())+COS(@GPSLat/180*PI())*COS(@Lat/180*PI())*COS((@GPSLng-@Lng)/180*PI()))
计算公式 可查球面弧长计算公式的详细解释
分页获取
//获取列表
string sql = string.Format(@"select * from (
SELECT *, dbo.GetDistance({0},{1},Latitude,Longitude)AS dis FROM cards
) B
where IsEnable=1 and dis>0
order by dis asc
OFFSET {2} ROW FETCH NEXT 20 ROW ONLY", lat, lng, skip);
List list = _cards.Context.Database.SqlQuery(sql).ToList();
原文出处：https://blog.csdn.net/ecdyf1989/article/details/8077943/