nc数据转换为二维

clear all;close all;clc
ncdisp('e:\05_0819_0822.nc') % display the information of a netCDF file
slp1=ncread('e:\05_0819_0822.nc','slp');
[dk di dj]=size(slp1); %获取该变量的维数本例为3维
fid=fopen('e:\05_0819_0822.dat','w'); % 打开一个文件,用于存储变量数据
for kk=1:dk %该循环用于将原来3维的数据销至2维,写到05_0819_0822.dat上.
slp2=squeeze(slp1(kk,:,:))+119765;
%-------------------------------------------------------------------------
%注意，有的nc数据为方便存储，故在实际数据的基础上进行了一些换算(不是所有的nc都有这种换算)
%数据还原换算公式是: 实际数据=nc数据*scale_factor+add_offset.
%scale_factor和add_offset的数值可以从ncbrower得到（具体见readme_nc.txt).
%或在命令窗口中输入ncdump，打开nc文件，据matlab窗口显示可知;
%本数据的scale_factor和add_offset分别为1和119765.
%-------------------------------------------------------------------------
for ii=1:di
fprintf(fid,'%6d ',slp2(ii,:));
fprintf(fid,'\n'); %换行
end
fprintf(fid,'\n'); %不同的kk之间用一个空行隔开,如无需空行，可删去此行
end
fclose(fid)
fclose all
以上可以转换为二维数据
的nc文件数据是全球数据，我通过这个程序可以读出这个文件中的数据并保存在文本文档中，但是全球数据太多。

我不想要这么多。

我只要某个经纬度之间的数据就OK，请问如何改动这个程序让读出的数据位于某个经纬度之间呢？
程序如下：
% 程序的主要目的是将一个海平面气压数据05_0819_0822.nc转成ascii数据文件及绘图，% 写成的ascii通常都比较大，如果我们并不需要这个ascii文件，而只想对数据进行处理或绘图，
% 其他nc文件可类似处理
clear all;close all;clc
%--------------------------------------------------------------------------
%--读取数据
f=netcdf('cma_cgcm3_mit.run5.monthly.tas_a1_commit_5_cgcm3.1_t47_20 51_2100.nc');% open the nc file., 注意有的nc文件后缀不是为.nc而是.cdf, 若要读取多个nc 文件,可在filename上设循环.
tas1=f{'tas'}(:); % 获取要从nc文件中提取的变量vas
lon1=f{'lon'}(:);
lat1=f{'lat'}(:);
%scale_factor=f{'vas'}.scale_factor;
%add_offset=f{'vas'}.add_offset;
%vas_real = vas1*scale_factor(1)+add_offset(1);
vas_real = tas1;
close(f)
[dk di dj]=size(tas1); %获取该变量的维数本例为3维
disp(dk);
disp(di);
disp(dj);
%-------------------------------------------------------------------------
%注意，有的nc数据为方便存储，故在实际数据的基础上进行了一些换算(不是所有的nc都有这种换算)
%数据还原换算公式是: 实际数据=nc数据*scale_factor+add_offset.
%scale_factor和add_offset的数值可以从ncbrower得到（具体见readme_nc.txt).
%或在命令窗口中输入ncdump，打开nc文件，据matlab窗口显示可知;
%本数据的scale_factor和add_offset分别为1和119765.
%-------------------------------------------------------------------------
%--------------------------------------------------------------------------
%-绘图
%for i = 1:dk
% figure
%vas = squeeze(vas_real(i,:,:));
%contourf(lon,lat,vas);
%colorbar
%end
%--------------------------------------------------------------------------
%--存储成ascii文件
fid=fopen('tas.txt','w'); % 打开一个文件,用于存储变量数据
for kk=1:dk %该循环用于将原来3维的数据销至2维,写到05_0819_0822.dat上.
tas_write = squeeze(tas_real(kk,:,:));
for ii=1:di
fprintf(fid,'%6d ',tas_write(ii,:));
fprintf(fid,'\n'); %换行
end
fprintf(fid,'\n'); %不同的kk之间用一个空行隔开,如无需空行，可删去此行
end
fclose(fid)
fclose all
ncdisp Display contents of NetCDF file in Command Window
ncinfo Return information about NetCDF file
ncread Read data from variable in NetCDF file
ncreadatt Read attribute value from NetCDF file
ncwrite Write data to NetCDF file
ncwriteatt Write attribute to NetCDF file
ncwriteschema Add NetCDF schema definitions to NetCDF file
创建NetCDF文件变量NCCREATE
在命令窗口中显示的内容ncdisp NetCDF文件
ncinfo NetCDF文件返回的信息
ncread读取数据NetCDF文件变量
ncreadatt读取属性值从NetCDF文件
ncwrite NetCDF文件写数据
ncwriteatt写NetCDF文件属性
ncwriteschema添加NetCDF模式定义为NetCDF文件。