R语言栅格迭代赋值错误：NA值被替换至错误行

我编写了一段代码，用于在栅格上绘制坐标，当提取的栅格数值大于4时，将对应行从数据框中移除。目前遇到的问题是：首次处理大于4的值时，子集化功能正常，但后续迭代中，sly_df$Eltelt_ido[value] = i * 16这行代码会把值替换到错误的行，导致匹配混乱。

dnbr_reclass_r是存储栅格文件名的数组，代码如下：

for (i in 1:length(dnbr_reclass_r)){
  DNBR = raster(paste0('E:/Teruletek/Solymos/dNBR_reclass/', dnbr_reclass_r[i]))
  plot(DNBR)
  points(sly_coor)
  
  extracted_values = extract(DNBR, sly_coor)
  
  for(value in 1:length(extracted_values)){
    if(extracted_values[value] > 4){
      sly_df$Eltelt_ido[value] = i * 16
    }
  }
  sly_subset = subset(sly_df, is.na(sly_df$Eltelt_ido))
  sly_coor = data.frame(x = sly_subset$x, y = sly_subset$y)
}

附sly_df数据框结构：

sly_df
   ID        x       y Eltelt_ido
1   1 250947.8 5126356         NA
2   2 250961.2 5126359         NA
3   3 251202.3 5125685         NA
4   4 251127.2 5125503         NA
5   5 251189.7 5125774         NA
6   6 251305.8 5126352         NA
7   7 250935.7 5126549         NA
8   8 250856.0 5126148         NA
9   9 251205.2 5126735         NA
10 10 250955.7 5126356         NA

我曾尝试用ID关联的方案但未成功，以下是可行的解决思路：

问题根源 #

迭代过程中，sly_coor是原数据框的子集，它的行索引和原sly_df的行索引不再一一对应。内层循环用value（sly_coor的行号）直接索引原sly_df的行，必然导致错位。

解决方法：用唯一ID关联行 #

利用sly_df中的ID列（唯一标识每行），替代行号进行匹配，彻底避免索引错位问题。

修改后代码（高效批量更新版） #

for (i in 1:length(dnbr_reclass_r)){
  DNBR = raster(paste0('E:/Teruletek/Solymos/dNBR_reclass/', dnbr_reclass_r[i]))
  plot(DNBR)
  points(sly_coor)
  
  extracted_values = extract(DNBR, sly_coor)
  # 筛选当前子集里提取值>4的ID
  update_ids = sly_subset$ID[extracted_values > 4]
  # 通过ID批量更新原数据框的对应行
  sly_df$Eltelt_ido[sly_df$ID %in% update_ids] = i * 16
  
  # 更新子集和坐标（保留ID用于下一次匹配）
  sly_subset = subset(sly_df, is.na(sly_df$Eltelt_ido))
  sly_coor = data.frame(x = sly_subset$x, y = sly_subset$y)
}

代码说明 #

1. 每次迭代时，sly_subset保留了原数据框的ID，通过update_ids定位需要更新的行，完全依赖唯一ID而非行号，避免错位。
2. 批量更新的写法替代内层循环，既提升效率，又减少出错概率。

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