def test_roundness():
Y, X = -24 + np.indices((48, 48)).astype(float)
r = roundness((Y**2. + X**2.) < 4**2.)
assert r > 0
r2 = roundness((Y**2. + 2 * X**2.) < 4**2.)
assert r2 > 0
assert r2 < r
def test_roundness():
Y,X = -24 + np.indices((48,48)).astype(float)
r = roundness( (Y ** 2. + X**2.) < 4**2. )
assert r > 0
r2 = roundness( (Y ** 2. + 2* X**2.) < 4**2. )
assert r2 > 0
assert r2 < r
def test_roundness():
Y,X = -240 + np.indices((480,480)).astype(float)
#X[0][0] = X[0][0] - 1
#print(X[0][0])
#print((Y ** 2. + X**2.) <= 100**2.)
r = roundness( (Y ** 2. + X**2.) <= 100**2.)
print(r)
#print(X)
#print(Y)
#print((Y ** 2. + X**2.) <= 150**2.)
r = roundness( ((Y+1) ** 2. + X**2.) <= 100**2. )
assert r > 0
print(r)
def test_zeros():
assert roundness(np.zeros((10, 10))) == 0
assert eccentricity(np.zeros((10, 10))) == 0
I = np.zeros((16, 16))
I[8:4:12] = 1
assert eccentricity(I) == 0
def test_zeros():
assert roundness(np.zeros((10,10))) == 0
assert eccentricity(np.zeros((10,10))) == 0
I = np.zeros((16,16))
I[8:4:12] = 1
assert eccentricity(I) == 0
