python中用于绘制各种图形、标注文本_在python中的图形上绘制常量文本-问答-阿里云开发者社区-阿里云...

我确定我可能在问一个愚蠢的问题，但找不到与我相同的问题。

我的朋友帮我写了一段代码，分析给出的数据并用趋势线将其绘制出来，我想在图的右上方添加一行文本，并在图上打印出其他内容，以表明它是什么文件（在代码的其他地方以常量形式编写）。

import numpy as np

import matplotlib.pyplot as plt

from scipy.optimize import curve_fit

from numpy import exp, loadtxt, pi, sqrt, random, linspace

from lmfit import Model

import glob, os

## Define gaussian

def gaussian(x, amp, cen, wid):

"""1-d gaussian: gaussian(x, amp, cen, wid)"""

return (amp / (sqrt(2\*i) * wid)) * exp(-(x-cen)\*2 / (2\*id\*2))

## Define exponential decay

def expdecay(x, t, A):

return A\*xp(-x/t)

## Define constants

fileToRun = 'Run15'

folderId = '\\'

baseFolder = 'C:'+folderId+'Users'+folderId+'ControlRoom6'+folderId+'Documents'+folderId+'PhD'+folderId+'Ubuntu-Analysis-DCF'+folderId+'DCF-an-b+decay'+folderId+'dcp-ap-27Al'+folderId+''

prefix = 'DECAY_COINC'

stderrThreshold = 10

minimumAmplitude = 0.1

approxcen = 780

MaestroT = 18

## Define paramaters

amps = []; ampserr = []; ts = []

folderToAnalyze = baseFolder + fileToRun + '\\'

## Gets number of files

files = []

os.chdir(folderToAnalyze)

for file in glob.glob(prefix + "\*Spe"):

files.append(file)

numfiles = len(files)

if numfiles<=1:

print('numfiles is {0}, minimum of 2 is required'.format(numfiles))

raise SystemExit(0)

## Generate the time array

for n in range(0, numfiles):

## Print progress

print('\rFile {0} / {1}'.format(n+1, numfiles), end='')

## Load text file

x = np.linspace(0, 8191, 8192)

fullprefix = folderToAnalyze + prefix + str(n).zfill(3)

y = loadtxt(fullprefix + ".Spe", skiprows= 12, max_rows = 8192)

## Make figure

fig, ax = plt.subplots(figsize=(15,8))

fig.suptitle('Coincidence Detections', fontsize=20)

plt.xlabel('Bins', fontsize=14)

plt.ylabel('Counts', fontsize=14)

## Plot data

ax.plot(x, y, 'bo')

ax.set_xlim(600,1000)

## Fit data to Gaussian

gmodel = Model(gaussian)

result = gmodel.fit(y, x=x, amp=8, cen=approxcen, wid=1)

## Plot results and save figure

ax.plot(x, result.best_fit, 'r-', label='best fit')

ax.legend(loc='best')

texttoplot = result.fit_report()

ax.text(0.02, 0.5, texttoplot, transform=ax.transAxes)

plt.close()

fig.savefig(fullprefix + ".png", pad_inches='0.5')

## Print progress

if n==numfiles-1:

print('\rDone')

## Append to list if error in amplitude and amplitude itself is within reasonable bounds

if result.params['amp'].stderr < stderrThreshold and result.params['amp'] > minimumAmplitude:

amps.append(result.params['amp'].value)

ampserr.append(result.params['amp'].stderr)

ts.append(MaestroT\*)

## Plot decay curve

fig, ax = plt.subplots()

ax.errorbar(ts, amps, yerr= 2\*p.array(ampserr), fmt="ko-", capsize = 5, capthick= 2, elinewidth=3, markersize=5)

plt.xlabel('Time', fontsize=14)

plt.ylabel('Peak amplitude', fontsize=14)

## Fit decay curve

emodel = Model(expdecay)

decayresult = emodel.fit(amps, x=ts, weights=1/np.array(ampserr), t=150, A=140)

ax.plot(ts, decayresult.best_fit, 'r-', label='best fit')

## Add text to plot

plottext = '{filetoRun}\n'

