高斯函数在图像处理中的应用资源-CSDN文库

高斯函数

4星 · 超过85%的资源需积分: 50 60 浏览量 2010-06-18 14:55:48 上传评论 5 收藏 390KB PDF 举报

资源详情

资源评论

3. The Gaussian kernel

"Everybody believes in the exponential law of errors: the experimenters, because

they think it can be proved by mathematics; and the mathematicians, because they

believe it has been established by observation" (Lippman in [Whittaker1967, p. 179]).

3.1 The Gaussian kernel

The Gaussian (better Gaußian) kernel is named after Carl Friedrich Gauß (1777-1855), a brilliant German

mathematician. This chapter discusses many of the nice and peculiar properties of the Gaussian kernel.

FEVinit‘;

FEVFunctions‘;

Show

Import

"Gauss10DM.gif"

;

Figure 3.1. The Gaussian kernel is apparent on every German banknote of DM 10,- where it is depicted

next to its famous inventor when he was 55 years old. The new Euro replaces these banknotes.

The Gaussian kernel is defined in 1-D, 2D and N-D respectively as

þþþþþþþþþþþþþþþþþ

! !!! !!!!

p s

þþþþþþþþþþþþþ

, G

x, y;

þþþþþþþþþþþþ

þþþþþþþþþþþþþþþþþþ

, G

;

þþþþþþþþþþþþþþþþþþþþþþþ

!! !!! !!!

p s



¸»



þþþþþþþþþþþþþ

The

determines the width of the Gaussian kernel. In statistics, when we consider the Gaussian probability

density function it is called the standard deviation, and the square of it,

, the variance. In the rest of this

book, when we consider the Gaussian as an aperture function of some observation, we will refer to

as the

inner scale or shortly scale. In the whole of this book the scale can only take positive values,

0. In the

process of observation

can never become zero. For, this would imply making an observation through an

infinitesimally small aperture, which is impossible. The factor of 2 in the exponent is a matter of convention,

because we then have a 'cleaner' formula for the diffusion equation, as we will see later on. The semicolon

between the spatial and scale parameters is conventionally put there to make the difference between these

parameters explicit. The scale-dimension is not just another spatial dimension, as we will thoroughly discuss

in the remainder of this book.

来源网址是：http://www.stat.wisc.edu/%7Emchung/teaching/MIA/reading/

3.2 Normalization

The term

þþþþþþþþþþþþþþþþþ

!!!!!! !!

p s

in front of the one-dimensional Gaussian kernel is the normalization constant. It comes from

the fact that the integral over the exponential function is not unity:

-





!!!!! !!!

p s

. With the

normalization constant this Gaussian kernel is a normalized kernel, i.e. its integral over its full domain is unity

for every

. This means that increasing the

of the kernel reduces the amplitude substantially. Let us look at

the graphs of the normalized kernels for

0.3,

1 and

2 plotted on the same axes:

Unprotect

gauss

; gauss

x_,

þþþþþþþþþþþþþþþþþ

!!!!!!!

Exp

þþþþþþþþþþ

;

Block

$DisplayFunction

Identity

p1, p2, p3

Plot

gauss

5, 5

, PlotRange

0, 1.4



.3, 1, 2

;

Show

GraphicsArray

p1, p2, p3

, ImageSize

450, 130

;

2 4

0.2

0.4

0.6

0.8

1.2

2 4

0.2

0.4

0.6

0.8

1.2

2 4

0.2

0.4

0.6

0.8

1.2

Figure 3.2. The Gaussian function at scales

.3,

1 and

2. The kernel is normalized, so the area

under the curve is always unity.

The normalization ensures that the average greylevel of the image remains the same when we blur the image

with this kernel. This is known as average grey level invariance.

3.3 Cascade property

The shape of the kernel remains the same, irrespective of the

. When we convolve two Gaussian kernels we

get a new wider Gaussian with a variance

which is the sum of the variances of the constituting Gaussians:

new

;

.; FullSimplify

-



gauss

0, Im

þþþþþþþþþþþþþþþþþþþþþþþþþ

þþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþþ

!!!!!!!!

"################

This phenomenon, i.e. that a new function emerges that is similar to the constituting functions, is called self-

similarity. The Gaussian is a self-similar function. Convolution with a Gaussian is a linear operation, so a

convolution with a Gaussian kernel followed by a convolution with again a Gaussian kernel is equivalent to

convolution with the broader kernel. Note that the squares of

add, not the

's themselves. Of course we can

concatenate as many blurring steps as we want to create a larger blurring step. With analogy to a cascade of

waterfalls spanning the same height as the total waterfall, this phenomenon is also known as the cascade

smoothing property.

2 03Gaussiankernel.nb

3.4 The scale parameter

In order to avoid the summing of squares, one often uses the following parametrization: 2



t, so the

Gaussian kernel get a particular short form. In N dimensions:G

, t

þþþþþþþþþþþþþþþþ



þþþþþþþþ

. It is this t that emerges

