在将坐标转换为UTM时，R会导出意外结果

Question

我想使用R将纬度/经度坐标转换为UTM（区域47）系统。我终于使用the last post on the bottom设法得到没有错误消息的结果。但结果并不像预期的那样。 这是一个例子：

    LongLatToUTM<-function(x,y,zone){
  xy <- data.frame(ID = 1:length(x), X = x, Y = y)
  coordinates(xy) <- c("X", "Y")
  proj4string(xy) <- CRS("+proj=longlat +datum=WGS84")  ## for example
  res <- spTransform(xy, CRS(paste("+proj=utm +zone=",zone," ellps=WGS84",sep='')))
  return(as.data.frame(res))
}

# Execute
x<-100.4605
y=21.568216
res.1<-LongLatToUTM(x,y,47)

我得到以下结果：x = 646172.7 y = 2876501带R，但我预计x = 651218 y = 2385744。我做错了什么？