selfuse

selfuse

path = '/media/hp/4T/datasets/total001-300/PCAA/'
for category in ["HGG","LGG"]:
    categoryPath = os.path.join(path, category)
    for folder in os.listdir(categoryPath):

        label = 0 if category== "LGG" else 1
        featuresDict['id'] = folder
        featuresDict['label'] = label
        print(folder,label)
        folderPath = os.path.join(categoryPath,folder)
        segPath = os.path.join(folderPath,'pca_'+folder+'_'+'seg.nii')
        print(segPath)
        for mode in ['img']:
            originPath = os.path.join(folderPath,'pca_'+folder+'_'+mode+'.nii')
            print(originPath)

027 1
/media/hp/4T/datasets/total001-300/PCAA/HGG/027/pca_027_seg.nii
/media/hp/4T/datasets/total001-300/PCAA/HGG/027/pca_027_img.nii
063 1
/media/hp/4T/datasets/total001-300/PCAA/HGG/063/pca_063_seg.nii
/media/hp/4T/datasets/total001-300/PCAA/HGG/063/pca_063_img.nii
034 1
/media/hp/4T/datasets/total001-300/PCAA/HGG/034/pca_034_seg.nii
/media/hp/4T/datasets/total001-300/PCAA/HGG/034/pca_034_img.nii
044 0
/media/hp/4T/datasets/total001-300/PCAA/LGG/044/pca_044_seg.nii
/media/hp/4T/datasets/total001-300/PCAA/LGG/044/pca_044_img.nii
132 0
/media/hp/4T/datasets/total001-300/PCAA/LGG/132/pca_132_seg.nii
/media/hp/4T/datasets/total001-300/PCAA/LGG/132/pca_132_img.nii
105 0


```bash
import numpy as np
import pandas as pd
import matplotlib as mpl
import matplotlib.pyplot as plt
import seaborn as sns
df = pd.read_csv('./features/Features_Resnet.csv')

# Draw Plot
i=2
for i in range(500):
    

    plt.figure(figsize=(13,10), dpi= 80)
    sns.violinplot(x='label', y=df.columns[i], data=df, scale='width', inner='quartile')

    # Decoration
    plt.title(str(i), fontsize=22)
    plt.show()

在这里插入图片描述

def merge_nii(folder):
    arrs = []
    for file in sorted(os.listdir(root_dir)):
        if "img" not in file:
            image = sitk.ReadImage(os.path.join(root_dir, file))
            arrs.append(sitk.GetArrayFromImage(image))

    new_arr = np.zeros(arrs[0].shape, dtype=np.uint8)
    for i in range(len(arrs)):
        new_arr[arrs[i] > 0] = 1
    
    new_image = sitk.GetImageFromArray(new_arr)
    new_image.SetSpacing(image.GetSpacing())
    new_image.SetOrigin(image.GetOrigin())
    new_image.SetDirection(image.GetDirection())
    
    sitk.WriteImage(new_image, os.path.join(root_dir, "merged.nii"))
plt.figure(figsize=(16,16))
for i in range(1,6,2):
    IMG = train_dataset[50]
    img = IMG['img']
    label = IMG['label']
    
    img = np.transpose(img, (1,2,0))
    img = img*0.5 + 0.5
    plt.subplot(3,2,i),plt.imshow(img,'gray'),plt.title('img'),plt.xticks([]),plt.yticks([])
    plt.subplot(3,2,i+1),plt.imshow(label,'gray'),plt.title('label'),plt.xticks([]),plt.yticks([])
    plt.show

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值