- 积分
- 28
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-8-30
- 最后登录
- 1970-1-1
|
发表于 2016-10-6 16:34:17
|
显示全部楼层
本帖最后由 ssy8stephy 于 2016-12-8 15:42 编辑
netCDF4簡明教程
http://blog.csdn.net/ssy8stephy/article/details/52563383
#讀nc文件需要用到netCDF4
import netCDF4 as nc
#打開nc文件
nc_obj=nc.Dataset("C:\\Users\\kong\\PycharmProjects\\rep\\meow")
#看看NC文件中有甚麼
print(nc_obj)
#<span style="background-color: rgb(240, 240, 240);">讀取U的數據值</span>
UU=nc_obj.variables['U'][:]
print(UU)
#關閉文件
nc_obj.close()
#新建文件並寫入數據
da=nc.Dataset("950hPa","w",format="NETCDF4")
da.createDimension("latsize",usize[0])
da.createDimension("lonsize",usize[1])
da.createVariable("u","f8",("latsize","lonsize"))
da.createVariable("v","f8",("latsize","lonsize"))
da.createVariable("z","f8",("latsize","lonsize"))
da.createVariable("lon","f8",("latsize","lonsize"))
da.createVariable("lat","f8",("latsize","lonsize"))
da.createVariable("num","f8",("latsize","lonsize"))
da.variables["u"][:]=output_u
da.variables["v"][:]=output_v
da.variables["z"][:]=output_z
da.variables["lat"][:]=lat
da.variables["lon"][:]=lon
da.variables["num"][:]=output_num
da.description="950"
da.author="kong"
da.createdate="2016-10-05"
da.close() |
|