我正在尝试使用 opencv 函数获得拉普拉斯金字塔pyrUp
的一层:和pyrDown
.
在这本书的文档和更详细的内容中,我发现第i 个拉普拉斯层应该通过以下表达式获得:
Li = Gi - pyrDown(Gi+1)
其中 Gi 是高斯金字塔的第 i 层。
我试过的是:
def get_laplacian_pyramid_layer(img, n):
gi = img
for i in range(n):
gi_prev = gi
gi = cv2.pyrDown(gi_prev)
pyrup = cv2.pyrUp(gi)
return cv2.addWeighted(gi_prev, 1.5, pyrup, -0.5, 0)
但是我得到了减法中涉及的不同大小的图像。我不明白,因为pyrUp
它应该反转高斯金字塔的过程,即 pyrDown (当然会丢失信息,但这不应该影响大小,对吧?)。
更新
我将我的代码重构为:
def get_laplacian_pyramid_layer(img, n):
'''Returns the n-th layer of the laplacian pyramid'''
currImg, i = img, 0
while i < n: # and currImg.size > max_level (83)
down, up = new_empty_img(img.shape), new_empty_img(img.shape)
down = cv2.pyrDown(img)
up = cv2.pyrUp(down, dstsize=currImg.shape)
lap = currImg - up
currImg = down
i += 1
return lap
如您所见,我通过函数的参数强制目标图像的大小与源图像的大小dstsize
相同pyrUp
。
pyrUp
但是,在执行函数时,这段代码也会给我一个错误。错误消息是:
OpenCV 错误:断言失败 (std::abs(dsize.width - ssize.width*2) == dsize.width % 2 && std::abs(dsize.height - ssize.height*2) == dsize.height % 2)在pyrUp_中,
在调试模式下,我检查了断言的表达式:
up.shape[1]-down.shape[1]*2 == up.shape[1] %2 and up.shape[0]-down.shape[0]*2 == up.shape[0] %2
它很满意。
所以,我不知道发生了什么。