matlab读取hdf的ndvi值，代码应该怎么写，数据来自MOD13A2

MATLAB 中可以使用 `hdf5read` 函数来读取 HDF 格式的数据，如 MODIS 的 HDF 文件（例如 MOD13A2）。NDVI（Normalized Difference Vegetation Index，归一化植被指数）通常是从红光波段和近红外波段的差值计算出来的，你可以从 HDF 文件中提取这两个波段的数据，然后计算 NDVI。

假设你需要从名为 `MOD13A2.hdf` 的文件中读取波段 1 和 2（红色和近红外），然后计算 NDVI，以下是一个简单的示例代码：

% 加载HDF5文件
filename = 'MOD13A2.hdf';
data = h5read(filename, '/Reflectance/Mean_1km');

% 提取红光波段（通常对应Band 1 或者 R1）
band1_idx = 'B1'; % 替换为实际波段编号或名称
red_data = data{find(strcmp(band1_idx, data.VariableName))};

% 提取近红外波段（通常对应Band 2 或者 NIR1）
band2_idx = 'B2'; % 同上
nir_data = data{find(strcmp(band2_idx, data.VariableName))};

% 计算NDVI
red = red_data.Value; % 提取数值数据部分
nir = nir_data.Value;
ndvi = (nir - red) ./ (nir + red); % 除以零可能会有问题，需加上处理异常的条件

% 如果需要保存为新矩阵
ndvi_matrix = reshape(ndvi, [size(data.Data, 1), size(data.Data, 2)]);

% 显示第一个像素的NDVI值
fprintf('First pixel NDVI value: %.2f\n', ndvi_matrix(1,1));

%