我有宇宙射线探测器的能谱.光谱遵循指数曲线,但它会有宽的(也可能是非常轻微的)肿块.显然,数据包含噪声元素.

我正在尝试平滑数据,然后绘制其渐变.

到目前为止,我一直在使用scipy sline函数来平滑它,然后使用np.gradient().

从图中可以看出,梯度函数的方法是找出每个点之间的差异,并且它不会非常清楚地显示肿块.

我基本上需要一个平滑的梯度图.任何帮助都会很棒!

我尝试了2个样条方法:

def smooth_data(y,x,factor):

print "smoothing data by interpolation..."

xnew=np.linspace(min(x),max(x),factor*len(x))

smoothy=spline(x,y,xnew)

return smoothy,xnew

def smooth2_data(y,factor):

xnew=np.linspace(min(x),factor*len(x))

f=interpolate.UnivariateSpline(x,y)

g=interpolate.interp1d(x,y)

return g(xnew),xnew

编辑:尝试数值区分:

def smooth_data(y,xnew

def minim(u,f,k):

""""functional to be minimised to find optimum u. f is original,u is approx"""

integral1=abs(np.gradient(u))

part1=simps(integral1)

part2=simps(u)

integral2=abs(part2-f)**2.

part3=simps(integral2)

F=k*part1+part3

return F

def fit(data_x,data_y,denoising,smooth_fac):

smy,xnew=smooth_data(data_y,data_x,smooth_fac)

y0,xnnew=smooth_data(smy,xnew,1./smooth_fac)

y0=list(y0)

data_y=list(data_y)

data_fit=fmin(minim,y0,args=(data_y,denoising),maxiter=1000,maxfun=1000)

return data_fit

但是,它只是再次返回相同的图形!

python渐变_嘈杂的数据中的渐变,python相关推荐

  1. python读取idx_通过Python从.idx3-ubyte文件或GZIP中提取图像 - python

    我已经使用OpenCV中的facerecognizer创建了一个简单的人脸识别功能.它可以很好地处理人的图像. 现在,我想通过使用手写字符而不是人来进行测试.我遇到了MNIST数据集,但它们将图像存储 ...

  2. python运用在大数据中精准生活_在大数据中“精准”生活 阅读答案

    在大数据中"精准"生活 阅读答案 在大数据中"精准"生活 ①万物皆互联,无处不计算.因为互联网.手机.无线传感器的普及,实时监测.远程协作.SOHO工作.数据管 ...

  3. python运用在大数据中精准生活_在大数据中精准生活阅读理解答案

    在大数据中精准生活阅读理解答案 ①万物皆互联,无处不计算.因为互联网.手机.无线传感器的普及,实时监测.远程协作.SOHO工作.数据管理已成为平常之事,信息像水电一样通过网络供应汩汩传输,计算机上有形 ...

  4. python运用在大数据中精准生活_在大数据中“精准”生活 阅读答案(2014滨州中考试题)...

    在大数据中"精准"生活 ①万物皆互联,无处不计算.因为互联网.手机.无线传感器的普及,实时监测.远程协作.SOHO工作.数据管理已成为平常之事,信息像水电一样通过网络供应汩汩传输, ...

  5. python查看数据大小_科多大数据带你看Python可以列为最值得学习的编程语言

    原标题:科多大数据带你看Python可以列为最值得学习的编程语言 不知道从什么时候开始,这句话开始流行.不过也从侧面反映出 Python 语言的特点:简单.高效. 从近期代表技术趋势的业界报告以及编程 ...

  6. python替代hadoop_Python连接Hadoop数据中遇到的各种坑(汇总)

    最近准备使用Python+Hadoop+Pandas进行一些深度的分析与机器学习相关工作.(当然随着学习过程的进展,现在准备使用Python+Spark+Hadoop这样一套体系来搭建后续的工作环境) ...

  7. python的顶级库_三大用于数据科学的顶级Python库

    Python有许多吸引力,如效率,代码可读性和速度,使其成为数据科学爱好者的首选编程语言.Python通常是希望升级其应用程序功能的数据科学家和机器学习专家的首选. 由于其广泛的用途,Python拥有 ...

  8. python大数据培训机构_学大数据开发需要学习python吗

    大数据现在互联网火热的一个名词,而和大数据关键词最紧密的相信就是Java和python了,在一年以前,Java大数据可能是很多培训机构的宣传标语.而到了2018年,python大数据则成为了潮流,无论 ...

  9. python 脑洞_从说韩语到写Python,这个数据媛的脑洞有毒吧

    你好,Hello,안녕하세요, こんにちは,我是会说四门语言,但是日语可以忽略,韩语也差不多已经忘记,转而写R和Python的语言学迷妹聂大哥. 我要说的脑洞有毒的数据媛就是我计己,哈哈.正式入坑数据 ...

最新文章

  1. 1080 Graduate Admission
  2. 【已解决】关于SQL2008 “不允许保存更改。您所做的更改要求删除并重新创建以下表。您对无法重新创建的标进行了更改或者启用了‘阻止保存要求重新创建表的更改’” 解决方案
  3. php self this parent
  4. 树莓派学习笔记(7):利用bypy实现树莓派NAS同步百度云
  5. edward_mj退役经验帖
  6. GCC如何产生core dump
  7. unity改变图片像素大小_类动森像素画反向转换 (xBR)
  8. linux 文件乱码_RedHatlinux系统虚机启动后/etc/fstab中的文件系统未能挂载上
  9. 记录 PHP 缓存区ob
  10. JavaWeb之分页代码
  11. RobotFramework中实现接口上传文件
  12. Markdown:插入图片和图注
  13. 大恒相机Python版Qt界面(一)
  14. Packet Tracer 5.0建构CCNA实验攻略——帧中继Frame Relay
  15. android ui 切图工具,APP切图标注教程:UI设计切图标注的小工具实用技巧
  16. C#通过函数名字符串执行相应的函数
  17. 论文解读 X-CLIP : Expanding Language-Image Pretrained Models for General Video Recognition
  18. 农家女靠养花赚钱,年收入几十万
  19. MySQL 编译安装错误记录
  20. 人工智能原理复习 | 课程背景

热门文章

  1. 9.3 国际化处理标签
  2. WPF中利用DynamicResource动态资源和资源字典实现软件中英文切换
  3. LVGL学习(一)PC模拟器
  4. AT指令集 CREG、CGREG、CEREG
  5. 计算机图形数学知识点,初三数学知识点归纳整理
  6. 基于C语言开发菜鸟驿站管理系统【100010822】
  7. 模糊理论在机器人传感器中的应用_超声波传感器和激光雷达传感器在机器人避障中的应用...
  8. BIM建模助手【获奖名单公布】2000大奖新鲜出炉
  9. 基于CentOS的Linux云服务器搭建
  10. 压力大睡不着怎么办?