使用Python轻松实现科研绘图

当撰写在学术期刊上发表的文章时,图表的布局和风格应符合预定义的格式要求。这样可以确保该出版物的所有文章都具有一致的风格,并且任何包含的图表在打印时都是高质量的。

Python在科学界广泛使用,并提供了创建科学绘图的好方法。然而,当使用Python中最流行的绘图库之一matplotlib时,默认的图表质量较差,需要进行调整以确保满足要求。更改matplotlib图形的样式可能非常耗时,这就是scienceplots库派上用场的地方。只需几行代码,就可以立即改变图形的外观,而无需花费太多时间来研究如何更改图形的不同部分。

scienceplots该库允许用户创建简单、信息丰富的图表,类似于学术期刊和研究论文中的图表。不仅如此,它还将所需的DPI设置为600(对于某些样式),这通常是出版物的要求,以确保打印出高质量的印刷图形。scienceplots库包含多种样式,包括对多种语言的支持,包括中文和日文等。可以通过下面的链接探索scienceplots库中的全部样式。

本文将探讨如何将一些基本和常见的数据可视化转换为可以包含在科学出版物中的形式。

1.设置scienceplots

在使用scienceplots库创建绘图之前,需要确保计算机上安装了LaTeX。LaTeX是一种排版系统,专为创建技术和科学文档而设计。如果正在使用Google Colab,则可以在单元格中运行以下代码来安装LaTeX。

!sudo apt-get install dvipng texlive-latex-extra texlive-fonts-recommended texlive-latex-recommended cm-super

在设置好LaTeX之后,可以使用pip安装scienceplots库:

pip install SciencePlots

在选择的平台上安装了库和LaTeX之后,就可以导入scienceplots库以及matplotlib

import scienceplots
import matplotlib.pyplot as plt

2.创建用于绘图的虚拟数据

在生成一些图表之前,首先需要创建一些样本数据。我们将在本文的后面部分展示scienceplots库如何处理现实世界的数据。

对于本文的这一部分,现在将使用np.linspace创建一些线性间隔的值,然后对该数据进行一些随机的数学计算。

# 生成x值
x = np.linspace(0, 10, 20)# 生成带有随机噪声的y值
y = np.sin(x)
y2 = np.cos(x)
y3 = y2 * 1.5

一旦创建了数据(或者如果从csv文件中加载数据,则将其加载到pandas中),就可以开始创建绘图了。

3.使用Matplotlib创建带有标记的折线图

我们将首先处理的是折线图。可以通过使用matplotlib.plot()函数并传入xy参数的必需数据来轻松创建此图表。由于我们处理的是从方程中派生的变量,因此有时将这些变量包含在图例中对于读者理解它们可能很方便。

Matplotlib的一个好处是我们可以使用LaTeX方程作为标签。我们只需要用美元符号($)将方程括起来即可。如下是使用scienceplotsmatplotlib创建带有标签的折线图的示例代码:

plt.figure(figsize = (6,6))
plt.plot(x, y, marker='o', label='$y=sin(x)$')
plt.plot(x, y2, marker='o', label='$y=cos(x)$')
plt.plot(x, y3, marker='o', label='$y=y2*1.5$')plt.xlabel('X')
plt.ylabel('Y')
plt.legend()
plt.show()

当运行上述代码时,会得到以下带有标准颜色的基本Matplotlib图像:

图片

应用scienceplots之前的基本Matplotlib折线图

尽管上图看起来可用,但它的质量(dpi和大小)以及样式可能并不完全适合在期刊中发表。

4.应用SciencePlots样式到线型图

要立即转换图像,可以添加一行代码:使用with语句,调用Matplotlib的style.context函数,并允许传入SciencePlots提供的众多可用样式之一。

with plt.style.context(['science', 'high-vis']):plt.figure(figsize = (6,6))plt.plot(x, y, marker='o', label='$y=sin(x)$')plt.plot(x, y2, marker='o', label='$y=cos(x)$')plt.plot(x, y3, marker='o', label='$y=y2*1.5$')plt.xlabel('X Variable (mm)')plt.ylabel('Y Variable')plt.legend()plt.show()

当运行上述代码时,将得到以下更适合于期刊发表的绘图。

应用scienceplots样式后的Matplotlib折线图

该图形简洁(即没有多余的图表元素),并且很容易区分不同的线条。此外,在Jupyter Notebook中查看此图形时,它可能会显示得非常大,即使我们设置了相对较小的图形尺寸。这是因为图形的DPI被设置为600,这通常是许多出版物的要求,并确保图形尽可能清晰。

