Os dados do Paraná são do Prof. Paulo Justiniano. O que te relatei fiquei sabendo quando fiz a disciplina de Geoestatística que ele ministrou na Esalq.
É possível fazer a conversão no R, usando {sp} e {rgdal}. Segue um exemplo para duas coordenadas tomadas no fuso 22 e calculadas para o fuso 23...
sapply(c('rgdal','sp'), require, character.only=T)
fuso22 <- data.frame(X=c(765755.252, 806567.125), Y=c(7726003.670, 7657549.738)); fuso22
coordinates(fuso22) <- ~X+Y
proj4string(fuso22) <- CRS("+proj=utm +zone=22 +south +datum=WGS84 +units=m +no_defs")
fuso22
fuso23 <- spTransform(fuso22, CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"))
fuso23