ifft 重建后的图像类型

信息处理 matlab 图片 IFFT 重建
2022-01-26 20:23:07

我使用 ifft 和平方根方法从 MRI k 空间重建图像。

imRef = ifftshift(ifftshift(ifft(ifft( ifftshift(ifftshift( kspacedata,1),2),[],1),[],2),1),2); 
imSOS    = squeeze(sqrt(sum( abs(angle(imRef)).^2, 3))); 
imagesc(abs(imSOS))

k = abs(imSOS); 
disp(class(k)); %double

图像矩阵k的类别显示为double

图像矩阵k是索引图像还是强度图像?

重建图像的样本值k如下所示。

4.2753    4.9807    4.5435    5.6548    6.1303    5.0229
3.3805    4.6260    5.1594    5.2692    4.1187    4.4885
5.8990    5.5275    4.3493    5.6182    6.7237    4.3071
6.4885    4.6861    4.4086    3.5034    5.2378    4.7466
6.1099    4.6995    4.1673    4.7408    3.8915    5.7531
5.4006    3.1289    5.6541    5.8782    4.6568    3.8166

通过引用http://in.mathworks.com/help/matlab/creating_plots/image-types.html,我了解到对于索引图像,数据矩阵由整数表示,整数是颜色图的索引。如果是强度图像,我觉得数据范围介于 0 和 1 之间。如果k表示索引图像,我怎么知道与之关联的颜色图?

谢谢。

2个回答

FFT 和 IFFT 是线性算子,因此,结果仅在线性强度空间中才有意义,而不是在索引到非线性映射空间时才有意义。

在 MRI 中,没有内在的“彩色图”。与计算机断层扫描相反,没有类似于 Hounsfield 单位。因此,您只有相对强度信息。

通常,如果您自己重建图像,您将在任意范围内获得双精度值。如果您将其与 MR 机器中的 DICOM 数据进行比较,您会注意到 DICOM 图像被重新缩放为整数值,比如说在 0 到 4095 之间。重新缩放因子只是一个数字,用于将信号拟合在预定义的范围内 -这里再次没有额外的定量信息。