用python在全球地图上画岩石圈厚度热力图

要画岩石圈厚度热力图，可以使用Python的Basemap库和matplotlib库。以下是一个简单的例子：

首先，需要安装Basemap和matplotlib库，在命令行中输入以下命令：

pip install basemap
pip install matplotlib

然后，导入所需要的库：

import numpy as np
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt

接下来，需要准备数据。这里使用了一个全球岩石圈厚度的数据集，可以从这里下载：https://www.ngdc.noaa.gov/mgg/fliers/01mgg04.html。

数据集包含一个名为"agegrid.xy"的文本文件，其中包含了全球各个点的经纬度坐标和岩石圈厚度数据。因此，可以读取该文件并将数据分别存储到三个数组中：

# 读取数据
data = np.loadtxt("agegrid.xy")

# 分别存储经度、纬度和厚度数据
lon = data[:,0]
lat = data[:,1]
thickness = data[:,2]

接下来，需要使用Basemap库创建一个地图，并将数据绘制到地图上：

# 创建地图
m = Basemap(projection='robin',lon_0=0,resolution='c')
m.drawcoastlines(linewidth=0.5)
m.drawcountries(linewidth=0.5)

# 绘制热力图
x, y = m(lon, lat)
m.hexbin(x, y, C=thickness, gridsize=100, cmap=plt.cm.jet, edgecolors='none')
plt.colorbar()

# 显示地图
plt.show()

这里选择了Robinson投影，但可以根据需要选择其他投影。gridsize参数控制热力图的精度，cmap参数指定了颜色映射，edgecolors参数指定了边界线颜色。

运行代码，就可以得到一个全球岩石圈厚度热力图了。注意，绘制地图可能需要一些时间，取决于计算机性能和数据集大小。