经纬度到 UTM 坐标答案

【问题标题】：Latitude and Longitude to UTM coordinates经纬度到 UTM 坐标
【发布时间】：2013-02-14 21:29:02
【问题描述】：

我有一个来自 UTM-zone 北 33 的大纬度经度信息。

我尝试了以下命令将此地理信息转换为 UTM 坐标（我的数据集对象最初称为 S3km）：

library(rgdal)
UTM33N<-"+proj=utm+zone=33+north"
UTM33N<-paste(UTM33N,"+ellps=WGS84",sep="")
UTM33N<-paste(UTM33N,"+datum=WGS84",sep="")
UTM33N<-paste(UTM33N,"+units=m+no_defs",sep="")
coord.UTM33N<-project(as.matrix(S3km[,c("Longitude","Latitude")]),UTM33N)

我收到以下错误消息：

Error in project(as.matrix(S3km[,c("Longitude","Latitude")]),UTM33N):
no arguments in initialization list.

有谁知道问题出在哪里？我下载了最新的 R 版本（即 R 2.15.2），rgdal-package 也是新下载的。

【问题讨论】：

您使用什么语言/库？请在您的问题中添加相关标签。

标签： r utm rgdal

【解决方案1】：

您使用的投影信息格式似乎不正确。这可能会导致函数无法识别投影字符串中的参数。如 rgdal 中所述，投影信息必须符合 PROJ.4 文档（即，+= 和 空格分隔参数 之间没有空格。例如：“+proj=lcc +lat_1=48 +lat_2 =33 +lon_0=-100" 将您的粘贴函数参数更改为 sep=" " 可能会解决此问题。

【讨论】：

【解决方案2】：

您的代码似乎至少有几个问题：

正如 Lucas 指出的，PROJ4 字符串在参数之间需要空格，因此请使用 sep = " "（paste() 的默认值）而不是 sep = ""。
此外，sp 和 rgdal 包中的函数需要将 proj4strings 包装在对 CRS() 实用程序函数的调用中。

这是一个你应该能够适应你的情况的工作示例：

library(rgdal)

## Create an example SpatialPoints object
pts <- SpatialPoints(cbind(-120:-121, 39:40), 
                     proj4string = CRS("+proj=longlat +datum=NAD27"))

## Construct a proper proj4string
UTM11N <- "+proj=utm +zone=11 +datum=NAD83 +units=m +no_defs"
UTM11N <- paste(UTM11N, "+ellps=GRS80 +towgs84=0,0,0")
UTM11N <-  CRS(UTM11N)

## Project your points
ptsUTM <- spTransform(pts, UTM11N)

## Check that it worked
ptsUTM
# SpatialPoints:
#      coords.x1 coords.x2
# [1,]  240111.6   4321052
# [2,]  158420.9   4435418
# Coordinate Reference System (CRS) arguments: +proj=utm +zone=11
# +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

【讨论】：

非常感谢你们两个！！ :) 通过适当的间距和 Josh 改编的示例，我终于设法转换了我的坐标。