遥感影像解译不确定性的评价与表达

《遥感影像解译不确定性的评价与表达》由会员分享，可在线阅读，更多相关《遥感影像解译不确定性的评价与表达（13页珍藏版）》请在金锄头文库上搜索。

1、遥感影像解译不确定性的评估与表达摘自遥感数据的不确定性问题承继成 郭华东 史文中等编著遥感数据的精度评估研究是从1975年开始的(1973年发射第一个遥感卫星 )。最早Hord和Brooner(1976)， Van Genderen 和 Lock(1977)及 Ginevan(1979) 曾提出了建 立测试评估地图的标准和技术的建议。Roslnfield(1982)，Congalton(1983) ，Aronoff(1985) 对遥感数据精度的评估标 准和技术进行了较深入的研究，以后又有更多的人参与了该项研究工作。误差矩阵是主要的方法，它能很好地表达专 题图的精度，已经成为普遍采用的方法。一、

2、遥感影像解译不确定性评估综述遥感解译有人工目视判读和计算机自动分类处理。在本 章中我们主要指计算机自动分类。造成遥感影像解译不确定 性的原因有遥感数据固有的不确定性 和遥感数据获取、处 理、传输、分类过程造成的误差。因此遥感数据解译过程中 的不确定性是客观存在、不可避免的。任何解译的成果图件 在不同程度上都存在着一定的不确定性，符合“任何人工模 拟产品与客观真实世界之间总是存在一定差异”的原理。遥感影像数据的不确定性是普遍存在的。一些遥感影 像的分辨率很低，经过各种处理影像分类的可信度尽管有所提高但仍然存在不确定性(表1)，一些地物的可信度仍很低。表1遥感影像分类的可信度()(据吴连喜，20X

3、X)地类城镇建筑农村居民点裸地大棚耕地园地林 地水体道路TM影像 Marr融合影像Brovey融合影像HIS融合影像7510PCA融合影像5839 5487遥感数据分类的不确定性度量方法通常用误差矩阵来度量。从误差矩阵中可以计算出分类精度的指标，如“正确分类比”。另一种指标是Cohen提出来的Kappa系数，后来经Foody(1992) 修正后称为Tau系数。遥感数据分类的专题不确定性是指专题值与其真值的接近程度，其度量随专题数据类型的不同而不同。专题数据的类型有两种：分类专题数据(categorical thematic data)和连续专题数据(continuous thematic da

4、ta), 也有将其分为定性数据 (qualitative data)和定量数据的(quantitative data) 。连续数据的不确定性度量指标与位置不确定性的度量指标相类似，如方差等 (Lanter and Veregin,1992 ； Heuvelink,1993 ； Goodchild et al,1992)。遥感数据不确定性的度量一般采用基于像元的分类结果评估，其不确定性度量评估流程如图图1基于像元的遥感数据不确定性评估流程图(据Lunetta et al,1991)二、基于采样的检验方法总结现有的文献主要有三种基于实验的检验方法：(1) 对于某一类或全集正确量测的百分比(Rose

5、nfield ， 1986);(2) 某一置信水平下某一类或全集正确量测的百分比(Aronoff,1985 ； Hord and Brooner,1976) ；(3) 基于某些参数的某一类或全集正确量测的百分比(Greenland et al,1985； Rosenfield andFitzpatrick-Lins,1986) 。以上三种方法适用于各种非连续属性值的精度评估。非连续属性数据的评估可以通过对一组分类结果的评价得以实现。地面真实数据有时也称之为参考数据，通过将量测数 据与参考数据的比较我们可以建立一个误差矩阵。该矩阵可 以描述某一类别的分类精度或整体分类精度。基于这一矩阵 可以对分

6、类精度进行进一步地讨论。本节将集中讨论遥感分类影像的精度评估问题。1. 参考数据的采样在基于采样数据的属性不确定性评估方法中，采样数据作为误差矩阵或其他统计分析的输入部分，在这个过程中选择适当的采样数据是非常重要的。有两个因素影响着采样数据的选择：采样样本的大小及采样数据分布模式。在评估分类遥感影像的分类精度时，采样点的数量是十 分重要的。获取地面采样数据是昂贵的，因此采样点的个数 应尽可能减少。另一方面，为了在统计上有意义，采样点的数 目应尽可能大，至少大于某一给定数目，例如30个。人们在 采样点个数方面进行了许多探讨(Van Genderen and Lock,1977 ； Rosenfi

7、eld ，1982; Congalton ，1988; Fukunaga and Hayes ， 1989)。另一个重要因素是采样模型。选择适当的采样点分布方 式,使所选择的样本可以代表全部分类的影像是十分重要 的。一个较差的采样可能导致精度评估的偏重，使得对精度的估计过高或过低。通常使用的采样模式有五种：简单随机 采样(Cochran ，1977);集群采样(Kish，1965);分类随机 采样(Barrett and Nutt,1979) ;系统采样 (Barreff and Nutt, 1979)以及分类系统非一致采样 (Berry and Baker,1968)。Congalton(1

8、988)对不同地区进行了采样模拟，并总结出在各种情况下简单随机采样与分类随机采样提供最佳采 样结果的规律。当采样个数与模式确定之后，即可实施采样，进而生成误差矩阵，依此进一步进行有关的属性精度评估。2. 误差矩阵 误差矩阵，有时也称为混淆矩阵，是一个用于表示分为某一类别的像素个数与地面检验为该类别数的比较阵列。一个误差矩阵的实例见表2。表2中的列通常表示参考数据，而行表示遥感分类的 结果。误差矩阵通常用于表示分类的精度，因为它可用于指 出某一类的或整体分类的精度。此外，用误差矩阵还可以表 示出包含与丢失两种误差。在表2中，A B、C是三类待分类的类别名称。第一行的数目，如总数“52”表示为A的

