A hybrid method based on the combination of generalized forward backward method (GFBM) and Green's function for the grounded dielectric slab together with the acceleration of the combination via a discrete Fourier transform (DFT) based algorithm is developed for the efficient and accurate analysis of electromagnetic radiation/scattering from electrically large, irregularly contoured two-dimensional arrays consisting of finite number of probe-fed microstrip patches. In this method, unknown current coefficients corresponding to a single patch are first solved by a conventional Galerkin type hybrid method of moments (MoM)/Green's function technique that uses the grounded dielectric slab's Green's function. Because the current distribution on the microstrip patch can be expanded using an arbitrary number of subsectional basis functions, the patch can have any shape. The solution for the array currents is then found through GFBM, where it sweeps the current computation element by element. The computational complexity of this method, which is originally O (N-tot(2)) (N-tot being the total number of unknowns) for each iteration, is reduced to O(N-tot) using a DFT based acceleration algorithm making use of the fact that array elements are identical and the array is periodic. Numerical results in the form of array current distribution are given for various sized arrays of probe-fed microstrip patches with elliptical and/or circular boundaries, and are compared with the conventional MoM results to illustrate the efficiency and accuracy of the method.