- 积分
- 67
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-11-15
- 最后登录
- 1970-1-1
|
发表于 2014-1-19 10:59:44
|
显示全部楼层
我用IDL写的NC格式,Grads下用SDFopen d:\a1.cn 提示SDF file has no discernable Y coordinate,请问什么原因,现将IDL代码和tif输入文件和nc写入文件附在下面
pro write_nc
varname='Temperature'
data=read_tiff(d:\a2.tif)
ID_write=NCDF_CREATE(d:\a1.nc, /CLOBBER)
size=size(data,/DIMENSIONS)
print,size[0],size[1]
xdimID=NCDF_DIMDEF(ID_write,'lon',size[0])
ydimID=NCDF_DIMDEF(ID_write,'lat',size[1])
tdimID=NCDF_DIMDEF(ID_write,'time',1)
varID=NCDF_VARDEF(ID_write,varname,[xdimID,ydimID,tdimID],/float)
NCDF_ATTPUT,ID_write,varID,'long_name','z-value'
varID_02=NCDF_VARDEF(ID_write,'lon',xdimID,/double)
varID_03=NCDF_VARDEF(ID_write,'lat',ydimID,/double)
varID_04=NCDF_VARDEF(ID_write,'time',tdimID,/double)
NCDF_CONTROL,ID_write, /ENDEF
NCDF_VARPUT,ID_write,varID,data
NCDF_CLOSE,ID_write
end |
|