In [1]:
%matplotlib inline
import math
import numpy
from matplotlib import pyplot
from allensdk.core.mouse_connectivity_cache import MouseConnectivityCache

# The manifest file is a simple JSON file that keeps track of all of
# the data that has already been downloaded onto the hard drives.
# If you supply a relative path, it is assumed to be relative to your
# current working directory.
mcc = MouseConnectivityCache(manifest_file='data/manifest.json')
# grab the Ontology instance
ontology = mcc.get_ontology()

In [2]:
ontology.df


Out[2]:
acronym atlas_id color_hex_triplet depth failed failed_facet graph_id graph_order hemisphere_id id ... neuro_name_structure_id neuro_name_structure_id_path ontology_id parent_structure_id safe_name sphinx_id st_level structure_id_path structure_name_facet weight
id
997 root -1 FFFFFF 0 False 734881840 1 0 3 997 ... NaN NaN 1 NaN root 1 NaN /997/ 385153371 8690
8 grey 0 BFDAE3 1 False 734881840 1 1 3 8 ... NaN NaN 1 997 Basic cell groups and regions 2 NaN /997/8/ 2244697386 8690
567 CH 70 B0F0FF 2 False 734881840 1 2 3 567 ... NaN NaN 1 8 Cerebrum 3 NaN /997/8/567/ 2878815794 8690
688 CTX 85 B0FFB8 3 False 734881840 1 3 3 688 ... NaN NaN 1 567 Cerebral cortex 4 NaN /997/8/567/688/ 3591311804 8690
695 CTXpl 86 70FF70 4 False 734881840 1 4 3 695 ... NaN NaN 1 688 Cortical plate 5 NaN /997/8/567/688/695/ 3945900931 8690
315 Isocortex 746 70FF71 5 False 734881840 1 5 3 315 ... NaN NaN 1 695 Isocortex 6 NaN /997/8/567/688/695/315/ 2323732626 8690
184 FRP 871 268F45 6 False 734881840 1 6 3 184 ... NaN NaN 1 315 Frontal pole cerebral cortex 7 NaN /997/8/567/688/695/315/184/ 2565719060 8690
68 FRP1 998 268F45 7 False 734881840 1 7 3 68 ... NaN NaN 1 184 Frontal pole layer 1 8 NaN /997/8/567/688/695/315/184/68/ 1397862467 8690
667 FRP2/3 1073 268F45 7 False 734881840 1 8 3 667 ... NaN NaN 1 184 Frontal pole layer 2/3 9 NaN /997/8/567/688/695/315/184/667/ 4268100038 8690
500 MO 203 1F9D5A 6 False 734881840 1 9 3 500 ... NaN NaN 1 315 Somatomotor areas 10 NaN /997/8/567/688/695/315/500/ 356591023 8690
107 MO1 1003 1F9D5A 7 False 734881840 1 10 3 107 ... NaN NaN 1 500 Somatomotor areas Layer 1 11 NaN /997/8/567/688/695/315/500/107/ 900612548 8690
219 MO2/3 1017 1F9D5A 7 False 734881840 1 11 3 219 ... NaN NaN 1 500 Somatomotor areas Layer 2/3 12 NaN /997/8/567/688/695/315/500/219/ 1075749695 8690
299 MO5 1027 1F9D5A 7 False 734881840 1 12 3 299 ... NaN NaN 1 500 Somatomotor areas Layer 5 13 NaN /997/8/567/688/695/315/500/299/ 851674589 8690
644 MO6a 787 1F9D5A 7 False 734881840 1 13 3 644 ... NaN NaN 1 500 Somatomotor areas Layer 6a 14 NaN /997/8/567/688/695/315/500/644/ 1003126892 8690
947 MO6b 825 1F9D5A 7 False 734881840 1 14 3 947 ... NaN NaN 1 500 Somatomotor areas Layer 6b 15 NaN /997/8/567/688/695/315/500/947/ 2730742230 8690
985 MOp 830 1F9D5A 7 False 734881840 1 15 3 985 ... NaN NaN 1 500 Primary motor area 16 NaN /997/8/567/688/695/315/500/985/ 1852742012 8690
320 MOp1 888 1F9D5A 8 False 734881840 1 16 3 320 ... NaN NaN 1 985 Primary motor area Layer 1 17 NaN /997/8/567/688/695/315/500/985/320/ 571569106 8690
943 MOp2/3 966 1F9D5A 8 False 734881840 1 17 3 943 ... NaN NaN 1 985 Primary motor area Layer 2/3 18 NaN /997/8/567/688/695/315/500/985/943/ 2488357079 8690
648 MOp5 929 1F9D5A 8 False 734881840 1 18 3 648 ... NaN NaN 1 985 Primary motor area Layer 5 19 NaN /997/8/567/688/695/315/500/985/648/ 628930507 8690
844 MOp6a 1095 1F9D5A 8 False 734881840 1 19 3 844 ... NaN NaN 1 985 Primary motor area Layer 6a 20 NaN /997/8/567/688/695/315/500/985/844/ 3473508879 8690
882 MOp6b 1100 1F9D5A 8 False 734881840 1 20 3 882 ... NaN NaN 1 985 Primary motor area Layer 6b 21 NaN /997/8/567/688/695/315/500/985/882/ 1442896821 8690
993 MOs 831 1F9D5A 7 False 734881840 1 21 3 993 ... NaN NaN 1 500 Secondary motor area 22 NaN /997/8/567/688/695/315/500/993/ 1043755260 8690
656 MOs1 930 1F9D5A 8 False 734881840 1 22 3 656 ... NaN NaN 1 993 Secondary motor area layer 1 23 NaN /997/8/567/688/695/315/500/993/656/ 4032915713 8690
962 MOs2/3 1110 1F9D5A 8 False 734881840 1 23 3 962 ... NaN NaN 1 993 Secondary motor area layer 2/3 24 NaN /997/8/567/688/695/315/500/993/962/ 3274108161 8690
767 MOs5 944 1F9D5A 8 False 734881840 1 24 3 767 ... NaN NaN 1 993 Secondary motor area layer 5 25 NaN /997/8/567/688/695/315/500/993/767/ 4144803096 8690
1021 MOs6a 1117 1F9D5A 8 False 734881840 1 25 3 1021 ... NaN NaN 1 993 Secondary motor area layer 6a 26 NaN /997/8/567/688/695/315/500/993/1021/ 3489757563 8690
1085 MOs6b 1125 1F9D5A 8 False 734881840 1 26 3 1085 ... NaN NaN 1 993 Secondary motor area layer 6b 27 NaN /997/8/567/688/695/315/500/993/1085/ 1225271489 8690
453 SS 339 188064 6 False 734881840 1 27 3 453 ... NaN NaN 1 315 Somatosensory areas 28 NaN /997/8/567/688/695/315/453/ 3016301862 8690
12993 SS1 NaN 188064 7 False 734881840 1 28 3 12993 ... NaN NaN 1 453 Somatosensory areas layer 1 29 NaN /997/8/567/688/695/315/453/12993/ 544207135 8690
12994 SS2/3 NaN 188064 7 False 734881840 1 29 3 12994 ... NaN NaN 1 453 Somatosensory areas layer 2/3 30 NaN /997/8/567/688/695/315/453/12994/ 3920828838 8690
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
116 chfl 721 AAAAAA 3 False 734881840 1 1204 3 116 ... NaN NaN 1 81 choroid fissure 1205 NaN /997/73/81/116/ 3192898784 8690
124 IVF 722 AAAAAA 2 False 734881840 1 1205 3 124 ... NaN NaN 1 73 interventricular foramen 1206 NaN /997/73/124/ 1644273448 8690
129 V3 723 AAAAAA 2 False 734881840 1 1206 3 129 ... NaN NaN 1 73 third ventricle 1207 NaN /997/73/129/ 562407244 8690
140 AQ 724 AAAAAA 2 False 734881840 1 1207 3 140 ... NaN NaN 1 73 cerebral aqueduct 1208 NaN /997/73/140/ 3509002335 8690
145 V4 725 AAAAAA 2 False 734881840 1 1208 3 145 ... NaN NaN 1 73 fourth ventricle 1209 NaN /997/73/145/ 2366263365 8690
153 V4r 726 AAAAAA 3 False 734881840 1 1209 3 153 ... NaN NaN 1 145 lateral recess 1210 NaN /997/73/145/153/ 4006455201 8690
164 c 727 AAAAAA 2 False 734881840 1 1210 3 164 ... NaN NaN 1 73 central canal spinal cord/medulla 1211 NaN /997/73/164/ 3758121442 8690
1024 grv 693 AAAAAA 1 False 734881840 1 1211 3 1024 ... NaN NaN 1 997 grooves 1212 NaN /997/1024/ 1816833424 8690
1032 grv of CTX 694 AAAAAA 2 False 734881840 1 1212 3 1032 ... NaN NaN 1 1024 grooves of the cerebral cortex 1213 NaN /997/1024/1032/ 4010583074 8690
1055 eg 697 AAAAAA 3 False 734881840 1 1213 3 1055 ... NaN NaN 1 1032 endorhinal groove 1214 NaN /997/1024/1032/1055/ 2757546022 8690
1063 hf 698 AAAAAA 3 False 734881840 1 1214 3 1063 ... NaN NaN 1 1032 hippocampal fissure 1215 NaN /997/1024/1032/1063/ 1845828307 8690
1071 rf 699 AAAAAA 3 False 734881840 1 1215 3 1071 ... NaN NaN 1 1032 rhinal fissure 1216 NaN /997/1024/1032/1071/ 915997995 8690
1078 ri 700 AAAAAA 3 False 734881840 1 1216 3 1078 ... NaN NaN 1 1032 rhinal incisure 1217 NaN /997/1024/1032/1078/ 1408661066 8690
1040 grv of CBX 695 AAAAAA 2 False 734881840 1 1217 3 1040 ... NaN NaN 1 1024 grooves of the cerebellar cortex 1218 NaN /997/1024/1040/ 306524468 8690
1087 pce 701 AAAAAA 3 False 734881840 1 1218 3 1087 ... NaN NaN 1 1040 precentral fissure 1219 NaN /997/1024/1040/1087/ 2084828473 8690
1095 pcf 702 AAAAAA 3 False 734881840 1 1219 3 1095 ... NaN NaN 1 1040 preculminate fissure 1220 NaN /997/1024/1040/1095/ 3439430666 8690
1103 pri 703 AAAAAA 3 False 734881840 1 1220 3 1103 ... NaN NaN 1 1040 primary fissure 1221 NaN /997/1024/1040/1103/ 1359711822 8690
1112 psf 704 AAAAAA 3 False 734881840 1 1221 3 1112 ... NaN NaN 1 1040 posterior superior fissure 1222 NaN /997/1024/1040/1112/ 706594763 8690
1119 ppf 705 AAAAAA 3 False 734881840 1 1222 3 1119 ... NaN NaN 1 1040 prepyramidal fissure 1223 NaN /997/1024/1040/1119/ 90069120 8690
3 sec 707 AAAAAA 3 False 734881840 1 1223 3 3 ... NaN NaN 1 1040 secondary fissure 1224 NaN /997/1024/1040/3/ 1603026106 8690
11 plf 708 AAAAAA 3 False 734881840 1 1224 3 11 ... NaN NaN 1 1040 posterolateral fissure 1225 NaN /997/1024/1040/11/ 3968698314 8690
18 nf 709 AAAAAA 3 False 734881840 1 1225 3 18 ... NaN NaN 1 1040 nodular fissure 1226 NaN /997/1024/1040/18/ 3778172686 8690
25 sif 710 AAAAAA 3 False 734881840 1 1226 3 25 ... NaN NaN 1 1040 simple fissure 1227 NaN /997/1024/1040/25/ 3504067713 8690
34 icf 711 AAAAAA 3 False 734881840 1 1227 3 34 ... NaN NaN 1 1040 intercrural fissure 1228 NaN /997/1024/1040/34/ 1415991612 8690
43 apf 712 AAAAAA 3 False 734881840 1 1228 3 43 ... NaN NaN 1 1040 ansoparamedian fissure 1229 NaN /997/1024/1040/43/ 2452402730 8690
49 ipf 713 AAAAAA 3 False 734881840 1 1229 3 49 ... NaN NaN 1 1040 intraparafloccular fissure 1230 NaN /997/1024/1040/49/ 2654107150 8690
57 pms 714 AAAAAA 3 False 734881840 1 1230 3 57 ... NaN NaN 1 1040 paramedian sulcus 1231 NaN /997/1024/1040/57/ 3972977495 8690
65 pfs 715 AAAAAA 3 False 734881840 1 1231 3 65 ... NaN NaN 1 1040 parafloccular sulcus 1232 NaN /997/1024/1040/65/ 771629690 8690
624 IPF 926 AAAAAA 2 False 734881840 1 1232 3 624 ... NaN NaN 1 1024 Interpeduncular fossa 1233 NaN /997/1024/624/ 1476705011 8690
304325711 retina NaN 7F2E7E 1 False 734881840 1 1233 3 304325711 ... NaN NaN 1 997 retina 1234 NaN /997/304325711/ 3295290839 8690

