基于OMEGA影像火星北極冰蓋季節性變化監測

張鼎凱， 劉召芹， 邸凱昌， 岳宗玉， 劉峰， 芶盛

(1.山東科技大學測繪科學與工程學院，青島　266590; 2.中國科學院遙感與數字地球研究所遙感科學國家重點實驗室，北京　100101)

基于OMEGA影像火星北極冰蓋季節性變化監測

張鼎凱1,2， 劉召芹2， 邸凱昌2， 岳宗玉2， 劉峰1， 芶盛2

(1.山東科技大學測繪科學與工程學院，青島266590; 2.中國科學院遙感與數字地球研究所遙感科學國家重點實驗室，北京100101)

摘要：火星兩極冰蓋每年季節性消融的空間范圍變化是火星全球及區域氣候變化最直觀的反映，針對火星冰蓋監測提出了一種利用歐空局火星快車搭載的OMEGA高光譜成像數據提取季節性冰蓋消融線的方法。基于28至29火星年覆蓋火星北半球6個時段的OMEGA紅外波段影像，利用監督分類方法辨別冰層與裸地,以提取季節性冰蓋的邊界，計算了火星北半球季節性冰蓋消融速度，并利用高分辨率HiRISE 數據驗證了本文研究方法的正確性。研究結果表明，火星北極季節性冰蓋的消融速度是每隔太陽經度(LS)10°冰層消融106 km2； 通過對比MOLA激光高度計數據和地形數據，發現局部區域冰蓋消融異常原因主要為撞擊坑引起的地形變化所致。

關鍵詞：火星北半球； OMEGA影像； HiRISE影像； 季節性冰蓋； 消融曲線

0引言

在火星高緯度地區，存在永久性冰蓋和季節性冰蓋，其中永久性冰蓋全年都存在，季節性冰蓋隨火星季節的變化每年周期性地生成和消失。關于火星冰蓋的探測，國內外學者利用現有的遙感影像數據對季節性冰蓋的變化規律進行了研究，發現季節性冰蓋主要是由大氣中的CO2和水汽凝結而成[1-2]。火星北極的永久性冰蓋主要是由水冰組成[3]，在水冰上面覆蓋的干冰厚度大約為1m左右，在夏季時完全消融； 南極永久性冰蓋由干冰和水冰覆蓋而成[1,4]，干冰冰層厚度在8m左右，夏季依然存在。由于季節性冰蓋的主要組成成分是干冰和水冰，隨著火星北半球季節性冰蓋的消融，大量CO2氣體進入大氣中，北半球氣壓升高，與此同時火星南半球進入冬季，CO2凝結，南北半球形成氣壓差，形成火星大氣循環，因此季節性冰層的揮發與凝結作為火星大氣循環的動力之一發揮著重要作用[5]。火星兩極冰蓋每年季節性隨時間和空間的變化是火星全球及區域氣候變化最直觀的反映[6]，通過對兩極冰蓋消融線的自動提取和制圖，計算冰蓋消融的面積和范圍，結合利用火星軌道激光高度計(Marsorbiterlaseraltimeter，MOLA)地形變化測量[7]，可以有效地估算進入大氣的CO2數量，為研究火星氣候變化提供基礎數據。

對火星季節性冰蓋的觀測方式主要分為地基觀測、太空觀測和遙感衛星觀測3個階段： 早期Antoniadi[8]和Fischbacher等[9]采用地基觀測手段在地球上進行了許多對火星極冠的觀測； 太空觀測階段Cantor等[10]利用哈勃太空望遠鏡對火星北極季節性冰蓋消融進行過數次觀測，進一步了解了火星季節性冰蓋的變化規律；Capen等[11]和Iwasaki等[12-14]由地基觀測得到的數據發現在火星北極冰蓋的消融過程中有一個停滯期； 對火星北極進行最早的遙感衛星觀測是Benson和James等[15-17]通過利用水手9號和海盜號進行的，發現大約數周時間冰蓋的邊界停在65°N；James等[18]還利用火星全球探勘者(Marsglobalsurveyor，MGS)攜帶的火星軌道相機(Marsorbitercamera，MOC)對火星北極冰蓋進行了周期性的觀測，從影像上利用人工識別的方法得到了季節性冰蓋的消融規律，并且提出了火星北半球太陽經度(solarlongitude，LS)與季節性冰蓋消融線緯度位置的函數，概略地闡述了火星塵暴對冰蓋消融的影響。

