python - 使用 K-Nearest 聚类颜色的有效方法

Question

我正在尝试将图像上的颜色聚类到预定义的类别(黑色、白色、蓝色、绿色、红色)。我正在使用以下代码:

import numpy as np
import cv2

src = cv2.imread('objects.png')

colors = np.array([[0x00, 0x00, 0x00],
                   [0xff, 0xff, 0xff],
                   [0xff, 0x00, 0x00],
                   [0x00, 0xff, 0x00],
                   [0x00, 0x00, 0xff]], dtype=np.float32)
classes = np.array([[0], [1], [2], [3], [4]], np.float32)
dst = np.zeros(src.shape, np.float32)

knn = cv2.KNearest()
knn.train(colors, classes)

# This loop is very inefficient!
for i in range(0, src.shape[0]):
    for j in range(0, src.shape[1]):
        sample = np.reshape(src[i,j], (-1,3)).astype(np.float32)
        retval, result, neighbors, dist = knn.find_nearest(sample, 1)
        dst[i,j] = colors[result[0,0]]

cv2.imshow('src', src)
cv2.imshow('dst', dst)
cv2.waitKey()

代码运行良好，结果如下所示。左边的图像是输入，右边的图像是输出。

然而，上面的循环非常低效，使转换变慢。替换上面循环的最有效的 Numpy 操作是什么？

Answer 1

如果您想要一个简单的平方差度量(“这是欧几里德最接近的数字”)，这会起作用。

计算差异

diff = ((src[:,:,:,None] - colors.T)**2).sum(axis=2)

(假设 src 的形状是 y,x,3)

选择最接近的颜色索引:

index = diff.argmin(axis=2)

新图片:

out = colors[index]

如果你的颜色真的要有 0 或 0xff 的分量值，你可以使用类似的东西

out = np.where(src>0x88, 0xff, 0)