1234 rows × 21 columns


In [3]:
tree_levels = [[] for i in range(0, 12)]
childs_map = dict()

for i in range(0, 12):
    level_ids = []
    for a in ontology.df.values:
        if (i == 0 and math.isnan(a[14])) or (i > 0 and a[14] in tree_levels[i - 1]):
            level_ids.append(a[9])
            childs_map[a[9]] = []
    tree_levels[i] = level_ids
    
for i in range(11, 0, -1):
    for a in ontology.df.values:
        if a[14] in tree_levels[i - 1]:
            childs_map[a[14]].append(a[9])
            childs_map[a[14]].extend(childs_map[a[9]])

parents = set([a[14] for a in ontology.df.values if not math.isnan(a[14])])
print "PARENTS:", len(parents)
leafs = [(a[0], a[9]) for a in ontology.df.values if a[9] not in parents]
leaf_ids = [l[1] for l in leafs]
print "LEAFS:", len(leafs), " of ", len(ontology.df.values), "MIN:", min(leaf_ids), "MAX:", max(leaf_ids)

def find_matching_tree_level(node_id):
    for i, arr_ids in enumerate(tree_levels):
        if node_id in arr_ids:
            return i+1
    return -1


PARENTS: 278
LEAFS: 956  of  1234 MIN: 1 MAX: 312782652

