以下是实现该功能的 Python 代码，使用 `geopandas` 库进行地理数据处理： ```python import os import geopandas as gpd # 定义基础路径 base_path = r'D:\A LINSHI\SHUJV\RENKOU\WGS - 84' output_base_path = r'D:\A LINSHI\SHUJV\RENKOU\CHULIHOU' boundary_path = r'D:\A LINSHI\SHUJV\XINGZHENGQU\HB BIANJIE\13,0000.shp' # 读取行政区划边界文件 boundary = gpd.read_file(boundary_path) # 遍历 2000 到 2020 年 for year in range(2000, 2021): # 生成文件地址 file_path = os.path.join(base_path, str(year), f'data_city_wgs84_{year}.shp') # 检查文件是否存在 if os.path.exists(file_path): # 读取数据 data = gpd.read_file(file_path) # 进行裁剪 clipped_data = gpd.clip(data, boundary) # 创建保存路径 output_folder = os.path.join(output_base_path, str(year)) if not os.path.exists(output_folder): os.makedirs(output_folder) # 保存裁剪后的数据 output_path = os.path.join(output_folder, f'data_city_wgs84_{year}_clipped.shp') clipped_data.to_file(output_path) ``` ### 代码解释： 1. **导入必要的库**：导入 `os` 用于文件路径操作，`geopandas` 用于地理数据处理。 2. **定义基础路径**：定义输入数据的基础路径、输出数据的基础路径和行政区划边界文件的路径。 3. **读取行政区划边界文件**：使用 `gpd.read_file` 读取行政区划边界文件。 4. **遍历年份**：使用 `range` 函数遍历 2000 到 2020 年。 5. **生成文件地址**：根据年份生成对应的文件地址。 6. **检查文件是否存在**：使用 `os.path.exists` 检查文件是否存在。 7. **读取数据并进行裁剪**：使用 `gpd.read_file` 读取数据，使用 `gpd.clip` 进行裁剪。 8. **创建保存路径**：使用 `os.makedirs` 创建保存裁剪后数据的文件夹。 9. **保存裁剪后的数据**：使用 `to_file` 方法将裁剪后的数据保存为 Shapefile 文件。