利用历史卫星影像重建南极Byrd冰川流速图

1、利用历史卫星影像重建南极Byrd冰川流速图【摘要】随着全球气候变暖，南极冰川流速研究成为海平面上升模型的一个关键参数。本文利用20世纪80年代的卫星影像，通过影像预处理及密集匹配，重建了南极Byrd冰川区域的历史流速图，为监测该区域的冰流速以及物质平衡提供了根底。【关键词】南极;冰川流速;Landsat中图分类号：P228.3;P223.4文献标识码：A文章编号：2095-2457202111-0001-002【Abstract】Withtheglobalwarming，velocityofglaciersinAntarctichasbecomeavitalparameteroftheseal

2、evelrisemodel.SatelliteimageshavebeenusedinthispapertorebuildthehistoricaliceflowmapofByrdGlacierinAntarcticbyimagepreprocessinganddensematching，whichprovidesafoundationformonitoringtheicevelocityandmassbalanceinthisarea.【Keywords】Antarctic;Glaciervelocity;Landsat0引言南极作为地球的冷源之一，是全球气候变化研究中最为关键和敏感的区域，

3、其物质平衡变化对全球海平面和气候变化具有重要的影响【1】。政府间气候变化专门委员会第五次评估报告指出，全球平均海平面在19012021年间上升的平均速率为1.7mm/y。南极冰层融化与全球海平面变化的定量关系及其影响机理，以及未来海平面变化的趋势仍然存在相当大的不确定性【2】。而冰川流速作为南极冰盖稳定性的一个指示器，是影响南极冰川物质平衡的一个重要因素，因此测定南极冰川流速在全球气候变化精密定量化中扮演了非常重要的角色。目前，南极冰流速监测主要分为两类：野外实地测量和遥感监测法。野外实地测量利用花杆、GPS以及全站仪、经纬仪和激光测距仪等测绘仪器3-6，通过在实地直接测量地面点的位移得到该区

4、域的冰流速。野外实地测量的突出优点是精度很高，但由于南极恶劣的气候条件以及特殊的地质构造大大限制了实地测量的应用，无法进行大范围、快速、持续的测量。自从20世纪70年代第一颗Landsat卫星发射成功以来，学者们便开始系统地对南极冰流速进行研究【7】。Rignot等综合使用了1996年至2021年的SAR数据，首次绘制了全南极的外表流速图8。由于南极历史冰流速数据的缺失，导致监测南极冰流速变化困难重重。本文利用80年代的Landsat-4及Landsat-5影像对Byrd冰川进行遥感监测，并获取了Byrd区域80年代的冰流速数据，使得对南极进行长时间的冰流速监测成为可能。1研究区域及数据1.1

5、研究区域Byrd冰川是南极洲的一个重要的冰川，长约136km，宽约24km，平均流速为650m/y9。同时Byrd冰川位于南极洲最大的集水区，每年流入Ross冰架的物质体积为23.62km310。1.2卫星影像数据本文使用Landsat-4以及Landsat-5的光学影像数据的获取时间、分辨率、波段以及波长如表1所示。2研究方法2.1高斯高通濾波由于冰流区底部基岩会引起冰流外表起伏，进而影响匹配结果，所以本文采用高斯高通滤波将这种由底部基岩引起的起伏剔除，同时增强影像外表的纹理。高斯高通滤波首先对输入的影像进行中心变化并计算图像的傅里叶变换DFT，即Fu，v，然后用滤波器传递函数Hu，v乘以F

6、u，v并计算傅里叶反变换，最后把图像由频率域返回到空间域中。其中截止频率处在距频率中心距离为D0的高斯高通滤波器GHPF的传递函数由下式给出：D0为截止频率，Du，v是点u，v到傅里叶频率域原点的距离，Du，v=u-M/22+v-N/221/2，n为滤波器的阶数。2.2归一化互相关匹配通过光学影像进行冰川流速的监测主要采用特征跟踪法，而提取完影像上的特征点之后，那么需要利用互相关的方法来衡量两组特征点之间的相似度。归一化互相关是一种根底且有效的检测影像同名点的方法，归一化互相关算法在视觉比照度较差的区域、薄云区域以及高视觉比照度区域中表现良好11。对于一对连续影像，将时间较早的图像作为参考影像

