如何获得平滑的边缘contourf情节规模对数散点图吗?

19日视图(30天)

显示旧的评论

欧文Brandriss 2023年5月11日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1961629-how-do-i-get-smooth-edges-for-a-contourf-plot-on-a-log-log-scale-scatterplot

评论道: 马修•诺伊 2023年6月12日在七51

答:接受马修•诺伊

你好,

我有一些麻烦与使用contourf函数对数尺度阴谋。我9点在二维彩色第三变量的散点图。我使用的代码绘制的数据以及情节包括如下。

一个=重塑(mtot_1 1 []);%转换矩阵行向量

b =重塑(MFR_1 1 []); %转换矩阵行向量

c =重塑(SN_maxes_1 1 []); %转换矩阵行向量

图(4)

clf

持有在

散射(b, a, [], c, “填充” )

集(gca), “xscale” , “日志” )

集(gca), “yscale” , “日志” )

colorbar

包含( “生产商” )

ylabel ( 总质量流量的 )

正如你所看到的,在对数尺度,datapoits形成一种“扭曲四边形”形状的边缘看起来“直”当绘制对数。我想创建一个等高线图从这些9分,但当我做的,它看起来像下面的情节因为生成等高线之间的线性插值方法创建直线上的点正常线性轴,然后看扭曲或弯曲时绘制对数刻度。我也包括我使用的代码。

图(5)

clf

持有在

contourf (MFR_1 mtot_1 SN_maxes_1, 100年, “线型” , “没有” )

散射(b, a, [], c, “填充” )

集(gca), “xscale” , “日志” )

集(gca), “yscale” , “日志” )

d = colorbar;

d.Label。字符串= “漩涡”。 ;

包含( “生产商” )

ylabel ( 总质量流量的 )

我想让等高线图的“直”边缘外点之间绘制对数刻度,这样等高线图本质上表现为一个四边形与直边的重对数坐标图而不是奇怪的弯曲的形状在上面的等高线图。有人能给我一些建议,如何实现这一目标?提前谢谢!

9日评论
显示8年长的评论隐藏8年长的评论

欧文Brandriss 2023年5月12日

太好了,谢谢!这是我使用的数据图:

首先是轴数据(它是一个矩阵):

MFR_1 = [0.93016, 0.13933, 0.04154;4.75072、0.96454、0.27638;16.1767,3.35929,1.03684)

然后y轴数据(矩阵):

mtot_1 = [0.00087393, 0.001293, 0.00161739;0.00146412、0.00182395、0.00211802;0.00195069,0.00228598,0.002528465)

然后“z”的数据(如果你会称呼它)。这是确定点的颜色。这也是一个矩阵:

SN_maxes_1 = [1.678801, 1.627564, 1.521288;1.535838、1.848008、1.7666569;1.419559,1.818278,1.963394)

我目前使用contourf函数绘制等高线与这行代码(你可以看到它在原文中的代码):

contourf (MFR_1 mtot_1 SN_maxes_1, 100年,“线型”,“没有一个”)

情节只是我使用的点(没有轮廓)散射函数。这个函数需要的输入是一个向量,所以我使用以下代码行改变矩阵向量和情节:

一个=重塑(MFR_1 1 []);%转换矩阵行向量

b =重塑(mtot_1 1 []);%转换矩阵行向量

c =重塑(SN_maxes_1 1 []); %转换矩阵行向量

散射(a, b, [], c,“填满”)

设置轴规模对数,我使用这些代码行数:

集(gca), “xscale” , “日志” )

集(gca), “yscale” , “日志” )

这应该只是覆盖了一切。让我知道如果有任何其他信息我可以提供帮助。提前再次感谢你的帮助!

马修•诺伊 2023年5月17日

像往常一样,这是我的荣幸!

登录置评。

在回答这个问题。

接受的答案

马修•诺伊 2023年5月15日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1961629-how-do-i-get-smooth-edges-for-a-contourf-plot-on-a-log-log-scale-scatterplot answer_1236894

移动:马修•诺伊 2023年5月16日

你好

我可以建议,虽然我认为这需要一些改进

MFR_1 = [0.93016, 0.13933, 0.04154;4.75072、0.96454、0.27638;16.1767,3.35929,1.03684);

%的轴数据(矩阵):

