我对使用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数据的值范围可能无效]
感谢您的任何帮助!
我对使用R刚起步,目前正在努力将上传到R中的经/纬度坐标从xlsx文档转换为要在其上绘制坐标的栅格的CRS。我是...
coords格式为coords = c("x", "y")
,所以只需反转纬度和经度,就可以使用。