Inventario imagenes kmc400 en colivri

MPRAGE

Basado en la lista enviada por pablo



In [1]:

    
asistieron_husi={9,9,15,19,25,29,31,35,44,51,54,56,64,65,69,71,72,73,75,83,90,93,95,107,108,113,119,121,123,124,125,128,129,134,138,141,143,144,145,149,151,153,154,156,157,165,173,175,176,177,182,185,186,195,197,198,201,202,205,208,216,219,221,225,227,230,231,232,235,237,253,256,259,261,263,264,266,266,277,292,293,300,301,307,310,310,312,313,314,319,320,322,327,328,331,332,333,344,346,348,353,355,356,356,357,358,364,369,371,374,381,390,396,399,409,413,416,417,423,424,426,427,429,431,432,440,447,452,456,458,464,469,472,478,480,483,484,485,491,496,499,504,517,526,532,535,536,537,542,544,545,547,548,549,552,566,576,579,580,592,593,595,599,600,602,610,611,615,616,619,623,625,630,631,645,651,662,665,670,675,684,686,689,694,696,712,715,734,754,761,765,769,783,784,786,789,790,791,804,806,815,818,821,829,840,841,848,850,861,863,868,869,874,876,877,878,879,884,891,892,893,894,898,905,906,912,913,918,928,934,935,939,940,942,953,954,965,966,971,982,984,992,994,1005,1006,1021,1026,1039,1049,1077,1212,1218,1221,1232,1237,1242,1244,1253,1260,1265,1267,1320,1322,1326,1333,1338,1340,1357,}
len(asistieron_husi)









    Out[1]:





264



In [2]:

    
mp_rage_colivri={292 ,153 ,593 ,874 ,353 ,552 ,322 ,331 ,333 ,154 ,483 ,348 ,141 ,429 ,761 ,966 ,216 ,176 ,548 ,54 ,884 ,143 ,1338 ,156 ,1021 ,431 ,611 ,850 ,458 ,327 ,29 ,1005 ,623 ,128 ,547 ,861 ,313 ,332 ,399 ,1326 ,129  ,452 ,151 ,307 ,670 ,256 ,356 ,396 ,992 ,464 ,754 ,545 ,675 ,230 ,478 ,123 ,1049 ,235 ,686 ,65 ,237 ,526 ,965 ,804 ,266 ,734 ,427 ,894 ,898 ,205 ,549 ,765 ,71 ,201 ,201 ,566 ,182 ,730 ,535 ,928 ,868 ,121 ,355 ,253 ,1221 ,31 ,818 ,472 ,35 ,301 ,277 ,381 ,829 ,954 ,416 ,344 ,942 ,878 ,197 ,863 ,409 ,19 ,783 ,1026 ,177 ,544 ,1212 ,918 ,293 ,310 ,195 ,579 ,432 ,417 ,537 ,469 ,20 ,95 ,95 ,95 ,689 ,119 ,789 ,769 ,1253 ,610 ,559 ,484 ,145 ,599 ,953 ,1333 ,982 ,592 ,231 ,202 ,532 ,1242 ,25 ,263 ,125 ,712 ,504 ,456 ,456 ,491 ,499 ,784 ,72 ,984 ,1006 ,320 ,869 ,906 ,157 ,165 ,364 ,64 ,1265 ,113 ,631 ,536 ,1357 ,595 ,616 ,879 ,715 ,651 ,198 ,225 ,500 ,56 ,630 ,1260 ,892 ,371 ,221 ,893 ,1232 ,264 ,357 ,219 ,791 ,600 ,346 ,44 ,940 ,1320 ,369 ,423 ,934 ,480 ,841 ,51 ,645 ,580 ,440 ,319 ,1039 ,542 ,876 ,905 ,426 ,314 ,15 ,227 ,815 ,625 ,877 ,877 ,891 ,9 ,576 ,994 ,138 ,186 ,517 ,790 ,358 ,107 ,684 ,175 ,390 ,173 ,665 ,939 ,1077 ,124 ,840 ,185 ,806 ,485 ,413 ,619 ,662 ,108 ,300 ,935 ,696 ,964 ,971 ,912 ,69 ,232 ,786 ,821 ,1340 ,144 ,848 ,602 ,83 ,}
len(mp_rage_colivri)









    Out[2]:





246



