土地利用覆盖变化地信息提取

1、实用标准土地利用/覆盖变化信息提取实验报告1.实验目的利用tm/etm3个时相卫星数据，应用envi软件进行土地利用/覆盖分类，在此基础上进一步分析其动态变化特征。2.实验内容金华市土地利用/覆被变化信息的提取。采用决策树分类法提取土地利用/覆被信息，它通过分析地物光谱特征和其他图像特征，充分利用高程、坡度等地理辅助信息可以有效地提高分类精度，比较适合于江南丘陵地形破碎、地物分布复杂的地区。和传统的监督分类法相比，它可以消除园地和林地、建设用地和裸地光谱相似所带来的影响。（1）tm影像数据的预处理。本文的遥感数据处理主要包括大气校正、几何校正和图像增强，并利用行政边界矢量图对影像进行裁剪。（2

2、）土地利用变化信息提取。首先对其中的一期影像(2003年)分别采用最大似然法、决策分类树法进行分类，提取土地利用/覆被信息，并对二者的提取精度进行比较，选择精度最高者作为最终的提取方法，进而提取19882003年金华市土地利用/土地覆被信息。（3）利用空间叠加获取土地利用/覆被变化的面积转移矩阵，进而通过面积转移矩阵分析土地利用/土地覆被的数量变化、空间结构变化和土地利用程度。3.实验方案4.数据预处理4.1数据源文案大全实用标准本文所采用的数据包括：两景金华市的landsattm和一景landsatetm陆地卫星影像，一景半spot全色影像；该地区1：50000地形图；该地区81m*81m分

3、辨率的数字高程模型（dem）；1：100万中国行政边界矢量图等。具体的见表4-1和4-2所示。表4-1研究区遥感影像数据获取时间2003年3月9日2003年3月26日1996年9月6日1988年12月5日传感器类型spot-5全色landsatetm+1-8波段landsattm1-7波段landsattm1-7波段数量（景）11/2111空间分辨率（m）5515m（全色）30m(多光谱)3030表4-2研究区其他资料及应用说明数据类型大比例尺地形图野外调查资料土地利用现状图应用说明最新时相的1：50000地形图，用于进行卫星遥感资料的几何校正野外控制点的测量，土地利用/覆盖分类训练样本区的调

4、查，建立判读标志，进行分类及信息提取精度检验等工作对比土地利用/覆盖动态变化及遥感影像分类精度参考4.2图像预处理数据预处理部分主要包括：对遥感影像进行大气校正、几何纠正、以及对研究区进行边界裁剪和图像增强。主要工作流程如下（图4-2）：文案大全实用标准图4-2数据预处理主要技术流程具体方法如下：1）大气校正。本文的大气纠正在pci软件的atcor2模块中完成，以去除薄云和大气对影像光谱的影响，尤其是96年的图像，经过大气校正后，图像质量得到了很大提高。2）几何纠正。影像几何纠正就是将所研究影像纳入到一个地面坐标系中，方法是利用地面控制点对各种因素引起的遥感影像的几何畸变进行纠正，以便确定影像

5、上每个像元在地面的坐标，其过程就是把目标由一个空间向另一个空间转换的过程。3）边界裁剪。对遥感影像进行上述处理后，利用所给的金华县行政边界矢量图边界对影像进行裁剪，分别得到1988年和1996年的tm图像以及2003年的spot和etm研究区影像。4）图像增强。本文采用最佳指数因子分析方法对3幅多光谱影像三波段组合方案进行评价最后得出最佳波段组合。最佳指数因子的计算公式为：（2-1）文案大全式中，分离性越高；实用标准为i波段图像的亮度标准差，其值越大，说明数据的离散度越大，所包含的信息量越大，可为三波段中任意两波段之间的相关系数，其值越小，表明图像数据独立性越高，信息的冗余度越小。oif越大，

6、组合图像的信息量越大，组合方案越佳。利用最佳指数因子分析方法计算1988、1996、2003年三幅多光谱影像的三波段组合值如下表（表4-3）：表4-3影像最佳波段组合信息含量表影像421432532oif组合方案54375274375488年tm96年tm03年etm26.33832.15416.24029.18234.86419.48325.52116.84521.12133.60536.02523.24027.79316.08920.88033.98336.16322.59233.91025.59820.427从上表中可得知，1988年和1996年两幅tm影像的743组合值都是最大的，20

