计算原始图像的拉普拉斯金字塔（LoG)

LoG（Laplacian of Gaussian）是一种图像增强算法，可以用于边缘检测、纹理分析等。下面是使用Python计算原始图像的LoG金字塔的示例代码：

import cv2
import numpy as np

# 读取原始图像
img = cv2.imread('lena.jpg', cv2.IMREAD_GRAYSCALE)

# 定义高斯金字塔参数
levels = 5
ksize = 5
sigma = 1.0

# 计算高斯金字塔
g = img.copy()
gp = [g]
for i in range(levels):
    g = cv2.GaussianBlur(g, (ksize, ksize), sigma)
    gp.append(g)

# 定义拉普拉斯金字塔参数
lp = [gp[levels-1]]
for i in range(levels-1, 0, -1):
    GE = cv2.pyrUp(gp[i])
    L = cv2.subtract(gp[i-1], GE)
    lp.append(L)

# 显示结果
for i in range(levels):
    cv2.imshow('Gaussian Pyramid ' + str(i+1), gp[i])
for i in range(levels):
    cv2.imshow('Laplacian Pyramid ' + str(i+1), lp[i])
cv2.waitKey(0)
cv2.destroyAllWindows()

在这个示例中，我们首先读取了一张灰度图像，并定义了高斯金字塔的参数。接着，我们计算了高斯金字塔，然后利用高斯金字塔计算了拉普拉斯金字塔。最后，我们将结果显示出来。

请注意，在计算拉普拉斯金字塔时，我们先从高斯金字塔的最高层开始，利用cv2.subtract()函数计算当前层和上一层金字塔的高斯图像的差值，得到当前层的拉普拉斯图像。然后将拉普拉斯图像存储到列表中，作为下一层金字塔的输入。

通过计算拉普拉斯金字塔，我们可以检测到图像中的边缘和纹理等特征。