使用surfer创建一个圆形mask

Surfer是Python的一个图像处理库，主要用于地理空间数据的可视化。如果你想在Surfer中创建一个圆形mask，这通常是在地图上应用某个范围限制时的做法。以下是一个简单的步骤：

1. 首先，你需要安装`matplotlib`和`cartopy`库，这两个库是Surfer的依赖，用于地图绘制。你可以通过命令行运行`pip install matplotlib cartopy`来安装。

2. 导入所需的库和模块：

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from surfer import GeoAxesSubplot

3. 定义圆的中心经纬度和半径：

center_lon = -90.0  # 圆心经度
center_lat = 40.0   # 圆心纬度
radius_km = 100      # 半径，单位公里

4. 创建GeoAxes并转换坐标系统以便计算圆：

fig = plt.figure(figsize=(10, 10))
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())

5. 计算圆的边界点，并创建mask：

# 将半径从km转换到地图投影坐标
radius_map = radius_km * 1e3  # 因为Surfer默认是米

# 创建圆的边界点列表
theta = np.linspace(0, 2*np.pi, 180)
lon, lat = center_lon + radius_map * np.cos(theta), center_lat + radius_map * np.sin(theta)

# 创建圆形mask数组
mask = ~np.isnan(lon) & ~np.isnan(lat)

6. 绘制地图和mask：

ax.coastlines(resolution='10m')
ax.set_global()
ax.scatter(lon[mask], lat[mask], transform=ccrs.Geodetic(), color='red', alpha=0.5)  # 标记圆内点
ax.imshow(mask.reshape(ax.projection.grid_shape), extent=ax.get_extent(), cmap='binary', interpolation='nearest', zorder=10)
plt.show()