python（Arcpy）批量nc文件转TIF

python（Arcpy）批量nc文件转TIF

# -*- coding: utf-8 -*-
"""
@Time    : 2022/3/19 15:55
@Author  : FJC
@File    : nc转tif.py
@Software: win10  python3.7
"""
# Name: MakeNetCDFRasterLayer_Ex_02.py
# Description: Create a raster layer from a netCDF file.
# Requirements: None

# Import system modules
import arcpy
import os

# Set local variables
arcpy.CheckOutExtension("ImageAnalyst")  # 检查许可
arcpy.CheckOutExtension("spatial")  # 检查许可

path = r"G:\中国气象模拟数据\全国降雨数据\降雨数据"
out_file = r"G:\中国气象模拟数据\全国降雨数据\降雨数据"

variable = "pre"
XDimension = "lon"
YDimension = "lat"
bandDimmension = "time"

# Execute MakeNetCDFRasterLayer
dirs = os.listdir(path)
print(dirs)
dir_path = []
for dir in dirs:
    path_1 = path + "//" + dir
    dir_path.append(path_1)
print(dir_path)
i = 0
for nc in dir_path:
    print(nc)
    nc_file=dirs[i].strip(".nc")
    arcpy.MakeNetCDFRasterLayer_md(nc, variable, XDimension, YDimension, nc_file, bandDimmension)
    outNC = nc.strip(".nc")
    out_path = outNC + ".tif"
    copy_raster = arcpy.CopyRaster_management(nc_file, out_path, background_value="", nodata_value="")
    print(out_path)
    i += 1

print("over")