关于批量作图~

小婕婕

在这里和大家分享下批量作图，当然我这里的批量也不算很多，由于数据的缘故，每周一个nc文件数据，目前我按一个月一个月来做，就拿南海sla为例，
for i=0:3
f=20030108+7*i;
f1=num2str(f);
fname=['dt_ref_global_merged_msla_h_',num2str(f),'_',num2str(f),'_20100503.nc'];
ncload(fname);
lon=NbLongitudes(307:381);%102-126.6667
lat=NbLatitudes(456:530);%-0.7576-23.249
sla=Grid_0001(307:381,456:530);
sla(sla>=10000)=nan;
   m_proj('lambert','lon',[102 127],'lat',[-1 24]);
   m_contourf(lon,lat,sla',20);
    hold on
       hc=colorbar;
     m_coast('patch',[0.6825 0.6825 0.6825],'edgecolor','b');
     m_grid('box','fancy','tickdir','out');
shading flat
print('-dtiff','-r200',f1);
end
这样就可以做出一个月的南海sla图像，大家可以回去试试，还有一个小问题就是虽然这样图可以做出来，但是会有很多的warning，这是为什么呢？

Aires

你把shading flat去掉，用pcolor才用shading flat，contour和contourf不用，你去掉后看还有没有waring

windprofiler

批量读入的话，
可以用filelist=dir（‘文件地址’）得到目录，
然后fopen（filelist(i).name）循环读入，
这样感觉更简单点

Aires

确实dir是比较方便的，婕婕可以随便造个数据试一下

小婕婕

Aires 发表于 2012-11-5 19:01

                               
登录/注册后可看大图

确实dir是比较方便的，婕婕可以随便造个数据试一下

嗷嗷，是的，我又弄了dir的，就是把所有数据文件放在一个文件夹里就行，然后就可以顺序读取啦，这次想弄多少个就能弄多少个，不会限制只能一个月的了，O(∩_∩)O~

Aires

小婕婕 发表于 2012-11-13 16:26

                               
登录/注册后可看大图

嗷嗷，是的，我又弄了dir的，就是把所有数据文件放在一个文件夹里就行，然后就可以顺序读取啦，这次想弄多 ...

嗯，一般批量处理都这么做

zjgi756

楼主可以再发下用dir后的程序怎么写么？谢谢了

gtouc

谢谢楼主 很受用

骆驼

很不错哈, 学习

maiyatang

顶！d=====(￣▽￣*)b，学习学习