nrrd to nifti

def nrrd2nii(in_fn, out_fn=None):
    if not in_fn.endswith('.nrrd'):
        logger.info('skipped')
    else:
        try:
            image = sitk.ReadImage(in_fn)
            if out_fn is None:
                out_fn = in_fn.replace('nrrd', 'nii.gz')
            sitk.WriteImage(image, out_fn)
        except:

            _nrrd = nrrd.read(in_fn)
            data = _nrrd[0]
            header = _nrrd[1]
            print('header', header)
            orient_list = header['space'].split('-')
            # print(orient_list)
            orient = np.array([1, 1, 1, 1])
            if orient_list[0] == 'right':
                orient[0] = 1
            elif orient_list[0] == 'left':
                orient[0] = -1
            if orient_list[1] == 'anterior':
                orient[1] = 1
            elif orient_list[1] == 'posterior':
                orient[1] = -1
            if orient_list[2] == 'superior':
                orient[2] = 1
            elif orient_list[1] == 'inferior':
                orient[2] = -1
            space = np.diag(orient)
            # affine[:3, :3] = header['space directions']
            # affine[:3, 3] = header['space origin'].T

            affine = np.eye(4)
            affine[:3, :3] = header['space directions'] * header['spacings']
            affine[:3, 3] = np.array(list(map(float, header['origin'][1:-1].strip().split()))).T

            affine = np.dot(space, affine)
            # print(affine)
            img = nib.Nifti1Image(data, affine)
            if out_fn is None:
                out_fn = in_fn.replace('nrrd', 'nii.gz')
            # print(out_fn)
            nib.save(img, out_fn)
            print(orient_list, nib.aff2axcodes(img.affine))


        # print(header['origin'])
        # print(temp)

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值