Plotten einer Formdatei mit dem Fehler ggplot2

8

In den letzten Tagen habe ich den Kopf für die Verwendung von R für die Kartierung ziemlich stark verändert. Ich habe R ausgiebig zum Modellieren usw. benutzt, aber nicht diese Art von Arbeit vorher. Ich habe einige Fragen und Probleme bezüglich Shapefiles, wie sie gelesen werden und so weiter.

Ich habe die Shape-Dateien vom Australian Bureau of heruntergeladen Statistiken gibt es zahlreiche Dateien mit Staatsgrenzen, Postleitzahlen, Städten und so weiter. Die Formdateien sind massiv, die australischen Staatsgrenzen haben etwa 1,8 Millionen Koordinatenpunkte, die andere Datei, die ich ausprobiert habe, war der statistische Bereich, der über 8 Millionen enthält. Ich habe mit dieser Datei nichts gemacht, da sie einfach zu groß für meine R-Einrichtung ist.

Ich habe die Shape-Datei mit readShapePoly eingelesen und so konvertiert

%Vor%

Nachdem ich die State-Borders-Shape-Datei von SpatialPolygonsDataFrame in einen normalen Daten-Frame konvertiert hatte, zeichnete ich sie erfolgreich auf, aber es dauerte ewig und das Detail war zu groß. Ich dachte, thinnedSpatialPoly zu verwenden, um es zu vereinfachen, aber es gibt den Fehler:

%Vor%

Welcher Google kann mir nicht helfen.

Meine nächste Strategie bestand darin, sie in SAS einzulesen und proc greduce zu verwenden, die die Datei aufnimmt und ein Dichtefeld erstellt, und Sie können wählen, wie dicht die Polygone sind.

%Vor%

SAS hat Mist-Grafiken und konnte das Ding nicht einmal für mich grafisch darstellen, also exportierte ich das Dataset und lese es in R mit dem neuen Dichtefeld, das ich hoffte, den Dataframe zu unterteilen und in meinen Plots zu verwenden.

Mein Problem ist jetzt, dass wenn ich in R in Plot gehe, ich diese

bekomme %Vor%

Ich vermute, es liegt daran, dass die Polygone nicht in Ordnung sind oder gebrochen sind? Ich habe diese Funktion verwendet, um zu versuchen, meinen Polygonen wieder beizutreten, aber immer noch kein Glück

%Vor%

Also, endlich meine Fragen! Wie gehen Menschen mit großen übermäßig detaillierten Formdateien um? Warum wurde das SpatialPoly-Working nicht ausgedünnt (ich möchte SAS möglichst vermeiden)? Wie kann ich meine Handlung nicht wie Mist aussehen lassen?

Endlich meine R-Spezifikationen:

%Vor%     
user1414259 11.07.2012, 08:04
quelle

1 Antwort

8

Wenn du zuerst mit dem Kopf tauchst, gehe nicht in das flache Ende.

Mein ziemlich alter PC kann die digitalen Grenzen im Shapefile-Format lesen, kein Problem:

%Vor%

aber die Karte ist deutlich zu detailliert. Ich habe es auch in Quantum GIS geladen, damit ich einen guten alten Zoom haben und herumschwenken kann, und jede kleine Insel ist da. Ich denke, es ist eines der detailliertesten Länderkarten, die ich je gesehen habe. Zweitens sollten Sie versuchen, eine leicht vereinfachte Karte der Staaten zu finden (siehe www.gadm.org für mögliche Fälle).

Also schauen wir ob gSimplify aus dem Paket: rgeos hilft:

%Vor%

das entfernt viele der kleinen Inseln, aber leider für mich (und ein großer Teil der Bevölkerung) entfernt es auch New South Wales. Nicht gut. Wenn ich die Toleranz schließlich herabsetze, kann ich etwas bekommen, das NSW hält:

%Vor%

Aber es gibt eindeutig ein Problem mit gSimplify oder den Shapefile-Daten selbst. Wie auch immer, wenn ich aus2 in ein Shapefile zurückspare, gibt es eine massive Verkleinerung der Größe, die .shp ist 180k statt 29 Megabyte.

Ich bleibe auch beim Plotten mit Basisgrafiken.

    
Spacedman 11.07.2012, 09:39
quelle

Tags und Links