Как правильно прочитать изображение в режиме YCbCr?

Question

Как узнать, правильно ли я прочитал изображение PNG в режиме YCbCr? Я получаю разные значения пикселей, что сбивает с толку.

def convert_rgb_to_ycbcr(img):
    y = 16. + (64.738 * img[:, :, 0] + 129.057 * img[:, :, 1] + 25.064 * img[:, :, 2]) / 255.
    cb = 128. + (-37.945 * img[:, :, 0] - 74.494 * img[:, :, 1] + 112.439 * img[:, :, 2]) / 255.
    cr = 128. + (112.439 * img[:, :, 0] - 94.154 * img[:, :, 1] - 18.285 * img[:, :, 2]) / 255.
    return np.array([y, cb, cr]).transpose([1, 2, 0])


# method 1 - read as YCbCr directly
img = scipy.misc.imread(path, mode='YCbCr').astype(np.float)
print(img[0, :5, 0]) 
# returns [32. 45. 68. 78. 92.]

# method 2 - read as RGB and convert RGB to YCbCr
img = scipy.misc.imread(path, mode='RGB').astype(np.float)
img = convert_rgb_to_ycbcr(img)
print(img[0, :5, 0]) 
# returns[44.0082902  55.04281961 75.1105098  83.57022745 95.44837255]

Я хочу использовать метод 1, поскольку scipy уже позаботился о преобразовании за меня, но мне не удалось найти для него исходный код. Поэтому я сам определил функцию преобразования, но получаю разные значения пикселей.

Answer 1

В последних версиях scipyimreadустарела. Однако он используетImage.convert от PIL преобразовать режимы.

Детали:

https://pillow.readthedocs.io/en/3.1.x/reference/Image.html?highlight=convert

https://pillow.readthedocs.io/en/3.1.x/handbook/concepts.html

https://github.com/scipy/scipy/blob/v0.18.0/scipy/misc/pilutil.py

Я изменил твою convert_rgb_to_ycbcr(img) функция и дает тот же результат.

Используемая реализация: https://docs.microsoft.com/en-us/openspecs/windows_protocols/ms-rdprfx/b550d1b5-f7d9-4a0c-9141-b3dca9d7f525?redirectedfrom=MSDN

Формула преобразования из RGB в YCbCr

import scipy.misc # scipy 1.1.0
import numpy as np

def convert_rgb_to_ycbcr(im):
    xform = np.array([[.299, .587, .114], [-.1687, -.3313, .5], [.5, -.4187, -.0813]])
    ycbcr = im.dot(xform.T)
    ycbcr[:,:,[1,2]] += 128
    return np.uint8(ycbcr)


# method 1 - read as YCbCr directly
img = scipy.misc.imread('test.jpg', mode='YCbCr').astype(np.float)
print(img[0, :5, 0]) 
# returns [32. 45. 68. 78. 92.]

# method 2 - read as RGB and convert RGB to YCbCr
img = scipy.misc.imread('test.jpg', mode='RGB').astype(np.float)
img = convert_rgb_to_ycbcr(img)
print(img[0, :5, 0]) 

[165. 165. 165. 166. 167.]
[165 165 165 166 167]