- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我刚刚从这个博客学习了如何从头开始将高斯滤波器应用于 python 中的灰度图像:http://www.adeveloperdiary.com/data-science/computer-vision/applying-gaussian-smoothing-to-an-image-using-python-from-scratch/
现在我想将高斯滤波器应用于 3 channel (RGB)图像。
为此,我实现了代码,但我得到的输出是亮度非常低的模糊图像。此外,图像的边缘没有正确模糊。
这是我的代码:
# import libraries
import cv2
import numpy as np
import os
import matplotlib.pyplot as plt
%matplotlib inline
import cv2
# loading image
img_orig = cv2.imread('home.jpg')
# convert GBR image to RGB image
img_orig = cv2.cvtColor(img_orig, cv2.COLOR_BGR2RGB)
# Gaussian function
def dnorm(x, mu, sd):
return 1 / (np.sqrt(2 * np.pi) * sd) * np.exp(-((x-mu)/sd)** 2 / 2)
# function for making gaussian kernel
def gaussian_kernel(kernel_size, mu = 0):
# initializing mu and SD
sd = np.sqrt(kernel_size)
# creating 1D kernel
kernel_1D = np.linspace(-(kernel_size // 2), kernel_size // 2, kernel_size)
# normalizing 1D kernel
for i in range(kernel_size):
kernel_1D[i] = dnorm(kernel_1D[i], mu, sd)
# creating 2D kernel
kernel_2D = np.outer(kernel_1D, kernel_1D)
kernel_2D /= kernel_2D.max()
return kernel_2D
# Covolution function with zero padding
def convolution(image, kernel):
# find row and column of 3 channel (RGB) image
img_row, img_col, img_channel = image.shape
kernel_size = kernel.shape[0]
padding_width = (kernel_size - 1) // 2
#initialize output image
output = np.zeros(image.shape, dtype = np.uint8)
# initialize padded image with zeros
padded_img = np.zeros((img_row + 2*padding_width, img_col + 2*padding_width, img_channel), dtype = np.uint8)
# copy orignal image inside padded image
padded_img[padding_width : padding_width + img_row, padding_width : padding_width + img_col] = image
# average pixel values using gaussian kernel
for i in range(img_row):
for j in range(img_col):
# average each pixel's R channel value
output[i, j, 0] = np.sum(padded_img[i : i+kernel_size, j : j+kernel_size, 0] * kernel) // (kernel_size * kernel_size)
# average each pixel's G channel value
output[i, j, 1] = np.sum(padded_img[i : i+kernel_size, j : j+kernel_size, 1] * kernel) // (kernel_size * kernel_size)
# average each pixel's B channel value
output[i, j, 2] = np.sum(padded_img[i : i+kernel_size, j : j+kernel_size, 2] * kernel) // (kernel_size * kernel_size)
return output
def gaussian_filter(image, kernel_size = 3):
# initialize mu
mu = 0
# create gaussian kernel
kernel = gaussian_kernel(kernel_size, mu)
# apply convolution to image
conv_img = convolution(image, kernel)
# return blurred image
return conv_img
plt.figure(figsize = (7, 5))
print('orignal image')
plt.imshow(img_orig)
plt.show()
plt.figure(figsize = (7, 5))
print('blurred image')
plt.imshow(gaussian_filter(img_orig, 11))
plt.show()
print('openCV blurred image')
plt.imshow(cv2.GaussianBlur(img_orig, (11,11), 0))
plt.show()
最佳答案
有两件事是错误的,导致图像强度无法保留:首先通过除以最大值来归一化内核,然后在卷积中除以内核中的样本数。当您创建内核时,通过除以内核值的总和,而不是这两个规范化,只规范化一次。这使得内核权重的总和等于 1,并使卷积保持平均图像强度。请注意,卷积计算的是局部加权平均值;在加权平均中,我们需要将权重加到 1 以避免偏差。
暗边是由填充引起的:你用零(黑色)填充,它与卷积中图像边缘的值混合。 OpenCV 可能使用不同的边界条件或图像填充。典型的选项包括镜像值,或者只是将边缘值扩展。
最后,您的代码运行缓慢的主要原因是您在 Python 中使用了循环。 Python 是一种解释型语言,因此速度很慢。你可以使用 Numba 来加速循环(它是 Python 的即时编译器),或者简单地使用 NumPy 中的卷积,它是用编译语言实现的。
你的代码很慢的另一个原因(在你修复第一个之前这并不重要)是你没有利用高斯的可分离性。您可以通过将两个 1D 高斯相乘来构建 2D 高斯,但是您可以依次应用两个 1D 卷积。对于您的 11x11 内核示例,计算成本从 11*11=121 次乘法和加法减少到 11+11=22 次乘法和加法。内核越大,速度增益越好。
关于python - 为什么应用高斯滤镜后图像的亮度会降低?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/60987086/
尝试使用集成到 QTCreator 的表单编辑器,但即使我将插件放入 QtCreator.app/Contents/MacOS/designer 也不会显示。不过,相同的 dylib 文件确实适用于独
在此代码示例中。 “this.method2();”之后会读到什么?在返回returnedValue之前会跳转到method2()吗? public int method1(int returnedV
我的项目有通过gradle配置的依赖项。我想添加以下依赖项: compile group: 'org.restlet.jse', name: 'org.restlet.ext.apispark', v
我将把我们基于 Windows 的客户管理软件移植到基于 Web 的软件。我发现 polymer 可能是一种选择。 但是,对于我们的使用,我们找不到 polymer 组件具有表格 View 、下拉菜单
我的项目文件夹 Project 中有一个文件夹,比如 ED 文件夹,当我在 Eclipse 中指定在哪里查找我写入的文件时 File file = new File("ED/text.txt"); e
这是奇怪的事情,这个有效: $('#box').css({"backgroundPosition": "0px 250px"}); 但这不起作用,它只是不改变位置: $('#box').animate
这个问题在这里已经有了答案: Why does OR 0 round numbers in Javascript? (3 个答案) 关闭 5 年前。 Mozilla JavaScript Guide
这个问题在这里已经有了答案: Is the function strcmpi in the C standard libary of ISO? (3 个答案) 关闭 8 年前。 我有一个问题,为什么
我目前使用的是共享主机方案,我不确定它使用的是哪个版本的 MySQL,但它似乎不支持 DATETIMEOFFSET 类型。 是否存在支持 DATETIMEOFFSET 的 MySQL 版本?或者有计划
研究 Seam 3,我发现 Seam Solder 允许将 @Named 注释应用于包 - 在这种情况下,该包中的所有 bean 都将自动命名,就好像它们符合条件一样@Named 他们自己。我没有看到
我知道 .append 偶尔会增加数组的容量并形成数组的新副本,但 .removeLast 会逆转这种情况并减少容量通过复制到一个新的更小的数组来改变数组? 最佳答案 否(或者至少如果是,则它是一个错
很难说出这里要问什么。这个问题模棱两可、含糊不清、不完整、过于宽泛或夸夸其谈,无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开,visit the help center . 关闭 1
noexcept 函数说明符是否旨在 boost 性能,因为生成的对象中可能没有记录异常的代码,因此应尽可能将其添加到函数声明和定义中?我首先想到了可调用对象的包装器,其中 noexcept 可能会产
我正在使用 Angularjs 1.3.7,刚刚发现 Promise.all 在成功响应后不会更新 angularjs View ,而 $q.all 会。由于 Promises 包含在 native
我最近发现了这段JavaScript代码: Math.random() * 0x1000000 10.12345 10.12345 >> 0 10 > 10.12345 >>> 0 10 我使用
我正在编写一个玩具(物理)矢量库,并且遇到了 GHC 坚持认为函数应该具有 Integer 的问题。是他们的类型。我希望向量乘以向量以及标量(仅使用 * ),虽然这可以通过仅使用 Vector 来实现
PHP 的 mail() 函数发送邮件正常,但 Swiftmailer 的 Swift_MailTransport 不起作用! 这有效: mail('user@example.com', 'test
我尝试通过 php 脚本转储我的数据,但没有命令行。所以我用 this script 创建了我的 .sql 文件然后我尝试使用我的脚本: $link = mysql_connect($host, $u
使用 python 2.6.4 中的 sqlite3 标准库,以下查询在 sqlite3 命令行上运行良好: select segmentid, node_t, start, number,title
我最近发现了这段JavaScript代码: Math.random() * 0x1000000 10.12345 10.12345 >> 0 10 > 10.12345 >>> 0 10 我使用
我是一名优秀的程序员,十分优秀!