批量TIF转NC并进行MK检验
这里主要记录一个批量进行tif文件转nc，并且将长序列数据进行mk检验的python代码。有问题随时联系：jia5678912。
import os
import numpy as np
import xarray as xr
from osgeo import gdal, osrdef Search_File(dirname,suffix):'''This function can search all files with the specified suffix in this dir.:param dirname: string, the path need to be searched:param suffix: string, the specified suffix need to be seached:return:filter_list: list, the path list need to be searched.'''filter = [suffix] # 设置过滤后的文件类型 当然可以设置多个类型filter_list = []for maindir, subdir, file_name_list in os.walk(dirname):#print(maindir) #当前主目录for filename in file_name_list:apath = os.path.join(maindir, filename)#合并成一个完整路径portion = os.path.splitext(apath)ext = portion[1] # 获取文件后缀 [0]获取的是除了文件名以外的内容if ext in filter:newname = portion[0] + suffixfilter_list.append((newname,portion[0].split("\\")[-1]))# print(filter_list)return filter_listdef combin():'''数据合并，沿时间轴:return:'''import xarray as xrimport pandas as pdimport os# 假设你的NC文件存放在‘data_folder’目录下data_folder = r"E:\530\VQI\vqi" # 请替换为你的NC文件所在的文件夹路径file_n = Search_File(data_folder,".nc"
)
