土地利用转移矩阵library(dplyr)
library(tidyr)
library(stringr)
#在excel表中选取全部数据区域,Ctr+C
data
data %>% mutate(GRID05 = str_sub(GRIDCODE, start = 1, end = 1), GRID10 = str_sub(GRIDCODE_1, start = 1, end = 1)) %>%
filter(GRID05 != 0 & GRID10 != 0) %>%
mutate(GRID05 = factor(as.numeric(GRID05), levels = 1:6, labels = c("林地", "草地", "湿地", "耕地", "人工表面", "其他"))) %>%
mutate(GRID10 = factor(as.numeric(GRID10), levels = 1:6, labels = c("林地", "草地", "湿地", "耕地", "人工表面", "其他"))) %>%
group_by(GRID05, GRID10) %>% summarise(Shape_Area = sum(Shape_Area)) %>%
spread(key = GRID10, value = Shape_Area)