Union of touching geometries is incorrect

[asinghvi17]

MWE:

using GeometryBasics
import GeometryOps as GO, GeoInterface as GI
p1 = GI.Polygon([decompose(Point2f, Rect2f(0,0,1,1))]) |> GO.fix
p2 = GI.Polygon([decompose(Point2f, Rect2f(1,1,1,1))]) |> GO.fix

julia> GO.intersects(p1, p2)
true

julia> GO.touches(p1, p2)
true

julia> GO.union(p1, p2; target = GI.PolygonTrait)
2-element Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}:
 GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}(GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}[GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}([(0.0, 0.0), (1.0, 0.0), (0.0, 1.0), (1.0, 1.0), (0.0, 0.0)], nothing, nothing)], nothing, nothing)
 GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}(GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}[GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}([(1.0, 1.0), (2.0, 1.0), (1.0, 2.0), (2.0, 2.0), (1.0, 1.0)], nothing, nothing)], nothing, nothing)

This should be a single polygon, that represents a rectangle of height 2 and width 1.

[asinghvi17]

This is actually fine - the geometries were self intersecting and so incorrect.