24

我正在尝试使用 NumPy 和 PIL 将两个图像添加在一起。我在MATLAB中执行此操作的方式类似于：

>> M1 = imread('_1.jpg');
>> M2 = imread('_2.jpg');
>> resM = M1 + M2;
>> imwrite(resM, 'res.jpg');

我得到这样的东西：

替代文字 http://www.deadlink.cc/matlab.jpg

使用合成程序并添加图像，MATLAB 结果似乎是正确的。

在 Python 中，我试图做同样的事情：

from PIL import Image
from numpy import *

im1 = Image.open('/Users/rem7/Desktop/_1.jpg')
im2 = Image.open('/Users/rem7/Desktop/_2.jpg')

im1arr = asarray(im1)
im2arr = asarray(im2)

addition = im1arr + im2arr

resultImage = Image.fromarray(addition)
resultImage.save('/Users/rem7/Desktop/a.jpg')

我得到这样的东西：

替代文字 http://www.deadlink.cc/python.jpg

为什么我得到所有这些时髦的颜色？我也尝试过使用ImageMath.eval("a+b", a=im1, b=im2)，但我收到有关 RGB 不支持的错误。

我还看到有一个Image.blend()但需要 alpha。

实现我想要的东西的最佳方法是什么？

源图像（图像已被删除）：

替代文字 http://www.deadlink.cc/_1.jpg 替代文字 http://www.deadlink.cc/_2.jpg

嗯，好的，我使用添加图像图标添加了源图像，它们在我编辑帖子时出现，但由于某种原因，图像没有出现在帖子中。

（图片已被删除）2013 05 09

4

5 回答 5

34

正如每个人已经建议的那样，您观察到的奇怪颜色已经溢出。正如您在schnaader 回答的评论中指出的那样，如果您添加这样的图像，您仍然会溢出：

addition=(im1arr+im2arr)/2

此溢出的原因是您的 NumPy 数组 ( im1arr im2arr ) 属于uint8类型（即 8 位）。这意味着数组的每个元素最多只能保存 255 的值，因此当总和超过 255 时，它会在 0 附近循环：

>>>array([255,10,100],dtype='uint8') +  array([1,10,160],dtype='uint8')
array([ 0, 20,  4], dtype=uint8)

为避免溢出，您的数组应该能够包含超过 255 的值。例如，您需要将它们转换为浮点数，执行混合操作并将结果转换回 uint8：

im1arrF = im1arr.astype('float')
im2arrF = im2arr.astype('float')
additionF = (im1arrF+im2arrF)/2
addition = additionF.astype('uint8')

你不应该这样做：

addition = im1arr/2 + im2arr/2

当您丢失信息时，通过在执行混合信息之前压缩图像的动态（您实际上使图像为 7 位）。

MATLAB 注意：您在 MATLAB 中看不到此问题的原因可能是因为 MATLAB 在其函数之一中隐式处理了溢出。

于 2009-02-08T17:05:45.870 回答
20

使用 alpha 值为 0.5 的 PIL 的 blend() 将等效于 (im1arr + im2arr)/2。Blend 不要求图像具有 Alpha 层。

试试这个：

from PIL import Image
im1 = Image.open('/Users/rem7/Desktop/_1.jpg')
im2 = Image.open('/Users/rem7/Desktop/_2.jpg')
Image.blend(im1,im2,0.5).save('/Users/rem7/Desktop/a.jpg')
于 2009-02-08T03:58:33.760 回答
2

您发布的代码似乎只是总结了大于 256 的值和值溢出。你想要“（a + b）/ 2”或“min（a + b，256）”之类的东西。后者似乎是您的 Matlab 示例的方式。

于 2009-02-08T01:28:27.573 回答
2

要钳制 numpy 数组值：

>>> c = a + b
>>> c[c > 256] = 256
于 2009-02-08T09:54:53.260 回答
0

您的示例图像没有显示在我的表格中，所以我要做一些猜测。

我不记得确切的 numpy 到 pil 转换是如何工作的，但有两种可能的情况。我 95% 确定它是 1，但我给出 2 以防万一我错了。1) 1 im1Arr 是一个 MxN 整数数组 (ARGB)，当您将 im1arr 和 im2arr 加在一起时，如果组件 b1+b2>255，您将从一个通道溢出到下一个通道。我猜 matlab 将他们的图像表示为 MxNx3 数组，因此每个颜色通道都是独立的。您可以通过拆分 PIL 图像通道然后制作 numpy 数组来解决此问题

2) 1 im1Arr 是一个 MxNx3 字节数组，当您将 im1arr 和 im2arr 添加在一起时，您正在包装组件。

在显示之前，您还必须将范围重新调整到 0-255 之间。您的选择是除以 2、缩放 255/array.max() 或进行剪辑。我不知道matlab是做什么的

于 2009-02-08T01:32:40.043 回答