c语言峰值算法_关于算法:将高斯曲线拟合到c ++中的某个直方图峰值-程序员宅基地

技术标签: c语言峰值算法  

关于如何将高斯曲线拟合到柱状图峰值,我有两个问题。我的第一个问题是一个非常基本的问题:

如何将高斯曲线拟合到整个柱状图?这仅仅意味着我必须找出并计算平均值(μ)和偏差(?)并将其放入高斯曲线的公式中?

下面的例子对吗?:假设(举个例子)我有一个带有5个颜色值的图像直方图。在x轴上有5个颜色值,在y轴上有每个颜色值的频率。即。:

值1:1次

值2:4次

值3:7次

值4:3次

值5:2次

现在平均值(μ)为3(μ=3)。

偏差(?)会是0.9(?=0.9)。公式:

现在我用密度函数公式中的这些值来计算高斯曲线?

对吗?不幸的是,我对数学背景有点不太了解。

我的第二个问题有点棘手:这次,我有一个有几个峰值的柱状图,但是我只想把高斯曲线拟合到最高的峰值。所以,我用一个简单的for循环遍历柱状图中的所有bin,在x轴上找到一个强度值(它包含图像的强度)。最大频率(显示在Y轴上)。这将是最高的山峰。但我该如何找出偏差呢?尤其是,因为我没有知道我的计算中应该包括哪些强度值。据我所知,高斯曲线的转折点在μ+?还有μ-?那有助于解决问题吗?

我很抱歉这个问题有点数学化,但我找不到更好的地方问它。我也读了一些类似的话题,但不幸的是,他们最终没有回答我的问题。

谢谢你的帮助!

尊敬的马克

我在这里没有看到编程问题。

当然,这是一个编程问题。这类计算通常是通过计算机编程来解决的。这就是为什么他们被称为计算机的原因。也许你只是看到一个你无法回答的编程问题;-)

对于单个峰值,您的简单方法应该是可以的。[顺便说一句,我认为这是一个与图像相关的编程问题。]

