-
大小: 2KB文件類型: .py金幣: 1下載: 0 次發(fā)布日期: 2021-06-15
- 語言: Python
- 標簽: python2.7??arcpy??GRACE_netCDF??
資源簡介
python_批量讀取netCDF4文件并輸出為tif
python_批量讀取 GRACE grid 數(shù)據(jù) netCDF4文件并輸出為tif
代碼片段和文件信息
#encoding=utf8
import?arcpy
from?arcpy.sa?import?*
import?os
def?extractAllNetCDF():
????
????variable?=?“l(fā)we_thickness“
????x_dimension?=?“l(fā)on“
????y_dimension?=?“l(fā)at“
????band_dimension?=?““
????dimension?=?“time“
????valueSelectionMethod?=?“BY_VALUE“
????#先找好或建好輸出路徑,然后在這里更改輸出路徑
????outLoc?=?r“G:\GSM_GRACE\GSR_grid\data“
????inNetCDF?=?r“G:\GSM_GRACE\GSR_grid\CSR_GRACE.nc“
????#Veriable
????variable?=?“l(fā)we_thickness“
????x_dimension?=?“l(fā)on“
????y_dimension?=?“l(fā)at“
????band_dimension?=?““
????dimension?=?“time“
????valueSelectionMethod?=?“BY_VALUE“
????nc_FP?=?arcpy.NetCDFFileProperties(inNetCDF)
????nc_Dim?=?nc_FP.getDimensions()
????for?dimension?in?nc_Dim:
????????top?=?nc_FP.getDimensionSize(dimension)
????????for?i?in?range(0?top):
????????????if?dimension?==?“time“:
????????????????dimension_values?=?nc_FP.getDimensionValue(dimension?i)
????????????????print?dimension_values[0:9].replace(‘?‘‘‘)
????????????????nowFile?=str(dimension_values[0:9].replace(‘?‘‘‘))
????????????????nowFile?=?nowFile.translate(None?‘/‘)
##????????????????print?nowFile
????????????????if?int(nowFile
評論
共有 條評論