【发布时间】:2020-05-28 21:21:23
【问题描述】:
我对在 R 中工作还很陌生,目前我正在努力将我上传到 R 中的经纬度坐标从 xlsx 文档转换为我想要绘制坐标的栅格的 CRS。我正在使用 RStudio 版本 1.2.5033。这里没有链接坐标文件,而是我上传的纬度/经度坐标(先经纬度): 1 阿勒夸什 (46.03323 -89.60874) 2 巴塞特 (44.99034 -93.27385) 3 奇佩瓦 (44.93223 -91.38599) 4 科莫 (45.10457 -91.53066) 5 弗朗博 (45.41387 -91.21639) 6 明尼哈哈 (44.918 -93.22536) 7 红雪松 (44.88263 -91.93599) 8 瓦 (45.03322 -93.28645) 9 鳟鱼 (45.09951 -91.65159)
我会先发布我的代码,然后再讨论我遇到了什么问题
site_coordinates<-read_excel(coordinates,sheet=1) #Proj=longlat datum=WGS84
rm(coordinates)
site_coordinates_crs<-CRS("+proj=longlat +datum=WGS84")
coordinates<-st_as_sf(site_coordinates,coords=c("Latitude","Longitude"),crs=site_coordinates_crs)
st_crs(coordinates)
coordinates_aes<-st_transform(coordinates,crs(nlcd))
# nlcd is the raster I'm plotting the coordinates on top of, uploaded online nlcd 2016 land use raster
我遇到的问题: 1. 当我在上传坐标后检查我的坐标,然后使用 st_as_sf 为它们分配 crs,我收到以下警告消息,我不知道为什么: 在 st_is_longlat(x) 中: 边界框可能对 longlat 数据的值范围无效
- 然后,当我使用 st_transform 转换我的坐标数据以匹配栅格数据时,我只转换了第一个站点(Allequash,甚至此转换似乎也不正确),而其他站点坐标保持为空。
感谢任何帮助!
【问题讨论】:
标签: r coordinates data-conversion