Add BinomialProcess

Author: juliohm

Project.toml

@@ -4,10 +4,12 @@ authors = ["Júlio Hoffimann"]
 version = "0.1.0"
 
 [deps]
+Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 GeoStatsBase = "323cb8eb-fbf6-51c0-afd0-f8fba70507b2"
 
 [compat]
-GeoStatsBase = "^0.7"
+Distributions = "0.21"
+GeoStatsBase = "0.7"
 julia = "1"
 
 [extras]

src/PointPatterns.jl

@@ -4,6 +4,15 @@
 
 module PointPatterns
 
-greet() = print("Hello World!")
+using GeoStatsBase
+
+using Distributions
+
+# point processes
+include("processes.jl")
+
+export
+  PointProcess,
+  BinomialProcess
 
 end # module

src/processes.jl

@@ -0,0 +1,24 @@
+# ------------------------------------------------------------------
+# Licensed under the ISC License. See LICENSE in the project root.
+# ------------------------------------------------------------------
+
+"""
+    PointProcess
+
+A spatial point process.
+"""
+abstract type PointProcess end
+
+"""
+    rand(p, r, n=1)
+
+Generate `n` realizations of spatial point process `p`
+inside spatial region `r`.
+"""
+Base.rand(p::PointProcess, r::AbstractRegion, n::Int=1) =
+  [rand_single(p, r) for i in 1:n]
+
+#-----------------
+# IMPLEMENTATIONS
+#-----------------
+include("processes/binomial.jl")

src/processes/binomial.jl

@@ -0,0 +1,22 @@
+# ------------------------------------------------------------------
+# Licensed under the ISC License. See LICENSE in the project root.
+# ------------------------------------------------------------------
+
+struct BinomialProcess <: PointProcess
+  npoints::Int
+end
+
+function rand_single(p::BinomialProcess, r::RectangleRegion{T,N}) where {N,T}
+  # region configuration
+  lo = lowerleft(r)
+  up = upperright(r)
+
+  # number of points
+  n  = p.npoints
+
+  # product of uniform distributions
+  U = product_distribution([Uniform(lo[i], up[i]) for i in 1:N])
+
+  # return point pattern
+  PointSet(rand(U, n))
+end