完整版遥感影像预处理.docx

1、完整版遥感影像预处理遥感影像预处理预处理是遥感应用的第一步， 也是非常重要的一步。 目前的技术也非常成熟， 大多数的商业化软件都具备这方面的功能。 预处理的大致流程在各个行业中有点差异，而且注重点也各有不同。本小节包括以下内容：数据预处理一般流程介绍预处理常见名词解释ENVI 中的数据预处理1、数据预处理一般流程数据预处理的过程包括几何精校正、配准、图像镶嵌与裁剪、去云及阴影处理和光谱归一化几个环节，具体流程图如图所示。图 1 数据预处理一般流程各个行业应用会有所不同， 比如在精细农业方面， 在大气校正方面要求会高点， 因为它需要反演；在测绘方面，对几何校正的精度要求会很高。2、数据预处理的各

2、个流程介绍（一）几何精校正与影像配准引起影像几何变形一般分为两大类： 系统性和非系统性。 系统性一般有传感 器本身引起的， 有规律可循和可预测性， 可以用传感器模型来校正； 非系统性几 何变形是不规律的， 它可以是传感器平台本身的高度、 姿态等不稳定， 也可以是 地球曲率及空气折射的变化以及地形的变化等。在做几何校正前，先要知道几个概念：地理编码：把图像矫正到一种统一标准的坐标系。地理参照：借助一组控制点，对一幅图像进行地理坐标的校正。图像配准：同一区域里一幅图像（基准图像）对另一幅图像校准影像几何精校正，一般步骤如下，（1）GCP（地面控制点）的选取这是几何校正中最重要的一步。可以从地形图（

3、 DRG）为参考进行控制选 点，也可以野外 GPS 测量获得，或者从校正好的影像中获取。选取得控制点有 以下特征：1、GCP 在图像上有明显的、清晰的点位标志，如道路交叉点、河流交叉点 等；2、地面控制点上的地物不随时间而变化。GCP均匀分布在整幅影像内，且要有一定的数量保证，不同纠正模型对控制点 个数的需求不相同。卫星提供的辅助数据可建立严密的物理模型，该模型只需 9 个控制点即可；对于有理多项式模型，一般每景要求不少于 30 个控制点，困难 地区适当增加点位； 几何多项式模型将根据地形情况确定， 它要求控制点个数多 于上述几种模型，通常每景要求在 30-50 个左右，尤其对于山区应适当增加

4、控 制点。（2）建立几何校正模型 地面点确定之后，要在图像与图像或地图上分别读出各个控制点在图像上的像元坐标（ x，y）及其参考图像或地图上的坐标（ X，Y），这叫需要选择一个 合理的坐标变换函数式 （即数据校正模型） ，然后用公式计算每个地面控制点的 均方根误差（ RMS ）根据公式计算出每个控制点几何校正的精度， 计算出累积的总体均方差误差，也叫残余误差，一般控制在一个像元之内，即 RMS1 。3）图像重采样重新定位后的像元在原图像中分布是不均匀的 ,即输出图像像元点在输入图 像中的行列号不是或不全是正数关系。 因此需要根据输出图像上的各像元在输入 图像中的位置， 对原始图像按一定规则重新

5、采样， 进行亮度值的插值计算， 建立 新的图像矩阵。常用的内插方法包括：1、最邻近法是将最邻近的像元值赋予新像元。 该方法的优点是输出图 像仍然保持原来的像元值， 简单， 处理速度快。 但这种方法最大可产生半个像元 的位置偏移，可能造成输出图像中某些地物的不连贯。2、双线性内插法是使用邻近 4 个点的像元值，按照其距内插点的距 离赋予不同的权重， 进行线性内插。 该方法具有平均化的滤波效果， 边缘受到平 滑作用，而产生一个比较连贯的输出图像，其缺点是破坏了原来的像元值。3、三次卷积内插法较为复杂，它使用内插点周围的 16 个像元值，用 三次卷积函数进行内插。 这种方法对边缘有所增强， 并具有均