In [4]:
template, template_info = mcc.get_template_volume()
annot, annot_info = mcc.get_annotation_volume()

# show that slice of all volumes side-by-side
f, ((ax1, ax2)) = pyplot.subplots(1, 2, figsize=(15, 6))
# pick a slice to show
slice_idx = 264

ax1.imshow(template[slice_idx,:,:], cmap='gray', aspect='equal', vmin=template.min(), vmax=template.max())
ax1.set_title("registration template")

ax2.imshow(annot[slice_idx,:,:], cmap='gray', aspect='equal', vmin=0, vmax=2000)
ax2.set_title("annotation volume")


Out[4]:
<matplotlib.text.Text at 0x105de1c90>

In [5]:
print template.min(), template.max(), annot.min(), annot.max()


0 516 0 312782652

In [6]:
annot_info


Out[6]:
{u'dimension': 3,
 u'encoding': 'gzip',
 u'endian': 'little',
 u'keyvaluepairs': {},
 u'kinds': ['domain', 'domain', 'domain'],
 u'sizes': [528, 320, 456],
 u'space': 'left-posterior-superior',
 u'space directions': [['25', '0', '0'], ['0', '25', '0'], ['0', '0', '25']],
 u'space origin': ['0', '0', '0'],
 u'type': 'unsigned int'}

In [7]:
found_ann = numpy.unique(annot)
print len(numpy.unique(template)), "distinct annotations", len(found_ann)


517 distinct annotations 589

In [8]:
all_region_idx = [a[9] for a in ontology.df.values]
background_idx = [ann for ann in found_ann if ann not in all_region_idx]
print "Background Value(s)", background_idx


Background Value(s) [0]

In [9]:
conn2voxels = {a[9]: (find_matching_tree_level(a[9]), a[0], (annot == a[9]).sum())
               for a in ontology.df.values}

In [10]:
mapped_regions = []
not_mapped_regions= []
for l in range(1, 11):
    ## ID, LABEL, VOXELS_NO
    level_data = [(i, a[1], a[2]) for i, a in conn2voxels.iteritems() if a[0]==l and a[2] >0]
    no_voxels = set([a[1] for i, a in conn2voxels.iteritems() if a[0]==l and a[2] ==0])
    voxels = set([i for i, a in conn2voxels.iteritems() if a[0]==l and a[2] >0])
    mapped_regions.extend(voxels)
    not_mapped_regions.extend(no_voxels)
    print "================================"
    print "LEVEL", l 
    print "Regions with no associated voxels:", len(no_voxels), no_voxels
    print ""
    print "Regions with at least one voxel associated:", len(level_data)
    print level_data
    
print "Total no of regions with associated voxels:", len(mapped_regions), "total not associated", len(not_mapped_regions)


================================
LEVEL 1
Regions with no associated voxels: 0 set([])

Regions with at least one voxel associated: 1
[(997, 'root', 40)]
================================
LEVEL 2
Regions with no associated voxels: 3 set(['grv', 'retina', 'grey'])

Regions with at least one voxel associated: 2
[(73, 'VS', 18385), (1009, 'fiber tracts', 236648)]
================================
LEVEL 3
Regions with no associated voxels: 11 set(['IVF', 'CH', 'cbf', 'cm', 'eps', 'mfbs', 'BS', 'IPF', 'grv of CBX', 'lfbs', 'grv of CTX'])

Regions with at least one voxel associated: 6
[(81, 'VL', 159548), (129, 'V3', 82850), (140, 'AQ', 28584), (512, 'CB', 1), (145, 'V4', 50656), (164, 'c', 80)]
================================
LEVEL 4
Regions with no associated voxels: 50 set(['vsp', 'cbc', 'ipf', 'ppf', 'pms', 'eg', 'CTX', 'SEZ', 'cbp', 'sec', 'epsc', 'IIIn', 'HB', 'Vn', 'apf', 'mfbc', 'vrt', 'icf', 'nf', 'tn', 'chpl', 'rf', 'XIIn', 'sif', 'RC', 'pcf', 'ri', 'pce', 'CNU', 'chfl', 'sst', 'pri', 'mfsbshy', 'VIn', 'plf', 'pfs', 'rst', 'Xn', 'IB', 'hf', 'psf', 'IXn', 'CBN', 'drt', 'XIn', 'lfbst', 'ctb', 'In', 'VIIIn', 'CBX'])

Regions with at least one voxel associated: 11
[(728, 'arb', 451669), (313, 'MB', 353356), (776, 'cc', 652685), (784, 'cst', 4959), (798, 'VIIn', 6612), (848, 'IIn', 3786), (877, 'tsp', 10009), (949, 'von', 17), (863, 'rust', 22757), (153, 'V4r', 53228), (911, 'IVn', 451)]
================================
LEVEL 5
Regions with no associated voxels: 71 set(['cett', 'fp', 'mfb', 'MBsta', 'ee', 'ptf', 'MBmot', 'ttp', 'pmx', 'iVIIn', 'CBXmo', 'tspc', 'sop', 'CBXgr', 'tp', 'cbt', 'tspd', 'rrt', 'nst', 'poc', 'apd', 'pap', 'CBXpu', 'scrt', 'cte', 'pvbh', 'mfbse', 'csc', 'sup', 'rstm', 'rstl', 'cstu', 'lotg', 'tct', 'IVd', 'fxs', 'cstc', 'cVIIIn', 'ntt', 'lab', 'sV', 'mfbsma', 'MBsen', 'step', 'aot', 'snp', 'HEM', 'php', 'pyd', 'vlt', 'bsc', 'shp', 'mfbst', 'stf', 'crt', 'fa', 'cvb', 'sttv', 'VERM', 'ccs', 'ccr', 'gVIIn', 'sttl', 'srp', 'stp', 'CTXpl', 'mfbsm', 'mtc', 'ccg', 'cpt', 'svp'])

Regions with at least one voxel associated: 34
[(6, 'int', 155764), (62, 'mlf', 67923), (78, 'mcp', 59418), (91, 'IP', 55207), (93, 'moV', 3952), (117, 'och', 18309), (125, 'opt', 43331), (158, 'pc', 3399), (14, 'im', 106), (190, 'py', 34577), (237, 'ts', 652), (301, 'st', 45648), (326, 'scp', 28914), (341, 'smd', 1000), (354, 'MY', 575501), (397, 'vtd', 4164), (413, 'vVIIIn', 15391), (477, 'STR', 175862), (549, 'TH', 295591), (703, 'CTXsp', 105638), (771, 'P', 384669), (803, 'PAL', 111434), (846, 'DN', 21468), (884, 'amc', 4933), (900, 'aco', 56240), (908, 'act', 13189), (924, 'cpd', 55637), (940, 'cing', 66890), (989, 'FN', 32651), (1016, 'onl', 251153), (1060, 'dtd', 2632), (1092, 'em', 8312), (1097, 'HY', 232031), (1123, 'icp', 47841)]
================================
LEVEL 6
Regions with no associated voxels: 64 set(['SCm', 'vc', 'jrb', 'hht', 'fxprg', 'DORpm', 'dl', 'ec', 'ab', 'per', 'P-sat', 'DT', 'PALv', 'SCs', 'PVR', 'dlf', 'aolt', 'sctv', 'PVZ', 'dtt', 'sctd', 'EP', 'P-sen', 'PRT', 'fpr', 'PIS', 'LZ', 'LSX', 'lotd', 'pvbt', 'phpl', '6b', 'dcm', 'fxpo', 'SNl', 'STRv', 'dscp', 'dc', 'DORsm', 'das', 'hc', 'BLA', 'PALc', 'STRd', 'ME', 'phpm', 'uf', 'phpd', 'cct', 'supv', 'PALd', 'MT', 'P-mot', 'PALm', 'ias', 'MY-mot', 'bct', 'supa', 'phpv', 'MY-sat', 'supd', 'MY-sen', 'RAmb', 'MEZ'])

