#install.packages("dplyr")
library(dplyr)

#该套代码用来将各个年份的每天数据合并到各个年份上

fun <- function(input,output){
  bb <- dir(input)
  cd <- paste(input,bb,sep="")
  
  m <- NULL
  
  for(i in 1:length(cd)){
    #读取各套数据模拟的PM2.5
    moni_data <- read.csv(cd[i],header = T) %>%dplyr::select(PAC,NAME,pm25_mean)
    
    #用来修改区县代码
    moni_data$PAC[moni_data$NAME == "苏州工业园区"] <- "320571"   
    
    
    #写入时间和合并模拟数据
    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)
  }
  
  #删除NAME列
  m <- m[,-2]
  ##写入文件
  f <- substring(w_time,1,4)
  filename = paste(output,f,".csv",sep = "") 
  write.csv(m,filename,row.names = F)
}

#TAP
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2013/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2014/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2015/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2016/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2017/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2018/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2019/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_PM2.5_Result/2901个_county_Result/2020/","G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/")

########################################################
#o3 

fun <- function(input,output){
  bb <- dir(input)
  cd <- paste(input,bb,sep="")
  
  m <- NULL
  
  for(i in 1:length(cd)){
    #读取各套数据模拟的PM2.5
    moni_data <- read.csv(cd[i],header = T) %>%dplyr::select(PAC,NAME,O3_mean)
    
    #用来修改区县代码
    moni_data$PAC[moni_data$NAME == "苏州工业园区"] <- "320571"   
    
    
    #写入时间和合并模拟数据
    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)
  }
  
  #删除NAME列
  m <- m[,-2]
  ##写入文件
  f <- substring(w_time,1,4)
  filename = paste(output,f,".csv",sep = "") 
  write.csv(m,filename,row.names = F)
}
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2013/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2014/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2015/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2016/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2017/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2018/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2019/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")
fun("G:/空气污染清理/TAP_O3_Result_xin/2901个county_Result/2020/","G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/")


###########################################################

#将2013-2020各个单独年份合并成一个
fun <- function(input,output){
  bb <- dir(input)
  cd <- paste(input,bb,sep="")
  
  m <- NULL
  
  for (i in 1:length(cd)){
    r_data <- read.csv(cd[i],header = T)  
    
    m <- rbind(m,r_data)
    print(i)  
  }
  
  ##写入文件
  filename = paste(output,sep = "")     
  write.csv(m,filename,row.names = F)
}

fun('G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/各个年份/',"G:/空气污染清理/TAP_pm25_2013-2020_苏州工业园区修改/TAP_pm25_2013-2020.csv")
fun('G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/各个年份/',"G:/空气污染清理/TAP_o3_2013-2020_苏州工业园区修改/TAP_o3_2013-2020.csv")