如何利用圆形缓冲和自定义函数提取栅格值？

要使用圆形缓冲和自定义函数，我们可以利用Rasterio和Shapely库。首先，加载栅格和点数据，并使用Shapely创建一个圆形缓冲区。然后，使用栅格的读取器方法和圆形缓冲区来获得其内部像素的所有值。最后，计算自定义函数并将其返回。以下是示例代码：

import rasterio
from shapely.geometry import Point
from shapely.geometry import Polygon

# 加载栅格数据
with rasterio.open("path/to/raster.tif") as raster:
    data = raster.read(1)

# 创建点并生成圆形缓冲区
point = Point(10, 10)
buffer = point.buffer(10)

# 获取缓冲区内部像素的值
window = raster.window(*buffer.bounds)
clip = raster.read(1, window=window)
transform = raster.window_transform(window)
out_shape = clip.shape
geometry = buffer.intersection(transform * rasterio.affine.Affine.from_gdal(*raster.profile["transform"]))
new_transform, width, height = rasterio.warp.transform(
    raster.crs, 'EPSG:4326', geometry=geometry, width=out_shape[1], height=out_shape[0])
bounded = rasterio.features.geometry_mask(
    [geometry], out_shape=out_shape, transform=new_transform,
    invert=True, fill_value=1)
data = clip[bounded]

# 自定义函数
def custom_function(arr):
    return sum(arr) / len(arr)

# 计算自定义函数并返回
result = custom_function(data)
print(result)