pie上将栅格数据分割成小区域分别导出，代码

在处理地理空间数据，如栅格数据（例如遥感影像），将它们分割成小区域并逐个导出，通常涉及到GIS（地理信息系统）库的操作，比如Python中的Rasterio和GeoPandas。以下是一个简单的例子，假设你已经安装了这些库：

import rasterio
import geopandas as gpd

# 打开栅格文件
with rasterio.open('input_raster.tif') as src:
    # 获取栅格的数据和元数据
    data = src.read()
    transform = src.transform
    crs = src.crs

# 创建一个方形网格，可以根据需要调整cell_size大小
cell_size = 50  # 每个小区域的边长，单位与栅格数据一致
rows, cols = data.shape
grid = (range(0, rows, cell_size), range(0, cols, cell_size))

# 将每个小区域的坐标转换为GeoDataFrame
gdf_list = []
for i in grid[0]:
    for j in grid[1]:
        window = ((i, i + cell_size), (j, j + cell_size))
        sub_data = src.read(window=window)
        sub_dataset = {
            'data': sub_data,
            'transform': transform,
            'crs': crs,
            'geometry': [src.window_to_box(window)]
        }
        sub_gdf = gpd.GeoDataFrame.from_dict(sub_dataset)
        gdf_list.append(sub_gdf)

# 对于每个子GeoDataFrame，可以单独保存到文件
for idx, gdf in enumerate(gdf_list):
    output_filename = f'sub_area_{idx}.tif'
    with rasterio.open(output_filename, 'w', driver='GTiff',
                       height=gdf['geometry'].iloc[0].bounds.height,
                       width=gdf['geometry'].iloc[0].bounds.width,
                       transform=gdf['transform'],
                       dtype=data.dtype,
                       crs=crs) as dst:
        dst.write(gdf['data'], indexes=1)