城市表层土壤重金属污染分析污染元素热力图代码

以下是一个示例代码，用于制作城市表层土壤重金属污染元素热力图：

import pandas as pd
import geopandas as gpd
import matplotlib.pyplot as plt

# 读取重金属含量数据，数据格式示例：采样点位经度、纬度、重金属含量
data = pd.read_csv('soil_data.csv')

# 创建地理空间数据框
gdf = gpd.GeoDataFrame(data, geometry=gpd.points_from_xy(data['经度'], data['纬度']))

# 选择需要制作热力图的重金属元素列
element = '铅'

# 绘制热力图
fig, ax = plt.subplots(figsize=(10, 8))
gdf.plot(column=element, cmap='Reds', ax=ax, legend=True)

# 设置图像标题和坐标轴标签
ax.set_title('城市表层土壤重金属污染热力图')
ax.set_xlabel('经度')
ax.set_ylabel('纬度')

# 显示图像
plt.show()

在上述代码中，需要将重金属含量数据保存为名为 'soil_data.csv' 的 CSV 文件，并确保数据中包含经度和纬度信息。可以根据实际情况修改数据文件的路径和列名。

代码中使用了 pandas 库和 geopandas 库来处理和可视化地理空间数据，同时使用了 matplotlib 库来绘制热力图。通过选择需要制作热力图的重金属元素列，并设置相应的颜色映射（这里使用了 'Reds'），可以生成相应元素的热力图。

请注意，上述代码仅为示例，实际应用时需要根据自己的数据和需求进行适当的修改和调整。