Modern GMT Series：Slice in 3D View (三维切片图)

Fancy的版本见九天学者的个人博客，关注文集博士干点啥或者微信公众号九天学者及时获取连载更新。

相信理工科方面的科研人员经常会遇到做三维切片图的时候，以地球物理为例，比如重力、磁异常三维反演，地震速度剖面+地形图，热液循环三维数值模拟（CFD范畴），地球动力学三维模拟等。当然了做三维切片图也有商业软件(比如Tecplot 360)和开源软件（比如Paraview），前者缺点是收费收费收费，后者导出的三维场景图效果很差（不是矢量图，colorbar位置错乱等问题），其优点是做现场演示。还有我最喜欢的编程语言python也可以做这个工作，但是我觉得它对坐标轴的label处理的不是那么完美，但也是个不错的选择。我的原则就是用开源的程序或者软件做出最漂亮最简洁并且高质量(本人追求矢量的eps、pdf、ps格式)的学术论文图。经过我四年的不断折腾和筛选，发现到目前为止，GMT 6.0最新版（我是直接从github克隆下来的开发版）对三维图支持已经非常棒了。不仅可以做三维地形图（上面标注文字，数据点，线，legend等），还可以做切片的等值线图或者image。

一贯的风格：先来一段回顾和概述，然后提出问题，再解决问题，最后点评一下存在的问题或可选方案

提出问题

这种图在什么情况下使用？直接上图感受一下

Olive, et al., 2015, Science

Schmid F. and Schlindwein V., 2016, G3

Zhao, et al., 2013, G3

GMT如何实现三维切片图

gmt新版推出的三维成图主要通过-p和-JZ实现的。前者进行旋转，后者给定垂向轴相关信息。gmt官网上有个用户给出的例子，但是没有讲清楚这个切片图到底如何准确实现，但是还是有借鉴意义的。因为本人的博士课题有部分数据涉密，不方面放上来，就以上面提到的这个例子里面的函数造一个数据进行成图。

-p参数

这个参数是新版gmt里面一个很通用的参数，在basemap, grdimage, grdcontour, grdview等都有这个参数，下面以basemap为例介绍-p参数的用法。具体请参见官网英文解释

• 参数解释
  -p后面可以跟x,y,z（这三个字母表示你的此图的纵轴是哪个），然后再跟方位角(azimuth)/仰角(elevation)，或者直接跟方位角/仰角。方位角的范围是