我有一组包含纬度和经度的数据集,我想将其转换为伊利诺伊州东部的州平面坐标,使用EPSG 2790 (http://spatialreference.org/ref/epsg/2790/)或者ESRI 102672 (http://spatialreference.org/ref/esri/102672/)。
这个问题之前肯定已经被问过了,我的代码基于这里的答案 ("Non Finite Transformation Detected" in spTransform in rgdal R Package 和 http://r-sig-geo.2731867.n2.nabble.com/Converting-State-Plane-Coordinates-td5457204.html)。
但出于某种原因,我无法让它运行:
library(rgdal)
library(sp)
data = data.frame(long=c(41.20,40.05), lat=c(-86.14,-88.15))
coordinates(data) <- ~ long + lat
proj4string(data) <- CRS("+init=epsg:4326") # latitude/longitude
data.proj <- spTransform(data, CRS("+init=epsg:2790")) # illinois east
提供:
non finite transformation detected:
long lat
41.20 -86.14 Inf Inf
Error in spTransform(data, CRS("+init=epsg:2790")) : failure in points 1
In addition: Warning message:
In spTransform(data, CRS("+init=epsg:2790")) :
2 projected point(s) not finite