def count_elements_in_all_bounding_regions(track_name, genome='testgenome', allow_overlaps=False):
bounding_regions = BoundingRegionShelve(genome, track_name, allow_overlaps).getAllBoundingRegions()
track = PlainTrack(track_name)
num_elements = 0
for tv in [track.getTrackView(region) for region in bounding_regions]:
num_elements += count_elements(tv)
return num_elements
def get_track_view(track_name, genome_region):
track = PlainTrack(track_name)
return track.getTrackView(genome_region)
