% 读取校准后土壤体积含水量数据
load('vwc_up.mat');
% 土壤体积含水量数据提取
date=datetime(vwc_up.date_up(24:2689,:)); %提取日期
water_up=vwc_up.calibration_up(24:2689,:); % 选择一定时间段的数据
%绘图代码
date=datetime(vwc_up.date_up(24:2689,:)); %提取日期
Soil_depth=[10 20 40 100 150 200]; %土壤深度
date_1=datenum(date);
[date_C,SD]=meshgrid(date_1,Soil_depth);
contourf(date_C,SD,(water_up.*100)'); %绘制颜色图
c=flipud(jet);colormap(c);%热图颜色
set(gca,'YDir','reverse'); %设定y轴朝下
set(gca,'xaxislocation','top'); %坐标轴位置设定上部显示
datetick('x','mmm-dd') %将坐标轴设置为日期格式
set(gca,'FontSize',10,'fontname','Times','FontWeight','bold');%设置y轴刻度标签字体
h=colorbar;
% h=colorbar('position',[0.78 0.2 0.05 0.5]); %colorbar([a b c d]),其中a，b为坐标，c为宽度，d为高度
set(h,'ylim',[0,8]); %设置colorbar范围
xlim([738693, 7.3880e+05]); %定义横坐标的范围
set(gca,'looseInset');
caxis([0 8]);%设置当前绘图区的颜色范围
% xlabel('GD_1_3','fontsize',12,'fontname','Times','FontWeight','bold')%设置x轴字体
ylabel('Soil depth (cm)','fontsize',12,'fontname','Times','FontWeight','bold');%设置y轴字体
title('(a1) VWC(%): GW13','fontsize',12,'fontname','Times','FontWeight','bold');
ax = gca;
ax.TitleHorizontalAlignment = 'left';%标题显示在左边
% set(gca,'TickDir','out','xticklabel',[]); %隐藏x轴刻度线标签
结果展示：
转载本文请联系原作者获取授权，同时请注明本文来自柏延文科学网博客。
链接地址：https://m.sciencenet.cn/blog-3402238-1382582.html?mobile=1
收藏
