The problem of modeling and analyzing point referencedbinary spatial data is addressed. We formulate a hierarchical model introducingspatial eects at the second stage. Rather than capturing thesecond stage spatial association using a Markov random eld specication,we employ a second order stationary Gaussian spatial process. In fact,we introduce a convenient latent Gaussian spatial process making our observeddata a realization of an indicator process on this latent process.Convenient...