读取数据文件\n df <- fread("E:/trainingdata/wheatOUTfile/75csv/Hebei_bazhou.csv")\n\n # 提取起始年份和结束年份\n start_year <- 1981\n end_year <- 2019\n\n # 循环计算1981-2019年的结果\n for (year in start_year:end_year) {\n # 提取播种日期 (day1) 和成熟日期 (day2)\n day1 <- df %%>%\n filter(year == year, StageName == "sowing") %%>%\n select(simulation_days)\n \n day2 <- df %%>%\n filter(year == year+1, StageName == "maturity") %%>%\n select(simulation_days)\n \n avg_temperature <- df %%>%\n filter(simulation_days >= day1[[1]] && simulation_days <= day2[[1]]) %%>%\n summarize(Avg_Temperature = mean(MaxT, na.rm = TRUE))\n \n total_rainfall <- df %%>%\n filter(simulation_days >= day1[[1]] && simulation_days <= day2[[1]]) %%>%\n summarize(Total_Rainfall = sum(Rain, na.rm = TRUE))\n \n gdd <- df %%>%\n filter(simulation_days >= day1[[1]] && simulation_days <= day2[[1]]) %%>%\n summarize(gdd = sum((MaxT+MinT)/2, na.rm = TRUE))\n \n # 提取成熟时的产量\n yield_at_maturity <- df %%>%\n filter(year == year+1, StageName == "maturity") %%>%\n select(yield)\n \n # 输出结果\n print(paste("Year:", year))\n print(paste("Average Temperature:", avg_temperature$Avg_Temperature))\n print(paste("Total Rainfall:", total_rainfall$Total_Rainfall))\n print(paste("GDD:", gdd$gdd))\n print(paste("Yield at Maturity:", yield_at_maturity$yield))\n print("-----------------------------------------------------")\n }