Python球体磁异常正演

球体磁异常正演是一种模拟方法，用于计算球体磁源对地表磁场的影响。在Python中，可以利用物理学公式和数值计算方法来实现球体磁异常正演。

以下是一个简单的示例代码，展示了如何计算球体磁异常正演：

import numpy as np

def magnetic_anomaly(x, y, z, x0, y0, z0, radius, magnetization):
    """
    计算球体磁异常正演
    :param x: 采样点x坐标
    :param y: 采样点y坐标
    :param z: 采样点z坐标
    :param x0: 磁源x坐标
    :param y0: 磁源y坐标
    :param z0: 磁源z坐标
    :param radius: 球体半径
    :param magnetization: 球体磁化强度
    :return: 磁异常值
    """
    r = np.sqrt((x - x0) ** 2 + (y - y0) ** 2 + (z - z0) ** 2)
    B = (magnetization * (radius ** 3)) / (4 * np.pi * r ** 3)
    return B

# 示例参数
x_sample = np.linspace(-100, 100, 100)
y_sample = np.linspace(-100, 100, 100)
z_sample = np.zeros_like(x_sample)
x0 = 0
y0 = 0
z0 = 100
radius = 50
magnetization = 1

# 计算磁异常值
magnetic_anomaly_values = magnetic_anomaly(x_sample, y_sample, z_sample, x0, y0, z0, radius, magnetization)

在上述示例代码中，我们定义了一个 `magnetic_anomaly` 函数来计算球体磁异常值。函数根据传入的采样点坐标、磁源坐标、球体半径和磁化强度，使用物理学公式计算球体磁异常值。