Regions with at least one voxel associated: 60
[(4, 'IC', 357773), (35, 'III', 2303), (66, 'LT', 382), (115, 'IV', 466), (128, 'MRN', 311620), (131, 'LA', 41076), (214, 'RN', 41425), (231, 'AT', 998), (246, 'RR', 27478), (271, 'SAG', 5980), (278, 'sAMY', 15406), (315, 'Isocortex', 6), (319, 'BMA', 69327), (374, 'SNc', 13715), (381, 'SNr', 107010), (530, 'df', 4736), (580, 'NB', 8983), (583, 'CLA', 31625), (603, 'fi', 90210), (616, 'CUN', 18893), (633, 'cic', 1282), (658, 'll', 60621), (673, 'mp', 1377), (681, 'mtg', 5319), (690, 'mtt', 10515), (697, 'ml', 34820), (698, 'OLF', 275171), (460, 'MEV', 2226), (749, 'VTA', 51455), (753, 'pm', 59), (757, 'VTN', 2406), (780, 'PA', 20241), (794, 'sptV', 101999), (795, 'PAG', 283394), (802, 'sm', 17327), (841, 'tb', 17735), (874, 'PBG', 5430), (912, 'LING', 7942), (920, 'CENT', 250358), (928, 'CUL', 448510), (936, 'DEC', 239591), (944, 'FOTU', 58946), (951, 'PYR', 58832), (957, 'UVU', 140142), (968, 'NOD', 80661), (975, 'EW', 2053), (1007, 'SIM', 351345), (1017, 'AN', 733048), (1025, 'PRM', 315291), (1033, 'COPY', 155348), (1049, 'FL', 84339), (1052, 'PPN', 29476), (1089, 'HPF', 171733), (1041, 'PFL', 356892), (705, 'mtV', 840), (595, 'fr', 10990), (482, 'bic', 21795), (611, 'hbc', 3617), (466, 'alv', 36066), (665, 'lot', 81071)]
================================
LEVEL 7
Regions with no associated voxels: 108 set(['SIMmo', 'PRNv', 'FLpu', 'EPI', 'sct', 'UVUpu', 'LINGpu', 'DCN', 'PYRgr', 'VNC', 'AOB', 'NODgr', 'PERI', 'LINGgr', 'CUL4', 'CUL5', 'MDRN', 'oct', 'SPF', 'GU', 'DECmo', 'db', 'NODpu', 'COPYgr', 'PFLpu', 'AHA', 'RHP', 'VIS', 'PHY', 'LAT', 'VENT', 'PRMgr', 'CUL4, 5', 'ACVI', 'PYRmo', 'AMB', 'NODmo', 'PGRN', 'PFLgr', 'AI', 'COPYpu', 'DECgr', 'PL', 'FOTUmo', 'CENT3', 'CENT2', 'SCzo', 'EV', 'INV', 'PRMmo', 'ECO', 'SSN', 'FLgr', 'ACA', 'SIMpu', 'ILA', 'ATN', 'TEa', 'ILM', 'ANcr2', 'ANcr1', 'COA', 'UVUmo', 'PMR', 'NLOT', 'rct', 'iaf', 'AUD', 'MRNp', 'FRP', 'DECpu', 'MRNm', 'mct', 'LTN', 'PFLmo', 'MED', 'VISC', 'SS', 'FOTUgr', 'UVUgr', 'PSCH', 'BMAp', 'z', 'ORB', 'MO', 'LINGmo', 'BMAa', 'FLmo', 'PYRpu', 'PTLp', 'MSC', 'CN', 'HIP', 'SOC', 'MBO', 'SIMgr', 'LS', 'GENv', 'PRMpu', 'MTN', 'GENd', 'FOTUpu', 'ICc', 'ICe', 'ICd', 'MRNmg', 'COPYmo', 'RSP'])

