광섬유 가공 후 데이터를 자동으로 분석하는 프로그램

#!/usr/bin/env python

# -*- coding: CP949 -*-

#

# generated by wxGlade 0.7.2 on Tue Apr 05 01:49:01 2016

#

# This is a converting program from thorlabs powermeter datafile taken from thorlab official data acquisition program to radius vs power after fiber taper process.

# IF you want to use this program for your own purpose except learning python and programming, you have to ask me licence policy and copyright.

# snowall@kaist.ac.kr / snowall@gmail.com / Keehwan Nam, Physics department, KAIST. 2016.

#

import wx

# begin wxGlade: dependencies

import gettext

# end wxGlade

# begin wxGlade: extracode

# end wxGlade

import numpy as np

import os

import re

import sys

import struct

f= open(‘sample.png’, ‘rb’)

w, h = struct.unpack(‘>LL’, f.read()[16:24])

displaysizeX = int(w)

displaysizeY = int(h)

f.close()

reload(sys)

sys.setdefaultencoding(‘cp949’)

#~ displaysizeX=1280

#~ displaysizeY=1024

x,y=0,0

a=0

b=0

c=0

sec=0

minute=0

hours=0

msec=0

class MyDialog(wx.Dialog):

x00=0

y00=0

def __init__(self, parent, id, title):

wx.Dialog.__init__(self, parent, id, title, size=(displaysizeX,displaysizeY))

self.canvas = wx.StaticBitmap(self, -1, wx.Bitmap(“sample.png”, wx.BITMAP_TYPE_ANY))

self.canvas.Bind(wx.EVT_LEFT_DOWN, self.positionnow)

def positionnow(self, event):

self.x00=wx.MouseEvent.GetX(event)

self.y00=wx.MouseEvent.GetY(event)

print((wx.MouseEvent.GetX(event),wx.MouseEvent.GetY(event)))

def drawTick(self, pl):

self.mydc = wx.MemoryDC()

self.canvas2 = wx.Bitmap(“sample.png”, wx.BITMAP_TYPE_ANY)

self.mydc.SelectObject(self.canvas2)

self.mypen=wx.Pen(wx.Colour(255,0,0),5)

self.mypen.SetWidth(1)

self.mydc.SetPen(self.mypen)

self.mydc.DrawLines(pl)

self.mydc.SelectObject(wx.NullBitmap)

self.canvas.SetBitmap(self.canvas2)

self.mydc = None

self.canvas = wx.StaticBitmap(self, -1, self.canvas2)

return

def graphroutine(self, myf):

self.convertedData=[]

self.pointlist=[]

src=open(myf, “r”)

src.readline()

d=src.readline()

a=d.split(” “)

self.mini=np.inf

self.maxi=0.

while True:

d=src.readline()

if d==””:

break

a=d.split(” “)

msec=float(((a[1]).split(“.”))[1])*0.001

sec=float(a[2].split(“:”)[2].split(“\t”)[0])

minute=float((a[2]).split(“:”)[1])*60.

hours=float((a[2]).split(“:”)[0])*3600.

self.time=msec+sec+minute+hours

self.convertedData+=[[self.time,float(a[2].split(“:”)[2].split(“\t”)[1])]]

if self.convertedData[-1][1] self.mini=self.convertedData[-1][1]

if self.convertedData[-1][1]>self.maxi:

self.maxi=self.convertedData[-1][1]

self.time0=self.convertedData[0][0]

self.time1=self.convertedData[-1][0]

self.timescale=displaysizeX/(self.time1-self.time0)

self.vertscale=displaysizeY/(self.maxi-self.mini)

for cd in self.convertedData:

self.pointlist+=[[np.ceil((cd[0]-self.time0)*self.timescale),displaysizeY-np.ceil((cd[1]-self.mini)*self.vertscale)]]

self.drawTick(self.pointlist)

src.flush()

src.close()

return

def findV(self,no):

logfile=open(“log.csv”,”r”)

try:

while True:

d=logfile.readline()

if d==””:

break

a=d.split(“,”)

if a[0]==str(no):

logfile.close()

return a

except:

logfile.close()

return “0”

def inverse_ratio(self,logdata):

return np.sqrt(((2.*float(logdata[15]))-float(logdata[14]))/((2.*float(logdata[15]))+float(logdata[14])))

def converting(self,f):

self.time0=(self.x00/self.timescale)+self.time0

tar=open(f[:-4]+”_processed.txt”, “w”)

errorlist=open(“errorlist.txt”,”a”)

plot=open(f[:-4]+”.gpl”, “w”)

shotnumber=re.search(“\d+”,f.split(“\\”)[-1]).group()

expcondition=self.findV(shotnumber)

shotlog=self.inverse_ratio(expcondition)

try:

for cd in self.convertedData:

tar.writelines(str(cd[0])+”\t”+str(cd[1])+”\t”+str(float(expcondition[9])*np.power(shotlog,(cd[0]-self.time0)/(float(expcondition[11])/float(expcondition[15]))))+”\t”+str(cd[1])+”\n”)

plot.writelines(“set logscale x\nplot \””+(f.split(“\\”)[-1])[:-4]+”_processed.txt\” using 3:4 w l “)

except:

errorlist.writelines(f+”\n”)

print “error”

finally:

tar.flush()

tar.close()

plot.flush()

plot.close()

class MainPlotFrame(wx.Frame):

def __init__(self, *args, **kwds):

# begin wxGlade: MainPlotFrame.__init__

kwds[“style”] = wx.DEFAULT_FRAME_STYLE

wx.Frame.__init__(self, *args, **kwds)

self.dig = wx.FileDialog(None, message=”Choose data files”, style = wx.FD_MULTIPLE)

self.dia = MyDialog(self, -1, ‘Close after click!’)

self.__set_properties()

self.__do_layout()

if self.dig.ShowModal() == wx.ID_OK:

for fn in self.dig.GetPaths():

self.dia.graphroutine(fn)

self.dia.ShowModal()

self.dia.converting(fn)

self.Close()

# end wxGlade

def __set_properties(self):

# begin wxGlade: MainPlotFrame.__set_properties

self.SetTitle((“frame_1”))

self.SetSize((displaysizeX, displaysizeY))

# end wxGlade

def __do_layout(self):

# begin wxGlade: MainPlotFrame.__do_layout

sizer_1 = wx.BoxSizer(wx.VERTICAL)

self.SetSizer(sizer_1)

self.Layout()

# end wxGlade

# end of class MainPlotFrame

if __name__ == “__main__”:

app = wx.App(0)

#~ wx.InitAllImageHandlers()

frame_main = MainPlotFrame(None, -1, “”)

frame_main.Show()

app.MainLoop()

댓글 남기기

이메일은 공개되지 않습니다.

This site uses Akismet to reduce spam. Learn how your comment data is processed.

Proudly powered by WordPress | Theme: Baskerville 2 by Anders Noren.

Up ↑

%d 블로거가 이것을 좋아합니다: