NC文件批量裁剪导出为CSV

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

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值