How to interpolate sf LINESTRING to another pre-defined grid?

I came up with this solution and it has been tested to interpolate to a finer grid. Maybe someone will find it useful :)

#A function to interpolate  lines to another coordination

lines_interpolation<-function(2Darray,raw_line){

 require(sf)


 DIMNAMES<-dimnames(2Darray)

 lon<-as.numeric(DIMNAMES[[1]])
 lat<-as.numeric(DIMNAMES[[2]])


 interpolated_LINE<-raw_line


 raw_line_crs<-st_coordinates(raw_line)

 interpolated_crs<-raw_line_crs

 for (k in 1:nrow(raw_line_crs)) {
    
  lon_div<-abs(lon-raw_line_crs[k,1])
    
  near_lon_index<-which(lon_div ==min(lon_div))
    
  interpolated_crs[k,1]<-lon[near_lon_index]
    
  lat_div<-abs(lat-raw_line_crs[k,2])
    
  near_lat_index<-which(lat_div ==min(lat_div))
    
  interpolated_crs[k,2]<-lat[near_lat_index]
    
 }

 interpolated_line_crs<-unique(interpolated_crs)
  
 interpolated_line<-st_linestring(x = as.matrix(interpolated_line_crs),
  dim= "XYM")
  

 interpolated_line<-st_as_sf(interpolated_line)  
 #Convert data frame to sf object and assign CRS

 st_crs(interpolated_line)<-"+proj=longlat +datum=WGS84" 


 return(interpolated_line)
}