### 使用Python按省份读取NetCDF（.nc）文件 为了实现按照特定省份来读取 NetCDF 文件中的数据，可以采用地理空间分析库 `geopandas` 和 `rioxarray` 来配合操作。具体来说，通过加载中国各省级行政区划的边界形状文件（Shapefile），并将其与 NetCDF 数据集相结合，从而筛选出指定省域内的气象或其他环境变量的数据。 #### 准备工作 安装必要的 Python 库： ```bash pip install netCDF4 geopandas rioxarray shapely fiona ``` #### 加载所需的模块和定义函数 ```python import geopandas as gpd from netCDF4 import Dataset import rioxarray as rxr import matplotlib.pyplot as plt ``` #### 获取中国的省级行政区域 Shapefile 并选择目标省份 假设已经有一个包含所有中国大陆省份边界的 shapefile 文件路径为 `china_provinces.shp`: ```python # 读入中国各省的矢量图层 provinces = gpd.read_file('path/to/china_provinces.shp') # 打印查看可用列名以便选取所需字段, 如 'NAME' 或者其他表示名称的属性 print(provinces.columns) # 基于某列值过滤得到单个省份的空间对象 selected_province = provinces[provinces['NAME'] == "广东省"] # 将 "广东省" 替换为目标省份的名字 ``` #### 处理 NetCDF 文件并与选定省份匹配 下面展示了一个简单的例子，说明如何打开一个 NetCDF 文件，并应用上述创建好的掩码来进行裁剪： ```python # 打开 .nc 文件作为 xarray.Dataset 对象 dataset = rxr.open_rasterio('example_data.nc', masked=True).squeeze() # 设置坐标系参数以确保两者一致 if not dataset.rio.crs: dataset = dataset.rio.write_crs("EPSG:4326") # EPSG代码取决于实际使用的投影系统 # 裁切至所选省份范围之内 clipped_dataset = dataset.rio.clip(selected_province.geometry.apply(mapping), selected_province.crs) # 可视化结果 fig, ax = plt.subplots(figsize=(8, 6)) clipped_dataset.plot(ax=ax) ax.set_title(f'Data clipped to {selected_province["NAME"].values[0]}') plt.show() ``` 以上过程实现了基于给定的 NetCDF 文件以及对应的省级行政区划信息对中国某一特定省份内数据的选择性提取[^1]。