Regions with at least one voxel associated: 143
[(7, 'PSV', 70385), (10, 'SCig', 133547), (12, 'IF', 2156), (17, 'SCiw', 134245), (23, 'AAA', 37040), (26, 'SCdg', 78038), (30, 'PVa', 816), (38, 'PVH', 15889), (42, 'SCdw', 21888), (50, 'PRC', 5291), (56, 'ACB', 263008), (63, 'PVHd', 3051), (67, 'INC', 3513), (72, 'ADP', 2188), (83, 'IO', 33071), (88, 'AHN', 49955), (100, 'IPN', 22723), (106, 'ISN', 456), (118, 'PVi', 6785), (126, 'PVp', 8772), (133, 'PVpo', 5380), (136, 'IRN', 153949), (146, 'PRNr', 164040), (159, 'AON', 240263), (162, 'LDT', 7472), (173, 'RCH', 16634), (194, 'LHA', 218492), (197, 'RL', 710), (203, 'LIN', 1644), (206, 'RM', 8135), (207, 'AP', 2893), (215, 'APN', 73669), (689, 'VLPO', 1968), (222, 'RO', 969), (223, 'ARH', 15693), (226, 'LPO', 65423), (230, 'RPA', 2242), (235, 'LRN', 31634), (238, 'RPO', 2777), (262, 'RT', 87419), (263, 'AVP', 4181), (272, 'AVPV', 4375), (280, 'B', 2272), (286, 'SCH', 6273), (287, 'BAC', 478), (292, 'BA', 4769), (298, 'MA', 21620), (303, 'BLAa', 46578), (307, 'MARN', 41122), (310, 'SF', 16876), (311, 'BLAp', 40334), (318, 'SG', 1128), (332, 'ASO', 100), (333, 'SH', 2365), (338, 'SFO', 1485), (342, 'SI', 178368), (347, 'SBPV', 9556), (350, 'SLC', 1954), (351, 'BST', 95464), (356, 'PST', 912), (358, 'SLD', 4961), (364, 'PSTN', 3901), (372, 'ICB', 624), (380, 'cuf', 1592), (388, 'grf', 190), (390, 'SO', 3713), (403, 'MEA', 100225), (429, 'SPVC', 84534), (437, 'SPVI', 75635), (445, 'SPVO', 32813), (451, 'BLAv', 16609), (452, 'MEPO', 10218), (470, 'STN', 11971), (507, 'MOB', 1028777), (515, 'MPN', 13213), (523, 'MPO', 53336), (531, 'MPT', 7298), (534, 'SUT', 10567), (536, 'CEA', 78676), (566, 'TR', 66078), (574, 'TRN', 45575), (576, 'ACVII', 350), (581, 'TRS', 13836), (587, 'ND', 1935), (589, 'TT', 98588), (591, 'CLI', 5341), (604, 'NI', 8536), (609, 'SPA', 3375), (443, 'dhc', 60607), (612, 'NLL', 41788), (614, 'TU', 30356), (621, 'V', 24252), (628, 'NOT', 5492), (634, 'NPC', 15510), (642, 'NTB', 6415), (651, 'NTS', 31890), (653, 'VI', 1290), (661, 'VII', 51799), (672, 'CP', 1627409), (679, 'CS', 43388), (693, 'VMH', 29872), (706, 'OP', 5117), (754, 'OT', 242728), (763, 'OV', 755), (765, 'x', 3979), (773, 'XII', 22521), (781, 'y', 1672), (788, 'PAA', 45973), (797, 'ZI', 88068), (814, 'DP', 15028), (830, 'DMH', 31008), (839, 'DMX', 5665), (842, 'SCsg', 97230), (851, 'SCop', 41274), (852, 'PARN', 118338), (859, 'PAS', 860), (867, 'PB', 55843), (872, 'DR', 8056), (880, 'DTN', 5586), (895, 'ECT', 101), (898, 'PCG', 38868), (903, 'ECU', 5822), (914, 'PD', 188), (931, 'PG', 60851), (946, 'PH', 47189), (952, 'EPd', 116160), (961, 'PIR', 822843), (966, 'EPv', 29142), (980, 'PMd', 5281), (998, 'FS', 23171), (1004, 'PMv', 5630), (1022, 'GPe', 103116), (1031, 'GPi', 27122), (1044, 'PP', 2799), (1048, 'GRN', 121816), (1061, 'PPT', 4278), (1069, 'PPY', 2505), (1093, 'PRNc', 198025), (1105, 'IA', 12448), (1109, 'PS', 1720), (449, 'vhc', 7495), (147, 'LC', 1042), (436, 'fx', 17309)]
================================
LEVEL 8
Regions with no associated voxels: 184 set(['DMHa', 'ACA2/3', 'ORB1', 'PVHlp', 'ORB5', 'LRNm', 'VMHdm', 'NTSm', 'NTSl', 'DMHp', 'VISpm', 'VISpl', 'PIR1-3', 'NLOT2', 'SPVOvl', 'NTSco', 'PVHmpv', 'BSTa', 'PIN', 'NTSce', 'MPNc', 'CUL4, 5gr', 'MOBopl', 'BSTp', 'ILA2', 'CUL4pu', 'DP6a', 'PVHm', 'ANcr2gr', 'ANcr2pu', 'PVHf', 'CUL4gr', 'CENT2pu', 'ORBv', 'ORBm', 'ORBl', 'SS2/3', 'PVHp', 'NLOT1-3', 'PAA1-3', 'ANcr1gr', 'CENT2gr', 'SPVOmdmv', 'DMHv', 'CUL4, 5mo', 'NIS', 'VIS2/3', 'PAA2', 'PAA3', 'PAA1', 'PTLp1', 'ORB6a', 'ORB6b', 'SS5', 'SS4', 'MO1', 'VISal', 'VISam', 'MO5', 'PTLp2/3', 'VIS4', 'AONl', 'AONm', 'MOs', 'AONd', 'AONe', 'DP2/3', 'ORBvl', 'TR1', 'TR3', 'TR2', 'CUL5gr', 'ACA1', 'AON2', 'AON1', 'SSp', 'SSs', 'SPVOmdmd', 'ENT', 'SPVOcdm', 'PTLp4', 'CNspg', 'RSPv', 'PIR3', 'PIR2', 'PIR1', 'OT1-3', 'CSm', 'CSl', 'SPVOrdm', 'MO2/3', 'isl', 'PBm', 'PBl', 'CUL5pu', 'PTLp5', 'PVHdp', 'PL2', 'ANcr1pu', 'VMHc', 'VMHa', 'PPYd', 'PPYs', 'CUL4mo', 'AUDv', 'ACAd', 'AUDp', 'CUL5mo', 'NC', 'MOBgr', 'AUDd', 'RSPd', 'ACAv', 'MOBgl', 'CENT2mo', 'NTSge', 'CENT3gr', 'CNlam', 'CENT3pu', 'MOBipl', 'CEAm', 'CEAl', 'CEAc', 'SS1', 'CA', 'MO6b', 'VMHvl', 'TTv', 'VISpor', 'AONpv', 'KF', 'AHNa', 'AHNc', 'AHNd', 'ACA5', 'PTLp6a', 'PTLp6b', 'AHNp', 'VISa', 'VISrl', 'DG', 'AUDpo', 'OT1', 'OT2', 'OT3', 'VISp', 'MOBmi', 'A13', 'AIv', 'AIp', 'MEApd', 'LRNp', 'TM', 'islm', 'CENT3mo', 'TTd', 'MEApv', 'VISl', 'VISli', 'VIS6a', 'VIS6b', 'NLLv', 'ANcr1mo', 'TR1-3', 'MEAav', 'VP', 'NLLh', 'NLLd', 'MO6a', 'MEAad', 'DP5', 'VIS1', 'DP2', 'DP1', 'VIS5', 'MPNm', 'MPNl', 'ORB2/3', 'ANcr2mo', 'ACA6b', 'CUL4, 5pu', 'ACA6a', 'SS6a', 'SS6b'])

