Invert polygons for fun and new functionality
Today’s story is inspired by a feature request opened for BigQuery GIS. The request, as I understand it, is to provide aggregate version of
ST_Intersection. BigQuery GIS has
ST_Union_Agg, but no
ST_Intersection_Agg. Frankly, I still wonder about a use case when this would be needed, but let’s create such a function anyway.
The reason this is interesting is we’ll use complementary polygons here. Usually they cause troubles when a polygon with incorrect orientation is loaded to BigQuery, but here we’ll use them to our advantage.
Short refresher: on a sphere, every polygon has a complementary polygon. For example, a polygon that describes the Earth’s continents would have a complementary polygon that describes the Earth’s oceans. They have the same boundary, but all inner points of the first polygon are outer points of the second polygon, and vice versa.
How can we make a complementary polygon for a polygon? A simple solution in BigQuery is to subtract it from the whole sphere. What’s left is complementary polygon. We can define a function for this:
create temp function ST_Invert(g GEOGRAPHY) AS (
Note that for polygons
ST_Invert function follows all the logic rules of set complement, and
ST_Invert(ST_Invert(foo)) equals to
foo. Be careful, that this only holds true for polygon interior, and does not work for points or lines: the
ST_Invert of any 0 or 1 dimensional shape returns full globe.
How can we use that to create common intersection, i.e.
Intersection = A ∩ B ∩ C ∩ …. ?
Let’s invert this equality, and use logic rules to simplify the result
¬Intersection = ¬(A ∩ B ∩ C) = (¬A) ∪ (¬B) ∪ (¬C) ∪ ….
So complement of intersection is union of complements. And we already have a function to compute union,
Another way to arrive to this formula is graphically.
The points that belong to intersection in the middle are points that belong to every polygon in the set. It’s complement, or points outside of the common intersection are points that are outside of at least one of the input polygon. I.e. complement of the common intersection is the union of complements:
Inverting this formula again we get
Intersection = ¬( (¬A) ∪ (¬B) ∪ (¬C) ∪ …. )
We can now define a formula for intersection:
Again, a reminder that this only works correctly for polygons interiors, don’t use it with points or lines input, or when the intersection result is a point or a line. E.g. for polygons that intersect on a single point this function returns empty geometry.