Pokud vše, co chcete dělat s SqlGeography, je sledovat body a využívat prostorové indexy SQL Server 2008, můžete, jak poznamenali jiní, skrýt svůj sloupec prostorových dat z Linq do SQL a používat UDF nebo uložené procedury. Předpokládejme, že máte tabulku AddressFields, která obsahuje pole Zeměpisná šířka a Zeměpisná délka. Přidejte tuto tabulku do svého souboru DBML a napište libovolný požadovaný kód, který nastaví pole Zeměpisná šířka a Zeměpisná délka. Potom kód SQL níže přidá do této tabulky pole Geo geogarphy a vytvoří spouštěč v databázi, který automaticky nastaví pole Geo na základě polí Zeměpisná šířka a Zeměpisná délka. Mezitím níže uvedený kód také vytváří další užitečné UDF a uložené procedury:DistanceBetween2 (Již jsem měl DistanceBetween) vrací vzdálenost mezi adresou reprezentovanou v AddressField a zadaným párem zeměpisná šířka/délka; DistanceWithin vrací různá pole ze všech polí AddressFields v rámci zadané mílové vzdálenosti; UDFDistanceWithin dělá totéž jako uživatelsky definovaná funkce (užitečné, pokud ji chcete vložit do většího dotazu); a UDFNearestNeighbors vrátí pole z AddressField odpovídající zadanému počtu sousedů nejblíže konkrétnímu bodu. (Jedním z důvodů pro použití UDFNearestNeighbors je to, že SQL Server 2008 nebude optimalizovat využití prostorového indexu, pokud zavoláte objednávku pouze voláním DistanceBetween2.)
Budete to muset upravit změnou AddressFields na vaši tabulku a přizpůsobením polí z této tabulky, která chcete vrátit (podívejte se na kód kolem odkazů na AddressFieldID). Poté to můžete spustit ve vaší databázi a zkopírovat výsledné uložené procedury a UDF do vašeho DBML a poté je můžete použít v dotazech. Celkově vám to umožňuje poměrně snadno využívat prostorový index bodů.
-----------------------------------------------------------------------------------------
--[1]
--INITIAL AUDIT
select * from dbo.AddressFields
GO
--ADD COLUMN GEO
IF EXISTS (SELECT name FROM sysindexes WHERE name = 'SIndx_AddressFields_geo')
DROP INDEX SIndx_AddressFields_geo ON AddressFields
GO
IF EXISTS (SELECT b.name FROM sysobjects a, syscolumns b
WHERE a.id = b.id and a.name = 'AddressFields' and b.name ='Geo' and a.type ='U' )
ALTER TABLE AddressFields DROP COLUMN Geo
GO
alter table AddressFields add Geo geography
--[2]
--SET GEO VALUE
GO
UPDATE AddressFields
SET Geo = geography::STPointFromText('POINT(' + CAST([Longitude] AS VARCHAR(20)) + ' ' +
CAST([Latitude] AS VARCHAR(20)) + ')', 4326)
--[3] VYTVOŘIT INDEX
IF EXISTS (SELECT name FROM sysindexes WHERE name = 'SIndx_AddressFields_geo')
DROP INDEX SIndx_AddressFields_geo ON AddressFields
GO
CREATE SPATIAL INDEX SIndx_AddressFields_geo
ON AddressFields(geo)
--UPDATE STATS
UPDATE STATISTICS AddressFields
--AUDIT
GO
select * from dbo.AddressFields
--[4] POSTUP VYTVOŘENÍ USP_SET_GEO_VALUE PARA 1 LATITUDE 2 LOGITUDE
IF EXISTS (SELECT name FROM sysobjects WHERE name = 'USPSetGEOValue' AND type = 'P')
DROP PROC USPSetGEOValue
GO
GO
CREATE PROC USPSetGEOValue @latitude decimal(18,8), @longitude decimal(18,8)
AS
UPDATE AddressFields
SET Geo = geography::STPointFromText('POINT(' + CAST(@longitude AS VARCHAR(20)) + ' ' +
CAST(@latitude AS VARCHAR(20)) + ')', 4326)
WHERE [Longitude] [email protected] and [Latitude] = @latitude
GO
--TEST
EXEC USPSetGEOValue 38.87350500,-76.97627500
GO
--[5] VYTVOŘIT SPUŠTĚNÍ PŘI ZMĚNĚ/VLOŽENÍ HODNOTY ŠÍŘKY/DLOUHÉ ---> NASTAVIT GEOKÓD
IF EXISTS (SELECT name FROM sysobjects WHERE name = 'TRGSetGEOCode' AND type = 'TR')
DROP TRIGGER TRGSetGEOCode
GO
CREATE TRIGGER TRGSetGEOCode
ON AddressFields
AFTER INSERT,UPDATE
AS
DECLARE @latitude decimal(18,8), @longitude decimal(18,8)
IF ( UPDATE (Latitude) OR UPDATE (Longitude) )
BEGIN
SELECT @latitude = latitude ,@longitude = longitude from inserted
UPDATE AddressFields
SET Geo = geography::STPointFromText('POINT(' + CAST(@longitude AS VARCHAR(20)) + ' ' +
CAST(@latitude AS VARCHAR(20)) + ')', 4326)
WHERE [Longitude] [email protected] and [Latitude] = @latitude
END
ELSE
BEGIN
SELECT @latitude = latitude ,@longitude = longitude from inserted
UPDATE AddressFields
SET Geo = geography::STPointFromText('POINT(' + CAST(@longitude AS VARCHAR(20)) + ' ' +
CAST(@latitude AS VARCHAR(20)) + ')', 4326)
WHERE [Longitude] [email protected] and [Latitude] = @latitude
END
GO
--[6] CREATE PROC USP_SET_GEO_VALUE_INITIAL_LOAD ----> POUZE JEDNO SPUŠTĚNÍ
IF EXISTS (SELECT name FROM sysobjects WHERE name = 'USPSetAllGeo' AND type = 'P')
DROP PROC USPSetAllGeo
GO
CREATE PROC USPSetAllGeo
AS
UPDATE AddressFields
SET Geo = geography::STPointFromText('POINT(' + CAST([Longitude] AS VARCHAR(20)) + ' ' +
CAST([Latitude] AS VARCHAR(20)) + ')', 4326)
GO
--[7] STÁVAJÍCÍ PROC DistanceBetween, který vrací vzdálenost mezi dvěma zadanými body
--podle dvojice souřadnic zeměpisné šířky/délky. --ALTER PROC VzdálenostMezi2
IF EXISTS (SELECT name FROM sysobjects WHERE name = 'DistanceBetween2' AND type = 'FN')
DROP FUNCTION DistanceBetween2
GO
CREATE FUNCTION [dbo].[DistanceBetween2]
(@AddressFieldID as int, @Lat1 as real,@Long1 as real)
RETURNS real
AS
BEGIN
DECLARE @KMperNM float = 1.0/1.852;
DECLARE @nwi geography =(select geo from addressfields where AddressFieldID = @AddressFieldID)
DECLARE @edi geography = geography::STPointFromText('POINT(' + CAST(@Long1 AS VARCHAR(20)) + ' ' +
CAST(@Lat1 AS VARCHAR(20)) + ')', 4326)
DECLARE @dDistance as real = (SELECT (@nwi.STDistance(@edi)/1000.0) * @KMperNM)
return (@dDistance);
END
GO--TEST
Vzdálenost mezi 2 12159,40,75889600,-73,99228900
--[8] VYTVOŘIT POSTUP USPDistanceWithin
-- VRÁTÍ SEZNAM ADRES Z tabulky AddressFields
POKUD EXISTUJE (VYBERTE název Z sysobjects WHERE name ='USPDistanceWithin' AND type ='P') PROCEDURA DROP USPDistanceWithin
GO
CREATE PROCEDURE [dbo].USPDistanceWithin
(@lat as real,@long as real, @distance as float)
AS
BEGIN
DECLARE @edi geography = geography::STPointFromText('POINT(' + CAST(@Long AS VARCHAR(20)) + ' ' +
CAST(@Lat AS VARCHAR(20)) + ')', 4326)
SET @distance = @distance * 1609.344 -- convert distance into meter
select
AddressFieldID
,FieldID
,AddressString
,Latitude
,Longitude
,LastGeocode
,Status
--,Geo
from
AddressFields a WITH(INDEX(SIndx_AddressFields_geo))
where
a.geo.STDistance(@edi) < = @Distance
END
GO
--TEST
--do 3 milUSPDistanceWithin 38,90606200,-76,92943500,3GO--do 5 milUSPDistanceWithin 38,90606200,-76,92943500,5GO--do 10 mil,29-02943500,01W,029-043PDistanceW.01
--[9] VYTVOŘIT FUNKCI FNDistanceWithin
-- VRÁTÍ SEZNAM ADRES Z tabulky AddressFields
POKUD EXISTUJE (VYBERTE název Z sysobjects WHERE name ='UDFDistanceWithin' AND type ='TF') DROP FUNCTION UDFDistanceWithin
GO
CREATE FUNCTION UDFDistanceWithin
(@lat as real,@long as real, @distance as real)
RETURNS @AddressIdsToReturn TABLE
(
AddressFieldID INT
,FieldID INT
)
AS
BEGIN
DECLARE @edi geography = geography::STPointFromText('POINT(' + CAST(@Long AS VARCHAR(20)) + ' ' +
CAST(@Lat AS VARCHAR(20)) + ')', 4326)
SET @distance = @distance * 1609.344 -- convert distance into meter
INSERT INTO @AddressIdsToReturn
select
AddressFieldID
,FieldID
from
AddressFields a WITH(INDEX(SIndx_AddressFields_geo))
where
a.geo.STDistance(@edi) < = @Distance
RETURN
END
GO
--TEST
--within 3 milesselect * from UDFDistanceWithin(38.90606200,-76.92943500,3)GO--within 5 milesselect * from UDFDistanceWithin( 38.90606200,-76.92943500,5)U70W6th-DistanceWithin from 0.0.
--[9] VYTVOŘIT FUNKCI UDFNearestNeighbors
-- VRÁTÍ SEZNAM ADRES Z tabulky AddressFields
POKUD EXISTUJE (VYBERTE název Z sysobjects WHERE name ='UDFNearestNeighbors' AND type ='TF') FUNKCE DOP FUNCTION UDFNearestNeighbors
GO
POKUD EXISTUJE (VYBERTE název ZE sysobjects WHERE name ='čísla' AND xtype ='u') DROP TABLE numbers
GO
-- First, create a Numbers table that we will use below.
SELECT TOP 100000 IDENTITY(int,1,1) AS n INTO numbers FROM MASTER..spt_values a, MASTER..spt_values b CREATE UNIQUE CLUSTERED INDEX idx_1 ON numbers(n)
GO
CREATE FUNCTION UDFNearestNeighbors
(@lat as real,@long as real, @neighbors as int)
RETURNS @AddressIdsToReturn TABLE
(
AddressFieldID INT
,FieldID INT
)
AS
BEGIN
DECLARE @edi geography = geography::STPointFromText('POINT(' + CAST(@Long AS VARCHAR(20)) + ' ' +
CAST(@Lat AS VARCHAR(20)) + ')', 4326)
DECLARE @start FLOAT = 1000;
WITH NearestPoints AS
(
SELECT TOP(@neighbors) WITH TIES *, AddressFields.geo.STDistance(@edi) AS dist
FROM Numbers JOIN AddressFields WITH(INDEX(SIndx_AddressFields_geo))
ON AddressFields.geo.STDistance(@edi) < @start*POWER(2,Numbers.n)
ORDER BY n
)
INSERT INTO @AddressIdsToReturn
SELECT TOP(@neighbors)
AddressFieldID
,FieldID
FROM NearestPoints
ORDER BY n DESC, dist
RETURN
END
GO
--TEST
--50 sousedů vyberte * z UDFNearestNeighbors(38,90606200,-76,92943500,50)GO--200 sousedů vyberte * z UDFNearestNeighbors( 38,90606200,-76,92943500200)GOP>GOP>