Plotten SHP-Datei R

Question

Ich versuche zu lesen und zeichnen Sie eine SHP-Datei aus Argentinien, um eine benutzerdefinierte Choropleth zu erstellen. Ich folge genau allen Codebeispielen, um diese Karte zu laden und zu plotten, aber nichts funktioniert.

Das rote Zeichen bleibt eingeschaltet, da R Informationen verarbeitet, aber nie etwas anzeigt oder nicht mehr funktioniert.

Dies ist mein Bildschirm, wenn ich auf "Plot" laufen und dann sagt Rstudio funktioniert nicht mehr. Irgendeine Idee, warum es nicht funktionieren könnte?

Answer 1

Zuerst ggplot2 scheint in der Lage zu sein, dies zu handhaben ziemlich gut:

library(rgdal) 
library(ggplot2) 

arg_shp <- readOGR("../data/arg_adm/ARG_adm1.shp", "ARG_adm1") 

arg_shp_map <- fortify(arg_shp) 

gg <- ggplot() 
gg <- gg + geom_map(data=arg_shp_map, map=arg_shp_map, 
        aes(long, lat, map_id=id), 
        color="#2b2b2b", size=0.15, fill=NA) 
gg <- gg + coord_map() 
gg <- gg + ggthemes::theme_map() 
gg 

Da Sie eine GADM Datei verwenden, brauchen Sie nicht wirklich eine externe Shape-Datei tho verwenden verwenden, nur getData():

library(maptools) 
library(raster) 

arg_adm <- getData('GADM', country='ARG', level=1) 
arg_adm_map <- fortify(arg_adm) 

(gleichen ggplot2 Code von oben funktioniert, nur aus dem Variablennamen sub).

Ich kann nicht das Basis-Plotting-System, um dies (entweder auf einem Grafikgerät) zu plotten. Ich denke, es gibt einige komplexe Polygone in einem Teil der Küste. Wir können damit umgehen ziemlich leicht tho:

library(maptools) 
library(raster) 
library(rgeos) 
library(magrittr) 
library(sp) 

arg_adm <- getData('GADM', country='ARG', level=1) 

gSimplify(arg_adm, 0.01, topologyPreserve=TRUE) %>% 
    SpatialPolygonsDataFrame(dat=arg_adm@data) -> new_arg 

plot(new_arg, lwd=0.25) 

(Grundstück Standardränder und die Art und Weise es mit Projektionen in dieser unterschiedlich großen Karte zur Folge hat)

Sofern Sie die Grenze brauchen Präzise, ​​vereinfachte Shapefiles können viel Zeit beim Plotten sparen. Sie können die optimierte Form in ggplot2 als auch verwenden, zusammen mit einem Argentinien zugeschnittene Kartenprojektion:

new_arg_map <- fortify(new_arg) 

gg <- ggplot() 
gg <- gg + geom_map(data=new_arg_map, map=new_arg_map, 
        aes(long, lat, map_id=id), 
        color="#2b2b2b", size=0.15, fill=NA) 
gg <- gg + ggalt::coord_proj("+proj=aeqd +lat_0=-37.869859624840764 +lon_0=-66.533203125") 
gg <- gg + ggthemes::theme_map() 
gg 

(Obwohl, coord_map() ‚s Standard-Mercator-Projektion pbly gd genug für die meisten Anwendungen ist)

Answer 2

Soweit ich weiß, es ist posible ein shp mit MapTools zu lesen, aber ihre Unterstützung ist begrenzt. Haben Sie versucht, es mit rgdal::readOGR(dir, layer) zu lesen?

Das Environment-Fenster zeigt, dass Sie ein "Large" -Objekt haben, aber wie groß ist es? wenn die Formen die nationalen oder staatlichen Grenzen sind, sollte es kein Problem geben, aber wenn Sie alle Straßen oder alle Blöcke des Landes in der Form haben, ändern sich die Dinge.

Versuchen Sie es mit einem Komponententest: Versuchen Sie, die internationale Grenze oder die Staatsgrenzen zu plotten. Wenn Sie das tun können, wäre das Problem die Größe der Formen.

PD: reebot Ihren Computer und versuchen Sie es erneut, es könnte helfen, wenn Rstudio es das einzige laufende Programm ist.