Le Blog SQL Server d'ElSüket » Systèmes d'information géographiques

Convertir des coordonnées sexagésimales en degrés et en radians

elsuket — Tue, 27 Oct 2009 17:14:11 +0000

Souvent, les coordonnées géographiques sont exprimées en degrés, minutes, secondes.
Mais nos systèmes sont bien plus à l’aise avec des nombres décimaux, et il est donc nécessaire de convertir des coordonnées sexagésimales en degrés ou bien en radians.
Voici une petite fonction pour le faire …

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

-------------------------------
-- Nicolas SOUQUET - 28/10/2009
-------------------------------
ALTER FUNCTION Fn_CoordonneesSexagesimalesToDegresRadians
(
@degresLat TINYINT,
@minutesLat TINYINT,
@secondesLat TINYINT,
--------------------
@degresLong TINYINT,
@minutesLong TINYINT,
@secondesLong TINYINT
)
RETURNS TABLE
AS
RETURN
(
SELECT latitudeDeg,
longitudeDeg,
RADIANS(latitudeDeg) AS latitudeRad,
RADIANS(longitudeDeg) AS longitudeRad
FROM
(
SELECT @degresLat + @minutesLat / 60.0 + @secondesLat / 3600.0 AS latitudeDeg,
@degresLong + @minutesLong / 60.0 + @secondesLong / 3600.0 AS longitudeDeg
) AS TMP
)

ElSuket

Vérifier l’unicité d’une position avec le type GEOGRAPHY sous SQL Server 2008

elsuket — Thu, 29 Oct 2009 22:33:34 +0000

SQL Server 2008 a introduit de nombreux nouveaux types de données, dont le type de données géographiques GEOGRAPHY.
Ce n’est pas un type habituel, puisque c’est un type CLR.NET intégré à SQL Server.

L’avantage présenté par l’intégration de ce type est l’ensemble des méthodes standard « livrées » avec ce type, qui permettent d’extraire très simplement une latitude (attribut Lat), une longitude (attribut Long), ou encore de connaître la distance entre deux points géographiques avec la méthode STDistance.

Mais il devient alors plus complexe de garantir l’unicité de positions dans une table où l’on stocke la celle de plusieurs villes.
En effet, si une colonne de ce type est spécifiée comme clé d’une contrainte d’unicité, le moteur de base de données SQL Server lève une exception.

Est-il possible de contourner ce problème ?

Si nous tentons de créer la table suivante :

1
2
3
4
5
6

CREATE TABLE dbo.TbCoordonneeVille
(
IDCoordonneeVille INT NOT NULL IDENTITY CONSTRAINT PK_TbCoordonneeVille PRIMARY KEY,
nomVille VARCHAR(45) NOT NULL CONSTRAINT UQ_TbCoordonneeVille_nomVille UNIQUE,
geoPositionVille GEOGRAPHY NOT NULL CONSTRAINT UQ_TbCoordonneeVille_geoPositionVille UNIQUE
)

Nous obtenons l’erreur :

Msg 1919, Niveau 16, État 1, Ligne 1
La colonne ‘geoPositionVille’ dans la table ‘TbCoordonneeVille’ n’est pas d’un type valide lui permettant d’être utilisée en tant que colonne clé dans un index.
Msg 1750, Niveau 16, État 0, Ligne 1
Impossible de créer la contrainte. Voir les erreurs précédentes.

Nous sommes donc contraints de créer la table sans la contrainte d’unicité sur la colonne geoPositionVille :

1
2
3
4
5
6

Comment faire pour garantir l’unicité des positions des villes sans perdre les avantages offerts par le type GEOGRAPHY ?
Il est possible de créer une contrainte d’unicité sur des colonnes calculées persistantes.

Ajoutons donc deux colonnes calculées, une pour la latitude, et une autre pour la longitude.
Attention, la casse sur les méthodes CLR est importante.

1
2
3
4
5
6

