如何在R中绘制球体上的点

Question

我正在开展一个项目，我需要在球体上绘制一些点。我所拥有的是每个点的经度和纬度，以及每个点的值。例如，

longitude latitude value
123       23       1.3
75        -34      2.8
190       48       2.1
...      

我想要一个如下的情节。我还希望点的大小与数据集中的值成比例。所以2.8将有一个更大的点，1.3将有一个更小的点，等等。世界地图是可选的。

我想知道R中是否有任何可以完成这项工作的软件包？任何建议或代码都非常感谢！

Answer 1

使用base绘图和包sp，rgdal和maptools的简单示例：

library(sp)
library(maptools)
library(rgdal)
xy <- data.frame(lon=c(-130,110,3,45),lat=c(60,-10,50,30)) #Some coordinates
value <- data.frame(value=c(1.5,0.8,2.3,2)) #Some values for the point size
df <- SpatialPointsDataFrame(xy,value,proj4string=CRS("+proj=lonlat"))
dfMoll <- spTransform(df, CRS("+proj=moll")) #Mollweide projection of the data
data(wrld_simpl) # A base world map
wrld_moll <- spTransform(wrld_simpl, CRS("+proj=moll")) # ... that we projects as well
plot(wrld_moll) #... and plot
points(dfMoll, cex=dfMoll$value, pch=20, col="red") #...with our data points

Answer 2

也许可以通过ggplot2：

来实现

ggplot() +
  geom_polygon(data=world,aes(x=long, y=lat, group=group), fill=NA,colour="black") +
  geom_point(data=d, aes(x=longitude, y=latitude, size=value), color="red") +
  coord_map("mollweide")

这给出了以下地图，不幸的是有一些故障：

顺便说一下，您可以使用coord_map而无需绘制地图：

ggplot() +
  geom_point(data=d, aes(x=longitude, y=latitude, size=value), color="red") +
  coord_map("mollweide")

Answer 3

看看mapproj包。它会为您做各种地图投影。它不会进行绘图，而是为您投射点，然后传递给您选择的函数。