in the diffusion equation



þþþþþþþ



þþþþþþþþþ



þþþþþþþþþ



þþþþþþþþþ



. It is often referred to as ’scale’ (like in: differentiation to

scale,



þþþþþþþ



), but a better name is variance.

To make the self-similarity of the Gaussian kernel explicit, we can introduce a new dimensionless spatial

parameter, x

þþþþþþþþþþþþþ

!!!!

. We say that we have reparametrized the x-axis. Now the Gaussian kernel becomes:

;

þþþþþþþþþþþþþþþþþ

!!!!!! !!

, or g

; t

þþþþþþþþþþþþþþþþ



. In other words: if we walk along the spatial axis in

footsteps expressed in scale-units (

’s), all kernels are of equal size or ’width’ (but due to the normalization

constraint not necessarily of the same amplitude). We now have a ’natural’ size of footstep to walk over the

spatial coordinate: a unit step in x is now

!!!!!

2 , so in more blurred images we make bigger steps. We call

this basic Gaussian kernel the natural Gaussian kernel g

;

. The new coordinate x

þþþþþþþþþþþþþ

!!!!

is called the

natural coordinate. It eliminates the scale factor

from the spatial coordinates, i.e. it makes the Gaussian

kernels similar, despite their different inner scales. We will encounter natural coordinates many times hereafter.

The spatial extent of the Gaussian kernel ranges from -



to +



, but in practice it has negligeable values for x

larger then a few (say 5)

. The numerical value at x=5

, and the area under the curve from x=5

to infinity

(recall that the total area is 1):

gauss

5, 1

D

Integrate

gauss

x, 1

x, 5, Infinity

<D 

1.48672



2.86652



The larger we make the standard deviation

, the more the image gets blurred. In the limit to infinity, the

image becomes homogenous in intensity. The final intensity is the average intensity of the image. This is true

for an image with infinite extent, which in practice will never occur, of course. The boundary has to be taken

into account. Actually, one can take many choices what to do at the boundary, it is a matter of concensus.

Boundaries are discussed in detail in chapter 5, where practical issues of computer implementation are

discussed.

3.5 Relation to generalized functions

The Gaussian kernel is the physical equivalent of the mathematical point. It is not strictly local, like the

mathematical point, but semi-local. It has a Gaussian weighted extent, indicated by its inner scale

. Because

scale-space theory is revolving around the Gaussian function and its derivatives as a physical differential

operator (in more detail explained in the next chapter), we will focus here on some mathematical notions that

are directly related, i.e. the mathematical notions underlying sampling of values from functions and their

derivatives at selected points (i.e. that is why it is referred to as sampling). The mathematical functions

involved are the generalized functions, i.e. the Delta-Dirac function, the Heavyside function and the error

function. We study in the next section these functions in more detail.

剩余12页未读，继续阅读

评论收藏

内容反馈

zyq211985

2013-03-20

很不错嘛，楼主的渊博的知识令小弟折服。

高斯函数在图像处理中的应用

评论12

最新资源

高斯函数在图像处理中的应用

评论12

最新资源

相关推荐

高斯滤波在图像滤波中的应用

用于数字图像处理的Matlab高斯滤波

Matlab 生成高斯函数图像

基于高斯函数假设的图像频谱恢复特性分析方法.pdf

Gauss_CMAC-master.rar_base3bl_cmac_图像处理_高斯函数

统计信号处理：非高斯信号处理及其应用

高斯-瑞利混合模型在遥感图像分割中的应用

数字图像处理2_lena_matlab_energyzwg_中值滤波_图像处理_

图像处理程序集

图像处理的matlab程序

MATLAB的图像处理基础的资料

小波分析及其在数字图像处理中的应用 [朱希安，曹林 编著] 2012年版

matlab辅助图像处理源代码

Matlab图像处理程序集

基于核回归 Kernel regression 的图像处理作者Takeda的matlab代码 应用于去噪 超分 去模糊等等

2012 Canny 边缘检测算子研究及其在医学图像中的应用.pdf

matlab-图像处理算法

老师布置的数字图像处理实验 c#（主要参考c#数字图像算法典型实例）

高斯伪谱法matlab代码

图像处理（第二版）章毓晋

python大作业 含爬虫、数据可视化、地图、报告、及源码（整和为一个文件）（2014-2020全国各地区原油加工量）.rar

仿真电路以及操作方法

【纯干货啊】华为IPD流程管理(完整版).pptx

可编程语言标准IEC61131-3中文版.pdf

OFDM完整仿真过程与教程.zip

信号与系统——保研复习资料.pdf

Landsat_WRS2.zip

小波分析及其在数字图像处理中的应用 [朱希安，曹林编著] 2012年版

基于核回归 Kernel regression 的图像处理作者Takeda的matlab代码应用于去噪超分去模糊等等

python大作业含爬虫、数据可视化、地图、报告、及源码（整和为一个文件）（2014-2020全国各地区原油加工量）.rar