Carte Mali Afrique occidentale

[aliou balde]

Bonjour,

Je débute sur R, et j'aimerais faire une carte du Mali avec les subdivisions administratives mais R me propose une carte vide et j'ai téléchargé (URL ci-dessous) un fichier shapefile et des fichiers rds mais j'ai du mal à les exploiter sous R . Pouvez-vous m'aider?
http://www.gadm.org/country
http://www.diva-gis.org/datadown

Code : Tout sélectionner

shapefile <- read.shapefile ("C:/Users/aliou/Downloads/mali_administrative/mali_administrative")


xx0 <- readRDS("C:/Users/aliou/Downloads/GADM/MLI_adm0.rds")
xx1 <- readRDS("C:/Users/aliou/Downloads/GADM/MLI_adm1.rds")
xx2 <- readRDS("C:/Users/aliou/Downloads/GADM/MLI_adm2.rds")
xx3 <- readRDS("C:/Users/aliou/Downloads/GADM/MLI_adm3.rds")
xx4 <- readRDS("C:/Users/aliou/Downloads/GADM/MLI_adm4.rds")

Merci
Amitiés

[dicko ahmadou]

Salut,

Votre question n'est pas très claire, qu'entendez vous par "J'ai du mal à les exploiter" ou par carte vide ?

Quel est le contenu de:

Code : Tout sélectionner

str(xx1)

Pouvez vous dire simplement ce qui ne marche pas pour que quelqu'un puisse vous aider.

Merci

[aliou balde]

Salut,

Code : Tout sélectionner

mapImage <- map(database="world",region="Mali")

quand j'utilise ce code il me renvoi l'image suivante (carte du Mali sans subdivisions administratives)
ensuite j'utilise le code suivant :

Code : Tout sélectionner

area <- readShapePoly("C:/Users/aliou/Documents/MLI_adm/MLI_adm1.shp")
area.points <- fortify(area)

  ggmap(mapImage) +
  geom_polygon(aes(x = long,
                   y = lat,
                   group = group),
               data = area.points,
               color = colors[9],
               fill = colors[6],
               alpha = 0.5) +
  labs(x = "Longitude",
       y = "Latitude")

qui ne marche pas
Voici un aperçu du contenu xx1 est :

Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
..@ data :'data.frame': 9 obs. of 13 variables:
.. ..$ OBJECTID : int [1:9] 1 2 3 4 5 6 7 8 9
.. ..$ ID_0 : int [1:9] 138 138 138 138 138 138 138 138 138
.. ..$ ISO : chr [1:9] "MLI" "MLI" "MLI" "MLI" ...
.. ..$ NAME_0 : chr [1:9] "Mali" "Mali" "Mali" "Mali" ...
.. ..$ ID_1 : int [1:9] 1 2 3 4 5 6 7 8 9
.. ..$ NAME_1 : chr [1:9] "Bamako" "Gao" "Kayes" "Kidal" ...
.. ..$ HASC_1 : chr [1:9] "ML.BA" "ML.GA" "ML.KY" "ML.KD" ...
.. ..$ CCN_1 : int [1:9] NA NA NA NA NA NA NA NA NA
.. ..$ CCA_1 : chr [1:9] "" "" "" "" ...
.. ..$ TYPE_1 : chr [1:9] "District" "Région" "Région" "Région" ...
.. ..$ ENGTYPE_1: chr [1:9] "District" "Region" "Region" "Region" ...
.. ..$ NL_NAME_1: chr [1:9] "" "" "" "" ...
.. ..$ VARNAME_1: chr [1:9] "" "" "" "" ...
..@ polygons :List of 9
.. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
.. .. .. ..@ Polygons :List of 1
.. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
.. .. .. .. .. .. ..@ labpt : num [1:2] -7.98 12.61
.. .. .. .. .. .. ..@ area : num 0.0203
.. .. .. .. .. .. ..@ hole : logi FALSE
.. .. .. .. .. .. ..@ ringDir: int 1
.. .. .. .. .. .. ..@ coords : num [1:182, 1:2] -7.93 -7.93 -7.93 -7.93 -7.93 ...
.. .. .. ..@ plotOrder: int 1
.. .. .. ..@ labpt : num [1:2] -7.98 12.61
.. .. .. ..@ ID : chr "1"
.. .. .. ..@ area : num 0.0203
.. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
.. .. .. ..@ Polygons :List of 1
.. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
.. .. .. .. .. .. ..@ labpt : num [1:2] 1.35 16.79
.. .. .. .. .. .. ..@ area : num 14.9
.. .. .. .. .. .. ..@ hole : logi FALSE
.. .. .. .. .. .. ..@ ringDir: int 1
.. .. .. .. .. .. ..@ coords : num [1:930, 1:2] 0.444 0.44 0.431 0.426 0.419 ...
.. .. .. ..@ plotOrder: int 1
.. .. .. ..@ labpt : num [1:2] 1.35 16.79
.. .. .. ..@ ID : chr "2"
.. .. .. ..@ area : num 14.9
.. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots

Merci

[dicko ahmadou]

Je ne comprends pas le problème car l'objet xx1 est bien un SpatialPolygon donc tu peux faire:

Code : Tout sélectionner

## Faire la carte
plot(xx1)

A mon humble avis, il faut d'abord revoir les bases du logiciel R et comprendre comment on manipule les objets.
Sinon ce code marche plutôt bien, tu peux t'en inspirer pour faire marcher le script que tu utilises.

Code : Tout sélectionner

library(raster)
library(ggplot2)

tmp <- tempdir()
mli_reg <- getData("GADM", country = "MLI", level = 1, path = tmp)
plot(mli_reg)

###
mli_reg_df <- fortify(mli_reg)

ggplot(mli_reg_df, aes(long, lat, group = group)) +
  geom_polygon(fill = "transparent", colour = "black") +
  coord_equal()

Du courage

[aliou balde]

Merci

Je débute vraiment et je me retrouve à faire une carte donc j'ai beaucoup à apprendre j'avoue

Merci encore

[Pierre-Yves Berrard]

Si c'est juste une carte, pourquoi ne pas utiliser un logiciel comme QGis ?

[aliou balde]

Je ne connais pas QGis, mais la finalité c'est d'arriver à avoir une carte avec les différentes aires de santé par régions voire cercles (département) et le nombre d'individus suivis par aire de santé par exemple