2

ImageJ HyperStacks 具有数据类型(8-32 位)、宽度、高度、通道数、切片数和帧数。它们代表 5D (XYCZT) 数据集。ImageJ 将它们存储为多页 tiff 文件,其中存储了 2D (XY) 图像的通道数乘以切片数乘以帧数。第一张图片似乎有两个自定义标签,ID 为 50838 和 50839。

我想创建包含来自 Matlab 的 5D 数据的 tif 文件,ImageJ 可以将其作为有效的 5D HyperStack 读取。

我可以使用 Matlab 将多个 2D 图像存储在多页 tiff 文件中,imwrite(matrix, file, 'WriteMode','append')但 ImageJ 只会将其读取为 3D (XYZ) 图像堆栈。不包含有关通道、切片和帧的信息。

我想我可以查看 ImageJ 的源代码以找出它们存储这些缺失信息的位置,然后使用 Matlab 的 LibTIFF 包装器来重新创建 ImageJ 的元信息。但是,如果您已经知道该怎么做,或者是否有其他选择,我想听听。

4

3 回答 3

3

我还认为 Bio-Fomats 是唯一的方法。但后来我意识到斐济的 ImageJ-MATLAB 允许您将 MATLAB 数据传递给 ImageJ。然后,您可以通过 ImageJ 保存图像,以便 ImageJ 可以肯定地打开它。

问题是人们(或者只是一个人?)创建了 ImageJ-MATLAB 在创建从 MATLAB 到 ImageJ 的路线方面做得很好,但他们似乎并没有费心优化行为。ImageJ-MATLABIJM.show('name')有很多隐藏或未记录的限制(我大量更改了文档,所以现在应该更清楚了)。最终,我编写了一个包装函数

ijmshow

https://github.com/kouichi-c-nakamura/ijmshow

示例 MATLAB 代码

你可以用 4 行来完成。

addpath '/Applications/Fiji.app/scripts' % depends your Fiji installation
ImageJ

imp = ijmshow(I,'YXCZT') % I is 5D array of uint8 or uint16
% the second argument determines which dimension represents what
% imp is a 5D hyperstack

ij.IJ.saveAsTiff(imp, 'image1.tif');

您可能需要在保存前设置每个通道的显示范围。

添加于 2018 年 5 月 7 日

copytoImagePlus提供比ijmshow上面更好的解决方案(您可以使用相同的语法),因为copytoImagePlus不再依赖IJM基础工作区中的变量。

https://github.com/kouichi-c-nakamura/copytoImagePlus

addpath '/Applications/Fiji.app/scripts' % depends your Fiji installation
ImageJ

imp = copytoImagePlus(I,'YXCZT') % I is 5D array of uint8 or uint16
% the second argument determines which dimension represents what
% imp is a 5D hyperstack

ij.IJ.saveAsTiff(imp, 'image1.tif');

另请参阅copytoImgPlus哪个创建 ImageJ2ImgPlus对象。

https://github.com/fiji/fiji/blob/master/scripts/copytoImgPlus.m

于 2018-05-03T20:18:07.073 回答
2

ImageJ 使用 TIFF 的第一个 IFD 的 ImageDescription 标签来存储它的元信息。以下是 Mitosis 样本数据集的一个示例(文件 > 打开样本 > 有丝分裂):

ImageJ=1.51d
images=510
channels=2
slices=5
frames=51
hyperstack=true
mode=composite
unit=\u00B5m
finterval=0.14285714285714285
loop=false
min=1582.0
max=6440.0

MATLAB 的Tiff类可能为您提供足够细粒度的控制,以便将自定义 ImageDescription 写入第一个 IFD;我不知道。

可能更容易使用Bio-Formats库的bfsave.m函数来编写 OME-TIFF,它可以是 5 维的;ImageJ 可以通过Bio-Formats 插件读取这种格式。我推荐使用 ImageJ 的Fiji发行版,它预装了 Bio-Formats。

MATLAB Central 上也发布了同样的问题

于 2017-11-26T19:06:10.500 回答
1

使用 Matlab 的Tiff 类

我可以组装一个我认为看起来不错的 ImageDescription 标签,ImageJ 将它(File>Open)作为 5D Hyperstack 读取,但数据已部分损坏。似乎存在偏移问题(我可以用 Matlab 再次很好地读取数据)。

但是,如果使用File>Import>Bio-Formats读取 tiff 文件,如果在 Bio-Formats 导入对话框中标记了带有 Hyperstack 的视图堆栈,它将在 ImageJ 中变为 5D (XYCZT) 堆栈。感谢ctrueden的有益评论。

d = ones(100, 200, 10, 2, 3, 'single') * 3.765;
hyperstack_write('test.tif', d);

function hyperstack_write(file, hyperstack)
% Writes an (up to) 5D stack into a (XYCZT) HyperStack Tiff file
% readable by ImageJ
%
% hyperstack must have class single

% simple checks
assert(nargin == 2, 'Not enough arguments');
assert(isa(hyperstack, 'single'), 'hyperstack must be single');

% get all five dimensions
d = zeros(5, 1);
for i = 1 : 5
    d(i) = size(hyperstack, i);
end

% assemble image description
s = sprintf('ImageJ=1.51\nnimages=%d\nchannels=%d\nslices=%d\nframes=%d\nhyperstack=true\nmode=color\nloop=false\nmin=%.1f\nmax=%.1f\n', prod(d(3:5)), d(3), d(4), d(5), floor(min(hyperstack(:))*10)/10, ceil(max(hyperstack(:))*10)/10);

% open tif file for writing and set file tags
t = Tiff(file, 'w');

ts.ImageLength = d(1);
ts.ImageWidth = d(2);
ts.Photometric = Tiff.Photometric.MinIsBlack;
ts.Compression = Tiff.Compression.None;
ts.BitsPerSample = 32;
ts.SamplesPerPixel = 1;
ts.SampleFormat = Tiff.SampleFormat.IEEEFP;
ts.RowsPerStrip = 5;
ts.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky;
ts.Software = 'MATLAB';
ts.ImageDescription = s;

% loop over dimensions 3, 4, and 5
for k = 1 : d(5)
    for j = 1 : d(4)
        for i = 1 : d(3)
            frame = hyperstack(:, :, i, j, k);
            t.setTag(ts)            
            t.write(frame);
            t.writeDirectory();
        end
    end
end

% close tif file
t.close();

end

似乎 Bio-Formats 是将 >2D 数据从 Matlab 导出到 ImageJ 的唯一完全铺平且可用的方法,尽管在没有太多元数据可用且需要传输的情况下,此解决方案可能是一种轻量级的替代方案。

于 2017-11-28T12:30:16.317 回答