library(raster)
library(ggplot2)
library(rgdal)
library(sf)
library(viridis)
library(ggthemes)
library(ggpubr)
setwd('D:/research/5_AOD/results/Analysis/Slope/test/')
list <- list.files()
#list=list[12:23]
j=0
for(i in list){
j=j+1
# print(i)
str_name=i
raster=raster(str_name)
outname= strsplit(str_name,split