python2.7_GLDAS_nc数据转为tif_批量

python2.7_GLDAS_nc数据转为tif_批量,每年的文件的存放12月的nc数据 可以将NC数据使用Python转为tif或者grid数据，并将nc文件中的数据按照时间逐一提取

# -*- coding: utf-8 -*-
import arcpy
from arcpy.sa import *
from arcpy import env
import os
s=[]
   
arcpy.CheckOutExtension（“spatial“） #检查许可
variable = “SoilMoi40_100cm_inst“ #变量名，40-100cm soil moisture
x_dimension = “lon“ #经度
y_dimension = “lat“ #纬度
band_dimension = ““
dimension = ““
dimension_values=““
valueSelectionMethod = “BY_VALUE“ 
yearinfo=range（20002018） #循环2000-2017年gldas文件下 每年12月 nc 格式的gldas数据
for y in yearinfo:
    inputLoc=r‘G:\GLDAS\GLDAS_NOAH025_V.2.1\DATA\ncData‘+‘\\‘+str（y） # 输入路径
    
    outLoc = r“G:\GLDAS\GLDAS_NOAH025_V.2.1\Tif_data\SM40100“+‘\\‘+str（y） #输出路径
    os.mkdir（outLoc） #建立输出文件夹
    env.workspace=inputLoc #设置工作环境
    prin