import arcpy
from arcpy import env
from arcpy.sa import *
import os
from netCDF4 import Dataset
import csv
import math
import re
tmp_root = r"D:\dataset\tmp"
pre_root = r"D:\dataset\pre"
# 获得文件路径
def get_path(filenames_in):
path_list = []
for root, dirs, files in os.walk(filenames_in):
for file in files:
# print(dir)
filePath = os.path.join(filenames_in, file)
# print(pathDir)
path_list.append(filePath)
return path_list
#定义各项参数
variable = "tmp"
XDimension = "lon"
YDimension = "lat"
outRasterLayer = "pre_Layer"
bandDimmension = "time"
dimensionValues = ""
valueSelectionMethod = ""
cellRegistration = ""
tmp_pathlist = get_path(tmp_root)
pre_pathlist = get_path(pre_root)
for path in tmp_pathlist:
out_path = path.split('\\')[-1]
print(out_path)
print(r'D:\\dataset\\caijiannc\\tmp\\' + out_path[:-3])
print(1)
nc_Ra = arcpy.MakeNetCDFRasterLayer_md(path, variable, XDimension, YDimension,
outRasterLayer, bandDimmension, dimensionValues,
valueSelectionMethod, cellRegistration)
outExtractByMask = ExtractByMask('pre_Layer', r"D:\高程\高程图裁剪.tif")
print("D:\dataset\caijiannc\tmp" + out_path)
outExtractByMask.save(f'D:\\dataset\\{out_path[:-3]}' )
print(out_path+'已裁剪成功')
break
NC文件批量裁剪导出为CSV
最新推荐文章于 2024-08-26 15:56:48 发布