1. initialise vector Anm[1 . . MOMMAX] to 0
2. for each pixel P(x,y) do
3. i:=1
4. for each n: n=0..ORDMAX do
5. for each m: 0 <= m <= n, n-m even do
6. Anm[i]:=Anm[i]+P(x,y)*Rnm(x,y)*exp(-j*m*f)
7. i:=i+1
8. endfor
9. endfor
10. endfor