Python 的生醫影像分析
Stephen Bailey
Instructor
空間範圍 為下列兩者的乘積：
# Calculate volume per voxel d0, d1, d2 = vol.meta['sampling'] dvoxel = d0 * d1 * d2 # Count label voxels nvoxels=ndi.sum(1, label, index=1) # Calculate volume of label volume = nvoxels * dvoxel volume
# Calculate volume per voxel d0, d1, d2 = vol.meta['sampling'] dvoxel = d0 * d1 * d2
# Count label voxels nvoxels=ndi.sum(1, label, index=1)
# Calculate volume of label volume = nvoxels * dvoxel volume
1249023
歐幾里得距離
# Create a left ventricle mask mask=np.where(labels == 1, 1, 0) # In terms of voxels d=ndi.distance_transform_edt(mask) d.max()
# Create a left ventricle mask mask=np.where(labels == 1, 1, 0) # In terms of voxels d=ndi.distance_transform_edt(mask)
d.max()
12.3847
# In terms of space d=ndi.distance_transform_edt(mask, sampling=vol.meta['sampling']) d.max()
5.8038
com=ndi.center_of_mass(vol, labels, index=1) com
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()
Preparing Video For Download...
