Hi Louise Thank you for your response. While waiting for the best solution,
I tried the code below. It works still now. I use it to convert from NetCDF
files to rasters.
import arcpy
from arcpy import env
from arcpy.sa import *
Input_Temp_folder = "E:\Mekong\DATA\JMP_Deeplearning_data\Mekong_temp\Temp_2017"
Out_Temp_folder_Max
="E:\Mekong\DATA\JMP_Deeplearning_data\Mekong_temp\Temp_2017\Max"
Out_Temp_folder_Min
="E:\Mekong\DATA\JMP_Deeplearning_data\Mekong_temp\Temp_2017\Min"
arcpy.env.workspace = Input_Temp_folder
arcpy.env.overwriteOutput = True
arcpy.env.scratchWorkspace = Out_Temp_folder_Max # Create scratch workspace
arcpy.env.scratchWorkspace = Out_Temp_folder_Min # create Scratch workspace
variableTemp1 = "T2MMAX"
variableTemp2="T2MMIN"
XDimension = "lon"
YDimension = "lat"
time = "time"
NetCDFfiles = arcpy.ListFiles("*.nc4.nc")
print NetCDFfiles
for filename in NetCDFfiles:
print ("Processing: " + filename)
day = filename[27:35]# extract date in the filenames
outRasterLayer_Temp_Max = variableTemp1 + "_" + day # create a
time series for names
outRasterLayer_Temp_Min = variableTemp2 + "_" + day # create a
time series for names
InNetCDTemp = arcpy.env.workspace + "/" + filename # get a file
arcpy.MakeNetCDFRasterLayer_md(InNetCDTemp, variableTemp1,
XDimension, YDimension, outRasterLayer_Temp_Max,time)
arcpy.MakeNetCDFRasterLayer_md(InNetCDTemp, variableTemp2,
XDimension, YDimension, outRasterLayer_Temp_Min,time)
outRasterTemp_Folder = Out_Temp_folder_Min
arcpy.SplitRaster_management
(outRasterLayer_Temp_Max,Out_Temp_folder_Max,outRasterLayer_Temp_Max,"POLYGON_FEATURES","TIFF","NEAREST","","","","DEGREES","","","Mekong_Basin","","","")
arcpy.SplitRaster_management
(outRasterLayer_Temp_Min,Out_Temp_folder_Min,outRasterLayer_Temp_Min,"POLYGON_FEATURES","TIFF","NEAREST","","","","DEGREES","","","Mekong_Basin","","","")