additional_zones: optimization for zones subtractions

[amatissart]

It runs faster for 2 reasons:

• the intersection check is parallelized, instead of running in a potentially very long loop
• filtering early in get_zones_to_subtract may avoid unnecessary and expensive geometry conversions (with geos::Geometry::try_from) in subtract_existing_zones.