Python 绘制地形与散点数据分布图
博客内容参考以下:
本文将介绍如何使用 Python 中的 Matplotlib、Cartopy、GDAL 库,读取 ETOPO1 海洋地形数据并绘制散点数据的分布图。该方法可以用于可视化地理数据,直观展示数据的空间分布特性。
功能概述
该代码实现以下功能:
- 读取地形数据(ETOPO1),提取目标区域水深信息。
- 生成随机散点数据,并叠加到地图上。
- 自定义地图样式,包括地形、散点、色标、边框和经纬度设置。
代码实现
1. 导入必要的库
import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from osgeo import gdal
确保已安装相关库:
pip install numpy matplotlib cartopy gdal
注意在Window中直接pip install gdal会失败,需要下载安装包再进行pip install,需要安装包的可私信小编。
MacOS和Linux则没有这个问题
2. 设置绘图范围与随机散点数据
设置经纬度范围
lonmin, lonmax = 105, 124
latmin, latmax = 12, 25
np.random.seed(42)
num_points = 100
longitudes = np.random.uniform(lonmin+5, lonmax-4, num_points)
latitudes = np.random.uniform(latmin+4, latmax-4, num_points)
angles = np.random.uniform(