在正投影中绘制世界地图给出“非有限点”

我有一个世界国家的形状文件，从下载here http://thematicmapping.org/downloads/TM_WORLD_BORDERS-0.3.zip。我可以使用 R 在 R 中绘制它

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))

现在我想用正交投影（从外太空看到的地球）来绘制它，所以我正在尝试

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

我还使用了 x_0 和 y_0 参数（如上所述here http://www.remotesensing.org/geotiff/proj_list/orthographic.html），但我总是收到错误：

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

有时在第三个多边形，有时在第七个多边形。那些“Inf”从哪里来？我需要更改任何参数吗？我想像这样绘制地图

但以南美洲为中心。感谢您的帮助！

给maps打包试试。它会发出有关无法投影的点的警告，但不会给出错误并关闭进程。稍微摆弄一下，即设置海洋的填充颜色（this https://stackoverflow.com/a/8473810/4055438答案帮助解决了这个问题），我能够用几行模拟您附加的地图：

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)