9、像素中根据实地检查有45个被分为类别“ A，2个像素被分为类别“ B”，5个像素被分为“ C”。表2误差矩阵实例分类 数据 A B C列总和 A 45 3 2 50地面真实数据 B 2 63 7 70 行总数 C 5 5 70 80 52 71 79 200根据误差矩阵可以导出若干关于总体分类或对于某一类别分类的精度描述指标。分类的总体精度是用误差距阵内对角线元素之和除以总的采样个数来表示的。例如，在表2的例子中，该值为(45+63+70)/200 = 89%即总体分类精度为89%为描述对某一类别的分类精度，我们定义了用户精度和生产者精度 (Story and Congalton , 1986

10、)。对于类别A的生产者精度是用下列公式计算的，即类别A的正确分类个数除以对于类别 A的总采样个数，即 A的列 总和。例如在以上的例子中，生产者精度为：45/50= 90%。 该指标指出了一个地面采样点被正确地分类的概率。事实 上，它是对丢失误差的一个量度，该误差指出了该采样数据 中没有被正确分类的百分比。丢失误差是该类所在列中非对 角线元素之和除以该列总和而得。例如在以上的例子中，丢 失误差为：(3+2)/50=10%。因此，有：生产者精度+丢失误差=100%另一方面，类别 A的用户精度定义为：正确分类为A的个数除以分为类别 A的总和(即A所在行的总和)。在以上 例子中，该值为45/52=%。

11、事实上，该指标指出一个采样分 类点表示实际地面真实情况的百分比。用户精度表示了包含 误差，A类的包含误差用 A所在行的非对角线元素之和除以 该行的总和。在以上的实例中，该值为：(2+5)/52= %。用户精度与包含误差有以下关系：用户精度+包含误差=100%Chrisman(1986)指出，输入GIS中数据应附有一个误差 矩阵。这应以原始的误差矩阵形式表示，而非该矩阵导出的 一系列参数。只有这个原始矩阵才能表示出每一类别的各种 精度与误差，用户可根据其自己的要求从中导出新的参数。系数两个最常用的属性精度量测量是二维正态概率和Kappa系数。于二维正态概率是基于“正确百分率”，因而 不能统计出包

12、含与丢失误差。另一方面， Kappa系数提供两 幅图观测协议的不同量度，而协议是几率形成的 (Congalton and Mead， 1983 ； Cangalton et al ， 1983)。Kappa 系数定 义为：K=(P0-Pe(1 -Pe)式中，P0是观测精度估计，而 Pe是期望精度估计。一个的Kappa系数可以解释为该分类以 80%的程度优于 随机地给像素赋类别值。Kappa系数的优点在于它已经包含 了丢失误差和包含误差。一个条件Kappa系数可以表示对于某一类别的分类精度。建议Kappa参数成为表示总体属 性不确定性的一个标准指数，而条件Kappa参数则成为某一类别精度描述的指

13、标。误差矩阵是一个常用的遥感影像分类不确定性描述模 型。为使对分类精度的描述具有代表性，样本大小及采样模 式是两个重要的考虑因素。基于误差矩阵，一系列的误差指 标可以被导出，如用户精度、生产者精度、包含误差、丢失 误差等。这些参数可用于描述某一类或整体分类的精度。然 而，在某些情况下需要知道每一个像素的不确定性。以上的 误差描述指标不能满足此要求。4.内部与外部检验确定属性数据的统计质量方法有三种：即演绎法推论、内部检验和外部检验 (Kennedy Smith ， 1986)。通常用演绎 法推断属性数据质量是利用具有内部或外部检验导出的属 性质量的量测值。内部检验方法是通过比较若干相互独立的

14、重复观测量，其平均值被视为“真值”。在质量控制中，内 部检验的结果是准确性。另一方面，外部检验是通过将量测值与“真值”或可以 写作“真值”的量进行比较。外部检验的结果可以满足用户 的需求(Kennedy Smith ，1986)，但这种检验不能区分开各 种误差源或过程的误差影响，其结果包含了各种误差的影 响。在利用外部检验确定属性的数据质量时，首先要选定一 定检查点，可以选择随机抽样点。为确使每一类别内均有一 定的点被选中为检查点，人们通常建议使用分类随机采样方 式。Hay(1979)建议对总体至少应选择50100个采样点，而对于每一类至少应选择30个采样点。使用外部检验法确定属性数据质量的过

15、程描述如下：(1) 定义一置信水平(例如%),从正态分布表中查出此 表所对应的值，即Z a =。(2) 确定采样个数(N)，例如N = 200。(3) 计算检查点正确分类的百分比(P),例如P=89% 利用以下不等式确定检验精度(Drummond , 1991):X2+X-NP 0对于N = 200以及z a =,有： x 因此可以说在的置信水平、对200个点采样精度为89% 的情况下分类精度为 81% 95%很明显，如果减少检查点的 个数或升高置信水平，确定的分类区间宽度将加大。该方法的一个缺点是整体分类精度有可能被拒绝，尽管对某一类的分类精度是可以接受的。遥感分类结果对于某些应用应具有一个最小的正确分类百分比。在该种情况下假设检验最为适合。预先确定精度 的假设检验是一个经常使用的质量控制方法，接受性采样是 质量控制的一个重要分支。关于统计质量控制的详细讨论见Grant 与 Leaven Worth (1988) 论述。类别属性数据外的另一种数据是连续属性数据。以下讨论连续属性数据不确定性的处理方法。三、误差矩阵的内容与表达遥感影像解释成果的可信度或不确定性问题往往是采用误差矩阵方法进行检验。这是公认的科学方法。北京国土资源遥感