-------------------------------
-- Nicolas SOUQUET - 29/10/2009
-------------------------------
ALTER TABLE dbo.TbCoordonneeVille
ADD latitudeVille AS (geoPositionVille.Lat) PERSISTED NOT NULL,
longitudeVille AS (geoPositionVille.Long) PERSISTED NOT NULL

Il nous suffit ensuite de créer une contrainte d’unicité sur ces deux colonnes :

1
2
3
4
5
6

-------------------------------
-- Nicolas SOUQUET - 29/10/2009
-------------------------------
ALTER TABLE dbo.TbCoordonneeVille
ADD CONSTRAINT UQ_TbCoordonneeVille_latitudeVille_longitudeVille
UNIQUE (latitudeVille, longitudeVille)

Essayons maintenant d’insérer deux villes dont le nom est différent, mais dont les coordonnées sont identiques :

1
2
3
4
5
6
7

INSERT INTO dbo.TbCoordonneeVille (nomVille, geoPositionVille)
VALUES ('Toulouse', GEOGRAPHY::STPointFromText('POINT(43.604503 1.444026)', 4326))
GO

INSERT INTO dbo.TbCoordonneeVille (nomVille, geoPositionVille)
VALUES ('Lyon', GEOGRAPHY::STPointFromText('POINT(43.604503 1.444026)', 4326))
GO

La première instruction se passe sans problèmes.
En revanche la seconde échoue, et nous obtenons l’erreur suivante :

Msg 2627, Niveau 14, État 1, Ligne 1
Violation de la contrainte UNIQUE KEY ‘UQ_TbCoordonneeVille_latitudeVille_longitudeVille’. Impossible d’insérer une clé en double dans l’objet ‘dbo.TbCoordonneeVille’.
L’instruction a été arrêtée.

Vérifions le contenu de la table :

Une solution alternative est de créer une vue indexée, ce qui présentera un bon avantage si la table est souvent lue.
Supprimons la table puis recréons-la :

1
2
3
4
5
6
7
8
9

DROP TABLE dbo.TbCoordonneeVille
GO

CREATE TABLE dbo.TbCoordonneeVille
(
IDCoordonneeVille INT NOT NULL IDENTITY CONSTRAINT PK_TbCoordonneeVille PRIMARY KEY,
nomVille VARCHAR(45) NOT NULL CONSTRAINT UQ_TbCoordonneeVille_nomVille UNIQUE,
geoPositionVille GEOGRAPHY NOT NULL
)

Créons ensuite la vue :

1
2
3
4
5
6
7
8
9
10

-------------------------------
-- Nicolas SOUQUET - 29/10/2009
-------------------------------
CREATE VIEW dbo.VwCoordonneeVille
WITH SCHEMABINDING
AS
SELECT nomVille,
geoPositionVille.Lat AS latitudeVille,
geoPositionVille.Long AS longitudeVille
FROM dbo.TbCoordonneeVille

Puis indexons-la :

1
2
3
4
5

-------------------------------
-- Nicolas SOUQUET - 29/10/2009
-------------------------------
CREATE UNIQUE CLUSTERED INDEX IXUQC_VwCoordonneeVille_latitudeVille_longitudeVille
ON dbo.VwCoordonneeVille(latitudeVille, longitudeVille)

Et tentons les mêmes insertions que celles que nous avons faites pour tester la contrainte d’unicité sur les colonnes calculées persistantes :

1
2
3
4
5
6
7

Là encore, la première instruction se passe sans problème, mais la seconde échoue, provoquant l’erreur suivante :

Msg 2601, Niveau 14, État 1, Ligne 1
Impossible d’insérer une ligne de clé en double dans l’objet ‘dbo.VwCoordonneeVille’ avec un index unique ‘IXUQC_VwCoordonneeVille_latitudeVille_longitudeVille’.
L’instruction a été arrêtée.

Interrogeons néanmoins la vue :

ElSuket

Obtenir une représentation sexagésimale d’une coordonnée géographique en degrés

elsuket — Thu, 19 Nov 2009 18:19:49 +0000

