Análise de Imagens Biomédicas em Python
Stephen Bailey
Instructor

Extensão espacial é o produto de:

# Calcular volume por voxel d0, d1, d2 = vol.meta['sampling'] dvoxel = d0 * d1 * d2# Contar voxels do rótulo nvoxels=ndi.sum(1, label, index=1)# Calcular volume do rótulo volume = nvoxels * dvoxel volume
1249023

Distância euclidiana
# Criar uma máscara do ventrículo esquerdo mask=np.where(labels == 1, 1, 0) # Em termos de voxels d=ndi.distance_transform_edt(mask)d.max()
12.3847
# Em termos de espaço
d=ndi.distance_transform_edt(mask,
sampling=vol.meta['sampling'])
d.max()
5.8038
com=ndi.center_of_mass(vol, labels, index=1)com
(5.5235, 128.0590, 128.0993)
plt.imshow(vol[5], cmap='gray')
plt.scatter(com[2], com[1])
plt.show()

Análise de Imagens Biomédicas em Python