1,背景
本人做光学的。。。最近研究实验偏向通信
做WDM实验的时候发现能用实验室TCP/IP通信的光谱仪(OSA)都被用了,翻翻找找只在角落找到一台Agilent 68146B的OSA。。。就是下图这个货,
想要捞取光谱数据,你看到那个软盘接口了吗。。。已经2020年奔小康了,我去哪找个软盘去读取呢,就算有软盘我笔记本也没有软盘口了。。
软盘测试取光谱数据肯定是不可能了,而且在测试方案设计上,需要不停的扫描WDM auto sacn模式去取得信噪比OSNR和功率power,然后利用判据去反馈调节光网络其他参数。
好消息是OSA后面有个GPIB的接口,后面经过调研(百度+bing),发现GPIB还是能用的,然鹅我的电脑也没有GPIB的口子,找了半天找到一根满是蜘蛛网的GPIB-USB的转接线。望着这一堆破烂。。我
后面在百度上看到有个类似的文章,也是用python控制仪器,关键的库为pyvisa。由于我研一用过python去爬过小姐姐,那段时间身体消瘦的特别快咳咳,也用xlwings对光谱数据进行过处理,因此对这个比较感兴趣。
这里实名感谢一下 @Luyao Zou 这个大哥的文章,特别致谢!给了我很大的启发,也让我意识到用python做成这件事情的可能性~,爱你!原文我摘取一部分,用python你能干什么事情
我的科研工作大量依赖自己开发的 Python 程序和脚本。比如,用 Python 和 PyVisa 控制实验仪器,用 numpy 和 scipy 拟合光谱
看到成功先例之后,后面我自己理解了一下仪器remote control的逻辑,其实就是分三部分。
- 通信协议以及驱动,包括当前流行的TCP/IP,或者RS232或者COM口通信以及早期GPIB通信,用于提供管道进行信息交互
- 仪器本身的编程手册,即发送给仪器的命令(SCPI协议?),将仪器的命令通过通信协议(TDC/IP,GPIB等)发送给仪器
- 实验逻辑,这个需要自己去制定专业的实验方案,和python没什么关系了,用py进行控制
有了上述三步,后面的事情就简单了(并不)
首先是装驱动天坑。。。,我用来实验的笔记本win7 32的,大概花了两天才从那根GPIB-USB的adlink官网的角落下载到了32位的驱动,然后找了半天NI-VISA,找到了一个非常古老的NI-VISA504离线版本装上了。。。通信协议以及驱动部分花了大概三天就完成了
第二步是装pyvisa,这个相比之下是最easy的了,直接pip install pyvisa就搞定了。安装失败一般是pypi的网速问题,请用国内阿里的源或者清华的源即可。
第三步去bing上面找了一圈PDF,都是OSA的编程指令,花了三小时调试好了老爷机,下面直接上光谱仪的WDM_AutoScan部分的代码,定义了一个类方便在实验的中控py文件调用
# -*- coding: utf-8 -*-
#!/usr/bin/pythonimport pyvisa
import xlwings as xw
import os
import matplotlib
import visa
import time
import paramiko#author_jzk 2020-7-22
#class OSA,仪器的指令不大好记,获取WDM模式下的波长,功率以及OSNR值class OSA():"""docstring for OSA"""'''初始化OSA,连接'''def __init__(self):osa_GPIB='GPIB0::25::INSTR' #GPIB口self.rm=visa.ResourceManager() #open visaself.instrument=self.rm.open_resource(osa_GPIB) #connect agilent osaprint('----------------> the OSA is connected successful n')def close(self):self.instrument.close()print('the osa connect is closed......')def wdm_AutoScan(self):'''进入WDM扫描光谱模式'''self.instrument.write('inst:sel "WDM_Autoscan"') print('----------------> now the OSA is using sub-system ----WDM Spectrum n')print('wait the sub-systems end for 20s ')time.sleep(20)def wdm_sweep(self):'''WDm auto test operation'''self.instrument.write('INITiate:IMMediate')print('----------------> auto-scan is running now, wait for 15s n')time.sleep(14)def wdm_wav_get(self):'''get channel wavelength'''wve=[]try:print('----------------> wavelength is calculating---')str_wve=self.instrument.query('CALCulate:DATA:CWAVelengths?').split(',')for i in str_wve:wve.append(float(i))time.sleep(1)return wveexcept Exception as e:print('------->None data is calculated, please check the optical input')raise edef wdm_OSNR_get(self):'''get osnr '''osnr=[]try:print('----------------> OSNR is calculating---')str_osnr=self.instrument.query('CALCulate:DATA:CSNR?').split(',')for i in str_osnr:osnr.append(float(i))time.sleep(1)return osnrexcept Exception as e:print('------->None data is calculated, please check the optical input')raise e# 获取对应信号的功率def wdm_Power_get(self):power = []try:print('----------------> the power is calculated')str_power = self.instrument.query('CALCulate:DATA:CPOWers?').split(',')for i in str_power:power.append(float(i))time.sleep(0.5)return powerexcept Exception as e:print('------> None data is get, please check the optical input...')raise e#single sweep一遍'''def sweep(self):self.instrument.write('init:imm')print('take a sweep for osa,please wait for 2 seconds until it is ended')time.sleep(0.5)
总结,python很多人都是从很多看起来就很热的小姐姐爬虫程序开始,但是后面可以做什么还是看个人的想象力和使用工具方式,无论是前阵子用xlwings来画图或者干别的都行
大佬们 @Luyao Zou 那篇文章都是2017年的了,也是我爬很热的小姐姐最多的那阵子,三年过去了终于开始使用工具走上正路了,俺导师看到这个估计也会替我高兴吧hhh,后续的GUI还没需求,有需求准备再挖个坑写一篇hhh