怎么样在电脑上画三维图(Nature论文中的三维图原来是这么画的)

高水平科研结果必须要用高水平数据可视化方式展现给编辑、审稿人和读者！利用GMT可以容易的实现三维数据可视化，这里将我多年摸索得来的独家秘方全部奉上！用GMT绘制三维图，学这一个例子就够了。提供全套解决方案，包括修复三维绘图bug之后的独家GMT版本，还有教程中所需要的所有的运行环境和程序，全部打包为Docker镜像。再也不用担心有些命令无法在你电脑上运行了。Nature子刊作者手把手教学，附带三小时4K教学视频和全部代码，将三维绘图和一些常用gmt命令和数据处理一网打尽！如果想看全部的视频和绘图代码，请点击文末的【了解更多】按钮。

幻灯片预览

三维数据可视化教程幻灯片预览

教程内容

效果展示

先上一个效果展示图，符合你的预期再往下看

图1NatureComms.(10.1038/s41467-020-15062-w)论文中的图1用GMT绘制的效果

绘图过程简介

GMT绘图复杂组合图，建议采用先拆分后组合的策略！

Figure1(a,b)

这是两个二维图，完全没有难度，直接上图。

Figure1(a)

Figure1(b)

Figure1(c)

海底地形三维立体图+三维散点，这个必须使用我修复bug之后的gmt才能画出这种效果，所需要的所有程序都在docker镜像里面。只需简单的三个命令就可以使用，完全无压力！

Figure1(d)

海底地形三维立体图+三维散点+切片数据+切片投点和曲线，这个必须使用我修复bug之后的gmt才能画出这种效果，所需要的所有程序都在docker镜像里面。只需简单的三个命令就可以使用，完全无压力！

最终组合

图1的最主要的四个子图，每一个都能够单独绘制，然后将其进行组合布局。因为gmt不是一个所见即所得的绘图软件，所以为了更方便的设置每个子图的位置，可以通过一个A4纸大小的网格作为辅助底图进行布局调整。

三小时高清精讲视频和全部绘图代码

如果想看全部的视频和绘图代码，请点击文末的【了解更多】按钮。

图1的不完整绘图代码示例

总共两百行代码就可搞定上面的三维+二维组合图

functionpreset(){path_data=../../data/rawbathy=${path_data}/bathy/bathy_Longqi_cart.ncbasalt=${path_data}/sites/site_basalt_cart_zhao_bathy.txt#数据范围zmin=-6zmax=-1.5zrange=$zmin/$zmaxxmin=-18.35xmax=33.35ymin=-29.9826ymax=11.653range=$xmin/$xmax/$ymin/$ymaxrange_xyz=${range}/$zrangexmin_eq=0xmax_eq=20ymin_eq=-10range_xyz_eq=$xmin_eq/$xmax_eq/$ymin_eq/$ymax_eq/$zmin_eq/$zmax_eq#视角angle_h=-40angle_v=45angle_view=${angle_h}/$angle_v#图片大小width_x_fig=8}gmtbeginFigure1pdf#(a)move_x=0move_y=12#gmtbasemap-JM$figwidth_inset-R$range_inset-Bwsen--MAP_FRAME_TYPE=plain#gmtbasemapgmtcoast-Ggray-JL$lon0_inst/$lat0_inst/$lat_min_inset/$lgmtgrdcontour$bathy_AUV-C0.02-A+f3p,Helvetica+ggreen-W0.2p-GD500e-p#采样点symbol_sample=cgmtpsxy$path_data/sites/M-active.txt-S$symbol_sample$si#legend#(d)move_x=-1move_y=-4gmtbasemap-JX$width_x_fig_eq/$width_y_fig_eq-JZ$width_zgmtgrdview$bathy-JZ-p-C$bathy_cpt-N${zmin_eq}+g$color_profile-Qi-I0+a45-Wf0.5p,blackawk'{print$1/1000.0,$2/1000.0}'$earthquake|gmtgrdtrack-G$bathy|gmtpsxyz-JZ-p-Sc$size_eq-Ggreen-Wblack#东西方向gmtbasemap-JX$width_x_fig_eq/$width_z_fig_eq-JZ$width_y_#boundarybetweentitleandlegend#symbolspacelegendtypelength(-)thickness,colorspaceLegendName#N3H10pHelveticaLegendoffigure(c)N2G${gap_rows}cS${left_symbol}ca0.3c$color_active0.2p,$color_active${symbol_text}cActivehydrothermalfieldS1cc${size_symbol}c$color_basalt0gmtendrmtmp*gmt.historybox.eq