Le portail francophone de la géomatique

Bonjour,

Merci pour ta réponse.

Note finale est de type string à la base

str(region@data$note_finale)
 chr [1:22] NA NA NA NA "5" NA "5" "4" NA "5" "5" ...

Du coup j'ai réussi à faire quelque chose en bidouillant (je passe par des break) mais je ne trouve pas ça propre et ça ne me satisfait pas tout à fait: comment faire si j'ai de vrais attributs qualitatif de type "bon","moyen",... sans avoir a transformer ça en valeurs numériques. Ce devrait être faisable simplement?

L'autre problème c'est qu'il y a bien le bon nombre de régions colorées, les bonnes notes...
Mais...-alors que ma jointure entre mon shapefile des régions et mes données fonctionne bien (pas de décalage)- les notes sur la carte finale ne sont pas du tout attribuées aux bonnes régions et je me demande si ce n'est pas du aux N.A (que maptools ne prendrais pas en compte provoquant un décalage)

Elles devraient être attribuées aux régions suivantes sur la carte BOURGOGNE, CENTRE, CHAMPAGNE-ARDENNE, FRANCHE-COMTE, HAUTE-NORMANDIE, ILE-DE-FRANCE, NORD-PAS-DE-CALAIS, PICARDIE et si on regarde la carte ce n'est pas du tout le cas...

[img]http://imageshack.com/a/img910/7830/jUqEIB.jpg[/img]

region <-readShapePoly("F:\\regions_metropole.shp" )

requete<- "SELECT * from qualification_repartition_rarete_all where cd_ref_demande='79734' order by nom_region asc"
colnames(region@data)

espece <- dbGetQuery(con, requete);espece 
 
nom_espece <- "Acer campestre L., 1753"

region@data <- merge(region@data,espece, by.x="NOM_REG",by.y="nom_region", all.x=TRUE) ;region2

pal<- c('#F41322','#F49E13','#E5D310','#F4F009', '#2DF70E','#0A9B0F')

region@data$test <- as.character(cut(as.numeric(region@data$note_finale),
breaks = c(0,1,2,3,4,5,100),
labels = pal,
include.lowest = TRUE,
right = FALSE))

legende2<-c('0','1','2','3','4','5')

x<-plot(region, col=region@data$test, border = "dark grey"  )

title(main=nom_espece) 
legend("bottomleft",legend =legende2,fill=pal,border = "dark grey",title = "Note finale",bty = "n")

Voilà. Je reste dans l'attente de bonnes idées

Petit ajout par rapport au message précédent : on peut bien sûr définir directement des vecteurs de valeurs dans la fonction mapvalues, sans passer par la création d'un tableau qui fasse office de dictionnaire.

region$testcol <- mapvalues(region$test, from = c(1, 2, 3), to = c("red", "green", "blue"))

Merci, ça a super bien fonctionné!

Une autre petite question (je tente ma chance).
Je cherche à creer des cartes en boucle et j'uilise normalement

liste<-list("espece1","espece2","espece3",...)
for (i in liste) {...}

la liste qui en résulte est de ce type

[[1]]
[1] "espece1"

[[2]]
[1] "espece2"

[[3]]
[1] "espece3"

ce qui me permet de faire varier le nom de l'espèce dans ma requête sql qui récupère les notes pour une espèce.

Comme j'ai plusieurs centaines d'espèces, plutôt que de mettre en dur la liste dans le code je cherchais à la générer à partir d'un vecteur (un champ d'une table sous postgres.

J'ai tenté le as.liste(vecteur) comme indiqué ici

http://stackoverflow.com/questions/1636 … ent-in-the

(j'ai vérifié l'exemple fonctionne bien)

requete<- "SELECT distinct espece from matable"
code_espece <- dbGetQuery(con, requete)

as.list(code_espece)

mais cela me renvoi

[1] "espece1"    "espece2"      "espece3" 

et une erreur a l'execution de ma boucle de type:

Could not create executeSELECT * from matable2 where nom_espece=c("espece1", "espece2", "espece3" [... truncated]

Je ne comprend pas que l'expression as.list fonctionne avec un vecteur de type c("1","2","3") et pas avec un vecteur issu d'une table de données: pourquoi ne lit-il pas chaque ligne comme un élément de liste?

Merci d'avance pour l'aide que vous pourrez m'apporter

Bonjour,

re-moi....

Savez vous comment on peut afficher un dataframe directement en dessous d'une carte?

Ne trouvant pas de doc a ce sujet j'ai tenté de modifier le:

text(x=coordx,y=coordy,labels="TEXTE", font=1, cex=0.5)

en faisant un:

text(x=coordx,y=coordy,labels=tableau.espece, font=1, cex=0.5)

Mais cela ne fonctionne pas.

Merci d'avance