在 Matlab 中绘制类似于使用 Gnuplot 获得的曲面图？

Question

更新：

链接到我的测试数据集

我有一些这样的数据：

10.0 11.0 0.5
10.0 12.0 0.5
10.0 14.5 0.5
10.0 16.0 0.5
10.0 18.5 0.5
10.0 19.0 0.5
10.0 19.5 0.5
10.0 20.0 0.5

使用Gnuplot，我可以绘制下图：

使用我编写的以下脚本：

set data style lines
set surface
set dgrid3d 30,30
set hidden3d
set ticslevel 0.8
set isosample 40,40

set view 60, 30, 1, 1
set contour base

splot "dataset" using 2:1:($3) notitle

我正在尝试对此图进行一些分析，所以想要将其输入 Matlab，但我不知道如何绘制它。我尝试将数据集导入 Matlab，这就是我所能达到的程度。我正在加载文件并使用以下 Matlab 语句进行绘图：

load data
rangeY = floor(min(data(:,2))):.5:ceil(max(data(:,2)))
rangeX = floor(min(data(:,1))):.5:ceil(max(data(:,1)))
[X,Y] = meshgrid(rangeX,rangeY)
Z = griddata(data(:,1),data(:,2),data(:,3),X,Y, 'cubic')
surf(X,Y,Z)

图片看起来与我截然不同（轴范围、形状等）有人可以告诉我如何获得类似于我通过 Gnuplot 获得的绘图吗？

Answer 1

surf 无法立即工作的原因是数据非常嘈杂，如果您调用，您会看到

plot3(data(:,1),data(:,2),data(:,3),'.')

因此，您首先需要平滑数据（我猜这就是gnuplot 函数也可以）。以下是平滑数据的一种方法

%# set bin edges (maybe you want to do 5:10:205)
xEdge = 10:10:200;
yEdge = xEdge;

%# bin the data
xId = sum(bsxfun(@ge,data(:,1),xEdge(2:end-1)),2)+1;
yId = sum(bsxfun(@ge,data(:,2),yEdge(2:end-1)),2)+1;

%# average z within every bin
smoothData = accumarray([xId,yId],data(:,3),[length(xCenter),length(yCenter)],@mean);

%# remove bins with no data
smoothData(smoothData==0) = NaN;

%# bin centers
xCenter = (xEdge(1:end-1)+xEdge(2:end))/2;
yCenter = xCenter;

%# plot
meshc(xCenter,yCenter,smoothData)

Answer 2

我很惊讶你的代码不起作用。我在我的一个数据集上尝试了一下，效果很好。作为第一个答案，我建议您发布数据集的链接（如果不敏感），并尝试其他插值方法（从默认的线性插值方法开始）。

然而，我注意到 griddata 在 XY 域的边界产生了相当多的 NaN，尽管我的数据集非常密集并且在该域中分布良好。我检查了文档，其中提到：

有时，griddata 可能会以 NaN 形式返回数据凸包上或非常靠近凸包的点。这是因为计算中的舍入有时会导致难以确定边界附近的点是否在凸包内。

我不认为这是问题的原因，但我仍然怀疑这个函数在非常稀疏的数据集上表现不佳。您应该尝试使用 TriScatteredInterp 来代替，这无论如何都会更快。