PythonPIL.Image.convert文件不是用最近的调色板替换颜色。

Question

这是一个后续问题：Convert image to specific palette using PIL without dithering

我也想创建一个脚本，它可以将图像转换为一组特定的颜色，而不必进行抖动。在

我已经实现了“自定义量化”功能，作为问题的答案。除了一个大问题外，大多数脚本都运行良好。在

浅绿色的RGB（130190,40）被浅棕色的RGB（166、141、95）代替。（见鬃毛左上角的浅绿色。）

from PIL import Image

def customConvert(silf, palette, dither=False):
    ''' Convert an RGB or L mode image to use a given P image's palette.
        PIL.Image.quantize() forces dither = 1. 
        This custom quantize function will force it to 0.
        https://stackoverflow.com/questions/29433243/convert-image-to-specific-palette-using-pil-without-dithering
    '''

    silf.load()

    # use palette from reference image made below
    palette.load()
    im = silf.im.convert("P", 0, palette.im)
    # the 0 above means turn OFF dithering making solid colors
    return silf._new(im)

palette = [ 
    0,0,0,
    0,0,255,
    15,29,15,
    26,141,52,
    41,41,41,
    65,105,225,
    85,11,18,
    128,0,128,
    135,206,236,
    144,238,144,
    159,30,81,
    165,42,42,
    166,141,95,
    169,169,169,
    173,216,230,
    211,211,211,
    230,208,122,
    245,245,220,
    247,214,193,
    255,0,0,
    255,165,0,
    255,192,203,
    255,255,0,
    255,255,255
    ] + [0,] * 232 * 3


# a palette image to use for quant
paletteImage = Image.new('P', (1, 1), 0)
paletteImage.putpalette(palette)


# open the source image
imageOrginal = Image.open('lion.png').convert('RGB')

# convert it using our palette image
imageCustomConvert = customConvert(imageOrginal, paletteImage, dither=False).convert('RGB')

CIE76三角洲-E：

当前：RGB（130190,40）-->；RGB（166，141，95）=57.5522

预期值：RGB（130190,40）-->；RGB（144238144）=31.5623

---

有人能解释一下如果我写的代码不正确或建议如何让它工作。在

Answer 1

我试着计算每个像素的CIE76 Delta-E函数，以得到最接近的颜色。Python不是我最好的语言，所以如果它能像您期望的那样工作，那么您可能需要再问一个问题来优化代码。在

我基本上把输入图像和调色板转换成Lab颜色空间，然后计算每个像素到每个调色板条目的CIE76 Delta-E值的平方，并取最近的一个。在

#!/usr/bin/env python3

import numpy as np
from PIL import Image
from skimage import color

def CIE76DeltaE2(Lab1,Lab2):
    """Returns the square of the CIE76 Delta-E colour distance between 2 lab colours"""
    return (Lab2[0]-Lab1[0])*(Lab2[0]-Lab1[0]) + (Lab2[1]-Lab1[1])*(Lab2[1]-Lab1[1]) + (Lab2[2]-Lab1[2])*(Lab2[2]-Lab1[2])

def NearestPaletteIndex(Lab,palLab):
    """Return index of entry in palette that is nearest the given colour"""
    NearestIndex = 0
    NearestDist   = CIE76DeltaE2(Lab,palLab[0,0])
    for e in range(1,palLab.shape[0]):
        dist = CIE76DeltaE2(Lab,palLab[e,0])
        if dist < NearestDist:
            NearestDist = dist
            NearestIndex = e
    return NearestIndex

palette = [ 
    0,0,0,
    0,0,255,
    15,29,15,
    26,141,52,
    41,41,41,
    65,105,225,
    85,11,18,
    128,0,128,
    135,206,236,
    144,238,144,
    159,30,81,
    165,42,42,
    166,141,95,
    169,169,169,
    173,216,230,
    211,211,211,
    230,208,122,
    245,245,220,
    247,214,193,
    255,0,0,
    255,165,0,
    255,192,203,
    255,255,0,
    255,255,255
    ] + [0,] * 232 * 3


# Load the source image as numpy array and convert to Lab colorspace
imnp = np.array(Image.open('lion.png').convert('RGB'))
imLab = color.rgb2lab(imnp) 
h,w = imLab.shape[:2]

# Load palette as numpy array, truncate unused palette entries, and convert to Lab colourspace
palnp = np.array(palette,dtype=np.uint8).reshape(256,1,3)[:24,:]
palLab = color.rgb2lab(palnp)

# Make numpy array for output image
resnp = np.empty((h,w), dtype=np.uint8)

# Iterate over pixels, replacing each with the nearest palette entry
for y in range(0, h):
    for x in range(0, w):
        resnp[y, x] = NearestPaletteIndex(imLab[y,x], palLab)

# Create output image from indices, whack a palette in and save
resim = Image.fromarray(resnp, mode='P')
resim.putpalette(palette)
resim.save('result.png')

我明白了：

---

使用scipy.spatial.distance的cdist()函数似乎更快、更简洁：

^{pr2}$

Answer 2

如果速度是问题所在，ImageMagick可以更快地完成这项工作。它安装在大多数Linux发行版上，可用于macOS和Windows。在

基本上，您将创建一个24x1图像，名为"map.png"，调色板中每种颜色都有一个像素，然后告诉ImageMagick将狮子图像重新映射到Lab颜色空间中的颜色映射，而不会抖动。因此，终端/命令提示符中的命令是：

magick lion.png +dither -quantize Lab -remap map.png result.png

不到0.3秒。如果您想从Python中执行此操作，可以这样执行：

^{pr2}$