目前對火星冰蓋季節性變化的研究主要基于火星快車獲取的遙感數據。Giuranna等[19]利用火星快車的行星傅里葉分光計的短波通道研究了火星北半球春季冰蓋的組成成分和厚度，分析了不同區域的干冰粒徑大小和水冰豐富度，利用反照率模型推算出不同緯度區域的干冰層的單位質量，進而推算出干冰的厚度；Appéré等[20]利用火星快車的光學與紅外礦物光譜儀(ObservatoirepourlaMinéralogie,l’Eau,lesGlacesetl’Activité，OMEGA)影像對火星北半球自27火星年(Martianyear,MY)的冬季到28火星年的春季進行了一系列觀測(火星年由Clancy[21]定義，地球公元1955年春分作為火星第1年的起始)，分別利用1.429μm和1.500μm波段對干冰和水冰反照率的閾值進行分析，得到了2種冰層的消融界線隨時間的變化規律，做出了季節性冰層的變化分布曲線圖。但是，由于不同時相影像的光照條件、大氣密度和傳感器高度不統一等因素，很難得到統一的大范圍、多時相的冰蓋反照率的閾值邊界，進而依據單波段閾值來確定冰蓋的消融邊界具有一定的不足。針對這些不足，本文提出利用北半球OMEGA影像，選取其中的128個近紅外波段，采用地物監督分類方法，對實驗區域的每幅影像選取感興趣區進行監督分類，這樣就避免了因多幅影像使用同一閾值提取冰蓋邊界造成的不準確性。通過監督分類方法提取冰蓋的消融曲線，分析季節性冰蓋的消融趨勢，擬合出冰蓋消融線隨時間變化的曲線，得到火星北極冰蓋變化的時空規律； 另外，由于火星表面的季節性沉積主要是水冰、干冰以及灰塵，春季升華時會在不同的空間尺度引起不尋常的、動態的消融現象，如深色斑點、扇形沉積、蛛網狀地形、裂縫、邊坡條紋、沙雪崩裂和明暗條紋等，在可見光范圍內可以直觀地觀察到火星季節性冰蓋變化引起的一些升華活動[15,18]。通過火星偵察軌道器攜帶的0.25m高分辨率成像科學實驗相機(highresolutionimagingscienceexperiment,HiRISE)能夠直接觀測到北半球季節性冰蓋的消融過程，驗證使用OMEGA影像提取冰蓋消融曲線的正確性。

1數據源及其預處理

1.1數據源

火星快車發射于2003年6月，12月進入火星軌道[22]，OMEGA是火星快車的一個重要載荷，主要用于采集可見光和近紅外數據，共352個波段。在0.35～1.0μm，0.93～2.65μm和2.51～5.1μm的光譜范圍采樣間隔分別是7nm，13nm和20nm，具體參數見表 1[23]。

依據干冰和水冰的波譜反射特征選取1～128波段用于冰蓋消融邊界的提取。為了檢驗分類結果的正確性，使用2005年8月HiRISE相機獲取的高分辨率影像進行驗證。HiRISE影像的空間分辨率達0.25m，其3個波段的波長范圍為： 藍—綠光波段(400～600nm)、紅光波段(550～850nm)與近紅外波段(800～1 000nm)，從影像上可以直接觀測到不同時期的火星冰蓋消融變化。

1.2數據預處理

從美國國家航空航天局(NASA)PDS網站提供的數據中選擇火星北半球LS在340°～110°范圍的原始數據。利用IDL對數據進行讀取，提取需要的信息，包括灰度值(idat)、輻射亮度值(jdat)、火星距日心處的太陽光譜值(specmars)和經緯度(geocube)。將數據信息導入ENVI中，得到需要的數據類型，選取需要的波段存儲為影像格式。

由于火星大氣組成比較單一，對太陽輻射產生影響的主要是CO2，因此使用OMEGA團隊提出的“經驗傳輸函數法”(ETF算法)進行大氣校正，其原理是利用火星上的奧林帕斯火山頂部和底部相近地物的輻射率(I/F)差異，對頂部和底部的光譜做比值，得到火星大氣輻射傳輸因子(S)，又稱為大氣光譜。然后將觀測像元光譜除以S，基本消除大氣對輻射傳輸的影響[24]。

1.3光譜曲線特征分析

