tbl = Table[
  f(Sqrt[x^2+y^2],ArcTan[x,y]),
  {x,-1,1,0.1},{y,-1,1,0.1}
]

NMinimize[Total[(Flatten[Fourier[tbl]]-Flatten[dataset])^2],{$Z_1,Z_2,...$}]