arcgis中轉換netCDF為柵格數據

最近有個同學詢問我一個問題，使用arcpy把netcdf轉化成柵格文件，忙活了兩個小時才搞定，其實主要代碼非常簡單，只不過要對arcgis 的功能比較熟悉（其實多思考和查考它的幫助文章，無聊）

 

# -*- coding: utf-8 -*-
import arcpy
from arcpy.sa import *


def NetCDF2Tiff(inNetCDFFile,toTIFFFile):

    variable = "Qs_GDS0_SFC"
    XDimension = "g0_lon_1"
    YDimension = "g0_lat_0"
    outRasterLayer = "gldasTest"
    bandDimmension = ""
    dimensionValues = ""
    valueSelectionMethod = ""

    print "start"
    arcpy.MakeNetCDFRasterLayer_md(inNetCDFFile, variable, XDimension, YDimension,
                                   outRasterLayer, bandDimmension, dimensionValues, 
                                   valueSelectionMethod)
    arcpy.CopyRaster_management(outRasterLayer, toTIFFFile,format="TIFF")

    print "finish"


## 

inNetCDFFile = "C:/Users/Administrator/Documents/GLDAS_NOAH025SUBP_3H.A2014058.0600.001.2017040053742.pss.nc"
toTIFFFile =  "C:/Users/Administrator/Documents/test.tif"

NetCDF2Tiff(inNetCDFFile,toTIFFFile)

還有，值得注意的是，在使用arcpy的時候，使用ipython和notebook的時候，總是運行程序的時候出現錯誤