红外谱图数据分析系统

[align=center][b]个人简介[/b][/align] [img=,640,477]http://ng1.17img.cn/bbsfiles/images/2017/08/201708091048_01_2984502_3.jpg[/img] 曾仲大，男，博士，现任大连达硕信息技术有限公司总经理。 曾博士师承梁逸曾教授，2006年获得工学博士学位，主要从事化学计量学基础算法研究，以及色、质、光谱等分析技术在制药、烟草和代谢组学等复杂体系分析中的应用及其数据分析挖掘等。近年来在大数据的分析与应用方面亦有涉猎。 曾博士先后工作于香港理工大学、澳洲RMIT大学、Monash大学，以及中国科学院大连化学物理研究所。迄今已发表SCI论文40余篇，在2013-2016近三年时间里，以第一作者或合作者在美国分析化学杂志发表7篇研究论文，同时获邀为TrAC等权威期刊撰写化学计量学及化学数据分析处理方面的综述。 曾博士曾获得中国科学院大连化学物理研究所“所百人”引进人才计划，大连“海创工程”计划、高层次人才创新创业支持计划、新兴技术创新成长计划，以及国家人社部高层次海归人才创业计划的支持。公司主要提供复杂化学与生物数据分析服务，数据挖掘软件产品开发，以及个性化数据应用的整体解决方案。[b]人生格言：[/b]有志者，事竟成。[align=center]~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~[/align][b][color=#7030A0] 随感：[/color][color=#7030A0]“[/color][color=#7030A0]我与近红外的故事[/color][color=#7030A0]”[/color][color=#7030A0]征文近一年了，看过许多老师情真意切的表达，真是把乐趣融入到了近红外的研究与应用之中，也更加深切地感受到同行们对国内近红外发展的使命感和责任感。而自己与近红外的故事，几次动笔却都没能写下几个字。时间肯定不是借口，惰性真是害人啊。好在拖到春节，总算能静下心来了。就像与近红外的相遇相知，既是机缘巧合，更是某种必然吧。[/color][/b] 初识近红外，都是博士毕业一年以后的事了。那时已经在香港理工大学周福添教授课题组从事博士后研究一年多了，主要方向还是老本行-化学计量学基础算法研究，解决中药和代谢组学等复杂体系分析中的数据处理问题，从GC-MS，[url=https://insevent.instrument.com.cn/t/Yp][color=#3333ff]LC-MS[/color][/url]到中药指纹与药物活性关系。一次Daniel MOK博士找到我，询问是否有意愿到陈新滋院士课题组从事中药质量分析与鉴别方面的工作，陈院士那时是理大副校长(后任香港浸会大学校长，现受聘中山大学教授、学委会主任)，研究组的条件与学术水准自不必说，就这样幸运地开始了近二年的近红外数据分析之旅。 对香港熟悉的朋友一定对其大街小巷的名贵中药材印象深刻，尤其是弥墩道，应该是内地赴港旅游人士的必经之地吧，一是去旺角购买电子产品的旅游大巴必定经过这里，另一方面则是这条大道两旁大大小小的中药材店。记得第一次见到时，很是疑惑哪来的那么多冬虫夏草、燕窝和野生人参？说回到陈院士负责的这个研究课题，由香港赛马会中药研究院提供500万研究经费，对包括上述中药，以及石斛、灵芝、阿胶等在内的30味名贵中药材进行质量鉴别分析和研究，目的是帮助那些大街小巷的药材经销店铺，中间批发商，甚至普通消费者，以快速、经济、简便的方法识别药材真假，甚至质量等级。这些药材大多价格不菲，若能够有效识别真假，其商用价值可想而知！顺便一提，香港赛马会中药研究院很多年前已经解散，个中原因无法深究，但在目前国家大力践行中医药研究开发与应用的今天，这也算是一件憾事吧，包括设想中的香港国际中医药中心。 说到这里，近红外分析可以派上用场了！无论是十年前，还是十年后的今天，应没有什么分析技术比近红外更适合完成这项使命，综合考虑时间效率、分析成本，亦或是平衡多重因素影响下定性定量分析结果的准确性！记得当时我们使用的是FOSS公司的XDS快速含量分析仪(Type XM 1100Series)，以及Polychromix手持式近红外分析仪(Model:1600-2400)。由于项目定位于实际应用，需要适应不同场合下的快速分析，对数据分析本身的要求同样也是比较高的，比如涉及模型传递，尽可能简化数据分析的过程及对使用者的要求，亦确保结果的准确可靠性。基于此编写了功能完备的近红外数据分析软件系统，一站式地完成近红外数据分析的完整流程，从各种各样的预处理方法到特征选择，再到定性定量模型的构建、评价与验证预测，以及模型传递等。[img=,587,242]http://ng1.17img.cn/bbsfiles/images/2017/08/201708091049_01_2984502_3.jpg[/img] 说实在的，那时对化学计量学的多元校正方法并不是特别熟悉，我的整个硕士和博士研究，都是多元分辨方向，也就是如何从中药和烟草等复杂体系分析的联用仪器数据中，发展“数学分离”的方法，获取化学纯组分的定性定量信息，即纯组分的光谱和色谱信息。幸运的是，得益于在梁逸曾教授研究组六年时间里耳濡目染的学习，比如许青松教授对统计分析的讲解，杜一平教授的QSAR研究等等，使得我无论对复杂数据的理解，还是化学计量学方法的应用与发展，都有足够基础支持我去解决近红外数据分析中遇到的各种问题。在香港的几年时间里，梁教授每年也都会利用假期去香港一段时间，与香港同行合作交流化学计量学及其应用方面的成果，更是继续指导我解决研究中遇到的实际难题。每每想到这些，总会浮现与恩师相处过程中的点点滴滴。至于上面提到的中药质量分析研究项目，我们对包括阿胶、珍珠、川贝母、藏红花、黄连在内的多味中药进行了深入分析研究，获得了非常不错的结果，陈院士对此也给予了很高的评价。很清楚地记得因此第一次上了电视新闻，是香港亚洲卫视针对我们使用近红外分析技术，如何快速识别真假中药，及其质量等级的采访报道。当然，这些研究很多也是和理工大学的同事，以及杨大坚教授(现任重庆市中药研究院院长)、董玮玮博士等一起完成的，我主要负责数据分析，以及数据软件产品开发与实现方面的工作。[img=,574,238]http://ng1.17img.cn/bbsfiles/images/2017/08/201708091050_01_2984502_3.jpg[/img] 离开香港后，很长一段时间内都没有与近红外分析有直接的关联。先是在Philip Marriott 教授课题组做research fellow，从事全二维色谱数据分析方面的工作，主要方向是全二维分离的模拟、预测，以及化学计量学新方法的发展。2012年回国后则作为引进人才，在中科院大连化物所许国旺教授研究组，从事代谢组学数据分析与高分辨[url=https://insevent.instrument.com.cn/t/Yp][color=#3333ff]LC-MS[/color][/url][i][sup]n[/sup][/i]数据处理新算法的研究等。看似这些工作与近红外分析不怎么挨着边，但老实说，同其他研究一样，数据分析也是一通百通的事！数据来源与数据结构可能不一样，数据背景与数据分析结果，以及数据处理方法亦可能存在差别，但数据分析的本质却是高度一致的，无论是色谱分离的模拟，亦或是代谢小分子标志物的发现！从这个意义上来说，也算是一直在这个圈子吧。 近红外技术的发展，面临非常多的机会，无论从国内快检还是工业智能化的需要来看，还是从国外近红外发展的轨迹来看。然而近红外分析更广阔的应用，仍有一系列需要解决的难题，这其中当然包括仪器硬件的小型化、便携式，以及智能化与场景化。但从数据及数据分析的角度来说，快速、准确的模型构建，模型的通用性、更新及转换等仍是需要加以研究的内容。基于此，离开化物所后创办的大连达硕信息技术有限公司，第一个数据产品“魔力”，便专注近红外数据的分析，这也算是真正走在了近红外技术与数据分析的商业应用之路上。希望能够以智慧化、便捷化的方式，分析挖掘科学研究与工业应用中的海量数据。无论对于近红外分析的初入者，还是有了相当经验的人员，一旦采集到数据，便能快速得到好用的模型及结果，这也是目前非常欠缺的，主要原因就在于近红外数据分析的过程长，可变因素多，涉及的算法也很多，传统上要快速得到一个好用的模型并不容易。尽管大多数研究者并没有把数据分析提升到特别核心的位置，但其价值显而易见，甚至在某些方面可与硬件本身相得益彰，弥补硬件的物理劣势！ 另一方面，近红外分析以其简单方便的前处理，加上非常快速的数据采集方式，使得数据的获取，甚至大数据的积累顺理成章。然而即使对同一组数据，不同的研究者亦极有可能得到完全不同，甚至相反的分析结果或结论，即使在固定分析方法的情况下！这是一个容易被忽视，却又至关重要的问题，否则不管如何将近红外分析的硬件评价，以及实验测试全过程标准化，也无法得到可相互比较的结果。数据“横看成岭侧成峰”的魅力，不应是由于数据分析方法或人员的不同导致，而是数据背景的属性差异或者数据分析目的的不同产生。基于此，我们也正采用近红外数据分析的通用准则，使用粒子群等最优化的方法，开发全新的近红外数据分析软件产品，自动优选数据分析算法，以及方法的使用顺序，并全局优化方法的参数。这样我们获得数据后，只需按照标准化的流程一步一步走，便可获得最优的数据分析模型与模型结果。从而使得近红外数据的分析，如同实验分析一样，结果的重现性与可比性也就不再是个问题。避免像现在这样，往往是漫无目的的数据探索，耗费漫长时间也不一定能得到合适好用的模型！这无论在研究中，还是在工业生产中，都是需要花大力气迎接的挑战。在这一过程中，得到了袁洪福教授、吴海龙教授、邵学广教授、杜一平教授、褚小立教授、闵顺耕教授等诸多老师的大力支持与帮助。从老师们关切的眼神中，能读懂那份殷殷之情，也唯有努力做点事情，为国内近红外的发展做些有益的工作，方不负此情。近红外分析能做的事情很多，近红外数据分析如是，尤其站在移动互联时代，站在大数据分析挖掘的视角与高度。近红外有其自身特有的巨大优势-本身就是物联网中的一个 绝佳传感器！从这个意义上来说，近红外分析代表着某种未来，只是通往未来的路上，还需要我辈站在前辈的肩膀上，不断付出智慧和汗水。 “师者也，教之以事而喻诸德也。”，数据分析之路上，深深地烙上了梁逸曾教授的影响。亦师亦友者，感恩、深切缅怀您。

