二维高斯函数的两种形式

第一种形式很常见

多元正态分布

多元正态分布(Multivariate Normal Distribution),也称为多变量正态分布或多维正态分布,是统计学中一种重要的概率分布,用于描述多个随机变量的联合分布。

假设有 n n n 个随机变量 X 1 , X 2 , … , X n X_1, X_2, \ldots, X_n X1,X2,,Xn,它们服从多元正态分布,则这 n n n 个随机变量的联合概率密度函数(PDF)可以表示为:

f ( x 1 , x 2 , … , x n ) = 1 ( 2 π ) n ∣ Σ ∣ exp ⁡ ( − 1 2 ( x − μ ) ⊤ Σ − 1 ( x − μ ) ) f(x_1, x_2, \ldots, x_n) = \frac{1}{\sqrt{(2\pi)^n |\Sigma|}} \exp\left(-\frac{1}{2} (x - \mu)^\top \Sigma^{-1} (x - \mu)\right) f(x1,x2,,xn)=(2π)n∣Σ∣ 1exp(21(xμ)Σ1(xμ))

其中:

  • x = ( x 1 , x 2 , … , x n ) ⊤ x = (x_1, x_2, \ldots, x_n)^\top x=(x1,x2,,xn) 是一个 n × 1 n \times 1 n×1 的列向量,表示 n n n 个随机变量的观测值。
  • μ = ( μ 1 , μ 2 , … , μ n ) ⊤ \mu = (\mu_1, \mu_2, \ldots, \mu_n)^\top μ=(μ1,μ2,,μn) 是一个 n × 1 n \times 1 n×1 的列向量,表示 n n n 个随机变量的均值。
  • Σ \Sigma Σ 是一个 n × n n \times n n×n 的对称正定矩阵,称为协方差矩阵,其元素 σ i j \sigma_{ij} σij 表示 X i X_i Xi X j X_j Xj 之间的协方差。
  • ∣ Σ ∣ |\Sigma| ∣Σ∣ 表示协方差矩阵 Σ \Sigma Σ 的行列式。
  • Σ − 1 \Sigma^{-1} Σ1 是协方差矩阵 Σ \Sigma Σ 的逆矩阵。
  • ( x − μ ) ⊤ (x - \mu)^\top (xμ) ( x − μ ) (x - \mu) (xμ) 的转置。

二维高斯函数(也称为二维正态分布函数)是统计学和物理学中常用的一个函数,用于描述二维空间中的随机变量的分布。这里我们讨论两种形式的二维高斯函数,其中第二种形式通过几何变换(如旋转和平移)从第一种形式推导出来。

二元正态分布

在二元正态分布中,如果已知两个变量 X X X Y Y Y的相关系数 ρ \rho ρ以及它们的标准差 σ X \sigma_X σX σ Y \sigma_Y σY,那么协方差矩阵 Σ \Sigma Σ可以表示为:

Σ = ( σ X 2 ρ σ X σ Y ρ σ X σ Y σ Y 2 ) \Sigma = \begin{pmatrix} \sigma_X^2 & \rho \sigma_X \sigma_Y \\ \rho \sigma_X \sigma_Y & \sigma_Y^2 \end{pmatrix} Σ=(σX2ρσXσYρσXσYσY2)

这里, σ X 2 \sigma_X^2 σX2 σ Y 2 \sigma_Y^2 σY2分别是 X X X Y Y Y的方差,而 ρ σ X σ Y \rho \sigma_X \sigma_Y ρσXσY是它们的协方差。

接下来,需要计算这个协方差矩阵的逆矩阵 Σ − 1 \Sigma^{-1} Σ1。首先,计算协方差矩阵的行列式 ∣ Σ ∣ |\Sigma| ∣Σ∣

∣ Σ ∣ = σ X 2 σ Y 2 − ( ρ σ X σ Y ) 2 = σ X 2 σ Y 2 ( 1 − ρ 2 ) |\Sigma| = \sigma_X^2 \sigma_Y^2 - (\rho \sigma_X \sigma_Y)^2 = \sigma_X^2 \sigma_Y^2 (1 - \rho^2) ∣Σ∣=σX2σY2(ρσXσY)2=σX2σY2(1ρ2)

然后,利用行列式的值和协方差矩阵的元素来构造逆矩阵。逆矩阵 Σ − 1 \Sigma^{-1} Σ1的元素可以表示为:

