如何对基于毫米的 3D 图像进行平滑处理？

Question

我正在模拟一篇论文，图像平滑的比例是基于毫米的。 3D 图像的格式是 DICOM。例如，通常使用窗口 size 和缩放 sX 进行平滑处理/em>，正在做如下：

f1 = fspecial('gaussian',[size,size],s);
Smooth1 = imfilter(X,f1);

有谁知道我怎样才能平滑它的刻度以毫米为单位？

我应该根据毫米而不是像素来更改窗口的大小吗？我怎样才能做到这一点？

Answer 1

如果您希望以毫米为单位定义高斯滤波器的形状，您可以使用 DICOM 元数据中的PixelSpacing 属性将毫米转换为像素。 PixelSpacing的单位是mm/pixel

info = dicominfo(filename);
X = double(dicomread(info));

scale_mm = 3;                                           % Signma of the filter in mm

scale_px = scale_mm / double(info.PixelSpacing(1));     % Converted to pixels

% Create a Gaussian filter using this sigma value
f1 = fspecial('gaussian', [height, width], scale_px);

% Apply the smoothing
smooth1 = imfilter(X, f1);