raster und ggplot map nicht ganz in reihenfolge R

Question

Ich versuche ein räumliches Raster mit ggplot2 zu plotten.

require(raster) 
require(ggplot2) 

Daten herunterladen, Last als Raster des raster Paket. Weitere Details zu diesem Datenprodukt finden Sie unter here. Dann wandle das Raster in Punkte um, so dass es gut mit ggplot spielt.

system('wget https://www.dropbox.com/s/7jsdqgc9wjcdznv/NADP_wet_deposition_nh4_0.5x0.5_grid_annual_R1.txt') 
layer<- raster("path/to/raster/NADP_wet_deposition_nh4_0.5x0.5_grid_annual_R1.txt") #you need to specify your own path here, wherever the downloaded file is saved. 
raster.points <- rasterToPoints(layer) 
raster.points <- data.frame(raster.points) 
colnames(raster.points) <-c('x','y','layer') 

Jetzt ggplot2 verwenden, um eine Karte zu machen, und legen Sie das Raster über.

mp <- NULL 
#grab US map and choose colors 
map.US <- borders("usa", colour='white',fill='black', lwd=0.4) 
mp <- ggplot(data=raster.points, aes(y=y, x=x)) 
mp <- mp + map.US 
mp <- mp + geom_raster(aes(fill=layer)) 
mp <- mp + theme(axis.text.y=element_blank(), 
       axis.text.x=element_blank(), 
       axis.title.y=element_blank(), 
       axis.title.x=element_blank(), 
       axis.ticks=element_blank(), 
       panel.background = element_rect(fill='black'), 
       plot.background = element_rect(fill='black'), 
       panel.grid.major=element_blank(), 
       panel.grid.minor=element_blank()) 
mp 

Die Ausgabe sieht wie folgt aus:

Wie Sie sehen können, die Dinge fast Linie, aber nicht ganz. alles ist leicht nach rechts verschoben. Was kann das verursachen und wie kann ich es beheben?

Answer 1

Nach der ORNL-Dokumentation ist die Grenze des Ndep-Diagramms tatsächlich mit den unteren linken Ecken des Gitters ausgerichtet. Um die x- und y-Positionen mit den mittleren Punkten auszurichten (Standard in ggplot), müssen Sie die x-Positionen um 1 Gitterintervall verschieben. Da das Gitterintervall in diesem Fall 0,5 Grad beträgt, subtrahierte ich einen halben Grad von meinem x-Koordinatenvektor.

Die Lösung für dieses Problem wurde von @ 42 in den Kommentaren vorgeschlagen.

So, wie vor, laden Sie die Daten:

system('wget https://www.dropbox.com/s/7jsdqgc9wjcdznv/NADP_wet_deposition_nh4_0.5x0.5_grid_annual_R1.txt') 
layer<- raster("path/to/raster/NADP_wet_deposition_nh4_0.5x0.5_grid_annual_R1.txt") #you need to specify your own path here, wherever the downloaded file is saved. 
raster.points <- rasterToPoints(layer) 
raster.points <- data.frame(raster.points) 
colnames(raster.points) <-c('x','y','layer') 

Und entscheidend ist, einen halben Grad von der x-Vektor der Koordinaten subtrahieren.

raster.points$x <- raster.points$x - 0.5 

Nun gehen Sie vor und Grundstück:

mp <- NULL 
#grab US map and choose colors 
map.US <- borders("usa", colour='white',fill='black', lwd=0.4) 
mp <- ggplot(data=raster.points, aes(y=y, x=x)) 
mp <- mp + map.US 
mp <- mp + geom_raster(aes(fill=layer)) 
mp <- mp + theme(axis.text.y=element_blank(), 
       axis.text.x=element_blank(), 
       axis.title.y=element_blank(), 
       axis.title.x=element_blank(), 
       axis.ticks=element_blank(), 
       panel.background = element_rect(fill='black'), 
       plot.background = element_rect(fill='black'), 
       panel.grid.major=element_blank(), 
       panel.grid.minor=element_blank()) 
mp 

alle aufgereiht!

Answer 2

Basierend auf was 42 und Colin sagen, bietet das ORNL eine Datei, die nicht korrekt ist. Und Sie sollten ihnen wahrscheinlich davon erzählen. Die Datei hat:

xllcorner -124 
yllcorner 25 

Wo es sollte offensichtlich sein:

xllcorner -124.5 
yllcorner 25 

Wenn ja, das bedeutet, dass die Datei zur Zeit der x-Koordinaten als Zellenkanten angibt und die y-Koordinaten als Zellzentren. Das ist in keinster Weise gut. Es gibt Formate, in denen Sie sowohl x als auch y Zellkanten anstelle von Zellzentren darstellen können, aber nicht eine der beiden; In beiden Fällen ist dies im verwendeten Dateiformat (arc-ascii) nicht erlaubt.

Eine gute Möglichkeit, diesen Fehler zu korrigieren, ist shift zu verwenden, nachdem die RasterLayer erstellt:

layer < raster("NADP_wet_deposition_nh4_0.5x0.5_grid_annual_R1.txt") 
layer <- shift(layer, -0.5, 0) 

und dann fortzufahren.Dies ist eine allgemeinere Lösung, da die Rasterdaten für andere Zwecke als ggplotting korrekt verwendet werden können.