Davide Cittaro avatar 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
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.