中科院1km土地利用数据处理流程

2022-02-09 13:54:01 浏览数 (1)

数据下载

代码语言:javascript复制
下载地址:https://www.resdc.cn/Default.aspx(需注册)

下载年份为:1995,2000年的吉林省,辽宁省共4份文件(非VIP会员受限制,一天仅可下载5个文件)

数据处理

处理流程如下:

  • ①数据拼接
  • ②裁剪
  • ③空白值处理
  • ④叠置(计算土地流转面积)
①数据拼接(以1995年为例)

工具:【镶嵌至新栅格】,拼接1995年吉林省与辽宁省

②按研究区裁剪图层

(若两图层投影不一致,裁剪前需先投影,此处投影已同),操作时在[环境]可设置extent,精度更高(如下图对比)

提取研究区范围

裁剪范围对比

未在环境设置extent

在环境中已设置extent

目标图层范围

注:在操作前,因先看空值个数,再判断是否需要赋值,本例由于空白值数目少,对结果影响甚微,故不做赋值也可

③空白值处理

为显示空白值,此处已对图层符号化,如图所示

左侧可知,所有属性值都赋黑色,为什么还会有空白?

点击属性查看空白值,为NoData,同时可观察空白有呈“带状”

敲黑板! 进行操作前,首先检查数据!一起来倒带 添加数据,并放大查看数据

放大局部图

实际上,空白值为边界未接合导致! 理论上,空白处应有值,根据地理学第一定律(地理环境越相似,地理特征越相近),故采用栅格计算器进行赋值

操作后

  • 可多次采用3*3像元大小迭代,直至nodata个数在一定阈值
  • 该阈值自行设置,迭代多次后,若发现nodata个数不发生变化,可尝试5*5像元大小,7*7像元大小
  • 迭代次数:不建议多次,随着像元设置越大,该像元取值失真
  • 在处理nodata值中,边界处也被赋值,故需裁剪(注:此处并没有采用设置环境处理范围,因像元会发生些许偏移)

裁剪如下

裁剪因赋值nodata而导致的边界扩张

对比空白值处理前后个数

nodata像元个数无法通过属性表得到,需设置为0,通过查看0值个数间接得nodata个数,对nodata进行赋0值

计算去空值图层剩余空值个数

由于赋值是对包含整个研究区的外接矩形进行,故需裁剪(需设置extent),裁剪后对比图如下

空值对比图

注:在操作前,因先看空值个数,再判断是否需要赋值,本例由于数据少,故不做赋值也可

④图层叠置

2000年数据经与上操作后(尤其注意裁剪设置范围extent,否者叠置会出现像元不完全重合),与1995年数据进行叠置,如下图

叠置计算流转

叠置结果

从结果可知,除了两者共同位置剩余5个nodata未叠置,其余像元地类代码均已匹配

新建字段,拆分地类代码

字段计算器对count值拆分地类代码

计算结果

将属性表导出即可计算不同地类之间转换

导出属性表

使用excel打开文件

使用EXCEL打开叠置属性表

删除属性为0值,并使用数据透视表计算流转情况

使用数据透视表查看土地流转情况

使用透视表时,注意行列选择

使用数据透视表查看土地流转情况

小结:

  1. 操作前先检查数据(如投影一致?拼接无缝?同一区域count值同?)
  2. 裁剪注意设置环境extent
  3. 空值,先判断是否需要进行处理,若是,需要迭代的次数需自行把握
  4. 叠置最容易出现像元偏移而导致地类代码在进行栅格计算时未能重合,叠置前需检验像元是否重合,count值是否相同。

0 人点赞