然后尝试应用另一种样式。这次将使用电气和电子工程师学会(IEEE)的样式。要做到这一点,只需要将high-vis替换为ieee即可改变样式。

with plt.style.context(['science', 'ieee']):plt.figure(figsize = (6,6))plt.plot(x, y, marker='o', label='$y=sin(x)$')plt.plot(x, y2, marker='o', label='$y=cos(x)$')plt.plot(x, y3, marker='o', label='$y=y2*1.5$')plt.xlabel('X')plt.ylabel('Y')plt.legend()plt.show()

当运行上述代码时,将得到以下使用IEEE推荐样式的绘图:

图片

应用scienceplots IEEE风格后的Matplotlib折线图

5.使用SciencePlots绘制直方图

在前面的示例中探讨了如何将样式应用于线型图,也可以将相同的样式应用于其他类型的绘图,接下来尝试如何将这种样式应用于直方图。

首先,使用以下代码创建一个Matplotlib图形,使用一些射线(地质形态的自然放射性测量)数据。为了显示第二个数据集,本文将同样的数据向右调整了20个API单位。

plt.figure(figsize = (6,6))
plt.hist(df['GR'], bins=100, label='GR1', alpha =0.5)
plt.hist(df['GR']+20, bins=100, label='GR2', alpha=0.5)
plt.xlim(0, 150)
plt.xlabel('Gamma Ray')
plt.ylabel('Frequency')
plt.legend()
plt.show()

当运行上面的代码时,将得到以下图形:

图片

射线测量的简单Matplotlib直方图

可以注意到它使用了Matplotlib的标准样式,并且看起来非常基础,两个数据集互相重叠,这导致一些信息被隐藏。现在查看一下IEEE样式如何改变这种情况:

with plt.style.context(['science', 'ieee']):plt.figure(figsize = (6,6))plt.hist(df['GR'], bins=100, label='GR1')plt.hist(df['GR']+20, bins=100, label='GR2')plt.xlim(0, 150)plt.xlabel('Gamma Ray')plt.ylabel('Frequency')plt.legend()plt.show()

当运行上述代码时,将得到以下使用IEEE样式的图形。然而,第二个射线数据集仍然遮挡了第一个数据集。

应用scienceplots IEEE样式后的射线测量的Matplotlib直方图

我们对SciencePlots库能够处理任何重叠并自动应用透明度抱有过高的期望,为此需要付出一点努力,需要为每个数据集添加alpha参数即可。

with plt.style.context(['science', 'ieee']):plt.figure(figsize = (6,6))plt.hist(df['GR'], bins=100, label='GR1', alpha=0.5)plt.hist(df['GR']+20, bins=100, label='GR2', alpha=0.5)plt.xlim(0, 150)plt.xlabel('Gamma Ray')plt.ylabel('Frequency')plt.legend()plt.show()

当运行带有alpha更改的上述代码时,将得到以下图形:

图片

应用scienceplots并为数据集添加透明度后的射线测量值的直方图

在这个图形中,可以看到两个数据集的条形的变化,建议检查所需发表的期刊的样式指南,以确保使用透明度是可接受的。

6.应用Science Plots到Seaborn图形

不仅可以将scienceplots库中的样式应用于matplotlib图形,还可以将其应用于Seaborn图形。这是因为Seaborn基于matplotlib代码开发

在创建图形时,有时Seaborn提供了比matplotlib更简单的绘图方式。例如,当有一个基于文本的分类变量时,我们希望能够绘制该变量,而不必为每个类别添加单独的散点图。

在这个例子中,有一些中子孔隙度和堆积密度数据,这是常见的测井测量。对于每个测量,还有一个岩性类别,这个数据集来源于Force 2020 Xeek机器学习竞赛数据集。

首先,需要在笔记本中导入seaborn:

import seaborn as sns

导入seaborn库后,可以使用以下代码创建散点图。

plt.figure(figsize=(6, 6))
sns.scatterplot(data=df, x='NPHI', y='RHOB', hue='LITH', s=10)plt.ylabel('Bulk Density (RHOB) - g/cc')
plt.xlabel('Neutron Porosity (NPHI) - dec')
plt.ylim(2.8, 1.4)
plt.show()

当运行上面的代码时,会得到以下数据按不同岩性着色的散点图:

图片

使用Seaborn生成的基本中子-密度交叉图