7、03年的etm543组合值最大，743次之，整个金华市的土地利用格局以耕地和林地等农用地为主，整个地区有较高的植被覆盖，而743组合更接近植被的真彩色，有利于植被的分类，所以三幅影像都采用743组合方案参与分类。5土地利用/土地覆被分类5.1监督分类法此法的关键在于训练区的选择。训练区的选取应与分类地区的特点和分类系统相适应。对训练区的统计特征应进行详细的分析，以选择最有效的参数变量(谱段)参与后续的分类。此外，应对训练区特征指标的外延性进行评估（赵英时，2003）。监督分类法中具体方法包括最小距离分类法、多级切割分类法、特征曲线窗口法和最大似然分类法，其中最大似然分类法用的最多。最大似然分类

8、法（maximumlikelihoodclassifier）在多类别分类时，常常采用统计学方法建立起一个判别函数集，然后根据这个判别函数集计算各待分像元的归属概率。这里，归属概率是指：对于待分像元x，它从属于分类类别的（后验）概率。设从类别中观测到的条件概率为，则归属概率可表示为如下形式的判别函数：（3-1）式中，为待分像元，为类别的先验概率，它可以通过训练区来决定。此外，由于上式中分母和类别无关，在类别间比较的时候可以忽略。1)训练样区的选取与纯化文案大全实用标准本文中采用的方法为试分类混淆矩阵分析法和j-m距离法来对所选取的训练样本纯度进行评价，通过对所得混淆矩阵进行分析，即可得到训练样本

9、占原来各类个体总数的百分比，以确定其分类的正确率，从而也检验了训练的纯度。在envi4.2中对纯化前和纯化后的训练样本区进行了分类，并以训练样本自身对分类结果进行精度检验，得到混淆矩阵。纯化前训练样本的训练区分类混淆矩阵显示：overallaccuracy=92.0142%，kappacoefficient=0.9165。对训练样本进行纯化后，训练样本的训练区分类混淆矩阵显示：overallaccuracy=96.3045%，kappacoefficient=0.9500。在envi4.2中未纯化前和纯化后的训练样本区进行统计其j-m距离，结果显示：纯化前训练样本的j-m距离，最小的是园地和林

10、地之间的j-m距离，只有1.3208。另外园地和耕地之间、居民点及工矿用地和未利用地之间的j-m距离也比较小，但均大于1.5。经纯化以后，只有园地和林地之间的j-m距离仍然小于1.5。其他均在1.8以上。训练样本之间的可分性明显增大。2)分类精度及结果分析在对影像做了图像预处理后，根据影像的光谱特征，选择适当的训练样区，将土地利用分为若干不同的类别。结合研究区的实际情况，根据目视判读选择训练样本，采用最大似然分类法(maximumlikelihoodclassification,mlc)，将金华市2003年的土地利用分为5类，即耕地、园地、林地、居民点及工矿用地、水域。（a）1988年（b）1

11、996年（c）2003年图5-1金华市最大似然法土地利用/土地覆被分类图利用上面确定的方法和已有的样本数据对研究区分类后图像进行精度估计，得到最终最大似然法的误差矩阵和各项统计指标（以2003年为例）。表5-1最大似然法分类精度矩阵未分居民点及工类别耕地园地林地水体未利用地总计类矿用地未分类00000000文案大全实用标准耕地园地林地0016531814101994000000000733697居民点及工0002140016矿用地水体0000010010未利0004801224用地合计16923119221012256表5-2最大似然法分类精度评价（）类别未分类耕地园地林地居民点及工矿参考总计

12、1692311922分类总计073369716正确分类数065149414生产精度-94.20%60.87%78.99%63.64%用户精度-89.04%38.89%96.91%87.50%用地水体未利101210241012100.00%100.00%100.00%50.00%用地合计256256209总精度209/25681.64表5-1和5-2显示了在最大似然分类中，园地和未利用地的分类精度比较低，分别为：38.89%和50.00%，这是因为园地、耕地、林地之间和未利用地和居民点及工矿用地之间的光谱特征很多比较相似（图5-2），从而导致误判。尤其是林地和园地错分现象比较多，36个园地样点

13、中，有19个被误分为林地，3个被误分为耕地。5.2决策树分类1）典型地物光谱分析为了获取研究区内各种地物类型光谱特征知识，对研究区内典型地物类型进行采样并加以统计，统计结果见表5-3（以03年为例）。其中03年耕地按耕地l(主要为水田)和耕地2(主要为旱地)两类采样，水体分为水体1(河流)和水体2(水库和坑塘水面)两类采样。88年图像上由于存在大量山体阴影，因此阴影也列为单独一类进行采样。表5-303年典型地物样本亮度值统计表文案大全实用标准b1b2b3b4b5b7地物波段最小值最大值均值均方差最小值最大值均值均方差最小值最大值均值均方差最小值最大值均值均方差最小值最大值均值均方差最小值最大值