mtot_1 = [0.00087393, 0.001293, 0.00161739;0.00146412、0.00182395、0.00211802;0.00195069,0.00228598,0.002528465);

%的“z”数据(如果你会称呼它)。这是确定点的颜色。这也是一个矩阵:

SN_maxes_1 = [1.678801, 1.627564, 1.521288;1.535838、1.848008、1.7666569;1.419559,1.818278,1.963394);

一个=重塑(mtot_1 1 []);%转换矩阵行向量

b =重塑(MFR_1 1 []);%转换矩阵行向量

c =重塑(SN_maxes_1 1 []);%转换矩阵行向量

图(4)

clf

持有在

散射(b, a, [], c,“填充”)

%用常数日志创建一个网格间距,并找到里面的点

%一个多边形(凸包)

提单= log10 (b);

艾尔= log10 (');

cl = log10 (c ');

k =边界(提单,艾尔,1);%定义外部船体

十五=提单(k);

青年志愿= al (k);

% (10。^十五10。^青年志愿,- r)

x = linspace (min(十五),max(十五),200);

y = linspace (min(青年志愿),max(青年志愿),200);

(X, Y) = meshgrid (X, Y);

x = x (:);

y = y (:);

在= inpolygon (x, y,十五,青年志愿);

鑫= 10 ^ (X ());

阴= 10。^ (Y ());

vq = griddata(提单,艾尔c log10(新),log10(阴));

散射(鑫,阴,[],矢量量化,“填充”)

%的阴谋(鑫,阴,“r”)

集(gca),“xscale”,“日志”)

集(gca),“yscale”,“日志”)

colorbar

包含(“生产商”)

ylabel (总质量流量的)

持有从

4评论
显示3年长的评论隐藏3年长的评论

欧文Brandriss 2023年6月8日18:54

嗨 @Mathieu诺 :

我有一个后续问题上面的情节。我试图让相同类型的情节,但是具有不同的值在轴上。的值轴和colorbar是相同的。我用下面的代码:

                                   图(4)
                                  
                                   clf
                                  
                                   持有在
                                  
                                   一个=重塑(mtot_1 1 []);%转换矩阵行向量
                                  
                                   a2 =重塑(MF_tot_1 1 []);%转换矩阵行向量
                                  
                                   b =重塑(MFR_1 1 []);%转换矩阵行向量
                                  
                                   c =重塑(SN_maxes_1 1 []);%转换矩阵行向量
                                  
                                   %这情节是一样的我们之前讨论的
                                  
                                   次要情节(1、2、1)
                                  
                                   %轮廓策划- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                                  
                                   提单= log10 (b);
                                  
                                   艾尔= log10 (');
                                  
                                   cl = log10 (c ');
                                  
                                   k =边界(提单,艾尔,1);%定义外部船体
                                  
                                   十五=提单(k);
                                  
                                   青年志愿= al (k);
                                  
                                   % (10。^十五10。^青年志愿,- r)
                                  
                                   x = linspace (min(十五),max(十五),1000);
                                  
                                   y = linspace (min(青年志愿),max(青年志愿),1000);
                                  
                                   (X, Y) = meshgrid (X, Y);
                                  
                                   x = x (:);
                                  
                                   y = y (:);
                                  
                                   在= inpolygon (x, y,十五,青年志愿);
                                  
                                   鑫= 10 ^ (X ());
                                  
                                   阴= 10。^ (Y ());
                                  
                                   vq = griddata(提单,艾尔c log10(新),log10(阴));
                                  
                                   散射(鑫,阴,[],矢量量化,“填充”)
                                  
                                   % - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                                  
                                   集(gca),“xscale”,“日志”)
                                  
                                   集(gca),“yscale”,“日志”)
                                  
                                   colorbar
                                  
                                   d = colorbar;
                                  
                                   d.Label。字符串=“漩涡”。;
                                  
                                   包含(“生产商”)
                                  
                                   ylabel (总质量流量的)
                                  
                                   盒子在
                                  
                                   网格在
                                  
                                   %这是新的情节(有问题)
                                  
                                   次要情节(1、2、2)
                                  
                                   %轮廓策划- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                                  
                                   提单= log10 (b);
                                  
                                   艾尔= log10 (a2);%使用a2代替轴值
                                  
                                   cl = log10 (c ');
                                  
                                   k =边界(提单,艾尔,1);%定义外部船体
                                  
                                   十五=提单(k);
                                  
                                   青年志愿= al (k);
                                  
                                   % (10。^十五10。^青年志愿,- r)
                                  
                                   x = linspace (min(十五),max(十五),1000);
                                  
                                   y = linspace (min(青年志愿),max(青年志愿),1000);
                                  
                                   (X, Y) = meshgrid (X, Y);
                                  
                                   x = x (:);
                                  
                                   y = y (:);
                                  
                                   在= inpolygon (x, y,十五,青年志愿);
                                  
                                   鑫= 10 ^ (X ());
                                  
                                   阴= 10。^ (Y ());
                                  
                                   vq = griddata(提单,艾尔c log10(新),log10(阴));
                                  
                                   散射(鑫,阴,[],矢量量化,“填充”)
                                  
                                   % - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
                                  
                                   集(gca),“xscale”,“日志”)
                                  
                                   集(gca),“yscale”,“日志”)
                                  
                                   colorbar
                                  
                                   d = colorbar;
                                  
                                   d.Label。字符串=“漩涡”。;
                                  
                                   包含(“生产商”)
                                  
                                   ylabel (总动量通量的)
                                  
                                   盒子在
                                  
                                   网格在
                                  
                                   集(gcf,“颜色”,' w ');
                                  
                                   sgtitle (轮廓的漩涡。”)

