学习NC读取心得_nc读取总结

学习NC读取心得由刀豆文库小编整理，希望给你工作、学习、生活带来方便，猜你可能喜欢“nc读取总结”。

打开和查看NC数据

ncid = netcdf.open(‘example.nc’,'nc_write’);% 打开文件，放入内存，记录文件的指针ncid；

[ndims,nvars,ngatts,unlimdimid] = netcdf.inq(ncid);% 返回nc和cdf文件的信息，其中ndims为维度信息，nvars为变量个数信息，ngatts和unlimdimid为整体特征；

[dimname, dimlen] = netcdf.inqDim(ncid,0);% 查询ndims中维度的名称和资料长度，包含0~ndims-1的维度特征；

[varname,xtype,dimids,natts] = netcdf.inqVar(ncid,0);%查询nvars中名称和资料维度，包含0~nvars-1个变量特征；

varid = netcdf.inqVarID(ncid,’varname’);% 返回查询变量varname 在nc和cdf文件中的nvars的位置，是变量的代号；

var= netcdf.getVar(ncid,varid);% 获取变量；返回具体变量的具体内容（开始度。。结尾度）

2操作

1.Create Dimensions函数： dimid = netcdf.defDim(ncid,dimname,dimlen)

2.用ID返回Dimensions的名和长度的函数： [dimname, dimlen] = netcdf.inqDim(ncid,dimid)

3.用名字返回Dimensions的ID号： dimid = netcdf.inqDimID(ncid,dimname)

4.重命名Dimensions： netcdf.renameDim(ncid,dimid,newName)；

f=netcdf(‘c: wrfout_d01_2006-05-06_03’,‘now-

rite’);

lon_start=f{‘LON_LL_D’}(:);

lon_end=f{‘LON_LR_D’}(:);

lat_start=f{‘LAT_LL_D’}(:);

lat_end=f{‘LAT_UR_D’}(:);

nc_rainc=f{‘RAINC’}(:);%取RAINC值

[ny nx nt]=size(nc_rainc);%取南北格点数、东西

向格点数、时间个数

dx=(lon_end-lon_start)/nx;%东西向格距

dy=(lat_end-lat_start)/ny;%南北向格距

使用函数m_contourf()绘制填色等值线图,结果见

图4。代码如下:

m_contourf(lon, lat, nc_rainc);

图4 使用NetCDF数据绘制的填色等值线图