Преобразование данных широты / долготы в UTM с данными из 2 разных зон UTM
Я пытаюсь преобразовать данные широты и долготы из Италии (UTM 32 и 33) в UTM. Я немного борюсь с этим.
Код, который я пытаюсь использовать, тот, который я нашел в предыдущем посте, я пытался изменить его, чтобы решить проблему нескольких UTM, но, похоже, не работает.
## Sample Data frame. x and y columns represent long and lat respectively.
## df can be any data frame that you want to reproject
df <- data.frame(x = Data_group_malesCLEAN$Latitude, y = Data_group_malesCLEAN$Longitude, ID = Data_group_malesCLEAN$Group.1)
long2UTM <- function(long) {
## Function to get the UTM zone for a given longitude
(floor((long + 180)/6) %% 60) + 1
}
LongLatToUTM <- function(df){
coordinates(df) <- ~x+y
proj4string(df) <- CRS("+proj=longlat +datum=WGS84")
zone <- long2UTM(df$x)
## Change CRS of the spatial data frame and convert to data frame
res1 <- spTransform(df, CRS(paste0("+proj=utm +zone=", zone, "+datum=WGS84")))
return(as.data.frame(res1))
}
res1
Большое спасибо! Лучшая Мара