||
% 读取校准后土壤体积含水量数据
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轴刻度线标签
结果展示:
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 09:27
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社