def test_pixel_alloverlap_corr(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 15, 20, 20, 0) Aover = _ellipse_overlap(pt1, pt2) * np.pi * 10**2 npt.assert_almost_equal(_ellipse_overlap(pt1, pt2, return_corr=True), Aover / (np.pi * 10 * 20), decimal=3)
def test_pixel_nooverlap_corr(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 31, 10, 10, 0) npt.assert_almost_equal(_ellipse_overlap(pt1, pt2, return_corr=True), 0.0, decimal=3)
def test_pixel_alloverlap(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 15, 20, 20, 0) npt.assert_almost_equal(_ellipse_overlap(pt1, pt2), 1.0, decimal=3)
def test_pixel_nooverlap(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 31, 10, 10, 0) npt.assert_almost_equal(_ellipse_overlap(pt1, pt2), 0.0, decimal=3)
def test_ellipse_overlap(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 15, 10, 10, 0) npt.assert_almost_equal(_circle_overlap(pt1, pt2), _ellipse_overlap(pt1, pt2), decimal=2)
def test_ellipse_corr(): pt1 = (10, 10, 10, 10, 0) pt2 = (10, 15, 10, 10, 0) npt.assert_almost_equal(_circle_overlap(pt1, pt2, return_corr=True), _ellipse_overlap(pt1, pt2, return_corr=True), decimal=2)
def test_pixel_corr_inside(): pt1 = (10, 10, 10, 8, 0) pt2 = (10, 10, 8, 6, 0) # Reduces to ratio of the radii for this case corr = np.sqrt(48 / 80.) npt.assert_almost_equal(_ellipse_overlap(pt1, pt2, return_corr=True), corr)
def test_pixel_corr_inside(): pt1 = (10, 10, 10, 8, 0) pt2 = (10, 10, 8, 6, 0) # Reduces to ratio of the radii for this case corr = np.sqrt(48/80.) npt.assert_almost_equal(_ellipse_overlap(pt1, pt2, return_corr=True), corr)