Faisant suite à ce billet, voici une petite fonction pour représenter sexagésimalement une coordonnée géographique exprimée en degrés :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

-------------------------------
-- Nicolas SOUQUET - 19/11/2009
-------------------------------
ALTER FUNCTION Fn_SIG_CoordonneesDegresToSexagesimales
(@_degresDecimaux FLOAT)
RETURNS VARCHAR(16)
WITH SCHEMABINDING
AS
BEGIN
-- Extraction des degrés
DECLARE @nDegresSexagesimaux TINYINT
SELECT @nDegresSexagesimaux = CAST(@_degresDecimaux AS TINYINT)

-- Extraction des minutes
DECLARE @fMinutesSexagesimaux FLOAT
SELECT @fMinutesSexagesimaux = (@_degresDecimaux - @nDegresSexagesimaux) * 60
DECLARE @nMinutesSexagesimaux TINYINT
SELECT @nMinutesSexagesimaux = CAST(@fMinutesSexagesimaux AS TINYINT)

-- Extraction des secondes
DECLARE @fSecondesSexagesimaux FLOAT
SELECT @fSecondesSexagesimaux = (@fMinutesSexagesimaux - @nMinutesSexagesimaux) * 60

-- Resultat
RETURN CAST(@nDegresSexagesimaux AS VARCHAR(2)) + '° ' +
CAST(@nMinutesSexagesimaux AS VARCHAR(2)) + '′ ' +
CAST(@fSecondesSexagesimaux AS VARCHAR(14)) + '"'
END

Ainsi l’exécution de :

SELECT dbo.Fn_SIG_CoordonneesDegresToSexagesimales (47.51027)

Nous retourne 47° 30′ 36.972″

ElSuket

Vérifier qu’un point est contenu dans un polygone de type GEOGRAPHY

elsuket — Thu, 10 Dec 2009 14:29:56 +0000

Voici comment vérifier qu’un point de type GEOGRAPHY, c’est à dire une position définie par sa latitude et sa longitude, appartient à la surface définie par un polygone, toujours sous le type GEOGRAPHY

Comme je suis toulousain, j’ai pris pour exemple des coordonnées de positions aux 4 coins de la Place du Capitole pour définir le polygone que voici :

1
2
3
4
5
6
7
8
9
10

DECLARE @gPlaceCapitole GEOGRAPHY

SET @gPlaceCapitole = GEOGRAPHY::STGeomFromText
(
'POLYGON
(
(43.604832 1.442787,43.60492 1.443754,43.603922 1.443987,43.603823 1.443028,43.604832 1.442787)
)',
4326
);

Il est nécessaire que :

– le premier point soit la même que la dernier point pour qu’on dispose effectivement d’un polygone,
– les points soient ordonnés dans le sens inverse des aiguilles d’une montre,
– les coordonnées soit données suivant un couple {latitude, longitude}

Prenons maintenant un point sur la Croix du Languedoc qui se situe sur la Place du Capitole, et un point dans l’Hôtel de Ville, qui jouxte la Place du Capitole :

1
2
3
4
5

DECLARE @pCroixLanguedoc GEOGRAPHY
SET @pCroixLanguedoc = GEOGRAPHY::Point(1.443304, 43.604386, 4326)

DECLARE @pHotelDeVille GEOGRAPHY
SET @pHotelDeVille = GEOGRAPHY::Point(43.604526,1.444323, 4326)

Ici, il est nécessaire les coordonnées soient données suivant un couple {longitude, latitude}.
Faites-y attention, on se fait facilement avoir

Vérifions maintenant que la Croix du Languedoc est bien sur la Place du Capitole :

SELECT @g.STIntersects(@pCroixLanguedoc )

L’exécution de ce code nous retourne un BIT ayant pour valeur 1.

Qu’en est-il de l’Hôtel de Ville ?

SELECT @g.STIntersects(@pHotelDeVille)

Ici nous obtenons 0.

Bons SIGs !

ElSuket.