读取数据文件

df <- fread("E:/trainingdata/wheatOUTfile/75csv/Hebei_bazhou.csv")

提取起始年份和结束年份

start_year <- 1981 end_year <- 2019

循环计算1981-2019年的结果

for (year in start_year:end_year) {

提取播种日期（day1）和成熟日期（day2）

day1 <- df %>% filter(year == year, StageName == "sowing") %>% select(simulation_days)

day2 <- df %>% filter(year == year+1, StageName == "maturity") %>% select(simulation_days)

avg_temperature <- df %>% filter(simulation_days >= day1[[1]] & simulation_days <= day2[[1]]) %>% summarize(Avg_Temperature = mean(MaxT, na.rm = TRUE))

total_rainfall <- df %>% filter(simulation_days >= day1[[1]] & simulation_days <= day2[[1]]) %>% summarize(Total_Rainfall = sum(Rain, na.rm = TRUE))

gdd <- df %>% filter(simulation_days >= day1[[1]] & simulation_days <= day2[[1]]) %>% summarize(gdd = sum((MaxT+MinT)/2, na.rm = TRUE))

提取成熟时的产量

yield_at_maturity <- df %>% filter(year == year+1, StageName == "maturity") %>% select(yield)

输出结果

print(paste("Year:", year)) print(paste("Average Temperature:", avg_temperature$Avg_Temperature)) print(paste("Total Rainfall:", total_rainfall$Total_Rainfall)) print(paste("GDD:", gdd$gdd)) print(paste("Yield at Maturity:", yield_at_maturity$yield)) print("-----------------------------------------------------")