SQL Server：計算Lat/Long的半徑？

Question

說我有一個城市的經緯度，我需要找出距離這個地點100英里範圍內的所有機場。我怎麼做到這一點？我的數據駐留在SQL Server中。 1張桌子有拉特和長的所有城市信息，另一張桌子有拉特和長的機場信息。

Answer 1

我在幾年前使用過/寫過，並且它足夠接近我所需要的。如果我沒有記錯的話，公式的一部分會考慮地球曲率，但它已經有一段時間了。我使用郵政編碼，但您可以輕鬆地適應城市 - 相同的邏輯。

ALTER PROCEDURE [dbo].[sp_StoresByZipArea] (@zip nvarchar(5), @Radius float) AS 

DECLARE @LatRange float 
DECLARE @LongRange float 
DECLARE @LowLatitude float 
DECLARE @HighLatitude float 
DECLARE @LowLongitude float 
DECLARE @HighLongitude float 

DECLARE @istartlat float 
DECLARE @istartlong float 

SELECT @iStartlat=Latitude, @iStartLong=Longitude from zipcodes where zipcode=@ZIP 

SELECT @LatRange = @Radius/((6076/5280) * 60) 
SELECT @LongRange = @Radius/(((cos((@iStartLat * 3.141592653589/180)) * 6076.)/5280.) * 60) 

SELECT @LowLatitude = @istartlat - @LatRange 
SELECT @HighLatitude = @istartlat + @LatRange 
SELECT @LowLongitude = @istartlong - @LongRange 
SELECT @HighLongitude = @istartlong + @LongRange 

/** Now you can create a SQL statement which limits the recordset of cities in this manner: **/ 

SELECT * FROM ZipCodes 
WHERE (Latitude <= @HighLatitude) AND (Latitude >= @LowLatitude) AND (Longitude >= @LowLongitude) AND (Longitude <= @HighLongitude) 

Answer 2

首先...轉換城市的數據點

DECLARE @point geography; 

SELECT geography::STPointFromText('POINT(' + CAST(@lat AS VARCHAR(20)) + ' ' + 
        CAST(@lon AS VARCHAR(20)) + ')', 4326) 

其中@lat和@lon是城市問題的經度和緯度。

然後你就可以查詢該表...

SELECT [column1],[column2],[etc] 
FROM [table] 
WHERE @point.STBuffer(160934.4).STIntersects(geography::STPointFromText(
    'POINT(' + CAST([lat] AS VARCHAR(20)) + ' ' + 
    CAST([lon] AS VARCHAR(20)) + ')', 4326)); 

，其中160934.4是米100英里數。

雖然這會很慢。如果你想做更多的空間工作，你可以添加一個持續的計算列（因爲經緯度點不會真的改變），然後使用空間索引。

ALTER TABLE [table] 
    ADD geo_point AS geography::STPointFromText('POINT(' + CAST([lat] AS VARCHAR(20)) 
    + ' ' + CAST([lon] AS VARCHAR(20)) + ')', 4326) PERSISTED; 

CREATE SPATIAL INDEX spix_table_geopt 
    ON table(geo_point) 
    WITH (BOUNDING_BOX = (0, 0, 500, 200)); --you'd have to know your data