【原创+整合版】matlab批量读取nc文件

an0621yan

感激万分，谢谢分享！

kongfeng0824

vendya 发表于 2014-9-30 23:19
楼主，请问，matlab读取nc数据的时候是怎么存储的？行是经度，列是维度吗？当我用ncread 读取出数据的时候 ...

需要对数据转置

njhanjing

谢谢楼主分享，解决一个大难题啊，哈哈哈

vendya

楼主，请问，matlab读取nc数据的时候是怎么存储的？行是经度，列是维度吗？当我用ncread 读取出数据的时候，画出数据的时候，发现数据是倒着的，我的代码如下，不知可否帮忙看一下？ 谢谢！
%define region
%%read lat and lon
lat=ncread('/Users/yhe/Documents/BMB/tasmax_day_CESM1-CAM5_rcp85_r1i1p1_20060101-20401231.nc','lat');
lon=ncread('/Users/yhe/Documents/BMB/tasmax_day_CESM1-CAM5_rcp85_r1i1p1_20060101-20401231.nc','lon');
%time=ncread('/Users/yhe/Documents/BMB/tasmax_day_CESM1-CAM5_rcp85_r1i1p1_20060101-20401231.nc','time');
%%find the id for region
i=find(lon>=277&lon<=283);
j=find(lat>=37&lat<=41);
%%read data for region
tasmax_wv=ncread('/Users/yhe/Documents/BMB/tasmax_day_CESM1-CAM5_rcp85_r1i1p1_20060101-20401231.nc','tasmax',[i(1),j(1),1],[i(end)-i(1)+1,j(end)-j(1)+1,365]);
%%display the map
latlim = [lat(j(1)) lat(j(end))];
lonlim = [lon(i(1))-360 lon(i(end))-360];
ax=usamap(latlim,lonlim);         
setm(gca,'MapProjection','lambert')
hold on
%Contour
[c,h]=contourfm(lat(j),lon(i),tasmax_wv(:,:,1));
clabelm(c,h)
%colorbar
getm(ax)
conus=load('conus')
plotm(conus.statelat,conus.statelon,'k')
alpha(0)

日不不不落

qingdu

楼主大神，支持

409696585

楼主，为什么我在执行第二条命令时就出错了呢？就是ncdisp('precip.mon.1981-2010.ltm.nc');这一条

笨小孩0410

好东西！谢谢楼主分享！

烟熏在你眼

正好今天在学习nc文件，好东西~~~

wangpang

一水天下 发表于 2014-8-17 13:40
先看下precip是几维的，然后在读取的时候增加一维s，PrecipData(:,:,:,s) = ncread(filename,'precip');
...

正解，多谢