Я новичок в работе с R, и в настоящее время я пытаюсь преобразовать координаты широты и долготы, которые я загрузил в R из документа xlsx, в CRS растра, на котором я хочу построить свои координаты. Я использую RStudio версии 1.2.5033. Вместо того, чтобы связывать файл координат, вот загруженные мной координаты широты и долготы (сначала широта, затем долгота): 1 Allequa sh (46,03323 -89,60874) 2 Bassett (44,99034 -93,27385) 3 Chippewa (44,93223 -91,38599) 4 Como (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)
* 1002 сначала опубликуйте мой код, а затем go, с какими проблемами я столкнулся
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. Когда я проверяю свои координаты после их загрузки, а затем назначаю им crs используя st_as_sf, я получаю следующее предупреждающее сообщение, и я не уверен, почему: В st_is_longlat (x): ограничивающая рамка потенциально имеет недопустимый диапазон значений для данных longlat
Когда я затем go преобразовываю свои данные координат в соответствии с данными растра с помощью st_transform, я затем конвертирую только первый сайт (Allequa sh, и даже это преобразование не кажется правильным), в то время как другое координаты сайта остаются пустыми.
Любая помощь приветствуется!