MATLAB——tiff文件数据读取,modis
读geotiff
[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中读取一张图像。
读，并显示信息
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); %像素样本格式