In [3]:

    
nii_in_colivri={1005, 1006, 1021, 1026, 1039, 1049, 107, 1077, 108, 113, 119, 121, 1212, 1221, 123, 1232, 124, 1242, 125, 1253, 1260, 1265, 128, 129, 1320, 1326, 1333, 1338, 1340, 1357, 138, 141, 143, 144, 145, 15, 151, 153, 154, 156, 157, 165, 173, 175, 176, 177, 182, 185, 186, 19, 195, 197, 198, 20, 201, 202, 205, 216, 219, 221, 225, 227, 230, 231, 232, 235, 237, 25, 253, 256, 263, 264, 266, 277, 29, 292, 293, 300, 301, 307, 31, 310, 312, 313, 314, 319, 320, 322, 327, 331, 332, 333, 344, 346, 348, 35, 353, 355, 356, 357, 358, 364, 369, 371, 374, 381, 390, 396, 399, 409, 413, 416, 417, 423, 426, 427, 429, 431, 432, 44, 440, 452, 456, 458, 464, 469, 472, 478, 480, 483, 484, 485, 491, 499, 500, 504, 51, 517, 526, 532, 535, 536, 537, 54, 542, 544, 545, 547, 548, 549, 552, 559, 56, 566, 576, 579, 580, 592, 593, 595, 599, 600, 602, 610, 611, 616, 619, 623, 625, 630, 631, 64, 645, 65, 651, 662, 665, 670, 675, 684, 686, 689, 69, 696, 71, 712, 715, 72, 730, 734, 754, 761, 765, 769, 783, 784, 786, 789, 790, 791, 804, 806, 815, 818, 821, 829, 83, 840, 841, 848, 850, 861, 863, 868, 869, 874, 876, 877, 878, 879, 884, 891, 892, 893, 894, 898, 9, 905, 906, 912, 918, 928, 934, 935, 939, 940, 942, 95, 953, 954, 964, 965, 966, 971, 982, 984, 992, 994, }
print len(nii_in_colivri)



In [4]:

    
sorted(nii_in_colivri-mp_rage_colivri)









    Out[4]:





[312, 374]



In [5]:

    
missing_mp_rage_in_colivri = asistieron_husi - mp_rage_colivri
print sorted(missing_mp_rage_in_colivri)









    



[73, 75, 90, 93, 134, 149, 208, 259, 261, 312, 328, 374, 424, 447, 496, 615, 694, 913, 1218, 1237, 1244, 1267, 1322]



In [6]:

    
missing_in_colivri = asistieron_husi - nii_in_colivri
print sorted(missing_in_colivri)









    



[73, 75, 90, 93, 134, 149, 208, 259, 261, 328, 424, 447, 496, 615, 694, 913, 1218, 1237, 1244, 1267, 1322]



In [7]:

    
pilotos={20,2,3,730}

Control de Calidad



In [8]:

    
moved_francoises = {31, 35, 72, 95, 123, 129, 399, 456, 547, 850, 994}
print len(moved_francoises)



In [9]:

    
uriza_star={125, 143, 195, 201, 219, 227, 230, 232, 235, 256, 263, 277, 292, 307, 331, 332, 357, 364, 390, 413, 416, 472, 478, 499, 500, 532, 595, 602, 684, 784, 790, 906, 984}
print len(uriza_star)



In [10]:

    
uriza_wrong ={31, 72, 95, 123, 129, 202, 300, 319, 322, 348, 355, 358, 381, 399, 456, 458, 464, 547, 645, 696, 821, 850, 898, 994, 994}
print len(uriza_wrong)



In [11]:

    
pablo_wrong={31, 95, 123, 186, 322, 348, 355, 399, 464, 499, 547, 611, 786, 898, 994, 662, 396, 532, 264, 821, }
len(pablo_wrong)









    Out[11]:





20



In [12]:

    
print len(moved_francoises.intersection(uriza_wrong))



In [13]:

    
print len(moved_francoises.intersection(pablo_wrong))



In [14]:

    
print len(pablo_wrong.intersection(uriza_wrong))



In [15]:

    
any_wrong = pablo_wrong.union(uriza_wrong).union(moved_francoises)
print len(any_wrong)



In [16]:

    
saved=mp_rage_colivri - any_wrong-pilotos
print len(saved)



In [65]:

    
print str(saved)









    



