利用MATLAB批量读写tif数据（Landsat7&8）

thegirl@chx

已于 2023-02-28 18:08:09 修改

阅读量1.3k

点赞数 3

分类专栏： MATLAB 文章标签： matlab 开发语言图像处理

于 2023-02-28 18:04:36 首次发布

本文链接：https://blog.youkuaiyun.com/weixin_43347033/article/details/129266813

版权

MATLAB 专栏收录该内容

2 篇文章

订阅专栏

文章讲述了使用MATLAB进行Landsat7和Landsat8的tif图像批量读取与处理。在读取Landsat7数据时因行列号差异导致导出错误，而Landsat8数据则能正常处理。主要步骤包括读取tif文件，计算NDVI，以及根据特定坐标系统导出图像。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

利用MATLAB批量读写tif数据（Landsat7&8）

代码来源于：批量读写
代码如下：

clc;clear
%% 读取NDVI-tif文件
file_path_N = 'E:\';% 图像文件夹路径 
img_path_list_N = dir(strcat(file_path_N,'*.tif'));%获取该文件夹中所有TIF格式的图像 
img_num_N = length(img_path_list_N);%获取图像总数量
II_N=cell(1,img_num_N);

if img_num_N > 0
    for jj = 1:img_num_N
        image_name_N = img_path_list_N(jj).name;% 图像名
        [image_N,geo] = geotiffread(strcat(file_path_N,image_name_N));
        II_N{jj}=image_N;
        fprintf(' %d %s\n',jj,strcat(file_path_N,image_name_N));% 显示正在处理的图像名
    end
end
%% 计算植被覆盖度
NDVImax = 0.8;
NDVImin = 0.2;
fc_1 = cell(1,img_num_L);
for i = 1:img_num_L
    fc_1{i} = ((II_N{i}-NDVImin)./(NDVImax-NDVImin)).^2;
end
 %% 输出TIF图像
 for j = 1:1:img_num_L
     CoordRefSysCode=32648; 
     outpath='E:\';
     outFile=strcat(outpath,img_path_list_L(j).name);
     geotiffwrite(outFile,fc_1{i},geo,'CoordRefSysCode',CoordRefSysCode);
 end