This function forces points of an object of class
[sf] located outside
the limits defined by an object of class sp::SpatialPolygons
to NA.
Details
The function try to apply sf::st_filter(). In case it fails,
it will try to rebuild y according to spherical geometry
(see sf::st_as_s2()) before filtering. If it still fail, it will return
x unfiltered.