7、，另一张影像那么作为待搜索影像。对于影像上面的特征点P，首先以P为中心在影像中确定一个矩形作为参考窗口，然后在影像中以P的同名点P为中心也开辟一个窗口作为搜索窗口，以公式2计算两者之间的相关系数12。两次匹配保存相关系数r最高的一对点并认为是候选同名点，只有当两次匹配结果一致时保存这对同名点。其中a，b是匹配窗口中心相对于参考窗口中心的坐标偏移量，fxi，yj和gxi，yj是参考窗口和匹配窗口xi，yj处的灰度值，f和g是参考窗口和匹配窗口的平均灰度值，m和n表示参考窗口和匹配窗口的长宽。3实验过程本实验中所采用的卫星影像来自于Landsat-4和Landsat-5，而这两个卫星所拍摄的影像定

8、位精度较低，需要进行正射校正。本文通过选取7个人工控制点利用Geomatica软件中的正射模块来对原始影像进行正射校正，经过处理所得到的正射影像的定位精度有了明显的提高，均小于一个像素表2。在正射影像的根底上，利用NCC匹配法进行影像匹配。在得到原始匹配点之后，根据相关系数以及冰流运动理论进行误匹配点的剔除，最终将剩余特征点插值得到Byrd冰川区域的冰流图。如图1所示，以Byrd冰川地区的LIMA影像为底图，展示了冰川的流速，主冰流区域的流速在200-870m/y。4结语本文利用南极区域的Landsat历史影像，通过归一化互相关匹配法以及高斯高通滤波得到同名像对的匹配点，并根据相关方法剔除误匹

9、配点，成功重建了东南极Byrd冰川80年代的流速图。该流速图为计算Byrd冰川区域的物质平衡提供数据根底，并为极地气候精密定量化研究提供了有力的方法以及数据支持。【参考文献】【1】张栋，孙波，柯长青，等.南极冰盖物质平衡与海平面变化研究新进展J.极地研究，2021，223：296-305.【2】IPCC.WorkingGroupIContributiontotheIPCCFifthAssessmentReport，ClimateChange2021：ThePhysicalScienceBasis=SummaryforPolieymakersR/OL.2021-10-28. ：/ climate

10、/images/uploads/WFIARS-SPM_Approved27Sep2021.pdf.【3】BuddWF，CorryMJ，JackaTH.ResultsFromTheAmeryIceShelfProjectJ.AnnalsofGlaciology，1982，3：36-41.【4】KingMA，ColemanR，MorganPJ，etal.VelocitychangeoftheAmeryIceShelf，EastAntarctica，duringtheperiod1968-1999J.JournalofGeophysicalResearch，2021，11

11、2F1：F01013.【5】UrbiniS，FrezzottiM，GandolfiS，etal.HistoricalbehaviourofDomeCandTalosDomeEastAntarcticaasinvestigatedbysnowaccumulationandicevelocitymeasurementsJ.Global&PlanetaryChange，2021，603：576-588.【6】AnM，WiensD，AnC，etal.AntarcticicevelocitiesfromGPSlocationsloggedbyseismicstationsJ.AntarcticS

12、cience，2021，272：210-222.【7】LucchittaBK，FergusonHM.Antarctica：measuringglaciervelocityfromsatelliteimagesJ.Science，1986，2344780：1105-8.8RignotE，ScheuchlB.IceflowoftheAntarcticicesheet.J.Science，2021，3336048：1427-30.9BrecherHH.PhotographicdeterminationofsurfacevelocitiesandelevationsonByrdGlacierJ.Antarct.J.US，1982，175：79-81.10RignotE，ThomasRH.MassbalanceofpolaricesheetsJ.Science，2021，2975586：1502-1506.11HeidT，KbA.evaluationofexistingimagematchingmethodsforderivingglaciersurfacedisplacementsgloballyfromopticalsatelliteimageryJ.RemoteSensingofEnvironment，2021，118118：339-355.12Scam