Skip to content

Commit

Permalink
fix issue using numpy factorial and remove abandoned nose test.
Browse files Browse the repository at this point in the history
  • Loading branch information
matthewtownson committed Sep 17, 2024
1 parent 0d3f031 commit 0aef577
Show file tree
Hide file tree
Showing 2 changed files with 16 additions and 4 deletions.
9 changes: 5 additions & 4 deletions aotools/functions/zernike.py
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
import numpy
import math
from . import circle

def phaseFromZernikes(zCoeffs, size, norm="noll", rot=0):
Expand Down Expand Up @@ -88,10 +89,10 @@ def zernikeRadialFunc(n, m, r):
for i in range(0, int((n - m) / 2) + 1):

R += numpy.array(r**(n - 2 * i) * (((-1)**(i)) *
numpy.math.factorial(n - i)) /
(numpy.math.factorial(i) *
numpy.math.factorial(int(0.5 * (n + m) - i)) *
numpy.math.factorial(int(0.5 * (n - m) - i))),
math.factorial(n - i)) /
(math.factorial(i) *
math.factorial(int(0.5 * (n + m) - i)) *
math.factorial(int(0.5 * (n - m) - i))),
dtype='float')
return R

Expand Down
11 changes: 11 additions & 0 deletions test/test_centroiders.py
Original file line number Diff line number Diff line change
@@ -1,5 +1,6 @@
from aotools import image_processing
import numpy
# from nose.tools import raises


def test_centre_of_gravity_single():
Expand Down Expand Up @@ -66,3 +67,13 @@ def test_correlation_many():
com = image_processing.correlation_centroid(im, ref, 0.3)
assert (com.shape[0] == 2)
assert(com.shape[1] == 5)


# @raises(ValueError)
# def test_correlation_error():
# im = numpy.random.random((10))
# ref = numpy.random.random((10, 10))
# com = image_processing.correlation_centroid(im, ref, threshold=0.3)



0 comments on commit 0aef577

Please sign in to comment.