©2020-2024 ioDraw All rights reserved
<>方法一、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); %像素样本格式