如何量化两个图像之间的差异?


178

这是我想做的:

我定期使用网络摄像头拍照。有点像延时的东西。但是,如果什么都没有真正改变,即图片看起来几乎相同,则我不想存储最新的快照。

我想有某种量化差异的方法,而我将不得不凭经验确定阈值。

我在寻找简单而不是完美。我正在使用python。


Answers:


268

大概的概念

选项1:将两个图像都加载为数组(scipy.misc.imread),然后计算逐个元素(逐像素)的差异。计算差异范数。

选项2:加载两个图像。计算每个特征向量的某些特征向量(如直方图)。计算特征向量之间的距离,而不是图像。

但是,首先要做出一些决定。

问题

您应该首先回答以下问题:

  • 图像的形状和尺寸是否相同?

    如果没有,您可能需要调整大小或裁剪它们。PIL库将帮助您使用Python做到这一点。

    如果使用相同的设置和相同的设备拍摄它们,则它们可能是相同的。

  • 图像是否对齐正确?

    如果不是,则可能要先运行互相关,然后首先找到最佳对齐方式。SciPy具有执行此功能的功能。

    如果相机和场景静止不动,则图像可能对齐良好。

  • 图像的曝光总是一样吗?(亮度/对比度是否相同?)

    如果不是,则可能要规范化图像。

    但是要小心,在某些情况下,这样做可能弊大于利。例如,深色背景上的单个明亮像素将使标准化图像非常不同。

  • 颜色信息重要吗?

    如果要注意颜色变化,则将具有每个点的颜色值向量,而不是灰度图像中的标量值。编写此类代码时,您需要更多注意。

  • 图像中是否有明显的边缘?他们可能会移动吗?

    如果是,则可以先应用边缘检测算法(例如,使用Sobel或Prewitt变换计算梯度,应用一些阈值),然后将第一个图像上的边缘与第二个图像上的边缘进行比较。

  • 图像中是否有噪点?

    所有传感器都会在一定程度上污染图像。低成本传感器噪声更大。您可能希望在比较图像之前进行一些降噪处理。在这里,模糊是最简单(但不是最好)的方法。

  • 您想注意哪些变化?

    这可能会影响要用于图像之间差异的标准的选择。

    考虑使用曼哈顿范数(绝对值的总和)或零范数(元素数量不等于零)来衡量图像已更改了多少。前者将告诉您图像偏移了多少,后者将仅告诉您有多少像素不同。

我认为您的图像对齐良好,大小和形状相同,可能具有不同的曝光度。为简单起见,即使它们是彩色(RGB)图像,我也将它们转换为灰度。

您将需要这些导入:

import sys

from scipy.misc import imread
from scipy.linalg import norm
from scipy import sum, average

主要功能,读取两张图像,转换为灰度,比较并打印结果:

def main():
    file1, file2 = sys.argv[1:1+2]
    # read images as 2D arrays (convert to grayscale for simplicity)
    img1 = to_grayscale(imread(file1).astype(float))
    img2 = to_grayscale(imread(file2).astype(float))
    # compare
    n_m, n_0 = compare_images(img1, img2)
    print "Manhattan norm:", n_m, "/ per pixel:", n_m/img1.size
    print "Zero norm:", n_0, "/ per pixel:", n_0*1.0/img1.size

如何比较。img1img2是2D SciPy的阵列,在这里:

def compare_images(img1, img2):
    # normalize to compensate for exposure difference, this may be unnecessary
    # consider disabling it
    img1 = normalize(img1)
    img2 = normalize(img2)
    # calculate the difference and its norms
    diff = img1 - img2  # elementwise for scipy arrays
    m_norm = sum(abs(diff))  # Manhattan norm
    z_norm = norm(diff.ravel(), 0)  # Zero norm
    return (m_norm, z_norm)

如果文件是彩色图像,则imread返回一个3D数组,平均RGB通道(最后一个数组轴)以获取强度。对于灰度图像(例如.pgm)无需这样做:

