遥感讲座——遥感影像预处理

遥感讲座——遥感影像预处理

据预处理是遥感应用的第一步，也是非常重要的一步。目前的技术也非常成熟，大多数的商业化软件都具备这方面的功能。预处理的大致流程在各个行业中有点差异，而且注重点也各有不同。下面是预处理中比较常见的流程。

1、数据预处理一般流程

数据预处理的过程包括几何精校正、配准、图像镶嵌与裁剪、去云及阴影处理和光谱归一化几个环节，具体流程图如图所示。

各个行业应用会有所不同，比如在精细农业方面，在大气校正方面要求会高点，因为它需要反演；在测绘方面，对几何校正的精度要求会很高。

2、数据预处理的各个流程介绍

（一）几何精校正与影像配准

引起影像几何变形一般分为两大类：系统性和非系统性。系统性一般有传感器本身引起的，有规律可循和可预测性，可以用传感器模型来校正；非系统性几何变形是不规律的，它可以是传感器平台本身的高度、姿态等不稳定，也可以是地球曲率及空气折射的变化以及地形的变化等。

在做几何校正前，先要知道几个概念：

地理编码：把图像矫正到一种统一标准的坐标系。

地理参照：借助一组控制点，对一幅图像进行地理坐标的校正。

图像配准：同一区域里一幅图像（基准图像）对另一幅图像校准

影像几何精校正，一般步骤如下，

（1）GCP（地面控制点）的选取

这是几何校正中最重要的一步。可以从地形图（DRG）为参考进行控制选点，也可以野外GPS测量获得，或者从校正好的影像中获取。选取得控制点有以下特征：

1、GCP在图像上有明显的、清晰的点位标志，如道路交叉点、河流交叉点等；

2、地面控制点上的地物不随时间而变化。

GCP均匀分布在整幅影像内，且要有一定的数量保证，不同纠正模型对控制点个数的需求不相同。卫星提供的辅助数据可建立严密的物理模型，该模型只需9个控制点即可；对于有理多项式模型，一般每景要求不少于30个控制点，困难地区适当增加点位；几何多项式模型将根据地形情况确定，它要求控制点个数多于上述几种模型，通常每景要求在30-50个左右，尤其对于山区应适当增加控制点。

（2）建立几何校正模型

地面点确定之后，要在图像与图像或地图上分别读出各个控制点在图像上的像元坐标（x，y）及其参考图像或地图上的坐标（X，Y），这叫需要选择一个合理的坐标变换函数式（即数据校正模型），然后用公式计算每个地面控制点的均方根误差（RMS）根据公式计算出每个控制点几何校正的精度，计算出累积的总体均方差误差，也叫残余误差，一般控制在一个像元之内，即RMS<1。

（3）图像重采样

重新定位后的像元在原图像中分布是不均匀的,即输出图像像元点在输入图像中的行列号不是或不全是正数关系。因此需要根据输出图像上的各像元在输入图像中的位置，对原始图像按一定规则重新采样，进行亮度值的插值计算，建立新的图像矩阵。常用的内插方法包括：

1、最邻近法是将最邻近的像元值赋予新像元。该方法的优点是输出图像仍然保持原来的像元值，简单，处理速度快。但这种方法最大可产生半个像元的位置偏移，可能造成输出图像中某些地物的不连贯。

2、双线性内插法是使用邻近4个点的像元值，按照其距内插点的距离赋予不同的权重，进行线性内插。该方法具有平均化的滤波效果，边缘受到平滑作用，而产生一个比较连贯的输出图像，其缺点是破坏了原来的像元值。

3、三次卷积内插法较为复杂，它使用内插点周围的16个像元值，用三次卷积函数进行内插。这种方法对边缘有所增强，并具有均衡化和清晰化的效果，当它仍然破坏了原来的像元值，且计算量大。

一般认为最邻近法有利于保持原始图像中的灰级，但对图像中的几何结构损坏较大。后两种方法虽然对像元值有所近似，但也在很大程度上保留图像原有的几何结构，如道路网、水系、地物边界等。

（二）数字图像镶嵌与裁剪

镶嵌

当研究区超出单幅遥感图像所覆盖的范围时，通常需要将两幅或多幅图像拼接起来形成一幅或一系列覆盖全区的较大的图像。

