why sf :: st_intersection returns duplicate features. I can't understand why st_intersection adds duplicate lines to my sf data frame. sf_points Simple feature collection with 156421 features and 9 fields geometry type: POINT dimension: XY bbox: xmin: 596640.7 ymin: 7317066 xmax: 607113.9 ymax: 7326187 epsg (SRID): 32722 proj4string: +proj=utm +zone=22 +south +datum=WGS84 +units=m +no_defs sf_polygons Simple feature collection with 225 features and 4 fields geometry type: POLYGON dimension: XY bbox: xmin: 543943.6 ymin: 7317069 xmax: 607110.9 ymax: 7342239 epsg (SRID): 32722 proj4string: +proj=utm +zone=22 +south +datum=WGS84 +units=m +no_defs nrow(sf_points[!duplicated(sf_points), ]) == nrow(sf_points) > TRUE st_crs(sf_points) == st_crs(parcelas) # test crs > TRUE sf_points_intersection <- st_intersection(sf_points, sf_polygon) sf_points_intersection Simple feature collection with 121716 features and 9 fields geometry type: POINT dimension: XY bbox: xmin: 596646.7 ymin: 7317070 xmax: 607109.9 ymax: 7326184 epsg (SRID): 32722 proj4string: +proj=utm +zone=22 +south +datum=WGS84 +units=m +no_defs nrow(sf_points_intersection[!duplicated(sf_points_intersection), ]) == nrow(sf_points_intersection) > FALSE nrow(sf_points_intersection[duplicated(as.data.frame(sf_points_intersection)), ]) > 37056 I appreciate some help in understanding this behavior.