7 H6 N6 V4 T/ ^' T9 q: ]" m+ W3 _
可快速关注我们 - e% R! R! U% W, a* j/ b
地理国情是从空间角度分析呈现国情信息,包括森林植被覆盖信息、人口及土地覆盖变化等。地理国情监测就是使用一些技术手段对监测区域进行观测,获取信息的现状并分析变化趋势。因此,国家的持续稳定发展离不开地理国情监测,监测结果对我国重大战略决策有着至关重要的作用。
6 W! x+ l4 L+ {8 R# o3 @, N 在早期的地理国情监测中,通常采用航天遥感数据作为主要数据源,以LANDSAT为代表的中低分辨率卫星,其影像分辨率低、辨析度较差,信息的判读和提取存在一定的困难,在一定程度上影响了地理国情监测的数据精度。为了满足日益增长的国情需要,高分辨率卫星(如高分一号卫星) 逐渐成为遥感数据来源的主流。 # C% A$ h3 C5 J! x
1.高分一号卫星
; f+ p5 d' b, D7 k# ]) T% W# w$ @ 高分一号卫星属于低轨卫星,于2013年4月成功发射。在单星上同时实现了大宽幅与高分辨率的结合:16m分辨率实现大于800km成像幅宽;2m高分辨率实现大于60km成像幅宽。适应多源遥感数据、多种高分辨率等综合需求。同时,实现了在无地面控制点的条件下达到50m图像的定位精度,及2×450 Mbps的数据传输能力,满足大数据量应用需求,大大提高了数据的传输能力,在国内同类卫星中属于最高水平。高分一号卫星是中国高分辨率遥感卫星技术历程的里程碑,标志着中国遥感卫星领域正式跨入了高分时代,推动了我国现代测绘技术的发展。 " t0 s3 {0 s. e q
2.原始遥感影像资料选取 ! W; h& J( x9 O$ m
本文所采用的高分卫星数据来自于我国自主研发设计的高分一号卫星,分辨率为8m4个波段的多光谱影像和2m的单波段全色影像。在影像处理前,需要对原始影像进行检查和筛选,筛选出的数据必须符合地理国情监测中土地覆盖变化研究的需要。
( Z; c) e" A- d3 U C 由于研究区域冬季时间长,冰雪天气较多,根据研究区域独特的气候特征,选取2013年7-10月的10景影像、2014年7-10月的10景影像作为原始数据。在数据准备阶段需对卫星影像数据进行初步检测,筛选出满足检测研究的影像数据。首先,对影像的重叠度进行检测,运用ENVI 5.1依次打开同年份的所有影像查看重叠度,根据监测研究的要求尽量选择重合度在5%-10%左右的影像数据,排除无重叠的影像;其次,对影像上云雪覆盖情况进行目视判别,尽管研究区内7-10月的气候条件较利于观测,但影像拍摄过程中偶尔还会存在云朵的影响,故要求研究影像中云朵覆盖率小于10%;最后,还需要对影像的质量进行粗略查看,检查影像数据是否存在噪声、黑点以及影像光谱信息的丰富度是否满足地理国情监测的条件。 ) [9 q3 O2 T0 ~( `
3.高分一号卫星数据处理 9 ~/ a8 C) M/ ~! M# P! r8 |
在对原始影像数据进行筛选检查结束后,筛选出的影像基本为可用数据。为了提供更加精准的变化监测数据,需要对这些数据进行相应的处理。在影像处理阶段使用ENVI 5.1对部分影像进行试验,制定出较完善的影像处理流程(如图1所示) 。
- _- l3 }. U$ S( R
, G& H6 z0 U& D9 w% F. i! o 2 w6 {& j4 d7 p% n& ~) h
图1 影像处理流程图3.1 辐射定标 9 g: S; F6 b0 Z/ [) q* S, V. j
辐射定标是遥感信息定量化的前提,是遥感数据可靠性及可应用性的基础保障。辐射定标分为绝对定标和相对定标,相对定标又称为传感器探测元件归一化,是为了校正传感器中各个探测元件响应度差异而对卫星传感器测量到的原始亮度值进行归一化的一种处理过程。 9 p. b L: S1 A; |& A/ L" {
因原始影像并未做过辐射定标处理,故使用ENVI 5.1对多光谱影像和全色影像进行相对辐射定标处理。图2为多光谱影像经过辐射定标的前后对比图。 6 V* T* t# Q5 [' V. n
) U" s; f2 r3 F. ?+ a$ _ ( S! h+ L. n# ?, t# K4 k; T1 T
图2 多光谱影像辐射定标前后对比图经过相对辐射定标后的影像亮度发生了变化。色调上经过辐射定标后影像对比度更强,地物辨析程度更高。通过查看波谱曲线,辐射定标后的数值主要集中在0-10范围内,图3为多光谱影像辐射定标后波谱曲线图,单位μW/( cm2×sr×nm)。
. i" |8 ]2 j+ |1 |- Z( X ; C1 V5 Y2 F' ~4 b6 ^1 y9 N6 P
5 S0 L$ D+ N! g* w 图3 多光谱影像辐射定标后波谱曲线图3.2 大气校正 ! }3 f( ?2 @6 t! o9 x4 a2 O
大气校正是遥感数据处理的重要部分,随着遥感技术的发展,利用遥感数据对土地覆盖、土地利用、气候变化等国情信息进行变化监测的要求也越来越高。电磁波透过大气层时,不仅改变了光线的方向,也会影响遥感图像的辐射特征。大气校正的目的在于通过图像处理减弱大气散射、吸收等引起的误差,使图像数据更加精确。 ( l& h4 h3 M e3 y
本研究所用影像均进行了FLAASH大气校正,FLAASH大气校正采用了MODTRAN 4+辐射传输模型。任何有关影像标准的MODTRAN大气校正模型和气溶胶类型都可以直接使用。FLAASH大气校正的优势在于:可以通过影像像素光谱上的特征来估计大气的属性,不依赖遥感成像时同步测量的大气参数数据。同时,有效去除气溶胶散射效应(水蒸气),并基于像素级校正目标像元和邻近像元交叉辐射的“邻近效应”,对由于人为抑止而导致的波谱噪声进行光谱平滑处理。 3 L8 s2 T$ p# l# L$ w$ b" o
对经过辐射定标后的多光谱影像进行FLAASH大气校正,其成像中心点经纬度FLAASH自动从影像中获取,传感器高度根据规定为645km,像元大小为8m,成像时间需要在真实时间基础上减去8h换算GMT时间,大气模型和气溶胶模型可根据影像区域类别进行选择。完成设置后存储文件输出,查看大气校正前后影像对比。图4为FLASSH校正前后对比图。 2 \6 H1 |/ l; q0 Q; W
) y& V. k6 a, x' o5 k: `! s * t& u# Z, M0 C! S! j
图4 FLAASH校正前后对比图从图4可以看出,经过大气校正后,图像的对比度增强了,同时去除了部分雾的干扰,视觉改善效果虽不显著,但像元值发生变化,校正后的水体波谱曲线更加接近真实水体波谱曲线。
9 r% @- b' h( o1 R( A% ?2 a2 j# A 3.3 正射校正
$ g- x3 m" M9 m- L9 y/ Y$ b 遥感影像在拍摄过程中,由于飞行器拍摄时的姿态、飞行速度、轨道高度以及地球自转等因素的影响,使得获取的图像相对于地面目标发生几何畸变,其表现为像元相对于地面目标的实际位置发生挤压、扭曲、拉伸和偏移等现象,通常要对这种几何畸变进行几何校正。几何校正是指通过数学模型改正和消除遥感影像成像时产生的原始图像上各地物的几何位置、形状、尺寸、方位等特征的变形。本研究主要采用的是正射校正。高分一号的L1A级包括了RPC文件,在经过辐射定标、大气校正等处理后,ENVI自动将RPC嵌入到处理结果中,进行基于无控制点的对多光谱和全色数据结果进行正射校正。图5为多光谱影像正射校正前后对比图。
& q+ B0 T1 X: K. J
3 |9 ?+ n; C9 ~ {5 Z6 J7 \ ' g7 p" N* V6 ?. X$ A
图5 多光谱影像正射校正前后对比图如图5所示,研究区域为山丘地带,地形起伏并不十分明显。从多光谱影像中能看出经过正射校正后影像的清晰度明显提高。
, s5 ^/ z7 k3 Y" E+ _ 3.4 影像融合
) [+ G: F! @+ W' j 高分一号卫星数据包括分辨率为8m的多光谱影像和分辨率为2m的全色影像,通过影像融合的方式生成高分辨率多光谱影像,新生成的高分辨率影像既保留了光谱特征,又具有高空间分辨率,满足了在影像上提取足够监测信息的需要。本研究使用Gram-schmidt算法的图像融合方法,该方法融合的影像保真度较好,且操作简单,通过改变多光谱影像存储顺序可以提高融合速率。为了验证此方法的可行性,进行了试验比较,首先不改变正射校正后影像的存储顺序进行融合,大约需要40min完成影像融合,分辨率为2m,融合影像质量较好,轮廓色调清晰; 再将同一幅经过正射校正后的多光谱影像存储为BIP模式后进行融合,大约需要15min完成影像融合,影像质量与未改变存储顺序融合的影像一致。 * M& m+ c1 c9 @( _
经融合后的影像在分辨率上明显优于融合前的多光谱影像,融合后的影像轮廓也变得更加清晰,色调明显优于融合前影像。地物特征更加鲜明,便于后期影像解译标志的建立及分类信息提取。
+ z. h8 _$ q/ Z9 a 3.5 影像镶嵌
* v, n8 o3 i! g 在卫星拍摄过程中,因相机宽幅有限,数据采集区域具有一定的局限性,因此使用的原始数据必须为具有一定重叠度的影像。本研究需要基于影像提取土地覆盖信息,为满足样本选择标准的统一性及土地覆盖变化研究的需要,对已有多景影像进行拼接。本研究使用ENVI 5.1对影像进行镶嵌,使用软件中无缝拼接工具Seamless Mosaic,对影像进行无缝拼接,并自动增强色调。
2 q& b; ~4 D6 Y) _ 在拼接过程中部分影像背景为黑色,需要对背景值赋值为0,从而消除在镶嵌过程中背景的影响。 7 _; }% g8 V1 ^
3.6 影像裁剪 . }/ a/ v: Z7 Q8 M. m3 b
在进行土地覆盖变化监测之前需要对研究区域影像进行裁剪,以保证数据范围一致、面积大小相等,进而保证结果的科学性和有效性。根据已有影像数据情况,在处理过程中按照行政边界对影像进行裁剪,获取符合监测需要的影像区域。
0 c* ~- }% G8 z% D1 [9 w% L 在裁剪过程中,首先将研究区域边界矢量化,由于高分一号卫星影像和研究区域边界矢量图的坐标系统不一致,不能直接进行裁剪,需先对边界矢量图进行投影转换,再根据转换后的矢量边界对影像进行裁剪,获取研究区域影像图。
- B+ D" M. F7 C. r% p; J 4.结束语
7 _% \; v- b9 I4 g 本文通过具体案例总结出高分一号卫星数据处理的关键要点,制定出高分一号卫星数据的处理流程,并获得符合地理国情监测中土地覆盖变化监测所需要的影像数据,进一步拓宽了高分一号卫星影像数据的应用范围,为其在地理测绘、海洋和气候气象观测、水利和林业资源监测、城市和交通精细化管理、疫情评估与公共卫生应急、地球系统科学研究等领域的应用提供了借鉴。 ) B4 n5 J* r6 q" \* m$ W
作者:于爱洁、李桂芬、张佳彬 9 A4 B, Y5 a4 K1 Y) u) |
来源:测绘与空间地理信息 6 n% u: ^' L9 Z- ]0 O
原标题:《行业应用|地理国情监测中高分一号卫星数据处理关键要点分析》
3 F/ r4 H" j6 G5 g. H' O- [' k K+ P T0 U6 u; o% z
; q8 J% {$ Q# n( W
! U% f$ T6 x$ R; J. B1 x4 r$ a
* \- Q1 @# t+ K" q' r, X, G A$ f5 E |