computation of the Zernike Invariants

[Anonymous]

I am studying your source code and I have a question about this method for the computation of the invariants in the Moments class:

    def CalcInvariants(self):
        momrange=self.GetMomentSize()+1
        self.invariants=[]
        for (i,j) in [(a,b) for a in range(momrange) for b in range(a+1)]:
            momtot=0
            for k in range(-b,b+1,1):
                mom=self.GetMoment(i,j,k)
                if mom<>"VOID":
                    momtot=mom.real**2+mom.imag**2
            self.invariants.append([i,j,momtot])

Shouldn't the line momtot=mom.real**2+mom.imag**2 be momtot+=mom.real**2+mom.imag**2 instead?

Thank you in advance,

Sebastian Daberdaku

[Scott Grandison]

Yes - it looks like you're right. It certainly wouldn't make sense as it is. Thanks for pointing it out.

[Anonymous]