Python PIL函数用于分割混合两个图像？

Question

编辑：多亏了Mark和zephyr，代码现在可以工作了。zephyr也有以下两种替代工作解决方案

我想用PIL分割并混合两个图像。我找到了ImageChops.multiply(image1, image2)，但找不到类似的divide(image, image2)函数

Divide Blend Mode Explained（我在这里使用前两个图像作为测试源。）

是否有我错过的内置divide blend函数（PIL或其他）

我下面的测试代码正在运行，并且正在接近我想要的。结果图像输出类似于此处的divide blend示例图像：Divide Blend Mode Explained

是否有一种更有效的方法来执行此分割混合操作（步骤更少，速度更快）？起初，我尝试在Image.eval和ImageMath.eval中使用lambda函数来检查黑色像素，并在除法过程中将它们翻转为白色，但我无法得到两者中的任何一个来产生正确的结果

编辑：修复了代码，由于马克和西风，代码缩短了。生成的图像输出与下面zephyr的numpy和scipy解决方案的输出相匹配

# PIL Divide Blend test

import Image, os, ImageMath

imgA = Image.open('01background.jpg')
imgA.load()
imgB = Image.open('02testgray.jpg')
imgB.load()

# split RGB images into 3 channels
rA, gA, bA = imgA.split()
rB, gB, bB = imgB.split()

# divide each channel (image1/image2)
rTmp = ImageMath.eval("int(a/((float(b)+1)/256))", a=rA, b=rB).convert('L')
gTmp = ImageMath.eval("int(a/((float(b)+1)/256))", a=gA, b=gB).convert('L')
bTmp = ImageMath.eval("int(a/((float(b)+1)/256))", a=bA, b=bB).convert('L')

# merge channels into RGB image
imgOut = Image.merge("RGB", (rTmp, gTmp, bTmp))

imgOut.save('PILdiv0.png', 'PNG')

os.system('start PILdiv0.png')

Answer 1

这里有一个除法函数的数学定义： http://www.linuxtopia.org/online_books/graphics_tools/gimp_advanced_guide/gimp_guide_node55_002.html

下面是一个使用scipy/matplotlib的实现：

import numpy as np
import scipy.misc as mpl

a = mpl.imread('01background.jpg')
b = mpl.imread('02testgray.jpg')

c = a/((b.astype('float')+1)/256)
d = c*(c < 255)+255*np.ones(np.shape(c))*(c > 255)

e = d.astype('uint8')

mpl.imshow(e)
mpl.imsave('output.png', e)

如果不想使用matplotlib，可以这样做（我假设您有numpy）：

imgA = Image.open('01background.jpg')
imgA.load()
imgB = Image.open('02testgray.jpg')
imgB.load()

a = asarray(imgA)
b = asarray(imgB)
c = a/((b.astype('float')+1)/256)
d = c*(c &lt 255)+255*ones(shape(c))*(c &gt 255)
e = d.astype('uint8')

imgOut = Image.fromarray(e)
imgOut.save('PILdiv0.png', 'PNG')