土地利用覆盖变化地信息提取
土地利用覆盖（LUCC）变化遥感解译.
⼟地利⽤覆盖（LUCC）变化遥感解译.⼟地利⽤/覆盖（LUCC）变化遥感解译1、资料的准备2、遥感影像的处理3、遥感影像的解译（初步解译）（1）解译的指标体系解译的指标⽤⼟地利⽤/覆盖分类系统，采⽤全国⼆级分类系统：⼀级分为6类，主要根据⼟地的⾃然⽣态和利⽤属性；⼆级分为25个类型，主要根据⼟地经营特点、利⽤⽅式和覆盖特征；耕地根据地形特征进⾏了三级划分，即进⼀步划分为平原、丘陵、⼭区和坡度⼤于25度的耕地。
分类系统及含义为：①耕地：指种植农作物的⼟地，包括熟耕地、新开荒地、休闲地、轮歇地、草⽥轮作地；以种植农作物为主的农果、农桑、农林⽤地；耕种三年以上的滩地和海涂。
11、⽔⽥：指有⽔源保证和灌溉设施，在⼀般年景能正常灌溉，⽤以种植⽔稻、莲藕等⽔⽣农作物的耕地，包括实⾏⽔稻和旱地作物轮种的耕地。
12、旱地：指⽆灌溉⽔源及设施，靠天然降⽔⽣长作物的耕地；有⽔源和浇灌设施，在⼀般年景下能正常灌溉的旱作物耕地；以种菜为主的耕地；正常轮作的休闲地和轮歇地。
②林地：指⽣长乔⽊、灌⽊、⽵类以及沿海红树林地等林业⽤地。
21、有林地：指郁闭度>30%的天然林和⼈⼯林。
包括⽤材林、经济林、防护林等成⽚林地。
22、灌⽊林地：指郁闭度>40%、⾼度在2⽶以下的矮林地和灌丛林地。
23、疏林地：指郁闭度为10-30%的稀疏林地。
24、其他林地：指未成林造林地、迹地、苗圃及各类园地（果园、桑园、茶园、热作林园等）。
③草地：指以⽣长草本植物为主、覆盖度在5%以上的各类草地，包括以牧为主的灌丛草地和郁闭度在10%以下的疏林草地。
31、⾼覆盖度草地：指覆盖度>50%的天然草地、改良草地和割草地。
此类草地⼀般⽔分条件较好，草被⽣长茂密。
32、中覆盖度草地：指覆盖度在20-50%的天然草地和改良草地，此类草地⼀般⽔分不⾜，草被较稀疏。
33、低覆盖度草地：指覆盖度在5-20%的天然草地，此类草地⽔分缺乏，草被稀疏，牧业利⽤条件差。
变化信息提取
波段替换法
变化信息提取方法
通常要求采用两种以上变化信息发现方法，有效防止变 化信息的遗漏。
光谱变异
主分量变化
监测结果பைடு நூலகம்
变化区域提取
➢ 阈值法 ➢ 区域生长法 ➢ 手工描绘法
阈值法
遥感图像中，各类地物都对应特定的灰度域。在变 化信息特征增强的图像上，变化区域的灰度值与其 它区域的灰度值一般是明显不同的。因此可以根据 直方图和影像特征，交互确定变化存在区灰度域的 上下限阈值；然后利用阈值将变化发生的区域从图 像中提取出来。
➢ 多波段主成份变换法 ➢ 主成份差异法
差异主成分法
前一时相融合影像
后一时相融合影像
影像相差取绝对值
主分量变换
变化信息 （第一分量）
差 异 主 成 分
多波段主成分变换法
前一时相融合影像 （n波段）
后一时相融合影像 （n波段）
两时相波段合成 （2n波段）
主分量变换 后几个分量波段组合
变化信息
多 波 段 主 成 分 变 换 法
手工描绘法
通过人机交互解译，从变化特征增强的图像中手 工描绘出变化范围。
解译的基本要素包括色调/颜色、大小、形状、纹 理、结构、高度、阴影、组合构型和所处的地理位 置等。
在PHOTOSHOP软件下，当变化图斑颜色单一时， 可用魔棒直接选择；反之则需进行勾绘，勾绘必须 在大于1：1的状态下，沿着融合图像的光谱与纹理 变异分界处以两个象素宽度向外勾绘图斑边界线。
差值法
分类后比较法
分类后比较法
分类后比较法是一种较为简单的变化发现 方法。首先运用统一的分类体系对每一时相 遥感影像进行分类,然后通过分类结果的比 较直接发现土地覆盖的变化。
基于遥感技术的土地利用／土地覆盖提取信息精度研究——以艾比湖湿地保护区为例
Ｅ ｖｒ ｍ ｎ ｌ ｒｔ ｔ ｎｏ ｉ ａｇ ｎｉ ｎ ｅｔ ｏ ｃ ｏ ｆ ￣ｉ ｏ ａＰ ｅｉ Ｘ ｎ
基 于 遥 感 技 术 的 土 地 利 用 ／ 地 覆 盖 提 取 信 息 精 度 研 究 土
— —
以艾 比湖 湿 地保 护 区为 例
４ 。７ 一 ５１ 总 面积 ２９６ ２ １９ ｋ ４３ ４ 。０ ， ５ ．７ ７ ｍ 。艾 比湖 湿
的评价就交付使用或入库 ， 缺乏较严格认真的精度评 估 。由于数据具有相应的复杂性噪声和误差 ， 数据 ＝（ 信息 ） 噪声 ＋（ 或误差） 而数据质量 又是数据的核心 ， ， 如果数据质 量 得不 到保 证 ， 数据 将 变得 毫 无 意义 … 。所 以 ， 据 数 的检验便是不可缺少 的一步 。数据是科学监测评价
遥感土地利用调查可为生态环境动态监测系统 提供本底数据 ， 而对所得数据 的验证是一个十分重要 却常被忽视的问题。很多情况下 ， 由遥感解译得到的
数据 ， 往 只经 过简单 的野外 实地 验证 进行 一些 定性 往
１ 研 究 区概 况 及数据 收集 １ １ 研 究 区概况 ．
区域范围是 在博尔塔拉蒙古 自治州 （ 以下简称 博州 ） 行政区域 内， 地理座标为 Ｅ８ 。Ｏ ～ ３５ Ｎ ２３ ８ 。０ ，
这 几种方法对该 区 Ｌ Ｃ Ｕ Ｃ的现 状 数 据 进 行 检 验 ， 果相 吻 合 。 结
关键 词 ： 感技 术 ；Ｕ Ｃ； 遥 Ｌ Ｃ 艾比湖 ； 精度 中圈分类号 ： ８ Ｘ７ 文献 标识码 ： Ａ 文章 编号 ：０ ８— ３ １ ２ １ ） ３— ０ ７— ４ １０ ２ ０ （ ０ １ ０ ０ ３ ０
土地利用变化遥感信息自动提取与自动整理
变化信息 和伪 变化信息等信息碎斑 ， 采用 图像收 缩与膨胀 运算
技 术能有效地进行 变化信息 的 自动整 理 ， 获得既符 合土地 利用 变化实际分布情况 又满足成 图精度要 求 的土地利用变 化 图像 ， 清晰地表现 出土地利 用变化的空 间分布和 强度 分布。 ［ 参考 文献 ］ ［ ］冯德俊 ， １ 李永树 ， 兰燕 ． 基于 主成 分变换 的动态监 测变化 信 息 自动发现 ［ ］ 计 算 机工 程与应 用 ，０ ４ ４ ３ ） １９— Ｊ． ２０ ，０（６ ：９
料 属于科技档案 管理 的范畴。在界定水土保持监测档 案收集范 围的基础上 ， 究 了监测档案案卷的主要 内容 与组成 、 研 档
案 工作 过 程 、 案 质 量 控 制 措 施 以及 档 案 的利 用 开 发 ， 出 了适 合 于水 土 保 持 监 测 与 信 息 系统 全 面运 行 后 日常 产 出 的 数 档 提
为 了实现 国家宏 观决 策 、 水土 流失 预测 预 报 、 水土 保持 科
研 以及水土保持工 作的信息化 ， 国家决定 分两期 建设 全 国水 土
保 持监测 网络 和信 息系统 ， 一期工 程主要 完成 了水利 部水 土保
持 监测中心 、 黄河流域 水 土保 持监测 中心站 、 江 流域 水土 保 长
持监 测 中心站和西部 １ １省及湖北 、 湖南等 省级监 测总 站 ， 以及
与监 测 总 站 对 应 的 １０个 监 测 分 站 的 建 设 ； 期 工 程 重 点 建 设 ０ 二
松辽 流域 、 淮河 流域 、 江流域 、 湖流域 ４个流 域机构 监测 中 珠 太
心站 ，８个省级 监测站及 ７ １ ５个监测 分站 ， 在全 国建设 ７ ８个 并 ３
地表覆盖知识点总结
第一章：绪论（★★★★）1.1地表覆盖与土地利用的基本概念1.1.1地表覆盖：是指自然营造物和人工建筑物所覆盖的地表诸多要素的综合体，包括地表植被、土壤、冰川、河流、湖泊、沼泽湿地及各种建筑物，主要侧重描述地球表面的自然属性，具有特定的时间和空间特性。
1.1.2土地利用：土地利用是土地利用方式、利用程度和利用效果的总称。
它包括的主要内容是：1、确定土地的用途；2、在国民经济各部门和各行业间合理分配土地资源；3、采取各种措施开发、整治、经营、保护土地资源，提高土地利用效果。
1.1.3地表覆盖与土地利用性质：1、地表覆盖反映了土地的自然属性（其性质主要取决于自然因素）2、影响土地利用的因素包括地表自然因素（气候、地形地貌、土壤、水文、地质条件）和社会经济因素（社会制度、政策、城市化与工业化、区位、交通条件、成本、效益、土地利用现状）（其性质主要取决于自然和经济因素）3、地表覆盖与土地利用有着密切的联系，其性质的变化构成了地表覆盖/土地利用变化4、自然力也会引起地表覆盖的变化5、地表覆盖具有特定的时间和空间属性，其形态和特征可在多种时空尺度上变化，而且产生变化的原因也具复杂多样性，并因此引起一系列的生态环境效应。
6、地表覆盖的空间分布反映着人类社会经济活动过程，决定着地表的水热和物质平衡，其变化直接影响到生物地球化学循环，改变着陆地-大气的水分、能量和碳循环，以至引起气候变化。
1.2地表覆盖与土地利用的研究内容1.2.1地表覆盖研究内容1、地表覆盖的特点；（与地理要素的关系）2、地表覆盖变化研究；3、地表覆盖分类体系研究；（满足地球系统模式需求；环境系统模式需求等）4、遥感测定地表覆盖及其变化；（数据尺度；时空分辨率；数据质量；数据处理；自动分类等）5、变化检测技术；6、地表过程模型化与动态模拟；举例：1、植被：利用遥感手段，可以在大范围内经济而有成效的查清植被资源和检测环境动态，从空间以不同尺度来研究地球植被层的空间结构和波动规律以及多种自然灾害和人类活动对生物圈的影响，并把植被遥感信息转换成图像和数据，供决策和管理者参考。
全国土地利用数据遥感信息提取土地利用数据说明
全国土地利用数据遥感信息提取土地利用数据说明北京揽宇方圆信息技术有限公司是中国科学院系统的遥感影像数据服务企业，专注于遥感影像数据一站式的基础卫星数据服务、卫星影像数据处理服务。
土地利用数据时间：1985年、1990年、1995年、2000年、2005年、2010年、2015年土地利用数据源：Landsat TM影像Landsat ETM影像土地利用数据遥感信息的提取：根据影像光谱特征，利用ARCGIS、易康软件、ENVI软件等，同时参照有关地理图件，对地物的几何形状，颜色特征、纹理特征和空间分布情况进行分析，提取土地利用信息。
土地利用/覆被变化信息的提取。
采用arcgis与易康结合，它通过分析地物光谱特征和其他图像特征，充分利用高程、坡度等地理辅助信息可以有效地提高分类精度，比较适合于地形破碎、地物分布复杂的地区。
基于Landsat TM遥感影像，采用全数字化人机交互遥感快速提取方法，同时参考国内外现有的土地利用/土地覆盖分类体系，以及遥感信息源的实际情况，将遥感影像进行解译并进行验证将土地利用数据类型划分为6个一级分类，24个二级分类以及部分三级分类的土地利用/土地覆盖数据产品，并结合本项目制定土地利用数据产品分类体系。
目视解译侧重于人的知识的参与，为了减少由于不同人员的主观差异性所造成的误差，提高遥感判读精度，因此建立统一解译标志是十分必要的。
根据影像光谱特征，结合野外实测资料，同时参照有关地理图件，对地物的几何形状，颜色特征、纹理特征和空间分布情况进行分析。
一、TM影像数据的预处理。
遥感数据处理主要包括大气校正、几何校正和图像增强，并利用行政边界矢量图对影像进行裁剪。
二、土地利用变化信息提取。
首先对其中的一期影像分别采用人工解译的方法，然后利用易康开始分类。
三、数据集成对数据形式特征（如格式、单位、分辨率、精度）等和内部特征（特征、属性、内容等）做出全部或部分的调整、转化、合成、分解等操作，形成充分兼容的数据库。
基于GIS的遥感影像土地利用／土地覆盖信息提取研究——以滇西北香格里拉县为例
似然 法 进 行 图 像 分 类 ，就会 导 致 分 类 精 度 降低 【 ” 。
虽然近 年来 ，人们 提 出了不少 研究遥 感 分类 的新 方
征 ，提 出多 步骤 遥感 影像 分类 方 法 （ １：基于 分 图 ） 辨 率融 合提 取地 物和 基 于 Ｎ Ｖ 计算 的 分 区 、分 层 Ｄ Ｉ
的思 想是通 过 分层和 分 区的途 径将 每一 次参 与监 督 分 类 的地 物类 别 尽 可 能 减 少 ， 以减 少 地 物 类 别 间 “ 物 同谱 、同谱 异 物 ” 的 干扰 ，在 此基 础 上再 利 异 用 Ｄ Ｍ 和坡 度数据 对 其进行 基 于知识 的影 像分 类 ， Ｅ 最后确 定 出像 元点 所属 的类 别 。
提取 精 度 达 ７ ％以 上 ，较 单 一的 信 息提 取 方 法 精 度有 很 大 提 高 。 ０
关 键 词 ：Ｇ Ｓ Ｉ ；土 地 利 用／ 地覆 盖 ；信 息提 取 ；香 格 里拉 县 土 中 图分 类 号 ：Ｘ ７ ８ 文献 标 识码 ：Ａ 文章 编 号 ：１０ ＿ ８ ２ （ｏ７ ２ ０ ９ — ５ ｏ １ ７ ５ ２ ０ ）０－ ０ ８ ０
图 １ Ｇ Ｓ支 持 下 的多 步 骤 遥 感 图 像分 类 流 程 图 Ｉ
Ｆ ｇ １ Ｆ ｏ ｉｇ ａ ｏ ｍｏ ｅ ｓ ｎ ｉ ｇｉ ｇ ｌｓ ｉ ｃ ｔ ｎ ｗ ｔ ｌ — ｔｇ ｙ ａ Ｇ Ｓ ｉ． ｌ ｗ ｄ ａ ｒ ｍ ｆ ＇ ｔ ｅ ｓｎ ｌ ｅ ｍａ ｅ ｃａ ｓ ａ ｏ ｉ ｍｕ ｔ ｓａ ｅ ｂ Ｉ ｉ ｆ ｉ ｈ ｉ
０ 引言
利 用遥 感 分 类 技术 进 行 土 地利 用／ 地覆 盖 分 土 类研 究一直 是遥感 应 用研究 的热点 。而 提 高计算 机
基于MODIS数据的朝鲜半岛土地利用变化信息提取
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
土地利用/覆盖变化信息提取实验报告1. 实验目的利用TM/ETM3个时相卫星数据，应用ENVI软件进行土地利用/覆盖分类，在此基础上进一步分析其动态变化特征。
2. 实验内容金华市土地利用/覆被变化信息的提取。
采用决策树分类法提取土地利用/覆被信息，它通过分析地物光谱特征和其他图像特征，充分利用高程、坡度等地理辅助信息可以有效地提高分类精度，比较适合于江南丘陵地形破碎、地物分布复杂的地区。
和传统的监督分类法相比，它可以消除园地和林地、建设用地和裸地光谱相似所带来的影响。
（1）TM影像数据的预处理。
本文的遥感数据处理主要包括大气校正、几何校正和图像增强，并利用行政边界矢量图对影像进行裁剪。
（2）土地利用变化信息提取。
首先对其中的一期影像(2003年)分别采用最大似然法、决策分类树法进行分类，提取土地利用/覆被信息，并对二者的提取精度进行比较，选择精度最高者作为最终的提取方法，进而提取1988～2003年金华市土地利用/土地覆被信息。
（3）利用空间叠加获取土地利用/覆被变化的面积转移矩阵，进而通过面积转移矩阵分析土地利用/土地覆被的数量变化、空间结构变化和土地利用程度。
3. 实验方案4. 数据预处理4.1 数据源本文所采用的数据包括：两景金华市的Landsat TM和一景Landsat ETM陆地卫星影像，一景半SPOT 全色影像；该地区1：50 000地形图；该地区81m*81m分辨率的数字高程模型（DEM）；1：100万中国行政边界矢量图等。
具体的见表4-1和4-2所示。
表4-1 研究区遥感影像数据获取时间传感器类型数量（景）空间分辨率（m）2003年3月9日SPOT-5全色15 1/252003年3月26日LandsatETM+1-8波段115m（全色）30m(多光谱)1996年9月6日LandsatTM1-7波段1301988年12月5日LandsatTM1-7波段130表4-2 研究区其他资料及应用说明数据类型应用说明大比例尺地形图最新时相的1：50000地形图，用于进行卫星遥感资料的几何校正野外调查资料野外控制点的测量，土地利用/覆盖分类训练样本区的调查，建立判读标志，进行分类及信息提取精度检验等工作土地利用现状图对比土地利用/覆盖动态变化及遥感影像分类精度参考4.2 图像预处理数据预处理部分主要包括：对遥感影像进行大气校正、几何纠正、以及对研究区进行边界裁剪和图像增强。
主要工作流程如下（图4-2）：图4-2数据预处理主要技术流程具体方法如下：1）大气校正。
本文的大气纠正在PCI软件的ATCOR2模块中完成，以去除薄云和大气对影像光谱的影响，尤其是96年的图像，经过大气校正后，图像质量得到了很大提高。
2）几何纠正。
影像几何纠正就是将所研究影像纳入到一个地面坐标系中，方法是利用地面控制点对各种因素引起的遥感影像的几何畸变进行纠正，以便确定影像上每个像元在地面的坐标，其过程就是把目标由一个空间向另一个空间转换的过程。
3）边界裁剪。
对遥感影像进行上述处理后，利用所给的金华县行政边界矢量图边界对影像进行裁剪，分别得到1988年和1996年的TM图像以及2003年的SPOT和ETM研究区影像。
4）图像增强。
本文采用最佳指数因子分析方法对3幅多光谱影像三波段组合方案进行评价最后得出最佳波段组合。
最佳指数因子的计算公式为：（2-1）式中，为i波段图像的亮度标准差，其值越大，说明数据的离散度越大，所包含的信息量越大，可分离性越高；为三波段中任意两波段之间的相关系数，其值越小，表明图像数据独立性越高，信息的冗余度越小。
OIF越大，组合图像的信息量越大，组合方案越佳。
利用最佳指数因子分析方法计算1988、1996、2003年三幅多光谱影像的三波段组合值如下表（表4-3）：表4-3影像最佳波段组合信息含量表影像OIF组合方案42143253254375274375488年TM26.33829.18225.52133.60527.79333.98333.91096年TM32.15434.86416.84536.02516.08936.16325.59803年ETM16.24019.48321.12123.24020.88022.59220.427从上表中可得知，1988年和1996年两幅TM影像的743组合值都是最大的，2003年的ETM543组合值最大，743次之，整个金华市的土地利用格局以耕地和林地等农用地为主，整个地区有较高的植被覆盖，而743组合更接近植被的真彩色，有利于植被的分类，所以三幅影像都采用743组合方案参与分类。
5 土地利用/土地覆被分类5.1 监督分类法此法的关键在于训练区的选择。
训练区的选取应与分类地区的特点和分类系统相适应。
对训练区的统计特征应进行详细的分析，以选择最有效的参数变量(谱段)参与后续的分类。
此外，应对训练区特征指标的外延性进行评估（赵英时，2003）。
监督分类法中具体方法包括最小距离分类法、多级切割分类法、特征曲线窗口法和最大似然分类法，其中最大似然分类法用的最多。
最大似然分类法（maximum likelihood classifier）在多类别分类时，常常采用统计学方法建立起一个判别函数集，然后根据这个判别函数集计算各待分像元的归属概率。
这里，归属概率是指：对于待分像元x，它从属于分类类别的（后验）概率。
设从类别中观测到的条件概率为，则归属概率可表示为如下形式的判别函数：（3-1）式中，为待分像元，为类别的先验概率，它可以通过训练区来决定。
此外，由于上式中分母和类别无关，在类别间比较的时候可以忽略。
1)训练样区的选取与纯化本文中采用的方法为试分类混淆矩阵分析法和J-M距离法来对所选取的训练样本纯度进行评价，通过对所得混淆矩阵进行分析，即可得到训练样本占原来各类个体总数的百分比，以确定其分类的正确率，从而也检验了训练的纯度。
在ENVI 4.2中对纯化前和纯化后的训练样本区进行了分类，并以训练样本自身对分类结果进行精度检验，得到混淆矩阵。
纯化前训练样本的训练区分类混淆矩阵显示：Overall Accuracy = 92.0142% ，Kappa Coefficient = 0.9165。
对训练样本进行纯化后，训练样本的训练区分类混淆矩阵显示：Overall Accuracy = 96.3045% ，Kappa Coefficient = 0.9500。
在ENVI 4.2中未纯化前和纯化后的训练样本区进行统计其J-M距离，结果显示：纯化前训练样本的J-M距离，最小的是园地和林地之间的J-M距离，只有1.3208。
另外园地和耕地之间、居民点及工矿用地和未利用地之间的J-M距离也比较小，但均大于1.5。
经纯化以后，只有园地和林地之间的J-M距离仍然小于1.5。
其他均在1.8以上。
训练样本之间的可分性明显增大。
2)分类精度及结果分析在对影像做了图像预处理后，根据影像的光谱特征，选择适当的训练样区，将土地利用分为若干不同的类别。
结合研究区的实际情况，根据目视判读选择训练样本，采用最大似然分类法(Maximum Likelihood Classification, MLC)，将金华市2003年的土地利用分为5类，即耕地、园地、林地、居民点及工矿用地、水域。
（a）1988年（b）1996年（c）2003年图5-1 金华市最大似然法土地利用/土地覆被分类图利用上面确定的方法和已有的样本数据对研究区分类后图像进行精度估计，得到最终最大似然法的误差矩阵和各项统计指标（以2003年为例）。
表5-1 最大似然法分类精度矩阵类别未分类耕地园地林地居民点及工矿用地水体未利用地总计未分类00000000耕地0658000073园地0314*******林地1119400097居民点及工0002140016矿用地水体0000010010未利0004801224用地合计16923119221012256表 5-2 最大似然法分类精度评价（％）类别参考总计分类总计正确分类数生产精度用户精度未分类100------耕地69736594.20%89.04%园地23361460.87%38.89%林地119979478.99%96.91%居民点及工矿22161463.64%87.50%用地水体101010100.00%100.00%未利122412100.00%50.00%用地合计256256209总精度＝209/256＝81.64％表5-1和5-2显示了在最大似然分类中，园地和未利用地的分类精度比较低，分别为：38.89%和50.00%，这是因为园地、耕地、林地之间和未利用地和居民点及工矿用地之间的光谱特征很多比较相似（图5-2），从而导致误判。
尤其是林地和园地错分现象比较多，36个园地样点中，有19个被误分为林地，3个被误分为耕地。
5.2 决策树分类1）典型地物光谱分析为了获取研究区内各种地物类型光谱特征知识，对研究区内典型地物类型进行采样并加以统计，统计结果见表5-3（以03年为例）。
其中03年耕地按耕地l(主要为水田)和耕地2(主要为旱地)两类采样，水体分为水体1(河流)和水体2(水库和坑塘水面)两类采样。
88年图像上由于存在大量山体阴影，因此阴影也列为单独一类进行采样。
表5-3 03年典型地物样本亮度值统计表地物波段耕地耕地2园地林地居民点及工矿用地水体1水体2未利用地B1最小值 4.0014.008.000.009.000.00 5.00 1.00最大值40.0041.0029.0024.0090.0047.0035.00101.00均值21.0624.2215.767.5628.7531.2514.2536.07均方差 4.03 3.58 2.88 2.85 5.93 4.29 4.499.45B2最小值10.0021.0013.000.0016.000.009.008.00最大值59.0062.0047.0032.00101.0079.0055.00118.00均值35.8735.5724.2115.2938.3548.3819.8757.00均方差 5.54 4.68 4.06 3.527.388.287.7112.52B3最小值10.0018.0011.000.0018.000.00 6.0013.00最大值76.0067.0057.0032.00117.0081.0063.00138.00均值35.2838.3623.3813.5245.3146.6415.3772.54均方差 6.67 6.13 5.79 3.569.199.437.9918.70B4最小值25.0019.0022.000.0017.000.000.0017.00最大值109.0090.0072.0067.0082.0065.0063.0095.00均值71.2042.1246.4341.8836.7813.18 5.4064.83均方差10.288.47 6.089.527.708.43 5.7510.54B5最小值30.0011.0020.000.0024.000.00 4.0026.00最大值117.00110.00102.0092.00172.0099.0095.00167.00均值75.4347.5860.5047.7567.1814.4310.41106.08均方差9.2617.839.7711.6212.1111.347.2118.14B7最小值16.007.0012.000.0019.000.00 1.0023.00最大值97.0089.0080.0059.00163.0075.0068.00155.00均值47.2431.0135.7424.5160.7110.597.5286.17通过03年典型地物样本亮度值统计可以得到03年典型地物样本波谱响应曲线：图5-2 03年典型地物样本波谱响应曲线2）决策树的构建对于03年图像：由于水体和阴影的低反射率，尤其是在长波部分更明显。