坐标上传和CRS转换

问题描述 投票:0回答:1

我对使用R刚起步,目前正在努力将上传到R中的经/纬度坐标从xlsx文档转换为要在其上绘制坐标的栅格的CRS。我正在使用RStudio版本1.2.5033。而不是链接坐标文件,这里是我上传的经/纬度坐标(先经纬度然后经纬度):1个Allequash(46.03323 -89.60874)2巴塞特(44.99034 -93.27385)3奇珀瓦(44.93223 -91.38599)4科莫(45.10457 -91.53066)5 Flambeau(45.41387 -91.21639)6 Minnehaha(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数据的值范围可能无效]

  1. [然后我使用st_transform转换坐标数据以匹配栅格的坐标数据时,我只转换了第一个站点(Allequash,甚至这种转换似乎都不正确),而其他站点坐标保持空白。] >
  2. 感谢您的任何帮助!

我对使用R刚起步,目前正在努力将上传到R中的经/纬度坐标从xlsx文档转换为要在其上绘制坐标的栅格的CRS。我是...

r coordinates data-conversion
1个回答
0
投票

coords格式为coords = c("x", "y"),所以只需反转纬度和经度,就可以使用。

© www.soinside.com 2019 - 2024. All rights reserved.