Hello,
thank you very much for this useful code!
I test it first using a Geopackage-file, but the path-row combinations in the result are east of Africa instead of Central Europe (same with geojson). Working with a bounding box leads to correct results.
Could it be, that in the function spatial_filter in the line
geo_json["coordinates"] = [map(reverse, geo_json["coordinates"][1])]
the option reverse is the reason for getting mismatching data?
Best regards,
Christian