def to_grayscale(arr):
    "If arr is a color image (3D array), convert it to grayscale (2D array)."
    if len(arr.shape) == 3:
        return average(arr, -1)  # average over the last axis (color channels)
    else:
        return arr

标准化很简单,您可以选择标准化为[0,1]而不是[0,255]。arr是一个SciPy数组,因此所有操作都是按元素进行的:

def normalize(arr):
    rng = arr.max()-arr.min()
    amin = arr.min()
    return (arr-amin)*255/rng

运行main函数:

if __name__ == "__main__":
    main()

现在,您可以将所有内容放入脚本中并针对两个图像运行。如果我们将图像与其自身进行比较,则没有区别:

$ python compare.py one.jpg one.jpg
Manhattan norm: 0.0 / per pixel: 0.0
Zero norm: 0 / per pixel: 0.0

如果我们模糊图像并与原始图像进行比较,则存在一些差异:

$ python compare.py one.jpg one-blurred.jpg 
Manhattan norm: 92605183.67 / per pixel: 13.4210411116
Zero norm: 6900000 / per pixel: 1.0

PS整个compare.py脚本。

更新:相关技术

由于问题是关于视频序列的,其中的帧可能几乎是相同的,并且您在寻找不寻常的东西,所以我想提及一些可能相关的替代方法:

  • 背景扣除和分割(以检测前景对象)
  • 稀疏光流(检测运动)
  • 比较直方图或其他统计信息而不是图像

我强烈建议您阅读“学习OpenCV”这本书,第9章(图像部分和分割)和第10章(跟踪和运动)。前者教导使用背景减法,后者给出有关光流方法的一些信息。所有方法都在OpenCV库中实现。如果使用Python,我建议使用OpenCV≥2.3及其cv2Python模块。

最简单的背景减法版本:

  • 了解背景每个像素的平均值μ和标准偏差σ
  • 比较当前像素值到(μ-2σ,μ+2σ)或(μ-σ,μ+σ)的范围

更高级的版本会考虑每个像素的时间序列,并处理非静态场景(例如移动的树木或草地)。

光流的概念是拍摄两个或更多帧,然后将速度矢量分配给每个像素(密集光流)或分配给其中一些像素(稀疏光流)。要估计稀疏的光流,可以使用Lucas-Kanade方法(它也在OpenCV中实现)。显然,如果有很多流量(在速度场的最大值上具有较高的平均值),则说明帧中正在移动某些内容,并且后续图像会有所不同。

比较直方图可能有助于检测连续帧之间的突然变化。Courbon等人,2010年使用了这种方法:

连续帧的相似性。测量两个连续帧之间的距离。如果过高,则意味着第二帧已损坏,因此图像被消除。所述的Kullback-Leibler距离,或互熵,对两帧的直方图:

$$ d(p,q)= \ sum_i p(i)\ log(p(i)/ q(i))$$

其中pq是帧的直方图。阈值固定为0.2。


RuntimeWarning: invalid value encountered in double_scalars在第44行(return (arr-amin)*255/rng)和ValueError: array must not contain infs or NaNs第30行(z_norm = norm(diff.ravel(), 0)
BioGeek

@BioGeek rng等于零。只需添加一张支票并进行设置rng = 1
haisi 18'Apr 24'18

76

一个简单的解决方案:

将图像编码为jpeg,并在filesize中寻找实质性的变化。

我已经用视频缩略图实现了类似的功能,并且取得了很多成功和可扩展性。


3
这是一个非常容易,简单的解决方案,比任何按像素进行的比较要好得多。如果网络摄像头的图像中有一点点噪点,或者图像偏移了一个像素,那么直接比较就会发现所有这些毫无意义的变化。一种更可靠的方法是计算离散余弦变换,然后在频域中比较图像。像这样使用JPEG压缩,将获得大多数收益,而无需深入研究傅里叶理论。
AndrewF

喜欢它。尽管其他解决方案也可以使用,但这在常见情况下具有很大的优势:如果您不想保存“基本”图像怎么办?只需将文件大小另存为哈希,然后仅将数字与减号进行比较即可。就我而言,我有4张图片,其中一张非常相似,其他3张完全不同。只需缩放到相同的尺寸,以jpg和减。非常好。
迭戈·安德烈斯·迪亚斯·埃斯皮诺萨

60

您可以使用PIL的功能比较两个图像。

import Image
import ImageChops

im1 = Image.open("splash.png")
im2 = Image.open("splash2.png")

diff = ImageChops.difference(im2, im1)

diff对象是一个图像,其中每个像素都是从第一图像减去第二图像中该像素的颜色值的结果。使用差异图像,您可以做几件事。最简单的一种是diff.getbbox()功能。它将告诉您包含两个图像之间所有更改的最小矩形。

您可能也可以使用PIL中的函数来实现此处提到的其他内容的近似值。


2
我要保存差异图像。表示保存图像差异的差异对象。ow,我要保存吗?
萨加尔2014年

2
@Anthony,您可以在diff对象上调用save()来指定图像名称。像这样:diff.save(“ diff.png”)它将为您保存差异图像。
萨加尔(Sagar)2015年

20

两种流行且相对简单的方法是:(a)已经建议的欧几里得距离,或(b)归一化互相关。与简单的互相关相比,归一化互相关趋向于在光照变化方面更加明显。维基百科给出了归一化互相关的公式。也存在更复杂的方法,但是它们需要大量工作。

使用类似numpy的语法,

dist_euclidean = sqrt(sum(((i1-i2)^ 2))/ i1.size

dist_manhattan = sum(abs(i1-i2))/ i1.size

dist_ncc = sum((i1-平均值(i1))*(i2-平均值(i2)))/(
  (i1.size-1)* stdev(i1)* stdev(i2))

假设i1i2是2D灰度图像阵列。


3
图像互相关函数内置于SciPy(docs.scipy.org/doc/scipy/reference/Generated/…)中,而使用FFT的快速版本在stsci python中可用(stsci.edu/resources/software_hardware/pyraf/ stsci_python
endolith

14

尝试的琐碎事情:

将两个图像重新采样为较小的缩略图(例如64 x 64),并将缩略图与特定阈值逐像素进行比较。如果原始图像几乎相同,则重新采样的缩略图将非常相似甚至完全相同。此方法要注意特别是在弱光场景中可能发生的噪声。如果使用灰度,可能会更好。


但是您如何比较像素?
航空母舰

有了缩略图后,您可以简单地将像素一一比较。如果您使用彩色,则将计算RGB值的“距离”,如果使用灰度,则将仅计算灰度之间的差异。
Ates Goral

1
“一一比较像素”。那是什么意思?如果每像素64 ^ 2像素测试之一失败,则测试是否会失败?
Federico A. Ramponi's

我的意思是“将缩略图与特定阈值逐像素进行比较”,是想出一种模糊算法来比较像素。如果计算出的差异(取决于模糊算法)超过某个阈值,则图像“不相同”。
Ates Goral

1
非常简单的例子,没有了“模糊算法”:通过每个像素并行循环(比较像素#Ñ图像#1到像素#Ñ图像#2),并且在值添加的差异的可变
MK12

7

我专门解决的问题是如何计算它们是否“足够不同”。我假设您可以弄清楚如何一一减去像素。

首先,我将拍摄一束没有任何变化的图像,并找出仅由于捕获的变化,成像系统中的噪声,JPEG压缩伪像以及照明的瞬时变化而导致的任何像素变化的最大量。 。也许您会发现,即使什么都没有发生,也将期待1或2位的差异。

然后,对于“真实”测试,您需要这样的条件:

  • 如果最多P个像素相差不超过E,则相同。

因此,也许,如果E = 0.02,P = 1000,则(大约)意味着,如果任何单个像素变化超过约5个单位(假设为8位图像),或者如果超过1000个,它将是“不同的”像素根本没有任何错误。

这主要是作为一种良好的“分类”技术,可以快速识别足够接近不需要进一步检查的图像。然后,“失败”的图像可能更多地是一种更为复杂/昂贵的技术,例如,如果相机晃动了一点,或者对照明的变化更鲁棒,则不会出现误报。

我运行一个开放源代码项目OpenImageIO,其中包含一个名为“ idiff”的实用程序,该实用程序将差异与类似的阈值进行比较(实际上更加精细)。即使您不想使用此软件,也可能要查看源代码来了解我们是如何做到的。它在商业上已经大量使用,并且开发了此阈值技术,以便我们可以拥有一个用于渲染和图像处理软件的测试套件,其中的“参考图像”可能与平台之间或平台之间的细微差异有所不同。 tha算法,因此我们希望进行“公差内匹配”操作。


6

我在工作中遇到了类似的问题,我正在重写图像转换端点,并且想检查新版本与旧版本产生的输出相同或几乎相同。所以我这样写:

https://github.com/nicolashahn/diffimg

它在相同大小的图像上并且在每个像素级别上运行,测量每个通道上的值差:R,G,B(,A),取这些通道的平均差,然后对所有像素,并返回比率。

例如,对于一个10x10的白色像素图像,而同一图像但一个像素变为红色,则该像素处的差异为1/3或0.33 ...(RGB 0,0,0与255,0,0 ),其他所有像素均为0。如果总共有100个像素,则0.33 ... / 100 =图像差异约为0.33%。

我相信这对于OP的项目将是完美的(我意识到这是一个非常老的帖子,但是为将来想要也比较python中的图像的StackOverflowers发布)。


5

给出的大多数答案都不会涉及照明水平。

在进行比较之前,我首先将图像标准化为标准亮度。


如果您要拍摄周期性图像并扩散相邻的图像对,则可以在有人打开灯后保留第一个图像。
walkytalky

5

衡量两个图像之间相似度的另一种不错的简单方法:

import sys
from skimage.measure import compare_ssim
from skimage.transform import resize
from scipy.ndimage import imread

# get two images - resize both to 1024 x 1024
img_a = resize(imread(sys.argv[1]), (2**10, 2**10))
img_b = resize(imread(sys.argv[2]), (2**10, 2**10))

# score: {-1:1} measure of the structural similarity between the images
score, diff = compare_ssim(img_a, img_b, full=True)
print(score)

如果其他人对比较图像相似性的更强大方法感兴趣,我会组合一个教程和Web 应用程序,以使用Tensorflow测量和可视化相似图像。


3
是的,skimage非常适合用于此应用程序。我用from skimage.measure import compare_ssim, compare_mse很多。skimage.measure文档
ximiki '18 -10-11

3

您是否看到过寻找相似图像算法问题?检查一下以查看建议。

我建议对您的帧进行小波变换(我使用Haar变换为此编写了C扩展);然后,比较两个图片之间最大(按比例)小波因子的索引,您应该得到数值相似度近似值。


2

如果为时已晚,我对此表示歉意,但是由于我一直在做类似的事情,所以我认为我可以有所作为。

也许使用OpenCV可以使用模板匹配。假设您正在使用网络摄像头,如您所说:

  1. 简化图像(可能是阈值吗?)
  2. 应用模板匹配并使用minMaxLoc检查max_val

提示:max_val(或min_val取决于所使用的方法)将为您提供数字,大数字。要获得百分比差异,请对同一张图片使用模板匹配-结果将是您的100%。

伪代码举例说明:

previous_screenshot = ...
current_screenshot = ...

# simplify both images somehow

# get the 100% corresponding value
res = matchTemplate(previous_screenshot, previous_screenshot, TM_CCOEFF)
_, hundred_p_val, _, _ = minMaxLoc(res)

# hundred_p_val is now the 100%

res = matchTemplate(previous_screenshot, current_screenshot, TM_CCOEFF)
_, max_val, _, _ = minMaxLoc(res)

difference_percentage = max_val / hundred_p_val

# the tolerance is now up to you

希望能帮助到你。


1

推土机距离可能正是您需要的。但是,实时实施可能会有点繁重。


我对这个答案的回答并不满意:“我在寻找简单而不是完美。我在使用python。”
-PilouPili

我认为随着这个问题线程的大量访问,吸引大多数观众的标题是关于如何量化两个图像之间的差异,因此在这里很有用。
Danoram

1

如何计算两个图像的曼哈顿距离。这样就可以得到n * n个值。然后,您可以执行类似行平均的操作以减少到n个值,并执行一个函数以获得一个单个值。


1

我在三脚架上用同一台相机拍摄的jpg图像的运气很不错,方法是:(1)大大简化(例如从3000像素宽变为100像素或更小)(2)将每个jpg数组展平为一个向量(3)使用简单的相关算法对序列图像进行成对相关以获取相关系数(4)对相关系数进行平方以得到r平方(即,一幅图像中的可变性分数由下一幅图像的变化解释)(5)如果r-square <0.9,则表示两个图像不同,并且两者之间发生了某些情况。

在我的实施中,这是强大且快速的(Mathematica 7)

值得一试的是,将您感兴趣的图像部分并通过将所有图像裁剪到该小区域来专注于该部分,否则将错过与相机遥远但重要的更改。

我不知道如何使用Python,但是可以确定它也具有相关性,不是吗?


1

您可以计算两个图像的直方图,然后计算Bhattacharyya系数,这是一种非常快速的算法,我用它来检测板球视频中的镜头变化(在C中使用openCV)


您可以计算图像本身的系数吗?
endlith 2011年

您将必须计算图像的直方图(根据要求使用直方图的bin大小)。
vishalv2050 2011年

1

查看isk-daemon如何实现Haar Wavelets 。您可以使用它的imgdb C ++代码即时计算图像之间的差异:

isk-daemon是一个开源数据库服务器,能够将基于内容的(可视)图像搜索添加到任何与图像相关的网站或软件中。

这项技术允许任何与图像相关的网站或软件的用户在小部件上草绘他们想要查找的图像,并让网站将最相似的图像回复给他们,或者在每个图像详细信息页面上简单地请求更多相似的照片。


1

我遇到了同样的问题,并编写了一个简单的python模块,该模块使用枕头的ImageChops比较两个相同尺寸的图像,以创建黑白差异图像,并对直方图值求和。

您可以直接获得此分数,也可以得到与全黑对比白对比的百分比值。

它还包含一个简单的is_equal函数,可以在图像通过(并包括)相等的情况下提供模糊阈值。

该方法不是很复杂,但是可能对于其他在同一问题上苦苦挣扎的人有用。

https://pypi.python.org/pypi/imgcompare/


1

一种更有原则的方法是使用全局描述符来比较图像,例如GIST或CENTRIST。散列函数,如所描述这里,还提供了类似的解决方案。


1
import os
from PIL import Image
from PIL import ImageFile
import imagehash
  
#just use to the size diferent picture
def compare_image(img_file1, img_file2):
    if img_file1 == img_file2:
        return True
    fp1 = open(img_file1, 'rb')
    fp2 = open(img_file2, 'rb')

    img1 = Image.open(fp1)
    img2 = Image.open(fp2)

    ImageFile.LOAD_TRUNCATED_IMAGES = True
    b = img1 == img2

    fp1.close()
    fp2.close()

    return b





#through picturu hash to compare
def get_hash_dict(dir):
    hash_dict = {}
    image_quantity = 0
    for _, _, files in os.walk(dir):
        for i, fileName in enumerate(files):
            with open(dir + fileName, 'rb') as fp:
                hash_dict[dir + fileName] = imagehash.average_hash(Image.open(fp))
                image_quantity += 1

    return hash_dict, image_quantity

def compare_image_with_hash(image_file_name_1, image_file_name_2, max_dif=0):
    """
    max_dif: The maximum hash difference is allowed, the smaller and more accurate, the minimum is 0.
    recommend to use
    """
    ImageFile.LOAD_TRUNCATED_IMAGES = True
    hash_1 = None
    hash_2 = None
    with open(image_file_name_1, 'rb') as fp:
        hash_1 = imagehash.average_hash(Image.open(fp))
    with open(image_file_name_2, 'rb') as fp:
        hash_2 = imagehash.average_hash(Image.open(fp))
    dif = hash_1 - hash_2
    if dif < 0:
        dif = -dif
    if dif <= max_dif:
        return True
    else:
        return False


def compare_image_dir_with_hash(dir_1, dir_2, max_dif=0):
    """
    max_dif: The maximum hash difference is allowed, the smaller and more accurate, the minimum is 0.

    """
    ImageFile.LOAD_TRUNCATED_IMAGES = True
    hash_dict_1, image_quantity_1 = get_hash_dict(dir_1)
    hash_dict_2, image_quantity_2 = get_hash_dict(dir_2)

    if image_quantity_1 > image_quantity_2:
        tmp = image_quantity_1
        image_quantity_1 = image_quantity_2
        image_quantity_2 = tmp

        tmp = hash_dict_1
        hash_dict_1 = hash_dict_2
        hash_dict_2 = tmp

    result_dict = {}

    for k in hash_dict_1.keys():
        result_dict[k] = None

    for dif_i in range(0, max_dif + 1):
        have_none = False

        for k_1 in result_dict.keys():
            if result_dict.get(k_1) is None:
                have_none = True

        if not have_none:
            return result_dict

        for k_1, v_1 in hash_dict_1.items():
            for k_2, v_2 in hash_dict_2.items():
                sub = (v_1 - v_2)
                if sub < 0:
                    sub = -sub
                if sub == dif_i and result_dict.get(k_1) is None:
                    result_dict[k_1] = k_2
                    break
    return result_dict


def main():
    print(compare_image('image1\\815.jpg', 'image2\\5.jpg'))
    print(compare_image_with_hash('image1\\815.jpg', 'image2\\5.jpg', 7))
    r = compare_image_dir_with_hash('image1\\', 'image2\\', 10)
    for k in r.keys():
        print(k, r.get(k))


if __name__ == '__main__':
    main()
  • 输出:



    image2 \ 5.jpg image1 \ 815.jpg
    image2 \ 6.jpg image1 \ 819.jpg
    image2 \ 7.jpg image1 \ 900.jpg
    image2 \ 8.jpg image1 \ 998.jpg
    image2 \ 9.jpg image1 \ 1012 .jpg

  • 示例图片:

    • 815.jpg
      815.jpg

    • 5.jpg
      5.jpg


0

我认为您可以简单地计算两个图像的亮度之间的欧几里德距离(即sqrt(差异平方和,逐个像素)),并且如果这低于某个经验阈值,则认为它们相等。而且,您最好将其包装为C函数。


0

有很多指标可用于评估两个图像的外观/外观。

在这里我将不讨论任何代码,因为我认为这应该是科学问题,而不是技术问题。

通常,问题与人对图像的感知有关,因此每种算法都对人的视觉系统特征有支持。

经典方法是:

可见差异预测器:一种评估图像保真度的算法(https://www.spiedigitallibrary.org/conference-proceedings-of-spie/1666/0000/Visible-differences-predictor--an-algorithm-for-the-评估/10.1117/12.135952.short?SSO=1

图像质量评估:从错误可见性到结构相似性(http://www.cns.nyu.edu/pub/lcv/wang03-reprint.pdf

FSIM:用于图像质量评估的特征相似性索引(https://www4.comp.polyu.edu.hk/~cslzhang/IQA/TIP_IQA_FSIM.pdf

其中,SSIM(图像质量评估:从错误可见性到结构相似性)是最容易计算的,并且开销也很小,如另一篇论文《基于梯度相似性的图像质量评估》(https://www.semanticscholar .org / paper /基于渐变效果的图像质量评估-刘林/ 2b819bef80c02d5d4cb56f27b202535e119df988)。

还有许多其他方法。如果您对艺术真正感兴趣,可以看看Google学术搜索并搜索诸如“视觉差异”,“图像质量评估”之类的内容。


0

有一种使用numpy的简单快速的解决方案,可通过计算均方误差来解决:

before = np.array(get_picture())
while True:
    now = np.array(get_picture())
    MSE = np.mean((now - before)**2)

    if  MSE > threshold:
        break

    before = now
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.