MATLAB读取文件

之前是写在OneNote里的，有自己的代码，也有参考别人的帖子总结的内容。

循环读取文件

path = ['G:\data']; 
dir_output=dir(fullfile(path,'*.dat'));
filename={dir_output.name};
for i=1:length(filename)
    file=[path,filename{i}];
    data=load(file);
end

读写xlsx文件

• 读取
分别将数值和字符存进DataNum和DataStr

name=20180331;
filename=['H:\EddyPro\output\excel\',num2str(name),'.csv'];
sheet=1;
[DataNum, DataStr] = xlsread(filename,sheet);
n=length(DataNum);

[Data,DataText,DataCell]=xlsread(‘new.xls’,‘Sheet1’)
DataCell是全部的cell cell 格式
Data带有NAN（代表字符）的全部元素 double 格式
DataText 只包含字符
[_,,DataCell]=xlsread(‘new.xls’,‘Sheet1’);

• 写入

%write
filename=[path,'MODIS_LST_LAC_T2019'];
data=[cellsta(time),num2cell(temp)];
xlswrite(filename,data);

读写二进制文件

• fopen打开文件

fid=fopen(文件名，读写方式)；
fileID = fopen(filename,permission)
fid：用于存储文件句柄值；小于0表示打开失败；大于0表示打开成功；
文件名：以字符串的形式书写（注意需要单引号）；
读写方式包括：
‘r’：只读方式打开文件（默认的方式），该文件必须已存在。
‘r+’：读写方式打开文件，打开后先读后写。该文件必须已存在。
‘w’：打开后写入数据。该文件已存在则更新；不存在则创建。
‘w+’：读写方式打开文件。先读后写。该文件已存在则更新；不存在则创建。
‘a’：在打开的文件末端添加数据。文件不存在则创建。
‘a+’：打开文件后，先读入数据再添加数据。文件不存在则创建。
注意：在打开方式后加’t’表示以文本方式打开；加‘b’表示以二进制数据方式打开；‘wb’‘wt’fopen默认以二进制方式打开；

• fread读取数据

[data,count]=fread(fid,size,数据类型)；
o output arguments
data：文件数据，作为列向量返回，如果指定sizeA参数，返回矩阵。默认double
count：读取字符数
o input arguments
fid：句柄值；
size：表示读取内容的大小；

数据类型precision

o 读取多种precision二进制文件 eg：
%# type and size in byte of the record fields
recordType = {‘double’ ‘double’ ‘int32’ ‘int8’ ‘char’};
recordLen = [8 8 4 1 1];
R = cell(1,numel(recordType));

%# read column-by-column
fid = fopen(‘file.bin’,‘rb’);
for i=1:numel(recordType)
%# seek to the first field of the first record
fseek(fid, sum(recordLen(1:i-1)), ‘bof’);

%# % read column with specified format, skipping required number of bytes
R{i} = fread(fid, Inf, [’*’ recordType{i}], sum(recordLen)-recordLen(i));
end
fclose(fid);

• fwrite写入数据
count=fwrite（fid，data，数据类型）；
fid：表示文件句柄值（打开的文件的句柄值）；
data：需要写入的数据
数据类型：double、float、int等等
count：返回所写入元素的个数（可省略）

• fclose（fid）关闭文件
在对文件进行完读写操作以后需要关闭文件，防止数据丢失或者损坏。

读写文本文件（ASCII）

• fscanf（≈fread）/txt也可以直接load（不适宜一行多个数据）
[A，COUNT]=fscanf(fid，format，size)
A用来存放读取的数据，
COUNT返回所读取的数据元素个数，可省略
fid为文件句柄
format用来控制读取的数据格式，由%加上格式符组成，常见的格式符有：d（整型）、f（浮点型）、s（字符串型）、c（字符型）等，在%与格式符之间还可以插入附加格式说明符，如数据宽度说明等。
size为可选项，决定矩阵A中数据的排列形式，它可以取下列值：N（读取N个元素到一个列向量）、inf（读取整个文件）、[M，N]（读数据到M×N的矩阵中，数据按列存放）。
• fprintf（≈fwrite）
fprintf（fid，format，A）
fid为文件句柄，指定要写入数据的文件
format是用来控制所写数据格式的格式符，与fscanf函数相同
eg：%6.3f’是指保存格式，6是指数据保存间隔，一般设置为有效数字位数+2，3是指小数点后几位，若是要换行就在后面加上\n
%d整数 %e实数：科学计数法 %f实数：小数 %g实数：自动 %s字符串
A是用来存放数据的矩阵。
eg: >> a=‘string’;

fid=fopen(‘d:\char1.txt’,‘w’);
fprintf(fid,’%s’,a);
fclose(fid);

数据定位

• fseek用于定位文件位置指针（从特定位置开始读取）
status=fseek(fid, offset, origin)
fid为文件句柄
offset表示位置指针相对移动的字节数

0 Move toward the end of the file.
= 0 Do not change position.
< 0 Move toward the beginning of the file.
origin表示位置指针移动的参照位置
‘bof’ or -1 Beginning of file
‘cof’ or 0 Current position in file
‘eof’ or 1 End of file。
若定位成功，status返回值为0，否则返回值为–1。
• ftell 返回文件指针的当前位置
position=ftell (fid)
返回值为从文件开始到指针当前位置的字节数。
若返回值为–1表示获取文件当前位置失败。

nc数据读写

%% 显示结构
% ncdisp(ncFilePath);%显示nc文件的所有结构，以便大概了解里面的内容
% ncdisp(ncFilePath,‘evap’);%显示指定变量的内容，注意一定要是变量variables才可以
% ncdisp(ncFilePath,’/’,‘min’);%简单显示结构以及定义
% ncdisp(ncFilePath,’/’,‘full’);%全部显示所有结构和定义信息

%% 读取变量值
% ncid = netcdf.open(ncFilePath,‘NOWRITE’); %打开nc文件返回索引ID
% [ndims,nvars,ngglobalatts,unlimdimid] = netcdf.inq(ncid);%获取维数，变量数，全局属性数量（此处ncid从0开始为第一个变量）
% [varname,xtype,dimids,natts] = netcdf.inqVar(ncid,0); %根据变量索引号获取变量的名称
%var=ncread(ncFilePath,‘varname’);%读取变量，ncread读出结果为真实值，即自动将fillvalue记为nan，vardata*scale+offset

%% 显示数据
% pcolor(lat,lon,evap1);
% [x,y]=meshgrid(lon,lat);%根据经纬度信息产生格网，284列（经度），164列（纬度）
% phandle=pcolor(x,y,evap1’);%显示一个矩阵，其中x,y,evap1的行列数必须一致。类似surface函数
% colorbar

%%time
在处理气象nc格式数据时，时间变量往往不是已YYYYMMDD格式给出日期，而是距某一天（如1900年1月1日）XXX小时或XXX日。
可由以下代码将nc的时间变量转为日期：
t0 = datetime(1900,1,1);
date_yyymmdd = t0 + double(time(😃)/24; %time为距1900年1月1日00时的小时数一维数组
%转回nc格式
time=time_matlab;
time_nc=(time-datenum(‘1900-01-01 00:00:00.0’))*24;

eg:

dat_filename=[path,'MOD11A1.006_1km_aid0001.nc'];
ncid = netcdf.open(dat_filename,'NOWRITE');
%ncdisp(dat_filename) %查看nc文件结构
[~,nvars,~,~] = netcdf.inq(ncid);%获取变量数
for i=1:nvars %vars
    varname{i}= netcdf.inqVar(ncid,i-1);%id从0开始
    eval([varname{i},'=ncread(dat_filename,varname{i});']);
end
netcdf.close(ncid)
clear i 

%%写入nc
% ---------------------------- DEFINE THE FILE --------------------------- %
ncid = netcdf.create(‘f:\ncl\ecco_ekman_mon_mean1.nc’,‘CLOBBER’);
%创建一个存放数据的nc文件

%-----------------------------define dimension-----------------------------%
%因为有个变量是三维的，所以定义了三个维度，其它一维二维可视具体情况定义
dimidx = netcdf.defDim(ncid,‘lat’,242);
dimidy = netcdf.defDim(ncid,‘dep’,50);
dimidz = netcdf.defDim(ncid,‘time’,12);

%----------------------------define new variables---------------------------------%
varid = netcdf.defVar(ncid,‘strf’,‘double’,[dimidx dimidy dimidz]);
varid2 = netcdf.defVar(ncid,‘lat’,‘double’,[dimidx]);
varid3 = netcdf.defVar(ncid,‘dep’,‘double’,[dimidy]);

%---------------------------define attributes of the new variables--------------%
netcdf.putAtt(ncid,varid,‘units’,‘Sv’); %单位信息和long_name，其它的信息可依此定义
netcdf.putAtt(ncid,varid2,‘units’,‘degress_north’);
netcdf.putAtt(ncid,varid3,‘units’,‘m’);

netcdf.putAtt(ncid,varid,'long_name','The Meridional streamfunction ');
netcdf.putAtt(ncid,varid2,'long_name','Latitude');
netcdf.putAtt(ncid,varid3,'long_name','Depth');
netcdf.endDef(ncid);

%--------------------------给新变量赋值-------------------------------------------%
netcdf.putVar(ncid,varid,fail_mon_mean);
netcdf.putVar(ncid,varid2,lat);netcdf.putVar(ncid,varid3,d);
netcdf.close(ncid);