6、衡化和清晰化的效 果，当它仍然破坏了原来的像元值，且计算量大。一般认为最邻近法有利于保持原始图像中的灰级，但对图像中的几何结构 损坏较大。 后两种方法虽然对像元值有所近似， 但也在很大程度上保留图像原有 的几何结构，如道路网、水系、地物边界等。（二）数字图像镶嵌与裁剪l 镶嵌当研究区超出单幅遥感图像所覆盖的范围时， 通常需要将两幅或多幅图像拼 接起来形成一幅或一系列覆盖全区的较大的图像。在进行图像的镶嵌时， 需要确定一幅参考影像， 参考图像将作为输出镶嵌图像的 基准，决定镶嵌图像的对比度匹配、 以及输出图像的像元大小和数据类型等。 镶 嵌得两幅或多幅图像选择相同或相近的成像时间， 使得图像的色

7、调保持一致。 但 接边色调相差太大时， 可以利用直方图均衡、 色彩平滑等使得接边尽量一致， 但 用于变化信息提取时，相邻影像的色调不允许平滑，避免信息变异。l 裁剪图像裁剪的目的是将研究之外的区域去除， 常用的是按照行政区划边界或自然区 划边界进行图像的分幅裁剪。它的过程可分为两步：矢量栅格化和掩膜计算 （Mask ）。矢量栅格化是将面状矢量数据转化成二值栅格图像文件，文件像元 大小与被裁剪图像一致；把二值图像中的裁剪区域的值设为 1 ，区域外取 0 值，与被裁剪图像做交集运算，计算所得图像就是图像裁剪结果。（三）大气校正遥感图像在获取过程中， 受到如大气吸收与散射、 传感器定标、 地形等因素

8、 的影响，且它们会随时间的不同而有所差异。因此，在多时相遥感图像中，除了 地物的变化会引起图像中辐射值的变化外， 不变的地物在不同时相图像中的辐射 值也会有差异。利用多时相遥感图像的光谱信息来检测地物变化状况的动态监测， 其重要前提是要消除不变地物的辐射值差异辐射校正是消除非地物变化所造成的图像辐射值改变的有效方法， 按照校正 后的结果可以分为 2 种，绝对辐射校正方法和相对辐射校正方法。 绝对辐射校正 方法是将遥感图像的 DN（Digital Number） 值转换为真实地表反射率的方法， 它 需要获取影像过境时的地表测量数据， 并考虑地形起伏等因素来校正大气和传感 器的影响， 因此这类方法

9、一般都很复杂， 目前大多数遥感图像都无法满足上述条 件。相对辐射校正是将一图像作为参考 （或基准）图像，调整另一图像的 DN 值， 使得两时相影像上同名的地物具有相同的 DN 值，这个过程也叫多时相遥感图像 的光谱归一化。这样我们就可以通过分析不同时相遥感图像上的辐射值差异来实 现变化监测。 因此，相对辐射校正就是要使相对稳定的同名地物的辐射值在不同 时相遥感图像上一致，从而完成地物动态变化的遥感动态监测。3、ENVI 中的数据预处理介绍（一）几何精校正与影像配准（1）选择几何校正模型ENVI 中支持有大多数商业化卫星的几何校正模型， 如 QuickBird 、Ikonos 、 Spot1-5

10、 、 P6、 WorldView-1 等，一般的校正模型包括二次多项式、仿射变换 和局部三角网图 2 何校正模型控制点选择方式可以是从影像上，也可以从矢量数据或者野外实测等。图 3 控制点选择方式选择控制点也非常的方便，包含了误差的结算重采样方式包含了三种方法。图 4 控制点选择图 5 重采样方式（二）数字图像镶嵌与裁剪l 镶嵌ENVI 支持有地理参照和没有地理参照影像数据的镶嵌，能够自动对镶嵌影 像进行颜色平衡， 并提供了多种影像增强和直方图匹配工具， 可以最大限度地消 除镶嵌影像间的色调和颜色差异多种色彩平衡方法图 6 颜色校正设置多种接边线编辑方式图 7 接边线镶嵌l 裁减在 ENVI 中做裁减的方法非常的多，提供多种方法进行图像的空间裁剪获得子区，包括：手动输入行列数、从图像中交互选择区域、输入地理坐标范围、和 另外图像文件的交集、使用滚动窗口中的图像和通过感兴趣区域图 8 影像的裁剪三）大气校正ENVI 的大气校正模块为 FLAASH 。详细情况参见帖子： ENVI 中的大气校正模块（ FLAASH ）的使用说明