def test_single_gray_hole(): one_gray_hole = add_holes((200, 200), hole_level=100, nholes=1) test_gray_hole = Projection(one_gray_hole, wcs=wcs.WCS()) test_bubble = BubbleFinder2D(test_gray_hole, beam=Beam(10), channel=0, sigma=40) test_bubble.multiscale_bubblefind(edge_find=False)
def test_random_gray_holes(): np.random.seed(375467546) gray_holes = add_holes((500, 500), hole_level=100, nholes=40, max_corr=0.1, rad_max=40) test_gray_holes = Projection(gray_holes, wcs=wcs.WCS()) scales = 3 * np.arange(1, 8, np.sqrt(2)) test_bubble = BubbleFinder2D(test_gray_holes, beam=Beam(10), sigma=10, channel=0, scales=scales) test_bubble.multiscale_bubblefind(edge_find=True, nsig=3, overlap_frac=0.6) test_bubble.visualize_regions(edges=True) # print(test_bubble.region_params) return test_bubble