遥感数据预处理

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

2、数据预处理的各个流程介绍
    (一)几何精校正与影像配准
    引起影像几何变形一般分为两大类:系统性和非系统性。系统性一般有传感器本身引起的,有规律可循和可预测性,可以用传感器模型来校正;非系统性几何变形是不规律的,它可以是传感器平台本身的高度、姿态角不稳定,也可以是地球曲率及空气折射的变化以及地形的变化等。
在做几何校正前,先要知道几个概念:
        地理编码:把图像矫正到一种统一标准的坐标系。
        地理参照:借助一组控制点,对一幅图像进行地理坐标的校正。
        图像配准:同一区域里一幅图像(基准图像)对另一幅图像校准

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

这是几何校正中最重要的一步。可以从地形图()为参考进行控制选点,也可以野外测量获得,或者从校正好的影像中获取。选取得控制点有以下特征:
在图像上有明显的、清晰的点位标志,如道路交叉点、河流交叉点等;

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

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

 

)建立几何校正模型
    地面点确定之后,要在图像与图像或地图上分别读出各个控制点在图像上的像元坐标()及其参考图像或地图上的坐标(),这叫需要选择一个合理的坐标变换函数式(即数据校正模型),然后用公式计算每个地面控制点的均方根误差(
      根据公式计算出每个控制点几何校正的精度,计算出累积的总体均方差误差,也叫残余误差,一般控制在一个像元之内,即
)图像重采样
    重新定位后的像元在原图像中分布是不均匀的即输出图像像元点在输入图像中的行列号不是或不全是正数关系。因此需要根据输出图像上的各像元在输入图像中的位置,对原始图像按一定规则重新采样,进行亮度值的插值计算,建立新的图像矩阵。常用的内插方法包括:
、最邻近法是将最邻近的像元值赋予新像元。该方法的优点是输出图像仍然保持原来的像元值,简单,处理速度快。但这种方法最大可产生半个像元的位置偏移,可能造成输出图像中某些地物的不连贯。
、双线性内插法是使用邻近个点的像元值,按照其距内插点的距离赋予不同的权重,进行线性内插。该方法具有平均化的滤波效果,边缘受到平滑作用,而产生一个比较连贯的输出图像,其缺点是破坏了原来的像元值。
、三次卷积内插法较为复杂,它使用内插点周围的个像元值,用三次卷积函数进行内插。这种方法对边缘有所增强,并具有均衡化和清晰化的效果,当它仍然破坏了原来的像元值,且计算量大。
     一般认为最邻近法有利于保持原始图像中的灰级,但对图像中的几何结构损坏较大。后两种方法虽然对像元值有所近似,但也在很大程度上保留图像原有的几何结构,如道路网、水系、地物边界等。

数字图像镶嵌与裁剪

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

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







### 多源遥感数据预处理的方法及工具 #### 数据输入与格式转换 多源遥感数据通常来自不同的传感器,因此第一步是对这些数据进行统一的格式化处理。ENVI软件支持多种常见的遥感影像格式,能够完成数据的读取和写入操作[^1]。通过这一过程,可以确保不同来源的数据能够在同一平台上被有效管理。 #### 几何校正 为了使多源遥感数据在同一地理坐标系下对齐,几何校正是必不可少的一个环节。ENVI提供了完整的几何校正模块,允许用户定义地面控制点并执行多项式变换或其他类型的校正算法。此外,在计算环境中,Google Earth Engine (GEE)也内置了几何校正的功能,可以直接调用已有的地图投影库来简化此流程[^2]。 #### 辐射定标与大气校正 辐射定标的目的是将原始DN值转化为物理单位(如反射率)。这一步骤对于消除成像过程中仪器特性变化的影响至关重要。而随后的大气校则进一步修正由于大气散射等因素引起的误差。在实际应用中,可以通过ENVI中的专用插件或者脚本来实现这两个步骤;而在GEE里,则可利用其预先加载的标准模型快速完成相应调整[^4]。 #### 图像配准与镶嵌 当研究区域跨越多个场景时,可能需要将若干幅单独拍摄的小图拼接起来形成一幅大图。此时就需要先做精确的空间位置匹配即图像配准工作,然后再实施无缝连接——也就是所谓的图像镶嵌。无论是传统桌面级应用程序还是现代端解决方案都能很好地满足这方面的需求[^3]。 #### 波段组合与增强显示 最后,在正式进入分类识别之前往往还需要做一些视觉上的优化以便更好地观察目标对象特征。比如根据不同波长特性的组合创造出假彩色合成效果等等。这类任务同样可以在上述提到的各种软硬件设施之上轻松达成目的。 ```python import ee ee.Initialize() # Example of loading Landsat data and applying basic preprocessing steps using GEE. dataset = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') \ .filterDate('2020-01-01', '2020-12-31') \ .map(lambda image: image.clip(ee.Geometry.Rectangle([79, 26, 80, 27]))) def applyRadiometricCalibration(image): calibratedImage = image.select(['B.*']).multiply(0.0001).rename(['calibrated_B.*']) return calibratedImage.copyProperties(image,['system:time_start']) preprocessedData = dataset.map(applyRadiometricCalibration) print(preprocessedData.first().getInfo()) ```
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值