Billet

Ah… 36 000 communes feront toujours la différence

Billet publié le 11/12/2011

Je viens de découvrir que l’IGN mettait maintenant gratuitement à disposition de tous le fichier shapefile GEOFLA communes, qui permet de produire des cartes de toutes les communes de France. Combinez cette carte avec des données en provenance de l’INSEE (par exemple sur la répartition des revenus), et hop :

Sur cette carte, plus c’est bleu, plus, en gros, “c’est riche”, et plus c’est rouge, plus, en gros, “c’est pauvre”. [l’indicateur est la médiane du revenu fiscal des ménages divisé par le nombre d’unité de consommation du ménage, en 2009]
Cette carte a été réalisé avec R, assez simplement : voir mon tutoriel sur la cartographie avec R pour en savoir plus, et en faire autant.

26 commentaires

Un commentaire par RjG (11/12/2011 à 22:22)

Tu mérites pleinement ton titre de géographe d’honneur, Baptiste!

Un commentaire par Ced (11/12/2011 à 22:22)

Comme on pourrait s’y attendre, les zones littorales présentent toutes une bande de revenus élevés, particulièrement au niveau de l’Europe rhénane. Cependant on observe un exception notable sur la façade méditerranéenne. Pourquoi cette exception ?

Un commentaire par Joël (11/12/2011 à 22:23)

J’ai assez peu utilisé les données de la DGI pour le moment, mais suffisamment pour m’interroger sur ce qu’elles nous disent exactement. Certes, c’est un indicateur de richesse, mais je me demande aussi dans quelle mesure ce n’est pas un indicateur d’ingéalités (à plus forte raison quand on prend la moyenne plutôt que la médiane)… Qu’en penses-tu ?
Autre question : comment as-tu fait avec les communes concernées par le secret statistique ?

Un commentaire par Baptiste Coulmont (11/12/2011 à 22:29)

> joel. je n’ai pas cherché loin. J’ai juste pris le fichier que j’ai trouve sur le site de l’insee. Le but etait de faire rapidement une carte…
> ced. Peut etre parce qu’il y a plus de retraités sur la cote mediterranéenne.

Un commentaire par N. Holzschuch (11/12/2011 à 23:07)

On voit vraiment bien la structure des villes, notamment en province : centre-ville pauvre, péri-urbain riche, et campagne pauvre à nouveau. Chaque petit rond bleu au centre de la mer de rouge correspond à une grande ville.

On remarque aussi bien “l’effet vignoble” : cote d’or, bordelais, beaujolais, Alsace et champagne…

Il ne reste que le gros rond bleu autour de Toulouse que je n’arrive pas à m’expliquer.

Un commentaire par Joël (11/12/2011 à 23:41)

@ced il y a bien une bande bleue sur le littoral méditerranéen de Menton jusqu’aux Saintes-Maries de la Mer, non comprises. Ca n’est que du côté des Pyrénées-Orientales et de l’Aude qu’on observe du rouge – on n’a pas dans ces départements l’attraction de cadres sup et de riches retraités qu’on a en PACA…

Un commentaire par Pierre M (12/12/2011 à 2:00)

J’aime.
Intéressant de comparer la couleur des iles. Entre Ouessant et Ré… par exemple.

Un commentaire par ageekculture (12/12/2011 à 20:02)

Sympa,
Il serait intéressant d’avoir les cartes électorale (rouge = gauche, bleu = droite) des communes et de les comparer ou de comparer cette carte avec d’autres indicateurs comme la taille de la population, l’age moyen etc… et de trouver un moyen sympa de visualiser les relations.

Un commentaire par Baptiste Coulmont (12/12/2011 à 20:46)

Concernant le vote, les données doivent se trouver quelque part, mais je ne sais pas trop où. Concernant l’âge moyen et autres données, on trouve beaucoup de choses sur le site du Recensement (INSEE)
http://www.recensement.insee.fr/home.action

Un commentaire par Tal Galili (13/12/2011 à 13:10)

Hello there, what a beautiful post!
May I invite you to join:
http://www.r-bloggers.com/lang/-/french
By submitting your blog here:
http://www.r-bloggers.com/lang/add-your-blog

With regards,
Tal

Un commentaire par DM (13/12/2011 à 14:22)

Félicitations!

Un commentaire par Create maps with maptools R package | Statisfaction (13/12/2011 à 20:27)

[…] Coulmont explains on his blog how to use the R package maptools. It is based on shapefile files, for example the ones offered by […]

Un commentaire par Joël (14/12/2011 à 20:32)

Pour les données électorales, c’était un peu compliqué de les obtenir jusque il y a peu. Elles sont désormais (une partie au moins) sur data.gouv.fr.

Un commentaire par Baptiste Coulmont (14/12/2011 à 20:47)

Et on peut faire de jolies cartes, comme cette carte des votes pour “Le Pen” au premier tour des présidentielles de 2007 : https://picasaweb.google.com/lh/photo/_GwJy_Ubcgrksp0dByd-wMa_X9CNRZARz6gXztIpSWQ?feat=directlink

