Add multipoints

src/GeometryOps.jl

@@ -17,10 +17,6 @@ const GB = GeometryBasics
 
 const TuplePoint = Tuple{Float64,Float64}
 const Edge = Tuple{TuplePoint,TuplePoint}
-const MultiGeomTrait = Union{
-    GI.MultiPointTrait, GI.MultiCurveTrait,
-    GI.MultiPolygonTrait, GI.GeometryCollectionTrait}
-
 
 include("primitives.jl")
 include("utils.jl")

src/methods/angles.jl

@@ -10,10 +10,11 @@ To provide an example, consider this rectangle:
 ```@example angles
 using GeometryOps
 using GeometryOps.GeometryBasics
-using Makie
+using Makie, CairoMakie
 
 rect = Polygon([Point(0,0), Point(0,1), Point(1,1), Point(1,0), Point(0, 0)])
 f, a, p = poly(rect; axis = (; aspect = DataAspect()))
+f
 ```
 This is clearly a rectangle, with angles of 90 degrees.
 ```@example angles
@@ -53,7 +54,7 @@ function angles(geom, ::Type{T} = Float64; threaded =false) where T <: AbstractF
 end
 
 # Points and single line segments have no angles
-_angles(::Type{T}, ::Union{GI.PointTrait, GI.LineTrait}, geom) where T = T[]
+_angles(::Type{T}, ::Union{GI.PointTrait, GI.MultiPointTrait, GI.LineTrait}, geom) where T = T[]
 
 #= The angles of a linestring are the angles formed by the line. If the first and last point
 are not explicitly repeated, the geom is not considered closed. The angles should all be on

test/methods/angles.jl

@@ -1,4 +1,5 @@
 pt1 = GI.Point((0.0, 0.0))
+mpt1 = GI.MultiPoint([pt1, pt1])
 l1 = GI.Line([(0.0, 0.0), (0.0, 1.0)])
 
 concave_coords = [(0.0, 0.0), (0.0, 1.0), (-1.0, 1.0), (-1.0, 2.0), (2.0, 2.0), (2.0, 0.0), (0.0, 0.0)]
@@ -20,6 +21,7 @@ c1 = GI.GeometryCollection([pt1, l2, p2])
 
 # Points and lines
 @test isempty(GO.angles(pt1))
+@test isempty(GO.angles(mpt1))
 @test isempty(GO.angles(l1))
 
 # LineStrings and Linear Rings