在进行图像的镶嵌时，需要确定一幅参考影像，参考图像将作为输出镶嵌图像的基准，决定镶嵌图像的对比度匹配、以及输出图像的像元大小和数据类型等。镶嵌得两幅或多幅图像选择相同或相近的成像时间，使得图像的色调保持一致。但接边色调相差太大时，可以利用直方图均衡、色彩平滑等使得接边尽量一致，但用于变化信息提取时，相邻影像的色调不允许平滑，避免信息变异。

裁剪

图像裁剪的目的是将研究之外的区域去除，常用的是按照行政区划边界或自然区划边界进行图像的分幅裁剪。它的过程可分为两步：矢量栅格化和掩膜计算（Mask）。矢量栅格化是将面状矢量数据转化成二值栅格图像文件，文件像元大小与被裁剪图像一致；把二值图像中的裁剪区域的值设为1，区域外取0值，与被裁剪图像做交集运算，计算所得图像就是图像裁剪结果。

（三）大气校正

遥感图像在获取过程中，受到如大气吸收与散射、传感器定标、地形等因素的影响，且它们会随时间的不同而有所差异。因此，在多时相遥感图像中，除了地物的变化会引起图像中辐射值的变化外，不变的地物在不同时相图像中的辐射值也会有差异。利用多时相遥感图像的光谱信息来检测地物变化状况的动态监测，其重要前提是要消除不变地物的辐射值差异。

辐射校正是消除非地物变化所造成的图像辐射值改变的有效方法，按照校正后的结果可以分为2种，绝对辐射校正方法和相对辐射校正方法。绝对辐射校正方法是将遥感图像的DN(Digital Number)值转换为真实地表反射率的方法，它需要获取影像过境时的地表测量数据，并考虑地形起伏等因素来校正大气和传感器的影响，因此这类方法一般都很复杂，目前大多数遥感图像都无法满足上述条件。相对辐射校正是将一图像作为参考（或基准）图像，调整另一图像的DN值，使得两时相影像上同名的地物具有相同的DN值，这个过程也叫多时相遥感图像的光谱归一化。这样我们就可以通过分析不同时相遥感图像上的辐射值差异来实现变化监测。因此，相对辐射校正就是要使相对稳定的同名地物的辐射值在不同时相遥感图像上一致，从而完成地物动态变化的遥感动态监测。

3、ERDAS中的数据预处理介绍

（一）几何精校正与影像配准

（1）选择几何校正模型

ERDAS有大多数商业化卫星的几何校正模型，如Landsat、QuickBird、Ikonos、Spot4、5、P5等。

控制点选择方式也有很多。

选择控制点也非常的方便，包含了误差的结算。

重采样方式包含了四种方法。

（二）数字图像镶嵌与裁剪 镶嵌

在ERDAS IMAIGEN V9.x中，镶嵌模块（mosaic）新增加了三种方式，除了之前的mosaic Tool,新添Mosaic pro、Mosaic Direct和Mosaic Wizard。

ERDAS的Mosaic模块功能比较全，包括了色彩平衡、接边线编辑、输出区域裁减等。操作也非常简单，流程化操作方式。

多种色彩平衡方法

多种接边线编辑方式。

裁减

在ERDAS中做裁减的方法非常的多，如Dice Image模块，Subset Image，Mask 模块，Mosaic下的裁减，Map Series Tool裁减。

如下图为Dice Image模块下的地图裁减。打开DatePre→Dice Image模块。输入分幅影像文件和输出的文件，设置划分间隔，如需要还可以设置重叠区宽度。

在Mosaic中也可以利用矢量文件对影像进行裁减。打开DatePre→Mosaic Images →Mosaic Tool模块。加载完分幅影像之后，在Set mode for Output Image下，点开Set Output Option dialog对话框，如图4。在method列表框下有6个选项，选择Polygon Vector File，之后选择一个分幅的矢量文件，同时还可以选择其中的一个字段作为分幅结果的文件名。

（三）大气校正

ERDAS的大气校正模块为ATCOR。详细情况参见帖子：https://www.360docs.net/doc/b93361416.html,/ESRI/viewthread.php?tid=11818&highlight=%B4% F3%C6%F8%D0%A3%D5%FD。