MATLAB读nc文件及创建写入nc

春风得意马蹄急

嗯嗯，初学者，分享了

春风得意马蹄急

嗯嗯，初学者，分享了

春风得意马蹄急

嗯嗯，初学者，分享了

LemonTree007

希望对我的毕设有用

LemonTree007

楼主你好，我按照你的分享结合自己的修改了一下，出现了错误，但我不会改，求指点。我是菜鸟，新新手。麻烦了。
代码如下，出错显示随后：

1. clear all
   
2. clc
   
3. %定义输入目录
   
4. inputdir = 'D:/matlabht/pdsi0505';
   
5. %定义输出目录
   
6. outputdir = 'D:/matlabht/pdsi';
   
7. %获取目录下所有nc文件的信息
   
8. inputfiles = dir([inputdir , '/', '*.nc']);
   
9. outfile = (['pdsi199101_201412.nc']);%定义输出的nc文件的文件名
   
10. outdirfile = [outputdir '/' outfile];%输出路径
    
11. outid = netcdf.create(outdirfile,'NC_NOCLOBBER');%创建nc文件
    
12. latdimID = netcdf.defDim(outid,'lat',360);%创建维度及维数
    
13. londimID = netcdf.defDim(outid,'lon',720);%同理
    
14. timedimID = netcdf.defDim(outid,'time',1380);%同
    
15. 
    
16. latid = netcdf.defVar(outid,'lat','double',latdimID);%定义变量，变量名等。
    
17. lonid = netcdf.defVar(outid,'lon','double',londimID);
    
18. timeid = netcdf.defVar(outid,'time','int',timedimID);
    
19. pdsiid = netcdf.defVar(outid,'pdsi','double',[londimID,latdimID,timedimID]);                 
    
20. netcdf.endDef(outid);%定义结束
    
21. for ifiles = 1:length(inputfiles) %length（inputfiles）是获取inputfiles的长度，就是nc文件的数目
    
22.     infile = inputfiles(ifiles).name;%获取nc文件名
    
23.     varName = infile(11:14);%把文件名的pdsi赋值给varname
    
24.     indirfile = [inputdir '/' inputfiles(ifiles).name];%获取nc文件的完整路径
    
25.     time = ncread(indirfile,'time');%读取时间变量
    
26.     Ti=size(time);%获取时间变量的尺度
    
27.     lon = ncread(indirfile,'lon');%读取经度
    
28.     lat = ncread(indirfile,'lat');%读取纬度
    
29.     pdsi = ncread(indirfile,'pdsi');
    
30.     pdsi = double(pdsi);
    
31. 
    
32.     netcdf.putVar(outid,latid,lat);%将值赋给变量
    
33.     netcdf.putVar(outid,lonid,lon);
    
34.     netcdf.putVar(outid,timeid,time);
    
35.     netcdf.putVar(outid,pdsiid,pdsi);
    
36.     netcdf.reDef(outid);%赋值结束
    
37. 
    
38.     netcdf.putAtt(outid,latid,'long_name','latitude');%对变量的说明
    
39.     netcdf.putAtt(outid,latid,'units','degrees_north');
    
40.     netcdf.putAtt(outid,latid,'standard_name','latitude');
    
41. 
    
42.     netcdf.putAtt(outid,lonid,'long_name','longitude');
    
43.     netcdf.putAtt(outid,lonid,'units','degrees_east');
    
44.     netcdf.putAtt(outid,lonid,'standard_name','longitude');
    
45. 
    
46.     netcdf.putAtt(outid,timeid,'long_name','time');
    
47.     netcdf.putAtt(outid,timeid,'units','hours since 1900-01-01 00:00:00');
    
48.     netcdf.putAtt(outid,timeid,'standard_name','time');
    
49.    
    
50.     netcdf.putAtt(outid,shumid,'long_name','Air specific humidity');
    
51.     netcdf.putAtt(outid,shumid,'units','kg kg-1');
    
52.     netcdf.putAtt(outid,shumid,'_FillValue','-9.99e+33');
    
53.     netcdf.putAtt(outid,shumid,'missing_value','-9.99e+33');
    
54.     netcdf.close(outid);
    
55. end

复制代码

错误命令行窗口
错误使用 netcdflib
无法使用指定的文件名创建文件。文件名可能包含不受支持的字符。

出错 netcdf.create (line 53)
        ncid = netcdflib('create', filename, mode);

出错 humidity (line 11)
outid = netcdf.create(outdirfile,'NC_NOCLOBBER');%创建nc文件

LemonTree007

之前的问的问题不够详细和准确问题，我自己发了帖子求助，求指教！
详见http://bbs.06climate.com/forum.p ... id=60398&extra=

woshishui321

实用性比较强大的

鱼啊鱼

谢谢分享，果断下载啦~

1314520tao

6“““好：：不”””

kingjames

谢谢楼主  十分感谢分享  666   雪中送炭