基于ERDAS变化信息提取方法-李英超..doc

变化信息提取 4.5.1 变化信息自动提取方法 土地利用变化监测（即变化信息自动发现）方法主要有影像相减法、植被指数相减法、变化矢量分析法、主分量分析法、光谱特征变异法、分类结果比较法等。 前5种方法只是检测出可能的变化，而并没有给出土地利用变化的定量信息（如面积）和变化中类型的转化信息（如地类属性）。 分类结果比较法的最终精度受到影像分类精度的限制，而且它对影像的全部范围都要进行分类计算而不管它们是否已经发生变化，增加了变化信息检测的计算量。即使对于同一地物，由于条件的不同，得到的影像灰度值也不太一样。因此对于影像相减法而言，单纯相减所得的变化模板中肯定会含有大量的假变化信息和噪声信息，要从这些信息中提取出真正的变化仍旧是个棘手问题。由于异物同谱现象的存在，许多真正的变化信息也会因为相减而被漏掉，从而影响了最终变化信息的获得。植被指数相减的方法同样有着本身的局限性，它对不同时相植被覆盖情况的变化敏感，而不能很好地发现其他类型变化。 4.5.2不同时相遥感影像变化信息发现 在没有土地利用基础图件的情况下，利用两个不同时相或序列不同时相的遥感影像进行变化信息的发现方法主要包括：光谱特征变异法、主成分分析法、假彩色合成法、图像差值法、分类后比较法、波段替换法及变化矢量分析法。下面对几种比较成熟的方法做简单介绍。 （一）光谱特征变异法 同一地物反映在SPOT影像上的信息与其反映在TM影像上的光谱信息是一一对应的。因此对同一时相的TM和SPOT影像进行融合后，地物光谱属性可以如实正确地表现出来。但如果同一地物在两者上的信息表现为不一致时，那么融合后影像中此地物的光谱就表现得与正常地物的光谱有所差别，此时称地物发生了光谱特征变异，我们也就可以根据发生变异的光谱特征确定变化信息的发生。 发生变化的区域在融合后的影像上表现为含有纹理特征的亮绿色，区别于周围植被地物的绿色特征和正常的居民地信息，判定为该区域新增了一块建设用地。 图4-5-1光谱特征变异法技术流程 图4-5-2光谱特征变异法 特点：物理意义明显、操作简捷的特点。但是发生光谱特征变异的地物在几何尺寸上要足够大才能被人工目视发现，而小于这一尺寸的地物变化就很难被目视解译出来。 假彩色合成法 由于地表的变化，相同传感器对同一地点所获取的不同时相的影像在灰度上有较大的区别。在进行变化信息的发现时，将前、后两时相的全色SPOT数据精确配准，再利用假彩色合成的方法，将后一时相的数据赋予红色通道，前一时相的数据赋予蓝色和绿色通道。利用三原色原理，形成假彩色影像。其中，地表未发生变化的区域，合成后影像灰度值接近，而土地利用发生了变化的区域则呈现出红色，即判定为变化信息的所在。 图4-5-3假彩色合成法 多波段主成分分析法 主成分分析法的数学基础是主成分变换，即基于变量之间的相互关系，在尽量不丢失 信息的前提下，利用线性变换的方法实现数据压缩和特征选取的目的。它的基本原理在第 三章中第二节“四融合方法的实现”中有了较为详细的介绍，这里不在重述。在土地利用 动态监测中，根据操作方法的不同，该方法又有以下几种方式： 主成分差异法:将两时相的融合影像经纠正、配准后，分别对它们进行 PC 变换。根据 PCA 的特性，变化结果的第一分量集中了影像的主要信息，其他分量则反映了影像的差异信息。因此将 两时相的第一主分量作差值取绝对值处理，生成的差值影像即为变化信息模板。 图4-5-4差异主成分技术流程 图4-5-5差异主成分法 多波段主成分分析法：由遥感理论得知，地物属性发生变化，必将导致其在影像某几个波段上的值发生变化， 所以只要找出两时相影像中对应波段上值的差别并确定这些差别的范围，便可发现土地利 用变化信息。在具体工作中将两时相的影像各波段进行组合成一个两倍于原影像波段数的新影像，对该影像作 PC 变换。由于变换结果前几个分量上集中了两个影像的主要信息，而后几个分量则反映出了两影像的差别信息， 因此可以试着抽取后几个分量进行波段组合 来产生变化信息。例如对以 TM5、4、3 波段组合并与 SPOT 融合的影像来说，把两时相的 影像进行波段合成，新产生的影像将有 6 个波段的信息。对这个影 像作 PC 变换的结果含有6个分量的影像值，这6个分量的前一、二分量便包含了原始影像的主要信息，而后几个 分量则反映了两时相影像的变化信息。一般来讲，在上述多波段主成分变 换之后，采用6、4、3 分量进行波段组合能较好地反映出前后时相影像的变化部分。在所示的例子中，变化信息在结果影像中呈现为蓝色。 多波段主成分分析方法的不足在于操作过程较复杂，数据量和计算量都较大，运行速度 较慢。而且，主分量选择的合理性难以确定，存在信息遗失的情况。 图4-5-6多波段主成分技术