在R中将shapefile转换为x，y，z(经纬度坐标)
在R中将shapefile转换为x，y，z(经纬度坐标)
提问于 2018-06-06 05:24:17
我有一个shapefile，在提取x，y，z分量时遇到问题。我还想将x和y从UTM转换为Lat和Lon。
最终我想要
plot(x, y)
map(add=T)并添加其他点。
对于ex，你如何理解这个投影：
library(PBSmapping)
r <- importShapefile("~/Desktop/BIO_CA_KelpPersistent/BIO_CA_KelpPersistent.shp")
str(r)
- attr(*, "prj")= chr "PROJCS[\"NAD_1983_California_Teale_Albers\",GEOGCS[\"GCS_North_American_1983\",DATUM[\"D_North_American_1983\",SPHEROID[\"GRS_1"| __truncated__更新：
require(sf)
nc <- st_read("BIO_CA_KelpPersistent.shp")
nc_points <- st_coordinates(st_geometry(st_centroid(nc)))
xy <- project(nc_points, paste("+proj=aea +lat_1=34 +lat_2=40.5 +lat_0=0 +lon_0=-120 +x_0=0 +y_0=-4000000 +datum=NAD83 +units=m +no_defs", sep = ""), inverse=T)
plot(xy, ylim=c(32, 33), xlim=c(-119, -117))
map(add=T, lwd=2)但是很多地方都没有了。如何正确地对其进行索引？
谢谢!
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/50709149
复制相关文章
相似问题