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 dans le package map
- france.png (28.77 KiB) Viewed 3561 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 par representation avec la fonction plot
- france-polygon.png (29.62 KiB) Viewed 3561 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 depuis la Db de RGeostats
- france-db.png (70.32 KiB) Viewed 3561 times
La difference entre les differentes representations graphiques est liee aux types de projections utilisees.
Merci de me dire ce qui ne fonctionne pas.