文章/答案/技术大牛

发布

社区首页 >问答首页 >纬度和经度到UTM坐标

问纬度和经度到UTM坐标
EN

Stack Overflow用户

提问于 2013-02-28 23:49:40

回答 2查看 2.8K关注 0票数 0

我有来自UTM北纬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也是新下载的。

utm

rgdal

回答 2

Stack Overflow用户

回答已采纳

发布于 2013-03-01 01:32:15

看起来你的代码至少有几个问题：

正如卢卡斯所指出的，PROJ4字符串在参数之间需要空格，因此使用sep = " " (paste()的默认值)而不是sep = ""。
此外，E19和E111包中的函数E210和E212期望proj4strings包装在对D13实用程序函数的调用中。H214

下面是一个工作示例，您应该能够适应您的情况：

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

票数 3

Stack Overflow用户

发布于 2013-03-01 01:09:24

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

票数 2

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/15139791

复制

相似问题

问纬度和经度到UTM坐标
EN

回答 2

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问纬度和经度到UTM坐标EN

回答 2

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问纬度和经度到UTM坐标
EN