我想将Mollweide投影应用于某些纬度/经度坐标。
set.seed(41)
n <- 200
s <- data.frame(lon = rnorm(n, 0, 180),
lat = rnorm(n, 0, 30))
coordinates(s) = ~lon + lat
projection(s) <- '+proj=longlat'
plot(s)
更改投影会引发错误:
s <- spTransform(s, CRS('+proj=moll'))
non finite transformation detected:
lon lat
Error in spTransform(s, CRS("+proj=moll")) : failure in points
Additionally: Warning:
In spTransform(s, CRS("+proj=moll")) : 1 projected point(s) not finite
我是地理数据和sp
软件包的新手,并不了解这里发生了什么以及如何解决它。