ICESat（冰、云和陸地高程衛(wèi)星）是NASA的一項(xiàng)衛(wèi)星任務(wù)，用于測(cè)量冰蓋?質(zhì)量平衡、云和氣溶膠高度，以及陸地地形和植被特征。ICESat 任務(wù)旨在提供確定冰蓋質(zhì)量平衡所需的高程數(shù)據(jù)以及云特性信息，特別是對(duì)于極地地區(qū)常見的平流層云。除了格陵蘭島和南極冰蓋的極地特定覆蓋范圍外，它還提供全球地形和植被數(shù)據(jù)。該衛(wèi)星被發(fā)現(xiàn)可用于評(píng)估重要的森林特征，包括樹木密度【1】。我之前寫過一個(gè)專欄，介紹了如何讀取ICESat數(shù)據(jù)ICESat-2數(shù)據(jù)讀取顯示 - 嗶哩嗶哩 (bilibili.com)。此處所用 ICESat 數(shù)據(jù)是由美國(guó)冰雪數(shù)據(jù)中心（National Snow and Ice DataCenter, NSIDC）發(fā)布的全球陸地表面高程數(shù)據(jù)（Global Land Surface Altimetry Data, GLA14）。

在這里我簡(jiǎn)要減少如何利用ICESat數(shù)據(jù)提取湖泊的水位，這里僅以青海湖為例。

【1】下載數(shù)據(jù)。進(jìn)入數(shù)據(jù)官網(wǎng)：https://openaltimetry.org/data/icesat/。進(jìn)入以下頁(yè)面

【2】選擇青海湖的區(qū)域，并保存為csv文件。

【3】打開下載的文件，前五列分別為經(jīng)度、緯度、高程、時(shí)間和軌道號(hào)。并保存至txt文件，利用matlab處理。

【4】數(shù)據(jù)讀取的代碼：

A = shaperead('qinghai.shp');

lon1 = A.X;

lat1 = A.Y;

plot(SD(1:302,2),SD(1:302,3))

box on

grid on

T = tabulate(SD(:,4));

in=inpolygon(SD(:,1),SD(:,2),lon1,lat1);

output_args=SD(in,:);

figure

hold on

line(lon1,lat1,'marker','.','markersize',1,'color','b');

scatter3(output_args(:,1),output_args(:,2),output_args(:,3),0.1,'r')

set(gca,'gridlinestyle',':','LineWidth',0.5,'GridAlpha',0.5)

box on

grid on

xlim([99.5 101]);

ylim([36.4 37.4]);

AX = tabulate(output_args(:,4))

% AX = output_args(:,4);

ind = find(AX(:,2)~=0);

參考文獻(xiàn)

【1】https://en.wikipedia.org/wiki/ICESat（2022/8/8）.