feat(masking vectors of layers with polygons with joint mask across l…

…ayers)
PoisotLab · Dec 16, 2024 · 7ce7bb5 · 7ce7bb5
1 parent 01445f8
commit 7ce7bb5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/polygons/polygons.jl b/src/polygons/polygons.jl
@@ -78,7 +78,7 @@ function _get_inclusion_from_polygon!(inclusion, layer, multipolygon::GeoJSON.Mu
 end 
 
 function SimpleSDMLayers.mask!(layers::Vector{SDMLayer}, multipolygon::GeoJSON.MultiPolygon) 
-    inclusion = zeros(eltype(first(layers).indices), size(first(layers)))
+    inclusion = .!reduce(.|, [l.indices for l in layers])
     _get_inclusion_from_polygon!(inclusion, first(layers), multipolygon)
     for layer in layers
         layer.indices .&= inclusion