发个自己写的逐日气温预报检验程序

schliezer

短期岗，每天跟气温较劲，跟中央台斗智斗勇，编个程序，看看中央台（excel里的sevpT）以及本省（XXT）对于实况（ObsT）的偏差情况吧。

其中，中央台(\micaps\cityfcst\)和本省的精细化预报(\micaps\....)以及最低温实况（\micaps\Surface\Tmin-p）文件的路径是本省的网络盘，最高温实况由于Surface\Tmax-p缺了很多站，是自己用业务平台生成的xls文件，其数据格式在包里已有示例。‘xx’均表示台站名，‘DDDDD’均表示台站号，2016_T_test.xls是结果，xlsappend.m是matlab官网上的xls写入函数，要和其余函数以及主m放一个路径下使用。

余辉

mofangbao

多谢分享，楼主可以更进一步直接出图显示

schliezer

mofangbao 发表于 2016-6-1 11:56
多谢分享，楼主可以更进一步直接出图显示

好的谢谢，我在考虑下如何出图

schliezer

mofangbao 发表于 2016-6-1 11:56
多谢分享，楼主可以更进一步直接出图显示

加入画图的模块：
YLabel=c(16:27,3);
XLabel=cell(1,34);
XLabel=num2cell([1:34]);
% YLabel=cell(1,12);
% YLabel=num2cell([1:34]);
Z=cell2mat(c(16:end-1,4:end));
x0=[1:length(XLabel)];
y0=[1:length(YLabel)];
[xx,yy]=meshgrid(x0,y0);
clf
set(gcf,'Position',[100  50 1600 800], 'color','w')
pcolor(xx,yy,Z);
h=colorbar('Eastoutside');
t = get(h,'YTickLabel');
t = strcat(t,'℃');
set(h,'YTickLabel',t);
set(gca,'YTickLabel',YLabel');
set(gca,'XTick',x0);
set(gca,'XTickLabel',XLabel);
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 14 5])
print(gcf,'-djpeg',[INPUTDIR_Obs,st,'T_test_Mosaic.jpg'],'-r300');

图片中，横轴数字代表站点（就不具体显示是哪些站了哈），马赛克代表预报对实况相差的度数：

realy

感谢楼主分享~

sunflower

给楼主点赞

xopowomimi

看看，呵呵，还要贡献啊