

In [39]:

    
matD = np.load('reference/reads/matrices/ref.depth20.seed0_0001.reads.npz')
plt.figure(figsize=(20,20))
for i,j in matD.iteritems():
    print i
    imshow(j,interpolation='none',cmap=cm.gray)
    savefig(i+'.png')









    



S1_775
S1_774
S1_478
S1_479
S1_771
S1_770
S1_773
S1_772
S1_472
S1_473
S1_470
S1_471
S1_476
S1_477
S1_474
S1_475
S1_924
S1_925
S1_926
S1_776
S1_920
S1_921
S1_922
S1_923
S1_928
S1_929
S1_487
S1_486
S1_485
S1_484
S1_483
S1_482
S1_481
S1_480
S1_847
S1_846
S1_845
S1_844
S1_843
S1_842
S1_489
S1_488
S1_270
S1_271
S1_272
S1_273
S1_274
S1_275
S1_276
S1_277
S1_278
S1_279
S1_373
S1_372
S1_375
S1_374
S1_377
S1_376
S1_506
S1_507
S1_504
S1_505
S1_502
S1_503
S1_500
S1_501
S1_779
S1_508
S1_509
S1_248
S1_979
S1_978
S1_573
S1_572
S1_571
S1_570
S1_577
S1_576
S1_575
S1_574
S1_579
S1_578
S1_322
S1_687
S1_686
S1_971
S1_970
S1_690
S1_323
S1_498
S1_691
S1_898
S1_899
S1_894
S1_895
S1_896
S1_897
S1_890
S1_891
S1_892
S1_893
S1_281
S1_280
S1_283
S1_282
S1_285
S1_284
S1_287
S1_286
S1_289
S1_288
S1_111
S1_110
S1_113
S1_112
S1_115
S1_114
S1_117
S1_116
S1_119
S1_118
S1_848
S1_57
S1_199
S1_56
S1_51
S1_50
S1_739
S1_738
S1_53
S1_731
S1_730
S1_733
S1_732
S1_735
S1_734
S1_737
S1_736
S1_46
S1_47
S1_44
S1_45
S1_42
S1_43
S1_805
S1_804
S1_809
S1_840
S1_48
S1_49
S1_234
S1_235
S1_236
S1_237
S1_230
S1_231
S1_232
S1_233
S1_59
S1_238
S1_239
S1_58
S1_389
S1_609
S1_608
S1_605
S1_604
S1_607
S1_606
S1_601
S1_600
S1_603
S1_602
S1_425
S1_424
S1_427
S1_426
S1_421
S1_420
S1_423
S1_422
S1_542
S1_543
S1_540
S1_541
S1_429
S1_428
S1_544
S1_545
S1_371
S1_370
S1_876
S1_877
S1_874
S1_875
S1_872
S1_873
S1_870
S1_871
S1_977
S1_976
S1_975
S1_974
S1_973
S1_972
S1_878
S1_879
S1_300
S1_301
S1_302
S1_303
S1_304
S1_305
S1_306
S1_307
S1_308
S1_309
S1_69
S1_35
S1_698
S1_34
S1_318
S1_62
S1_63
S1_201
S1_39
S1_200
S1_740
S1_623
S1_185
S1_203
S1_622
S1_202
S1_777
S1_205
S1_150
S1_204
S1_207
S1_741
S1_250
S1_206
S1_369
S1_629
S1_628
S1_641
S1_640
S1_643
S1_642
S1_645
S1_644
S1_647
S1_646
S1_649
S1_648
S1_538
S1_434
S1_91
S1_90
S1_93
S1_92
S1_95
S1_94
S1_97
S1_96
S1_99
S1_98
S1_162
S1_163
S1_164
S1_165
S1_166
S1_167
S1_469
S1_468
S1_742
S1_743
S1_744
S1_745
S1_746
S1_747
S1_461
S1_460
S1_463
S1_462
S1_465






    



---------------------------------------------------------------------------
KeyboardInterrupt                         Traceback (most recent call last)
<ipython-input-39-4deeb00b24e5> in <module>()
      4     print i
      5     imshow(j,interpolation='none',cmap=cm.gray)
----> 6     savefig(i+'.png')

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/pyplot.pyc in savefig(*args, **kwargs)
    686 def savefig(*args, **kwargs):
    687     fig = gcf()
--> 688     res = fig.savefig(*args, **kwargs)
    689     fig.canvas.draw_idle()   # need this if 'transparent=True' to reset colors
    690     return res

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/figure.pyc in savefig(self, *args, **kwargs)
   1563             self.set_frameon(frameon)
   1564 
-> 1565         self.canvas.print_figure(*args, **kwargs)
   1566 
   1567         if frameon:

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/backend_bases.pyc in print_figure(self, filename, dpi, facecolor, edgecolor, orientation, format, **kwargs)
   2230                 orientation=orientation,
   2231                 bbox_inches_restore=_bbox_inches_restore,
-> 2232                 **kwargs)
   2233         finally:
   2234             if bbox_inches and restore_bbox:

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc in print_png(self, filename_or_obj, *args, **kwargs)
    525 
    526     def print_png(self, filename_or_obj, *args, **kwargs):
--> 527         FigureCanvasAgg.draw(self)
    528         renderer = self.get_renderer()
    529         original_dpi = renderer.dpi

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc in draw(self)
    472 
    473         try:
--> 474             self.figure.draw(self.renderer)
    475         finally:
    476             RendererAgg.lock.release()

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/artist.pyc in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/figure.pyc in draw(self, renderer)
   1157         dsu.sort(key=itemgetter(0))
   1158         for zorder, a, func, args in dsu:
-> 1159             func(*args)
   1160 
   1161         renderer.close_group('figure')

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/artist.pyc in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/axes/_base.pyc in draw(self, renderer, inframe)
   2322 
   2323         for zorder, a in dsu:
-> 2324             a.draw(renderer)
   2325 
   2326         renderer.close_group('axes')

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/artist.pyc in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/image.pyc in draw(self, renderer, *args, **kwargs)
    381 
    382         if self._check_unsampled_image(renderer):
--> 383             self._draw_unsampled_image(renderer, gc)
    384         else:
    385             if self._image_skew_coordinate is not None:

/Users/skasinat/anaconda/lib/python2.7/site-packages/matplotlib/image.pyc in _draw_unsampled_image(self, renderer, gc)
    358 
    359         renderer.draw_image(gc, xmin, ymin, im, dxintv, dyintv,
--> 360                             trans_ic_to_canvas)
    361 
    362     def _check_unsampled_image(self, renderer):

KeyboardInterrupt: