栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > Python

nii文件python可视化并转化为gif

Python 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

nii文件python可视化并转化为gif

今儿刚好用的医学图像处理,看了一些博客,就封装了一下,把三维nii ,转化为 gif,懒得废话了,直接看代码的注释吧:

import cv2
import nibabel as nib
import numpy as np
from PIL import Image


def show_nii_gif(nii_file, save_path="array.gif", show_frame=True):
    """
    @Brife:
        将 shape=(arr_len, h, w) 的 nii 文件转化为 gif
    @Param:
        nii_file   : nii 文件的路径
        show_frame : 是否展示这是第几帧, 默认 True
        save_path  : 默认为当前路径的 array.gif 文件
    """
    
    # 调包加载 nii 文件
    img=nib.load(nii_file)
    # 转化成 numpy.ndarray 的格式
    img_arr = img.get_fdata()
    img_arr = np.squeeze(img_arr)
    
    # 找到最大最小值, 方便之后归一化
    img_max, img_min = img_arr.max(), img_arr.min()
    
    # 归一化
    img_arr = (img_arr - img_min) / (img_max - img_min) * 255
    img_arr = img_arr.astype(np.uint8)
    
    # 将单通道的灰度图转化为RGB三通道的灰色图, (不转化这一步, 没法写字)
    img_RGB_list = [cv2.cvtColor(img, cv2.COLOR_GRAY2BGR) for img in img_arr]
    
    assert 3==len(img_arr.shape) # 如果是别的shape的, 那处理不了
    
    arr_len, h, w = img_arr.shape
    
    if show_frame:
        # 在每一帧上打上这是第几帧
        for i in range(arr_len):
            
            img_RGB_list[i] = cv2.putText(img_RGB_list[i], 
                                     '{:>03d}/{:>03d}'.format(i+1, arr_len),
                                     (50, 50),
                                     cv2.FONT_HERSHEY_COMPLEX,
                                     1,
                                     (255, 255, 255),
                                     2)
    
    # 将所有的 ndarray 转化为 PIL.Image 格式
    imgs = [Image.fromarray(img) for img in img_RGB_list]
    
    # 保存
    # duration is the number of milliseconds between frames; this is 40 frames per second
    imgs[0].save(save_path, save_all=True, append_images=imgs[1:], duration=50, loop=0)


show_nii_gif(nii_file='IMG_0002.nii')

展示一哈:

参考了:

Python读取nii文件
https://blog.csdn.net/PKnotnull/article/details/80841410

用openCV在 array 上写字:
https://blog.csdn.net/GAN_player/article/details/78155283

灰度图直接转化为三通道灰度图:
https://blog.csdn.net/llh_1178/article/details/77833447

Python字符串格式化,"{:>03d}" 死活记不住hhh [向右对齐,并用0补充]
https://www.runoob.com/python/att-string-format.html

将ndarray的照片生成git
https://stackoverflow.com/questions/63047707/how-do-i-convert-a-numpy-array-to-a-gif

Matplotlib 保存gif (这次没有用上,几年前第一次用,就感觉是个老麻烦的东西,hhh)
https://blog.csdn.net/qq_28888837/article/details/85778395

Python skimage 库安装:
https://blog.csdn.net/j_starry/article/details/100192206

转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/757459.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号