[b]作者：[/b]丁延伟，[color=#2d374b]中国科学技术大学理化科学实验中心副主任。[/color]众所周知，通过热分析/红外光谱联用技术，在得到样品在温度变化过程中质量或热效应的信息的同时还可以得到逸出气体的官能团随温度变化的信息。此处所指的热分析部分主要指热重法（TG）、热重-差热分析（TG-DTA）和热重-差示扫描量热法（DSC），这部分的数据分析将在之后的相关文章中进行介绍，限于篇幅本部分不做进一步展开。在许多专著和教材中，对于红外光谱的理论、实验及数据分析已有更为详尽的介绍，本文中对于这些内容也不做过多的介绍。本文中重点介绍由热分析/红外光谱联用技术得到的数据分析方法，为了叙述方便，将重点介绍由热重/红外光谱（TG/IR）联用实验所得到的数据分析过程。概括来说，由TG/IR可以得到样品在不同温度或时间下得到的所有红外光谱的[color=#000000]三维图（通常称为逸出气体光谱图，简称EGS），还可以得到[/color][color=#000000]逸出气体剖面图（即在实验期间得到的由样品释放出的气体总量为时间或温度的函数曲线，通常通过应用Gram-Schmidt重建算法根据记录得到的红外光谱数据计算出来的剖面图，简称EGP)[/color][color=#000000]、官能团剖面图（即用来表示在热分析实验过程中释放的气体量中特定的化学官能团随测量时间或的变化关系，通常通过对实验过程中所选光谱区域上的红外光谱数据的吸光值积[/color]分来得到该剖面图，简称FGP）、特定气体剖面图（即当选定的光谱区域中含有特定气体（如氨或一氧化碳）时，由软件所得到的一种特殊类型的官能团剖面图，简称SGP）。在TG/IR分析中，通常需要得到EGS、EGP和FGP图。由于在之前的《热分析/红外光谱联用数据分析时三维图的作图方法》一文中已经介绍了EGS的作图法，因此在本部分内容中以美国Perkin-Elmer公司的TG/IR仪所附带的Timebase软件为例介绍EGP和FGP的分析方法。首先，在Timebase软件中打开需要分析的数据的原始文件（图1）。[align=center] [img=,560,301]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170951414224_4507_3224499_3.jpg!w560x301.jpg[/img] [/align][align=center]图1 打开原始文件[/align]图1中左上部分的曲线为任一时刻得到的实时红外光谱谱图，右上部分蓝色的三维图为实验过程中不同时刻/温度下的所有的红外光谱图，下半部分即为EGP曲线。EGP曲线的导出方法为：（1）选中该文件（图2）；[align=center][img=,188,103]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170951540346_4853_3224499_3.jpg!w188x103.jpg[/img][/align][align=center]图2 选中GS曲线[/align]（2）选中File菜单下的Save As选项，会弹出图4所示的窗口，点击保存即得到Excel格式的GS曲线的数据文件（图4）。导出的文件在Origin、Excel等文件中可以直接作图。[align=center][img=,176,258]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952031126_9770_3224499_3.jpg!w176x258.jpg[/img][/align][align=center]图3[/align][align=center][img=,560,431]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952153848_188_3224499_3.jpg!w560x431.jpg[/img][/align][align=center]图4[/align]（3）如需直接复制软件中的曲线图，点击Edit菜单下的Copy选项即可（图5）。[align=center][img=,144,117]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952337376_1353_3224499_3.jpg!w144x117.jpg[/img][/align][align=center]图5[/align]下面介绍FGP曲线的得到方法。用鼠标拖动GS曲线中的黑色竖线即可用来查看在不同时刻或温度下的红外光谱曲线（图6），此时可以看到图6中左上部分的红外光谱图在发生变化。在左上图中移动红色的竖线至感兴趣的官能团所对应的波数位置，例如CO[sub]2[/sub]的特征吸收位置2357cm[sup]-1[/sup]处，右击鼠标会弹出如图7所示的菜单，选中View Profile选项，即可得到该波数随时间或温度变化的FGP曲线，如图8中绿色曲线所示。该曲线在软件中的导出方法同EGP曲线。[align=center][img=,560,295]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952410486_6248_3224499_3.jpg!w560x295.jpg[/img][/align][align=center]图6[/align][align=center][img=,164,152]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952494022_9207_3224499_3.jpg!w164x152.jpg[/img][/align][align=center]图7[/align][align=center][img=,560,261]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170952589752_4178_3224499_3.jpg!w560x261.jpg[/img][/align][align=center]图8[/align]如需要导出不同时刻/温度下的实时红外光谱图，如图8中1499.1s处的红外光谱图，则需选中左上图下方的文件（图9），然后按照EGP曲线的导出方法即可得到该时刻的红外光谱图。[align=center][img=,560,274]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170953092976_4639_3224499_3.jpg!w560x274.jpg[/img][/align][align=center]图9[/align]另外，如需在软件中将EGP和FGP曲线的横坐标由时间转化为温度，则可以选中Process菜单中的Convert Time Axis选项（图10），弹出如图11所示的对话框。下方的选项为直接导入热重实验的文件的时间-温度曲线，上方的选项为定义两个指定的温度所对应的时间，由此会生成一个时间-温度的函数曲线，也可以在Excel文件中进行转换。[align=center] [img=,244,327]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170953226375_1876_3224499_3.jpg!w244x327.jpg[/img][/align][align=center]图10[/align][align=center][img=,404,344]https://ng1.17img.cn/bbsfiles/images/2019/10/201910170953340117_8499_3224499_3.jpg!w404x344.jpg[/img][/align][align=center]图11[/align]

[align=center][b]热分析/红外光谱联用数据分析时三维图的作图方法 [/b][/align][align=center]丁延伟[/align][align=center]中国科学技术大学理化科学实验中心[/align][align=center]（安徽省合肥市金寨路96号）[/align]使用热分析/红外光谱联用技术，在得到样品在温度变化过程中质量的变化信息的同时还可以得到逸出气体随温度变化的信息。通过红外光谱测定，可以判定未知样品中存在哪些有机官能团，这为最终确定未知物的化学结构奠定了基础。在通过热分析/红外光谱联用技术得到数据后，为了便于分析通常会在作图软件如Origin软件中进行重新作图，许多人在得到由仪器的分析软件导出的数据后，不知如何做出如图1所示的三维图。当然通过热分析/红外光谱联用技术得到的信息不仅仅局限于三维图，还可以得到总吸光值、不同官能团的吸光值等随温度和/或时间的变化关系曲线，据此可以分析样品随温度变化过程中结构的变化过程，这些内容将在以后的文章中进行介绍。在本文中将介绍如何通过PE公司的热重/红外光谱联用仪得到的导出的ASCII文件（通常为txt或者excel格式）在Origin软件中得到三维图的方法。[align=center] [img=,634,359]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051615551824_6299_3237657_3.jpg!w634x359.jpg[/img] 图1仪器分析软件得到的三维FTIR图[/align]由于热分析时间较长，对于逸出气体在使用红外光谱仪进行实时监测时，通常采用较快的分析时间（8波数分辨率时，DTGS检测器，平均1s得到一张谱图）。因此，对于一个普通的热分析实验来说（例如，20℃/min加热速率，室温-800摄氏度的试验温度范围），一次实验平均得到2000张以上的红外光谱图。由此得到的excel文件通常为20M左右，在excel软件中直接打开时往往无法显示所有的谱图数据。因此，通常在Origin软件中直接导入该文件，打开方式如图2所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051615109434_2866_3237657_3.png!w690x388.jpg[/img][/align][align=center]图2 Origin软件中导入Excel文件的方法示意图[/align]点击图中选项后会弹出以下窗口（图3）：[align=center][img=,690,387]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051616122214_3758_3237657_3.png!w690x387.jpg[/img][/align][align=center]图3[/align]点击open即可将文件在origin软件中打开，如图4所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051616271004_2412_3237657_3.png!w690x388.jpg[/img][/align][align=center]图4[/align]图4中，黄色区域为实验时的时间列，由于软件设置的问题显示在了标题栏，不过这不会影响后期的分析，第一列为波数。点击左上方空白栏全选图中所有数据（注意图中黑色箭头），如图5所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051616538527_4468_3237657_3.png!w690x388.jpg[/img][/align][align=center]图5[/align]将选中的数据复制在新建的一个窗口中，并删除第一列波数数据，如图6所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051617164670_4524_3237657_3.png!w690x388.jpg[/img][/align][align=center]图6[/align]然后按照图7所示将表中所有格式转换为矩阵形式。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051617346634_4957_3237657_3.png!w690x388.jpg[/img][/align][align=center]图7[/align]点击open dialog选项，弹出以下对话窗口（图8）。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051617510244_1980_3237657_3.png!w690x388.jpg[/img][/align][align=center]图8[/align]继续点击ok选项，在新的窗口中生成了矩阵形式的数据（图9）。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051618066250_8843_3237657_3.png!w690x388.jpg[/img][/align][align=center]图9[/align]按照图10的方法设定X和Y坐标的信息，例如X轴为波数，Y轴为温度。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051618226084_1618_3237657_3.png!w690x388.jpg[/img][/align][align=center]图10[/align]在弹出的窗口中输入X轴和Y轴的起止范围，如图11所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051618442934_5220_3237657_3.png!w690x388.jpg[/img][/align][align=center]图11[/align]点击ok，接下来就可以作图了。点击图下方的图标，例如以图12的方式生成三维图。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051619074444_3843_3237657_3.png!w690x388.jpg[/img][/align][align=center]图12[/align]点击选项后，生成以下三维图（图13）。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051620062364_7396_3237657_3.png!w690x388.jpg[/img][/align][align=center]图13[/align]可以点击图中曲线改变曲线的颜色，点击XYZ轴的图注信息编辑不同坐标的信息，编辑后如图14所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051620315994_1435_3237657_3.png!w690x388.jpg[/img][/align][align=center]图14[/align]还可以点击图中三维图的显示方式来达到更好的效果，如图15所示。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051620499033_7940_3237657_3.png!w690x388.jpg[/img][/align][align=center]图15[/align]顺时针旋转后，变为图16。[align=center][img=,690,388]https://ng1.17img.cn/bbsfiles/images/2018/09/201809051621037004_7930_3237657_3.png!w690x388.jpg[/img][/align][align=center]图16[/align]还有其他的形式，可以根据具体需要来进行调整，不再详述。

在繁忙的检测工作中，数据分析是我们每位仪器分析人员在做检测时必做的事，如同家常便饭，下面以同一样品溶液中测铅，同一次检测，不同波长其数据相差大时的分析过程为例，望与各位版友交流和讨论. 铅的分析谱线选择220nm和283nm,数据结果分别为:Pb(220):280mg/kg,Pb(283)308mg/kg,两波长的数据相差接近10%，在标准文件或内部作业指导书没有明确规定釆用那条谱线的数据为结果时，需要仪器分析人员去分析并做出判断。1.从谱图上看，两个波长都没有受到干拢。见下图：http://ng1.17img.cn/bbsfiles/images/2013/04/201304011118_433179_1827064_3.jpghttp://ng1.17img.cn/bbsfiles/images/2013/04/201304011118_433180_1827064_3.jpg2.查看各波长曲线的相关线性,两条波长的线性都大于0.99999，见下图：http://ng1.17img.cn/bbsfiles/images/2013/04/201304011120_433183_1827064_3.jpghttp://ng1.17img.cn/bbsfiles/images/2013/04/201304011120_433184_1827064_3.jpg3.查看标样空白、样品空白和样品的基线强度。各波长标样空白、样品空白和样品的基线强度都很接近，见下图： 标样空白图：http://ng1.17img.cn/bbsfiles/images/2013/04/201304011121_433185_1827064_3.jpghttp://ng1.17img.cn/bbsfiles/images/2013/04/201304011121_433186_1827064_3.jpg样品空白图http://ng1.17img.cn/bbsfiles/images/2013/04/201304011121_433187_1827064_3.jpghttp://ng1.17img.cn/bbsfiles/images/2013/04/201304011122_433188_1827064_3.jpg样品的基线见1中的样品谱图。4.查看测样前后QC质控，QC质控理论浓度为0.5mg/L，见下图：测样前，各波长QC质控值：http://ng1.17img.cn/bbsfiles/images/2013/04/201304011123_433189_1827064_3.jpg测样后，各波长QC质控值：波长220nm和283nm的质控测样前后都有增大(Pb(220):+1.1%,Pb(283):2.8%)，但可接受。综上所述，我会以波长283nm结果出据结果。希望各位版友就数据分析中遇到的问题，发表提问或看法。我所想到与各位版友交流的有以下问题：a.同一元素，不同波长，同一次仪器读数中，数据相差大的原因？b.不同波长测的数据相差大时，版友们是怎么做出判断的？c.日常分析中，是用那些依据来进行数据分析呢？

说明：本文首发于“热分析与吸附”公众号（[url=http://mp.weixin.qq.com/s?__biz=MzI5MjUzMzQ0OA==&mid=2247484524&idx=1&sn=0471e13d67419c99e627f97c71c549a7&chksm=ec7ea1cbdb0928ddad4c1d2da2c2feee6c54624c0aa1a3df9f70aeead5cd180d5f8546c1d475&token=790168635&lang=zh_CN#rd]链接[/url]），欢迎关注公众号了解更多的热分析与吸附相关的内容。[b][font=华文楷体][size=14.0pt]在本部分内容中将介绍与红外光谱部分的数据处理与作图相关的内容。[/size][/font][font=华文楷体][size=14.0pt]1. [/size][/font][font=华文楷体][size=14.0pt]红外光谱数据的导入[/size][/font][font=华文楷体][size=14.0pt]点击红外光谱软件Timebase图标（图1），软件打开后的界面如图2所示。点击Flie菜单下的Open选项，按照数据的保存路径找到需要分析的质谱原始数据文件（图3）。点击图3中的Viewtemplate框中的不同选项（图4），可以在窗口中同时显示三维图、GS曲线以及不同时刻的红外光谱图等信息。例如，当选择GS Profile，Spectrum，Stack plot选项时，软件中的窗口变为图5所示。图5中，上半部分左侧窗口为不同时刻所对应的红外光谱图。右侧为实验范围内所得到的所有的红外光谱图的堆积图，该图的横坐标为波数，纵坐标为时间（可以换算为温度），Z轴为吸光值。图5的下半部分为GS曲线，可以反映实验过程中气体的逸出情况。在每一个窗口中点击鼠标右键，可以改变曲线的形状和颜色，进行进一步的处理（图6）。[/size][/font][/b][align=center][img=,206,247]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220717089227_1446_3132575_3.png!w206x247.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图1[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,425]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220717282510_137_3132575_3.png!w558x425.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图2[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,292]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220717385373_7173_3132575_3.png!w558x292.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图3[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,452,359]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220717520698_3790_3132575_3.png!w452x359.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图4[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,294]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220718040758_7245_3132575_3.png!w558x294.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图5[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,348,312]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220718193144_2314_3132575_3.png!w348x312.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图6[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]2. GS[/size][/font][font=华文楷体][size=14.0pt]曲线的数据分析与导出[/size][/font][font=华文楷体][size=14.0pt]GS[/size][/font][font=华文楷体][size=14.0pt]曲线的导出方法为：[/size][/font][font=华文楷体][size=14.0pt]（1）选中该文件（图7）；[/size][/font][/b][align=center][img=,348,446]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220718318097_8039_3132575_3.png!w348x446.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图7 选中GS曲线[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]（2）选中File菜单下的Save As选项，会弹出图9所示的窗口，点击保存即得到Excel格式的GS曲线的数据文件（图9）。导出的文件在Origin、Excel等文件中可以直接作图。[/size][/font][/b][align=center][img=,348,518]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220718454632_6002_3132575_3.png!w348x518.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图8 [/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,424]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220718599155_5299_3132575_3.png!w558x424.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图9[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]（3）如需直接复制软件中的曲线图，点击Edit菜单下的Copy选项即可（图10）。[/size][/font][/b][align=center][img=,354,291]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220719102311_511_3132575_3.png!w354x291.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图10[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]3.[/size][/font][font=华文楷体][size=14.0pt]官能团剖面图曲线（FGP曲线）的数据分析与导出[/size][/font][font=华文楷体][size=14.0pt]下面介绍FGP曲线的得到方法。[/size][/font][font=华文楷体][size=14.0pt]用鼠标拖动GS曲线中的黑色竖线即可用来查看在不同时刻或温度下的红外光谱曲线（图11），此时可以看到图11中左上部分的红外光谱图在发生变化。在左上图中移动红色的竖线至感兴趣的官能团所对应的波数位置，例如CO[sub]2[/sub]的特征吸收位置2357cm[sup]-1[/sup]处，右击鼠标会弹出如图12所示的菜单，选中View Profile选项，即可得到该波数随时间或温度变化的FGP曲线，如图13中绿色曲线所示。该曲线在软件中的导出方法同GS曲线。[/size][/font][/b][align=center][img=,510,291]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220719241444_1083_3132575_3.png!w510x291.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图11[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,270,370]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220719398126_9542_3132575_3.png!w270x370.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图12[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,611,370]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220719532513_5809_3132575_3.png!w611x370.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图13[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]4. [/size][/font][font=华文楷体][size=14.0pt]特定时刻/温度下的实时红外光谱图的数据分析与导出[/size][/font][font=华文楷体][size=14.0pt]如需要导出特定时刻/温度下的实时红外光谱图，如图14中2202.6s处的红外光谱图，则需选中左上图下方的文件（图14），然后按照GS曲线的导出方法即可得到该时刻的红外光谱图。[/size][/font][/b][align=center][img=,611,370]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220720058247_1677_3132575_3.png!w611x370.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图14[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]5. [/size][/font][font=华文楷体][size=14.0pt]实验范围内所有时刻/温度下的实时红外光谱图的数据分析与导出[/size][/font][font=华文楷体][size=14.0pt]数据分析时，需导出实验范围内所有时刻/温度下的实时红外光谱图时，点击Flie 菜单下的Save Time Resolved Data选项（图15），在弹出的窗口中将文件格式改为.csv格式（图16），即可得到所有的红外光谱图。图17为导出的所有红外光谱曲线的Excel文件，由图可见实验过程中总共1079张红外光谱图。图17中表格的第5行为实验中的不同时刻，A列为波数，其他列为不同时刻得到的红外光谱图的吸光值。在本系列第7部分内容中，将详细介绍在Origin 软件中如何利用这些数据得到三维图的方法。[/size][/font][/b][align=center][img=,295,469]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220720264578_3892_3132575_3.png!w295x469.jpg[/img][/align][align=center][b][font=华文楷体][size=14.0pt]图15[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,469]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220720415303_1476_3132575_3.png!w558x469.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图16[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,558,580]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220720545421_7937_3132575_3.png!w558x580.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图17[/size][/font][/b][/align][b][font=华文楷体][size=14.0pt]6. [/size][/font][font=华文楷体][size=14.0pt]横坐标由时间转化为温度的方法[/size][/font][font=华文楷体][size=14.0pt]如需在软件中将EGP和FGP曲线的横坐标由时间转化为温度，则可以选中Process菜单中的Convert Time Axis选项（图18），弹出如图19所示的对话框。下方的选项为直接导入热重实验的文件的时间-温度曲线，上方的选项为定义两个指定的温度所对应的时间，由此会生成一个时间-温度的函数曲线，也可以在Excel文件中进行转换。[/size][/font][/b][align=center][b][font=华文楷体][size=14.0pt] [img=,369,489]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220721088436_9061_3132575_3.png!w369x489.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图18[/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt][img=,403,345]https://ng1.17img.cn/bbsfiles/images/2020/06/202006220721252497_2222_3132575_3.png!w403x345.jpg[/img][/size][/font][/b][/align][align=center][b][font=华文楷体][size=14.0pt]图19[/size][/font][/b][/align]

[size=18px][b] 学习了文献计量学后，我好像对大数据分析有更深的领悟了。林肯说过一句富有哲理的话：你可以在一个时期蒙骗所有人，也可以长时期蒙骗一些人，但你无法在长时期蒙骗所有人。[/b][/size][size=18px][b] 我引申为对大数据分析的理解。你可以一时蒙骗过大数据，但只要你长时间这么做，肯定蒙骗不了大数据。[/b][/size][b][size=18px] [/size][/b]

求助数据分析，钴酸锂电池爆炸产物，谢谢，weiyongxing999@tom.com 请高手指教[~187404~][~187405~][img]http://ng1.17img.cn/bbsfiles/images/2009/12/200912020937_187406_1635498_3.jpg[/img][img]http://ng1.17img.cn/bbsfiles/images/2009/12/200912020938_187407_1635498_3.jpg[/img]