Regions with at least one voxel associated: 120
[(15, 'PT', 17129), (19, 'IG', 1772), (27, 'IGL', 4556), (36, 'GU1', 28255), (59, 'IMD', 9498), (64, 'AD', 10074), (68, 'FRP1', 41158), (84, 'PL6a', 16388), (96, 'DCO', 37734), (97, 'TEa1', 28872), (101, 'VCO', 48367), (104, 'AId', 8), (105, 'SOCm', 9435), (114, 'SOCl', 12408), (122, 'POR', 8213), (127, 'AM', 21509), (132, 'PL6b', 1333), (143, 'AMBv', 1406), (148, 'GU4', 12954), (149, 'PVT', 25631), (155, 'LD', 68148), (169, 'PRP', 10012), (170, 'LGd', 46211), (171, 'PL1', 30765), (177, 'NR', 1389), (178, 'LGv', 25851), (180, 'GU2/3', 43435), (181, 'RE', 24576), (186, 'LH', 20631), (187, 'GU5', 52850), (188, 'AOBgl', 9454), (189, 'RH', 4099), (196, 'AOBgr', 13589), (202, 'MV', 86910), (204, 'AOBmi', 17753), (209, 'LAV', 9057), (210, 'LM', 2062), (218, 'LP', 69350), (225, 'SPIV', 39468), (234, 'TEa4', 14555), (250, 'LSc', 34325), (255, 'AV', 26109), (258, 'LSr', 138535), (260, 'NLOT1', 10457), (266, 'LSv', 7287), (289, 'TEa5', 58112), (304, 'PL2/3', 34208), (321, 'SubG', 1726), (325, 'SGN', 2009), (335, 'PERI6a', 16967), (362, 'MD', 86351), (363, 'PL5', 28852), (366, 'SMT', 11317), (368, 'PERI6b', 5078), (414, 'SPFm', 5070), (422, 'SPFp', 16827), (427, 'ECT2/3', 40003), (475, 'MG', 30179), (483, 'MH', 19039), (491, 'MM', 36227), (494, 'SCig-a', 8243), (502, 'SUB', 289037), (503, 'SCig-b', 26045), (511, 'SCig-c', 4918), (525, 'SUM', 12409), (540, 'PERI1', 29487), (556, 'ILA2/3', 33218), (564, 'MS', 39040), (575, 'CL', 15427), (596, 'NDB', 67827), (599, 'CM', 10754), (638, 'GU6a', 47156), (639, 'COAa', 36755), (647, 'COAp', 124500), (662, 'GU6b', 4219), (667, 'FRP2/3', 7768), (685, 'VM', 47254), (692, 'PERI5', 41974), (707, 'ILA1', 17329), (711, 'CU', 19078), (729, 'TEa6a', 30318), (786, 'TEa6b', 8572), (804, 'FF', 3662), (827, 'ILA5', 29789), (836, 'ECT1', 25197), (843, 'PAR', 59900), (849, 'VISC6b', 2660), (857, 'VISC6a', 30680), (888, 'PERI2/3', 46144), (894, 'RSPagl', 3439), (897, 'VISC1', 18373), (907, 'PCN', 16836), (930, 'PF', 26914), (939, 'AMBd', 1210), (970, 'PGRNd', 9881), (977, 'ECT6a', 19220), (978, 'PGRNl', 42251), (982, 'FC', 4025), (985, 'MOp', 34027), (988, 'ECT5', 40220), (1010, 'VISC4', 12560), (1020, 'PO', 75867), (1029, 'POL', 6759), (1037, 'POST', 53374), (1039, 'GR', 4320), (1045, 'ECT6b', 6388), (1054, 'ILA6a', 22853), (1058, 'VISC5', 30944), (1077, 'PR', 1164), (1081, 'ILA6b', 6116), (1084, 'PRE', 67752), (1098, 'MDRNd', 22557), (1106, 'VISC2/3', 26981), (1107, 'MDRNv', 48209), (1113, 'IAD', 5646), (1120, 'IAM', 2981), (1127, 'TEa2/3', 42070), (1139, 'NLOT3', 14544), (217, 'SUV', 22710), (629, 'VAL', 46277)]
================================
LEVEL 9
Regions with no associated voxels: 89 set(['DGlb', 'BSTam', 'BSTal', 'RSPv2', 'SSp-tr', 'BSTif', 'BSTd', 'POST3', 'POST2', 'POST1', 'BSTv', 'PBlc', 'DGmb', 'BSTfu', 'SSp-un', 'PVHmpd', 'SSp-ul', 'TTv2', 'DG-sgz', 'PAR3', 'PAR2', 'PAR1', 'BSTov', 'SSp2/3', 'AMd', 'AMv', 'MEApd-c', 'MEApd-b', 'MEApd-a', 'BSTpr', 'PBlv', 'MGd', 'COApl', 'COApm', 'Mmme', 'TTv1', 'SSp-ll', 'PVHam', 'ORBm2', 'MGm', 'MGv', 'PVHap', 'SUBv', 'TTd1', 'TTd2', 'TTd3', 'TTd4', 'PVHpv', 'SSp6a', 'SUBd', 'SSp6b', 'COAa1', 'COAa3', 'COAa2', 'MDl', 'MDm', 'MDc', 'PBld', 'PBle', 'PBls', 'PRE1', 'TTv3', 'PRE3', 'PRE2', 'BSTju', 'SUMm', 'SUMl', 'LGvm', 'LGvl', 'BSTrh', 'AIp6b', 'BSTse', 'SSp-m', 'SSp-n', 'TTv1-3', 'BSTtr', 'PVHpm', 'PBmv', 'BSTdm', 'PBme', 'PBmm', 'SSp4', 'SSp5', 'SSp1', 'PVHmm', 'BSTmg', 'SSp-bfd', 'DGcr', 'TTd1-4'])

