estimating overdispersed spatial Poisson process

Dear R users,
I have a "spatial" (1D) grid of points; the hypothesis is that each is generating data according to an overdispersed Poisson prcoess with unknown parameters.
My question is twofold:
- Ignoring the spatial aspect, how can I estimate the two parameters independently for each grid point in R (using glm perhaps ? But how exactly)
- How would I include spatial correlations into this model ?

Thanks a lot,

Markus