新手快被逼疯了，matlab一直在报错

zyw199818

close all;
clc;
clear;
cd ('E:\TUXIANG\') ; %待处理的数据文件所在文件夹
Allname=struct2cell (dir) ; %得到上述文件夹下的所有文件名
[m, n]=size (Allname) ;
data= [];

for i=1:n %前两个不是文件名
name=Allname{1,1}
filename=strcat(name);
ncidsst=netcdf.open('TU.nc','NOWRITE'); %打开HDF文件
lonID=netcdf. inqVarID(ncidsst,'lon') ;%返回经度变量ID
lon=netcdf.getVar(ncidsst, lonID); %读入经度值
latID=netcdf. inqVarID(ncidsst,'lat');%返回维度变量ID
lat=netcdf. getVar (ncidsst, latID) ; %读入纬度值
sstID=netcdf. inqVarID(ncidsst,'sst') ;%返回sst变量ID
sst=netcdf.getVar(ncidsst,sstID); %读入sst值
netcdf. close (ncidsst) ; %关闭HDF文件

longl=(0+0)*4;
long2=(10+0)*4;
lati1=(50+90)*4;
lati2=(60+90)*4;
就是这后面一句一直显示说不是整数或者逻辑型变量
RSsst=sst(longl: long2,lati1:lati2)  %截取sst
RSsst(RSsst==-999)=0;%陆地部分为零
RSsstnew=double(RSsst)*0.01;%变成正确的温度
RSsstnew=imrotate(RSsstnew,90) ;%把图像旋转摆正
%imshow(RSsstnew) ;%显示看下是否颠倒
%下步用图像处理滑动邻域操作
fun=inline('mean (mean (x))') ; %做矩阵平均
%B=nlfilter(RSsstnew, [2 2], fun);
B=blkproc (RSsstnew,[4 4],fun) ; %分块矩阵平均
Bl=B(:)';%把矩阵变为一行， 是按照一列一列的方式变成一行
data(i,:)=Bl;%有多少天都放在一起
end
xlswrirte ('D: \englangweather.xls', data) ;

阳光正好

long1等于0肯定是不行的，至少等于1

zyw199818

我去试试，谢谢您的帮助

zyw199818

大佬然后想问您一下，我用什么函数画这个nc文件的图?有高大上一点的嘛

zyw199818

阳光正好 发表于 2020-2-14 20:13
long1等于0肯定是不行的，至少等于1

大佬然后想问您一下，我用什么函数画这个nc文件的图?有高大上一点的嘛