```r
pnts<-read.csv(paste(path,"study_points.csv",sep="/"))

#coordinates(pnts) = ~X+Y  # Set coordinates
#pnts<-st_as_sf(pnts)  # Make sf object
#st_crs(pnts) <-27700 ## Set CRS


pnts<-st_as_sf(pnts, coords = c("X", "Y"), crs = 27700) ## Better way
```