14、均值均方差耕地4.0040.0021.064.0310.0059.0035.875.5410.0076.0035.286.6725.00109.0071.2010.2830.00117.0075.439.2616.0097.0047.249.46耕地214.0041.0024.223.5821.0062.0035.574.6818.0067.0038.366.1319.0090.0042.128.4711.00110.0047.5817.837.0089.0031.0113.23园地8.0029.0015.762.8813.0047.0024.214.0611.0057.0023.385.79

15、22.0072.0046.436.0820.00102.0060.509.7712.0080.0035.748.59林地0.0024.007.562.850.0032.0015.293.520.0032.0013.523.560.0067.0041.889.520.0092.0047.7511.620.0059.0024.517.13居民点及工矿用地9.0090.0028.755.9316.00101.0038.357.3818.00117.0045.319.1917.0082.0036.787.7024.00172.0067.1812.1119.00163.0060.7111.64水体10.

16、0047.0031.254.290.0079.0048.388.280.0081.0046.649.430.0065.0013.188.430.0099.0014.4311.340.0075.0010.598.38水体25.0035.0014.254.499.0055.0019.877.716.0063.0015.377.990.0063.005.405.754.0095.0010.417.211.0068.007.525.26未利用地1.00101.0036.079.458.00118.0057.0012.5213.00138.0072.5418.7017.0095.0064.8310.54

17、26.00167.00106.0818.1423.00155.0086.1716.75通过03年典型地物样本亮度值统计可以得到03年典型地物样本波谱响应曲线：图5-203年典型地物样本波谱响应曲线文案大全实用标准2）决策树的构建对于03年图像：由于水体和阴影的低反射率，尤其是在长波部分更明显。因此可以利用tm4/tm3的比值来区分大部分林地。然后通过缨帽变换的亮度分量可以区分大部分耕地，剩下的耕地、园地和林地可以通过地学辅助知识如高程和坡度信息来加以区分。最后通过分析得到03年的分类决策树：图4-503年分类决策树依据同样的原理，通过分析分别得到88年和96年的分类决策树：图5-696年分类决

18、策树文案大全实用标准图5-788年分类决策树3）分类精度及结果分析通过决策树分类以及分类后的合并处理最后得到三个时期的土地利用/覆被分类图：（a）1988年（b）1996年（c）2003年图5-8金华市决策树分类法土地利用/土地覆盖分类图通过分层采样和已有的样本数据对研究区分类后图像进行精度估计，得到最终决策树的误差矩阵和各项统计指标（以2003年为例）。表5-6决策树分类精度矩阵未分居民点及工类别耕地园地林地水体未利用地总计类矿用地未分类00000000文案大全实用标准耕地园地林地00077163112301050010000008312114居民点及工0200211024矿用地水体1000

19、012013未利000020810用地合计1861610823139256表5-7决策树分类精度评定（）类别未分类耕地园地林地居民点及工矿参考总计1861610823分类总计0831211424正确分类数0771110521生产精度-89.53%68.75%97.22%91.30%用户精度-92.77%91.67%92.11%87.50%用地水体未利139131012892.31%88.89%92.31%80.00%用地合计256256234总精度209/25691.41表5-6和5-7显示了决策树分类中，园地和未利用地的分类精度明显有所提高，分别为：91.67%和80.00%，林地和园地错分

20、现象明显减少，12个园地样点中，只有1个被误分为耕地，同样，在10个未利用地样点中，只有2个样点被误分为居民点及工矿用地。这是因为在决策树分类过程中，充分利用了地物的各种特征信息，除了利用园地、耕地、林地之间和未利用地和居民点及工矿用地之间的光谱特征之外，还有效使用了其他的特征信息，如纹理信息，通过ndvi可以很好的初步分离耕地、园地和林地，利用缨帽变换的亮度、湿度和绿度信息可以很好的区分未利用地和居民点及工矿用地，在初步分类结果的基础上，通过gis辅助数据，像高程和坡度等信息，可以进一步分离混分部分，使分类精度大大提高。5.3两种分类精度评定及结果分析对两种分类方法的分类结果分别进行精度评定

21、最后得到了两种分类方法三个时期总体分类精度：表5-8两种分类方法的分类精度比较年份最大似然法决策树分类法文案大全实用标准方法1988年1996年2003年总体精度84.776.9581.64kappa系数0.78460.69490.7461总体精度87.8987.8991.41kappa系数0.82880.82900.8751从上述表格两种方法分类结果的比较可以看出：（1）决策树分类方法从整体上优越于最大似然法，而且整体上比较平均，该方法在各单项地物类型的分类精度之间相差比较小，分类结果较为满意，总精度达到了85以上，均达到最低允许判别精度0.7的要求（lucasifj，1994）。（2）基于

