Преобразование данных широты / долготы в 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

Большое спасибо! Лучшая Мара

0 ответов

Другие вопросы по тегам