火星北極季節性冰蓋的組成成分是干冰和水冰的混合物，Langevin等[25]在季節性冰蓋覆蓋的區域得到了包含水冰、干冰以及干冰水冰混合物的光譜曲線。對實驗區影像預處理后進行反照率計算，使用OMEGA數據文件中的jdat與specmars做比值運算，然后選取影像1～128波段范圍數據，得到地表的光譜曲線(圖 1)。圖1中實線所示光譜為LS=31°時位于(55.37°N，14.76°E)區域的裸地光譜曲線； 虛線所表示的是LS=31°時位于(79.45°N，2.16°E)區域的冰蓋光譜曲線。在LS=31°時冰蓋的消融線大約位于65°N左右的位置[15]，因此選擇上述2條曲線作為選取感興趣區的參考光譜曲線是可靠的，以此選取感興趣區對實驗區域進行監督分類，提取季節性冰蓋的消融曲線。

(OMEGA影像序列號： ORB5285)

利用不同地物的光譜特征，采用監督分類的馬氏距離法把地物分為裸地和冰蓋2大類。為了驗證分類方法的正確性，選擇南半球LS=195°E時相的影像進行分類并與其他分類結果對比(圖 2)。將本文分類后結果與Schmidt分類結果[6]疊加對比發現，兩者的冰蓋消融分界線位置基本一致，都位于60°S附近。Schmidt的實驗結果是依據由不同物質的波段深度得到的經驗公式進行閾值計算后得到的。由此說明了本文提取季節性冰蓋消融曲線的方法是可行的，且方法更為簡單直觀。

(a)本文實驗結果 (b)Schmidt分類結果

圖2南極區域消融對比

Fig.2Recessioncomparisonofsouthpolarcap

2實驗結果與驗證

2.1實驗結果

選取在MY28年春季和MY29年春季火星北半球的部分OMEGA影像為實驗數據，共95景影像。依據影像覆蓋范圍分為6個時間段，采用馬氏距離法分別對OMEGA影像進行監督分類，結果如圖 3所示。

(a)LS范圍： 10°～14°(b) LS范圍： 20°～24° (c) LS范圍： 30°～34°

(d) LS范圍： 67°～72°(e) LS范圍： 80°～84° (f) LS范圍： 90°～94°

圖3北極地區冰蓋提取結果

Fig.3Extractionresultsofnorthpolarcap

圖3中綠色部分表示冰蓋分布區域，紅色部分表示裸地區域，以紅綠之間的分界線(白色實線)提取出不同時間段的冰蓋消融邊界。

從圖 3中冰蓋消融變化可以看出，在(195°E，71°N)處的撞擊坑附近，出現消融異常的現象。因此，對撞擊坑附近的消融異常區進行重點分析。圖4(a)綠色方框為撞擊坑的位置； 圖4(b)為使用MGS的MOLA數據生成的暈渲圖，可以觀察到撞擊坑周圍的地形比較粗糙； 圖4(c)和(d)展示了撞擊坑在2個時期的冰蓋消融變化現象，發現撞擊坑附近消融速度明顯慢于同緯度的平坦區域。通過對其他撞擊坑的統計發現，撞擊坑上的冰蓋消融速度普遍比相同緯度的平坦地區緩慢，因此推斷可能是在撞擊坑附近，由于坑沿地形粗糙度較大，產生的陰影致使太陽輻射減少，使冰蓋的消融速度降低，形成消融異常區。

(a)撞擊坑位置(b) 暈渲圖(c) LS=64.9°(d) LS=66.2°

圖4消融異常區域

Fig.4Abnormalrecessionalarea

通過對火星北半球得到的分類結果分析發現，在春季開始時冰蓋消融曲線基本沿著同一緯線圈，不同區域的地形對冰蓋消融的影響不大。隨著時間推移，消融曲線呈現出不規則的形狀，隨著緯度的升高，太陽高度角變小，地形起伏產生的陰影成為影響冰蓋消融的主要因素，因此推斷隨時間的推移地形對冰蓋消融起到了主要作用。對每個時間段的冰蓋面積進行統計，得到的結果疊加到火星北極極區底圖上，如圖 5 所示。不同時間段的冰蓋面積統計見表 2。對冰蓋面積和太陽經度進行曲線擬合，得到圖 6所示擬合曲線，冰蓋面積與太陽經度近似呈線性關系。

圖5　冰蓋面積統計

LS/(°)面積/km212848506522743864932644994770273015482100901892851446