Regions with at least one voxel associated: 169
[(1, 'TMv', 1884), (33, 'VISp6a', 82014), (41, 'VISpm2/3', 19917), (74, 'VISl6a', 12778), (120, 'AIp1', 31189), (121, 'VISl6b', 1610), (156, 'AUDd6a', 52134), (163, 'AIp2/3', 33257), (211, 'ACAd2/3', 50656), (233, 'VISal5', 14082), (243, 'AUDd6b', 10217), (249, 'AUDpo6a', 1307), (251, 'AUDp2/3', 39182), (252, 'AUDd5', 75964), (257, 'VISpm6a', 11585), (269, 'VISpl2/3', 9689), (274, 'RSPd6a', 46322), (279, 'RSPagl6b', 158), (281, 'VISam1', 8083), (288, 'ORBvl2/3', 35411), (296, 'ACAv2/3', 46352), (305, 'VISp6b', 10335), (314, 'AIp6a', 13892), (320, 'MOp1', 149419), (328, 'AId2/3', 70598), (330, 'RSPd6b', 1223), (344, 'AIp5', 29191), (377, 'VISpl6a', 9301), (382, 'CA1', 555065), (393, 'VISpl6b', 126), (401, 'VISam4', 5000), (412, 'ORBl2/3', 44449), (421, 'VISl1', 9779), (423, 'CA2', 88966), (430, 'RSPv2/3', 88080), (433, 'VISam5', 13976), (434, 'RSPd2/3', 64270), (440, 'ORBl6a', 11953), (441, 'VISam6b', 1086), (442, 'RSPd1', 30893), (448, 'ORBl1', 17882), (456, 'AUDpo6b', 270), (463, 'CA3', 354791), (10704, 'DG-po', 51187), (469, 'VISpm6b', 2275), (484, 'ORBm1', 19533), (488, 'ORBl6b', 19), (501, 'VISpm4', 7169), (520, 'AUDv6a', 30840), (527, 'AUDd1', 29577), (542, 'RSPv1', 35900), (545, 'RSPd4', 976), (565, 'VISpm5', 16877), (573, 'VISl4', 7752), (582, 'ORBm2/3', 23368), (588, 'ACAv1', 33890), (590, 'RSPv6a', 59776), (593, 'VISp1', 56148), (598, 'AUDv6b', 4998), (600, 'AUDd2/3', 60142), (601, 'VISal6a', 7715), (608, 'ORBvl6a', 9998), (610, 'RSPd5', 76561), (613, 'VISl5', 21668), (620, 'ORBm5', 15388), (622, 'RSPv6b', 5614), (630, 'ORBl5', 30371), (632, 'DG-sg', 101404), (643, 'AUDpo2/3', 3137), (648, 'MOp5', 344358), (649, 'VISal6b', 1441), (656, 'MOs1', 121223), (671, 'RSPagl1', 19519), (675, 'AIv6a', 6165), (678, 'AUDd4', 25950), (680, 'ORBvl6b', 76), (687, 'RSPv5', 118143), (694, 'AIv2/3', 39670), (696, 'AUDpo1', 1728), (699, 'AIv6b', 8), (704, 'AIv1', 30945), (718, 'VPL', 40066), (721, 'VISp4', 61357), (725, 'VPLpc', 2273), (733, 'VPM', 102689), (735, 'AUDp1', 16294), (741, 'VPMpc', 8633), (750, 'VISpl1', 4462), (755, 'AUDv2/3', 33060), (759, 'AUDpo4', 1356), (767, 'MOs5', 191770), (772, 'ACAv5', 37081), (774, 'RSPagl5', 24606), (778, 'VISp5', 132000), (783, 'AId6a', 19039), (791, 'AUDpo5', 3728), (800, 'AIv5', 25696), (805, 'VISpm1', 11155), (806, 'SSs2/3', 105435), (810, 'ACAv6a', 18507), (816, 'AUDp4', 16563), (819, 'ACAv6b', 4443), (821, 'VISp2/3', 132154), (831, 'AId6b', 418), (844, 'MOp6a', 244017), (847, 'AUDp5', 46593), (862, 'SSs6a', 68575), (869, 'VISpl4', 3003), (873, 'SSs1', 49609), (882, 'MOp6b', 43771), (893, 'SSs6b', 8102), (902, 'VISpl5', 14368), (905, 'VISal2/3', 11661), (906, 'RSPagl6a', 14114), (910, 'ORBm6a', 6524), (918, 'ENTl', 526878), (919, 'ACAd6a', 22529), (926, 'ENTm', 191405), (927, 'ACAd6b', 3005), (934, 'ENTmv', 16552), (935, 'ACAd1', 34249), (943, 'MOp2/3', 330988), (954, 'AUDp6a', 21820), (959, 'AUDv1', 18862), (962, 'MOs2/3', 208905), (965, 'RSPagl2/3', 24510), (969, 'ORBvl1', 25156), (973, 'VISl23', 20753), (990, 'AUDv4', 16821), (996, 'AId1', 40001), (1005, 'AUDp6b', 6090), (1015, 'ACAd5', 43701), (1021, 'MOs6a', 79750), (1023, 'AUDv5', 47654), (1035, 'SSs4', 50294), (1046, 'VISam6a', 9415), (1066, 'VISam2/3', 13266), (1074, 'VISal1', 5504), (1085, 'MOs6b', 5942), (1090, 'SSs5', 93754), (1101, 'AId5', 41939), (1114, 'VISal4', 4975), (1125, 'ORBvl5', 15049), (1126, 'TMd', 936), (312782550, 'VISa1', 5549), (312782554, 'VISa2/3', 10105), (312782558, 'VISa4', 4552), (312782562, 'VISa5', 10867), (312782566, 'VISa6a', 6155), (312782570, 'VISa6b', 899), (312782578, 'VISli1', 4320), (312782582, 'VISli2/3', 7962), (312782586, 'VISli4', 2933), (312782590, 'VISli5', 9603), (312782594, 'VISli6a', 5424), (312782598, 'VISli6b', 892), (312782604, 'VISrl1', 7720), (312782608, 'VISrl2/3', 16931), (312782612, 'VISrl4', 7218), (312782616, 'VISrl5', 19519), (312782620, 'VISrl6a', 10905), (312782624, 'VISrl6b', 2699), (312782632, 'VISpor1', 9528), (312782636, 'VISpor2/3', 24180), (312782640, 'VISpor4', 2354), (312782644, 'VISpor5', 23387), (312782648, 'VISpor6a', 8268), (312782652, 'VISpor6b', 979), (10703, 'DG-mo', 257103)]
================================
LEVEL 10
Regions with no associated voxels: 66 set(['ENTl6a', 'DGcr-po', 'ENTl6b', 'COApl1-3', 'COApl1-2', 'ENTl4', 'ENTl5', 'ENTl2', 'ENTl3', 'DGlb-mo', 'SUBd-sp', 'DGmb-mo', 'DGcr-mo', 'COApm1', 'COApm2', 'COApm3', 'ENTmv4', 'SUBd-sr', 'ENTmv1', 'ENTmv2', 'ENTmv3', 'CA3sr', 'CA3sp', 'SUBd-m', 'SUBv-m', 'ENTm5/6', 'CA2sp', 'CA1slm', 'ENTm2a', 'ENTm2b', 'CA3so', 'DGmb-sg', 'ENTl2b', 'DGcr-sg', 'CA3slm', 'COApm1-2', 'COApm1-3', 'ENTm5', 'ENTm4', 'CA2slm', 'ENTm6', 'ENTm1', 'CA2so', 'PVHpmm', 'ENTm2', 'COApl3', 'COApl2', 'COApl1', 'ENTm3', 'PVHpml', 'CA3slu', 'DGlb-sg', 'CA2sr', 'SUBv-sr', 'SUBv-sp', 'ENTl2/3', 'ENTl2a', 'DGlb-po', 'DGmb-po', 'CA1sp', 'CA1sr', 'ENTmv5/6', 'CA1so', 'ENTl4/5', 'ENTl1', 'ENTl5/6'])

Regions with at least one voxel associated: 42
[(2, 'SSp-m6b', 7487), (9, 'SSp-tr6a', 10674), (113, 'SSp-ll2/3', 33461), (201, 'SSp-bfd2/3', 112563), (450, 'SSp-ul1', 20915), (461, 'SSp-tr6b', 2045), (478, 'SSp-ll6a', 24511), (510, 'SSp-ll6b', 6193), (558, 'SSp-n1', 21892), (577, 'SSp-ul4', 30178), (625, 'SSp-ul5', 61493), (654, 'SSp-n4', 32103), (657, 'SSp-m2/3', 126863), (670, 'SSp-tr2/3', 22678), (702, 'SSp-n5', 55951), (838, 'SSp-n2/3', 63473), (854, 'SSp-ul2/3', 60653), (878, 'SSp-m1', 50816), (889, 'SSp-n6a', 52755), (182305693, 'SSp-un1', 8938), (929, 'SSp-n6b', 5198), (182305701, 'SSp-un4', 12962), (182305705, 'SSp-un5', 27361), (182305709, 'SSp-un6a', 21825), (182305713, 'SSp-un6b', 4331), (950, 'SSp-m4', 57011), (974, 'SSp-m5', 103521), (981, 'SSp-bfd1', 35823), (1006, 'SSp-tr1', 8145), (1026, 'SSp-ul6b', 10803), (1030, 'SSp-ll1', 13610), (1038, 'SSp-bfd6a', 78454), (1047, 'SSp-bfd4', 51691), (1062, 'SSp-bfd6b', 14293), (1070, 'SSp-bfd5', 103351), (1086, 'SSp-tr4', 6378), (1094, 'SSp-ll4', 14951), (1102, 'SSp-m6a', 79118), (1111, 'SSp-tr5', 25067), (1128, 'SSp-ll5', 33760), (182305697, 'SSp-un2/3', 29066), (945, 'SSp-ul6a', 46906)]
Total no of regions with associated voxels: 588 total not associated 646

