体位转换

最新推荐文章于 2024-01-09 15:09:39 发布

shchojj

最新推荐文章于 2024-01-09 15:09:39 发布

阅读量417

点赞数

分类专栏： python

本文链接：https://blog.csdn.net/fanre/article/details/116533235

版权

python 专栏收录该内容

16 篇文章 2 订阅

订阅专栏

该博客介绍了如何使用Euler3DTransform和Resample函数进行医学影像的中心定位和方向校正。通过get_center函数获取图像的物理中心点，然后利用rotation设置旋转矩阵，实现图像与参考图像的配准。convert_to_hfs函数则用于将图像转换到同一坐标系下，确保对比分析的准确性。

摘要由CSDN通过智能技术生成

def get_center(img):
    """
    This function returns the physical center point of a 3d sitk image
    :param img: The sitk image we are trying to find the center of
    :return: The physical center point of the image
    """
    width, height, depth = img.GetSize()
    return img.TransformIndexToPhysicalPoint((int(np.ceil(width/2)),
                                              int(np.ceil(height/2)),
                                              int(np.ceil(depth/2))))
def convert_to_hfs(image,reference_image, interpolator=sitk.sitkLinear, direction=None):
    rangefilter = sitk.MinimumMaximumImageFilter()
    rangefilter.Execute(image)
    default_value = rangefilter.GetMinimum()
    if direction is None:
        direction = reference_image.GetDirection()
    image_center = get_center(image)
    rotation = sitk.Euler3DTransform()
    rotation.SetCenter(image_center)
    rotation.SetMatrix(direction)
    resampled_image = sitk.Resample(image, reference_image, rotation,interpolator , default_value)
    return resampled_image