set([1026, 427, 517, 9, 526, 15, 19, 535, 536, 537, 29, 542, 544, 545, 548, 549, 177, 552, 44, 559, 51, 1077, 566, 56, 576, 65, 579, 580, 69, 71, 314, 592, 593, 83, 599, 600, 356, 602, 1039, 869, 610, 616, 107, 108, 623, 625, 630, 631, 121, 124, 125, 277, 128, 619, 364, 138, 651, 141, 143, 144, 145, 25, 665, 154, 156, 157, 670, 1221, 675, 165, 113, 684, 173, 686, 175, 176, 689, 182, 185, 1212, 195, 197, 198, 712, 201, 715, 205, 1232, 216, 292, 1242, 219, 221, 734, 225, 227, 1253, 230, 231, 232, 235, 1260, 237, 1265, 754, 761, 253, 256, 769, 263, 266, 452, 783, 784, 472, 1049, 789, 790, 791, 804, 293, 806, 1320, 301, 1326, 815, 818, 307, 1333, 310, 313, 1338, 1340, 829, 320, 992, 54, 327, 840, 841, 331, 332, 1357, 848, 344, 346, 861, 863, 353, 868, 357, 874, 876, 877, 878, 879, 369, 371, 884, 891, 892, 893, 894, 64, 491, 390, 905, 906, 151, 912, 918, 153, 409, 413, 416, 417, 934, 423, 426, 939, 940, 429, 942, 431, 432, 440, 953, 954, 928, 964, 965, 966, 971, 333, 469, 982, 984, 478, 480, 483, 484, 485, 935, 1005, 1006, 765, 595, 500, 504, 119, 1021])



In [17]:

    
print "\n".join(map(str,sorted(saved)))

copy info to google drive



In [19]:

    
rows = index.split()



In [20]:

    
def has_mprage(s):
    try:
        s_i = int(s)
        if s_i in mp_rage_colivri:
            return "1"
    except Exception:
        pass
    return "0"
col_uniandes=[has_mprage(s) for s in rows]



In [21]:

    
print "\n".join(col_uniandes)



In [22]:

    
730 in mp_rage_colivri









    Out[22]:





True

Tabla mprage



In [23]:

    
index ="""9
9
15
19
31
35
44
51
56
56
64
65
69
72
72
73
83
95
107
108
119
121
123
124
125
128
129
138
141
143
144
145
151
153
156
157
165
173
175
176
177
182
185
186
195
198
201
202
205
216
219
221
225
227
230
231
232
235
237
256
259
263
292
293
300
307
310
310
320
322
327
331
332
333
344
346
348
355
356
357
358
364
369
371
381
390
399
409
413
416
417
423
426
427
429
440
447
452
456
458
464
469
472
478
480
483
485
491
499
500
504
517
526
535
539
542
544
545
547
548
549
559
566
576
579
592
593
595
600
610
611
616
619
623
625
630
631
645
645
651
665
675
684
686
696
712
730
734
754
761
765
769
783
784
786
791
815
829
848
850
863
868
869
874
878
879
891
892
893
898
928
935
939
942
953
954
964
971
971
984
992
994
1005
1021
1039
1049
1077
20pil
2pil
3pil
54
789
790
25
253
301
877
934
662
905
266
552
113
982
599
670
396
353
580
1026
532
966
940
374
71
314
1265
602
1232
312
861
1221
840
536
537
154
277
818
1320
689
1333
1006
432
431
484
1253
918
313
841
1267
319
804
356
1326
264
906
806
821
1338
894
1242
1212
197
884
965
876"""



In [24]:

    
subjs=index.split()



In [25]:

    
def get_approved(s):
    try:
        s_i = int(s)
        if s_i in saved:
            return "1"
    except ValueError:
        pass
    return "0"



In [26]:

    
approved_col=[get_approved(s) for s in subjs]



In [27]:

    
print "\n".join(approved_col)



In [28]:

    
i_subjs = [int(s) for s in subjs if s.isdigit()]



In [29]:

    
missing_saved = saved - set(i_subjs)



In [30]:

    
missing_bad = any_wrong - set(i_subjs)



In [31]:

    
print missing_saved









    



set([715, 1260, 1357, 912, 1340, 29])



In [32]:

    
print missing_bad









    



set([])



In [33]:

    
print "\n".join(map(str,missing_saved))



In [34]:

    
def get_star(s):
    try:
        s_i = int(s)
        if s_i in uriza_star:
            return "1"
    except ValueError:
        pass
    return "0"



In [35]:

    
starred = [get_star(s) for s in subjs]
print "\n".join(map(str,starred))



In [36]:

    
uriza_star - set(i_subjs)









    Out[36]:





set()

FreeSurfer Completed succesfully

We will check the existance of a aparc+aseg image



