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的时候,总是运行程序的时候出现错误