1 打开和查看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数据绘制的填色等值线图
