Laplace算子和Laplacian矩阵

1 Laplace算子的物理意义

Laplace算子的定义为梯度的散度。

在Cartesian坐标系下也可表示为:

或者,它是Hessian矩阵的迹:

以热传导方程为例,因为热流与温度的梯度成正比,那么温度的梯度的散度就是热量的损失率。

由此可见,Laplace算子可用于表现由于物质分布不均引起的物质输送。

2 Laplace算子的数学意义

现在,在一维空间中简单分析上面的式子:

也可以写作:

把分子第一项和第二项分别按泰勒展开:

可以看出Laplace算子实际上是一个使函数取平均的算子。多维空间相似。

3 Laplace方程

若Laplace算子右边为零,称为Laplace方程。Laplace方程的解称为调和函数。若右边是一个函数,称为泊松方程。

4 Laplace算子在图像处理的运用

图像处理是以像素作为基础离散化,如下:

5 Laplacian 矩阵

是一种用于表示图的矩阵。 它的维度是 |V|-by-|V| ( |V| 是节点的数目 )。 James Demmel提供了一种由Incidence matrix转化为Laplacian矩阵的方法。

In(G)是一个 |V|-by-|E| 矩阵( |E| 是边的数目 ), 设边e=(i,j),这一列除了第i行(为+1)和第j行(为-1)外都为零。 需要说明的是,根据这个定义,对于无向图 e=(i,j) 和 e=(j,i) 是等价的, 看似会生成很多不同的In图(根据每条边不同的取向)。但是实际上可以证明,无论边的方向怎么取,由In图生成的L图都是唯一的。 也就是说, e=(i,j) 和 e=(j,i) 怎么取是无关紧要的。 如何使用In图生成L图:

可得知Laplacian矩阵的两个重要性质:一是为对称阵。二是存在一个为零的特征值(秩为|V|-1)。三是一个半正定矩阵。 注意Laplace算子是负定的。

在求解含Laplacian矩阵的方程组时,常常要求为正定矩阵。观察发现这是因为Laplacian矩阵每列相加等于零。这时只需要手动更改第一行和第一列(比如第一个元素设为1,其余设为零),破坏其结构,令秩等于|V|就可以了。

对于非正定矩阵,左乘个transpose of the matrix, 推导如下:

Ax – b = 0

最小化 ||Ax – b||^2,展开后对x求导数:

可转化为正定方程组。

6 Laplace算子和Laplacian矩阵的关系

Laplace算子可以推广到多维情况计算。Laplacian矩阵主要用于三维以下的图形学计算,可以表现复杂的几何结构。而Lapace方程使用了Laplace算子来表示Laplacian矩阵。

时间: 2024-12-13 13:24:54

Laplace算子和Laplacian矩阵的相关文章

灰度图像--图像分割 Sobel算子,Prewitt算子和Scharr算子平滑能力比较

学习DIP第47天 转载请标明本文出处:http://blog.csdn.net/tonyshengtan ,出于尊重文章作者的劳动,转载请标明出处!文章代码已托管,欢迎共同开发: https://github.com/Tony-Tan/DIPpro 开篇废话 依然是废话,这篇主要想对比下Sobel,Prewitt和Scharr算子的平滑能力,由于一阶微分对噪声响应强,进行微分之前进行降噪是非常必要的,这里我们进行的实验是,以lena图作为实验原图,取其中一行数据作为无噪声的原始信号,分别加上不

高效的计算Laplacian 矩阵

