【2022.11.16】Lammps+Python+MATLAB在绘制维诺图时遇到的问题

写在前面

记录一下做毕设时候遇到的问题
大目标是绘制g6®的图片

相关文章：
[1] Zu M , Liu J , Hua T , et al. Density affects the nature of the hexatic-liquid transition in two-dimensional melting of core-softened systems[J]. Physical Review Letters, 2016, 117(8).
[2] V. Nosenko, S.K. Zhdanov, A.V. Ivlev, et al.2D Melting of Plasma Crystals: Equilibrium and Nonequilibrium Regimes[J]. Physical Review Letters, 2009, 103.

顺便附上第二篇论文的paper reading（第一篇直接一对一在老板这边讲的，忘记录了）

组会论文讲解实录：2D Melting of Plasma Crystals，PRL 2009

一些确定会写入源代码中的代码不会展示，因为鄙人搜了一下，发现CSDN上的代码也是会参与查重的，被查到怪麻烦的

绘制g6®执行步骤

1、 导入数据
2、 利用论文1中的方法找到最近邻
3、 利用r进行分类
4、 算出psi6®
5、 利用公式算出g6®

问题1：数据导入

用文本文档无法打开数据文件，转而在Python中进行阅读
代码略

问题2：利用Python选取想要的数据

报错：ValueError: 4 columns passed, passed data had 6 columns
通过仔细查看lamms中的模拟语句，以及在Python中进行校验，