matlab使用geotiffread函数、geotiffinfo函数和geotiffwrite函数输出多波段.tif格式文件中的每个波段

您可以使用以下代码来读取、获取信息和写入多波段 GeoTIFF 文件的每个波段：

% 读取多波段 GeoTIFF 文件
filename = 'multiband.tif';
[A, cmap, R] = geotiffread(filename);

% 获取 GeoTIFF 文件信息
info = geotiffinfo(filename);

% 获取每个波段的数据和元数据
for i = 1:size(A, 3)
    % 获取当前波段的数据
    data = A(:,:,i);
    
    % 获取当前波段的元数据
    metadata = info.GeoTIFFTags.GeoKeyDirectoryTag.GeoDoubleParams{i};
    
    % 写入当前波段为单波段 GeoTIFF 文件
    output_filename = sprintf('band%d.tif', i);
    geotiffwrite(output_filename, data, R, 'GeoKeyDirectoryTag', metadata);
end

其中，`geotiffread` 函数用于读取多波段 GeoTIFF 文件，返回一个包含所有波段数据的三维数组 `A`，一个色彩映射表 `cmap` 和一个参考对象 `R`。`geotiffinfo` 函数用于获取 GeoTIFF 文件的元数据信息，返回一个结构体 `info`。在循环中，我们使用 `i` 迭代每个波段，从 `A` 数组中提取当前波段的数据，并从 `info` 结构体中获取当前波段的元数据。最后，我们使用 `geotiffwrite` 函数将当前波段的数据和元数据写入单波段 GeoTIFF 文件中。文件名为 `band{i}.tif`，其中 `{i}` 为当前波段的索引。