我想计算一个未知函数轮廓内的面积。我的函数值保存在一个Numpy数组中,绘制时如下所示:

我想计算不同高度的横截面积,例如峰高的50%。在

此外,有时特征周围会有噪声,因此会出现其他峰值,这意味着特定高度的轮廓可能会导致创建多个区域。如果从上面看的话,我只对围绕“兴趣峰”的轮廓内的区域感兴趣。在

在上图中,“感兴趣的峰值”是在其顶点有红点的那个。我有一个算法可以识别感兴趣的峰值,所以我总是知道我对哪个特定的峰值感兴趣。在

有人能帮我吗?

非常感谢

我找到了this question,但无法运行给定的示例。我对vs = contour.get_paths()[0].vertices行有问题,列表索引超出范围。在

我还需要修改代码,以确保该区域包含了峰值位置

这是尝试查找区域的代码部分

水平=[0.5,0.75]calc_levels = (radial[x_peak,y_peak])*np.transpose(levels)

cs = plt.contour(Z,X_,radial, levels = calc_levels,colors='k')

plt.clabel(cs, inline=2, fontsize=10)

for loop in range(len(levels)):

vs = None

contour_ = None

contour_ = cs.collections[loop]

vs = contour_.get_paths()[0].vertices

ax.scatter(vs[:,0],vs[:,1],calc_levels[loop],c='red',s=20)

# Compute area enclosed by vertices

a = area(vs)

features['Radial_Area_' + str(levels[loop]*100)] = area(vs)

print("r = {}: a = {}".format(levels[loop],a))

在哪里

^{pr2}$

以这个问题的链接帖子为例。但是,it代码总是在contour_ = cs.collections[loop]行崩溃,显然列表索引超出范围。我已经尝试了我能想到的每一个代码修改来让它工作,但是我现在没有主意了。在

有趣的是,如果我一行一行地运行代码,那么它不会崩溃并提供输出,即使它是一个不太可能的小值

有趣的是,我最终让这段代码在一个简单的例子中工作,其中只包含了感兴趣的峰值。我不知道我改变了上面的代码让它工作。我仍然需要一些帮助,以确定在存在多个峰值的情况下计算哪个区域。如果你能帮我这个忙,请告诉我

python求峰面积_Python:计算contou的面积相关推荐

  1. python求零点_python – 计算大序列的过零点的结果不同

    这个问题源于查看 this关于计算 zero crossings数量的问题的答案.提供了几个解决问题的答案,但NumPy appproach在时间上摧毁了其他问题. 问题:为什么NumPy解决方案提供 ...

  2. python判断三边是否构成三角形并求面积_Python 计算三角形的面积

    参考方法: # 通过用户输入三角形三边长度,并计算三角形的面积 # 已知三角形三边a,b,c,则 # (海伦公式)(p=(a+b+c)/2) # S=sqrt[p(p-a)(p-b)(p-c)] # ...

  3. python求圆面积_python如何求圆的面积

    首先我们要知道圆的面积计算公式:S = πr²,公式中S为所求圆的面积,π为圆周率,r为圆的半径. 示例: # 定义一个方法来计算圆的面积 def findArea(r): PI = 3.142 re ...

  4. python求圆面积_python如何求圆的面积 python求圆的面积方法

    python如何求圆的面积?本篇文章小编给大家分享一下python求圆的面积方法,代码介绍的很详细,小编觉得挺不错的,现在分享给大家供大家参考,有需要的小伙伴们可以来看看. 首先我们要知道圆的面积计算 ...

  5. python求周长_计算三角形的周长和面积

    计算三角形的周长和面积 题目要求 写一段程序,让用户输入三角形的三条边长,如果三条边长不能构成三角形,则提示用户重新输入 如果可以构成三角形,则计算周长和面积 思路分析 对于用户的输入,首先要约定格式 ...

  6. python求对数_python中取对数

    技术 | Python从零开始系列连载(二十七) 我们接着上一期的Python,继续跟大家分享有关Python中常用的数值计算和正则表达式. 运算符 1)数值运算 +.-.*./ % (求余数) // ...

  7. 用python求兀_python求π

    <像计算机科学家一样思考Python(第2版)>--2.5 操作顺序 本节书摘来自异步社区<像计算机科学家一样思考Python(第2版)>一书中的第2章,第2.5节,作者[美] ...

  8. python求积分_python 求定积分和不定积分示例

    求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分 sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi import math import numpy a ...

  9. python求乘积_Python实现求笛卡尔乘积的方法

    Python实现求笛卡尔乘积的方法 本文实例讲述了Python实现求笛卡尔乘积的方法.分享给大家供大家参考,具体如下: 在数学中,两个集合X和Y的笛卡尓乘积(Cartesian product),又称 ...

最新文章

  1. Set Matrix Zeroes leetcode
  2. iOS 7 Searchbar右侧空白
  3. Hive的下载安装,以及配置mysql作为元数据库
  4. webpack的四个核心概念
  5. Scala集合实现WordCount代码实现
  6. # 20175213 2018-2019-2 《Java程序设计》第1周学习总结
  7. POJ1789-Truck History .
  8. Ubuntu下安装Chrome浏览器的两个方法
  9. insert mysql后加where,如何在MySQL Insert語句中添加where子句?
  10. 音视频开发(39)---语音增强
  11. php5.3降5.2,php5.3降级为5.2
  12. 扔掉JUnit,迎接Spock
  13. 阿里消息中间件ONS消息乱序问题(二)
  14. 378.有序矩阵中第K小的元素
  15. 不容错过 | 华为内部资料--成功的项目管理PPT(123页)
  16. 任意类型变量转换成char类型——sprintf函数使用方法
  17. SpringCloud:Gateway网关配置及使用
  18. Nuvoton I91200 语音MCU入门教程
  19. 季琦谈创业:三对矛盾和三个关口[转]
  20. 详解Docker的基本原理与实践操作

热门文章

  1. uniapp H5、app、小程序不同端之间跳转
  2. 2022 IEEE 全球机器人、人工智能与信息技术学术会议(GCRAIT2022)
  3. Android正方形View
  4. 编码规则-驼峰式拼写法
  5. 【修真院java小课堂】ArrayList浅析
  6. 调用华为云GES服务业务面API相关参数的获取
  7. 使用FeignClient注解,进行远程http第三方调用
  8. 【力扣】6. Z 字形变换
  9. 2014河北职称计算机考试大纲,2014职称计算机考试大纲
  10. 路由器与三层交换机有啥区别?