Performance-improvement: Speed up de9im pattern matching

[Croydon-Brixton]

I feel like there has to be a way to make this much more efficient, maybe through stuff like numba or Cython, but that's something to leave for later.

Originally posted by @herbiebradley in #28 (comment)

Further info:
Side note: when searching for what the relate function does we came to this: https://github.com/libgeos/geos/blob/186bbd32fbf07d8b5d419cdfd64c14e2270a418a/src/geom/IntersectionMatrix.cpp#L89 which are the functions which a lot of Shapely stuff calls.

Overall, we haven't concentrated on this part much bc the string comparison is much faster than the overlap computations and rtree query so the overall speedup in the function would be only around 5%.