Un commentaire par 36 000 communes, oui, mais donnons-leur un poids – Joël Gombin (15/12/2011 à 7:12)

[…] Coulmont a publié il y a quelques jours un billet dans lequel il présente une carte représentant, pour chaque commune, la médiane du revenu fiscal des ménages divisé par le nombre […]

Un commentaire par Joël (15/12/2011 à 7:13)

Je me suis permis de proposer une autre version de cette carte, en forme de cartogramme : http://www.joelgombin.fr/?p=577

Un commentaire par Alain (28/12/2011 à 2:20)

Erreur importation shapefile GEOFLA communes

Bonjour Mr Coulmont,

Suite à la lecture de vos documents (de très bonne qualité), j’ai tenté de créer une carte des communes avec R. Lors de l’importation du fichier GEOFLA des communes (téléchargé sur l’IGN version shapefile), j’obtiens une erreur dans le code R :

Erreur dans validityMethod(as(object, superClass)) :
Geographical CRS given to non-conformant data: 1242375 7110524

J’obtiens également la même erreur avec le shapefile des départements. Comment avez-vous importé dans R le fichier des communes ? Avez-vous eu la même erreur ? Existe t-il un fichier Rdata ?

J’utilise R2.13.0 64 bits sous Windows Seven.
Ci joint une partie du code utilisé :

library(ggplot2)
library(sp)
library(maptools)
library(gpclib)
gpclibPermit()

france<-readShapeSpatial("C:/Users/Afo/Downloads/GEOFLA/COMMUNES/COMMUNE.shp",proj4string=CRS"+proj=longlat"))
plot(france)

Merci de votre aide.

Cordialement,

AF.

Un commentaire par Baptiste Coulmont (28/12/2011 à 18:08)

Je pense que :
readShapeSpatial(“C:/Users/Afo/Downloads/GEOFLA/COMMUNES/COMMUNE.shp”)
suffit
car la projection utilisée par “geofla communes” est la projection “Lambert”, pas long/lat… ceci dit, je n’ai toujours pas compris comment fonctionnent les projections dans les SHP…

Un commentaire par Didi (06/02/2012 à 16:52)

Bonjour,

vous ne voulez pas publier votre code ?

On va pas le revendre sur ebay.

Un commentaire par Baptiste Coulmont (06/02/2012 à 16:59)

> Didi : mon “tutoriel” donne de nombreuses indications. J’aurai pu aussi vous l’envoyer, si vous m’aviez donné une véritable adresse e-mail.

Un commentaire par mathieu rajerison (12/03/2012 à 19:13)

Bonjour Mr Coulmont,

Suite au lien de Mr Gombin sur une représentation en anamorphose, je me permets de vous proposer ici une autre vision, cette fois, sous forme de paysage d’archipels de prospérité:
http://datagistips.blogspot.com/2012/03/archipels-de-prosperite.html

Tant qu’on y est… :)

Un commentaire par Les blogs de Data Geek : Mathieu Rajerison, la passion de la GéoData et des représentations | Data Publica (11/07/2012 à 13:52)

[…] mon blog. Cette carte faisait suite à une série de visualisations proposées par des bloggeurs: Baptiste Coulmont avec une choroplèthe, Joël Gombin avec une […]

Un commentaire par Alexandre (25/09/2012 à 11:56)

Bonjour,

Merci pour ce blog et toutes ces informations (assez rares sur maptools)
Effectivement, le code du tutoriel ne fonctionne pas avec chez moi non plus avec les fichiers Geofla

le code :
france<-readShapeSpatial("departement.shp",proj4string=CRS("+proj=longlat"))

renvoie l'erreur :
Geographical CRS given to non-conformant data: 1242375 7110524

Dès lors, il est possible par la suite d'afficher la carte mais pas de placer des points

Si vous avez trouvé une solution depuis, je suis serais curieux de la connaitre

Bien à vous,

Merci encore pour la diffusion de tous ces conseils

Un commentaire par Baptiste Coulmont (25/09/2012 à 13:18)

A mon avis, c’est parce que vous avez téléchargé un shapefile qui est dans une autre projection (Lambert).
Il y a, sur le site de l’IGN, 2 shapefiles avec des projections différentes, il faut télécharger le bon (celui qui est en Lon / lat)

Un commentaire par Maps with R, and polygon boundaries | Freakonometrics (12/02/2013 à 2:10)

[…] to draw maps. Let us start with something simple, like French regions. Baptiste mentioned on his blog that shapefiles can be downloaded from http://ign.fr/ website. Hence, if you extract […]

Un commentaire par Arnaud (03/10/2013 à 18:03)

Bonjour, je sais que ce post est un peu ancien, mais je suis très intéressé par le code. En effet, j’ai bien compris la maîtrise des cartes sous R (au passage un très grand merci pour votre tutorial), mais je bute, concernant les communes sur des questions d’affichage (limites des communes et autres).
Merci d’avance