R语言提取TIF文件值结果为NA的原因排查

你的提取结果全为NA，核心是代码里的几个关键逻辑错误导致的，逐个拆解：

1. 坐标顺序完全搞反 #

raster::extract()要求传入的坐标是先经度（lon）后纬度（lat），你写的是cbind(lat,lon)，顺序反了——这是返回NA的最主要原因，相当于给了栅格范围外的无效坐标。

2. TIF文件匹配与加载逻辑错误 #

你循环里直接用tif = tif_list[i]按索引取文件，但CSV里的filename和tif_list的顺序未必对应，应该用grep找到的匹配项加载栅格，而非硬套循环索引。另外grep没找到匹配会返回空向量，用same_file > 0判断会报错，正确做法是判断返回结果的长度。

3. 结果拼接逻辑错误 #

每次循环执行ext_data = cbind(data, ext)，会把整个原始数据和当前ext重复绑定，最后只有最后一行结果有效，前面全是NA。正确方式是给data新增一列，然后给当前行单独赋值。

4. 文件路径的潜在问题 #

list.files用了full.names = FALSE，如果工作目录和TIF目录不一致会找不到文件，建议用full.names = TRUE直接获取完整路径。

修正后的代码 #

# 加载依赖包
library(readr)
library(raster)
library(sf)

# 设置工作目录
setwd("D:\\Raju_V\\data\\ocean bottom\\out")

# 读取CSV数据
data <- read_csv("name.csv", show_col_types = FALSE)

# 获取所有TIF文件的完整路径（用\\.tif$精准匹配后缀）
tif_list <- list.files(path = getwd(), pattern = "\\.tif$", full.names = TRUE)

# 新增结果列，初始设为NA
data$ext_value <- NA

# 循环处理每一行
for (i in 1:nrow(data)) {
  filename <- data$filename[i]
  lat <- data$latitude[i]
  lon <- data$longitude[i]
  
  # 精准匹配文件名（fixed=TRUE避免正则匹配干扰）
  same_file <- grep(pattern = filename, x = tif_list, fixed = TRUE)
  
  if (length(same_file) > 0) {
    # 加载匹配到的第一个TIF文件（多匹配场景可自行调整）
    ras <- raster(tif_list[same_file[1]])
    # 修正坐标顺序：先lon后lat
    ext <- extract(ras, cbind(lon, lat), method = 'bilinear')
    # 给当前行赋值结果
    data$ext_value[i] <- ext
    print(ext)
  } else {
    message(paste("未找到匹配文件：", filename))
  }
}

# 保存结果（去掉row.names避免冗余列）
write.csv(data, 'result.csv', row.names = FALSE)
print(data)

额外提示：可以用raster::crs(ras)查看栅格的坐标参考系，确认你的经纬度坐标和栅格CRS是否一致，若不一致需要用sf包的st_transform做坐标转换。

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