手把手教你用Python绘制断层‘伍尔夫网’：从走向/倾角数据到极射赤平投影

手把手教你用Python绘制断层‘伍尔夫网’：从走向/倾角数据到极射赤平投影

地质构造分析中，断层的走向和倾角数据是理解地壳运动的关键参数。但将这些抽象的数字转化为直观的图形，一直是地质工作者和学生面临的挑战。极射赤平投影（伍尔夫网）作为经典的地质可视化工具，能将三维空间中的几何关系展现在二维平面上。本文将带你用Python的matplotlib和mplstereonet库，实现从原始数据到专业图形的全流程自动化处理。

1. 环境配置与数据准备

在开始之前，确保你的Python环境已安装以下库：

pip install matplotlib numpy mplstereonet

地质数据通常以走向（strike）和倾角（dip）的形式记录。走向指断层线与正北方向的夹角（0-360°），倾角则是断层面与水平面的夹角（0-90°）。我们用一个简单的CSV文件存储示例数据：

station,strike,dip
A,120,45
B,30,60
C,250,30

读取和处理这些数据的Python代码如下：

import pandas as pd

data = pd.read_csv('fault_data.csv')
strikes = data['strike'].values
dips = data['dip'].values

注意：地质学中走向的表示有不同惯例，本文采用"右旋走向"（右手定则），即从走向方向顺时针旋转90°得到倾向方向。

2. 创建基础伍尔夫网投影

mplstereonet库提供了极射赤平投影的核心功能。我们先创建一个基础投影网：

import matplotlib.pyplot as plt
import mplstere