该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
我的nc文件数据是全球数据,我通过这个程序可以读出这个文件中的数据并保存在文本文档中,但是全球数据太多。我不想要这么多。我只要某个经纬度之间的数据就OK,请问如何改动这个程序让读出的数据位于某个经纬度之间呢?
程序如下:
% 程序的主要目的是将一个海平面气压数据05_0819_0822.nc转成ascii数据文件及绘图,
% 写成的ascii通常都比较大,如果我们并不需要这个ascii文件,而只想对数据进行处理或绘图,
% 其他nc文件可类似处理
clear all;close all;clc
%--------------------------------------------------------------------------
%--读取数据
f=netcdf('pcmdi.ipcc4.cccma_cgcm3_1.commit.run5.monthly.tas_a1_commit_5_cgcm3.1_t47_2051_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 = va