精华内容
下载资源
问答
  • 介绍如何利用ARCGIS生成土地利用转移矩阵的两种方法
  • 如果你要做土地利用变化的研究,那就会用到土地转移矩阵。而土地利用转移矩阵常规有两个办法。 第一种办法就是利用ENVI里面的change detection statics工具。 第二种办法就是利用arcgis里面的相交工具,但这个是...

    如果你要做土地利用变化的研究,那就会用到土地转移矩阵。而土地利用转移矩阵常规有两个办法。

    第一种办法就是利用ENVI里面的change detection statics工具。

    第二种办法就是利用arcgis里面的相交工具,但这个是基于矢量的。

    这个博客我介绍第三种办法,使用arcgis自带的栅格计算器制作土地利用转移矩阵。

    该方法的原理简单介绍:将第一年的栅格土地利用类型乘以100,再加上第二年的栅格土地利用栅格数据。

    即:第一年*100+第二年。

    下面展示例子:我有某市2010年与2017年30m的土地利用栅格数据。其中每个代码的含义如下:

    注意一下,2010年多一个第11类的cloud云类别。

    某市2010年与2017年的土地利用图如下所示:

    在你已经有研究区两景影像时,

    第一步:打开栅格计算器,将第一年的乘以100,再加上第二年的,如下图所示。这样得到的每个像素的像素值就是土地转移矩阵。比如说在某个像素上,如果第一年类别为耕地1,第二年变为了林地2,那经过计算该值为102.

    第二步,打开arcgis的属性表,并将属性表导出,转为excel。我们需要注意到,value值就是转移类型,count就是像元素。

    3.在excel中整理土地转移矩阵。比如像下面一样:把数值重新整理,然后依次填入对应的空格。

     

    到此利用arcgis进行土地转移矩阵的计算完毕,该方法可直接利用网上的土地利用数据进行计算。而不需要额外的转为矢量,更为方便。希望能够帮到你。(作者:快多宝)

     

     

     

     

    展开全文
  • 土地覆被变化制图与土地利用转移矩阵土地利用数据来源自己解译哨兵数据易康面向对象SVM分类eCognition(易康)面向对象分类与分类结果处理网上各种土地覆被数据清华宫鹏课题组:http://data.ess.tsinghua.edu.cn/...

    7f5cc86d5687bf1b88b50243aec22d26.png

    土地覆被变化制图与土地利用转移矩阵

    土地利用数据来源

    1. 自己解译
    5b6e94bd0d050dc4b2323b4bc6e8abd7.png
    哨兵数据易康面向对象SVM分类
    2343f11b2188b960c1ad40594e8246fc.png
    eCognition(易康)面向对象分类与分类结果处理
    1. 网上各种土地覆被数据
    • 清华宫鹏课题组:http://data.ess.tsinghua.edu.cn/
    • 中科院资源环境科学与数据中心:http://www.resdc.cn/data.aspx?DATAID=99
    • GlobalLand30: http://www.webmap.cn/commres.do?method=globeIndex
    • MCD12Q1:https://e4ftl01.cr.usgs.gov/MOTA/MCD12Q1.006/
    • 国家地球系统科学数据中心土地覆被专题库:http://www.geodata.cn/data/index.html?categoryId=18

    数据源与制图目的

    已有数据

    已有2010,2018年两个年份的矢量土地覆被图

    252b538085cbe6b2b35af961f2a9dabc.png
    图1. 黄河口2010,2018主要植物群落和土地覆被分布

    分析结果

    根据上面的土地覆被图,绘制下面的侵蚀变化图,植物群落与地类变化图。制作土地利用转移矩阵。

    06bb8fb5e1cdfc3065cb82a45f349e15.png
    图2. 三角洲侵蚀变化图
    7f1695e5a859e48cd35b049c51c7d9d4.png
    图3. 黄河三角洲植物群落与地类变化图
    914a53a008311aae9373661261236959.png
    表1. 黄河口群落与土地利用转移矩阵

    三角洲侵蚀变化图制作

    总体思路

    • 提取水陆边界
    • 提取2010,2018水陆边界交集
    • 2010,2018土地覆被求并集
    • 土地覆被并集-水陆边界交集,求侵蚀变化情况

    具体做法

    1. 根据图1绘制水陆分界
    548a10592433b582727c93f9eb1b5ee1.png
    图4. 水陆分界图
    1. 将图4两个矢量求交集,求得2010,2018均为陆地的区域(输出要素类:OLI水陆分界_Intersect)
    2. 将图1两个矢量求交集(输出要素类:TM2010_Intersect1)
    3. 在TM2010_Intersect1减去OLI水陆分界_Intersect(输出要素类:TM2010_Intersect1_Erase)
    4. 新建一个字段,计算侵蚀淤积变化情况
    351e76d61bd3f64cd6c061c42f7cc050.png
    使用python进行计算

    需要注意的是,我前面用的是汉字,下面代码里面判断的是地类编码,由于ArcGIS中Python2对于汉字支持不太好,所以还是建议判断数字。

    def Cal(a,b):
      if a==b:
        return "无变化"
      elif a ==1 and b !=a:
        return "淤积"
      else:
        return "侵蚀"
     
     
    Cal(  !CODE! , !CODE_1! )
    1. 最后对Change字段进行分类符号化即可。
    99cf86dd96c5e01de0fe169f1d9cbd35.png
    变化制图

    黄河三角洲植物群落与地类变化图

    以水域为例,介绍地类变化图制作

    • 按属性选择,筛选2010,2018年均为水域的所有要素,并导出要素类
    • 新建一个Change字段(文本型),计算变化情况
    • 计算代码和上面侵蚀淤积代码类似
    • 转入转出分类符号化
    a84ab992a214e5e4f7a3e651f91d574c.png
    根据转入转出分类符号化

    土地利用转移矩阵

    • 将前面TM2010_Intersect1要素类属性表导出
    2554de3f01952ab7bb44998206f5c898.png
    属性表导出为文本型
    • NOTEPAD++改一下ANSI编码,防止EXCEL打开乱码
    • 数据透视表制作转移矩阵,行2010年地类,列2018年地类,求和项为面积
    f056dcf2f13f98571be8f66a24dcdc94.png
    制作数据透视表
    • 由于数据使用的是投影坐标系,单位为米,对应的面积单位为平方米,因此土地利用转移矩阵需要进一步进行单位换算,换算为所需要的单位。

    f2215886c998816c1b73153d171e3b50.png

    展开全文
  • 土地利用转移矩阵生成的几种方法土地利用转移矩阵生成的几种方法根据你的数据类型选用不同的数据生成方法若你的数据是Raster格式:则有如下方法1ErdasImagine----Interpreter---GisAnalysis---Matrix,输入两个时相...

    土地利用转移矩阵生成的几种方法

    土地利用转移矩阵生成的几种方法

    根据你的数据类型选用不同的数据生成方法

    若你的数据是

    Raster

    格式:则有如下方法

    1

    Erdas

    Imagine----Interpreter---Gis

    Analysis---Matrix

    ,输入两个时相的

    Raster

    数据即可

    做这一步之前记得先对两时相的数据进行重编码(

    nterpreter---Gis

    Analysis---Recode

    )

    一般运行如果出现错误

    肯定是重编码没做好,请继续查证。

    2

    先在

    Erdas

    中利用

    Modeler

    计算

    如下公式

    NC(I,J)=NC(I)*10+NC(J),(J>I)

    其中:

    NC(I,J)

    表示

    i

    j

    两年份的土地利用变化图;

    NC

    (

    i

    )表示

    i

    年份遥感

    分类影像;

    NC

    (

    j

    )表示

    j

    年份的遥感分类影像。

    在此计算的基础上,将以上变化影像图转化为

    BIL

    格式,再利用

    ARC/INFO

    GRID

    模块将影像转为

    GRID

    格式,然后利用

    GRID

    模块中的属性

    表(

    vat

    )查看命令对影像灰度值进行统计,最后得出土地利用转化举

    证。(注

    :

    此方法本人尚未实现过,不知可行否)。

    若数据是

    Vector

    格式

    1

    Erdas

    Imagine----Interpreter---Gis

    Analysis---Matrix

    ,输入两个时相的

    Vector

    数据即可

    此时注意

    输出栅格大小不应设的太小

    要不一运行就会提示你的空间不

    做这一步之前,请做好前期的地理编码。

    2 ArcView3.3

    加载

    spatial analysis

    模块

    把两时相的

    Vector

    图转成

    grid

    格式(当然中间有一些单位的设置根据你

    做的图的分辨率来设置即可)

    analysis---mapcaculate

    直接计算即可。

    3

    把两期解译完的

    Vector

    文件在

    arctoolbox——overlay——union

    中叠加,

    注意:两个文件不能用同一个字段名,比如一个用

    93Type

    ,另一个时相

    则用

    00Type

    叠加后的文件在

    Arcmap

    中打开,选中文件,然后点右键

    ——Property

    ——

    空间查询,输入条件语句,比如:

    93Type=‘1’And

    00Type=‘2’

    ;查询

    结果即为第一种类型转化为第二种类型的图形,可以另建一图层比如:

    12

    ,把查询结果复制到

    12

    图层上。统计出面积,依进行,就可以得到土

    地利用类型转移矩阵

    展开全文
  • 下面介绍如何在ArcGIS中建立转移矩阵一、数据准备(图1)准备两幅不同时相的土地利用现状图(shp格式),每幅图的属性表都要有一个表示土地利用类型的字段,并且要使用不同的名称加以区分,如Type1995,Type2000。土地...

    下面介绍如何在ArcGIS中建立转移矩阵

    一、数据准备(图1)

    准备两幅不同时相的土地利用现状图(shp格式),每幅图的属性表都要有一个表示土地利用类型的字段,并且要使用不同的名称加以区分,如Type1995,Type2000。土地利用类型名称必须统一,并且完整,如都使用“城镇用地”、“有林地”等。

    二、数据融合(图2)

    在ArcMap里分别打开两个时相的图层,打开ArcToolbox,选择Data Management Tools | Generalization | Dissolve工具。Input Feature选择要融合的图层,Output Feature Class选择输出结果存储的位置及名称,Dissolve Field(s)选择土地利用类型字段(如Type1995),然后勾选Creat multipart features选项,点击OK完成。重复此过程,对另一时相数据进行融合。此步骤使相同利用类型的记录融合为一个记录,以提高后面步骤的计算速度。

    三、叠置分析(图3)

    在ArcMap中打开两个时相融合后的数据,在ArcToolbox中选择Analysis Tools | Overlay | Intersect工具,Input Features选择两个时相的图层,Output Feature Class选择叠加结果存储的位置及名称,其余选项可以忽略,单击【OK】完成。

    四、计算面积并导出属性表(图4-6)

    在ArcMap中打开叠加后的图层数据,在该图层上右键打开属性表,选择Option | Add field… 新建一个字段,命名为NewArea。

    在Editer工具条中选择Editer | Start Editing,然后在属性表中NewArea字段上单击右键选择Calculate Geometry… ,在打开的Calculate Geometry对话框中,Property选择Area,Units选择要使用的面积单位,单击【OK】完成图斑面积计算。依次选择Editer | Save Edits / End Editing保存和退出编辑状态。

    在属性表中选择Option | Export… 将属性表保存为dbf文件。

    五、制作转移矩阵(图7-10)(以Excel2007为例)

    在Excel中打开上一步保存的dbf,另存为Excel格式并打开。在Excel中选中所有数据(不要点左上角,只选择有效数据),点击【插入】选项卡,选择【数据透视表】|【数据透视表】,点击【确定】。

    在打开的数据透视表中按图示将字段拖入相应区域。

    Excel自动计算矩阵,将该表稍事整饰就得到美观的土地利用转移矩阵。矩阵中r(I, j)就表示i类型向j类型转移的土地面积,空值表示i类型向j类型没有转移。

    展开全文
  • 匿名用户1级2018-01-10 回答ⅠENVI下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。1、准备数据两个时相的土地利用分类结果,它是单波段、专题类型的伪彩色图像(ENVI Classification)。2、计算...
  • 利用遥感和GIS来土地利用转移矩阵计算,非常方便,找了很久找到
  • %土地利用转移矩阵假设目前有两个不同年份的土地利用数据,分别是1990年和2000年的,要求不同土地利用类型间的相互转换关系,本文以matlab为例进行求解,非常的快速方面,代码如下[a,R]=geotiffread('D:\qixiang\...
  • 土地利用转移矩阵制作

    千次阅读 多人点赞 2019-07-11 18:03:00
    土地利用转移矩阵制作 1 马尔可夫模型(Markov) 1.1马尔可夫模型 马尔科夫模型(Markov)是俄国数学家安德烈·马尔科夫提出来的一种随机过程的研究,随机过程具有无后效性(即在系统状态转移过程中,系统将来的状态...
  • 数据基础:已有的土地利用转移矩阵数据处理平台:Excel目前得到的土地利用转移矩阵已经是一个很齐全的表格了,进行进一步分析的目的主要是为了计算新增建设用地占用农用地和新增建设用地占用耕地的情况,同时了解并...
  • 土地利用转移矩阵步骤

    千次阅读 2020-08-14 13:32:46
    一、2000年-2005年土地利用转移矩阵 1、转换工具-由栅格转出-栅格转面:grid转shp 2、合并之后,属性表添加字段:类型和面积(一定要带年份,方便后续处理) 3、Data Management Tools → Generalization → ...
  • 1.数据准备 30米全球地表覆盖数据GlobeLand30....再计算图斑面积 如何变化的-表转EXCEL 图面显示变化区域 1)图层属性表添加字段 2)计算 4)出图-2010-2020土地利用变化图 5)出对比图
  • Arcgis土地利用转移矩阵制作 数据准备** 土地利用数据是栅格格式,只有地类代表的数字,我们在栅格数据上添加一个NAME字段,把土地利用类型的名称加上去,比如1代表耕地,我们就在NAME字段中输入“耕地”。 然后...
  • arcgis制作土地利用转移矩阵

    千次阅读 多人点赞 2020-07-05 21:49:09
    1.首先下载土地利用的影像,在这里我用的是资源环境数据云平台的数据,分别为1980年与2015年的土地利用栅格影像,然后打开影像。 2.下载影像后进行重分类影像,因为影像原分类是二级分类,我这里需要的是一级分类...
  • 凤凰财经2019-07-09 01:301038ⅠENVI下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。1、准备数据两个时相的土地利用分类结果,它是单波段、专题类型的伪彩色图像(ENVI Classification)。2、...
  • 土地利用转移矩阵的几种实现方法

    万次阅读 多人点赞 2018-05-29 09:11:07
    土地利用转移矩阵生成的几种方法 查阅相关的资料,也没有得到土地利用类型转换矩阵确切的定义,我理解为不同时间段内同一区域内土地利用类型的相互转换关系,一般用二维表来表达,从二维表中可以快速查看各个地类间...
  • 这里是网上搜到的生成土地利用变化转移矩阵的几种方法,以飨来者:A 栅格数据做转移矩阵“一般习惯列为早期的数据,行为近期的数据,就如你说的早期的在上,晚期的在下。”B 矢量数据做转移矩阵根据你的数据类型选用...
  • 关于土地分类格式互转、土地利用转移矩阵计算详解 一、ENVI、ArcGIS格式互转 1.1 ENVI分类格式转TIF ENVI分类格式使用File下面的Save as另存为TIF格式: 对于部分5.3版本可能会出现报错: 那我们可以使用Classic...
  • 利用ArcGIS制作土地利用转移矩阵 经常遇到土地利用转移矩阵制作的问题,一段时间不接触,时间长了就容易忘记,而网上的教程参差不一,今天索性抽时间出一期教程,方便以后查询。01数据准备两期土地利用数据(栅格或...
  • 本文档将用图解操作详细说明利用ArcGIS对不同时期土地利用现状地图数据进行叠加统计分析,并利用Excel计算汇总,生成漂亮的转移矩阵
  • 假设有两期用地数据2012年和2017年,分成三种用地类型,首先得到这两期数据的栅格图像,每个栅格单元存储的值代表一种用地类型,然后以矩阵的方式读入内存中,再用numpy的flatten的方法将其一维化,最后调用python...
  • 转移矩阵操作

    2014-06-14 22:13:11
    土地利用变化常常利用转移矩阵进行衡量,转移矩阵有一下几种实现方式,已经亲身试验过,按照步骤操作即可。
  • nbsp手工制作ArcMap制作土地利用转移矩阵——图文详解版.pdf5页本文档一共被下载:次,您可全文免费在线阅读后下载本文档。 下载提示1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔...
  • 土地利用转移矩阵生成的几种方法 根据你的数据类型选用不同的数据生成方法 若你的数据是Raster格式:则有如下方法 1 Erdas Imagine----Interpreter---Gis Analysis---Matrix,输入两个时相的Raster数据即可 做这...
  • 深耕于土地利用变化的同学,可能经常在一些论文中看到土地模拟预测的方法,比如马尔科夫链、人工神经网络、clue-s模型、元胞自动机(CA)、最小二乘法、ANN-CA模型、PLS-PP模型等时间序列、马尔科夫(Markov)以及元胞...
  • 桑基图桑,基能量分流图,也叫桑基能量平衡图。...对于土地利用文献中大多使用转移矩阵的方式。本博客实现2000年,2005年,2010年,2015年,2019年土地利用之间的转移情况使用桑基图可视化。 首选.

空空如也

空空如也

1 2 3
收藏数 51
精华内容 20
关键字:

土地转移矩阵