Я очищаю файл sf с помощью пакета sf. Когда я записываю файл, используя write_sf, там есть вся геометрия. Но когда я читаю его обратно с использованием readOGR, он сообщает: «Удаление нулевой геометрии» и сбрасывает 430 строк. Мне нужно прочитать его обратно с OGR, потому что я затем привязываю его к другому объекту. То, как я делаю подмножество, довольно шаткое - это то, что заставляет его сбрасывать геометрию? Любая помощь приветствуется!
'fsvtrt <- read_sf("C:/Users/A02305035/Desktop/JamelaResearch/Data/GIS/FS/fsvtrt.shp") #raster package'
Добавить ObjectID
fsvtrt$OBJECTID <- seq(1:5464)
Удалить неполные процедуры
fsunwanted1 <- fsvtrt$DATE_COMPL[is.na(fsvtrt$DATE_COMPL)]
incomplete <- fsvtrt[fsvtrt$DATE_COMPL %in% fsunwanted1,]
salvagetx <- incomplete[na.omit(incomplete$STAGE_VALU),]
unique(incomplete$STAGE_VALU)
Найти и удалить запланированные будущие процедуры
planned <- fsvtrt$DATE_PLANN[fsvtrt$DATE_PLANN>as.Date("2020-02-13")]
Определение запланированных будущих обработок
fsvtrt <- fsvtrt[!fsvtrt$DATE_PLANN %in% planned,]
Показать обработанные процедуры, исключая запланированные
Удалить неполные обработки
incomplete <- fsvtrt$DATE_COMPL[is.na(fsvtrt$DATE_COMPL)]
fsvtrt <- fsvtrt[!fsvtrt$DATE_COMPL %in% incomplete,]
Удалить пожары
wildfires <- fsvtrt$ACTIVITY[fsvtrt$ACTIVITY == "Wildfire - Natural Ignition"]
fsvtrt <- fsvtrt[!fsvtrt$ACTIVITY %in% wildfires,]
wildfires1 <- fsvtrt$ACTIVITY[fsvtrt$ACTIVITY == "Wildland Fire Use"]
fsvtrt<-fsvtrt[!fsvtrt$ACTIVITY %in% wildfires1,]
wildfires2 <- "Wildfire - Fuels Benefit"
fsvtrt <- fsvtrt[!fsvtrt$ACTIVITY %in% wildfires2, ]
wildfires3 <- "Wildfire - Human Ignition"
fsvtrt <- fsvtrt[!fsvtrt$ACTIVITY %in% wildfires3, ]