两次插值构建三维速度模型
2023-03-04
在三维地震成像中通常需要构建三维初始模型，这些模型往往来自已有的公共速度模型或其他前人研究，尽管提供的模型多为规则网格，由于我们所需网格大小是自定义的，所以需要对其进行网格化。而网格化所用的插值算法多为线性插值，它存在一个弊端：如果我们的网格范围比原先的大，将会出现NaN结果，无法用于后续计算。
491 words
|
2 minutes
将Fortran中的数组与矩阵保存为numpy中的.npz文件
2022-02-01
我们之前讨论过如何用python在三维矩阵中切片，在Python中保存三维矩阵和它的三个坐标轴的方法有用numpy.savez将数组和矩阵用不同的键值对保存至一个.npz文件。这里将讨论如何用Fortran实现这一过程。
398 words
|
2 minutes
用GMT绘制中国大陆及邻区地质年代图
2021-05-14
用GMT绘制绘制地质图一直是一个难题，此前Po-Chin Tseng的博客详细介绍了制作和绘制地质图的过程，并成功绘制了东南亚的地质图。本文将介绍用相似的方法绘制中国大陆及邻区的地质图。
630 words
|
3 minutes
用GMT、PyVista和Paraview对速度模型、地形起伏、Moho面起伏进行三维可视化
2021-03-06
我们通过计算获得了地壳三维速度模型、Moho面起伏，为例对比壳内低速区与Moho面起伏的关系，可以将这些结果用Paraview进行三维可视化。由于数据格式的不同，我们需要借助GMT和PyVista辅助绘制。
1347 words
|
7 minutes
用GMT绘制中国及周边构造分区与GPS速度场
2020-01-30
GMT对地图的处理十分优秀，而且只需要短短几行代码。本文介绍了如何用等距圆锥投影绘制不同范围大小的地形图，以及如何在地形图上绘制断层、文字和矢量等。
2087 words
|
10 minutes
判断地理坐标点是否在多边形内的Python实现
2019-10-08
在地学问题中经常会遇到判断点是否在多边形内的问题，例如成像问题中射线是否穿过某个网格即为此类问题。在Python中我们对处理这样的问题需要满足以下条件：
302 words
|
2 minutes