22、知识的决策树分类方法是提高遥感影像土地利用/覆被自动分类精度的有效途径之一，其关键在于知识的获取与规则的创建，只有在对地类光谱特征、几何特征、地类间相互关系、地类变化规律、地学分布规律等知识认真分析的基础上建立合理的规则，通过有效地推理判断才能得到高质量的分类结果。本研究在对金华研究区地物分布规律、地物光谱特征分析的基础上所建立的规则用于分类取得了较好的效果，说明其规则的创建是合理而有效的。6土地利用/覆被变化分析6.1金华市土地利用类型的数量变化根据土地利用/覆被遥感图像分类的结果，1988年、1996年和2003年金华市各类土地利用/覆被类型的面积和比例的变化，如表4-1和4-2所示。表6

23、-11988年-1996年金华市土地利用总量变化分析表土地利用类1988年1996年面积增加比例比例变化率型面积/hm2面积/hm2/hm2耕地园地林地居民点及工103476.242442.9687942.335617.9850.191.1942.652.7267770.7215425.64102268.3512702.5132.877.4849.66.16-35705.5212982.6814326.027084.53-34.51531.4316.29126.10矿用地水体未利用地合计3111.033586.05206176.591.511.741005589.92419.47206176.

24、592.721.171002478.87-1166.580.0079.68-32.530.00文案大全实用标准表6-21996年-2003年金华市土地利用总量变化分析表土地利用类1996年2003年面积增加比例比例变化率型面积/hm2面积/hm2/hm2耕地园地林地居民点及工67770.7215425.64102268.3512702.5132.877.4849.66.1672926.738804.25101124.2715782.6735.374.2749.057.655156.01-6621.39-1144.083080.167.61-42.92-1.1224.25矿用地水体未利用地合计5

25、589.92419.47206176.592.721.171005588.461950.21206176.592.710.95100-1.44-469.260.00-0.03-19.400.006.2金华市土地利用类型转换变化单纯看面积增加或减少，难以反映土地利用的内部结构变化。为揭示各土地利用类型的内部转移特征，基于1988年、1996年和2003年金华市土地利用分类图，利用arcgis软件中空间分析模块的叠加功能，分别将每两个时期分类图像叠加，得到土地利用面积转移矩阵（表6-3、表6-4和表6-5）。表6-319881996金华市土地利用转换矩阵（unit：hm2）动态转移矩阵耕地园地19

26、96年居民点及林地工矿用地水域未利用地1988年合计耕地p/%园地p/%林地p/%61793.371423.7171.032287.6265.6013575.9632.57438.48745.2921.3716128.6338.70421.3821.0284457.88208.4519.6969.33.46274.777.882044.894.9119.440.97115.473.311720.444.1370.653.5264.081.84103476.24100.002442.96100.0087942.33100.001988年居民点及工矿用地651.24277.56651.243618

27、.36321.6697.925617.98p/%水域p/%未利用地p/%1996年合计32.57131.0424.731483.7447.1267770.7213.8868.1312.86318.4210.1115425.6432.56134.4625.38475.215.09102268.35161.3730.46369.6311.7412701.5216.092580.93501.6615.935589.94.9034.836.57437.42419.47100.003111.03100.003586.05100.00206176.59文案大全实用标准总的来说，耕地和未利用地一直在减少，耕

28、地减少的最多，是其他用地的主要来源；居民点及工矿用地的面积呈大幅度地增加，主要来自于耕地；林地、园地、水域面积有不同程度的增加。1）耕地转换分析耕地在前期（1988年1996年）向各种类型均有不同程度的转变，其中主要向园地、林地和居民点及工矿用地转变，转化量分别为13575.96hm2、16128.63hm2、8204.54hm2，分别占由耕地转换为其他用地总量的32.57、38.7.和19.69，在后一时期（1996年2003年）耕地有稍微的增加，主要由林地、园地、居民点及工矿用地转化而来，分别为8088.21hm2、8047.98hm2、3402.36hm2。总的来说，耕地面积一直在减少，从1988年到2003年，主要转化为园地、林地和居民点及工矿用地，分别达到了6732.99hm2、15919.29hm2、11062.26hm2，占用由耕地转为非耕地的18.14、42.88和29.80。2）园地转换分析园地的变化起伏比较大，先增加后减少。前期主要是大幅度地增加，转出主要是变为耕地和林地，增加的面积主要由耕地而来，达到了13575.96hm2，后期面积逐渐减少，到2003年时主要转向了耕地和园地，分别为8088.21hm2、327402hm2，占了由园地变为非园地总量的62.09和