In [37]:

    
all_aparc={292, 153, 593, 874, 353, 552, 322, 331, 333, 154, 483, 348, 141, 429, 761, 966, 216, 176, 548, 54, 884, 143, 1338, 156, 1021, 431, 611, 850, 458, 327, 29, 1005, 623, 128, 547, 861, 313, 332, 399, 1326, 129, 452, 151, 307, 670, 256, 396, 992, 464, 754, 545, 675, 230, 478, 123, 1049, 235, 686, 65, 237, 526, 965, 804, 266, 734, 427, 894, 898, 205, 549, 765, 71, 201, 566, 182, 730, 535, 928, 868, 121, 355, 253, 1221, 31, 818, 472, 35, 301, 277, 381, 73, 829, 954, 416, 344, 942, 878, 197, 863, 409, 19, 783, 1026, 177, 544, 1212, 374, 918, 293, 310, 195, 579, 432, 417, 537, 469, 20, 95, 689, 119, 789, 769, 1253, 610, 559, 484, 145, 599, 953, 1333, 982, 592, 231, 202, 532, 1242, 25, 263, 125, 712, 504, 456, 491, 499, 784, 72, 984, 1006, 320, 869, 906, 157, 165, 364, 64, 1265, 113, 631, 536, 1357, 595, 616, 879, 715, 651, 198, 225, 500, 56, 630, 1260, 892, 371, 221, 893, 1232, 264, 357, 219, 791, 600, 346, 44, 940, 1320, 369, 423, 934, 480, 841, 51, 645, 580, 440, 319, 1039, 542, 876, 905, 426, 314, 15, 227, 815, 625, 877, 891, 9, 576, 994, 138, 186, 517, 790, 358, 107, 684, 175, 390, 173, 665, 939, 1077, 124, 840, 185, 806, 485, 413, 619, 662, 108, 300, 935, 696, 964, 971, 912, 69, 232, 786, 821, 1340, 144, 848, 602, 83}
print len(all_aparc)



In [38]:

    
print saved - all_aparc









    



set([356])

IMPORTANTE, TOCA PROCESARLO

DTI



In [39]:

    
dti_colivri={ 292, 153, 593, 874, 353, 552, 322, 331, 333, 154, 483, 348, 141, 429, 761, 966, 216, 176, 548, 54, 884, 143, 1338, 156, 1021, 431, 611, 850, 458, 327, 29, 1005, 623, 128, 547, 861, 313, 332, 399, 1326, 129, 452, 151, 307, 670, 256, 356, 396, 992, 464, 754, 545, 675, 230, 478, 123, 1049, 235, 686, 65, 237, 526, 804, 266, 734, 427, 894, 898, 205, 549, 765, 71, 201, 566, 182, 535, 928, 868, 121, 355, 253, 1221, 31, 818, 472, 35, 277, 381, 829, 954, 416, 344, 942, 312, 878, 197, 863, 409, 19, 783, 1026, 177, 544, 1212, 374, 918, 293, 310, 195, 579, 432, 417, 537, 469, 95, 689, 119, 789, 769, 1253, 610, 484, 145, 599, 953, 1333, 982, 592, 231, 202, 532, 1242, 25, 263, 125, 712, 504, 456, 491, 499, 784, 72, 984, 1006, 320, 869, 906, 157, 165, 364, 64, 1265, 113, 631, 536, 1357, 595, 616, 879, 715, 651, 198, 225, 56, 630, 1260, 892, 371, 221, 893, 1232, 264, 357, 219, 791, 600, 346, 44, 940, 1320, 369, 423, 934, 480, 841, 51, 645, 580, 440, 319, 1039, 542, 876, 905, 426, 314, 15, 227, 815, 625, 877, 891, 576, 994, 138, 186, 517, 790, 358, 107, 684, 175, 390, 173, 665, 939, 1077, 124, 840, 185, 806, 485, 413, 619, 662, 108, 300, 935, 696, 964, 971, 912, 69, 232, 786, 821, 1340, 144, 848, 602, 83, }
print len(dti_colivri)



In [40]:

    
saved - dti_colivri









    Out[40]:





{9, 301, 500, 559, 965}



