- 积分
- 109
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-7
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
楼主近期才开始学习arcgis和python,因为有个netCDF文件想要转换成栅格进行处理,就用了arcgis的Make NetCDF Raster Layer工具,但由于在下需要各个time点上的栅格,于是试着用python写代码。代码如下:
inNetCDFfile=ncfile
>>> ncfp=arcpy.NetCDFFileProperties(inNetCDFfile)
>>> Variables=ncfp.getVariables()
>>> timevariable=Variables[2]
>>> lvariable=Variables[3]
>>> print 'varriable to be processded is ' +lvariable
>>> nsteps=ncfp.getDimensionSize(timevariable)
... print "total time steps to be precessed is :" + str(nsteps)
>>> for i in range(0,nsteps):
... nctime=ncfp.getDimensionValue(timevariable,i)
... dimensionvalue="time "+str(nctime)
... print "time is "+nctime
... outRasterLayer=nctime3 #set output raster
...
... arcpy.MakeNetCDFRasterLayer_md(inNetCDFfile, lwevariable, XDimension, YDimension,outRasterLayer,"",dimensionvalue)
然后上面的步骤单独运行都没有问题,就最后利用Make NetCDF Raster Layer工具出现问题。从报错的原因来看应该是dimensionvariable出错了。想请教下,dimensionvalue="time "+str(nctime)是否是有问题的。如有问题正确的应该怎么写?
|
|