注意: 此文章主要针对本人在工作中遇到的问题进行记录和经验上的总结,难免存在错误和遗漏,如有任何问题,请留言,我会尽快回复。
最近用Matlab绘制图像时,发现用pcolor和imagesc画图区别挺大的,一不小心就容易弄错,因此记录一下,免得以后再出问题。
废话不多说,先上代码。
% 创建一个5x5的矩阵
Data = [1,2,3,4,5;
6,7,8,9,10;
11,12,NaN,14,15;
16,17,18,19,20;
21,22,23,24,NaN];
% 给5x5的矩阵赋予X,Y坐标
X = [1,2,3,4,5];
Y = [2,4,6,8,10];
% 分别用imagesc和pcolor绘图,并查看区别
figure(1);
subplot(1, 2, 1);
imagesc(X,Y,Data);
colorbar;
title('imagesc 绘图');
axis equal tight;
grid on;
subplot(1, 2, 2);
pcolor(X,Y,Data);
colorbar;
title('pcolor 绘图');
% 调整坐标轴,防止拉伸,保持固定比例
axis equal tight;
grid on;
运行以上代码,可以得到如下结果:
根据结果可发现:
pcolor绘图会漏数据
这是因为pcolor绘图是以像素左下角为坐标基准点绘制图像的,比如,我们的X轴坐标为1-5,pcolorX轴画到5就停止了,而5所代表的像素并不会显现,Y轴同理。因此,用pcolor画图会漏掉最右侧和最顶端的数据。相反,imagesc是以像素中心点为基准绘制图像的,因此不会漏数据。而且,一般像GLDAS等大数据在提供lat和lon的时候,都是提供的像素中心点的坐标,而不是像素左下角的坐标。最重要的是,一般在研究中,我们用Matlab画图还会添加shapefile等文件凸显研究区域,用pcolor绘制时,可能会出现像素对不准的问题(除非手动调教坐标数据),所以,一般我画图,都用imagesc命令。
imagesc绘图时,Y轴默认是上小下大
当然,用imagesc绘图也会出现一些问题。例如,imagesc绘图时,Y轴的左边系统默认上小下大。这对于我们画纬度(北半球为正值,位于坐标轴上方,南半球为负,位于坐标轴下方)时,会出现图像上下颠倒的情况。而pcolor绘制图像时,Y轴默认设置跟imagesc是反的,所以不会出现图像颠倒的情况。
imagesc的Y轴颠问题可以用下面一排代码解决:
h = imagesc(X, Y, Data);
set(gca, 'YDir', 'normal'); % 设置 Y 轴方向为正常,即下方为负,上方为正
imagesc绘图时,NaN数据无法自动显示透明色
如上图所示,imagesc在绘制图像时,遇到数据为空集(NaN)的情况时,仍然会为像素上色。而pcolor绘制图像时,则会将数据为空集(NaN)的像素设置为透明色。
imagesc的空集(NaN)上色问题可以用以下代码解决:
alphaData = ~isnan(Data); % 查找无数据的像素
set(h, 'AlphaData', alphaData); %将图中无数据的区域设为透明颜色;