Error in the speedyweather.jl test #16
Description
Running the speedyweather.jl test errors with:
julia> A = ConservativeRegridding.intersection_areas(src_polys, dst_polys)
┌ Error: Intersection failed!
│ i1 = 2381
│ i2 = 2244
└ @ ConservativeRegridding ~/development/ConservativeRegridding.jl/src/regridder.jl:64
ERROR: ArgumentError: Length of array must be >= 3 for GeoInterface.Wrappers.LinearRing
Stacktrace:
[1] _length_error(T::Type, f::Function, x::Vector{Tuple{Float64, Float64}})
@ GeoInterface.Wrappers ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:356
[2] #_#32
@ ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:281 [inlined]
[3] LinearRing
@ ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:256 [inlined]
[4] LinearRing
@ ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:172 [inlined]
[5] _broadcast_getindex_evalf
@ ./broadcast.jl:709 [inlined]
[6] _broadcast_getindex
@ ./broadcast.jl:682 [inlined]
[7] getindex
@ ./broadcast.jl:636 [inlined]
[8] macro expansion
@ ./broadcast.jl:1004 [inlined]
[9] macro expansion
@ ./simdloop.jl:77 [inlined]
[10] copyto!
@ ./broadcast.jl:1003 [inlined]
[11] copyto!
@ ./broadcast.jl:956 [inlined]
[12] copy
@ ./broadcast.jl:928 [inlined]
[13] materialize
@ ./broadcast.jl:903 [inlined]
[14] (GeoInterface.Wrappers.Polygon{nothing, nothing})(geom::Vector{Vector{Tuple{Float64, Float64}}}; extent::Nothing, crs::Nothing)
@ GeoInterface.Wrappers ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:295
[15] Polygon
@ ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:256 [inlined]
[16] Polygon
@ ~/.julia/packages/GeoInterface/4tyIo/src/wrappers.jl:172 [inlined]
[17] _trace_polynodes(alg::GeometryOps.FosterHormannClipping{…}, ::Type{…}, a_list::Vector{…}, b_list::Vector{…}, a_idx_list::Vector{…}, f_step::typeof(GeometryOps._inter_step), poly_a::GeoInterface.Wrappers.Polygon{…}, poly_b::GeoInterface.Wrappers.Polygon{…})
@ GeometryOps ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/clipping_processor.jl:963
[18] _intersection(alg::GeometryOps.FosterHormannClipping{…}, ::TraitTarget{…}, ::Type{…}, ::PolygonTrait, poly_a::GeoInterface.Wrappers.Polygon{…}, ::PolygonTrait, poly_b::GeoInterface.Wrappers.Polygon{…}; exact::GeometryOpsCore.True, kwargs::@Kwargs{})
@ GeometryOps ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:87
[19] _intersection
@ ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:76 [inlined]
[20] #intersection#102
@ ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:45 [inlined]
[21] intersection (repeats 2 times)
@ ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:42 [inlined]
[22] #intersection#103
@ ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:56 [inlined]
[23] intersection (repeats 2 times)
@ ~/.julia/packages/GeometryOps/FkfF4/src/methods/clipping/intersection.jl:53 [inlined]
[24] (::ConservativeRegridding.var"#3#4"{SparseArrays.SparseMatrixCSC{…}, Vector{…}, Vector{…}})(i1::Int64, i2::Int64)
@ ConservativeRegridding ~/development/ConservativeRegridding.jl/src/regridder.jl:62
[25] macro expansion
@ ~/.julia/packages/GeometryOps/FkfF4/src/utils/LoopStateMachine/LoopStateMachine.jl:84 [inlined]
[26] do_dual_query(f::ConservativeRegridding.var"#3#4"{…}, predicate::typeof(Extents.intersects), node1::SortTileRecursiveTree.STRLeafNode{…}, node2::SortTileRecursiveTree.STRLeafNode{…})
@ GeometryOps.SpatialTreeInterface ~/.julia/packages/GeometryOps/FkfF4/src/utils/SpatialTreeInterface/SpatialTreeInterface.jl:171
[27] macro expansion
@ ~/.julia/packages/GeometryOps/FkfF4/src/utils/LoopStateMachine/LoopStateMachine.jl:84 [inlined]
[28] do_dual_query(f::ConservativeRegridding.var"#3#4"{…}, predicate::typeof(Extents.intersects), node1::SortTileRecursiveTree.STRNode{…}, node2::SortTileRecursiveTree.STRNode{…})
@ GeometryOps.SpatialTreeInterface ~/.julia/packages/GeometryOps/FkfF4/src/utils/SpatialTreeInterface/SpatialTreeInterface.jl:191
[29] macro expansion
@ ~/.julia/packages/GeometryOps/FkfF4/src/utils/LoopStateMachine/LoopStateMachine.jl:84 [inlined]
[30] do_dual_query(f::ConservativeRegridding.var"#3#4"{…}, predicate::typeof(Extents.intersects), node1::SortTileRecursiveTree.STRNode{…}, node2::SortTileRecursiveTree.STRNode{…})
@ GeometryOps.SpatialTreeInterface ~/.julia/packages/GeometryOps/FkfF4/src/utils/SpatialTreeInterface/SpatialTreeInterface.jl:191
[31] macro expansion
@ ~/.julia/packages/GeometryOps/FkfF4/src/utils/LoopStateMachine/LoopStateMachine.jl:84 [inlined]
[32] do_dual_query(f::ConservativeRegridding.var"#3#4"{…}, predicate::typeof(Extents.intersects), node1::SortTileRecursiveTree.STRNode{…}, node2::SortTileRecursiveTree.STRNode{…})
@ GeometryOps.SpatialTreeInterface ~/.julia/packages/GeometryOps/FkfF4/src/utils/SpatialTreeInterface/SpatialTreeInterface.jl:191
[33] macro expansion
@ ~/.julia/packages/GeometryOps/FkfF4/src/utils/LoopStateMachine/LoopStateMachine.jl:84 [inlined]
[34] do_dual_query(f::ConservativeRegridding.var"#3#4"{…}, predicate::typeof(Extents.intersects), node1::SortTileRecursiveTree.STRtree{…}, node2::SortTileRecursiveTree.STRtree{…})
@ GeometryOps.SpatialTreeInterface ~/.julia/packages/GeometryOps/FkfF4/src/utils/SpatialTreeInterface/SpatialTreeInterface.jl:191
[35] compute_intersection_areas!(areas::SparseArrays.SparseMatrixCSC{…}, grid1::Vector{…}, grid2::Vector{…}; manifold::Planar, nodecapacity1::Int64, nodecapacity2::Int64)
@ ConservativeRegridding ~/development/ConservativeRegridding.jl/src/regridder.jl:55
[36] compute_intersection_areas!
@ ~/development/ConservativeRegridding.jl/src/regridder.jl:38 [inlined]
[37] #intersection_areas#1
@ ~/development/ConservativeRegridding.jl/src/regridder.jl:35 [inlined]
[38] intersection_areas(src_field::Vector{GeoInterface.Wrappers.Polygon{…}}, dst_field::Vector{GeoInterface.Wrappers.Polygon{…}})
@ ConservativeRegridding ~/development/ConservativeRegridding.jl/src/regridder.jl:26
[39] top-level scope
@ REPL[27]:1
Some type information was truncated. Use `show(err)` to see complete types.