获取 MR 图像特定区域的强度

Getting the intensities of a certain region of an MR image

我有一张 3D MR 图像作为 NIfTI 文件 (.nii.gz)。我还有一个 'mask' 图像作为 NIfTI 文件,它只是一堆 0 和 1。这张掩码图像中的 1 代表我感兴趣的 3D MR 图像区域。

我想检索存在于掩模中的 3D MRI 图像中像素的强度(即在掩模图像文件中为 1)。我发现的唯一强度特征是 sitk.MinimumMaximumImageFilter 这不是很有用,因为它使用整个图像(而不是特定区域),并且只给出所述图像的最小值和最大值。

我也不认为 GetPixel() 函数在这种情况下对我有帮助,因为它输出的 'pixel value' 与我在 ITK-SNAP 查看器中观察到的强度不同。这是正确的吗?

在这种情况下我可以使用什么工具或功能来提供帮助?

您可能想使用 itk::Statistics::MaskedImageToHistogramFilter,然后是 min = histogram.Quantile(0, 0.0)max = histogram.Quantile(0, 1.0)。您可能需要使用比示例更多的垃圾箱。

使用itk::BinaryImageToStatisticsLabelMapFilter