rgdal封装横向/纵向->；UTM

有人能看出这有什么问题吗：

library(rgdal)

# Make a two-column matrix, col1 = long, col2 = lat
xy <- cbind(c(-107), c(26))
# Convert it to UTM coordinates (in units of meters)
project(xy, "+proj=utm +zone=51 ellps=WGS84")

---

谢谢

使用正数表示经度（107，而不是-107）。负数表示西半球

你是说（-107,26）的UTM坐标与（107,26）的UTM坐标相同吗？我的意思是（lon，lat）=（-107,26）按惯例是指西107度，北26度。UTM 51区位于菲律宾附近的某个地方的顶部。所以我假设你的坐标是东半球。您必须为坐标选择正确的UTM区域，如果该区域太过偏离，则无法像本例中那样进行投影。

xy我的问题与此相关：。因此，我必须为每个经度确定正确的区域。我知道有60个区域。是否没有更简单的方法将我的LAT/LONG传输到UTM？谢谢。我查了一下Mollweide或Aitoff的投影，它是一种地图投影，用相等的面积覆盖了整个地球。通过将球形地球展开为矩形，无法避免扭曲几何体。你保留一些属性（例如，角度、大小），根据这个家伙的说法，你最好保留面积。如果你有你最喜欢的教科书，你可以仔细查看其中地图的标题，看看是否有提到地图投影。为什么是-1？这是一个正确的问题。
In project(xy, "+proj=utm +zone=51 ellps=WGS84") :
  1 projected point(s) not finite