это относится к последнему вопросу, который я задал (все еще относительное новшество ie для R, с использованием RStudio, версия 1.2.5033). По сути, у меня есть растровый файл, который я скачал в Интернете (землепользование NLCD 2016), и я хочу построить шейп-файлы, которые очерчивают бассейны рек и координаты участков отбора проб поверх растровой карты. Я могу построить шейп-файлы поверх растра, но по какой-то причине я не могу получить координаты для построения поверх растра. RStudio вместо этого создает новую фигуру для координат. Может ли кто-нибудь помочь решить это? Спасибо!
Информация о растре (начинается в коде как nlcd, а затем переклассифицируется и обрезается до r c): class: Размеры RasterLayer: 7302, 11233, 82023366 (nrow, ncol, ncell ) разрешение: 30, 30 (x, y) протяженность: 175185, 512175, 2386365, 2605425 (xmin, xmax, ymin, ymax) crs: + proj = aea + lat_1 = 29,5 + lat_2 = 45,5 + lat_0 = 23 + lon_0 = -96 + x_0 = 0 + y_0 = 0 + ellps = WGS84 + towgs84 = 0,0,0, -0, -0, -0,0 + units = m + no_defs источник: / private / var / folder / f5 / 7wtm6gsn1qn70rxh1ny6m7m40000gn / T / RtmpRxBkDJ / raster / r_tmp_2020-05-28_090052_49820_66166.grd имена: NLCD_2016_Land_Cover_L48_20190424 значения: 1, 100 * 100 (* минимум, одна форма) информация о характеристиках: 1 объект и 1 тип геометрии поля: POLYGON Размер: XY bbox: xmin: 338204,4 ymin: 2456656 xmax: 350204,4 ymax: 2471206 CRS: + proj = aea + lat_1 = 29,5 + lat_2 = 45,5 + lat_0 = 23 + lon_0 = -96 + x_0 = 0 + y_0 = 0 + ellps = WGS84 + towgs84 = 0,0,0, -0, -0, -0,0 + units = m + no_defs Имя геометрии 1 WI2 0200528013216810000 POLYGON ((346454.4 2471206, ...
Информация о координатах: Простая коллекция объектов с 9 объектами и 1 типом геометрии поля: Размер ТОЧКИ: XY bbox: xmin: 213539 ymin: 2439512 xmax: 495172.2 ymax: 2577206 CRS: + proj = aea + lat_1 = 29,5 + lat_2 = 45,5 + lat_0 = 23 + lon_0 = -96 + x_0 = 0 + y_0 = 0 + ellps = WGS84 + towgs84 = 0,0,0, -0, -0 , -0,0 + units = m + no_defs
Таблица: 9 x 2
Геометрия имени 1 Allequa sh (495172.2 2577206) 2 Bassett (214667.1 2447725) 3 Chippewa (363545.1 2447006) 4 Комо (351253,3 2465617) 5 Фламбо (374201,5 2501172) 6 Миннехаха (218719,3 2439788) 7 Красный кедр (320474 2439512) 8 Гонт (213539 2452466) 9 Форель (341781,7 2464613)
All_Bas<-st_read("./Allequash/layers","globalwatershed")
Chip_Bas<-st_read("./Chippewa/layers","globalwatershed")
Bass_Bas<-st_read("./Bassett/layers","globalwatershed")
Como_Bas<-st_read("./Como/layers","globalwatershed")
Flam_Bas<-st_read("./Flambeau/layers","globalwatershed")
Minn_Bas<-st_read("./Minnehaha/layers","globalwatershed")
Red_Bas<-st_read("./Red Cedar/layers","globalwatershed")
Shin_Bas<-st_read("./Shingle/layers","globalwatershed")
Trout_Bas<-st_read("./Trout/layers","globalwatershed")
All_Bas_aea<-st_transform(All_Bas,crs(nlcd))
Chip_Bas_aea<-st_transform(Chip_Bas,crs(nlcd))
Bass_Bas_aea<-st_transform(Bass_Bas,crs(nlcd))
Como_Bas_aea<-st_transform(Como_Bas,crs(nlcd))
Flam_Bas_aea<-st_transform(Flam_Bas,crs(nlcd))
Minn_Bas_aea<-st_transform(Minn_Bas,crs(nlcd))
Red_Bas_aea<-st_transform(Red_Bas,crs(nlcd))
Shin_Bas_aea<-st_transform(Shin_Bas,crs(nlcd))
Trout_Bas_aea<-st_transform(Trout_Bas,crs(nlcd))
Basins<rbind(Chip_Bas_aea,Minn_Bas_aea,All_Bas_aea,Bass_Bas_aea,Como_Bas_aea,Flam_Bas_aea,Red_Bas_aea,Shin_Bas_aea,Trout_Bas_aea)
nlcd_crop<-crop(nlcd,Basins)
coordinates <-"Site Coordinates_degrees.xlsx"
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("Longitude","Latitude"),crs=site_coordinates_crs)
coordinates_aes<-st_transform(coordinates,crs(nlcd))
m<-c(0,10,0,11,11,1,12,20,0,21,24,2,25,40,0,41,43,3,44,80,0,81,87,4,88,89,0,90,95,5)
rclmat<-matrix(m,ncol=3,byrow=TRUE)
rc<-reclassify(nlcd_crop,rclmat)
pal<-colorRampPalette(c("gray90","gray47","cornflowerblue","gray90","lightgreen","gray90","gray90","gray90","red","lightgreen"))
jpeg("Simplifed Landuse code.jpeg",width=1000,height=750)
plot(rc,col=pal(50),legend=FALSE)
plot(All_Bas_aea,color="white",lwd=2,fill=NA,add=TRUE)
plot(Chip_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Minn_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Bass_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Como_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Flam_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Red_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Shin_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(Trout_Bas_aea,color="black",lwd=2,fill=NA,add=TRUE)
plot(coordinates_aes) # This is where I'm running into the problem of coordinates not plotting on top of raster map
dev.off()'''