圖6　冰蓋面積與太陽經度的關系

從圖 6中可以得出北半球季節性冰蓋的消融速度大約是106km2/10°。當太陽經度到達90°時，冰蓋面積保持穩定，由圖 5也可以觀察到冰蓋消融到80°N時基本不再消融，剩余的冰蓋到達夏季仍未消融，作為永久性冰蓋存在。

2.2基于HiRISE影像的冰蓋變化驗證

從HiRISE高分辨率相機影像上直觀得到的冰蓋隨時間的變化規律，可以驗證利用高光譜影像進行地物分類的正確性。(70°N，103°E)區域不同時間序列的HiRISE紅光波段影像如圖 7所示。

圖7　位于(70°N，103°E)處不同時間序列的HiRISE紅光波段影像

從圖7中可以看到LS在37°～48°時開始有裸露土地(紋狀線間的黑色斑塊)出現，白色物質與周圍地物對比明顯，隨著時間推移，白色物質越來越少，因此推斷白色物質即為冰蓋。在LS=55°時，冰蓋基本消失，只是在地表縫隙處存在較明顯的冰蓋，可能是在地表縫隙中由于太陽直射較少影響冰蓋揮發，在圖 3(c)中LS=30°時可以發現北緯70°附近為冰蓋，在圖 3(d)中LS=70°時北緯70°附近已經沒有冰蓋，與圖 7中冰蓋的變化一致。在圖 7所示區域LS=55°時只是在表面的大縫隙中存在少量冰蓋，但在圖 8和圖9不同地區的假彩色影像中LS=78°時在大的縫隙中仍存有冰蓋，而相同時間圖 7所示區域已經看不到冰蓋的存在。由此可以得出： 緯度越高，冰蓋消融的越慢； 同時發現沙丘坡度較緩的一側冰蓋消融速度要快于坡度較陡的一側； 冰蓋的消融與太陽光照有關。

圖8　位于(75°N，94°E)處不同時間序列的HiRISE假彩色影像

圖9　位于(76°N , 89°E)處不同時間序列的HiRISE假彩色影像

2.3對比分析

Cantor等[10]和James[26]利用MOC影像在可見光范圍描繪出的季節性冰蓋的變化規律圖像，與本文中的實驗結果對比，在LS=20°和LS=70°時對應的消融邊界基本吻合，同時發現在經度范圍E210°～270°之間的消融速度明顯快于E315°～30°范圍的冰蓋消融速度。這是由于在E210°～270°區域的海拔高于E315°～30°[27]，海拔越高氣壓越低，冰蓋的升華點隨壓強的減小而逐漸降低，即相同光照條件下，氣壓越低冰蓋越容易揮發[28]。

Appéré等[20]通過OMEGA影像的1.429μm和1.500μm波段得到的深度信息研究冰蓋的消融邊界。圖 10展示了利用MOC，OMEGA和MARCI影像得到的冰蓋消融邊界變化。通過圖 10中多元數據實驗結果對比分析，本文得到的北半球季節性冰蓋的邊界與使用其它數據和方法得到的結果基本一致，驗證了本文方法的可靠性。

(黑色擬合曲線為本文的實驗結果)

3結論

依據火星季節性冰蓋與裸地的不同光譜特征，選取火星快車的OMEGA高光譜影像，利用監督分類的方法提取季節性冰蓋的消融曲線，探討了冰蓋邊界線隨時間的變化關系，以此推斷出冰蓋消融速度為每隔太陽經度10°冰蓋面積變化106km2。由HiRISE高分辨率影像直觀地展示了火星北半球季節性冰蓋在春季消融的現象，對不同時刻的影像所展示的冰蓋消融現象與擬合的冰蓋消融曲線做對比分析，驗證了冰蓋消融曲線的正確性。通過研究季節性冰蓋消融規律，能夠更好地了解火星冰蓋變化的規律，為進一步探究火星大氣變化以及全球氣候變化對比研究奠定了基礎。如果結合激光高度計數據測量的冰蓋變化量，可以進一步更精確地計算消融的體積及質量，為火星氣候環境模擬提供數據支撐。

參考文獻(References)：

[1]BibringJP,LangevinY,PouletF,etal.PerennialwatericeidentifiedinthesouthpolarcapofMars[J].Nature,2004,428(6983):627-630.

