Question

我正在尝试使用pm3d为绘图表面着色，并希望使用来自另一列而不是z轴的值进行着色。

输入文件（test.file，制表符分隔）为：

atom_num    residue_name    X   Y   Z
288 1   45.3    36.6    79.3
301 1   38.9    197.4   72.5
314 1   118.2   53.8    76.5
327 1   58.2    139.1   78.5
353 1   1.9 14.4    71.9
366 1   156.9   180.0   72.1
379 1   183.2   5.4 69.5
392 1   71.7    155.4   75.8
457 1   83.4    11.8    74.8
613 1   97.1    180.7   77.5
626 1   145.2   160.3   71.7
678 2   73.1    76.3    81.0
704 3   30.3    46.5    79.3
717 2   216.0   130.7   85.5
743 2   55.0    137.2   74.4
756 2   23.4    67.3    78.3
769 2   46.9    156.1   77.3
821 2   145.4   143.9   80.7
990 2   7.8 119.3   79.8
1016    3   44.3    67.3    76.7
1042    3   12.8    44.4    74.3
1055    3   149.1   79.9    78.2
1068    3   100.8   35.8    76.1
1081    3   57.6    196.8   76.8
1094    3   214.7   122.8   79.5
1107    3   82.0    190.0   74.4
1120    3   150.9   39.4    71.3
1133    3   50.4    143.7   75.3
1146    1   42.9    104.7   74.3
1159    1   139.0   48.8    73.4
1172    1   66.8    165.3   71.5
1198    1   190.7   150.1   84.2
1211    1   92.1    5.1 75.8
1224    1   211.8   177.7   74.1
1237    1   131.6   0.2 73.6
1250    2   103.8   104.2   76.6
1276    2   132.4   5.0 70.0
1289    2   94.4    9.4 73.0
1302    2   72.6    33.7    74.3
1315    2   14.4    162.6   74.7
1406    2   171.4   143.6   86.1
1419    2   209.5   52.9    77.4
1445    2   11.6    14.7    72.3
1458    1   115.5   165.0   73.0
1549    1   147.1   45.5    76.1
1575    1   115.8   36.6    74.5
1588    1   35.8    37.3    76.2
1601    1   65.4    28.2    76.9
1614    1   13.4    199.9   76.5

我使用的命令是：

set dgrid3d 30,30
set hidden3d
set palette rgbformulae 33,13,10
splot "test.file" u 3:4:5 w pm3d

图像显示如下：

默认情况下，绘图基于Z轴值着色（第5列）。我使用剩余名称（第2列）的值对图进行着色，其范围为1-3。是否可以选择定义哪个颜色用于着色？理想情况下，我希望有相同的情节，但根据第2列进行着色，这样我就可以看到哪种“残留类型”在哪个轮廓中。

任何帮助都会非常有帮助。

Answer 1

由于您的residue是一个整数，因此不清楚是否要将其插入到网格中。

但是，如果这是您想要的，您可以使用Plotting 3D surface from scatter points and a png on the same 3D graph中的解决方案，但在编写表格时不要使用with pm3d。这是一个使用快速且有点脏的unix技巧来合并表的解决方案：

set terminal push #Save current terminal settings
set terminal unknown #dummy terminal
set table "surface.dat"
set dgrid3d
splot 'test.dat' using 3:4:5
set table "residue.dat"
splot 'test.dat' using 3:4:2
unset dgrid3d
unset table
set term pop #reset current terminal settings
!paste surface.dat residue.dat > test_grid.dat
splot "test_grid.dat" u 1:2:3:7 w pm3d

使用列值为pm3d表面着色

1 个答案: