经纬度转UTM（27区）失败求助：含缺失空间数据的R代码问题

你的代码存在几个关键错误，导致转换失败或出错：

1. 函数名大小写错误：定义的带号计算函数是long2utm，但调用时写成了long2UTM，R语言大小写敏感，会直接报错找不到函数
2. UTM带号逻辑错误：原函数用纬度计算带号，实际上UTM带号由经度决定；且你明确要转换到27区，无需自动计算带号
3. 坐标顺序颠倒：sp包中coordinates()要求先经度（x轴）后纬度（y轴），原代码写反成~lat+long，会导致坐标转换完全错误
4. 未处理NA值：包含NA的行无法被转换为空间对象，直接调用coordinates()会触发报错

先确保加载sp包：

library(sp)

使用以下修正后的函数，它会处理NA值、修正坐标顺序，并直接指定UTM27区：

LongLatToUTM <- function(df, utm_zone = 27) {
  # 拆分有完整坐标和缺失坐标的行
  has_coords <- complete.cases(df$lat, df$long)
  df_valid <- df[has_coords, ]
  df_missing <- df[!has_coords, ]
  
  if (nrow(df_valid) > 0) {
    # 创建空间数据框：注意经度(long)在前，纬度(lat)在后
    coordinates(df_valid) <- ~long + lat
    # 设置原始投影为WGS84经纬度
    proj4string(df_valid) <- CRS("+proj=longlat +datum=WGS84")
    
    # 转换到UTM 27区（南半球需加+south参数）
    utm_crs <- CRS(paste0("+proj=utm +zone=", utm_zone, " +south +datum=WGS84"))
    df_utm <- spTransform(df_valid, utm_crs)
    
    # 转回普通数据框，并重命名UTM坐标列
    df_utm <- as.data.frame(df_utm)
    names(df_utm)[names(df_utm) %in% c("long", "lat")] <- c("utm_easting", "utm_northing")
  } else {
    df_utm <- df_valid
  }
  
  # 合并数据并恢复原始行顺序
  result <- rbind(df_utm, df_missing)
  result <- result[match(df$site, result$site), ]
  
  return(result)
}

调用函数获取转换结果：

brazil_utm_data <- LongLatToUTM(brazil_legacy_data, utm_zone = 27)

• NA值处理：先拆分数据，只转换有完整坐标的行，避免报错
• 坐标顺序修正：~long + lat符合空间数据的x/y轴定义
• 指定UTM带号：直接使用你需要的27区，避免自动计算的逻辑错误
• 南半球参数：南纬区域必须加+south，否则UTM坐标会出现负值错误
• 行顺序保留：合并后通过match恢复原始数据的行顺序，保证数据对应关系正确

内容的提问来源于stack exchange，提问作者Eizy