renamed intersect to intersect points as not to conflict with R base intersect()