% 读取校准后土壤体积含水量数据

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轴刻度线标签

结果展示：