nc文件的读取作图

%%读取nc文件并作图
clear all
cd('C:\Users\FEI ZHU\Desktop\MATLAB课件\数据');%更改路径到nc文件的目录
domain=[70 160 10 60];%限定作图区域
ncid=netcdf.open('air.mon.mean.nc','nowrite');%打开数据
ncdisp('air.mon.mean.nc')%显示数据属性
%获取变量
level=netcdf.getVar(ncid,0,'double');
lat=netcdf.getVar(ncid,1,'double');
lon=netcdf.getVar(ncid,2,'double');
time=netcdf.getVar(ncid,3,'double');
air=netcdf.getVar(ncid,4,'double');
%或lat1=ncread('air.mon.mean.nc','lat');%读取纬度等变量这个方法好。
%寻找经度，纬度 时间的位置（find）
lonk1=find(lon==domain(1));
lonk2=find(lon==domain(2));
latk2=find(lat==domain(4));
latk1=find(lat==domain(3));
%日期转化成日期向量表示
%查看日期存储的方式，判断从哪一年开始以何种方式存储的。
%寻找day1（开始）和day2（结束）在time中的位置
day3=find(time==day1);
day4=find(time==day2);
%从air中提取温度数据   按照维数限定条件
air3=air(lonk1:lonk2,latk1:latk2,:,day2:day3);
%利用上述的限定条件提取数据  经度  纬度 层数  时间
%然后利用相关函数求统计量
%%作图  
%t提取所要范围的经纬度
ladata=lat(latk1:latk2);