users > sri24 tzo116 atlas -- missing labels
Showing 1-2 of 2 posts
Display:
Results per page:
Jan 27, 2012  03:01 PM | Yaroslav Halchenko
sri24 tzo116 atlas -- missing labels
It seems that some labels defined in SRI24-tzo116plus.txt are not defined in tzo116plus.nii.gz and vise versa. how come? ;)

here is the output of a little python snippet listed below

(array([83, 84]), array([98, 98]), array([55, 55]), array([0, 0]))
Labels present in the map but not in the labels file:  set([422, 424, 426, 428, 430, 432, 434, 476, 478])
Labels present in the labels file but not in the map:  set([227, 201, 203, 205, 413, 219, 221])

import nibabel as nib;
import numpy as np

atl_tzo116 = "/usr/share/data/sri24-atlas/tzo116plus.nii.gz"
atl_tzo116_labelsfile = '/usr/share/data/sri24-atlas/SRI24-tzo116plus.txt'
atl_tzo116_labels = np.loadtxt(atl_tzo116_labelsfile, dtype=str)
atl_tzo116_labels = dict([(int(i), l) for i,l in atl_tzo116_labels[:, :2]])

ni = nib.load(atl_tzo116)
m = ni.get_data()

mset = set(np.unique(m))
lset = set(np.unique(atl_tzo116_labels.keys()))

print np.where(m == 432)
# so we have voxels labels as 432 which is not present in the labels file

print "Labels present in the map but not in the labels file: ", mset - lset
print "Labels present in the labels file but not in the map: ", lset - mset
Jan 27, 2012  03:01 PM | Torsten Rohlfing - Google LLC
RE: sri24 tzo116 atlas -- missing labels
Created a tracker item for SRI24:

http://www.nitrc.org/tracker/index.php?f...