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_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
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
