如何利用R语言做一个土地利用转移矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

change0510.xlsx (30.86 KB)

2016-12-19 16:19:14 上传

转移矩阵.xls (28.5 KB)

2016-12-19 16:19:13 上传

原始数据如下：GRIDCODE和GRIDCODE_1是土地利用类型
101-112林地
21-24草地
31-37湿地
41-42耕地
51-54人工表面
61-69其他
根据GRIDCODE筛选，然后汇总Shape_Area

Excel制作结果：


我想知道如何用R编程实现啊，有没有大神指导一下，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：土地利用 R语言 idcode EXCEL Shape R语言

1. library(dplyr)
   
2. library(tidyr)
   
3. library(stringr)
   
4. #在excel表中选取全部数据区域，Ctr+C
   
5. data <- tbl_df(read.table("clipboard", header = TRUE))
   
6. data %>% mutate(GRID05 = str_sub(GRIDCODE, start = 1, end = 1), GRID10 =  str_sub(GRIDCODE_1, start = 1, end = 1)) %>%
   
7.     filter(GRID05 != 0 & GRID10 != 0) %>%
   
8.     mutate(GRID05 = factor(as.numeric(GRID05), levels = 1:6, labels = c("林地", "草地", "湿地", "耕地", "人工表面", "其他"))) %>%
   
9.     mutate(GRID10 = factor(as.numeric(GRID10), levels = 1:6, labels = c("林地", "草地", "湿地", "耕地", "人工表面", "其他"))) %>%
   
10.     group_by(GRID05, GRID10) %>% summarise(Shape_Area = sum(Shape_Area)) %>%
    
11.     spread(key = GRID10, value = Shape_Area)

复制代码

你好，请问一下你是怎么用excel做出转移矩阵的呢？求指教啊谢谢谢谢~！

Palette221 发表于 2018-4-18 21:41
你好，请问一下你是怎么用excel做出转移矩阵的呢？求指教啊谢谢谢谢~！

大致就是1楼说的，先根据地类汇总，汇总计算面积，挨个汇总，前后相减，大概是这个思路，具体操作我也不记得了，好久之前做的

Excel数据透视表