Open
Description
I was trying to do a geojoin involving United States census block level data, and ran into an error. I tracked it down to a specific geometry, where calling simplexify would introduce NaN into the vertices. I have created a minimum working example below.
julia> ring1 = Ring(
Point(0., 0.),
Point(0., 3.),
Point(2., 3.),
Point(2., 2.),
Point(3., 2.),
Point(3., 0.)
);
julia> ring2 = Ring(
Point(1., 1.),
Point(1., 2.),
Point(2., 2.),
Point(2., 1.)
);
julia> pa = PolyArea(ring1, ring2);
julia> simplexify(pa)
10 SimpleMesh{2,Float64}
10 vertices
├─ Point(0.0, 0.0)
├─ Point(0.0, 3.0)
├─ Point(2.0, 3.0)
├─ Point(NaN, NaN)
├─ Point(NaN, NaN)
├─ Point(1.0, 2.0)
├─ Point(1.0, 1.0)
├─ Point(2.0, 1.0)
├─ Point(3.0, 2.0)
└─ Point(3.0, 0.0)
10 elements
├─ Triangle(7, 1, 2)
├─ Triangle(6, 2, 3)
├─ Triangle(5, 6, 3)
├─ Triangle(3, 4, 5)
├─ Triangle(6, 7, 2)
├─ Triangle(10, 1, 7)
├─ Triangle(10, 7, 8)
├─ Triangle(10, 8, 5)
├─ Triangle(9, 10, 5)
└─ Triangle(5, 4, 9)Note the NaN's in the list of vertices above.
Here is a visualization of the two rings composing the PolyArea:
