如何在 Python 中将 RGB 图像转换为灰度?

2024-12-11 08:47:00
admin
原创
184
摘要:问题描述:我正在尝试matplotlib读取 RGB 图像并将其转换为灰度。在 matlab 中我使用这个:img = rgb2gray(imread('image.png')); 在matplotlib 教程中,他们没有介绍它。他们只是在图像中读取import matplotlib.image as mpi...

问题描述:

我正在尝试matplotlib读取 RGB 图像并将其转换为灰度。

在 matlab 中我使用这个:

img = rgb2gray(imread('image.png'));

在matplotlib 教程中,他们没有介绍它。他们只是在图像中读取

import matplotlib.image as mpimg
img = mpimg.imread('image.png')

然后他们对数组进行切片,但据我所知,这与将 RGB 转换为灰度不同。

lum_img = img[:,:,0]

我很难相信 numpy 或 matplotlib 没有内置函数来将 rgb 转换为灰色。这不是图像处理中的常见操作吗?

我编写了一个非常简单的函数,可以在 5 分钟内处理导入的图像imread。它的效率非常低,但这就是我希望内置专业实现的原因。

Sebastian 改进了我的功能,但我仍然希望找到内置的功能。

matlab(NTSC/PAL)的实现:

import numpy as np

def rgb2gray(rgb):

    r, g, b = rgb[:,:,0], rgb[:,:,1], rgb[:,:,2]
    gray = 0.2989 * r + 0.5870 * g + 0.1140 * b

    return gray

解决方案 1:

使用Pillow来做怎么样:

from PIL import Image
img = Image.open('image.png').convert('L')
img.save('greyscale.png')

如果输入图像中存在 alpha(透明度)通道且应保留,则使用模式LA

img = Image.open('image.png').convert('LA')

使用 matplotlib 和公式

Y' = 0.2989 R + 0.5870 G + 0.1140 B 

你可以这样做:

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.image as mpimg

def rgb2gray(rgb):
    return np.dot(rgb[...,:3], [0.2989, 0.5870, 0.1140])

img = mpimg.imread('image.png')     
gray = rgb2gray(img)    
plt.imshow(gray, cmap=plt.get_cmap('gray'), vmin=0, vmax=1)
plt.show()

解决方案 2:

您还可以使用scikit-image,它提供了一些将图像转换为 的函数ndarray,例如rgb2gray

from skimage import color
from skimage import io

img = color.rgb2gray(io.imread('image.png'))

:此转换中使用的权重已针对当代 CRT 荧光粉进行了校准:Y = 0.2125 R + 0.7154 G + 0.0721 B

或者,您可以通过以下方式读取灰度图像:

from skimage import io
img = io.imread('image.png', as_gray=True)

解决方案 3:

使用在 Ubuntu 16.04 LTS(带有 SSD 的 Xeon E5 2670)上的 Python 3.5 运行的 1000 个 RGBA PNG 图像(224 x 256 像素)对三种建议的方法的速度进行了测试。

平均运行时间

pil :1.037 秒

scipy:1.040 秒

sk :2.120 秒

PIL 和 SciPy 给出了相同的numpy数组(范围从 0 到 255)。SkImage 给出的数组范围从 0 到 1。此外,颜色转换略有不同,请参见 CUB-200 数据集中的示例。

SkImage: 图像处理软件

PIL : 太平船务

SciPy : 科学

Original: 原来的

Diff : 在此处输入图片描述

代码

  1. 表现

run_times = dict(sk=list(), pil=list(), scipy=list())
for t in range(100):
    start_time = time.time()
    for i in range(1000):
        z = random.choice(filenames_png)
        img = skimage.color.rgb2gray(skimage.io.imread(z))
    run_times['sk'].append(time.time() - start_time)

start_time = time.time()
for i in range(1000):

z = random.choice(filenames_png)
img = np.array(Image.open(z).convert('L'))

run_times['pil'].append(time.time() - start_time)

start_time = time.time()
for i in range(1000):

z = random.choice(filenames_png)
img = scipy.ndimage.imread(z, mode='L')

run_times['scipy'].append(time.time() - start_time)



for k, v in run_times.items():
    print('{:5}: {:0.3f} seconds'.format(k, sum(v) / len(v)))


  1. 输出

z = 'Cardinal_0007_3025810472.jpg'
img1 = skimage.color.rgb2gray(skimage.io.imread(z)) * 255
IPython.display.display(PIL.Image.fromarray(img1).convert('RGB'))
img2 = np.array(Image.open(z).convert('L'))
IPython.display.display(PIL.Image.fromarray(img2))
img3 = scipy.ndimage.imread(z, mode='L')
IPython.display.display(PIL.Image.fromarray(img3))
  1. 比较

img_diff = np.ndarray(shape=img1.shape, dtype='float32')
img_diff.fill(128)
img_diff += (img1 - img3)
img_diff -= img_diff.min()
img_diff *= (255/img_diff.max())
IPython.display.display(PIL.Image.fromarray(img_diff).convert('RGB'))
  1. 导入

import skimage.color
import skimage.io
import random
import time
from PIL import Image
import numpy as np
import scipy.ndimage
import IPython.display
  1. 版本

skimage.version
0.13.0
scipy.version
0.19.1
np.version
1.13.1

解决方案 4:

imread您可以从一开始就使用OpenCV将图像文件读取为灰度:

img = cv2.imread('messi5.jpg', 0)

此外,如果您想将图像读取为 RGB,请进行一些处理,然后转换为灰度,您可以使用cvtcolorOpenCV:

gray_image = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)

解决方案 5:

目前最快且最有效的方法是使用Pillow,通过 安装pip install Pillow

代码如下:

from PIL import Image
img = Image.open('input_file.jpg').convert('L')
img.save('output_file.jpg')

解决方案 6:

使用此公式

Y' = 0.299 R + 0.587 G + 0.114 B 

我们可以做到

import imageio
import numpy as np
import matplotlib.pyplot as plt

pic = imageio.imread('(image)')
gray = lambda rgb : np.dot(rgb[... , :3] , [0.299 , 0.587, 0.114]) 
gray = gray(pic)  
plt.imshow(gray, cmap = plt.get_cmap(name = 'gray'))

但是,GIMP将彩色图像转换为灰度图像的软件有三种算法来完成这项任务。

解决方案 7:

本教程有点作弊,因为它从用 RGB 编码的灰度图像开始,所以他们只是切分单个颜色通道并将其视为灰度。您需要执行的基本步骤是从 RGB 颜色空间转换为使用近似亮度/色度模型进行编码的颜色空间,例如 YUV/YIQ 或 HSL/HSV,然后切分类似亮度的通道并将其用作灰度图像。 matplotlib似乎没有提供转换为 YUV/YIQ 的机制,但它确实允许您转换为 HSV。

尝试使用matplotlib.colors.rgb_to_hsv(img)然后从数组中切出最后一个值 (V) 作为灰度。它与亮度值不太一样,但这意味着您可以在 中完成所有操作matplotlib

背景:

或者,您可以使用 PIL 或内置函数colorsys.rgb_to_yiq()转换为具有真实亮度值的色彩空间。您也可以全力以赴,推出自己的仅限亮度的转换器,尽管这可能有点过头了。

解决方案 8:

你可以这样做:

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.image as mpimg

def rgb_to_gray(img):
        grayImage = np.zeros(img.shape)
        R = np.array(img[:, :, 0])
        G = np.array(img[:, :, 1])
        B = np.array(img[:, :, 2])

        R = (R *.299)
        G = (G *.587)
        B = (B *.114)

        Avg = (R+G+B)
        grayImage = img.copy()

        for i in range(3):
           grayImage[:,:,i] = Avg
           
        return grayImage       

image = mpimg.imread("your_image.png")   
grayImage = rgb_to_gray(image)  
plt.imshow(grayImage)
plt.show()

解决方案 9:

如果你已经在使用 NumPy/SciPy,那么你也可以使用:

scipy.ndimage.imread(file_name, mode='L')

解决方案 10:

使用 OpenCV 很简单:

import cv2

im = cv2.imread("flower.jpg")

# To Grayscale
im = cv2.cvtColor(im, cv2.COLOR_BGR2GRAY)
cv2.imwrite("grayscale.jpg", im)

# To Black & White
im = cv2.threshold(im, 127, 255, cv2.THRESH_BINARY)[1]
cv2.imwrite("black-white.jpg", im)

在此处输入图片描述

解决方案 11:

使用img.Convert(),支持“L”、“RGB”和“CMYK”模式。

import numpy as np
from PIL import Image

img = Image.open("IMG/center_2018_02_03_00_34_32_784.jpg")
img.convert('L')

print np.array(img)

输出:

[[135 123 134 ...,  30   3  14]
 [137 130 137 ...,   9  20  13]
 [170 177 183 ...,  14  10 250]
 ..., 
 [112  99  91 ...,  90  88  80]
 [ 95 103 111 ..., 102  85 103]
 [112  96  86 ..., 182 148 114]]

解决方案 12:

我通过谷歌找到了这个问题,寻找一种将已加载的图像转换为灰度的方法。

以下是使用 SciPy 实现的方法:

import scipy.misc
import scipy.ndimage

# Load an example image
# Use scipy.ndimage.imread(file_name, mode='L') if you have your own
img = scipy.misc.face()

# Convert the image
R = img[:, :, 0]
G = img[:, :, 1]
B = img[:, :, 2]
img_gray = R * 299. / 1000 + G * 587. / 1000 + B * 114. / 1000

# Show the image
scipy.misc.imshow(img_gray)

解决方案 13:

当像素中所有 3 个颜色通道(RGB)的值相同时,该像素将始终为灰度格式。

将 RGB 图像转换为灰度的一个简单直观的方法是取每个像素中所有颜色通道的平均值,然后将该值分配回该像素。

import numpy as np
from PIL import Image

img=np.array(Image.open('sample.jpg')) #Input - Color image
gray_img=img.copy()

for clr in range(img.shape[2]):
    gray_img[:,:,clr]=img.mean(axis=2) #Take mean of all 3 color channels of each pixel and assign it back to that pixel(in copied image)

#plt.imshow(gray_img) #Result - Grayscale image

输入图像:
输入图像

输出图像:
输出图像

解决方案 14:

假设我的图像原始形式为 3 通道

my_image = cv2.imread("./5d10e5939c5101174c54bb98.png")
#greyscaling the image
image_sum = my_image.sum(axis=2)
new_image = image_sum/image_sum.max() 

new_image这是我的单通道灰度图像

plt.figure(figsize=[12,8])
plt.imshow(new_image, cmap=plt.cm.gray)
plt.show()

[这没有使用 cv2 的cv2.COLOR_BGR2GRAY参数或 PIL 的.convert('L')方法]

解决方案 15:

image=myCamera.getImage().crop(xx,xx,xx,xx).scale(xx,xx).greyscale()

可以greyscale()直接用来进行转换。

相关推荐
  政府信创国产化的10大政策解读一、信创国产化的背景与意义信创国产化,即信息技术应用创新国产化,是当前中国信息技术领域的一个重要发展方向。其核心在于通过自主研发和创新,实现信息技术应用的自主可控,减少对外部技术的依赖,并规避潜在的技术制裁和风险。随着全球信息技术竞争的加剧,以及某些国家对中国在科技领域的打压,信创国产化显...
工程项目管理   2500  
  为什么项目管理通常仍然耗时且低效?您是否还在反复更新电子表格、淹没在便利贴中并参加每周更新会议?这确实是耗费时间和精力。借助软件工具的帮助,您可以一目了然地全面了解您的项目。如今,国内外有足够多优秀的项目管理软件可以帮助您掌控每个项目。什么是项目管理软件?项目管理软件是广泛行业用于项目规划、资源分配和调度的软件。它使项...
项目管理软件   1541  
  PLM(产品生命周期管理)项目对于企业优化产品研发流程、提升产品质量以及增强市场竞争力具有至关重要的意义。然而,在项目推进过程中,范围蔓延是一个常见且棘手的问题,它可能导致项目进度延迟、成本超支以及质量下降等一系列不良后果。因此,有效避免PLM项目范围蔓延成为项目成功的关键因素之一。以下将详细阐述三大管控策略,助力企业...
plm系统   16  
  PLM(产品生命周期管理)项目管理在企业产品研发与管理过程中扮演着至关重要的角色。随着市场竞争的加剧和产品复杂度的提升,PLM项目面临着诸多风险。准确量化风险优先级并采取有效措施应对,是确保项目成功的关键。五维评估矩阵作为一种有效的风险评估工具,能帮助项目管理者全面、系统地评估风险,为决策提供有力支持。五维评估矩阵概述...
免费plm软件   23  
  引言PLM(产品生命周期管理)开发流程对于企业产品的全生命周期管控至关重要。它涵盖了从产品概念设计到退役的各个阶段,直接影响着产品质量、开发周期以及企业的市场竞争力。在当今快速发展的科技环境下,客户对产品质量的要求日益提高,市场竞争也愈发激烈,这就使得优化PLM开发流程成为企业的必然选择。缺陷管理工具和六西格玛方法作为...
plm产品全生命周期管理   26  
热门文章
项目管理软件有哪些?
曾咪二维码

扫码咨询,免费领取项目管理大礼包!

云禅道AD
禅道项目管理软件

云端的项目管理软件

尊享禅道项目软件收费版功能

无需维护,随时随地协同办公

内置subversion和git源码管理

每天备份,随时转为私有部署

免费试用