Matlab中如何基于非均匀散点数据绘制带插值的热力图？

问题1：生成10000×10000异常图像的原因及解决 #

你用meshgrid(D, all_depth)直接基于原始数据点生成网格，相当于把10000个日期点和10000个深度点做全组合，自然会得到10000×10000的超大矩阵，这完全没必要。正确的做法是基于原始数据的范围生成均匀采样的网格，既节省内存，又能得到平滑的热力图。

代码示例：

% 先将日期转为Excel数值格式
D = convertTo(all_dates, 'excel');

% 获取原始数据的范围
date_range = [min(D), max(D)];
depth_range = [min(all_depth), max(all_depth)];

% 自定义采样密度，比如x轴（时间）取100个点，y轴（深度）取100个点，可按需调整
num_time_points = 100;
num_depth_points = 100;

% 生成均匀采样的网格
Xi = linspace(date_range(1), date_range(2), num_time_points);
Yi = linspace(depth_range(1), depth_range(2), num_depth_points);
[Xi_grid, Yi_grid] = meshgrid(Xi, Yi);

% 执行插值
F = scatteredInterpolant(D, all_depth, example_data);
Ci = F(Xi_grid, Yi_grid);

问题2：还原日期轴显示 #

把插值后的数值型日期转回Matlab的datetime格式，直接传给imagesc即可自动识别日期轴，再配合刻度调整优化显示：

% 将数值日期转回datetime格式
Xi_datetime = datetime(Xi, 'ConvertFrom', 'excel');

% 绘制热力图
imagesc(Xi_datetime, Yi, Ci);
colorbar; % 显示Z值对应的颜色条

% 设置轴标签与标题
xlabel('时间');
ylabel('深度');
title('参数Z随时间-深度的热力图');

% 优化x轴显示：旋转标签避免重叠
xtickangle(45);

% 可选：反转y轴，让深度从上到下递增（符合常规观测习惯）
set(gca, 'YDir', 'reverse');

额外优化建议 #

• 插值方法：scatteredInterpolant默认用线性插值，你可以根据数据特性指定其他方法，比如F = scatteredInterpolant(D, all_depth, example_data, 'nearest');（最近邻插值）或'cubic'（三次插值）。
• 采样密度：如果觉得热力图不够平滑，可以增大num_time_points和num_depth_points的数值，但注意不要过大导致内存占用过高。

内容的提问来源于stack exchange，提问作者Donald Martocello