单文件绘图代码正常，循环处理海洋光谱数据时陷入无限执行？

我需要绘制海洋中可见光波长随深度的变化情况，单文件绘图代码可以正常运行，但尝试循环处理文件列表、逐个绘图并保存为独立文件时，代码陷入无限执行状态，输出异常数字而非图像。以下是可运行的单文件代码、尝试的循环代码，以及spec_bats_ng_normalised的数据结构示例，需要修改循环逻辑实现正确批量绘图并单独保存。

可运行的单文件代码 #

# Create an empty list to store the results 
spec_bats_ng_normalised_long <- list()

# Loop through each dataframe in the list 
for (i in seq_along(spec_bats_ng_normalised)) {
   # Reshape the current dataframe into long format   
   long_df <- spec_bats_ng_normalised[[i]] %>%
     pivot_longer(cols = -`depth(m)`, names_to = "Wavelength", values_to = "Intensity")
        
   long_df$Wavelength <- as.numeric(long_df$Wavelength)  
   # Append the long dataframe to the list   
   spec_bats_ng_normalised_long[[i]] <- long_df 
}

depth_intervals <- seq(0, 200, by = 10)  

# Create the plot outside the loop 
spec_bats_ng_plot <- ggplot() +
  scale_y_continuous(limits=c(0,1)) +
  scale_x_continuous(breaks = seq(400,700,by = 50)) +
  theme_minimal()  

# Loop over depth intervals and add layers to the plot 
for (depth in depth_intervals) {
   # Find the index of the closest depth to the one in depth_intervals   
   closest_index <- which.min(abs(spec_bats_ng_normalised_long[[1]]$`depth(m)` - depth))  
   # Extract the closest depth value   
   closest_depth <- spec_bats_ng_normalised_long[[1]]$`depth(m)`[closest_index]   
   # Filter data for the closest depth   
   filtered_data <- filter(spec_bats_ng_normalised_long[[1]], `depth(m)` == closest_depth)   
   # Add layers to the plot   
   spec_bats_ng_plot <- spec_bats_ng_plot +
     geom_line(data = filtered_data, aes(colour = rev(`depth(m)`), x = Wavelength, y = Intensity)) +
     labs(title = "Spectral Intensity for PAR with depth - BATS (No Gaussian)",
          x = "Wavelength (nm)", y = "Intensity (normalised to surface Irradiance)")
}

# Print the final plot 
print(spec_bats_ng_plot + labs(colour = "Depth (m)", breaks = depth_intervals))  
ggsave("SpecE_BATS_No_G.jpg",spec_bats_ng_plot,width = 12, height = 6, dpi = 300)

尝试的循环代码（存在问题） #

# Create an empty list to store the results 
spec_bats_ng_normalised_long <- list()

# Loop through each dataframe in the list 
for (i in seq_along(spec_bats_ng_normalised)) {
   # Reshape the current dataframe into long format   
   long_df <- spec_bats_ng_normalised[[i]] %>%
     pivot_longer(cols = -`depth(m)`, names_to = "Wavelength", values_to = "Intensity")
        
   long_df$Wavelength <- as.numeric(long_df$Wavelength)  
   # Append the long dataframe to the list   
   spec_bats_ng_normalised_long[[i]] <- long_df 
}

depth_intervals <- seq(0, 150, by = 10)  

# Create the plot outside the loop 
spec_bats_ng_plot <- ggplot() +
  scale_y_continuous(limits=c(0,1)) +
  scale_x_continuous(breaks = seq(400,700,by = 50)) +
  theme_minimal()  

# Loop over each spectral element within the list 
for (spec_data in spec_bats_ng_normalised_long){   
  # Loop over depth intervals and add layers to the plot   
  for (depth in depth_intervals) {     
    # Find the index of the closest depth to the one in depth_intervals     
    closest_index <- which.min(abs(spec_data$`depth(m)` - depth))          
    # Extract the closest depth value     
    closest_depth <- spec_data$`depth(m)`[closest_index]          
    # Filter data for the closest depth     
    filtered_data <- filter(spec_data, `depth(m)` == closest_depth)          
    # Add layers to the plot     
    spec_bats_ng_plot <- spec_bats_ng_plot +
      geom_line(data = filtered_data, aes(colour = `depth(m)`, x = Wavelength, y = Intensity)) +
      labs(title = paste("Spectral Intensity for PAR with depth - BATS (No Gaussian)",i),
           x = "Wavelength (nm)", y = "Intensity (normalised to surface Irradiance)")
  }  
}

