#######################################
#
# Räumliche Statistik SS 2011
# Blatt 1, Aufgabe 2
#
#####################################

# Nutze die bedingte Gleichverteilungseigenschaft 
# des Poisson-Prozesses, um eine Realisierung
# mit Parameter lambda auf einem rechteckigen
# Beobachtungsfenster zu erzeugen.
getPoissonPP<-function(window_x_lim, window_y_lim, lambda)
{
   windowSize<-abs(window_x_lim[2]-window_x_lim[1])*abs(window_y_lim[2]-window_y_lim[1])
   k<-rpois(1,lambda*windowSize)
   points<-c()

   for(i in 1:k)
   {
      x<-runif(1,min=window_x_lim[1], max=window_x_lim[2])
      y<-runif(1,min=window_y_lim[1], max=window_y_lim[2])

      points<-rbind(points, c(x,y))
   }

   return(points)
}


points<-getPoissonPP(c(0,100), c(0,100), 0.01)
#print(points)
plot(points)