能否提供一个实现三层拉普拉斯金字塔算法的详细代码示例？

三层拉普拉斯金字塔（Laplacian Pyramid）通常用于图像处理，特别是图像的分级表示和缩放。它通过将原始图像分解成一系列低分辨率版本和高频细节来构建。以下是Python中使用numpy库实现一个简单三层拉普拉斯金字塔的一个简略示例：

import numpy as np

def laplacian_pyramid(image, levels=3):
    # 创建一个空的金字塔列表
    pyramid = [image]
    
    for _ in range(1, levels + 1):
        # 对当前层进行低通滤波（这里假设使用的是简单的平均滤波）
        blurred = np.mean(pyramid[-1], axis=(0, 1), keepdims=True)
        
        # 计算差分，得到高斯金字塔
        detail = pyramid[-1] - blurred
        
        # 把结果添加到金字塔中
        pyramid.append(blurred)
        pyramid.append(detail)
    
    return pyramid[:levels+1]

# 示例用法
input_image = np.random.rand(500, 500, 3)  # 500x500 RGB图像
pyramid = laplacian_pyramid(input_image, levels=3)