Add plot recipe for PointPattern

Author: juliohm

Project.toml

@@ -6,11 +6,13 @@ version = "0.1.0"
 [deps]
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 GeoStatsBase = "323cb8eb-fbf6-51c0-afd0-f8fba70507b2"
+RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 
 [compat]
 Distributions = "0.21"
 GeoStatsBase = "0.7"
+RecipesBase = "0.7"
 StaticArrays = "0.12"
 julia = "1"
 

src/PointPatterns.jl

@@ -8,12 +8,16 @@ using GeoStatsBase
 
 using Distributions: Uniform, Poisson, product_distribution
 using StaticArrays: SVector, MVector
+using RecipesBase
 
 import GeoStatsBase: npoints, coordinates, coordinates!
 
 include("pattern.jl")
 include("processes.jl")
 
+# plot recipes
+include("plotrecipes/pattern.jl")
+
 export
   # point pattern
   PointPattern,

src/plotrecipes/pattern.jl

@@ -0,0 +1,30 @@
+# ------------------------------------------------------------------
+# Licensed under the ISC License. See LICENSE in the project root.
+# ------------------------------------------------------------------
+
+@recipe function f(pp::PointPattern{T,N}) where {N,T}
+  X = coordinates(pp)
+
+  seriestype --> :scatter
+  markersize --> 2
+  color --> :black
+  legend --> false
+
+  if N == 1
+    @series begin
+      X[1,:], fill(0, npoints(pp))
+    end
+  elseif N == 2
+    aspect_ratio --> :equal
+    @series begin
+      X[1,:], X[2,:]
+    end
+  elseif N == 3
+    aspect_ratio --> :equal
+    @series begin
+      X[1,:], X[2,:], X[3,:]
+    end
+  else
+    @error "cannot plot in more than 3 dimensions"
+  end
+end