非均匀间隔的数据点的二维等高线图

66(30天)

显示旧的评论

一个 2015年5月26日

2
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/218758-2-d-contour-plot-for-non-uniformly-spaced-data-points

评论道: 约翰D 'Errico 2020年7月1日

我从COMSOL数组数据非均匀网格:x, y节点位置和u为每个节点,节点细化的一个区域。我想创建一个二维等高线图的数据在MatLab在x - y图和颜色代表u。发现似乎涉及到的所有信息通过“meshgrid”创建一个矩形网格,然后让它插入数据到这个网。我在想如果有一个简单的方法,我可以画出原始数据所以我不排除任何价值观?

任何帮助深表感谢。

2的评论
显示1年长的评论藏1年长的评论

员工切赫Michalek 2020年7月1日

你好,

我也有同样的问题。我得到的数据通过LiveLink Comsol Matlab。我想创建一个contourf情节从x, y和U,我一个一个的U值x和y的价值。然而,你是行向量1 x12792双和x和y的值分布不均一,因为它们是三角形网格节点的坐标。我可以生成一个统一的网格从x, y使用(x, y) = meshgrid (x, y),但我怎么生成数组从U ?我想使用插值?

我的代码

                               模型= mphload (“缸flow.mph”);
                              
                               速度= mpheval(模型,“spf.U”);
                              
                               密度= mpheval(模型,“spf.rho”);
                              
                               x = Velocity.p(1,1:结束);
                              
                               y = Velocity.p(2, 1:结束);
                              
                               xindex = Velocity.t(1,1:结束);
                              
                               yindex = Velocity.t(2, 1:结束);
                              
                               zindex = Velocity.t(3、1:结束);
                              
                               meshindex = Velocity.ve(1:结束,1);
                              
                               U = Velocity.d1;
                              
                               ρ= Density.d1;
                              
                               (X, Y) = meshgrid (X, Y);
                              
                               图(1)%网显示,效果很好
                              
                               mphmesh(模型);持有在;
                              
                               图(2)%散点图,效果很好
                              
                               散射(x, y, U,“填充”);
                              
                               图(3)% 3 d trimesh图——没问题
                              
                               三=德劳内(x, y);
                              
                               trisurf(三,x, y, U);colorbar;
                              
                               图(4)% contourf图不工作
                              
                               contourf (X, Y, U);

约翰D 'Errico 2020年7月1日

@员工切赫Michalek:

Irecall,有两个工具一个轮廓的分散数据文件交换。也许两者都是calledd tricontour,所以要么应该工作。

登录置评。

在回答这个问题。

答案(2)

阿施施Gudla 2015年5月27日

1
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/218758-2-d-contour-plot-for-non-uniformly-spaced-data-points answer_180713

如果我理解正确的话,你正试图画出原始数据,而不是创建一个正则网格使用“meshgrid”。

您可以使用“等值”或“contourf”功能的解释文档:

//www.tianjin-qmedu.com/help/matlab/ref/contour.html

还“网格”功能可以间隔不规则向量输入。看到这了一个例子。

2的评论
显示1年长的评论藏1年长的评论

一个 2015年5月28日

轮廓只能当“u”(z)是一个矩阵x和y的点。

“如果X和Y是向量,那么长度(X)必须等于大小(Z, 2)和长度(Y)必须等于大小(Z, 1)。向量必须严格严格增加或减少,不能包含任何重复的价值观”。

我的数据是3 向量成千上万的点,但可能有更多比y = x, y = 0 1因为ComSol不是均匀的网格。所以我不能简单地重塑你成一个矩阵不丢失数据。这是我的问题。

登录置评。

嘉里蒂迈克 2015年5月27日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/218758-2-d-contour-plot-for-non-uniformly-spaced-data-points answer_180723

轮廓的命令可以处理曲线网格,但不是多分辨率网格。所以我可以这样做:

                             (θ,rad) = meshgrid (linspace(0,π,49),linspace (1、3、49));
                            
                             (x, y) = pol2cart(θ,rad);
                            
                             z =山峰;
                            
                             contourf (x, y, z)

得到一个轮廓的2 d网格包装在半环,但X, Y, Z的输入都是二维数组。

一种选择是将数据转换成一个三角形网格,然后使用达伦Engwirda tricontour文件交换的功能。

另一个是使用一些其他的轮廓。例如,如果您的网格足够好,你可以使用分散在每个位置画彩色圆圈:

散射(x, y, 120 c,“填充”)

会给你任何插值,就原始数据。

5个评论
显示4年长的评论隐藏4年长的评论

嘉里蒂迈克 2015年5月28日

看着你的照片,看起来有一些机会,你向量实际上是所谓的样本直线网格。如果是这样的话,那么以下的方法可能会奏效。

让我们通过一个例子。我没有你的数据,所以我把我自己的。

                                  % %直线网格上创建一些样例数据
                                 
                                  (x, y) = meshgrid (6 * sin (linspace(0,π,32)3);
                                 
                                  z =山峰(x, y);
                                 
                                  % %转化为向量
                                 
                                  x = x (:);
                                 
                                  y = y (:);
                                 
                                  z z = (,);

现在让我们看看我们是否能重建直线网格。我们会发现独特的X和Y值,然后看看有一个值为每个组合。

                                  gx =独特的(x);
                                 
                                  gy =独特(y);
                                 
                                  广州= 0(长度(gx), (gy));
                                 
                                  数量= 0(长度(gx), (gy));%错误检查
                                 
                                  为i = 1:长度(x)
                                 
                                  c =找到(gx = = x (i));
                                 
                                  r =找到(gy = = y (i));
                                 
                                  广州(r、c) = z(我);
                                 
                                  计数(r、c) =计数(r、c) + 1;
                                 
                                  结束
                                 
                                  disp ([“数=”num2str (min(计数(:))),…
                                 
                                  ”、“num2str (max(计数(:)))));

你可以看到它不太适合我,因为有些数不等于1。发生了什么事,有X和Y值非常接近对方。

但让我们来看看不管怎样,因为结果可能不够好:

contourf (gx gy,广州)

如果不工作,那么我认为你坚持我提到的其他方法之一。

                                  % %散射方法
                                 
                                  散射(120 x, y,, z,“填充”)

                                  % % Trimesh方法
                                 
                                  三=德劳内(x, y);
                                 
                                  trisurf(三,x, y, z)