Plotten Weltkarte in orthographische Projektion gibt "nicht endliche Punkte"

Question

Ich habe ein Shapefile von Ländern der Welt, heruntergeladen von here. Ich kann es in R plotten

countries <- readOGR("shp","TM_WORLD_BORDERS-0.3",encoding="UTF-8",stringsAsFactors=F) 
par(mar=c(0,0,0,0),bg=rgb(0.3,0.4,1)) 
plot(countries,col=rgb(1,0.8,0.4)) 

mit

Jetzt möchte ich es in Orthogonalprojektion plotten (Erde aus dem Weltraum zu sehen ist), so versuche ich,

countries <- spTransform(countries,CRS("+proj=ortho +lat_0=-10 +lon_0=-60")) 

ich auch mit dem x_0 und Y_0 gespielt Parameter (wie here angegeben), aber ich habe immer den Fehler:

non finite transformation detected: 
[1] 45.08332 39.76804  Inf  Inf 
Erro em .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, : 
    failure in Polygons 3 Polygon 1 points 1 
Além disso: Mensagens de aviso perdidas: 
In .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, : 
    108 projected point(s) not finite 

manchmal im 3. Polygon, manchmal im 7.. Woher kommen diese "Inf"? Ich muss einen Parameter ändern? Ich möchte, dass die Karte wie diese

aber über Südamerika zentriert plotten. Danke für Ihre Hilfe!

Answer 1

Geben Sie das maps Paket einen Versuch. Es gibt eine Warnung über die Punkte, die nicht projiziert werden können, aber es gibt keinen Fehler und schließt den Prozess. Mit ein wenig Gefummel, nämlich die Füllfarbe für den Ozean Einstellung (this Antwort geholfen, dieses Problem zu lösen), konnte ich die Karte, die Sie mit ein paar Zeilen an emulieren:

library(maps) 

## start plot & extract coordinates from orthographic map 
o <- c(-10,-60,0) # oreantation 
xy <- map("world",proj="orthographic", orientation=o, bg="black") 
xy <- na.omit(data.frame(do.call(cbind, xy[c("x","y")]))) 

## draw a circle around the points for coloring the ocean 
polygon(max(xy$x)*sin(seq(0,2*pi,length.out=100)),max(xy$y)*cos(seq(0,2*pi,length.out=100)), 
     col="blue4", border=rgb(1,1,1,0.5), lwd=2) 

## overlay world map 
colRamp <- colorRampPalette(c("lemonchiffon", "orangered")) 
map("world",proj="orthographic", orientation=o, 
    fill=TRUE, col=colRamp(5), add=TRUE)