def test_centered_peak():
img = np.zeros((40, 40))
# Peak centered at 14,24
img[12:17, 22:27] = 1
img[13:16, 23:26] = 2
img[14:15, 24:25] = 3
expected = peak_local_max(img)
assert_equal(expected[0], (14, 24))
#result_gaussian = subpixel_peak(img, method='gaussian', algorithm=algorithm)
#assert_array_almost_equal(result_gaussian, expected)
result_quadrature = subpixel_peak(img, method='quadrature', algorithm=algorithm)
assert_array_almost_equal(result_quadrature, expected)
result_biparabolic = subpixel_peak(img, method='biparabolic', algorithm=algorithm)
assert_array_almost_equal(result_biparabolic, expected)
def test_centered_peak():
img = np.zeros((40, 40))
# Peak centered at 20,24
img[18:23, 22:27] = 1
img[19:22, 23:26] = 2
img[20:21, 24:25] = 3
expected = peak_local_max(img)[0]
assert_equal(expected, (20, 24))
result_gaussian = subpixel_peak(img, method='gaussian', algorithm=algorithm)
assert_array_almost_equal(result_gaussian, expected)
result_quadrature = subpixel_peak(img, method='quadrature', algorithm=algorithm)
assert_array_almost_equal(result_quadrature, expected)
result_biparabolic = subpixel_peak(img, method='biparabolic', algorithm=algorithm)
assert_array_almost_equal(result_biparabolic, expected)
def test_offcentered_peak_biparabolic():
img = np.zeros((40, 40))
# Peak centered at 14,24
img[17:22, 21:26] = 1
img[19:22, 23:26] = 2
img[20:21, 24:25] = 4
expected = peak_local_max(img, min_distance=3)[0]
assert_equal(expected, (20, 24))
result_biparabolic = subpixel_peak(img, method='biparabolic', algorithm=algorithm)
assert_array_less(result_biparabolic, expected)
assert_array_less(expected - 2, result_biparabolic)
def test_offcentered_peak_quadrature():
img = np.zeros((40, 40))
# Peak centered at 20,24
img[17:22, 21:26] = 1
img[19:22, 23:26] = 2
img[20:21, 24:25] = 4
expected = peak_local_max(img, min_distance=3)[0]
assert_equal(expected, (20, 24))
result_quadrature = subpixel_peak(img, method='quadrature', algorithm=algorithm)
assert_array_less(result_quadrature, expected + 1e-3) # + epsilon
assert_array_less(expected - 2, result_quadrature)