求助，这样能不能读取网格中的数据

喃喃楠楠喃喃

clc;
clear;
datadir='F:\海洋数据\Jason2\';%指定批量数据所在的文件夹
filelist=dir([datadir,'*.nc']);%指定批量数据的类型
a=filelist(1).name;%查看要读取的文件的编号
b=filelist(2).name;
k=length(filelist);
for s1=1:k
     filename=[datadir,filelist(s1).name];
     ncid=netcdf.open(filename,'NC_NOWRITE');
     ncdisp('F:\海洋数据\Jason2\JA2_GPN_2PdP219_010_20140613_002440_20140613_012053.nc');
     ncid=netcdf.open('F:\海洋数据\Jason2\JA2_GPN_2PdP219_010_20140613_002440_20140613_012053.nc','NOWRITE');
     ncdisp('F:\海洋数据\Jason2\JA2_GPN_2PdP219_010_20140613_002440_20140613_012053.nc');
     LonData{s1}=ncread(filename,'lon');
     LatData{s1}=ncread(filename,'lat');
     Swh_kuData{s1}=ncread(filename,'swh_ku');
     B=find(LonData{s1})>=180;
     LonData{s1}(B)=LonData{s1}(B)-360;
     C{s1}=[LonData{s1},LatData{s1},Swh_kuData{s1}];
     netcdf.close(ncid);
     plot(LonData{s1},LatData{s1},'b.');
     hold on
end     
D=reshape(C,s1,1);
E=cell2mat(D);%将D中的所有数据汇总到一个表格中
id=find( E(:,3)>11 | isnan(E(:,3))==1);%选出swh_ku中大于11 和出现nan数据的项
E(id,:)=[];
n2=length(E);
figure(2)
plot(E(:,1),E(:,2),'g.')
grid on
set(gca,'xtick',[-180:3:180],'ytick',[-90:3:90]);
for i=1:60
    ymin(i)=3*i-90-3;ymax(i)=3*i-90;%y分度值为3
    for j=1:120
    xmin(j)=3*j-180-3;xmax(j)=3*j-180;%x的分度值为3
    x0{i,j}=xmin(j)+1.5;y0{i,j}=ymin(i)+1.5;
    ind{i,j}=find(E(:,2)>=ymin(i) & E(:,2)<=ymax(i) & E(:,1)>=xmin(j) & E(:,1)<=xmax(j));%找到的就是在ymin<=E(:,2）<=ymax和xmin<=E(:,1）<=xmax之间的点
    F{i,j}=E(ind{i,j},:);
    x{i,j}=F{i,j}(:,1);y{i,j}=F{i,j}(:,2);z{i,j}=F{i,j}(:,3);
    H{i,j}=var(z{i,j});
    end
end
我想要的是将网格3*3进行划分，然后读取网格化之后的数据

二爷名声在外

1.能不能，自己试一试就知道了，只看代码，我也不知道能不能，只有matlab知道
2.看楼主的程序，建议还是找师姐指导一下吧，我感觉里面元胞和数组之间的关系很乱，最好还是先了解熟悉一下二者的特性再进行使用。

喃喃楠楠喃喃

二爷名声在外 发表于 2015-11-2 19:23
1.能不能，自己试一试就知道了，只看代码，我也不知道能不能，只有matlab知道
2.看楼主的程序，建议还是找 ...

我试过，就是不知道对不对。嗯呢，我还是再看看元胞与数组之间的关系

Carol_CHH

你好，我最近在做一个类似的呢，目前是读取了nc文件存储在元胞里啦，但是里面有一些异常数据需要处理，然后还需要进行计算、插值。这种情况也是需要把所有的矩阵合并为一个表格吗？