Σ − 1 = 1 ∣ Σ ∣ ( σ Y 2 − ρ σ X σ Y − ρ σ X σ Y σ X 2 ) \Sigma^{-1} = \frac{1}{|\Sigma|} \begin{pmatrix} \sigma_Y^2 & -\rho \sigma_X \sigma_Y \\ -\rho \sigma_X \sigma_Y & \sigma_X^2 \end{pmatrix} Σ1=∣Σ∣1(σY2ρσXσYρσXσYσX2)

将行列式的值代入上式,得到:

Σ − 1 = 1 σ X 2 σ Y 2 ( 1 − ρ 2 ) ( σ Y 2 − ρ σ X σ Y − ρ σ X σ Y σ X 2 ) \Sigma^{-1} = \frac{1}{\sigma_X^2 \sigma_Y^2 (1 - \rho^2)} \begin{pmatrix} \sigma_Y^2 & -\rho \sigma_X \sigma_Y \\ -\rho \sigma_X \sigma_Y & \sigma_X^2 \end{pmatrix} Σ1=σX2σY2(1ρ2)1(σY2ρσXσYρσXσYσX2)

进一步化简,得到:

Σ − 1 = 1 1 − ρ 2 ( 1 σ X 2 − ρ σ X σ Y − ρ σ X σ Y 1 σ Y 2 ) \Sigma^{-1} = \frac{1}{1 - \rho^2} \begin{pmatrix} \frac{1}{\sigma_X^2} & -\frac{\rho}{\sigma_X \sigma_Y} \\ -\frac{\rho}{\sigma_X \sigma_Y} & \frac{1}{\sigma_Y^2} \end{pmatrix} Σ1=1ρ21(σX21σXσYρσXσYρσY21)

这个逆矩阵代入多元正态分布就是第一种形式。

第一种形式:标准二维高斯函数

标准二维高斯函数(也称为二维正态分布的概率密度函数)的表达式为:

f ( x , y ) = 1 2 π σ x σ y 1 − ρ 2 exp ⁡ ( − 1 2 ( 1 − ρ 2 ) [ ( x − μ x σ x ) 2 − 2 ρ ( x − μ x σ x ) ( y − μ y σ y ) + ( y − μ y σ y ) 2 ] ) f(x, y) = \frac{1}{2\pi\sigma_x\sigma_y\sqrt{1-\rho^2}} \exp\left(-\frac{1}{2(1-\rho^2)}\left[\left(\frac{x-\mu_x}{\sigma_x}\right)^2 - 2\rho\left(\frac{x-\mu_x}{\sigma_x}\right)\left(\frac{y-\mu_y}{\sigma_y}\right) + \left(\frac{y-\mu_y}{\sigma_y}\right)^2\right]\right) f(x,y)=2πσxσy1ρ2 1exp(2(1ρ2)1[(σxxμx)22ρ(σxxμx)(σyyμy)+(σyyμy)2])

其中,

  • μ x \mu_x μx μ y \mu_y μy 分别是 x x x y y y 的均值(位置参数)。
  • σ x \sigma_x σx σ y \sigma_y σy 分别是 x x x y y y 的标准差(尺度参数)。
  • ρ \rho ρ x x x y y y 的相关系数(形状参数)。

第二种出现在Gabor滤波器中,第二种又有多种变形,认清它是高斯函数。

Gabor函数

Gabor滤波器是一种基于Gabor函数的特定频率和方向选择性滤波器,是一种用于图像纹理分析和特征提取的线性滤波器,由Dennis Gabor于1946年提出,广泛应用于图像处理和计算机视觉领域。在空间域中,一个二维Gabor滤波器可以看作是一个正弦平面波和高斯核函数的乘积。

