1. i:=1
2. for each n: n=0 . . ORDMAX do
3. for each m: 0 <= m <= n , n-m even do
4. Anm[i]:=0
5. for each pixel P(x,y) do
6. Anm[i]:=Anm[i]+P(x,y)*Rnm(x,y)*exp(-j*m*f)
7. endfor
8. i:=i+1
9. endfor
10. endfor
Nota: înmultirea fiecarui Anm[i] cu factorul (n+1)/p nu este aratata aici.