地理信息系统基础与实践Geographical Information System Theory and Practice 第六章第六章 空间分析空间分析—— —— 叠加分析叠加分析 旅游与地理学院旅游与地理学院 陈世发陈世发2024/8/311韶关学院旅游与地理学院 陈世发 一、空间叠置概念一、空间叠置概念 空间叠置分析空间叠置分析(Spatial Overlay Analysis)是指在统一空间参照是指在统一空间参照系统条件下，每次将同一地区系统条件下，每次将同一地区两个地理对象的图层进行叠置，两个地理对象的图层进行叠置，以产生空间区域的多重属性特以产生空间区域的多重属性特征，或建立地理对象之间的空征，或建立地理对象之间的空间对应关系间对应关系空间叠置分析空间叠置分析目的：寻找和确定同时具有几种地目的：寻找和确定同时具有几种地目的：寻找和确定同时具有几种地目的：寻找和确定同时具有几种地理属性的地理要素的分布，或是按理属性的地理要素的分布，或是按理属性的地理要素的分布，或是按理属性的地理要素的分布，或是按照确定的地理指标，对叠加后产生照确定的地理指标，对叠加后产生照确定的地理指标，对叠加后产生照确定的地理指标，对叠加后产生的具有不同属性级的多边形进行分的具有不同属性级的多边形进行分的具有不同属性级的多边形进行分的具有不同属性级的多边形进行分类或分级类或分级类或分级类或分级 2024/8/312韶关学院旅游与地理学院 陈世发 二、叠置分析的分类二、叠置分析的分类§将将将将同同同同一一一一地地地地区区区区，，，，同同同同一一一一比比比比例例例例尺尺尺尺的的的的两两两两组组组组或或或或更更更更多多多多的的的的多多多多边边边边形形形形要要要要素素素素的的的的数数数数据据据据层层层层进进进进行行行行叠叠叠叠置置置置，，，，会会会会出出出出现现现现两两两两种种种种成成成成果果果果形形形形式式式式，，，，一一一一种种种种是是是是根根根根据据据据两两两两组组组组多多多多边边边边形形形形边边边边界界界界的的的的交交交交点点点点来来来来建建建建立立立立具具具具有有有有多多多多重重重重属属属属性性性性的的的的多多多多边边边边形形形形，，，，称称称称为为为为合合合合成成成成叠叠叠叠置置置置；；；；另另另另一一一一种种种种是进行多边形范围的属性特征的统计分析，称为是进行多边形范围的属性特征的统计分析，称为是进行多边形范围的属性特征的统计分析，称为是进行多边形范围的属性特征的统计分析，称为统计叠置统计叠置统计叠置统计叠置。
合成叠置得到一张新的叠置图，合成叠置得到一张新的叠置图，合成叠置得到一张新的叠置图，合成叠置得到一张新的叠置图，产生了许多新的多边形，每个多产生了许多新的多边形，每个多产生了许多新的多边形，每个多产生了许多新的多边形，每个多边形内都边形内都边形内都边形内都具有两种以上的属性具有两种以上的属性具有两种以上的属性具有两种以上的属性，，，，通过区域多重属性的模拟，寻找通过区域多重属性的模拟，寻找通过区域多重属性的模拟，寻找通过区域多重属性的模拟，寻找和确定同时具有几种地理属性的和确定同时具有几种地理属性的和确定同时具有几种地理属性的和确定同时具有几种地理属性的分布区域分布区域分布区域分布区域§ §统计叠置得到的是统计叠置得到的是统计叠置得到的是统计叠置得到的是新的属性数据新的属性数据新的属性数据新的属性数据，，，，得到的新属性值可能与原数据层的属得到的新属性值可能与原数据层的属得到的新属性值可能与原数据层的属得到的新属性值可能与原数据层的属性意义不一致，因为得到的结果可能性意义不一致，因为得到的结果可能性意义不一致，因为得到的结果可能性意义不一致，因为得到的结果可能是名称、分类值，也可能是分析后的是名称、分类值，也可能是分析后的是名称、分类值，也可能是分析后的是名称、分类值，也可能是分析后的结果。
结果2024/8/313韶关学院旅游与地理学院 陈世发 叠置分析的叠置分析的数学基础数学基础（空间逻辑运算）（空间逻辑运算）叠加过程往往是对空间信息和对应的属性信息作集合的叠加过程往往是对空间信息和对应的属性信息作集合的叠加过程往往是对空间信息和对应的属性信息作集合的叠加过程往往是对空间信息和对应的属性信息作集合的交、并、交、并、交、并、交、并、差、余差、余差、余差、余运算，也可再进一步对属性作其他的数学运算运算，也可再进一步对属性作其他的数学运算运算，也可再进一步对属性作其他的数学运算运算，也可再进一步对属性作其他的数学运算 为讨论方便将空间图层为讨论方便将空间图层A，，B，，C定义为定义为二值图象二值图象 3、空间逻辑差运算；、空间逻辑差运算；A - B =X X∈∈A 且且 X∈∈B 1、空间逻辑并（、空间逻辑并（或或）运算；）运算；A∪∪B =X X∈∈A 或或 X∈∈ B 2、空间逻辑交（、空间逻辑交（与与）运算；）运算；A∩B = X X∈∈A 且且 X∈∈B 4、空间包含；、空间包含；A B2024/8/314韶关学院旅游与地理学院 陈世发 þþ例：土壤厚度（大于例：土壤厚度（大于例：土壤厚度（大于例：土壤厚度（大于5050厘米）和土壤类型（红壤和其他类型）厘米）和土壤类型（红壤和其他类型）厘米）和土壤类型（红壤和其他类型）厘米）和土壤类型（红壤和其他类型）两个二值化图层，不同的逻辑运算结果如下：两个二值化图层，不同的逻辑运算结果如下：两个二值化图层，不同的逻辑运算结果如下：两个二值化图层，不同的逻辑运算结果如下：ððANDAND关系：结果是将土层厚度大于关系：结果是将土层厚度大于关系：结果是将土层厚度大于关系：结果是将土层厚度大于5050厘米，且土壤为红壤的厘米，且土壤为红壤的厘米，且土壤为红壤的厘米，且土壤为红壤的土壤单元显示出来；土壤单元显示出来；土壤单元显示出来；土壤单元显示出来；ððOROR关系：结果将土层厚度大于关系：结果将土层厚度大于关系：结果将土层厚度大于关系：结果将土层厚度大于5050厘米，或者土壤为红壤的土厘米，或者土壤为红壤的土厘米，或者土壤为红壤的土厘米，或者土壤为红壤的土壤单元显示出来；壤单元显示出来；壤单元显示出来；壤单元显示出来；ððXORXOR：：：：结果将土层厚度小于结果将土层厚度小于结果将土层厚度小于结果将土层厚度小于5050厘米，或者土壤不是红壤的土厘米，或者土壤不是红壤的土厘米，或者土壤不是红壤的土厘米，或者土壤不是红壤的土壤单元显示出来；壤单元显示出来；壤单元显示出来；壤单元显示出来；ððNOT:NOT:如结果是将土层厚度大于如结果是将土层厚度大于如结果是将土层厚度大于如结果是将土层厚度大于5050厘米，但土壤不是红壤的土厘米，但土壤不是红壤的土厘米，但土壤不是红壤的土厘米，但土壤不是红壤的土壤单元显示出来；壤单元显示出来；壤单元显示出来；壤单元显示出来；ü逻辑关系运算逻辑关系运算例例2024/8/315韶关学院旅游与地理学院 陈世发 栅格数据的算术运算 1)1)算术运算算术运算 指两层以上的对应网格值经加、减运算，而得到新的栅指两层以上的对应网格值经加、减运算，而得到新的栅格数据系统的方法。
格数据系统的方法 （（4 4）叠加分析的数学运算叠加法）叠加分析的数学运算叠加法2024/8/316韶关学院旅游与地理学院 陈世发 算术运算算术运算算术运算算术运算————以灰度以灰度以灰度以灰度( ( ( (亮度亮度亮度亮度) ) ) )赋值栅格为例（赋值栅格为例（赋值栅格为例（赋值栅格为例（乘除运算乘除运算乘除运算乘除运算））））1010505050501001002020100100100100200200×2＝＝（图像反差增大）（图像反差增大）（图像反差增大）（图像反差增大）阳阳阳阳光光光光砂岩砂岩砂岩砂岩页岩页岩页岩页岩花岗岩花岗岩花岗岩花岗岩阳坡阳坡阳坡阳坡阴坡阴坡阴坡阴坡波波波波 段段段段MSSMSS4 4MSSMSS5 5MSSMSS4 4/ / / /MSSMSS5 5阳阳阳阳 坡坡坡坡343451510.670.67阴阴阴阴 坡坡坡坡222234340.650.65砂岩对阳光的反射情况砂岩对阳光的反射情况砂岩对阳光的反射情况砂岩对阳光的反射情况乘法乘法乘法乘法案例案例案例案例除法案例除法案例除法案例除法案例0.670.670.650.652024/8/317韶关学院旅游与地理学院 陈世发 无数据置换无数据置换无数据置换无数据置换①①①①算术运算算术运算算术运算算术运算————赋值栅格的算术条件赋值栅格的算术条件赋值栅格的算术条件赋值栅格的算术条件叠加叠加叠加叠加比较后取最大值比较后取最大值比较后取最大值比较后取最大值比较后取最小值比较后取最小值比较后取最小值比较后取最小值2024/8/318韶关学院旅游与地理学院 陈世发 ①①①①算术运算算术运算算术运算算术运算————赋值栅格的算术条件赋值栅格的算术条件赋值栅格的算术条件赋值栅格的算术条件叠加叠加叠加叠加数据置换数据置换数据置换数据置换2024/8/319韶关学院旅游与地理学院 陈世发 2)2)函数运算：指两个以上层面的栅格数据系统以某种函数运算：指两个以上层面的栅格数据系统以某种函数运算：指两个以上层面的栅格数据系统以某种函数运算：指两个以上层面的栅格数据系统以某种函数关函数关函数关函数关系系系系作为复合分析的依据进行逐网格运算，从而得到新的栅作为复合分析的依据进行逐网格运算，从而得到新的栅作为复合分析的依据进行逐网格运算，从而得到新的栅作为复合分析的依据进行逐网格运算，从而得到新的栅格数据系统的过程。
格数据系统的过程格数据系统的过程格数据系统的过程 这种复合叠置分析方法被广泛地应用到地学综合分析、环这种复合叠置分析方法被广泛地应用到地学综合分析、环这种复合叠置分析方法被广泛地应用到地学综合分析、环这种复合叠置分析方法被广泛地应用到地学综合分析、环境质量评价、遥感数字图像处理等领域中境质量评价、遥感数字图像处理等领域中境质量评价、遥感数字图像处理等领域中境质量评价、遥感数字图像处理等领域中 利用土壤侵蚀通用方程式利用土壤侵蚀通用方程式利用土壤侵蚀通用方程式利用土壤侵蚀通用方程式计算土壤侵蚀量时，就可利计算土壤侵蚀量时，就可利计算土壤侵蚀量时，就可利计算土壤侵蚀量时，就可利用多层面栅格数据的函数运用多层面栅格数据的函数运用多层面栅格数据的函数运用多层面栅格数据的函数运算复合分析法进行自动处理算复合分析法进行自动处理算复合分析法进行自动处理算复合分析法进行自动处理一个地区土壤侵蚀量的大小一个地区土壤侵蚀量的大小一个地区土壤侵蚀量的大小一个地区土壤侵蚀量的大小是降雨是降雨是降雨是降雨(R)(R)、植被覆度、植被覆度、植被覆度、植被覆度(C)(C)、、、、坡度坡度坡度坡度(S)(S)、坡长（、坡长（、坡长（、坡长（L L）、土壤）、土壤）、土壤）、土壤抗蚀性抗蚀性抗蚀性抗蚀性(SR)(SR)等因素的函数等因素的函数等因素的函数等因素的函数 在在Arc View中，使用中，使用Map Calculator可以很方便地实现可以很方便地实现栅格图层的复合（叠加）运算栅格图层的复合（叠加）运算土壤侵蚀多因子函数运算复合分析示意图 函函数数运运算算露点露点2024/8/3110韶关学院旅游与地理学院 陈世发 n n1 1、视觉信息叠加：、视觉信息叠加：、视觉信息叠加：、视觉信息叠加：将不同专题的内容叠加显示在结果图将不同专题的内容叠加显示在结果图将不同专题的内容叠加显示在结果图将不同专题的内容叠加显示在结果图件上件上件上件上 ，，，， 视觉信息叠加之后，参加叠加的平面之间没发生任视觉信息叠加之后，参加叠加的平面之间没发生任视觉信息叠加之后，参加叠加的平面之间没发生任视觉信息叠加之后，参加叠加的平面之间没发生任何逻辑关系，仍保留原来的数据结构。
何逻辑关系，仍保留原来的数据结构何逻辑关系，仍保留原来的数据结构何逻辑关系，仍保留原来的数据结构 面状图、线状图和点状图之间的叠加；面状图、线状图和点状图之间的叠加；面状图、线状图和点状图之间的叠加；面状图、线状图和点状图之间的叠加；üü面状图区域边界之间或一个面状图与其他专题区域边界之间面状图区域边界之间或一个面状图与其他专题区域边界之间面状图区域边界之间或一个面状图与其他专题区域边界之间面状图区域边界之间或一个面状图与其他专题区域边界之间的叠加；的叠加；的叠加；的叠加；üü遥感影像与专题地图的叠加；遥感影像与专题地图的叠加；遥感影像与专题地图的叠加；遥感影像与专题地图的叠加；üü专题地图与数字高程模型叠加显示立体专题图；专题地图与数字高程模型叠加显示立体专题图；专题地图与数字高程模型叠加显示立体专题图；专题地图与数字高程模型叠加显示立体专题图；üü遥感影像与遥感影像与遥感影像与遥感影像与DEMDEM叠加生成真三维地物景观叠加生成真三维地物景观叠加生成真三维地物景观叠加生成真三维地物景观 叠加分析分为以下五类：视觉信息叠加；点与多边形叠加分析分为以下五类：视觉信息叠加；点与多边形叠加分析分为以下五类：视觉信息叠加；点与多边形叠加分析分为以下五类：视觉信息叠加；点与多边形叠加；叠加；叠加；叠加； 线与多边形叠加线与多边形叠加线与多边形叠加线与多边形叠加 ；多边形叠加；；多边形叠加；；多边形叠加；；多边形叠加； 栅格图层叠栅格图层叠栅格图层叠栅格图层叠加加加加 面状图、面状图、线状图和线状图和点状图之点状图之间的叠加间的叠加面状图与专题区域边界面状图与专题区域边界之间的叠加之间的叠加遥感影像与专题地图的叠加遥感影像与专题地图的叠加遥感影像与专题地图的叠加遥感影像与专题地图的叠加专题地图与数字高程模型叠加专题地图与数字高程模型叠加遥感影像与遥感影像与DTM叠加叠加2024/8/3111韶关学院旅游与地理学院 陈世发 2、点与多边形的叠加üü叠加图层：叠加图层：叠加图层：叠加图层：将一个含有点的图层（目标图层）叠加在另一个含有将一个含有点的图层（目标图层）叠加在另一个含有将一个含有点的图层（目标图层）叠加在另一个含有将一个含有点的图层（目标图层）叠加在另一个含有多边形的图层（操作图层）上，以确定每个点落在哪个区域内。
多边形的图层（操作图层）上，以确定每个点落在哪个区域内多边形的图层（操作图层）上，以确定每个点落在哪个区域内多边形的图层（操作图层）上，以确定每个点落在哪个区域内 üü例如：例如：例如：例如：一个中国政区图（多边形）和一个全国矿产分布图（点），一个中国政区图（多边形）和一个全国矿产分布图（点），一个中国政区图（多边形）和一个全国矿产分布图（点），一个中国政区图（多边形）和一个全国矿产分布图（点），二者经叠加分析后，并且将政区图多边形有关的属性信息加到矿二者经叠加分析后，并且将政区图多边形有关的属性信息加到矿二者经叠加分析后，并且将政区图多边形有关的属性信息加到矿二者经叠加分析后，并且将政区图多边形有关的属性信息加到矿产的属性数据表中，然后通过属性查询，可以查询指定省有多少产的属性数据表中，然后通过属性查询，可以查询指定省有多少产的属性数据表中，然后通过属性查询，可以查询指定省有多少产的属性数据表中，然后通过属性查询，可以查询指定省有多少种矿产，产量有多少；而且可以查询指定类型的矿产在哪些省里种矿产，产量有多少；而且可以查询指定类型的矿产在哪些省里种矿产，产量有多少；而且可以查询指定类型的矿产在哪些省里种矿产，产量有多少；而且可以查询指定类型的矿产在哪些省里有分布等信息。
有分布等信息有分布等信息有分布等信息 ü例如将水井与规划区图层相叠置，可确定每口井所属的规划区范例如将水井与规划区图层相叠置，可确定每口井所属的规划区范围1 *3 *4 *2A DB C点 属性1234多边形 属性1 属性2ABCD点 多边形 点属性 面属性1 面属性21 A2 D3 C4 B2024/8/3112韶关学院旅游与地理学院 陈世发 3、线与多边形的叠加、线与多边形的叠加üü叠加图层叠加图层叠加图层叠加图层：将线的图层（目标图层）叠加在多边形的图层：将线的图层（目标图层）叠加在多边形的图层：将线的图层（目标图层）叠加在多边形的图层：将线的图层（目标图层）叠加在多边形的图层（操作图层）上，以确定一条线落在哪个多边形内操作图层）上，以确定一条线落在哪个多边形内操作图层）上，以确定一条线落在哪个多边形内操作图层）上，以确定一条线落在哪个多边形内。
üü叠加原理叠加原理叠加原理叠加原理：与前面不同的是，往往一个线目标跨越多个多边：与前面不同的是，往往一个线目标跨越多个多边：与前面不同的是，往往一个线目标跨越多个多边：与前面不同的是，往往一个线目标跨越多个多边形，这时需要先进行线与多边形的求交，并将线目标进行切形，这时需要先进行线与多边形的求交，并将线目标进行切形，这时需要先进行线与多边形的求交，并将线目标进行切形，这时需要先进行线与多边形的求交，并将线目标进行切割，形成一个新的空间目标（新的线目标）的结果集割，形成一个新的空间目标（新的线目标）的结果集割，形成一个新的空间目标（新的线目标）的结果集割，形成一个新的空间目标（新的线目标）的结果集üü例如：例如：例如：例如：当线状图层为河流，叠加的结果是多边形将穿过它的当线状图层为河流，叠加的结果是多边形将穿过它的当线状图层为河流，叠加的结果是多边形将穿过它的当线状图层为河流，叠加的结果是多边形将穿过它的所有河流打断成弧段，可以查询任意多边形内的河流长度，所有河流打断成弧段，可以查询任意多边形内的河流长度，所有河流打断成弧段，可以查询任意多边形内的河流长度，所有河流打断成弧段，可以查询任意多边形内的河流长度，进而计算它的河流密度等；如果线状图层为道路网，叠加的进而计算它的河流密度等；如果线状图层为道路网，叠加的进而计算它的河流密度等；如果线状图层为道路网，叠加的进而计算它的河流密度等；如果线状图层为道路网，叠加的结果可以得到每个多边形内的道路网密度，内部的交通流量，结果可以得到每个多边形内的道路网密度，内部的交通流量，结果可以得到每个多边形内的道路网密度，内部的交通流量，结果可以得到每个多边形内的道路网密度，内部的交通流量，进入、离开各个多边形的交通量，相邻多边形之间的相互交进入、离开各个多边形的交通量，相邻多边形之间的相互交进入、离开各个多边形的交通量，相邻多边形之间的相互交进入、离开各个多边形的交通量，相邻多边形之间的相互交通量。
通量 2024/8/3113韶关学院旅游与地理学院 陈世发 例如：确定某一行政区内各种等级道路的里程数例如：确定某一行政区内各种等级道路的里程数线与多边形的叠合的目的是确定某一线状图层上的弧段落与多边形的叠合的目的是确定某一线状图层上的弧段落在另一多变性图层上的哪个多边形内，以便为图层的每条弧段另一多变性图层上的哪个多边形内，以便为图层的每条弧段建立新的属性建立新的属性路网密度？路网密度？路网密度？路网密度？————叠加！叠加！叠加！叠加！•各种道路长=830KM•区域面积=2797KM2•结果？2024/8/3114韶关学院旅游与地理学院 陈世发 4、多边形与多边形的叠加、多边形与多边形的叠加 多边形与多边形的叠置是指将两个不同图层的多边形要素多边形与多边形的叠置是指将两个不同图层的多边形要素相叠加，根据两组多边形的交点来建立多重属性的多边形或相叠加，根据两组多边形的交点来建立多重属性的多边形或进行多边形范围内的属性特征的统计分析进行多边形范围内的属性特征的统计分析 多边形叠置过程分几何求交过程和属性确定过程，算法的多边形叠置过程分几何求交过程和属性确定过程，算法的核心是多边形求交。
核心是多边形求交①①对两个多边形进行边界求交和弧段分割运算，并以新弧段对两个多边形进行边界求交和弧段分割运算，并以新弧段为单位重建拓扑关系；为单位重建拓扑关系；②②判断重建多边形落在原始多边形层的哪个多边形内，从而判断重建多边形落在原始多边形层的哪个多边形内，从而建立新叠置多边形与原始多边形的关系，并抽取属性建立新叠置多边形与原始多边形的关系，并抽取属性2024/8/3115韶关学院旅游与地理学院 陈世发 多边形的不同叠加方式多边形的不同叠加方式取本底多边形取上叠多边形两个多边形相交？两个多边形是包含关系？多边形求交，重建拓扑关系属性抽取，建立关联连接上叠多边形结束？本底多边形结束？ynnn建立包含关系属性抽取yyy结束多多边边形形的的叠叠置置流流程程图图2024/8/3116韶关学院旅游与地理学院 陈世发 新多边形的属性多边形之间的叠多边形之间的叠置置IDID属性属性属性属性101101A AIDID属性属性属性属性1 1X X2 2Y Y3 3Z Z新多边新多边新多边新多边形形形形IDID层层层层1 1多边多边多边多边形属性形属性形属性形属性层层层层2 2多边形多边形多边形多边形属性属性属性属性1 10 0A A2 20 0X X3 3A AX X4 40 0X X5 5A A0 06 6A A0 07 7A AY Y8 8A AZ Z9 90 0Y Y10100 0Z Z11110 0Y Y1212Z Z0 0AXYZ37910145611 2层层1层层2新层新层12坡度图层坡度图层坡度图层坡度图层肥力图层肥力图层肥力图层肥力图层叠加？叠加？叠加？叠加？2024/8/3117韶关学院旅游与地理学院 陈世发 多边形叠多边形叠加产生碎加产生碎屑多边形屑多边形由于矢量有限精度原因，几何对象由于矢量有限精度原因，几何对象不可能完全匹配，出现碎屑多边形，不可能完全匹配，出现碎屑多边形，通常设模糊容限消除通常设模糊容限消除2024/8/3118韶关学院旅游与地理学院 陈世发 5、基于栅格数据的、基于栅格数据的叠置叠置分析分析 特点：栅格数据的特点：栅格数据的叠置叠置算法，虽然数据存贮量比较大，但运算法，虽然数据存贮量比较大，但运算过程比较简单。
算过程比较简单变换方法：变换方法：（（1）点变换）点变换（（2）区域变换方法）区域变换方法（（3）邻域变换方法）邻域变换方法基于栅格数据的基于栅格数据的叠置叠置分析是参与分析的两个图层的要素分析是参与分析的两个图层的要素均为栅格数据均为栅格数据栅格叠置的作用包括以下几种：栅格叠置的作用包括以下几种：栅格叠置的作用包括以下几种：栅格叠置的作用包括以下几种： ①①①① 类型叠置；类型叠置；类型叠置；类型叠置；②②②② 数量统计数量统计数量统计数量统计 ③③③③ 动态分析动态分析动态分析动态分析 ④④④④ 益本分析益本分析益本分析益本分析 ⑤⑤⑤⑤ 几何提取几何提取几何提取几何提取 2024/8/3119韶关学院旅游与地理学院 陈世发 （（1）点变换）点变换 点变换只依据参与叠置图层相应点的属性值进行新的运算，点变换只依据参与叠置图层相应点的属性值进行新的运算，既与各图层的邻域点的属性无关，也不受区域内一般特征的既与各图层的邻域点的属性无关，也不受区域内一般特征的影响运算方法包括：算术运算，指数运算，三角函数运算运算方法包括：算术运算，指数运算，三角函数运算等等2 2 4 4 4 4 6 63 3 3 3 4 4 6 63 3 3 3 7 7 7 72 2 2 25 5 4 42 2 3 36 6 4 43 3 3 36 6 6 6-=0 0-2-21 1-2-2-1-10 02 2-2-20 00 0-1-1-1-1点变换示意点变换示意结果特征：运算后得到的新属性值可能与原图层的属性意结果特征：运算后得到的新属性值可能与原图层的属性意义完全不同。
义完全不同2024/8/3120韶关学院旅游与地理学院 陈世发 (2)区域变换方法区域变换方法 指在计算新图层相应的属性值时，不仅与原图层对应的指在计算新图层相应的属性值时，不仅与原图层对应的栅格的属性值有关，而且要顾及原图层所在区域的集合特栅格的属性值有关，而且要顾及原图层所在区域的集合特征（区域长度、面积、周长等）征（区域长度、面积、周长等）3)邻域变换方法邻域变换方法 指在计算新图层相应的属性值时，不仅考虑原图层对应指在计算新图层相应的属性值时，不仅考虑原图层对应的栅格及其属性，而且还应顾及与该栅格相关联的邻域或的栅格及其属性，而且还应顾及与该栅格相关联的邻域或者影响半径内的栅格属性值的影响者影响半径内的栅格属性值的影响2024/8/3121韶关学院旅游与地理学院 陈世发 叠加分析与视觉信息叠加的区别叠加分析与视觉信息叠加的区别§ §叠叠叠叠置置置置分分分分析析析析是是是是把把把把参参参参与与与与叠叠叠叠加加加加的的的的各各各各数数数数据据据据层层层层经经经经过过过过算算算算术术术术的的的的、、、、几几几几何何何何的的的的、、、、逻逻逻逻辑辑辑辑的的的的运运运运算算算算生生生生成成成成新新新新的的的的数数数数据据据据层层层层；；；；而而而而视视视视觉觉觉觉信信信信息息息息叠叠叠叠加加加加仅仅仅仅是是是是把把把把不不不不同同同同的的的的数数数数据据据据层层层层同同同同时时时时显显显显示示示示到到到到一一一一个个个个屏屏屏屏幕幕幕幕上上上上或或或或绘绘绘绘制制制制在在在在一一一一幅幅幅幅图图图图形形形形里里里里，，，，各各各各数数数数据据据据层层层层间间间间不不不不存存存存在在在在逻逻逻逻辑辑辑辑关关关关系系系系，，，，没没没没有有有有生生生生成成成成新新新新的的的的数数数数据据据据文文文文件件件件，，，，在在在在同同同同一区域内，地图投影与比例尺相同便可进行视觉信息复合。
一区域内，地图投影与比例尺相同便可进行视觉信息复合一区域内，地图投影与比例尺相同便可进行视觉信息复合一区域内，地图投影与比例尺相同便可进行视觉信息复合叠置分析和空间集合分析的区别叠置分析和空间集合分析的区别§ §叠置分析和空间集合分析虽然都可以在基于矢量的地理信叠置分析和空间集合分析虽然都可以在基于矢量的地理信叠置分析和空间集合分析虽然都可以在基于矢量的地理信叠置分析和空间集合分析虽然都可以在基于矢量的地理信息系统和基于栅格的地理信息系统中完成，但由于基于栅息系统和基于栅格的地理信息系统中完成，但由于基于栅息系统和基于栅格的地理信息系统中完成，但由于基于栅息系统和基于栅格的地理信息系统中完成，但由于基于栅格的系统的叠加分析要容易和快得多，而矢量数据在叠置格的系统的叠加分析要容易和快得多，而矢量数据在叠置格的系统的叠加分析要容易和快得多，而矢量数据在叠置格的系统的叠加分析要容易和快得多，而矢量数据在叠置操作时容易产生多边形误差且边界搜索也很困难，运算速操作时容易产生多边形误差且边界搜索也很困难，运算速操作时容易产生多边形误差且边界搜索也很困难，运算速操作时容易产生多边形误差且边界搜索也很困难，运算速度又慢等原因，所以一般用栅格数据进行空间集合分析。
度又慢等原因，所以一般用栅格数据进行空间集合分析度又慢等原因，所以一般用栅格数据进行空间集合分析度又慢等原因，所以一般用栅格数据进行空间集合分析2024/8/3122韶关学院旅游与地理学院 陈世发 101010106060606010101010202060603030303030306060101010106060606010101010303060603030303030306060-=0 00 00 00 00 00 010100 00 00 00 00 0空间叠加分析实例空间叠加分析实例1：土地利用变化区域探测：土地利用变化区域探测80年遥感影像年遥感影像90年遥感影像年遥感影像点变换后影像点变换后影像点变换后影像分析点变换后影像分析通过通过80和和90年两期影像的相减运算后得到变年两期影像的相减运算后得到变换影像，如果：换影像，如果：变换影像值变换影像值 ＝＝ 0；说明该区未发生变化；说明该区未发生变化变换影像值变换影像值 ≠ 0；说明该区已发生变化；说明该区已发生变化10 耕地耕地20 居民点居民点30 水域水域40 草地草地50 未利用地未利用地60 林地林地Legend注意：注意：此处的遥感影像可以是此处的遥感影像可以是分类结果分类结果，也可以，也可以是是原始的遥感影像原始的遥感影像。
在一般应用中，多使用原始在一般应用中，多使用原始的遥感影像，可提高变化探测速度的遥感影像，可提高变化探测速度2024/8/3123韶关学院旅游与地理学院 陈世发 ü信息获取信息获取ü研究方法研究方法ü变化过程分析变化过程分析 ü机制分析机制分析 实例实例2 2：：土地利用变化分析土地利用变化分析更新土地利用图社经统计数据动态变化研究变化信息的解译与制图解译与制图变化数据管理原始遥感信息专题图件数据预处理空间操作与分析图像预处理影像分类变化信息探测结果比较变化信息分类高分辨率图像融合 土地利用/覆盖动态变化研究框架2024/8/3124韶关学院旅游与地理学院 陈世发 信信息息获获取取1999、2003和2007年遥感影像图2024/8/3125韶关学院旅游与地理学院 陈世发 面积统计面积统计年份年份 类型类型耕地耕地园地园地林地林地城镇村城镇村与工与工矿用矿用地地交通用交通用地地水域水域未利用未利用地地19991135.37258.32678.22192.1230.276.4722.2420031112.65266.32580.97194.31105.6286.4946.462007976.25319.772684.80197.20105.4686.3823.082024/8/3126韶关学院旅游与地理学院 陈世发 土地类型转移矩阵土地类型转移矩阵:根据两个不同时间（根据两个不同时间（t和和t+△△t））的土地利用图计算从一种类型的土地利用图计算从一种类型到另一种类型的转换概率，来分析土地利用变化过程。
它依靠到另一种类型的转换概率，来分析土地利用变化过程它依靠GIS技术，将两个不同时间的土地利用图进行栅格化处理，计技术，将两个不同时间的土地利用图进行栅格化处理，计算算t时刻上时刻上A类有多少格网点转换成在类有多少格网点转换成在t+△△t时刻的时刻的B，，C，，D等等类型，转换点数占该类型总数的百分比可称为转移概率类型，转换点数占该类型总数的百分比可称为转移概率ü研究方法研究方法使用景观生态学方法，计算土地类型转移矩阵使用景观生态学方法，计算土地类型转移矩阵 12┄nai+1a11a12┄a1na1+2a21a22┄a2na2+┆┆┆┆ Nan1a32┄annan+a+ja+1a+2┄a+na注：表中注：表中 ai+ 代表土地类型代表土地类型 I 在转化前的总量；在转化前的总量；a+j 代表代表土地类型土地类型 I在转化后的总量；在转化后的总量；a 是土地总面积是土地总面积 类型转移类型转移距阵示例距阵示例 2024/8/3127韶关学院旅游与地理学院 陈世发 朱溪河小流域朱溪河小流域1999-2007土地利用类型转移矩阵土地利用类型转移矩阵 hm2土地利用转移矩土地利用转移矩阵阵耕地耕地园地园地林地林地城镇村及工城镇村及工矿用地矿用地交通交通用地用地水域水域未利未利用地用地1999年年合计合计1999年转年转移量移量耕地耕地624.5159.42295.474.0245.5129.376.891135.37510.86园地园地37.73161.5343.2811.846.691.340.12258.396.77林地林地245.0573.12285.3131.4129.969.210.162678.22392.91城镇村及工矿用城镇村及工矿用地地36.0417.2754.8470.2813.280.330.05192.12121.84交通用地交通用地6.563.230.4517.68.010.42030.221.19水域水域22.650.025.401.242.82076.4733.65未利用地未利用地6.1744.22.231.982.157.5122.2414.732007年合计年合计976.25319.572684.68197.2105.4686.3823.084392.92--2007年转移量年转移量 351.74158.04399.37126.9297.4543.5615.57--1191.952024/8/3128韶关学院旅游与地理学院 陈世发 n n分析目的和评价准则：分析目的和评价准则：分析目的和评价准则：分析目的和评价准则：n n1 1）））） 估计住宅用地被洪水淹设而造成的损失；估计住宅用地被洪水淹设而造成的损失；估计住宅用地被洪水淹设而造成的损失；估计住宅用地被洪水淹设而造成的损失；n n2 2）））） 洪水水位的相对高程为洪水水位的相对高程为洪水水位的相对高程为洪水水位的相对高程为500500米；米；米；米；n n3 3）））） 损失的大小和居民的财产、地基的稳定性有关；损失的大小和居民的财产、地基的稳定性有关；损失的大小和居民的财产、地基的稳定性有关；损失的大小和居民的财产、地基的稳定性有关；应用实例应用实例3：洪水淹没损失分析：洪水淹没损失分析获取数据：获取数据：获取数据：获取数据：1 1）））） 数字化的地块数字化的地块数字化的地块数字化的地块多边形地图多边形地图多边形地图多边形地图。
每个地块均有土地使用、可遭每个地块均有土地使用、可遭每个地块均有土地使用、可遭每个地块均有土地使用、可遭损失的财产状况（简称估计财产）、不同地基类型等属性损失的财产状况（简称估计财产）、不同地基类型等属性损失的财产状况（简称估计财产）、不同地基类型等属性损失的财产状况（简称估计财产）、不同地基类型等属性2 2）））） 地块多边形属性表中有地均财产这一项，地均财产地块多边形属性表中有地均财产这一项，地均财产地块多边形属性表中有地均财产这一项，地均财产地块多边形属性表中有地均财产这一项，地均财产= =估计估计估计估计财产财产财产财产/ /地块面积地块面积地块面积地块面积3 3）））） 对每一类地基，可估计其稳定性，并估计房屋倒坍的可能对每一类地基，可估计其稳定性，并估计房屋倒坍的可能对每一类地基，可估计其稳定性，并估计房屋倒坍的可能对每一类地基，可估计其稳定性，并估计房屋倒坍的可能性，称损失系数（见表）性，称损失系数（见表）性，称损失系数（见表）性，称损失系数（见表）4 4）））） 数字化的等高线地形图数字化的等高线地形图数字化的等高线地形图数字化的等高线地形图2024/8/3129韶关学院旅游与地理学院 陈世发 地块多边形图地块多边形图地块属性表地块属性表地形等高线及其组成的多边形地形高程表2024/8/3130韶关学院旅游与地理学院 陈世发 空间分析过程空间分析过程n n将将将将地地地地块块块块多多多多边边边边形形形形和和和和高高高高程程程程多多多多边边边边形形形形叠叠叠叠合合合合，，，，产产产产生生生生地地地地块块块块——高高高高程程程程多多多多边边边边形形形形地图和对应的属性表。
地图和对应的属性表地图和对应的属性表地图和对应的属性表n n在在在在地地地地块块块块——高高高高程程程程属属属属性性性性表表表表中中中中选选选选择择择择高高高高程程程程小小小小于于于于等等等等于于于于500500，，，，土土土土地地地地使使使使用用用用性性性性质质质质为为为为住住住住宅宅宅宅的的的的记记记记录录录录和和和和地地地地基基基基——损损损损失失失失系系系系数数数数对对对对照照照照表表表表连连连连接接接接，，，，获获获获得得得得新新新新的的的的地地地地块块块块——高高高高程程程程属属属属性性性性表表表表估估估估计计计计损损损损失失失失= =面面面面积积积积× ×地地地地均均均均财财财财产产产产× ×损损损损失失失失系系系系数n n从从从从表表表表可可可可知知知知，，，，当当当当洪洪洪洪水水水水淹淹淹淹没没没没了了了了500500米米米米以以以以下下下下的的的的地地地地区区区区时时时时，，，，每每每每个个个个地地地地块块块块财产的大致损失状况财产的大致损失状况财产的大致损失状况财产的大致损失状况n n对对对对地地地地块块块块————高高高高程程程程图图图图按按按按对对对对应应应应属属属属性性性性进进进进行行行行分分分分类类类类，，，，得得得得到到到到洪洪洪洪水水水水淹淹淹淹没没没没损损损损失分布图。
失分布图失分布图失分布图2024/8/3131韶关学院旅游与地理学院 陈世发 叠合后地块-高程属性表的数据项（内容略）地基类型地基类型——损失系数对照表损失系数对照表地块和地形叠合后的多边形地块和地形叠合后的多边形 洪水淹没损失分布洪水淹没损失分布 l将地块多边形和高程多边形叠合，产生地块—高程多边形地图和对应的属性表l估计损失=面积×地均财产×损失系数损失估计表 l从表可知，当洪水淹没了500米以下的地区时，每个地块财产的大致损失状况分析结论表 2024/8/3132韶关学院旅游与地理学院 陈世发 思考与讨论思考与讨论n n1 1、什么是空间叠置分析？、什么是空间叠置分析？、什么是空间叠置分析？、什么是空间叠置分析？GISGIS叠加分析可以分哪些类别？叠加分析可以分哪些类别？叠加分析可以分哪些类别？叠加分析可以分哪些类别？n n2 2、多边形叠加分析是如何实现的？、多边形叠加分析是如何实现的？、多边形叠加分析是如何实现的？、多边形叠加分析是如何实现的？n n3 3、在一个区域内，欲为一公园选址，要求在人口密度较、在一个区域内，欲为一公园选址，要求在人口密度较、在一个区域内，欲为一公园选址，要求在人口密度较、在一个区域内，欲为一公园选址，要求在人口密度较大的地区，且距离已有公园较远，在主干路旁，河流流经，大的地区，且距离已有公园较远，在主干路旁，河流流经，大的地区，且距离已有公园较远，在主干路旁，河流流经，大的地区，且距离已有公园较远，在主干路旁，河流流经，有林地分布。
需要的数据有哪些？怎样制作选址图？说明有林地分布需要的数据有哪些？怎样制作选址图？说明有林地分布需要的数据有哪些？怎样制作选址图？说明有林地分布需要的数据有哪些？怎样制作选址图？说明其步骤，并画出操作流程图其步骤，并画出操作流程图其步骤，并画出操作流程图其步骤，并画出操作流程图 ？？？？n n4 4、由于城市的发展，韶关需要增加一个垃圾填埋场，要、由于城市的发展，韶关需要增加一个垃圾填埋场，要、由于城市的发展，韶关需要增加一个垃圾填埋场，要、由于城市的发展，韶关需要增加一个垃圾填埋场，要求该填埋场避开良田，交通便捷，地形坡度小于求该填埋场避开良田，交通便捷，地形坡度小于求该填埋场避开良田，交通便捷，地形坡度小于求该填埋场避开良田，交通便捷，地形坡度小于5°5°，地质，地质，地质，地质条件安全，并避开居民区，试以条件安全，并避开居民区，试以条件安全，并避开居民区，试以条件安全，并避开居民区，试以GISGIS方法，用框图表示其方法，用框图表示其方法，用框图表示其方法，用框图表示其运行过程，并说明需要的相关数据源、数据采集与组织方运行过程，并说明需要的相关数据源、数据采集与组织方运行过程，并说明需要的相关数据源、数据采集与组织方运行过程，并说明需要的相关数据源、数据采集与组织方法，说出空间分析过程？法，说出空间分析过程？法，说出空间分析过程？法，说出空间分析过程？2024/8/3133韶关学院旅游与地理学院 陈世发 。