libspatialjoin integretion

[patrickbr]

This PR aims to add our libspatialjoin to the spatialjoin SERVICE. Currently highly experimental and largely untested, but should work.

Supported geometry types are POINT, MULTIPOINT, LINESTRING, MULTILINESTRING, POLYGON, MULTIPOLYGON. GEOMETRYCOLLECTIONS are not supported atm, but this is mainly because I have to first outfactor the parsing code into its own function in libspatialjoin to re-use it here.

For example, the following query dynamically calculates a self-join on <touches> on all geometries contained in Liechtenstein. To do this, set the qlss:algorithm to <libspatialjoin> and set the qlss:joinType to <touches>. Supported join types are <intersects>, <contains>,<covers>, <crosses>, <touches>, <equals> and <overlaps>.

PREFIX osmkey: <https://www.openstreetmap.org/wiki/Key:>
PREFIX osmrel: <https://www.openstreetmap.org/relation/>
PREFIX ogc: <http://www.opengis.net/rdf#>
PREFIX geo: <http://www.opengis.net/ont/geosparql#>
PREFIX qlss: <https://qlever.cs.uni-freiburg.de/spatialSearch/>

SELECT ?a ?b
WHERE {
  osmrel:1155955 ogc:sfContains ?a .
  ?a geo:hasGeometry/geo:asWKT ?a_location .
  SERVICE qlss: {
    _:config qlss:algorithm <libspatialjoin> ;
             qlss:left ?a_location ;
             qlss:right ?b_location ;
             qlss:payload ?b ;
             qlss:joinType <touches> .
    {
      osmrel:1155955 ogc:sfContains ?b .
      ?b geo:hasGeometry/geo:asWKT ?b_location .
    }
  }
}

Running time to generate the 59,982 rows of the query above (full self-join on everything in Liechtenstein):

2025-01-31 12:34:50.778 - INFO: Done processing query and sending result, total time was 2,155 ms

[sparql-conformance]

Conformance check passed ✅

No test result changes.

Details: https://qlever.cs.uni-freiburg.de/sparql-conformance-ui?cur=226740b8c46076d2caf6a41dab8a67ddd5ee1843&prev=a307781592842f82bdfef78fc47fc832fd37368d