plottext = 'N: {0} / {1}\n'.format(len(ts), numfiles)

plottext += 't: {0:.2f} ± {1:.2f}\n'.format(decayresult.params['t'].value, decayresult.params['t'].stderr)

plottext += 'A: {0:.2f} ± {1:.2f}\n'.format(decayresult.params['A'].value, decayresult.params['A'].stderr)

plottext += 'Reduced $χ^2$: {0:.2f}\n'.format(decayresult.redchi)

ax.text(0.5, 0.55, plottext, transform=ax.transAxes, fontsize=14)

plt.show()

## Save figure

fig.savefig(folderToAnalyze + "A_" + prefix + "_decayplot.pdf", pad_inches='0.5')

我想要它（在这种情况下，Run15显示在该图中显示“ N = 28/50”的上方）。我已经尝试了方括号和plt.text的各种组合，但是说实话，我真的不知道自己在做什么，这对我来说是全新的。它不会出现这样的错误，只输出没有所需文本的图形

问题来源：stackoverflow

python中用于绘制各种图形、标注文本_在python中的图形上绘制常量文本-问答-阿里云开发者社区-阿里云...相关推荐

python合法标识符 40xl_Python基础测验（答案篇）-问答-阿里云开发者社区-阿里云
一.填空题 Python使用符号 # 标示单行注释:以缩进对齐划分语句块. Python序列类型包括字符串 . 列表 . 元组三种: _字典_是Python中唯一的映射类型. Python序列 ...
python中str用法_python中的str()不能直接用吗 -问答-阿里云开发者社区-阿里云
str函数是Python的内置函数,它将参数转换成字符串类型,即人适合阅读的形式. 其语法格式为 1 str(object) 返回值: 返回object的字符串形式使用示例无参调用当str()函 ...
python如何实时捕捉cmd显示_如何从Python脚本中捕获Python解释器和/或CMD.EXE的输出？ -问答-阿里云开发者社区-阿里云...
如果您正在谈论python解释器或CMD.exe,它是您脚本的"父",那么不可能.在每个类似POSIX的系统中(现在你正在运行Windows,看起来可能有一些我不知道的怪癖,YMM ...
python语言打小数点_如何在python中打小数点-问答-阿里云开发者社区-阿里云
python提供了三种浮点值:内置的float与complex类型,以及标准库的decimal.Decimal类型. float类型存放双精度的浮点数,具体取值范围依赖于构建python的c编译器,由 ...
python键盘怎么输入双引号_python中怎么输入引号 -问答-阿里云开发者社区-阿里云...
Python中的引号可分为单引号.双引号和三引号. 在Python中我们都知道单引号和双引号都可以用来表示一个字符串,比如 str1 = 'python' str2 = "python&qu ...
python中如何输出中文_python中怎么输出中文-问答-阿里云开发者社区-阿里云
方法一: 用encode和decode 如: ? 1 2 3 4 5 6 7 8 9 10 11 import os.path import xlrd,sys Filename='/home/tom/ ...
python if else用法同一行_在Python的同一行中使用if else for和del吗？-问答-阿里云开发者社区-阿里云...
我有一个列表,其项目是可变长度的列表. 如果这些可变长度的列表项超过此长度,则需要将其截断为特定长度(x). 我做了这个小的功能. def truncateList(batch_, trim_len) ...
python中右对齐_python中如何右对齐-问答-阿里云开发者社区-阿里云
例如,有一个字典如下: dic = { "name": "botoo", "url": "http://www.123.com&q ...
python中tmp什么意思_python中temp是什么意思-问答-阿里云开发者社区-阿里云
==tempfile 模块== [Example 2-6 #eg-2-6] 中展示的 tempfile 模块允许你快速地创建名称唯一的临时文件供使用. ====Example 2-6. 使用 temp ...

python中用于绘制各种图形、标注文本_在python中的图形上绘制常量文本-问答-阿里云开发者社区-阿里云...

python中用于绘制各种图形、标注文本_在python中的图形上绘制常量文本-问答-阿里云开发者社区-阿里云...相关推荐

最新文章

热门文章