Erste Fehler bei der Verwendung von fortify-Funktion in R (ggplot)

Ich bin versucht, die Funktion zu stärken, die in R zum zeichnen der Karte. Aber ich bin immer ein Fehler während der Verarbeitung der Daten. Ich lese die shape-Datei mit dem folgenden code:

shape <- readShapeSpatial("Lond_City/lond_city.shp")

Der shape-Datei finden Sie auf der folgenden dropbox-Speicherort:
https://www.dropbox.com/sh/d4w6saailydtu1r/5oIa56xV6S

Dann ich versuchte, verwenden Sie den folgenden Befehl ein :

shape1 <- fortify(shape,region="data")

In dem obigen code, was soll ich anstelle von "Daten". Dies ist der Ort, das scheint verwirrend. Ich kann öffnen Sie die Daten oder shapefile in R. Aber, wenn ich das fortify Zeile bekomme ich die folgende Fehlermeldung:

Fehler in '[.Daten.Rahmen'(attr, region) : undefined Spalten ausgewählt

Kann mir jemand sagen, wie man richtig zu bestärken, eine shape-Datei mit Beispiel ? Den Grund möchte ich stärken ist, dass ich möchten, kombinieren Sie die shape-Datei mit den Daten Punkte.

Danke Ihnen so sehr.

Jdbaba

Es sieht aus wie Sie sich mit etwas aus der Ordnance Survey? Wenn Sie einen link zu der shape-Datei, so dass andere es downloaden können, oder zu erklären, genau Woher hast du die Daten aus, dann denke ich, wir können wahrscheinlich helfen.
Slowlearner, ich danke Ihnen für Ihre Antwort. Ich habe den link Hinzugefügt von dem shapefile in die Frage. Bitte haben Sie einen Blick. danke.
Danke für, die. Sie wird fast immer eine schnellere und bessere Reaktion auf StackOverflow, wenn Sie uns die Daten als auch code.

InformationsquelleAutor Jd Baba | 2013-01-17

ggplot2 r

Dein problem ist, dass Sie brauchen, um zu geben fortify eine Spalte, die tatsächlich existiert. Wenn Sie str(lon.shape) von meinem Beispiel-code unten, werden Sie sehen, dass es keine region Spalte:

> str(lon.shape)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame': 1 obs. of  7 variables:
  .. ..$ ons_label: chr "00AA"
  .. ..$ name     : chr "City of London"
  .. ..$ label    : chr "02AA"
  .. ..$ X_max    : num 533843
  .. ..$ y_max    : num 182198
  .. ..$ X_min    : num 0
  .. ..$ y_min    : num 0
  ..@ polygons   :List of 1
  .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
  .. .. .. ..@ Polygons :List of 1
  .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
  .. .. .. .. .. .. ..@ labpt  : num [1:2] 532464 181220
  .. .. .. .. .. .. ..@ area   : num 3151465
  .. .. .. .. .. .. ..@ hole   : logi FALSE
  .. .. .. .. .. .. ..@ ringDir: int 1
  .. .. .. .. .. .. ..@ coords : num [1:771, 1:2] 531027 531029 531036 531074 531107 ...
  .. .. .. ..@ plotOrder: int 1
  .. .. .. ..@ labpt    : num [1:2] 532464 181220
  .. .. .. ..@ ID       : chr "0"
  .. .. .. ..@ area     : num 3151465
  ..@ plotOrder  : int 1
  ..@ bbox       : num [1:2, 1:2] 530967 180404 533843 182198
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "x" "y"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
  .. .. ..@ projargs: chr "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs"

Stattdessen versuchen Sie es mit name wie diese:

lon.df <- fortify(lon.shape, region = "name")

Denen diese Daten erzeugt frame:

> head(lon.df)
      long      lat order  hole piece            group             id