In [11]:
#Find conflicting mappings
for m in mapped_regions:
    my_children = childs_map[m]
    children_mapped = set(my_children).intersection(mapped_regions)
    children_mapped = [conn2voxels[c] for c in children_mapped]
    if children_mapped:
        print "================================"
        print "Conflict between parent associations for node:", m, conn2voxels[m]
        if len(children_mapped) > 20:
            print "and ", len(children_mapped), "child nodes (all with associated voxels as well)"
        else:
            print "and children with voxels", children_mapped


================================
Conflict between parent associations for node: 997 (1, 'root', 40)
and  587 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 73 (2, 'VS', 18385)
and children with voxels [(3, 'V3', 82850), (3, 'c', 80), (3, 'V4', 50656), (3, 'AQ', 28584), (3, 'VL', 159548), (4, 'V4r', 53228)]
================================
Conflict between parent associations for node: 1009 (2, 'fiber tracts', 236648)
and  56 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 512 (3, 'CB', 1)
and children with voxels [(6, 'CUL', 448510), (6, 'PRM', 315291), (6, 'FOTU', 58946), (6, 'SIM', 351345), (6, 'DEC', 239591), (6, 'COPY', 155348), (5, 'DN', 21468), (6, 'UVU', 140142), (6, 'LING', 7942), (6, 'PFL', 356892), (6, 'AN', 733048), (6, 'NOD', 80661), (6, 'PYR', 58832), (6, 'CENT', 250358), (6, 'FL', 84339), (5, 'IP', 55207), (5, 'FN', 32651)]
================================
Conflict between parent associations for node: 145 (3, 'V4', 50656)
and children with voxels [(4, 'V4r', 53228)]
================================
Conflict between parent associations for node: 784 (4, 'cst', 4959)
and children with voxels [(5, 'cpd', 55637), (5, 'py', 34577), (5, 'int', 155764)]
================================
Conflict between parent associations for node: 877 (4, 'tsp', 10009)
and children with voxels [(5, 'dtd', 2632)]
================================
Conflict between parent associations for node: 848 (4, 'IIn', 3786)
and children with voxels [(5, 'och', 18309), (5, 'opt', 43331)]
================================
Conflict between parent associations for node: 313 (4, 'MB', 353356)
and  43 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 863 (4, 'rust', 22757)
and children with voxels [(5, 'vtd', 4164)]
================================
Conflict between parent associations for node: 771 (5, 'P', 384669)
and  23 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 803 (5, 'PAL', 111434)
and children with voxels [(7, 'TRS', 13836), (7, 'GPi', 27122), (7, 'MA', 21620), (8, 'MS', 39040), (7, 'SI', 178368), (8, 'NDB', 67827), (7, 'BAC', 478), (7, 'GPe', 103116), (7, 'BST', 95464)]
================================
Conflict between parent associations for node: 549 (5, 'TH', 295591)
and  40 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 703 (5, 'CTXsp', 105638)
and children with voxels [(6, 'LA', 41076), (7, 'EPv', 29142), (6, 'CLA', 31625), (6, 'PA', 20241), (7, 'BLAa', 46578), (7, 'BLAv', 16609), (7, 'BLAp', 40334), (7, 'EPd', 116160), (6, 'BMA', 69327)]
================================
Conflict between parent associations for node: 1097 (5, 'HY', 232031)
and  42 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 477 (5, 'STR', 175862)
and children with voxels [(7, 'CP', 1627409), (8, 'LSr', 138535), (7, 'BA', 4769), (7, 'SF', 16876), (7, 'FS', 23171), (7, 'IA', 12448), (8, 'LSv', 7287), (7, 'SH', 2365), (7, 'ACB', 263008), (7, 'OT', 242728), (7, 'MEA', 100225), (6, 'sAMY', 15406), (7, 'AAA', 37040), (7, 'CEA', 78676), (8, 'LSc', 34325)]
================================
Conflict between parent associations for node: 354 (5, 'MY', 575501)
and  44 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 278 (6, 'sAMY', 15406)
and children with voxels [(7, 'CEA', 78676), (7, 'IA', 12448), (7, 'MEA', 100225), (7, 'BA', 4769), (7, 'AAA', 37040)]
================================
Conflict between parent associations for node: 795 (6, 'PAG', 283394)
and children with voxels [(7, 'ND', 1935), (7, 'PRC', 5291), (7, 'INC', 3513)]
================================
Conflict between parent associations for node: 698 (6, 'OLF', 275171)
and children with voxels [(7, 'PIR', 822843), (8, 'AOBgr', 13589), (8, 'COAp', 124500), (8, 'COAa', 36755), (8, 'AOBmi', 17753), (7, 'TT', 98588), (7, 'DP', 15028), (8, 'NLOT3', 14544), (7, 'PAA', 45973), (7, 'TR', 66078), (8, 'NLOT1', 10457), (7, 'MOB', 1028777), (8, 'AOBgl', 9454), (7, 'AON', 240263)]
================================
Conflict between parent associations for node: 315 (6, 'Isocortex', 6)
and  240 child nodes (all with associated voxels as well)
================================
Conflict between parent associations for node: 1089 (6, 'HPF', 171733)
and children with voxels [(8, 'FC', 4025), (9, 'ENTmv', 16552), (9, 'CA2', 88966), (8, 'PAR', 59900), (9, 'ENTl', 526878), (8, 'POST', 53374), (9, 'DG-mo', 257103), (9, 'DG-po', 51187), (8, 'IG', 1772), (9, 'ENTm', 191405), (8, 'SUB', 289037), (9, 'DG-sg', 101404), (9, 'CA3', 354791), (8, 'PRE', 67752), (9, 'CA1', 555065)]
================================
Conflict between parent associations for node: 10 (7, 'SCig', 133547)
and children with voxels [(8, 'SCig-c', 4918), (8, 'SCig-a', 8243), (8, 'SCig-b', 26045)]
================================
Conflict between parent associations for node: 797 (7, 'ZI', 88068)
and children with voxels [(8, 'FF', 3662)]
================================
Conflict between parent associations for node: 895 (7, 'ECT', 101)
and children with voxels [(8, 'ECT6a', 19220), (8, 'ECT5', 40220), (8, 'ECT2/3', 40003), (8, 'ECT1', 25197), (8, 'ECT6b', 6388)]
================================
Conflict between parent associations for node: 104 (8, 'AId', 8)
and children with voxels [(9, 'AId2/3', 70598), (9, 'AId6b', 418), (9, 'AId1', 40001), (9, 'AId5', 41939), (9, 'AId6a', 19039)]
================================
Conflict between parent associations for node: 894 (8, 'RSPagl', 3439)
and children with voxels [(9, 'RSPagl6a', 14114), (9, 'RSPagl6b', 158), (9, 'RSPagl2/3', 24510), (9, 'RSPagl5', 24606), (9, 'RSPagl1', 19519)]
================================
Conflict between parent associations for node: 985 (8, 'MOp', 34027)
and children with voxels [(9, 'MOp5', 344358), (9, 'MOp1', 149419), (9, 'MOp6b', 43771), (9, 'MOp6a', 244017), (9, 'MOp2/3', 330988)]

In [ ]: