Glas数据初步处理

glas数据是点数据，本文采用matlab来读取glas数据并对结果进行存储和展示。

FILE_NAME = 'D:\GLA14-05.10.21\GLAH14_634_2113_001_1275_0_01_0001.H5';
file_id = H5F.open (FILE_NAME, 'H5F_ACC_RDONLY', 'H5P_DEFAULT');
%Open the datasets
LATFIELD_NAME='Data_1HZ/Geolocation/d_lat';
lat_id=H5D.open(file_id, LATFIELD_NAME);
LONFIELD_NAME='Data_1HZ/Geolocation/d_lon';
lon_id=H5D.open(file_id, LONFIELD_NAME);
LEVFIELD_NAME='Data_1HZ/Atmosphere/d_Surface_temp';
temp_id=H5D.open(file_id, LEVFIELD_NAME);
% Read the datasets.
latitude=H5D.read(lat_id,'H5T_NATIVE_DOUBLE', 'H5S_ALL', 'H5S_ALL','H5P_DEFAULT');
lon=H5D.read(lon_id,'H5T_NATIVE_DOUBLE', 'H5S_ALL', 'H5S_ALL','H5P_DEFAULT');
temp=H5D.read(temp_id,'H5T_NATIVE_DOUBLE', 'H5S_ALL', 'H5S_ALL','H5P_DEFAULT');
data=[lon,latitude,temp]; %分别是经度、纬度和地表温度
xlswrite('属性.xlsx',data) %对结果进行输出

打开属性数据，添加字段名称，分别为x，y和temp。然后通过arcgis中add x y数据的方式加载数据，结果如下所示：

image.png

值得注意的是该数据集与传统的栅格数据不一致，该数据只是在特定的点上有，而并不是全部覆盖每个栅格。