[SOLVED] Unités des cordonnées spatiales

If you have problems using RGeostats package or if you simply need some help on an unlisted subject.

[SOLVED] Unités des cordonnées spatiales

Postby ngop » Thu Mar 28, 2019 4:07 pm

Bonjour,

Je ne comprends pas pourquoi quand j'affiche la carte de mes données, les données sont mal réparties (je devrais voir se dessiner la carte de France).

Mes coordonnées spatiales sont de la forme :
Long Lat
4,948042804 46,31373926
4,997248819 46,01967662
4,948480086 45,92523506

En faisant ma db.create j'ai pourtant bien défini les locators avec la longitude et la latitude. Avec Qgis je n'ai pas de problème de représentation (chargé en tant que EPSG 4326).

Merci de votre aide
ngop
 
Posts: 4
Joined: Thu Mar 28, 2019 4:02 pm

Re: Unités des cordonnées spatiales

Postby Didier Renard » Sun Apr 07, 2019 8:33 pm

Bonjour
Je ne comprends pas exactement le probleme rencontre. J'ai essaye de le mimer en faisant tourner le code suivant:

J'utilise la librairie maptools afin de visualiser les cotes francaises:

Code: Select all
library(maps)
library(maptools)
map("france")

france.png
France dans le package map
france.png (28.77 KiB) Viewed 3558 times


Je recupere les coordonnees des traits de cotes (et les limites de departements) de la France utilisant la fonction map mais avec l'option plot=FALSE

Code: Select all
x = map('france',plot=FALSE)$x
y = map('france',plot=FALSE)$y


Je peux verifier que x et y sont des vecteurs de taille 13466: les valeurs NA correspondent a un "lever de plume" dans la representation graphique du polygone. Cette option est compatible avec la fonction generique de trace plot:

Code: Select all
plot(x,y,type="l")

france-polygon.png
France par representation avec la fonction plot
france-polygon.png (29.62 KiB) Viewed 3558 times


Je peux maintenant creer la Db dans RGeostats. Le chargement direct de x dans une variable nommee x1 et de y dans la variable x2 assure que ces variables seront automatiquement attachees aux locators des coordonnees

Code: Select all
library(RGeostats)
db = db.create(x1=x,x2=y)
print(db)


On peut verifier que les coordonnees de points sont bien definis par les longitudes et les latitudes: la longitude varie de -5.14 a 9.56 alors que la latitude varie de 41.37 a 51.09.

Code: Select all
print(db,flag.stats=T,names="x*")


Il suffit maintenant d'en faire la representation graphique:
Code: Select all
plot(db,cex=0.1)

france-db.png
France depuis la Db de RGeostats
france-db.png (70.32 KiB) Viewed 3558 times


La difference entre les differentes representations graphiques est liee aux types de projections utilisees.

Merci de me dire ce qui ne fonctionne pas.
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm

Re: Unités des cordonnées spatiales

Postby ngop » Thu Apr 11, 2019 4:58 pm

Merci beaucoup pour ces explications.

Mon problème provenait du fait que mes coordonnées étaient au format factor et non numeric d'où problème à l'affichage.

Je vais réutiliser le fond de carte des départements que vous avez présenté.
ngop
 
Posts: 4
Joined: Thu Mar 28, 2019 4:02 pm

Re: Unités des cordonnées spatiales

Postby Fabien Ors » Fri Apr 12, 2019 8:08 am

Merci pour votre message.
Cette information pourra servir à bien d'autres utilisateurs de RGeostats !
Fabien Ors
Administrateur du site
 
Posts: 226
Joined: Thu Sep 20, 2012 1:07 pm


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 6 guests

cron