meteoinfo可以批量导出asc吗？

mrfeng1992

想把grads数据导出为asc，但时间点太多，每个时间点导出有点耗时间，有代码实现批量导出吗？

MeteoInfo

可以写脚本批量导出。

mrfeng1992

王老师可以提供一下脚本的代码吗？

mofangbao

mrfeng1992 发表于 2017-3-13 10:51
王老师可以提供一下脚本的代码吗？

那你要做什么

mrfeng1992

mrfeng1992 发表于 2017-3-13 10:51
王老师可以提供一下脚本的代码吗？

导入到arcgis做栅格数据集

MeteoInfo

mrfeng1992 发表于 2017-3-13 13:32
导入到arcgis做栅格数据集

参考下面的代码：

1. f = addfile('D:/Temp/GrADS/model.ctl')
   
2. for i in range(0, f.timenum()):
   
3.     t = f.gettime(i)
   
4.     print t.strftime('%Y-%M-%d %H')
   
5.     ps = f['PS'][i,'10:60','60:140']
   
6.     outfn = 'D:/Temp/GrADS/test_out_' + t.strftime('%Y%M%d%H') + '.dat'
   
7.     ps.savegrid(outfn)

f = addfile('D:/Temp/GrADS/model.ctl')<br /> for i in range(0, f.timenum()):<br />     t = f.gettime(i)<br />     print t.strftime('%Y-%M-%d %H')<br />     ps = f['PS'][i,'10:60','60:140']<br />     outfn = 'D:/Temp/GrADS/test_out_' + t.strftime('%Y%M%d%H') + '.dat'<br />     ps.savegrid(outfn)

mrfeng1992

MeteoInfo 发表于 2017-3-13 14:30
参考下面的代码：

请问ps = f['PS'][i,'10:60','60:140']中数字是什么意思？

MeteoInfo

mrfeng1992 发表于 2017-3-13 15:30
请问ps = f['PS']中数字是什么意思？

数据是全球的，这个是为了取中国区域数据，经纬度范围限定。

mrfeng1992

MeteoInfo 发表于 2017-3-13 15:41
数据是全球的，这个是为了取中国区域数据，经纬度范围限定。

   ps.savegrid(outfn)这一句找不到savegrid的函数

MeteoInfo

mrfeng1992 发表于 2017-3-13 16:07
ps.savegrid(outfn)这一句找不到savegrid的函数

错误提示是什么？