我正在尝试计算地图上两个位置之间的距离。 我存储了我的数据:经度,纬度,X POS,Y POS。
我以前使用过以下代码段。
DECLARE @orig_lat DECIMAL
DECLARE @orig_lng DECIMAL
SET @orig_lat=53.381538 set @orig_lng=-1.463526
SELECT *,
3956 * 2 * ASIN(
SQRT( POWER(SIN((@orig_lat - abs(dest.Latitude)) * pi()/180 / 2), 2)
+ COS(@orig_lng * pi()/180 ) * COS(abs(dest.Latitude) * pi()/180)
* POWER(SIN((@orig_lng - dest.Longitude) * pi()/180 / 2), 2) ))
AS distance
--INTO #includeDistances
FROM #orig dest
但我不相信这些数据,似乎会给出稍微不准确的结果。
如果您需要,可以提供一些示例数据
Latitude Longitude Distance
53.429108 -2.500953 85.2981833133896
任何人都可以帮我解决我的代码问题,我不介意你是否想要修复我已有的东西,如果你有一种新方法可以实现这一目标,那就太棒了。
请说明您的结果所在的计量单位。
答案 0 :(得分:110)
由于您使用的是SQL Server 2008,因此您可以使用geography
数据类型,该类型专为此类数据而设计:
DECLARE @source geography = 'POINT(0 51.5)'
DECLARE @target geography = 'POINT(-3 56)'
SELECT @source.STDistance(@target)
给出
----------------------
538404.100197555
(1 row(s) affected)
告诉我们距伦敦(近)爱丁堡(近)约538公里。
当然,首先要学习一些知识,但是一旦你知道它,它就比实施你自己的Haversine计算容易得多;另外你还有很多功能。
如果您想保留现有的数据结构,仍然可以使用STDistance
,使用Point
方法构建合适的geography
实例:
DECLARE @orig_lat DECIMAL(12, 9)
DECLARE @orig_lng DECIMAL(12, 9)
SET @orig_lat=53.381538 set @orig_lng=-1.463526
DECLARE @orig geography = geography::Point(@orig_lat, @orig_lng, 4326);
SELECT *,
@orig.STDistance(geography::Point(dest.Latitude, dest.Longitude, 4326))
AS distance
--INTO #includeDistances
FROM #orig dest
答案 1 :(得分:37)
以下函数以英里为单位给出两个地理坐标之间的 距离
create function [dbo].[fnCalcDistanceMiles] (@Lat1 decimal(8,4), @Long1 decimal(8,4), @Lat2 decimal(8,4), @Long2 decimal(8,4))
returns decimal (8,4) as
begin
declare @d decimal(28,10)
-- Convert to radians
set @Lat1 = @Lat1 / 57.2958
set @Long1 = @Long1 / 57.2958
set @Lat2 = @Lat2 / 57.2958
set @Long2 = @Long2 / 57.2958
-- Calc distance
set @d = (Sin(@Lat1) * Sin(@Lat2)) + (Cos(@Lat1) * Cos(@Lat2) * Cos(@Long2 - @Long1))
-- Convert to miles
if @d <> 0
begin
set @d = 3958.75 * Atan(Sqrt(1 - power(@d, 2)) / @d);
end
return @d
end
以下函数以千米为单位给出两个地理坐标之间的 距离
CREATE FUNCTION dbo.fnCalcDistanceKM(@lat1 FLOAT, @lat2 FLOAT, @lon1 FLOAT, @lon2 FLOAT)
RETURNS FLOAT
AS
BEGIN
RETURN ACOS(SIN(PI()*@lat1/180.0)*SIN(PI()*@lat2/180.0)+COS(PI()*@lat1/180.0)*COS(PI()*@lat2/180.0)*COS(PI()*@lon2/180.0-PI()*@lon1/180.0))*6371
END
以下函数以千米为单位给出两个地理坐标之间的 距离 使用sql server 2008中引入的 Geography 数据类型
DECLARE @g geography;
DECLARE @h geography;
SET @g = geography::STGeomFromText('LINESTRING(-122.360 47.656, -122.343 47.656)', 4326);
SET @h = geography::STGeomFromText('POINT(-122.34900 47.65100)', 4326);
SELECT @g.STDistance(@h);
<强>用法:强>
select [dbo].[fnCalcDistanceKM](13.077085,80.262675,13.065701,80.258916)
答案 2 :(得分:5)
看起来微软侵入了所有其他受访者的大脑,并使他们尽可能地编写复杂的解决方案。 这是没有任何附加函数/声明语句的最简单方法:
SELECT geography::Point(LATITUDE_1, LONGITUDE_1, 4326).STDistance(geography::Point(LATITUDE_2, LONGITUDE_2, 4326))
只需替换您的数据,而不是LATITUDE_1
,LONGITUDE_1
,LATITUDE_2
,LONGITUDE_2
,例如:
SELECT geography::Point(53.429108, -2.500953, 4326).STDistance(geography::Point(c.Latitude, c.Longitude, 4326))
from coordinates c
答案 3 :(得分:4)
当您使用SQL 2008或更高版本时,我建议您查看GEOGRAPHY数据类型。 SQL内置了对地理空间查询的支持。
e.g。你在GEOGRAPHY类型的表格中有一个列,其中将填充坐标的地理空间表示(请查看上面链接的MSDN参考示例)。此数据类型然后公开允许您执行一系列地理空间查询的方法(例如,找到2点之间的距离)
答案 4 :(得分:4)
Create Function [dbo].[DistanceKM]
(
@Lat1 Float(18),
@Lat2 Float(18),
@Long1 Float(18),
@Long2 Float(18)
)
Returns Float(18)
AS
Begin
Declare @R Float(8);
Declare @dLat Float(18);
Declare @dLon Float(18);
Declare @a Float(18);
Declare @c Float(18);
Declare @d Float(18);
Set @R = 6367.45
--Miles 3956.55
--Kilometers 6367.45
--Feet 20890584
--Meters 6367450
Set @dLat = Radians(@lat2 - @lat1);
Set @dLon = Radians(@long2 - @long1);
Set @a = Sin(@dLat / 2)
* Sin(@dLat / 2)
+ Cos(Radians(@lat1))
* Cos(Radians(@lat2))
* Sin(@dLon / 2)
* Sin(@dLon / 2);
Set @c = 2 * Asin(Min(Sqrt(@a)));
Set @d = @R * @c;
Return @d;
End
GO
用法:
选择dbo.DistanceKM(37.848832506474,37.848732506474,27.83935546875,27.83905546875)
输出:
0,02849639
您可以使用已注释的浮点数更改@R参数。
答案 5 :(得分:1)
除了之前的答案,这里有一种计算SELECT内部距离的方法:
CREATE FUNCTION Get_Distance
(
@La1 float , @Lo1 float , @La2 float, @Lo2 float
)
RETURNS TABLE
AS
RETURN
-- Distance in Meters
SELECT GEOGRAPHY::Point(@La1, @Lo1, 4326).STDistance(GEOGRAPHY::Point(@La2, @Lo2, 4326))
AS Distance
GO
用法:
select Distance
from Place P1,
Place P2,
outer apply dbo.Get_Distance(P1.latitude, P1.longitude, P2.latitude, P2.longitude)
标量函数也可以工作,但在计算大量数据时效率非常低。
我希望这可以帮助别人。