还需要确保样式适合预期的期刊,并且颜色对所有读者都可访问和使用。要应用scienceplots样式,可以使用与之前相同的语法:

with plt.style.context(['science', 'ieee']):plt.figure(figsize=(10, 8))sns.scatterplot(data=df, x='NPHI', y='RHOB', hue='LITH', s=10)plt.ylim(2.8, 1.4)plt.title('RHOB vs NPHI Crossplot')plt.show()

当运行上面的代码时,会得到以下具有改进样式的图,包括新的调色板。

图片

Seaborn散点图,使用scienceplots样式显示堆积密度与中子孔隙度,按岩性变化着色

选择图形的调色板可能会很棘手且耗时,然而它可以使图形对具有视觉问题的读者更易访问和理解。此外,一些颜色在黑白打印时可能不容易区分,可以考虑为不同类别分配不同的形状,这对于从实验室过程中获得的小型数据集中尤为重要。

综上,本文探讨了如何快速将基本的matplotlib图形转化为可以轻松添加到科学出版文章中的形式。这些图形可能仍然需要进一步调整,但通过使用scienceplots库,可以实现大部分需求。另外,建议查看所选择期刊的作者工具包,以确保创建的图形符合要求的标准。

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

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

相关文章

mysql中的各种日志文件redo log、undo log和binlog

mysql中的各种日志文件redo log、undo log和binlog mysql中的各种日志文件redo log、undo log和binlog1.MySQL日志文件类型2.redo log日志2.1 作用2.2工作原理:2.3详解 3.undo log日志4.binlog日志5.总结 mysql中的各种日志文件redo log、undo log和binlog 1.MySQL…

keras转onnx,TensorFlow转tf.keras.models.load_model,onnx精度转换

参考: https://blog.csdn.net/Deaohst/article/details/126864267 转onnx 别直接转onnx。 先转PB: import tensorflow as tfmodel_path ./models/model.h5 # 模型文件 model tf.keras.models.load_model(model_path) model.sa…

jbase虚拟M层的设计

对于只是自己产品内部使用的打印程序来说(比如打印收费单,打印结算单等),打印逻辑写在js,获取其他层都是没毛病的。但是对于类型检验报告这种打印来说,打印格式控制逻辑写在js层是百分百不行的。因为检验报…

合众汽车选用风河Wind River Linux系统

导读合众新能源汽车股份有限公司近日选择了Wind River Linux 用于开发合众智能安全汽车平台。 合众智能安全汽车平台(Hozon Automo-tive Intelligent Security Vehicle Plat-form)是一个面向高性能服务网关及车辆控制调度的硬件与软件框架,将于2024年中开始投入量产…

线性表的概念

目录 1.什么叫线性表2.区分线性表的题 1.什么叫线性表 线性表(linear list)是n个具有相同特性的数据元素的有限序列。 线性表是一种在实际中广泛使用的数据结构,常见的线性表:顺序表、链表、栈、队列、字符串… 线性表在逻辑上是…

【面试】测试/测开(未完成版)

1. 黑盒测试方法 黑盒测试:关注的是软件功能的实现,关注功能实现是否满足需求,测试对象是基于需求规格说明书。 1)等价类:有效等价类、无效等价类 2)边界值 3)因果图:不同的原因对应…

《硅基物语.AI写作高手:从零开始用ChatGPT学会写作》《从零开始读懂相对论》

文章目录 《硅基物语.AI写作高手:从零开始用ChatGPT学会写作》内容简介核心精华使用ChatGPT可以高效搞定写作的好处如下 《从零开始读懂相对论》内容简介关键点书摘最后 《硅基物语.AI写作高手:从零开始用ChatGPT学会写作》 内容简介 本书从写作与ChatG…

如何在Docker部署Draw.io绘图工具并远程访问

文章目录 前言1. 使用Docker本地部署Drawio2. 安装cpolar内网穿透工具3. 配置Draw.io公网访问地址4. 公网远程访问Draw.io 前言 提到流程图,大家第一时间可能会想到Visio,不可否认,VIsio确实是功能强大,但是软件为收费&#xff0…

matlab二维曲面散点图插值方法

在 MATLAB 中,你可以使用以下函数进行二维曲面散点插值: griddata: 该函数可以在散点数据上进行二维插值,生成平滑的曲面。它支持多种插值方法,包括三次样条插值、最近邻插值、线性插值和自然邻近法插值。 scatteredInterpolant:…

大数据基础设施搭建 - Hadoop

