matlab对tif文件的相关操作

前言：
本片博客记录了使用matlab读取tif文件的像素数组、标签信息，以及写tif文件。

问题前瞻：

1. 无法对已有的tif文件重新读取和写操作
2. 写操作成功，但显示跟原tif文件不一样，有时出现全黑的图片

代码如下：

clear all;

% 读取tif文件标签信息，为下面tag设置做参考
% Info = imfinfo('1.tif');
% 或者
t0 = Tiff('1.tif','r');

tags = ["ImageLength","ImageWidth","SampleFormat","BitsPerSample",...
    "SamplesPerPixel","Photometric","RowsPerStrip","PlanarConfiguration"];
for i = 1:length(tags)
    tag = getTag(t0, tags(i));
    tagValue(1,i) = tags(i);
    tagValue(2,i) = tag;
end

%t0Img = read(t0);   % bitspersample 为12位！！！！

% 待存储的图像数组，可以为任意二维数组，就以 1.tif 为例
data = imread('1.tif');
min1 = min(min(data));
max1 = max(max(data));
%data = data * 16;
data = uint16(data);
min2 = min(min(data));
max2 = max(max(data));
% imshow(data,[])	% 加[]使适应灰度值的显示

% 生成一个tif对象，为写入数据做准备
mytif = 'demo.tif';
t = Tiff(mytif,'w');

% 设置tif标签信息

tagstruct.ImageLength     = size(data,1);   % 长度
tagstruct.ImageWidth      = size(data,2);   % 宽度
tagstruct.SampleFormat    = 1;   %   对数据类型的解释(1-4)
tagstruct.BitsPerSample   = 16;     % 像素数值位数，single为单精度，32位系统为32.!!!!=12
tagstruct.SamplesPerPixel = 1;      % 像素波段个数，一般为1或3（rgb），遥感图像常大于3
tagstruct.Photometric     = 1;    %   颜色空间解释方式:最小值表示黑色
tagstruct.RowsPerStrip    = 12;     %
tagstruct.PlanarConfiguration = 1;  %
tagstruct.Software        = 'MATLAB';   %生成图片的软件

% 将上述tag信息写入
t.setTag(tagstruct)

% 开始写数据
t.write(data);
% 关闭图片
t.close;

% 检验矩阵是否相同，显示是否正常
img = imread('demo.tif');
imshow(img)
isequal(img,data)

问题解决：
1. 无法对tif图像进行读操作
tif源文件BitsPerSample属性为12位，无法使用matlab中的Tiff直接读取！这也是下面显示图片为全黑的一个原因。

2. 读写成功后的文件，用电脑查看显示为几乎全黑
以上述为例，待存储的图像灰度值数组范围为 0 ~ 4095，这个数字正好是上述12位能够存储的数据范围，即0 ~ 2^12 -1。但上面写tif文件时，设置BitsPerSample为16位（读写tif时都不能为12位，所以用imread读数据，用16位写tif），也就是说能够存储的灰度数据范围为0~2^16-1。

显示为几乎全黑的原因就来了，保存成功后，当用电脑自带的图片浏览器查看图片时，浏览器检测到BitsPerSample为16位，所以就以 0~ 2^16-1 的灰阶显示图片，但实际上，图片的灰度值范围只为 0~ 2^12-1。正常情况下，如果以0~ 2^12 -1灰阶范围来显示图片，最暗为0，最亮为2^12-1，此时为正常，但如果以0~ 2^16-1灰阶来显示数据量为2 ^12 -1的灰度值图像时，同样最暗为0，但此时2^12-1灰度值的像素此时就不是最亮的了，而是最亮的1/16左右，因此会显示为灰色，最大灰度值显示尚且如此，那么整张图可不就显示为几乎全黑了吗。

如果把原始的待存储的图像灰度值数组data*16，也就是使数组灰度值范围等比扩大到 致0~2^16-1范围，最后就会发现，写tif后显示的图片跟原图片几乎一样。

BTW，imshow()函数中，加入[]的意思就是以像素灰度值 min~max的范围显示图片。