这段代码将生成以下情节:

我的问题是这样的:我怎么能让它右边的情节有一个完整的轮廓(这样就不会失踪triagular块的轮廓在左下角)。左边的轮廓是我们前面所讨论的,是完美的。我只是想让右边的匹配方面,拥有一个完整的轮廓。

就像我说的,数据都是一样的,除了轴在第二个情节。的值MF_tot_1(用于轴值在第二个图)如下:

提前感谢你的帮助。

马修•诺伊 2023年6月9日11

你好

您只需要调整阈值边界得到正确的结果

所以对于第二个情节,必须降低阈值从1到0.5

这条线:

k =边界(提单,艾尔,1);%定义外部船体

现在是:

k =边界(提单,艾尔,0.5);%定义外部船体

完整的代码:

MFR_1 = [0.93016, 0.13933, 0.04154;4.75072、0.96454、0.27638;16.1767,3.35929,1.03684);

%的轴数据(矩阵):

mtot_1 = [0.00087393, 0.001293, 0.00161739;0.00146412、0.00182395、0.00211802;0.00195069,0.00228598,0.002528465);

%的“z”数据(如果你会称呼它)。这是确定点的颜色。这也是一个矩阵:

SN_maxes_1 = [1.678801, 1.627564, 1.521288;1.535838、1.848008、1.7666569;1.419559,1.818278,1.963394);

MF_tot_1 = [0.001616265021936 0.005243033092750 0.015546133692972;

0.005271449006048 0.007910496071212 0.017789070120675;

0.014974604975415 0.016927091212513 0.025444016339618);

一个=重塑(mtot_1 1 []);%转换矩阵行向量

a2 =重塑(MF_tot_1 1 []);%转换矩阵行向量

b =重塑(MFR_1 1 []);%转换矩阵行向量

c =重塑(SN_maxes_1 1 []);%转换矩阵行向量

%这情节是一样的我们之前讨论的

次要情节(1、2、1)

%轮廓策划- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

提单= log10 (b);

艾尔= log10 (');

cl = log10 (c ');

k =边界(提单,艾尔,1);%定义外部船体

十五=提单(k);

青年志愿= al (k);

x = linspace (min(十五),max(十五),1000);

y = linspace (min(青年志愿),max(青年志愿),1000);

(X, Y) = meshgrid (X, Y);

x = x (:);

y = y (:);

在= inpolygon (x, y,十五,青年志愿);

鑫= 10 ^ (X ());

阴= 10。^ (Y ());

vq = griddata(提单,艾尔c log10(新),log10(阴));

散射(鑫,阴,[],矢量量化,“填充”)

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

集(gca),“xscale”,“日志”)

集(gca),“yscale”,“日志”)

colorbar

d = colorbar;

d.Label。字符串=“漩涡”。;

包含(“生产商”)

ylabel (总质量流量的)

盒子在

网格在

%这是新的情节(有问题)

次要情节(1、2、2)