1 531026.9 181611.1     1 FALSE     1 City of London.1 City of London
2 531028.5 181611.2     2 FALSE     1 City of London.1 City of London
3 531036.1 181611.5     3 FALSE     1 City of London.1 City of London
4 531074.0 181610.3     4 FALSE     1 City of London.1 City of London
5 531107.0 181609.3     5 FALSE     1 City of London.1 City of London
6 531117.1 181608.9     6 FALSE     1 City of London.1 City of London

Hier ist ein Weg, um darüber zu gehen, von Anfang bis Ende:

library(rgdal)
library(ggplot2)
library(rgeos)

shape.dir <- "c:/test/london" # use your directory name here

lon.shape <- readOGR(shape.dir, layer = "lond_city")
lon.df <- fortify(lon.shape, region = "name")

ggplot(lon.df, aes(x = long, y = lat, group = group)) +
    geom_polygon(colour = "black", fill = "grey80", size = 1) +
    theme()

...was in den folgenden:

Erste Fehler bei der Verwendung von fortify-Funktion in R (ggplot)

SlowLearner, ich Danke Ihnen so sehr für Ihre wunderbare Lösung. Aber wie es scheint, ich kann nicht laufen, da bekomme ich die folgende Fehlermeldung "isTRUE(gpclibPermitStatus()) ist nicht WAHR. Wissen Sie, wie Sie dieses Problem beheben ? Habe ich versucht zu installieren "gpclibPermit", aber es sagt, dieses Paket ist nicht verfügbar für die version 2.15.2. Irgendwelche Vorschläge ?
Dies ist, weil SlowLearner ist mit readOGR eher als readShapeSpatial. Sie können dieses Problem beheben, durch ausführen gpclibPermit (), aber vorsichtig sein, da dies verboten ist, wenn Sie die software für einen kommerziellen Zweck. Sie können es vorziehen, um dieses Problem zu vermeiden, durch die Annahme der Ansatz in meiner Antwort. Alternativ können Sie die Installation der rgeos - Paket, das keine kommerzielle Einschränkung und sollte das problem auch lösen.
bewahrt die proj4-Informationen, so dass ich in der Regel lieber. Ich bearbeitet meine Antwort ein wenig zu zählen gpclib und die gpclibPermit() nennen.
Wenn Sie library(rgeos) Sie brauchen nicht zu laufen gpclibPermit() als gut, da rgeos ist eine alternative zu gpclib.
Ah, vielen Dank, denke, das ist, warum in einem anderen Stück code ich habe es ist - eine unerklärliche Aufruf rgeos! Bin ich verwirrt, wie, um die Bibliotheken zu liefern, was: sp, map, maptools, gpclib, rgdal, rgeos alle Durcheinander zusammen in meinem Kopf. Wird Bearbeiten.
Danke Jungs für das clearing Dinge. Ich habe eine andere Aufgabe zu tun. Ich muss hinzufügen, die Punkte und die Kontur auf das Grundstück und schneiden Sie die Kontur von der shape-Datei. Geht das in R ? Danke.
Auf StackOverflow, die wir gerne tun", eine Frage auf einmal" so bitte stellen Sie dieser eine separate Frage beinhalten und die Daten (oder fake-Daten, die in etwa den realen Daten).
Ja, es ist möglich, die Suche nach Antworten auf gis.stackexchange.com als auch StackOverflow.

InformationsquelleAutor SlowLearner

2

Sollten Sie label im Ort data.
```
require(maptools); require(ggplot2)
shape <- readShapeSpatial("Lond_City/lond_city.shp")
shape1 <- fortify(shape,region="label")
```
Um das herauszufinden, untersuchte ich die Daten element shape:
```
> shape@data
  ons_label           name label    X_max    y_max X_min y_min
0      00AA City of London  02AA 533842.7 182198.4     0     0
```
Dies zeigt ons_label oder name würde auch gearbeitet haben und möglicherweise besser geeignet für Ihre Zwecke.

Beachten, dass your shape-Datei ist etwas ungewöhnlich, weil es scheint nur eine Gruppe.

InformationsquelleAutor orizon

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.