#改代码是对TAP最新的PM2.5和臭氧的合并
fun <- function(input,output){
  bb <- dir(input)  
  cd <- paste(input,bb,sep="")
  
  m <- NULL
  
  for(i in 1:length(cd)){
    
    #moni_data <- read.csv(cd[i],header = T) %>%select(PAC,pm25_mean)  #pm25
    moni_data <- read.csv(cd[i],header = T) %>%select(PAC,O3_mean)  #o3
    
    
    w_time <- substring(cd[i],nchar(cd[i])-11)
    date1 <- substring(w_time,1,8)
    date <- cbind(moni_data,date1)
    m <- rbind(m,date)
    print(i)
  }
  
  f <- substring(w_time,1,4)
  filename = paste(output,f,".csv",sep = "") 
  write.csv(m,filename,row.names = F)
}

#pm25  PM2.5_mean PAC  pm25_mean
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2013/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2014/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2015/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2016/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2017/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2018/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2019/","C:/Users/zang/Desktop/result/pm25/")
#fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2020/","C:/Users/zang/Desktop/result/pm25/")

#O3  PAC  O3_mean
#fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2013/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2014/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2015/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2016/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2017/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2018/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2019/","C:/Users/zang/Desktop/result/O3/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2020/","C:/Users/zang/Desktop/result/O3/")