%轮廓策划- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

提单= log10 (b);

艾尔= log10 (a2);%使用a2代替轴值

cl = log10 (c ');

k =边界(提单,艾尔,0.5);%定义外部船体

十五=提单(k);

青年志愿= al (k);

x = linspace (min(十五),max(十五),1000);

y = linspace (min(青年志愿),max(青年志愿),1000);

(X, Y) = meshgrid (X, Y);

x = x (:);

y = y (:);

在= inpolygon (x, y,十五,青年志愿);

鑫= 10 ^ (X ());

阴= 10。^ (Y ());

vq = griddata(提单,艾尔c log10(新),log10(阴));

散射(鑫,阴,[],矢量量化,“填充”)

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

集(gca),“xscale”,“日志”)

集(gca),“yscale”,“日志”)

colorbar

d = colorbar;

d.Label。字符串=“漩涡”。;

包含(“生产商”)

ylabel (总动量通量的)

盒子在

网格在

集(gcf,“颜色”,' w ');

sgtitle (轮廓的漩涡。”)

马修•诺伊 2023年6月12日在七51

我的荣幸!

祝你有美好的一天

登录置评。

答案(1)

Sulaymon Eshkabilov 2023年5月12日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1961629-how-do-i-get-smooth-edges-for-a-contourf-plot-on-a-log-log-scale-scatterplot answer_1234009

这里是如何得到的直线边缘:

MFR_1 = [0.93016, 0.13933, 0.04154;4.75072、0.96454、0.27638;16.1767,3.35929,1.03684);

%的轴数据(矩阵):

mtot_1 = [0.00087393, 0.001293, 0.00161739;0.00146412、0.00182395、0.00211802;0.00195069,0.00228598,0.002528465);

%的“z”数据(如果你会称呼它)。这是确定点的颜色。这也是一个矩阵:

SN_maxes_1 = [1.678801, 1.627564, 1.521288;1.535838、1.848008、1.7666569;1.419559,1.818278,1.963394);

一个=重塑(mtot_1 1 []);%转换矩阵行向量

b =重塑(MFR_1 1 []);%转换矩阵行向量

c =重塑(SN_maxes_1 1 []);%转换矩阵行向量

图(4)

clf

持有在

散射(b, a, [], c,“填充”)

集(gca),“xscale”,“日志”)

集(gca),“yscale”,“日志”)

colorbar

包含(“生产商”)

ylabel (总质量流量的)

图(5)

clf

持有在

contourf(日志(MFR_1),日志(mtot_1),日志(SN_maxes_1), 100年,“线型”,“没有”)

散射(日志(b),日志(a),[],日志(c),“填充”)

集(gca),“xscale”,“线性”)

集(gca),“yscale”,“线性”)

d = colorbar;

d.Label。字符串=“漩涡”。;

包含(“生产商”)

ylabel (总质量流量的)

1评论
显示没有隐藏没有

欧文Brandriss 2023年5月15日

编辑:欧文Brandriss 2023年5月15日

@Sulaymon Eshkabilov 这是非常接近我想要的,但是我需要轴上的值匹配的数据。情节上的轴值显示的日志值据点,但我需要原始值。换句话说,我需要点一样在我提供的等高线图,只是直边缘轮廓。轴也应该有一个日志

登录置评。

在回答这个问题。

类别

MATLAB 图形二维和三维图表面,卷和多边形表面和网格图

找到更多的在表面和网格图在帮助中心和文件交换

s manbetx 845

MATLAB

释放

R2022a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

如何获得平滑的边缘contourf情节规模对数散点图吗?

9日评论
显示8年长的评论隐藏8年长的评论

接受的答案

4评论
显示3年长的评论隐藏3年长的评论

答案(1)

1评论
显示没有隐藏没有

另请参阅

类别

标签

s manbetx 845

释放

社区寻宝

如何获得平滑的边缘contourf情节规模对数散点图吗?

9日评论 显示8年长的评论隐藏8年长的评论

接受的答案

4评论 显示3年长的评论隐藏3年长的评论

答案(1)

1评论 显示没有隐藏没有

另请参阅

类别

标签

s manbetx 845

释放

社区寻宝

9日评论
显示8年长的评论隐藏8年长的评论

4评论
显示3年长的评论隐藏3年长的评论

1评论
显示没有隐藏没有