<>方法一、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); %像素样本格式

技术
今日推荐
下载桌面版
GitHub
百度网盘(提取码:draw)
Gitee
云服务器优惠
阿里云优惠券
腾讯云优惠券
华为云优惠券
站点信息
问题反馈
邮箱:[email protected]
QQ群:766591547
关注微信