- 积分
- 121
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-5-3
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
自己按照kongfeng0824大大提供的MATLAB读取nc程序,自己改了一下,无脑的自己改了一下。求问,问题出在哪里?
原地址:http://bbs.06climate.com/forum.php?mod=viewthread&tid=11870
程序:
clear all;close all;clc
datadir='F:\nc\'; %指定数据所在文件夹
ncid = netcdf.open('F:\nc\sst.ltm.1981-2010.nc','NOWRITE'); %打开nc文件
SstData = ncread('sst.ltm.1981-2010.nc,''sst'); %读入变量sst
TimeData = ncread('sst.ltm.1981-2010.nc.ltm.nc','time'); %读入变量time
LonData = ncread('sst.ltm.1981-2010.nc.ltm.nc','lon'); %读入变量lon
LatData = ncread('sst.ltm.1981-2010.nc.ltm.nc','lat'); %读入变量lat
Valid_yr_countData = ncread('sst.ltm.1981-2010.nc.nc','valid_yr_count'); %读入变量validprecip_yr_count
ClimatologyData = ncread('sst.ltm.1981-2010.nc','climatology_bounds'); %读入变量climatology_bounds
[X, Y] = meshgrid(LatData,LonData);
contourf(Y-180.0,X,PrecipData(:,:,9)); %画9月份等值面图
shading flat; %去掉等值线
colorbar('SouthOutside','Position',[0.142,0.03,0.75,0.04]); %添加颜色条 [左右,上下,长,宽]
load coast %加载全球海岸线,但是不显示出来
geoshow(lat,long); %显示出海岸线,lat和long是coast的属性
hold on;
set(gca,'LineWidth',1,'FontSize',10,'Ylim',[-90,90],'Xlim',[-180,180],'Position',[0.142,0.09,0.75,0.84]...
,'XTick',[-180:60:180],'XTicklabel',{'-180W','-120W','-60W','0','60E','120E','180E'}...
,'YTick',[-90:30:90],'YTicklabel',{'-90S','-60S','-30S','0','30N','60N','90N'}); %添加经纬度信息
hold off;
netcdf.close(ncid); %关闭nc文件
结果:
Error using internal.matlab.imagesci.nc/openToRead (line 1243)
Could not open sst.ltm.1981-2010.nc,'sst for reading
Error in internal.matlab.imagesci.nc (line 122)
this.openToRead();
Error in ncread (line 54)
ncObj = internal.matlab.imagesci.nc(ncFile);
出来的结果是这样的,求众位大大指教,问题出在哪里呢?
|
-
-
|