多峰比较困难。解决峰值的过程称为反褶积(http://en.wikipedia.org/wiki/decolution),通常没有唯一的答案。您可能需要决定有多少个峰值,或者一个峰值的最小方差是多少(否则您可以为每个bin创建一个峰值:-)。

我认为只有当所有高斯曲线都有相同的方差时,反褶积才能解决这个问题。

不,但这确实使它更容易。你没有说你的多模性是由什么引起的,但我想多峰大致是等变的。

我认为你估计单个高斯参数的方法是正确的。

对于多个高斯人,您可能希望查找混合模型或更具体地说是高斯人的混合。只需注意几点:

(你可能已经知道)组合中最强的高斯值不一定是最高峰值的高斯值。

如果您有两个具有不同平均值的高斯人的总和,柱状图中的峰值将与平均值不匹配。

使用离散直方图中单个峰值的高度可能不是一个好主意:实际峰值可能位于两个直方图箱之间。[说明:在这种情况下,两个存储箱中的计数可能明显低于以存储箱为中心的峰值。]而且,对于大多数类型的数据,柱状图将包含大量噪声。至少在柱状图上使用某种平均过滤。

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_39836943/article/details/111543378

智能推荐

WCE Windows hash抓取工具 教程_wce.exe -s aaa:win-9r7tfgsiqkf:0000000000000000000-程序员宅基地

文章浏览阅读6.9k次。WCE 下载地址:链接:https://share.weiyun.com/5MqXW47 密码:bdpqku工具界面_wce.exe -s aaa:win-9r7tfgsiqkf:00000000000000000000000000000000:a658974b892e

各种“网络地球仪”-程序员宅基地

文章浏览阅读4.5k次。Weather Globe(Mackiev)Google Earth(Google)Virtual Earth(Microsoft)World Wind(NASA)Skyline Globe(Skylinesoft)ArcGISExplorer(ESRI)国内LTEarth(灵图)、GeoGlobe(吉奥)、EV-Globe(国遥新天地) 软件名称: 3D Weather Globe(http:/_网络地球仪

程序员的办公桌上,都出现过哪些神奇的玩意儿 ~_程序员展示刀,产品经理展示枪-程序员宅基地

文章浏览阅读1.9w次,点赞113次,收藏57次。我要买这些东西,然后震惊整个办公室_程序员展示刀,产品经理展示枪

霍尔信号、编码器信号与电机转向-程序员宅基地

文章浏览阅读1.6w次,点赞7次,收藏63次。霍尔信号、编码器信号与电机转向从电机出轴方向看去,电机轴逆时针转动,霍尔信号的序列为编码器信号的序列为将霍尔信号按照H3 H2 H1的顺序组成三位二进制数,则霍尔信号翻译成状态为以120°放置霍尔为例如不给电机加电,使用示波器测量三个霍尔信号和电机三相反电动势,按照上面所说的方向用手转动电机得到下图① H1的上升沿对应电机q轴与H1位置电角度夹角为0°,..._霍尔信号

个人微信淘宝客返利机器人搭建教程_怎么自己制作返利机器人-程序员宅基地

文章浏览阅读7.1k次,点赞5次,收藏36次。个人微信淘宝客返利机器人搭建一篇教程全搞定天猫淘宝有优惠券和返利,仅天猫淘宝每年返利几十亿,你知道么?技巧分享:在天猫淘宝京东拼多多上挑选好产品后,按住标题文字后“复制链接”,把复制的淘口令或链接发给机器人,复制机器人返回优惠券口令或链接,再打开天猫或淘宝就能领取优惠券啦下面教你如何搭建一个类似阿可查券返利机器人搭建查券返利机器人前提条件1、注册微信公众号(订阅号、服务号皆可)2、开通阿里妈妈、京东联盟、拼多多联盟一、注册微信公众号https://mp.weixin.qq.com/cgi-b_怎么自己制作返利机器人

【团队技术知识分享 一】技术分享规范指南-程序员宅基地

文章浏览阅读2.1k次,点赞2次,收藏5次。技术分享时应秉持的基本原则:应有团队和个人、奉献者(统筹人)的概念,同时匹配团队激励、个人激励和最佳奉献者激励;团队应该打开工作内容边界,成员应该来自各内容方向;评分标准不应该过于模糊,否则没有意义,应由客观的基础分值以及分团队的主观综合结论得出。应有心愿单激励机制,促进大家共同聚焦到感兴趣的事情上;选题应有规范和框架,具体到某个小类,这样收获才有目标性,发布分享主题时大家才能快速判断是否是自己感兴趣的;流程和分享的模版应该有固定范式,避免随意的格式导致随意的内容,评分也应该部分参考于此;参会原则,应有_技术分享

随便推点

O2OA开源企业办公开发平台:使用Vue-CLI开发O2应用_vue2 oa-程序员宅基地

文章浏览阅读1k次。在模板中,我们使用了标签,将由o2-view组件负责渲染,给o2-view传入了两个参数:app="内容管理数据"和name="所有信息",我们将在o2-view组件中使用这两个参数,用于展现“内容管理数据”这个数据应用下的“所有信息”视图。在o2-view组件中,我们主要做的事是,在vue组件挂载后,将o2的视图组件,再挂载到o2-view组件的根Dom对象。当然,这里我们要在我们的O2服务器上创建好数据应用和视图,对应本例中,就是“内容管理数据”应用下的“所有信息”视图。..._vue2 oa

[Lua]table使用随笔-程序员宅基地

文章浏览阅读222次。table是lua中非常重要的一种类型,有必要对其多了解一些。

JAVA反射机制原理及应用和类加载详解-程序员宅基地

文章浏览阅读549次,点赞30次,收藏9次。我们前面学习都有一个概念,被private封装的资源只能类内部访问,外部是不行的,但这个规定被反射赤裸裸的打破了。反射就像一面镜子,它可以清楚看到类的完整结构信息,可以在运行时动态获取类的信息,创建对象以及调用对象的属性和方法。

Linux-LVM与磁盘配额-程序员宅基地

文章浏览阅读1.1k次,点赞35次,收藏12次。Logical Volume Manager,逻辑卷管理能够在保持现有数据不变的情况下动态调整磁盘容量,从而提高磁盘管理的灵活性/boot分区用于存放引导文件,不能基于LVM创建PV(物理卷):基于硬盘或分区设备创建而来,生成N多个PE,PE默认大小4M物理卷是LVM机制的基本存储设备,通常对应为一个普通分区或整个硬盘。创建物理卷时,会在分区或硬盘的头部创建一个保留区块,用于记录 LVM 的属性,并把存储空间分割成默认大小为 4MB 的基本单元(PE),从而构成物理卷。

车充产品UL2089安规测试项目介绍-程序员宅基地

文章浏览阅读379次,点赞7次,收藏10次。4、Dielecteic voltage-withstand test 介电耐压试验。1、Maximum output voltage test 输出电压试验。6、Resistance to crushing test 抗压碎试验。8、Push-back relief test 阻力缓解试验。7、Strain relief test 应变消除试验。2、Power input test 功率输入试验。3、Temperature test 高低温试验。5、Abnormal test 故障试验。

IMX6ULL系统移植篇-系统烧写原理说明_正点原子 imx6ull nand 烧录-程序员宅基地

文章浏览阅读535次。镜像烧写说明_正点原子 imx6ull nand 烧录