Значения XY исчезают в таблице атрибутов при использовании функций «координаты», «proj4string» и «writeOGR» в R для создания ESRI .shp

Question

У меня есть фрейм данных в R, который я хотел бы сохранить в качестве шейп-файла ESRI. Фрейм данных имеет X, Y и Z в качестве столбцов. Когда я использую следующий код

coordinates(df)=~X+Y
proj4string(df) <- CRS("+init=epsg:25832")
writeOGR(df, dsn=getwd(), "name", driver="ESRI Shapefile", overwrite_layer = TRUE)

, шейп-файл создается, но когда я открываю его, например, в ArcGIS, таблица атрибутов содержит только столбец Z, X и Y исчезли. Может кто-нибудь объяснить мне, почему это происходит, и еще лучше, как этого избежать? Это связано с выбранным CRS?

Заранее спасибо!