[2]LangevinY,PouletF,BibringJP,etal.SummerevolutionofthenorthpolarcapofMarsasobservedbyOMEGA/Marsexpress[J].Science,2005,307(5715):1581-1584.

[3]SchmittB,SchmidtF,DoutéS,etal.RecessionofthenorthernseasonalcondensatesonMarsbyOMEGA/Marsexpress[J].LPIContributions,2006,1323:8050.

[4]PaigeDA,HerkenhoffKE,MurrayBC.Mariner9observationsofthesouthpolarcapofMars:EvidenceforresidualCO2frost[J].JournalofGeophysicalResearch:SolidEarth,1990,95(B2):1319-1335.

[5]FeldmanWC,PrettymanTH,BoyntonWV,etal.CO2frostcapthicknessonMarsduringnorthernwinterandspring[J].JournalofGeophysicalResearch:Planets,2003,108(E9):5103.

[6]SchmidtF.DoutéS,SchmittB,etal.AutomaticdetectionofH2OandCO2icesinOMEGA/MEXimagesforthemonitoringofthesouthpolarcaprecession[C]//SchmidtF,DoutéS,SchmittB,etal.Proceedingsofthe37thannuallunarandplanetaryscienceconference.League,Texas:LPI,2006:1979.

[7]SmithDE,ZuberMT,NeumannGA.SeasonalvariationsofsnowdepthonMars[J].Science,2001,294(5549):2141-2146.

[8]AntoniadiEM.LaPlaneteMercureEtLaRotationDesSatellites;EtudeBaseeSurLesResultatsObtenusAvecLaGrandeLunetteDeL’observatoiredeMeudon[M].Paris:Gauthier-Villars,1934.

[9]FischbacherGE,MartinLJ,BaumWA,etal.Seasonalbehaviorofthemartianpolarcaps[J].PublicationsoftheAstronomicalSocietyofthePacific,1969,81:538.

[10]CantorBA,WolffMJ,JamesPB,etal.Recessionofmartiannorthpolarcap:1990—1997hubblespacetelescopeobservations[J].BulletinoftheAmericanAstronomicalSociety,1997,29:963.

[11]CapenCF,CapenVW.Martiannorthpolarcap,1962—1968[J].Icarus,1970,13(1):100-108.

[12]IwasakiK,SaitoY,AkabaneT.Behaviorofthemartiannorthpolarcap,1975—1978[J].JournalofGeophysicalResearch:SolidEarth(1978-2012),1979,84(B14):8311-8316.

[13]IwasakiK,SaitoY,AkabaneT.Martiannorthpolarcap1979—1980[J].JournalofGeophysicalResearch:SolidEarth(1978—2012),1982,87(B12):10265-10269.

[14]IwasakiK,ParkerDC,LarsonS,etal.Martiannorthpolarcap1996—1997[J].Icarus,1999,138(1):20-24.

[15]BensonJL,JamesPB.Yearlycomparisonsofthemartianpolarcaps:1999—2003Marsorbitercameraobservations[J].Icarus,2005,174(2):513-523.

[16]JamesPB.Recessionofmartiannorthpolarcap:1977—1978vikingobservations[J].JournalofGeophysicalResearch:SolidEarth(1978—2012),1979,84(B14):8332-8334.

[17]JamesPB,NorthGR.TheseasonalCO2cycleonMars:Anapplicationofanenergybalanceclimatemodel[J].JournalofGeophysicalResearch,1982,87(B):10271-10283.

[18]JamesPB,CantorBA.Martiannorthpolarcaprecession:2000Marsorbitercameraobservations[J].Icarus,2001,154(1):131-144.

[19]GiurannaM,HansenG,FormisanoV,etal.Spatialvariability,compositionandthicknessoftheseasonalnorthpolarcapofMarsinmid-spring[J].PlanetaryandSpaceScience,2007,55(10):1328-1345.

[20]AppéréT,SchmittB,LangevinY,etal.WinterandspringevolutionofnorthernseasonaldepositsonMarsfromOMEGAonMarsexpress[J].JournalofGeophysicalResearch,2011,116(E5):E05001.

[21]ClancyRT,SandorBJ,WolffMJ,etal.Anintercomparisonofground-basedmillimeter,mgstes,andvikingatmospherictemperaturemeasurements:SeasonalandinterannualvariabilityoftemperaturesanddustloadingintheglobalMarsatmosphere[J].JournalofGeophysicalResearch:Planets(1991—2012),2000,105(E4):9553-9571.

