Regarding reading diffrent files

2011 9 19
1 个回答
7 次查看（30 天）

0 个投票

I have some .nc files for 13 years (13 files) e.g GTE_bb_CH_1997.nc,GTE_bb_CH_1998.nc,...., GTE_bb_CH_2009.nc I want to read and plot all files same time. So how I can write loop so I can read separate file for every year ?

请先登录，再回答此问题。

回答（1 个）

Fangjun Jiang
Fangjun Jiang 2011-9-19

0 个投票

Files=dir('*.nc');
for k=1:length(Files)
FileName=Files(k).name;
%load file
end

7 个评论
显示 5更早的评论 隐藏 5更早的评论

Uday
Uday 2011-9-19
how to load file here?
Fangjun Jiang
Fangjun Jiang 2011-9-19
That depends. What is .nc file? There is no built-in function to load .nc file directly. If it's a text file, you have to understand its format and then load it use uiimport(),importdata(), textscan() or fopen(), fread().
Uday
Uday 2011-9-19
netcdf is special data file, generally its stores huge amount of data.
I have used following codes to read files , it works for single file , but I would like to run all once
path=('c:/GEFD/');
dir_list=dir('*.nc');
for i=1:length(dir_list)
FileName=dir_list(i).name;
ncid=netcdf.open(FileName,'NOWRITE');
latid=netcdf.inqVarID(ncid,'lat');
latitude=double(netcdf.getVar(ncid,latid));
lonid=netcdf.inqVarID(ncid,'lon');
long=double(netcdf.getVar(ncid,lonid));
longitude=rem((long+180),360)-180;% conversion from 0-360 to -180 to 180
emiss_bbid=netcdf.inqVarID(ncid,'emiss_bb');
bb=double(netcdf.getVar(ncid,emiss_bbid));
Uday
Uday 2011-9-19
but these codes are not reading all files
Fangjun Jiang
Fangjun Jiang 2011-9-19
Okay, Thank you! I learned netcdf() today! It is reading all files. But you are over-writing the data every time, right? You need to add an "end" line to complete the for-loop. To store data for every file, you need to declare an array before the for-loop. Something like latitude=zeros(length(dir_list),1). Then, inside the for-loop, use latitude(i)=double(netcdf.getVar(ncid,latid));
Uday
Uday 2011-9-19
its not over writing
dir_list=dir('*.nc');
for i=1:length(dir_list)
FileName=dir_list(i).name;
ncid=netcdf.open(FileName,'NOWRITE');
the above para does not read the data files those I have in that folder ( 13 different year files).
Sorry I just forgot to wrote end in the previous codes.
Fangjun Jiang
Fangjun Jiang 2011-9-19
Try this:
dir_list=dir('*.nc');
for i=1:length(dir_list)
FileName=dir_list(i).name
end
Do you see 13 file names appear in the Command Window?

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 HDF5 Files 的更多信息

标签

提问:

Uday
Uday
2011-9-19

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by