第七章-遥感在土地利用覆盖中的应用

1第七章遥感应用1—遥感在土地利用/土地覆盖中的应用2土地利用和土地覆盖土地利用是人类在生产活动中为达到一定的经济效益、社会效益和生态效益，对土地资源的开发、经营、使用方式的总称。土地覆盖的定义有：1)“国际地圈与生物圈计划”（IGBP）和“全球环境变化人文计划”（HDP）：地球陆地表层和近地面层的自然状态，是自然过程和人类活动共同作用的结果。2)美国“全球环境变化委员会”（USSGCR）：覆盖着地球表面的植被及其它特质。3概念（续）从两者的定义可以看出，土地利用和土地覆盖既有一定的联系又有差别。土地利用重点是表示与土地相结合的人类活动而产生的不同利用方式。土地覆盖主要是表示地球表面存在的不同类型的覆盖特征，强调的是土地的表面形状。4概念（续）我国土地利用分类的目的是查清各县各种土地分类的面积、分布和利用状况，为国家制定国民经济计划，农业区划和规划，指导农业生产和服务与土地统计登记制度等管理工作。因此土地利用分类的主要依据是土地用途、土地经营方式、土地利用方式和土地覆盖特征等。土地覆盖只是土地利用分类的一个依据，但遥感图像最能够直接反映的是土地覆盖。大多数遥感土地应用使用土地覆盖/利用的概念。5分类标准全国土地遥感动态监测的分类标准。该标准基本遵循了全国农业区划委员会1984年颁发的《土地利用现状调查规程》的分类体系，并考虑遥感图像自身的特点。土地覆盖/利用类别分类定为以下8大类49二级类别（表）。编码在矢量图形中为多边形的属性码，在栅格图像上为栅格值。6分类标准（续）7示例－福建土地利用动态监测8数据福州市遥感影像：1990年，1996年，4景LANDSAT-TM土地利用图：长乐县1993年，1：10万1：10万，1：5万地形图农历、耕作、农事、等文字资料9光谱分析福州市主要地物光谱曲线050100150200250123457波段灰度城市村庄建设中耕地（有作物）耕地（无作物）园林地荒地河流浅海深海10变化信息提取配准时相①TM影像地理辅助数据时相②TM影像波段分析假彩色合成训练区选择样本分析与校正最大似然方法分类结果（Ｉ90和Ｉ97）行政界线清绘扫描或数字化矢量图形属性投影变换高斯克吕格投影矢量数据（V）地图代数计算动态变化图象（dI）dI和V配准迭加初始动态图统计数据人机交互目视纠正11变化信息提取（续）时相①遥感影像时相②遥感影像配　准裁　剪和高克投影配准主成分分析法分析图像特征选择反映变化信息明显的分量分　类初始动态分类（有变化）确定动态变化的具体内容（由什么变成什么）初始动态变化图详查图数字化形式配　准迭加12变化信息提取（续）动态信息分类法三种变化信息提取的评述：13动态信息提取利用公式B1*10+B2计算得到一幅新的图像，例如某点像元在90年为耕地（灰度值为3），97年变为居民地（灰度值为0），则动态图像上该像元的灰度值为30，表示耕地变为居民地。14成果图件制作和整饰初始动态变化图象≥５个像元编码野外实证，记录标注各图斑变更情况表格修正的动态图矢量境界图真彩色合成图图框公里网格指北针图例排版整饰正式结果分析报告正式表格提交最终结果迭加15精度检验方法鉴于混淆矩阵成为评价遥感图像分类结果的标准方法，拟根据实地对变化图斑的考察记录，选取不少于考察图斑60%的图斑，对比遥感得到的结果，形成混淆矩阵，进行精度评价。利用混淆矩阵进行精度评价时，偶然的一致往往会发生。而总体精度并不能排除这种偶然一致性。因此在混淆矩阵基础上，加上kappa统计来进行精度评价。Kappa统计的意义是：如果Kappa统计为0.7，则表示所用的分类方法比随机赋予各点某一类别的方法优越70%。16精度检验方法（续）因为变化精度与图斑的大小有关，一般而言图斑越大，精度越高，反则反之。按图斑的大小进行分级，不同级别图斑精度要求不同。在经费允许的情况下，使用GPS将实地考察的变化图斑输入计算机中，进行精度检验。17分类结果18动态结果19结果120结果2213提高分类精度的方法改进由于遥感信息存在固有的“同物异谱，异物同谱”，严重制约了基于光谱特征的统计分类精度的提高。而遥感专题分类的精度决定了遥感信息的实际应用价值。增加遥感图像的分类精度作出了大量的工作。大致可以分为以下几类：223提高分类精度的方法改进3.1传统方法的改进3.2分层分区3.3图像空间信息分类3.4多源辅助数据综合分类233.1传统最大似然法的改进FabioMaselli等认为传统遥感分类往往得不到满意的结果，部分原因是算法内在的缺陷。基于多维正态分布假设的最大似然法，统计上稳定而强健，但是缺乏灵活性，在复杂或非均质的情况下，这一假设往往难以满足而不能得到正确的面积估计。从而根据无参数方法，从图像灰度直方图中抽取先验概率，加入到MLC中，使分类结果明显改善（kappa0.3749-0.5860）。243.1传统最大似然法的改进（续）C.Conese等（1993）认为经过主成分分析的第一主分量包含了亮度信息，和地形的影响密切相关，通过对图像训练样本进行主成分分析，去除第一分量的影响来修订最大似然分类法，以减少地形的影响。这一修订的方法在进行地形破碎的地区分类时改善了分类的精度。25提高分类精度的方法改进3.1传统方法的改进3.2分层分区3.3图像空间信息分类3.4多源辅助数据综合分类263.2分层分区分层分区则在对不同的类别利用不同的信息子集。在山区植被分类中对图像进行四叉树分割，利用距离判据来决定是否需要更多的信息参与分类。这样每一步尽可能少地利用信息，只有当判断模糊时，才询问辅助信息，以降低时间开销。但基于四叉树的图像分割不具有地理学意义，当图像中方格状图斑多时，效率较高；如果地物复杂，四叉树划分的效率不高。273.2分层分区（续）在研究草场植被类型分类中，鉴于单纯的监督分类不能有效的区分不同景观单元的草场植被类型，采用了先对图像按照地貌类型分区在在各个分区上分别进行监督分类。在使用多景TM图像进行流域范围尺度上地形覆盖类型分类时，在图像镶嵌后，将图像分为三个生态区，进行分类，克服了单景分类不具有任何生态学意义的缺点。283.2分层分区（续）Loveland（1991）在利用NOAA/AVHRR图像进行美国土地利用分类制图时，先根据GIS数据库中气候、植被、环境等数据将全国划分为几个小的区域，然后再对各个小的区域进行非监督分类，减少了误分的像元数。JesusS.M.A等（1997）利用TM和SPOT图像派生出比值图像主成分分析图像等，对于各种图像组合进行了比较研究，肯定了分层分区的优点。293提高分类精度的方法改进3.1传统方法的改进3.2分层分区3.3图像空间信息分类3.4多源辅助数据综合分类303.3图像空间信息的分类遥感信息综合了地物种类、组成、坡度、坡向、阴影、太阳高度角等信息，而传统的统计分类仅仅使用了图像的单像元的光谱特性，忽略了像元之间的相互空间关系。仅仅依赖于光谱信息的分类结果在实际应用中是不成功的。遥感图像的空间信息在区分地物时具有非常重要的作用，是目视解译的重要解译标志之一。将图像空间信息用于计算机自动分类是自然而然的想法。313.3图像空间信息的分类（续）早在1980年，Switzer将邻域像元的平均值输入到最大似然法中，是计算机图像分类运用空间信息的开始。纹理（Texture）是灰度在平面空间上的变化，是遥感信息的重要内容。要利用图像空间信息提高分类精度，合理而有效地度量纹理至关重要。323.3图像空间信息的分类（续）图像纹理的表示方法分为两类，一是结构方法，表示基本原始模式在空间的重复，如砖墙。结构方式中排列的方式描述十分复杂。二是统计方法。由于地物的组成、空间分布的复杂性和多样性，反映在遥感图像上的千差万别的纹理不会象砖墙那样有规则不变的局部模式和简单的周期性的重复。遥感图像的纹理特征和重复性往往只具有统计学上的意义。因此统计纹理分析方法是遥感图像纹理分析的主要方法333.3图像空间信息的分类（续）无论从历史发展还是从当前进展来看，统计的方法仍然占主导地位（王润生，1995）。统计方法又分为空间域和频率域方法。空间域方法基于统计图像像元灰度级的分布情况，如利用直方图。频率域方法首先将图像变换到傅立叶变换的频率域中，然后抽取相应的特征量，做高频或低频的压缩或扩展后，再变换回空间域处理纹理。频率域方法是遥感光学图像处理中常采用的纹理分析方法。343.3图像空间信息的分类（续）图像模式识别领域对纹理的研究和应用比较深入，描述纹理的方法归纳起来主要有：共生矩阵，纹理能量测量，Markov随机场模型，分形布朗模型等。353.3图像空间信息的分类（续）纹理是图像处理相关领域的表达图像空间信息的一个重要度量，遥感应用中更多地使用结构（也称上下文，context）。研究人员根据其应用目的，提出了一些用于图像分类的结构方法。出现频率。出现频率F(i,j,v)定义为灰度值为v的像元在以i,j为中心。边长为l的窗口内出现的频率。对单波段来讲，v是一个值，多波段时v为向量。由图像生成出现频率表，表项的数目跟灰度级(n)和波段数(m)有关，最大数目为nm。出现频率表较空间统计度量具有两个优点。一是比任何一个空间域度量要包含更多的信息。事实上常用的空间域度量可以从出现频率中直接求出。二是出现频率表可以快速地生成，无须占据硬盘空间。363.3图像空间信息的分类（续）C.Ricotta等（1996）在进行森林研究时，对NDVI图像使用了计算简洁的纹理绝对差值算法来利用空间信息。应用表明纹理差值在对于植被和非植被类型的区分时，能够抓住两者之间纹理特征的本质。显然本文研究的类别只有两种，是否对更细类别分类时同样实用，需要实验。373.3图像空间信息的分类（续）利用小波变换，产生的小波子图来利用图像空间信息。在连同多波段图像输入KOHONEN自组织神经网络分类器后的分类结果表明，小波子图在区分细致地物特别是边缘检测的优势是传统分类方法所无法比拟的。383.3图像空间信息的分类（续）－评论空间信息的使用主要通过“窗口”来进行的。窗口过大，将别的类别纳入，往往产生误差边缘；窗口过小，常常引入不相关或相关小的像元，产生不可靠的纹理，并引入随机纹理。纹理在图像分类的作用严重受到“窗口”的影响，对于分类的精度提高不是简单的提高。由于缺乏对遥感图像纹理机理的深刻认识，成功的机遇并不很大，带有一定的随即性、偶然性和片面性。393.3图像空间信息的分类（续）－评论要真正地利用遥感图像的空间信息，还必须模拟目视解译过程，把判读人员的经验加入到影响分类过程中，建立基于知识的纹理和线形特征识别以及邻域分析方法（ArgialsD.P等，1990）。Paul(1997)在这方面做了探讨，初步显示了基于知识的图像空间信息方法在提高分类精度的潜力。403提高分类精度的方法改进3.1传统方法的改进3.2分层分区3.3图像空间信息分类3.4多源辅助数据综合分类413.4多源辅助数据综合分类遥感图像信息的原型自然环境的综合体是非常复杂的。遥感获得的信息并非是自然综合体信息的全部，而是自然综合体历史演进过程中瞬时的、能够在二维平面上表现的那一部分。这就决定了遥感信息的单纯数学、物理处理结果具有不确定性或多解性。为了提高遥感专题信息识别结果的正确性和可靠性，必须将遥感未带回的信息补充进去，即加入地学相关信息。然而遥感信息是主要的，地理相关信息是辅助性的。423.4.1辅助数据的类型①地形图及其它等值线图。由地形图得到数字地面模型（DTM），再由DTM计算出绝对高程、相对高程、坡度和坡向。这些信息是遥感图像辅助信息的最常用的。其他的如降雨等值线图、温度等值线图等，可以采用类似DTM的处理方法。433.4.1辅助数据的类型（续）②专题地图。专题地图是指与某一具体专题研究有关的图件。如土地利用现状图、土壤分布图、植被类型分布图、地质图、地貌图等。从专题图上可以提