2017-12-20 26 views
-1

저는 뇌 세분화 프로젝트를 진행 중이며 (BRATS2015에서) MRI를 전처리해야합니다. 그리고 내 3D 사진의 막대 그래프를보고 싶습니다. 내 코드 :3D 이미지의 히스토그램을 얻는 방법

#-*- coding:utf-8 -*- 
import SimpleITK as sitk 
import numpy as np 
from matplotlib import pyplot 
import matplotlib.pyplot as plt 

url = r'H:\myfile\VSD.Brain.XX.O.MR_T1.54525.mha' 
image = sitk.ReadImage(url) 
result = sitk.GetArrayFromImage(image) 
print(type(result)) 
print(result.shape) 
plt.figure('historgram') 
result = result.flatten() 
n, bins, patches = plt.hist(result, bins=256, normed=0, facecolor='red', alpha=0.75) 
plt.show() 

그러나 일부 틀렸다는 것이다 이렇게 싶어 : enter image description here

답변

0

오, 나는 그것을 처리합니다.

#-*- coding:utf-8 -*- 
import SimpleITK as sitk 
import numpy as np 
from matplotlib import pyplot 
import matplotlib.pyplot as plt 

url = r'H:\deepfortest\brats2015_train\HGG\brats_2013_pat0003_1\VSD.Brain.XX.O.MR_T1.54525.mha' 
image = sitk.ReadImage(url) 
result = sitk.GetArrayFromImage(image) 
print(type(result)) 
print(result.shape) 
plt.figure('historgram') 
result = result.flatten() 
n, bins, patches = plt.hist(result, bins=256, range= (1,result.max()),normed=0, facecolor='red', alpha=0.75,histtype = 'step') 
plt.show()