如何在MATLAB中提取三维阵元的随机位移强度差？

Question

我正在对三维医学图像进行图像分割。在一篇相关的论文中，作者提取了一个我不知道该如何计算的特征，f(x, u) = I(x + u) − I(x)是x和mm^3中偏移量u之间的随机移动强度差。如果我们认为医学图像的pixel_spacing是0.36x0.36x0.7 mm^3，那么如何在内存和计算效率上为u=2 mm^3提取和获取整个卷的这一特征？谢谢你的帮助。

Answer 1

正如@user630317的评论所指出的那样，对于3D图像，u应该是一个3D向量。

此外，图像是‘体化’i，e，离散。因此，x是一个包含整数条目和单位为体素(而不是mm)的向量。U的元素可能是整数，也可能不是整数，但有单位作为体素。如果u有分数元，那么就必须对分数移位进行插值，从而丢失信息。

matlab函数可以同时完成整数移位和分数移位。

假设，我是一个三维数组，在matlab中表示三维体积.设，u是平移的三维矢量(单位为体素)。

那么翻译后的图像将是I2 = imtranslate(I，u)

这是一个具有3D数组的示例代码：

I = randn(3,3,3);
u = [1.5,1,1];
I2 = imtranslate(I,u);