文章目录 一、下载安装包二、上传压缩包三、解压压缩包四、配置环境变量五、测试Hadoop5.1 测试hadoop命令5.2 测试wordcount案例5.2.1 创建wordcount输入文本信息5.2.2 执行程序5.2.3 查看结果 六、分发压缩包到集群中其他机器6.1 分发压缩包6.2 解压压缩包6.3 配置环境变量 七…

Go ZIP压缩文件读写操作

创建zip文件 golang提供了archive/zip包来处理zip压缩文件,下面通过一个简单的示例来展示golang如何创建zip压缩文件: func createZip(filename string) {// 缓存压缩文件内容buf : new(bytes.Buffer)// 创建zipwriter : zip.NewWriter(buf)defer writ…

flowable消息事件

一&#xff1a;启动事件 定义消息。 引用消息。 <startEvent id"msgStart" name"消息启动事件" isInterrupting"true"><messageEventDefinition messageRef"myMsgStart"></messageEventDefinition> </startE…

【ROS导航Navigation】五 | 导航相关的消息 | 地图 | 里程计 | 坐标变换 | 定位 | 目标点和路径规划 | 激光雷达 | 相机

致谢&#xff1a;ROS赵虚左老师 Introduction Autolabor-ROS机器人入门课程《ROS理论与实践》零基础教程 参考赵虚左老师的实战教程 一、地图 nav_msgs/MapMetaData 地图元数据&#xff0c;包括地图的宽度、高度、分辨率等。 nav_msgs/OccupancyGrid 地图栅格数据&#…

reactive和effect,依赖收集触发依赖

通过上一篇文章已经初始化项目&#xff0c;集成了ts和jest。本篇实现Vue3中响应式模块里的reactive方法。 前置知识要求 如果你熟练掌握Map, Set, Proxy, Reflect&#xff0c;可直接跳过这部分。 Map Map是一种用于存储键值对的集合&#xff0c;并且能够记住键的原始插入顺…

【6】Spring Boot 3 集成组件:knift4j+springdoc+swagger3

目录 【6】Spring Boot 3 集成组件&#xff1a;knift4jspringdocswagger3OpenApi规范SpringFox Swagger3SpringFox工具&#xff08;不推荐&#xff09; Springdoc&#xff08;推荐&#xff09;从SpringFox迁移引入依赖配置jAVA Config 配置扩展配置&#xff1a;spring securit…

Linux安装DMETL5与卸载

Linux安装DMETL5与卸载 环境介绍1 DM8数据库配置1.1 DM8数据库安装1.2 初始化达梦数据库1.3 创建DMETL使用的数据库用户 2 配置DMETL52.1 解压DMETL5安装包2.2 安装调度器2.3 安装执行器2.4 安装管理器2.5 启动dmetl5 调度器2.6 启动dmetl5 执行器2.7 启动dmetl5 管理器2.8 查看…

计算机毕业设计选题推荐-一周穿搭推荐微信小程序/安卓APP-项目实战

✨作者主页&#xff1a;IT研究室✨ 个人简介&#xff1a;曾从事计算机专业培训教学&#xff0c;擅长Java、Python、微信小程序、Golang、安卓Android等项目实战。接项目定制开发、代码讲解、答辩教学、文档编写、降重等。 ☑文末获取源码☑ 精彩专栏推荐⬇⬇⬇ Java项目 Python…

系列二十六、idea安装javap -c

一、概述 javap -c是一个能够将.java文件反编译为.class文件的指令&#xff0c;例如我在idea中编写了一个Car.java文件&#xff0c;我想看看这个类被编译后长什么样的&#xff0c;就可以使用该指令进行查看。 二、配置 2.1、 Java Bytecode Decompiler File>Settings>Pl…

TMS320F28335使用多个串口时,SCIRXST Register出现错误

TMS320F28335使用多个串口时&#xff0c;SCIRXST Register出现错误 void ClearErrorState(void) {if((SciaRegs.SCIRXST.bit.FE 1)||(SciaRegs.SCIRXST.bit.BRKDT 1)){SciaRegs.SCICTL1.bit.SWRESET 0;SciaRegs.SCICTL1.bit.SWRESET 1;}if((ScibRegs.SCIRXST.bit.FE 1)||(S…

报错资源不足,k8s使用containerd运行容器修改挂载点根目录换成/home

运行k8s一段时间发现存储不足报错 发现这里用的是根路径的挂载&#xff0c;修改一下