vendredi 27 novembre 2015

Rentrer des points manuellement dans QGIS

(ou presque manuellement)

L'idée est : comment faire, à partir des coordonnées connues d'un point (par exemple, un point GPS noté manuellement et dont on ne possèderait pas le fichier .gpx), pour l'afficher dans QGIS ?

Il suffit de créer un fichier texte ou Coma Separated Value (*.txt ou *.csv) contenant les coordonnées du (ou des) point(s), comme par exemple :

N 43° 18' 53.1" ; E 0° 21' 56.6388"

Le contenu minimum de ce fichier est le suivant :

X_COORD;Y_COORD
-0 21 56.6388;43 18 53.1

Remarques :
Les coordonnées X sont les coordonnées de longitude : d'Ouest en Est, et les coordonnées Y sont les coordonnées de latitude : du Nord au Sud.

Les valeurs de X à l'Est du méridien de Greenwich (et sans doute aussi celles de Y dans l'hémisphère sud) doivent être notées en négatif !

Les symboles des degrés (º), minutes (') et secondes ('') doivent être omis et remplacés par des espaces.

Il est tout à fait possible de rajouter autant de colonnes que souhaitées pour stocker toutes les informations nécessaires.

Les colonnes peuvent être séparées par les symboles suivants : virgule (,) - tabulation (\t) - espace ( ) - colonne (|) - point virgule (;)

Si le symbole virgule (,) est utilisé comme séparateur de colonnes, il faut utiliser un point (.) comme séparateur décimal. Si la virgule (,) est utilisé comme séparateur décimal, il faut l'indiquer à QGIS en cochant la case correspondante.

Les intitulés des colonnes sont complètement libres !

Pour l'importer on aura recours au bouton "Virgule" (Ajouter une couche de texte délimité).

Par exemple le fichier contenant les lignes suivantes :

Projet;site;pt_GPS;X_Coord;Y_Coord
Test GPS;Univ. Pau;001;-0 21 56.6388;43 18 53.1

Peut être importé comme suit :

Une dernière fenêtre demandera quel est le système de coordonnées de références (SCR) à utiliser pour projeter le point.
S'il s'agit d'un point GPS, quel que soit le mode d'affichage degré, minute, seconde (DMS) retenu, on utilisera le SCR WGS 84 (nº de code EPSG : 4326).

Attention : si le SCR est mal choisi, on peut arriver à des erreurs d'affichage de l'ordre de 50 à 100 m d'écart !

Quelques infos pour aller plus loin :

Les modes d'affichages DMS peuvent être les suivants :
degrés décimaux [D,DDº]
degrés minutes décimales [Dº M,MM']
degrés minutes secondes décimales [Dº M' S,SS''] 

NB : si on enlève les décimales on perd en précision.

Rappels de conversion :

Un degré (º) est égal à 60 minutes (') et à 3600 seconds (") :

1º = 60'

1º = 3600"

1' = 60'' 

1' = (1/60)º = 0.01666667º

1'' = (1/60)' = 0.01666667'

1" = (1/3600)º = 2.77778e-4º = 0.000277778º

On note D la partie entière des degrés et DD les degrés décimaux.

Par exemple : 43,314750º est une valeur en degrés décimaux (DD), et 43º est une valeur en degrés entiers (D). "43" est la part entière de la valeur en degrés décimaux.

On converti de la manière suivante :

D,DDº -> Dº M' S''

DD = 43,314750
D = part entière de DD = 43

MM = (DD - D) × 60
MM = (43,314750 - 43) × 60
MM = 0,314750 × 60
MM = 18,885
M = 18

SS = (MM - M) × 60
SS = (18,885 - 18) × 60
SS = (0,885) × 60
SS = 53,1
S = 53

Donc 43,314750º = 43º 18' 53,1'' ~ 43º 18' 53''

Dº M' S'' -> D,DDº

DD = D + M ÷ 60 + S ÷ 3600

DD = 43º + (18' ÷ 60)º + (53,1'' ÷ 3600)º
DD = 43º + 0,3º + 0,01475º
DD = 43,31475º


Source : 
http://www.rapidtables.com/convert/number/degrees-minutes-seconds-to-degrees.htm