Masking layer with polygon hotfix (#595)

* fix(sdt): masking layers with polygons bugfix and simplification

* fix(sdt): simplify fix for tests

* fix(sdt): polygon simplification test has fcs

Author: gottacatchenall

CHANGELOG.md

@@ -5,6 +5,10 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## `v1.9.2` 
+
+- **fixed** a bug for masking layers with polygons with holes
+
 ## `v1.9.1` - 2026-04-09
 
 - **added** a `getindex` method for a vector of layers and a vector of occurrences

Project.toml

@@ -1,6 +1,6 @@
 name = "SpeciesDistributionToolkit"
 uuid = "72b53823-5c0b-4575-ad0e-8e97227ad13b"
-version = "1.9.1"
+version = "1.9.2"
 authors = ["Timothée Poisot <timothee.poisot@umontreal.ca>"]
 
 [deps]

src/boundingbox.jl

@@ -110,7 +110,7 @@ const _SDMPOLY_TYPES = Union{
 }
 
 function boundingbox(fc::T; kwargs...) where {T <: _SDMPOLY_TYPES}
-    return SimpleSDMPolygons.boundingbox(fc.geometry; kwargs...)
+    return SimpleSDMPolygons.boundingbox(fc; kwargs...)
 end
 
 function boundingbox(fc::SimpleSDMPolygons.FeatureCollection; kwargs...)

src/polygons/polygons.jl

@@ -85,27 +85,13 @@ function change_inclusion!(inclusion, layer, polygon::P) where {P}
     chunk_size = max(1, length(grid) ÷ (10 * Threads.nthreads()))
     data_chunks = Base.Iterators.partition(grid, chunk_size)
 
-    coords = SimpleSDMPolygons.GI.coordinates(transformed_polygon.geometry)
 
     tasks = map(data_chunks) do chunk
         Threads.@spawn begin
             for position in chunk
                 coord = (E[position[2]], N[position[1]])
-                val = false
-                if polygon isa MultiPolygon
-                    inpoly = [
-                        [PolygonOps.inpolygon(coord, ci) for ci in c] for
-                        c in coords
-                    ]
-                else
-                    inpoly = [PolygonOps.inpolygon(coord, c) for c in coords]
-                end
-                val = if eltype(inpoly) <: Vector
-                    any(isone, vcat(inpoly...))
-                else
-                    any(isone, inpoly)
-                end
-                inclusion[position] = val
+                inpoly = SimpleSDMPolygons.AG.contains(transformed_polygon.geometry, SimpleSDMPolygons.AG.createpoint(coord))
+                inclusion[position] = inpoly
             end
         end
     end