Skip to content

Instantly share code, notes, and snippets.

@sergpolly
Created February 7, 2023 20:06
Show Gist options
  • Select an option

  • Save sergpolly/c0243679c199f81c49ebfa84661b71d0 to your computer and use it in GitHub Desktop.

Select an option

Save sergpolly/c0243679c199f81c49ebfa84661b71d0 to your computer and use it in GitHub Desktop.
new cooler query engine
f = clr.open('r')
reader = cooler.core._rangequery.CSRReader(f['pixels'], f['indexes/bin1_offset'][:])
sel = cooler.core._rangequery.DirectRangeQuery2D(reader, "count", bbox=[0,1000,0,1000], chunksize=1_000_000)
sel_low = cooler.core._rangequery.FillLowerRangeQuery2D(reader, "count", bbox=[0,1000,0,1000], chunksize=1_000_000)
cooler.core._rangequery.FillLowerRangeQuery2D
plt.imshow(np.log(sel_low.to_array()))
sel.tasks
for g in sel:
print(g)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment