方法一、geotiffread（）geotiffwrite（）读写

[A,R] = geotiffread(filename)
从filename指定的GeoTIFF文件中读取一个地理参考的灰度、RGB或多光谱图像或数据网格到A中，并创建一个空间参考对象，R。
[X, cmap, R] = geotiffread(filename)读取一个索引的图像到X和相关的彩色地图到cmap，并创建一个空间引用对象，R。
[A, refmat, bbox] = geotiffread (filename)将地理参考的灰度、RGB或多光谱图像或数据网格读入A，将相应的参考矩阵读入refmat，将边界框读入bbox。
[X, cmap, refmat, bbox] = geotiffread(filename)读取一个索引图像到X，相关的彩色地图到cmap，参考矩阵到refmat，和边界框到bbox。引用矩阵必须由GeoTIFF文件明确地定义，否则它和边界框将返回空。= geotiffread(url)从url读取GeoTIFF图像。= geotiffread(，idx)从多图像GeoTIFF文件或URL中读取一张图像。

基础代码

[Data,R] = geotiffread(); 
[SS,R] = geotiffread(samplename);  %SS为tif数据的信息，是一个二维数组，R是tif数据的地理信息
info = geotiffinfo(samplename);  % 读取tif数据的地理信息，为后面导出为tif数据提供地理信息
[nrows,ncols] = size(SS);  %读取SS的行列数
%导出为tif
outname = ['##.tif'];  %存储位置和名字
geotiffwrite(outname,SS, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);

批量TIFF读取、年均值计算、批量写入（自定义输出路径）

在这里插入代码片

方法二 imfinfo、imread、imshow读取信息、数据，显示图像

读，并显示信息

filepath='D:\study\AOD\AVR\2021002.tif';                                         %%图像名称与路径
Info=imfinfo(filepath);                                      %%获取图片信息并判断是否为tif
 
tif='tif';
format=Info.Format;
if  (strcmp(format ,tif)==0)
    disp('载入的不是tif图像，请确认载入的数据');                %%确保载入的图像是tiff图像
end
 
Slice=size(Info,1);                                          %%获取图片z向帧数
Width=Info.Width;
Height=Info.Height;
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Image=zeros(Height,Width,Slice);
 
for i=1:Slice
    Image(:,:,i)=imread(filepath,i);                         %%一层一层的读入图像
end

读tiff

imageData = read(t) 从与 Tiff 对象 t 关联的 TIFF 文件的当前图像文件目录 (IFD) 中读取图像数据。
[Y,Cb,Cr] = read(t) 从 TIFF 文件的当前图像文件目录中读取 YCbCr 分量数据。仅对具有 YCbCr 光度解释的图像使用此语法。

根据 YCbCrSubSampling 标记的值，Cb 和 Cr 分量的大小可能与 Y 分量不同。

创建一个 Tiff 对象并从 TIFF 文件中读取数据。

t = Tiff('peppers_RGB_tiled.tif','r');
imageData = read(t);

读取 YCbCr TIFF 图像数据

t = Tiff('peppers_YCbCr_tiled.tif','r');
[Y,Cb,Cr] = read(t);

显示图像的 Y 分量。

imshow(Y); 
title('Peppers Image (Y Component)');

写，创建

1、创建tiff图像

MATLAB函数：obj = Tiff(filename, mode)

obj就是创建的Tiff对象
filename是文件名，要带上格式后缀吼，比如filename = ‘img.tiff’
mode是文件访问类型，默认是‘r’读，保存文件的话这里要修改。对于单张图像，改成 ‘w’ 就好啦，如果是多张的话，改成 ‘a’ 。‘r+’ 打开（而非创建）文件以便读写。

2、设置标记值

根据图像的布局，您必须设置以下额外标记：
条状布局 - 设置 RowsPerStrip 标记。
分块布局 - 设置 TileWidth 和 TileHeight 标记。

创建文件后，在写入数据前，要先设置以下标记：
ImageWidth：图像宽度
ImageLength：图像长度
BitsPerSample：比特数/一个像素样本。磕科研中常用的三种：8位，16位，32位。
SamplePerPixel：像素样本数/一个像素
PlanarConfiguration：存储配置。chunky是连续存储每个像素的分量值，seperate是分开存储每个通道。在科研中处理的大多数为单通道灰度图像，这个设置都可。
Photometric：图像数据颜色空间。有好多种可以选。科研处理的灰度图像，一般选’MinIsBlack’，即像素值为0时是黑色的。
SampleFormat：像素样本格式。uint8和uint16时选择’UInt’，32位single时选’IEEEFP’（IEEE浮点数据）
Compression：图像数据压缩方案，科研中选无压缩

t.setTag('ImageLength', size(img,1));
t.setTag('ImageWidth', size(img,2));
t.setTag('Photometric', Tiff.Photometric.MinIsBlack);   %图像数据的颜色空间
t.setTag('BitsPerSample', dtype);                       %数据位数
t.setTag('SamplesPerPixel', 1);
t.setTag('PlanarConfiguration', Tiff.PlanarConfiguration.Chunky);
t.setTag('Compression',Tiff.Compression.None);          %无压缩
t.setTag('SampleFormat',sf);                            %像素样本格式

MATLAB——tiff文件数据读取,modis