Expanding on Grzegorz Sapijaszko answer. As he points out, there are 60 UTM zones. They are 6 degrees wide, and the zone your data are in can be estimated like this for a longitude:
lon <- 5
ceiling((lon + 180) / 6)
#[1] 31
If that is your zone (have a look here, it is a bit more complex than the formula above), your coordinate reference system would be "+proj=utm +zone=31"
, and you can set that to your data (here with terra
)
library(terra)
crs <- "+proj=utm +zone=31"
p1 <- vect(df, geom=c("EASTING", "NORTHING"), crs=crs)
p1
# class : SpatVector
# geometry : points
# dimensions : 6, 0 (geometries, attributes)
# extent : 535234, 553310, 7528392, 7579475 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs
And then project to lon/lat and plot on top of data you have for the area to see if it correct
p2 <- project(p1, "+proj=longlat")
#plot(x)
#points(p2)