Commits

Davide Cittaro committed 653a6ce

bugfix

Comments (0)

Files changed (1)

 n = imsize - 1
 
 for chrom in csize.keys():
+  if '_' in chrom: continue
   sys.stderr.write("Processing chromosome %s\n" % chrom)
   
   chrom_length = csize[chrom]
   hc = np.zeros((imsize, imsize))
   data = Sample1.get_as_array(chrom, 0, chrom_length)
   data[np.isnan(data)] = 0
-  for i in xrange(0, len(chrom_length)):
+  for i in xrange(0, chrom_length):
     c, r = d2xy(n, int(round((n**2 - 1) * i / chrom_length)))
     hc[r][c] += data[i]
   # dump hilbert curve in image and array
   hc = np.zeros((imsize, imsize))
   data = Sample2.get_as_array(chrom, 0, chrom_length)
   data[np.isnan(data)] = 0
-  for i in xrange(0, len(chrom_length)):
+  for i in xrange(0, chrom_length):
     c, r = d2xy(n, int(round((n**2 - 1) * i / chrom_length)))
     hc[r][c] += data[i]
   # dump hilbert curve in image and array