In [41]:

    
tracula_disc = {1005, 1021, 1026, 1039, 1049, 107, 1077, 108, 113, 119, 121, 1212, 1221, 1232, 124, 125, 1253, 1260, 1265, 128, 129, 1320, 1326, 1333, 1338, 1340, 1357, 138, 141, 143, 144, 145, 15, 151, 153, 154, 156, 157, 165, 173, 175, 176, 177, 182, 185, 186, 195, 197, 198, 201, 202, 205, 216, 219, 221, 225, 227, 230, 231, 232, 235, 237, 25, 253, 256, 263, 264, 266, 277, 29, 292, 293, 300, 307, 31, 310, 313, 314, 319, 320, 322, 327, 331, 332, 333, 344, 369, 371, 374, 381, 390, 396, 399, 409, 413, 416, 417, 423, 426, 429, 431, 432, 44, 440, 452, 456, 458, 464, 469, 472, 478, 480, 483, 484, 485, 491, 499, 504, 51, 517, 526, 532, 535, 536, 537, 54, 542, 544, 545, 547, 548, 549, 552, 56, 566, 576, 579, 580, 592, 593, 595, 599, 600, 602, 610, 611, 616, 619, 623, 625, 630, 64, 645, 65, 651, 662, 665, 670, 675, 684, 686, 689, 69, 696, 71, 712, 715, 72, 734, 754, 761, 765, 769, 783, 784, 786, 789, 790, 791, 804, 806, 815, 818, 821, 829, 83, 840, 841, 848, 850, 861, 863, 868, 869, 874, 876, 877, 878, 879, 884, 891, 892, 893, 894, 898, 905, 906, 912, 918, 928, 934, 935, 939, 940, 942, 95, 953, 954, 964, 966, 971, 982, 984, 992, 994, }
print len(tracula_disc)



In [42]:

    
saved.intersection(dti_colivri - tracula_disc)









    Out[42]:





{19, 346, 353, 356, 357, 364, 427, 631, 1006, 1242}

19 done
346 Done
353 done
356 done
357 done
364 done
427 problems, doesnt converge
631 Done
1006 done
1242 No bvec/bval

    15480.pts-7.guitaca     (Detached)



In [43]:

    
camino_tracks={292, 153, 593, 874, 353, 552, 322, 331, 333, 154, 483, 348, 141, 429, 761, 966, 216, 176, 548, 54, 884, 143, 1338, 156, 1021, 431, 611, 850, 458, 327, 29, 623, 128, 547, 861, 313, 332, 399, 1326, 129, 452, 151, 307, 670, 256, 396, 992, 464, 754, 545, 675, 230, 478, 123, 235, 686, 65, 237, 526, 804, 266, 734, 427, 894, 898, 205, 549, 765, 71, 201, 566, 182, 535, 928, 868, 121, 355, 253, 1221, 31, 818, 472, 35, 277, 381, 829, 954, 416, 344, 942, 878, 197, 863, 409, 783, 1026, 177, 544, 1212, 374, 918, 293, 310, 195, 579, 432, 417, 537, 469, 95, 689, 119/130, 119/90, 119, 119/120, 119/110, 119/100, 789, 769, 1253, 610, 484, 145, 599, 953, 1333, 982, 592, 231, 202, 532, 25, 263, 125, 712, 504, 456, 491, 499, 784, 72, 984, 320, 869, 906, 157, 165, 364, 64, 1265, 113, 631, 536, 1357, 595, 616, 879, 715, 651, 198, 225, 56, 630, 1260, 892, 371, 221, 893, 1232, 264, 357, 219, 791, 600, 346, 44, 940, 1320, 369, 423, 934, 480, 841, 51, 645, 580, 440, 319, 1039, 542, 876, 905, 426, 314, 15, 227, 815, 625, 877, 891, 576, 138, 186, 517, 790, 358, 107, 684, 175, 390, 173, 665, 939, 124, 840, 185, 806, 485, 413, 619, 662, 108, 300, 935, 696, 964, 971, 912, 69, 232, 786, 821, 1340, 144, 848, 602, 83, }
print len(camino_tracks)



In [44]:

    
saved.intersection(dti_colivri-camino_tracks)









    Out[44]:





{19, 356, 1005, 1006, 1049, 1077, 1242}

19 : Done, problematic, no tracts
356 : Done
1005: Done
1006: Done
1049: Done
1077: Done
1242: no bval



In [45]:

    
19 in all_aparc









    Out[45]:





True



In [46]:

    
356 in all_aparc









    Out[46]:





False



In [47]:

    
models = { 1357, 198, 29, 364, 535, 675, 874, 1021, 1026, 1039, 1049, 107, 1077, 108, 113, 119, 121, 1212, 1221, 123, 1232, 124, 1242, 125, 1253, 1260, 1265, 128, 129, 1320, 1326, 1333, 1338, 1340, 292, 293, 300, 301, 307, 31, 310, 313, 314, 319, 320, 322, 327, 331, 332, 333, 344, 346, 348, 35, 353, 355, 357, 358, 138, 141, 143, 144, 145, 15, 151, 153, 154, 156, 157, 165, 173, 175, 176, 177, 182, 185, 186, 195, 197, 369, 371, 374, 381, 390, 396, 399, 409, 413, 416, 417, 423, 426, 427, 429, 431, 432, 44, 440, 452, 456, 458, 464, 469, 472, 478, 480, 483, 484, 485, 491, 499, 500, 504, 51, 517, 526, 532, 536, 537, 54, 542, 544, 545, 547, 548, 549, 552, 559, 56, 566, 576, 579, 580, 592, 593, 595, 599, 600, 602, 610, 611, 616, 619, 623, 625, 630, 631, 64, 645, 65, 651, 662, 665, 670, 684, 686, 689, 69, 696, 71, 712, 715, 72, 73, 730, 734, 754, 761, 765, 769, 783, 784, 786, 789, 790, 791, 804, 806, 815, 818, 821, 829, 83, 840, 841, 848, 850, 861, 863, 868, 869, 876, 877, 878, 879, 884, 891, 892, 893, 894, 898, 9, 905, 906, 912, 918, 928, 934, 935, 939, 940, 942, 95, 953, 954, 964, 965, 966, 971, 982, 984, 992, 994, 20, 201, 202, 205, 216, 219, 221, 225, 227, 230, 231, 232, 235, 237, 25, 253, 256, 263, 264, 266, 277, }
print len(models)



In [48]:

    
saved - models









    Out[48]:





{19, 356, 1005, 1006}



In [49]:

    
str(saved.intersection(dti_colivri))









    Out[49]:





'set([1026, 939, 517, 526, 15, 19, 535, 536, 537, 29, 542, 544, 545, 548, 549, 689, 552, 44, 861, 51, 1077, 566, 56, 576, 65, 579, 580, 69, 71, 592, 593, 83, 599, 600, 868, 602, 1039, 357, 610, 616, 619, 108, 623, 113, 630, 631, 121, 124, 125, 789, 128, 107, 876, 138, 651, 141, 143, 144, 145, 25, 665, 154, 156, 157, 670, 197, 675, 165, 625, 684, 173, 686, 175, 176, 177, 182, 185, 1212, 195, 1221, 198, 712, 201, 715, 205, 1232, 216, 804, 1242, 219, 221, 734, 225, 227, 1253, 230, 231, 232, 235, 1260, 237, 1265, 754, 761, 253, 256, 769, 263, 266, 783, 784, 984, 1049, 277, 790, 791, 964, 292, 293, 806, 1320, 1326, 815, 818, 307, 1333, 310, 313, 314, 1340, 829, 320, 480, 54, 327, 840, 841, 331, 332, 1357, 848, 344, 346, 1338, 863, 353, 356, 869, 874, 364, 877, 878, 879, 369, 371, 884, 891, 892, 893, 894, 64, 390, 905, 906, 151, 912, 918, 153, 409, 413, 416, 417, 934, 423, 426, 427, 940, 429, 942, 431, 432, 935, 440, 953, 954, 928, 452, 966, 971, 333, 469, 982, 472, 478, 992, 483, 484, 485, 491, 1005, 1006, 765, 595, 504, 119, 1021])'



In [50]:

    
old_freesurfer=[517,544,545,559,1039,992,891,892,953,964,765]



In [51]:

    
str(sorted((saved-set(old_freesurfer))))









    Out[51]:





'[9, 15, 19, 25, 29, 44, 51, 54, 56, 64, 65, 69, 71, 83, 107, 108, 113, 119, 121, 124, 125, 128, 138, 141, 143, 144, 145, 151, 153, 154, 156, 157, 165, 173, 175, 176, 177, 182, 185, 195, 197, 198, 201, 205, 216, 219, 221, 225, 227, 230, 231, 232, 235, 237, 253, 256, 263, 266, 277, 292, 293, 301, 307, 310, 313, 314, 320, 327, 331, 332, 333, 344, 346, 353, 356, 357, 364, 369, 371, 390, 409, 413, 416, 417, 423, 426, 427, 429, 431, 432, 440, 452, 469, 472, 478, 480, 483, 484, 485, 491, 500, 504, 526, 535, 536, 537, 542, 548, 549, 552, 566, 576, 579, 580, 592, 593, 595, 599, 600, 602, 610, 616, 619, 623, 625, 630, 631, 651, 665, 670, 675, 684, 686, 689, 712, 715, 734, 754, 761, 769, 783, 784, 789, 790, 791, 804, 806, 815, 818, 829, 840, 841, 848, 861, 863, 868, 869, 874, 876, 877, 878, 879, 884, 893, 894, 905, 906, 912, 918, 928, 934, 935, 939, 940, 942, 954, 965, 966, 971, 982, 984, 1005, 1006, 1021, 1026, 1049, 1077, 1212, 1221, 1232, 1242, 1253, 1260, 1265, 1320, 1326, 1333, 1338, 1340, 1357]'