g ( x , y ; λ , θ , ψ , σ , γ ) = exp ⁡ ( − x ′ 2 + γ 2 y ′ 2 2 σ 2 ) cos ⁡ ( 2 π x ′ λ + ψ ) g(x, y; \lambda, \theta, \psi, \sigma, \gamma) = \exp\left(-\frac{x'^2 + \gamma^2 y'^2}{2\sigma^2}\right) \cos\left(2\pi\frac{x'}{\lambda} + \psi\right) g(x,y;λ,θ,ψ,σ,γ)=exp(2σ2x′2+γ2y′2)cos(2πλx+ψ)

其中,

  • ( x , y ) (x, y) (x,y) 是图像中像素的坐标。
  • λ \lambda λ 是正弦函数的波长。
  • θ \theta θ 是Gabor核函数的方向,表示Gabor函数平行条纹的方向。
  • ψ \psi ψ 是相位偏移,通常设置为0或 π 2 \frac{\pi}{2} 2π
  • σ \sigma σ 是高斯函数的标准差,决定了Gabor滤波器的带宽。
  • γ \gamma γ 是空间纵横比,决定了Gabor滤波器的椭圆率。当 γ = 1 \gamma = 1 γ=1时,滤波器形状是圆的;当 γ < 1 \gamma < 1 γ<1时,滤波器形状沿平行于 θ \theta θ的方向拉长。

x ′ x' x y ′ y' y 是通过旋转坐标轴得到的坐标,它们与原始坐标 ( x , y ) (x, y) (x,y)的关系为:

( x ′ y ′ ) = ( cos ⁡ θ sin ⁡ θ − sin ⁡ θ cos ⁡ θ ) ( x y ) \begin{pmatrix} x' \\ y' \end{pmatrix} = \begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} (xy)=(cosθsinθsinθcosθ)(xy)

这个旋转矩阵用于将Gabor滤波器的方向调整到 θ \theta θ指定的角度。

Gabor滤波器的一个重要特性是它能够捕捉图像中特定频率和方向的信息。通过调整 λ \lambda λ θ \theta θ σ \sigma σ γ \gamma γ等参数,可以生成一系列具有不同频率和方向选择性的Gabor滤波器,从而用于图像的纹理分析和特征提取。

在实际应用中,通常会使用多个Gabor滤波器(具有不同的参数组合)对图像进行滤波,然后分析滤波后的图像以提取有用的特征信息。

第二种形式高斯函数

1

f ( x , y ) = A exp ⁡ ( − x ′ 2 + γ 2 y ′ 2 2 σ 2 ) f(x, y) =A\exp\left(-\frac{x'^2 + \gamma^2 y'^2}{2\sigma^2}\right) f(x,y)=Aexp(2σ2x′2+γ2y′2)是高斯函数,可由标准正态分布通过几何变换推出。

2

也可以写为
f ( x , y ) = A e − 1 2 ( x ˊ 2 σ x 2 + x ˊ 2 σ y 2 ) f(x, y) =A{\rm e} ^{-{ \tfrac {1}{2} \left( \tfrac { \acute{x} ^{2} } {\sigma_x ^{2}} + \tfrac { \acute{x} ^{2} } { \sigma_y ^{2}} \right) }} f(x,y)=Ae21(σx2xˊ2+σy2xˊ2)

σ x \sigma_x σx σ y \sigma_y σy x x x y y y方向上高斯函数的标准差, γ = σ x / σ y \gamma = \sigma_x /\sigma_y γ=σx/σy表示高斯函数的椭圆率。

滤波器在 x x x方向上的“宽度”可以由标准差 σ x \sigma_x σx控制,而在 y y y方向上的“宽度”可以由另一个标准差 σ y \sigma_y σy控制。但在实际应用中,我们并不直接设置这两个标准差,而是通过设置一个主标准差 σ \sigma σ和一个空间纵横比 γ \gamma γ来间接控制。

空间纵横比 γ \gamma γ定义了滤波器在 x x x y y y方向上的相对“宽度”。当 γ = 1 \gamma = 1 γ=1时,滤波器是圆形的;当 γ < 1 \gamma < 1 γ<1时,滤波器在 y y y方向上的“宽度”比 x x x方向小,形成椭圆形,且椭圆的长轴与 x x x轴平行;当 γ > 1 \gamma > 1 γ>1时,情况相反。

在Gabor滤波器的标准定义中,通常只关注一个标准差 σ \sigma σ和一个空间纵横比 γ \gamma γ。这两个参数共同决定了滤波器在不同方向上的形状和带宽。

3

f ( x , y ) f(x, y) f(x,y)写成关于 x x x y y y的函数,这是一个指数项有交叉项的形式。

z = A e − [ a ( x − μ x ) 2 + 2 b ( x − μ x ) ( y − μ y ) + c ( y − μ y ) 2 ] z = A e^{-\left[a(x - \mu_x)^2 + 2b(x - \mu_x)(y - \mu_y) + c(y - \mu_y)^2\right]} z=Ae[a(xμx)2+2b(xμx)(yμy)+c(yμy)2]

其中,

  • a = cos ⁡ 2 θ 2 δ x 2 + sin ⁡ 2 θ 2 δ y 2 a = \frac{\cos^2\theta}{2\delta_x^2} + \frac{\sin^2\theta}{2\delta_y^2} a=2δx2cos2θ+2δy2sin2θ
  • b = − sin ⁡ ( 2 θ ) 4 δ x 2 + sin ⁡ ( 2 θ ) 4 δ y 2 b = -\frac{\sin(2\theta)}{4\delta_x^2} + \frac{\sin(2\theta)}{4\delta_y^2} b=4δx2sin(2θ)+4δy2sin(2θ)
  • c = sin ⁡ 2 θ 2 δ x 2 + cos ⁡ 2 θ 2 δ y 2 c = \frac{\sin^2\theta}{2\delta_x^2} + \frac{\cos^2\theta}{2\delta_y^2} c=2δx2sin2θ+2δy2cos2θ
    在这里插入图片描述
    在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/53353.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【LabVIEW学习篇 - 18】:人机界面交互设计02

文章目录 错误处理函数简单错误处理器通用错误处理器清楚错误合并错误错误代码至错误簇转换查找第一个错误 鼠标指针 错误处理函数 在LabVIEW中&#xff0c;是通过错误输入簇和错误输出簇来传递错误信息&#xff0c;可以将底层错误信息传递到上层VI。设计人员需要对不同程度的…

2024 数学建模高教社杯 国赛(C题)| 农作物的种植策略 | 建模秘籍文章代码思路大全

铛铛&#xff01;小秘籍来咯&#xff01; 小秘籍团队独辟蹊径&#xff0c;运用等多目标规划等强大工具&#xff0c;构建了这一题的详细解答哦&#xff01; 为大家量身打造创新解决方案。小秘籍团队&#xff0c;始终引领着建模问题求解的风潮。 抓紧小秘籍&#xff0c;我们出发…

高并发内存池(一):项目介绍与定长内存池的实现

目录​​​​​​​ 项目介绍 池化技术 内存池 内存碎片 malloc工作原理 定长内存池 申请内存 释放内存 定位new VirtualAlloc函数 封装VirtualAlloc 定长内存池的最终代码 项目介绍 项目原型&#xff1a;goole的开源项目tcmalloc&#xff08;Thread-Caching Mal…

CIOE中国光博会&电巢科技即将联办“智能消费电子创新发展论坛”

在科技浪潮汹涌澎湃的当下&#xff0c;从通信领域的高速光传输&#xff0c;到消费电子中的高清显示与先进成像技术&#xff0c;光电技术的应用范围不断拓展且日益深化。而AIGC 凭借其丰富的内容供给与个性化反馈能力&#xff0c;正为新一代消费电子及智能穿戴产品开辟崭新的发展…

前端请求的路径baseURL怎么来的 ?nodejs解决cors问题的一种方法

背景&#xff1a;后端使用node.js搭建&#xff0c;用的是express 前端请求的路径baseURL怎么来的 &#xff1f; 前后端都在同一台电脑上运行&#xff0c;后端的域名就是localhost&#xff0c;如果使用的是http协议&#xff0c;后端监听的端口号为3000&#xff0c;那么前端请求…

S3C2440开发板:时钟,PWM定时器控制蜂鸣器发声

时钟 时钟和电源管理模块由三部分组成&#xff1a;时钟控制&#xff0c;USB 控制和电源控制。 S3C2440A 中的时钟控制逻辑可以产生必须的时钟信号&#xff0c;包括 CPU 的 FCLK&#xff0c;AHB 总线外设的 HCLK 以及 APB 总线外设的 PCLK。S3C2440A 包含两个锁相环&#xff08…

VBA进行excel坐标转换

在Excel里利用坐标绘图时&#xff0c;可以比较容易想到采用数据透视表&#xff0c;但是数据透视表生成的图不可更改&#xff0c;因此本案例采用VBA进行坐标变换而不改变原始值来转换图像&#xff0c;即实现图像的左右翻转和上下翻转&#xff0c;如下图所示&#xff0c;选择map的…

fastadmin 文件上传腾讯云

1-安装腾讯云SDK composer require qcloud/cos-sdk-v5 2-腾讯云配置 <?phpnamespace app\common\controller;use Qcloud\Cos\Client; use think\Controller; use think\Db;class Tencent extends Controller {/*** 上传文件* param $config* param $key* return array*/p…

Linux下快速判断当前终端使用的是bash or csh

在Linux下设置环境变量的时候&#xff0c;可能你也遇到过export: Command not found一类的错误。这是因为当前终端使用的不是bash&#xff0c;如何快速判断当前终端使用的是哪种类型的shell呢&#xff1f; echo $0判断shell类型 最简单的方法就是在终端输入echo $0&#xff0…

每日一题,零基础入门FPGA——工程师在线精讲来咯

传送门&#xff1a;zzfpga.com/StudentPlatform/Train/StudentArticleDetails?id149

M3U8工作原理以及key解密视频流详解

文章目录 前言一、M3U8是什么&#xff1f;二、HLS—M3U8的工作原理1.分段视频流2.生成播放列表3.客户端请求和解析4.片段下载和播放 三、.m3u8文件内部是什么样的&#xff1f;四、简单介绍下AES-128算法五、拿到KEY后如何去解密&#xff1f;1.手动解密.ts文件2.前人栽树&#x…

spring security 如何解决跨域的

一、什么是 CORS CORS(Cross-Origin Resource Sharing) 是由 W3C制定的一种跨域资源共享技术标准&#xff0c;其目就是为了解决前端的跨域请求。在JavaEE 开发中&#xff0c;最常见的前端跨域请求解决方案是早期的JSONP&#xff0c;但是JSONP 只支持 GET 请求&#xff0c;这是一…

深度学习从入门到精通——基于unet++算法实现细胞分割

模型定义 import torch from torch import nn__all__ [UNet, NestedUNet]class VGGBlock(nn.Module):def __init__(self, in_channels, middle_channels, out_channels):super().__init__()self.relu nn.ReLU(inplaceTrue)self.conv1 nn.Conv2d(in_channels, middle_channe…

FPGA速度优化

速度优化 文章目录 速度优化前言一、时序优化1.1 减少关键路径上的时序1.1.1 关键路径重组1.1.2 解决扇出问题1.1.3 路径上插入寄存器1.1.4 寄存器平衡1.1.5 并行结构1.1.6 消除代码优先级 总结 前言 速度优化&#xff0c;主要就是设计时序进行优化 吞吐量&#xff1a;每个时…

web渗透:RCE漏洞

RCE漏洞&#xff0c;即远程代码执行漏洞&#xff0c;是一种安全缺陷&#xff0c;它允许攻击者通过网络在目标系统上执行任意代码。一旦成功利用&#xff0c;攻击者可以完全控制受影响的系统&#xff0c;包括读取敏感数据、安装恶意软件、修改系统配置等。RCE漏洞通常发生在应用…

数据结构---双向链表---循环链表---栈

目录 一、双向链表 1.1.创建双向链表 1.2.头插法 1.3.尾插法 1.4.查询节点 1.5.修改节点 1.6.删除节点 1.7.打印节点 1.8.销毁链表 二、循环链表 2.1.单循环链表 2.2.双循环链表 三、栈 3.1.顺序栈 1.创建栈 2.判断栈是否满 3.判断栈是否为空 4.进栈 5.出栈…

SAP 生产订单工序删除状态撤回简介

SAP 生产订单工序删除状态撤回简介 一、业务场景二、处理办法三、系统控制一、业务场景 生产订单正常没有按工序分配物料,系统会自动会把物料分配到第一道工序中 生产订单中的0010工序中对应的组件的栏位被标识,表示有物料分配到了0010的工序中,正常情况下0010的工序被分配…

【微服务】springboot 自定义注解+反射+aop实现动态修改请求参数

目录 一、前言 二、动态修改接口请求参数的场景 2.1 动态修改请求参场景汇总 2.1.1 数据格式标准化 2.1.2 安全需要 2.1.3 参数校验与默认值设定 2.1.4 数据隐私保护 2.1.5 适配不同客户端 2.1.6 统计与监控 2.1.7 高级功能特性 三、springboot 使用过滤器和拦截器动…

Oracle rac模式下undo表空间爆满的解决

文章目录 前言一、确认对应实例的undo表空间二、确认对应实例undo的文件位置三、确认回滚段使用情况四、检查undo segment状态五、创建新的undo表空间并进行切换六、等待原undo表空间segment状态变更为offline七、删除原undo表空间以及数据文件 前言 一、确认对应实例的undo表空…

【云原生】Helm来管理Kubernetes集群的详细使用方法与综合应用实战

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全…