# Print the final plot 
print(spec_bats_ng_plot + labs(colour = "Depth (m)", breaks = depth_intervals))
ggsave("SpecE_BATS_No_G.jpg",spec_bats_ng_plot,width = 12, height = 6, dpi = 300)

数据结构示例 #

> dput(limited_list[[1]][1:5, 1:5])
structure(list(`depth(m)` = c(0, 0.8, 1.6, 2.4, 3.2), `400` = c(1, 
0.972531036836584, 0.945817414828785, 0.919838341546468, 0.894573597461042
), `405` = c(1, 0.974393076923276, 0.949442566850803, 0.925131625134213, 
0.901443839849326), `410` = c(1, 0.975436719339077, 0.951477439563053, 
0.928107292566113, 0.905311776650204), `415` = c(1, 0.976315609524299, 
0.953192772729279, 0.930618160888521, 0.908578761970095)), row.names = c(NA, 
5L), class = "data.frame")

问题分析 #

• 图层叠加溢出：所有文件的绘图图层都叠加到同一个spec_bats_ng_plot对象上，随着循环执行，图层数量急剧增加，导致程序卡顿甚至假死，最终输出混乱的叠加图像。
• 变量未定义：外层循环使用spec_data遍历时，标题里的i变量未被正确赋值，会引发错误或使用残留值。
• 重复标签覆盖：内层循环反复添加labs，会多次覆盖标题等标签，导致显示异常。
• 未批量保存：没有为每个文件生成独立的图像文件，最后仅保存了一个叠加所有数据的图。

修改后的批量绘图代码 #

# 转换为长格式数据（原逻辑保留）
spec_bats_ng_normalised_long <- list()

for (i in seq_along(spec_bats_ng_normalised)) {
  long_df <- spec_bats_ng_normalised[[i]] %>%
    pivot_longer(cols = -`depth(m)`, names_to = "Wavelength", values_to = "Intensity")
  
  long_df$Wavelength <- as.numeric(long_df$Wavelength)
  spec_bats_ng_normalised_long[[i]] <- long_df
}

depth_intervals <- seq(0, 150, by = 10)  

# 循环处理每个文件，生成独立图像并保存
for (i in seq_along(spec_bats_ng_normalised_long)) {
  # 为当前文件创建独立的ggplot对象，避免图层叠加
  current_plot <- ggplot() +
    scale_y_continuous(limits = c(0,1)) +
    scale_x_continuous(breaks = seq(400,700,by = 50)) +
    theme_minimal() +
    labs(title = paste("Spectral Intensity for PAR with depth - BATS (No Gaussian)", i),
         x = "Wavelength (nm)", 
         y = "Intensity (normalised to surface Irradiance)",
         colour = "Depth (m)")
  
  # 为当前文件添加各深度的曲线
  spec_data <- spec_bats_ng_normalised_long[[i]]
  for (depth in depth_intervals) {
    closest_index <- which.min(abs(spec_data$`depth(m)` - depth))
    closest_depth <- spec_data$`depth(m)`[closest_index]
    filtered_data <- filter(spec_data, `depth(m)` == closest_depth)
    
    current_plot <- current_plot +
      geom_line(data = filtered_data, aes(colour = factor(`depth(m)`), x = Wavelength, y = Intensity))
  }
  
  # 设置颜色刻度与深度区间匹配
  current_plot <- current_plot + scale_colour_discrete(breaks = depth_intervals)
  
  # 打印当前图像（可选，用于预览）
  print(current_plot)
  
  # 保存为唯一文件名，用索引区分不同文件
  ggsave(paste0("SpecE_BATS_No_G_", i, ".jpg"), 
         current_plot, 
         width = 12, height = 6, dpi = 300)
}

代码说明 #

• 为每个文件创建独立的current_plot对象，彻底避免图层叠加导致的性能问题和图像混乱。
• 将labs和主题设置放在图的初始化阶段，避免重复添加覆盖标签。
• 使用索引i作为文件名后缀，确保每个文件生成唯一的图像文件。
• 将depth(m)转换为因子类型映射颜色，确保深度区间的颜色对应准确，也可根据需求改为连续颜色映射。
• 内层循环仅针对当前文件的数据处理，逻辑清晰，避免数据交叉污染。

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