I’m getting bizarre results when I intersect some point data with polygons using the sf package in R, and while I think I’ve found a solution, I’m so baffled by the issue that I’d really appreciate it if someone could help me understand what’s going on and hence whether I’m missing any potential issues.