### 使用 Python 创建全球地震散点图 为了创建全球地震数据的散点图，可以利用 `matplotlib` 和 `geopandas` 库来处理地理空间数据并绘制地图。下面是一个简单的例子说明如何实现这一目标。 #### 导入必要的库 首先安装所需的包可以通过 pip 安装命令完成： ```bash pip install pandas geopandas matplotlib requests ``` 接着，在脚本中导入这些库： ```python import pandas as pd import geopandas as gpd from shapely.geometry import Point import matplotlib.pyplot as plt import requests ``` #### 获取地震数据 可以从美国地质调查局 (USGS) 的 API 下载最新的地震活动记录作为 JSON 文件格式的数据源[^1]。 ```python url = 'https://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/all_week.geojson' response = requests.get(url) data = response.json() features = data['features'] df_earthquakes = pd.DataFrame(features) # 提取经纬度和其他有用的信息到 DataFrame 中 geometry = [Point(xy) for xy in zip(df_earthquakes['geometry'].apply(lambda x: x['coordinates'][0]), df_earthquakes['geometry'].apply(lambda x: x['coordinates'][1]))] crs = {'init': 'epsg:4326'} gdf_earthquakes = gpd.GeoDataFrame(df_earthquakes, crs=crs, geometry=geometry) ``` #### 绘制世界地图和地震位置 加载自然地球的世界边界文件，并将其与地震地点一起绘制成图形。 ```python world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres')) fig, ax = plt.subplots(figsize=(15, 8)) world.plot(ax=ax, color='lightgray') # 添加地震事件的位置标记 magnitude = gdf_earthquakes.properties.apply(lambda x: x['mag']) gdf_earthquakes.plot(ax=ax, markersize=magnitude*2, alpha=0.7, edgecolor='k', marker="o", label="Earthquakes") plt.title("Global Earthquake Locations Over the Past Week") plt.legend() plt.show() ``` 通过上述代码片段，能够成功地获取最近一周内发生的地震信息，并在全球地图上以不同大小的圆圈表示各个地震的发生地点及其相对强度[^1]。