直接上代码,可以直接调用: function L1 = CP_Laplacian( M ) %UNTITLED Summary of this function goes here % Detailed explanation goes here % number of vertices (size of matrix) nverts = size(M.vertices,1); % assume manifold and faces oriented uniformly edges = [M.f

基于边缘的分割

边缘检测是图像分割的另一种重要方法,利用图像灰度级在边缘处的突变,找到目标物体的边缘,图像中边缘处像素的灰度值不连续,这种不连续性可通过求导数来检测到.对于阶跃状边缘,其位置对应一阶导数的极值点,对应二阶导数的过零点(零交叉点).因此常用微分算子进行边缘检测.常用的一阶微分算子有Roberts算子.Prewitt算子和Sobel算子,二阶微分算子有Laplace算子和Kirsh算子等.在实际中各种微分算子常用小区域模板来表示,微分运算是利用模板和图像卷积来实现.这些算子对噪声敏感,只适合于噪声较

边缘检测与图像分割

1图像分割原理 图像分割的研究多年来一直受到人们的高度重视,至今提出了各种类型的分割算法.Pal把图像分割算法分成了6类:阈值分割,像素分割.深度图像分割.彩色图像分割,边缘检测和基于模糊集的方法.但是,该方法中,各个类别的内容是有重叠的.为了涵盖不断涌现的新方法,有的研究者将图像分割算法分为以下六类:并行边界分割技术.串行边界分割技术.并行区域分割技术.串行区域分割技术.结合特定理论工具的分割技术和特殊图像分割技术.而在较近的一篇综述中,更有学者将图像分割简单的分割数据驱动的分割和模型驱动的分

OpenCV-跟我一起学数字图像处理之拉普拉斯算子

Laplace算子和Sobel算子一样,属于空间锐化滤波操作.起本质与前面的Spatial Filter操作大同小异,下面就通过Laplace算子来介绍一下空间锐化滤波,并对OpenCV中提供的Laplacian函数进行一些说明. 数学原理 离散函数导数 离散函数的导数退化成了差分,一维一阶差分公式和二阶差分公式分别为, Laplace算子的差分形式 分别对Laplace算子x,y两个方向的二阶导数进行差分就得到了离散函数的Laplace算子. 在一个二维函数f(x,y)中,x,y两个方向的二阶

【数字图像处理】边缘检测与图像分割

原文链接:边缘检测与图像分割 作者:HUSTLX 1图像分割原理 图像分割的研究多年来一直受到人们的高度重视,至今提出了各种类型的分割算法.Pal把图像分割算法分成了6类:阈值分割,像素分割.深度图像分割.彩色图像分割,边缘检测和基于模糊集的方法.但是,该方法中,各个类别的内容是有重叠的.为了涵盖不断涌现的新方法,有的研究者将图像分割算法分为以下六类:并行边界分割技术.串行边界分割技术.并行区域分割技术.串行区域分割技术.结合特定理论工具的分割技术和特殊图像分割技术.而在较近的一篇综述中,更有学

【OpenCV入门教程之十二】OpenCV边缘检测:Canny算子,Sobel算子,Laplace算子,Scharr滤波器合辑

本系列文章由@浅墨_毛星云 出品,转载请注明出处. 文章链接:http://blog.csdn.net/poem_qianmo/article/details/25560901 作者:毛星云(浅墨)    微博:http://weibo.com/u/1723155442 知乎:http://www.zhihu.com/people/mao-xing-yun 邮箱: [email protected] 写作当前博文时配套使用的OpenCV版本: 2.4.9 本篇文章中,我们将一起学习OpenCV中

学习 opencv---(11)OpenC 边缘检测:Canny算子,Sobel算子,Laplace算子,Scharr滤波器

本篇文章中,我们将一起学习OpenCV中边缘检测的各种算子和滤波器--Canny算子,Sobel算子,Laplace算子以及Scharr滤波器.文章中包含了五个浅墨为大家准备的详细注释的博文配套源代码.在介绍四块知识点的时候分别一个,以及最后的综合示例中的一个.文章末尾提供配套源代码的下载. **** 给大家分享一个OpenCv中写代码是节约时间的小常识.其实OpenCv中,不用nameWindow,直接imshow就可以显示出窗口.大家看下文的示例代码就可以发现,浅墨在写代码的时候并没有用na

paper 82:边缘检测的各种微分算子比较(Sobel,Robert,Prewitt,Laplacian,Canny)

不同图像灰度不同,边界处一般会有明显的边缘,利用此特征可以分割图像.需要说明的是:边缘和物体间的边界并不等同,边缘指的是图像中像素的值有突变的地方,而物体间的边界指的是现实场景中的存在于物体之间的边界.有可能有边缘的地方并非边界,也有可能边界的地方并无边缘,因为现实世界中的物体是三维的,而图像只具有二维信息,从三维到二维的投影成像不可避免的会丢失一部分信息:另外,成像过程中的光照和噪声也是不可避免的重要因素.正是因为这些原因,基于边缘的图像分割仍然是当前图像研究中的世界级难题,目前研究者正在试图