[22]ChicarroA,MartinP,TrautnerR.TheMarsExpressMission:AnOverview[M]//WilsonA.MarsExpress:TheScientificPayload.Noordwijk,Netherlands:ESAPublicationsDivision,2004,1240:3-13.

[23]BibringJP,SoufflotA,BerthéM,etal.OMEGA:ObservatoirepourlaMinéralogie,l’Eau,lesGlacesetl’Activité[M]//WilsonA.MarsExpress:TheScientificPayload.Noordwijk,Netherlands:ESAPublicationsDivision,2004,1240:37-49.

[24]胡全一,祝民強.火星高光譜遙感大氣校正方法[J].東華理工大學學報:自然科學版,2008,31(3):242-248.

HuQY,ZhuMQ.MethodsofatmosphericcorrectionforhyperspectralremotesensingonMars[J].JournalofEastChinaInstituteofTechnology:NaturalScienceEdition,2008,31(3):242-248.

[25]LangevinY,VincendonM,BibringJP,etal.Multi-bandprocessingofobservationsofthepolarregionsofMarswiththeOmegaimagingspectrometeronboardMarsexpress[C]//Proceedingsofthe2ndHyperspectralImageandSignalProcessing:EvolutioninRemoteSensing.Reykjavik:IEEE,2010:1-4.

[26]JamesPB,ClancyRT,LeeSW,etal.MonitoringMarswiththehubblespacetelescope:1990—1991observations[J].Icarus,1994,109(1):79-101.

[27]ZuberMT,SmithDE,SolomonSC,etal.ObservationsofthenorthpolarregionofMarsfromtheMarsorbiterlaseraltimeter[J].Science,1998,282(5396):2053-2060.

[28]FrayN,SchmittB.Sublimationoficesofastrophysicalinterest:Abibliographicreview[J].PlanetaryandSpaceScience,2009,57(14/15):2053-2080.

MonitoringofseasonalchangesofMartiannorthpolar

(責任編輯： 陳理)

icecapwithOMEGAimages

ZhangDingkai1,2,LIUZhaoqin2,DIKaichang2,YUEZongyu2,LIUFeng1,GOUSheng2

(1. Survey and Mapping Institute of Science and Technology， Shandong University of Science and Technology, Qingdao 266590,China; 2. State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China)

Abstract:The spatial variation extent of the annual seasonal melt of Martian polar ice caps is the most direct reflection of Mars global and regional climate changes. A method that utilizes hyperspectral images from OMEGA spectrometer on board ESA Mars Express for monitoring Mars ice cap by extracting seasonal ice cap ablation line is proposed in this study. Based on OMEGA infrared images from 6 periods of Martian year 28 and 29 that cover Martian northern hemisphere, the boundary of seasonal ice cap was extracted by supervised classification method which can distinguish between ice and bare land, and the melting rate of seasonal ice cap was also calculated and analyzed, with validation from high resolution HiRISE images. The results show that the melting rate of Martian northern polar ice cap is about 106 km2 every 10° solar longitude (LS). In addition, the comparison between the results and the terrain data from MOLA reveals that the regional abnormality of ice cap melting is mainly caused by the crater-induced topographic variation.

Keywords:Northern hemisphere of Mars; OMEGA images; HiRISE images; seasonal ice cap; melting curve

doi:10.6046/gtzyyg.2016.02.16

收稿日期：2014-11-26；

修訂日期：2015-03-16

基金項目：國家自然科學基金項目“火星壁壘撞擊坑成因機制分析及數值模擬”(編號： 41472303)資助。

中圖法分類號：TP 79

文獻標志碼：A

文章編號：1001-070X(2016)02-0099-07

第一作者簡介：張鼎凱(1987-),男,碩士生,主要從事行星遙感制圖方面的研究。Email： zhangdk123@sina.com。

通信作者：劉召芹(1973-),男,博士,副研究員,主要從事遙感制圖與探測車導航定位方面的研究。Email： liuzq@radi.ac.cn。

引用格式： 張鼎凱,劉召芹,邸凱昌,等.基于OMEGA影像火星北極冰蓋季節性變化監測[J].國土資源遙感,2016,28(2):99-105.(ZhangDK,LiuZQ,DiKC,etal.MonitoringofseasonalchangesofMartiannorthpolaricecapwithOMEGAimages[J].RemoteSensingforLandandResources,2016,28(2):99-105.)