Я нашел довольно простой пример того, как это сделать, но не могу заставить его работать на меня. Я довольно новичок в R

library(rgdal) 
xy <- cbind(c(118, 119), c(10, 50)) 
project(xy, "+proj=utm +zone=51 ellps=WGS84") 
          [,1]    [,2] 
[1,] -48636.65 1109577 
[2,] 213372.05 5546301

Но это с примерами чисел. У меня есть тысячи координат, которые я должен преобразовать, и я не могу понять, как получить их из моей таблицы в этот скрипт

В моем наборе данных есть 3 столбца: ID, X и Y. Как я могу преобразовать их, используя это уравнение? Я застрял на этом в течение нескольких недель