Traçage de Cartes Thématiques dans R en Utilisant les fichiers de formes et de Fichiers de Données à partir de Différentes Sources
Donné un shapefile, comment la forme et l'usage d'un fichier de données afin d'être en mesure de tracer des cartes thématiques à l'aide des identifiants qui correspondent à la forme des régions dans le shapefile?
#Download English Government Office Network Regions (GOR) from:
#http://www.sharegeo.ac.uk/handle/10672/50
tmp_dir = tempdir()
url_data = "http://www.sharegeo.ac.uk/download/10672/50/English%20Government%20Office%20Network%20Regions%20(GOR).zip"
zip_file = sprintf("%s/shpfile.zip", tmp_dir)
download.file(url_data, zip_file)
unzip(zip_file, exdir = tmp_dir)
library(maptools)
#Load in the data file (could this be done from the downloaded zip file directly?
gor=readShapeSpatial(sprintf('%s/Regions.shp', tmp_dir))
#I can plot the shapefile okay...
plot(gor)
#and I can use these commands to get a feel for the data...
summary(gor)
attributes(gor@data)
gor@data$NAME
#[1] North East North West
#[3] Greater London Authority West Midlands
#[5] Yorkshire and The Humber South West
#[7] East Midlands South East
#[9] East of England
#9 Levels: East Midlands East of England ... Yorkshire and The Humber
#download data from http://www.justice.gov.uk/downloads/publications/statistics-and-data/courts-and-sentencing/csq-q3-2011-insolvency-tables.csv
#insolvency<- read.csv("~/Downloads/csq-q3-2011-insolvency-tables.csv")
insolvency=read.csv("http://www.justice.gov.uk/downloads/publications/statistics-and-data/courts-and-sentencing/csq-q3-2011-insolvency-tables.csv")
insolvencygor.2011Q3=subset(insolvency,Time.Period=='2011 Q3' & Geography.Type=='Government office region')
#tidy the data
require(gdata)
insolvencygor.2011Q3=drop.levels(insolvencygor.2011Q3)
names(insolvencygor.2011Q3)
#[1] "Time.Period" "Geography"
#[3] "Geography.Type" "Company.Winding.up.Petition"
#[5] "Creditors.Petition" "Debtors.Petition"
levels(insolvencygor.2011Q3$Geography)
#[1] "East" "East Midlands"
#[3] "London" "North East"
#[5] "North West" "South East"
#[7] "South West" "Wales"
#[9] "West Midlands" "Yorkshire and the Humber"
#So what next? Ayant eu jusque-là, comment dois-je prendre la prochaine étape pour générer une thématique/choropleth carte, que les couleurs de chaque région en fonction de l'Débiteurs.Pétition de la valeur, par exemple?
(J'ai également remarqué une possible gotcha - il y a une incohérence dans la capitalisation GOR niveaux: "Yorkshire et Humber" et "Yorkshire et Humber" )
Il semble que ce DONC post (le premier sur la "Relative" de la liste sur le côté droit de la barre en) pourrait aider à obtenir le reste de la voie.
J'ai vu ça, mais il a manqué le correctif le premier couple de fois rond... Alors ce que je dois faire quelque chose comme: gor@data=merge(insolvencygor.2011T3,gor@data,par.x='Géographie',par.y='NOM') plot(gor,col=niveaux(du gr@data$Créanciers.Pétition)) mais avec une bonne couleur de la cartographie, et le bon nom de la région de cartes (je remarque que ce n'est pas seulement Yorks et Humberside que incompatibilité...)
Pourriez-vous réduire votre jeu de données pour un petit exemple qui illustre vous le problème? Vous pourriez économiser en utilisant la commande enregistrer et télécharger soit, ou à un serveur et poste le lien ici. Ce serait faire de la réduction de la grande quantité de code que vous avez beaucoup plus facile.
Je ne sais pas comment faire le shapefile les petits... je suis en train de créer le poste de: a) désireux de produire des cartes thématiques, b) sachant que les fichiers de formes existent pour générer des contours sur les cartes, et c) err, c'est ça... C'est peut-être aussi important de dire que la raison pour laquelle je suis à l'aide de R, dans ce cas, c'est parce qu'il me semble que R peut gérer les fichiers de formes...
Il y a d'autres outils qui peuvent gérer les fichiers de formes assez facilement, par exemple, QGis ou GVSig, qui sont à la fois open source, logiciels de SIG.
J'ai vu ça, mais il a manqué le correctif le premier couple de fois rond... Alors ce que je dois faire quelque chose comme: gor@data=merge(insolvencygor.2011T3,gor@data,par.x='Géographie',par.y='NOM') plot(gor,col=niveaux(du gr@data$Créanciers.Pétition)) mais avec une bonne couleur de la cartographie, et le bon nom de la région de cartes (je remarque que ce n'est pas seulement Yorks et Humberside que incompatibilité...)
Pourriez-vous réduire votre jeu de données pour un petit exemple qui illustre vous le problème? Vous pourriez économiser en utilisant la commande enregistrer et télécharger soit, ou à un serveur et poste le lien ici. Ce serait faire de la réduction de la grande quantité de code que vous avez beaucoup plus facile.
Je ne sais pas comment faire le shapefile les petits... je suis en train de créer le poste de: a) désireux de produire des cartes thématiques, b) sachant que les fichiers de formes existent pour générer des contours sur les cartes, et c) err, c'est ça... C'est peut-être aussi important de dire que la raison pour laquelle je suis à l'aide de R, dans ce cas, c'est parce qu'il me semble que R peut gérer les fichiers de formes...
Il y a d'autres outils qui peuvent gérer les fichiers de formes assez facilement, par exemple, QGis ou GVSig, qui sont à la fois open source, logiciels de SIG.
OriginalL'auteur psychemedia | 2012-01-12
Vous devez vous connecter pour publier un commentaire.
N'ayant pas vu la forêt pour les arbres, pour répondre à ma propre question, voici un moyen (code suivant à partir de code dans la question):
De sorte que cette approche n'est fusionner les données numériques dans le shapefile, puis tracer directement.
Cela dit, n'est-ce pas le moyen le plus propre serait de conserver le fichier de données et le fichier de formes distinctes? (Je ne suis toujours pas sûr de savoir comment le faire?)
Et je partage votre sentiment le code pourrait être plus court et plus au point, mais sans reproductible exemple, il est difficile à démontrer.
Reproductible dans quel sens? L'ajout de la réponse ci-dessus le code pour le code en question génère la carte thématique, certes le téléchargement étapes sont nécessaires comme un manuel d'intervention, plutôt que d'être manipulé par le script.
Je veux dire, pas reproductible dans le fait que je ne peut pas le coller dans R et obtenir ce que vous avez. Juste essayez de le rendre aussi facile que possible pour les gens à s'impliquer.
Bon, le point de prise de vue - mais, comme suggéré dans les commentaires à la question, je ne sais pas comment faire pour charger un fichier de forme à partir du fichier zip, il est inclus dans le shapegeo site web?;-)
OriginalL'auteur psychemedia