def test_join_no_overlap():
bub1 = Bubble3D.load_bubble("data/no_overlap_join_bubble_1.pkl")
bub2 = Bubble3D.load_bubble("data/no_overlap_join_bubble_2.pkl")
join_regions = join_bubbles([[bub1, bub2]])
# There should be one less region in the output
initial_num_regions = len(bub1.twoD_regions) + len(bub2.twoD_regions)
final_num_regions = len(join_regions[0])
assert final_num_regions == initial_num_regions
def test_join_single_overlap():
'''
Tests the joining function when there are multiple regions within a
channel
'''
bub1 = Bubble3D.load_bubble("data/single_overlap_join_bubble_1.pkl")
bub2 = Bubble3D.load_bubble("data/single_overlap_join_bubble_2.pkl")
join_regions = join_bubbles([[bub1, bub2]])
# There should be one less region in the output
initial_num_regions = len(bub1.twoD_regions) + len(bub2.twoD_regions)
final_num_regions = len(join_regions[0])
assert final_num_regions == initial_num_regions - 1