In [52]:

    
str(saved.intersection(dti_colivri)-set(old_freesurfer))









    Out[52]:





'set([1026, 427, 526, 15, 19, 535, 536, 537, 29, 542, 548, 549, 177, 552, 44, 51, 1077, 566, 56, 576, 65, 579, 580, 69, 71, 592, 593, 83, 599, 600, 356, 602, 869, 610, 616, 107, 108, 623, 625, 630, 631, 121, 124, 125, 277, 128, 619, 364, 138, 651, 141, 143, 144, 145, 25, 665, 154, 156, 157, 670, 1221, 675, 165, 113, 684, 173, 686, 175, 176, 689, 182, 185, 1212, 195, 197, 198, 712, 201, 715, 205, 1232, 216, 292, 1242, 219, 221, 734, 225, 227, 1253, 230, 231, 232, 235, 1260, 237, 1265, 754, 761, 1338, 253, 256, 769, 263, 266, 783, 784, 472, 1049, 789, 790, 791, 804, 293, 806, 1320, 1326, 815, 818, 307, 1333, 310, 313, 314, 1340, 829, 320, 54, 327, 840, 841, 331, 332, 1357, 848, 344, 346, 861, 863, 353, 868, 357, 874, 876, 877, 878, 879, 369, 371, 884, 893, 894, 64, 491, 390, 905, 906, 151, 912, 918, 153, 409, 413, 416, 417, 934, 423, 426, 939, 940, 429, 942, 431, 432, 440, 954, 928, 452, 966, 971, 333, 469, 982, 984, 478, 480, 483, 484, 485, 935, 1005, 1006, 595, 504, 119, 1021])'

fMRI



In [53]:

    
with_atencion_raw={ 1265,   176,   327 ,   485,   670,   804,   1021,   1039,   107,   1077,   108,   113,   119,   121,   1212,   1221,   123,   1253,   128,   1333,   1340,   138,   141,   143,   15,   151,   153,   156,   173,   175,   182,   185,   186,   19,   197,   198,   20,   202,   219,   221,   227,   230,   231,   235,   25,   256,   263,   266,   277,   29,   31,   313,   314,   331,   333,   344,   348,   35,   355,   356,   364,   371,   390,   396,   409,   416,   417,   423,   427,   431,   44,   452,   456,   458,   469,   484,   491,   500,   51,   517,   535,   536,   542,   547,   548,   552,   559,   56,   566,   579,   580,   595,   599,   600,   602,   610,   625,   630,   631,   64,   645,   65,   675,   686,   689,   69,   696,   712,   730,   734,   754,   761,   789,   790,   806,   821,   829,   83,   841,   848,   878,   884,   894,   898,   9,   912,   928,   934,   935,   939,   940,   942,   954,   984,   992,  }



In [54]:

    
len(with_atencion_raw)









    Out[54]:





135



In [55]:

    
with_atencion=(with_atencion_raw-pilotos).intersection(saved)



In [56]:

    
len(with_atencion)









    Out[56]:





118



In [57]:

    
for_matlab=','.join(map(lambda x:"'%s'"%x,sorted(with_atencion)))



In [58]:

    
print for_matlab









    



'9','15','19','25','29','44','51','56','64','65','69','83','107','108','113','119','121','128','138','141','143','151','153','156','173','175','176','182','185','197','198','219','221','227','230','231','235','256','263','266','277','313','314','327','331','333','344','356','364','371','390','409','416','417','423','427','431','452','469','484','485','491','500','517','535','536','542','548','552','559','566','579','580','595','599','600','602','610','625','630','631','670','675','686','689','712','734','754','761','789','790','804','806','829','841','848','878','884','894','912','928','934','935','939','940','942','954','984','992','1021','1039','1077','1212','1221','1253','1265','1333','1340'



In [59]:

    
with_coordinacion_raw = {1265 ,   176 ,   237 ,   327 ,    485 ,   576 ,   670 ,   804 ,   893 ,   1005 ,   1006 ,   1021 ,   1026 ,   1039 ,   1049 ,   107 ,   1077 ,   108 ,   113 ,   119 ,   121 ,   1212 ,   1221 ,   123 ,   1232 ,   124 ,   1242 ,   125 ,   1253 ,   1260 ,   128 ,   129 ,   1320 ,   1333 ,   1338 ,   1340 ,   1357 ,   138 ,   141 ,   143 ,   144 ,   145 ,   15 ,   151 ,   153 ,   154 ,   156 ,   157 ,   165 ,   173 ,   175 ,   177 ,   182 ,   185 ,   186 ,   19 ,   195 ,   197 ,   198 ,   20 ,   201 ,   202 ,   205 ,   216 ,   219 ,   221 ,   225 ,   227 ,   230 ,   231 ,   232 ,   235 ,   25 ,   253 ,   256 ,   263 ,   266 ,   277 ,   29 ,   292 ,   293 ,   300 ,   307 ,   31 ,   310 ,   313 ,   314 ,   319 ,   320 ,   322 ,   331 ,   332 ,   333 ,   344 ,   346 ,   348 ,   35 ,   353 ,   355 ,   356 ,   357 ,   358 ,   364 ,   369 ,   371 ,   374 ,   381 ,   390 ,   396 ,   399 ,   409 ,   413 ,   416 ,   417 ,   423 ,   427 ,   429 ,   431 ,   432 ,   44 ,   440 ,   452 ,   456 ,   458 ,   464 ,   469 ,   472 ,   478 ,   480 ,   483 ,   484 ,   491 ,   499 ,   500 ,   504 ,   51 ,   517 ,   535 ,   536 ,   537 ,   54 ,   542 ,   544 ,   545 ,   547 ,   548 ,   549 ,   552 ,   559 ,   56 ,   566 ,   579 ,   580 ,   592 ,   593 ,   595 ,   599 ,   600 ,   602 ,   610 ,   611 ,   616 ,   619 ,   623 ,   625 ,   630 ,   631 ,   64 ,   645 ,   65 ,   662 ,   665 ,   675 ,   684 ,   686 ,   689 ,   69 ,   696 ,   71 ,   712 ,   715 ,   72 ,   730 ,   734 ,   754 ,   761 ,   765 ,   769 ,   783 ,   784 ,   786 ,   789 ,   790 ,   791 ,   806 ,   815 ,   818 ,   821 ,   829 ,   83 ,   840 ,   841 ,   848 ,   850 ,   861 ,   863 ,   868 ,   869 ,   874 ,   876 ,   877 ,   878 ,   879 ,   884 ,   891 ,   892 ,   894 ,   898 ,   9 ,   905 ,   906 ,   912 ,   918 ,   928 ,   934 ,   935 ,   939 ,   940 ,   942 ,   95 ,   953 ,   954 ,   964 ,   965 ,   966 ,   971 ,   982 ,   984 ,   992 ,   994 ,  }
len(with_coordinacion_raw)









    Out[59]:





240



In [60]:

    
with_coordinacion = saved.intersection(with_coordinacion_raw)



In [61]:

    
len(with_coordinacion)









    Out[61]:





206



In [62]:

    
len(saved)









    Out[62]:





211



In [63]:

    
for_matlab=','.join(map(lambda x:"'%s'"%x,sorted(with_coordinacion)))



In [64]:

    
print for_matlab









    



'9','15','19','25','29','44','51','54','56','64','65','69','71','83','107','108','113','119','121','124','125','128','138','141','143','144','145','151','153','154','156','157','165','173','175','176','177','182','185','195','197','198','201','205','216','219','221','225','227','230','231','232','235','237','253','256','263','266','277','292','293','307','310','313','314','320','327','331','332','333','344','346','353','356','357','364','369','371','390','409','413','416','417','423','427','429','431','432','440','452','469','472','478','480','483','484','485','491','500','504','517','535','536','537','542','544','545','548','549','552','559','566','576','579','580','592','593','595','599','600','602','610','616','619','623','625','630','631','665','670','675','684','686','689','712','715','734','754','761','765','769','783','784','789','790','791','804','806','815','818','829','840','841','848','861','863','868','869','874','876','877','878','879','884','891','892','893','894','905','906','912','918','928','934','935','939','940','942','953','954','964','965','966','971','982','984','992','1005','1006','1021','1026','1039','1049','1077','1212','1221','1232','1242','1253','1260','1265','1320','1333','1338','1340','1357'



In [64]: