{"id":10795,"date":"2016-04-06T08:54:00","date_gmt":"2016-04-06T08:54:00","guid":{"rendered":"http:\/\/melotopia.net\/b\/?p=10795"},"modified":"2016-04-06T08:54:00","modified_gmt":"2016-04-06T08:54:00","slug":"%ea%b4%91%ec%84%ac%ec%9c%a0-%ea%b0%80%ea%b3%b5-%ed%9b%84-%eb%8d%b0%ec%9d%b4%ed%84%b0%eb%a5%bc-%ec%9e%90%eb%8f%99%ec%9c%bc%eb%a1%9c-%eb%b6%84%ec%84%9d%ed%95%98%eb%8a%94-%ed%94%84%eb%a1%9c%ea%b7%b8","status":"publish","type":"post","link":"http:\/\/melotopia.net\/b\/?p=10795","title":{"rendered":"\uad11\uc12c\uc720 \uac00\uacf5 \ud6c4 \ub370\uc774\ud130\ub97c \uc790\ub3d9\uc73c\ub85c \ubd84\uc11d\ud558\ub294 \ud504\ub85c\uadf8\ub7a8"},"content":{"rendered":"
\n
\n #!\/usr\/bin\/env python
\n
\n # -*- coding: CP949 -*-
\n
\n #
\n
\n # generated by wxGlade 0.7.2 on Tue Apr 05 01:49:01 2016
\n
\n #
\n
\n # This is a converting program from thorlabs powermeter datafile taken from thorlab official data acquisition program to radius vs power after fiber taper process.
\n
\n # 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.
\n
\n # snowall@kaist.ac.kr \/ snowall@gmail.com \/ Keehwan Nam, Physics department, KAIST. 2016.
\n
\n #
\n
\n import wx<\/p>\n

# begin wxGlade: dependencies
\n
\n import gettext
\n
\n # end wxGlade<\/p>\n

# begin wxGlade: extracode
\n
\n # end wxGlade
\n
\n import numpy as np
\n
\n import os
\n
\n import re
\n
\n import sys<\/p>\n

import struct<\/p>\n

f= open(‘sample.png’, ‘rb’)<\/p>\n

w, h = struct.unpack(‘>LL’, f.read()[16:24])
\n
\n displaysizeX = int(w)
\n
\n displaysizeY = int(h)
\n
\n f.close()
\n
\n reload(sys)
\n
\n sys.setdefaultencoding(‘cp949’)
\n
\n #~ displaysizeX=1280
\n
\n #~ displaysizeY=1024
\n
\n x,y=0,0
\n
\n a=0
\n
\n b=0
\n
\n c=0
\n
\n sec=0
\n
\n minute=0
\n
\n hours=0
\n
\n msec=0
\n
\n class MyDialog(wx.Dialog):
\n
\n x00=0
\n
\n y00=0
\n
\n def __init__(self, parent, id, title):
\n
\n wx.Dialog.__init__(self, parent, id, title, size=(displaysizeX,displaysizeY))
\n
\n self.canvas = wx.StaticBitmap(self, -1, wx.Bitmap(“sample.png”, wx.BITMAP_TYPE_ANY))
\n
\n self.canvas.Bind(wx.EVT_LEFT_DOWN, self.positionnow)<\/p>\n

def positionnow(self, event):
\n
\n self.x00=wx.MouseEvent.GetX(event)
\n
\n self.y00=wx.MouseEvent.GetY(event)
\n
\n print((wx.MouseEvent.GetX(event),wx.MouseEvent.GetY(event)))<\/p>\n

def drawTick(self, pl):
\n
\n self.mydc = wx.MemoryDC()
\n
\n self.canvas2 = wx.Bitmap(“sample.png”, wx.BITMAP_TYPE_ANY)
\n
\n self.mydc.SelectObject(self.canvas2)
\n
\n self.mypen=wx.Pen(wx.Colour(255,0,0),5)
\n
\n self.mypen.SetWidth(1)
\n
\n self.mydc.SetPen(self.mypen)
\n
\n self.mydc.DrawLines(pl)
\n
\n self.mydc.SelectObject(wx.NullBitmap)
\n
\n self.canvas.SetBitmap(self.canvas2)
\n
\n self.mydc = None
\n
\n self.canvas = wx.StaticBitmap(self, -1, self.canvas2)
\n
\n return<\/p>\n

def graphroutine(self, myf):
\n
\n self.convertedData=[]
\n
\n self.pointlist=[]
\n
\n src=open(myf, “r”)
\n
\n src.readline()
\n
\n d=src.readline()
\n
\n a=d.split(” “)
\n
\n self.mini=np.inf
\n
\n self.maxi=0.
\n
\n while True:
\n
\n d=src.readline()
\n
\n if d==””:
\n
\n break
\n
\n a=d.split(” “)
\n
\n msec=float(((a[1]).split(“.”))[1])*0.001
\n
\n sec=float(a[2].split(“:”)[2].split(“\\t”)[0])
\n
\n minute=float((a[2]).split(“:”)[1])*60.
\n
\n hours=float((a[2]).split(“:”)[0])*3600.
\n
\n self.time=msec+sec+minute+hours
\n
\n self.convertedData+=[[self.time,float(a[2].split(“:”)[2].split(“\\t”)[1])]]
\n
\n if self.convertedData[-1][1]\n self.mini=self.convertedData[-1][1]
\n
\n if self.convertedData[-1][1]>self.maxi:
\n
\n self.maxi=self.convertedData[-1][1]
\n
\n self.time0=self.convertedData[0][0]
\n
\n self.time1=self.convertedData[-1][0]
\n
\n self.timescale=displaysizeX\/(self.time1-self.time0)
\n
\n self.vertscale=displaysizeY\/(self.maxi-self.mini)
\n
\n for cd in self.convertedData:
\n
\n self.pointlist+=[[np.ceil((cd[0]-self.time0)*self.timescale),displaysizeY-np.ceil((cd[1]-self.mini)*self.vertscale)]]
\n
\n self.drawTick(self.pointlist)
\n
\n src.flush()
\n
\n src.close()
\n
\n return<\/p>\n

def findV(self,no):
\n
\n logfile=open(“log.csv”,”r”)
\n
\n try:
\n
\n while True:
\n
\n d=logfile.readline()
\n
\n if d==””:
\n
\n break
\n
\n a=d.split(“,”)
\n
\n if a[0]==str(no):
\n
\n logfile.close()
\n
\n return a
\n
\n except:
\n
\n logfile.close()
\n
\n return “0”<\/p>\n

def inverse_ratio(self,logdata):
\n
\n return np.sqrt(((2.*float(logdata[15]))-float(logdata[14]))\/((2.*float(logdata[15]))+float(logdata[14])))<\/p>\n

def converting(self,f):
\n
\n self.time0=(self.x00\/self.timescale)+self.time0
\n
\n tar=open(f[:-4]+”_processed.txt”, “w”)
\n
\n errorlist=open(“errorlist.txt”,”a”)
\n
\n plot=open(f[:-4]+”.gpl”, “w”)
\n
\n shotnumber=re.search(“\\d+”,f.split(“\\\\”)[-1]).group()
\n
\n expcondition=self.findV(shotnumber)
\n
\n shotlog=self.inverse_ratio(expcondition)
\n
\n try:
\n
\n for cd in self.convertedData:
\n
\n 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”)
\n
\n plot.writelines(“set logscale x\\nplot \\””+(f.split(“\\\\”)[-1])[:-4]+”_processed.txt\\” using 3:4 w l “)
\n
\n except:
\n
\n errorlist.writelines(f+”\\n”)
\n
\n print “error”
\n
\n finally:
\n
\n tar.flush()
\n
\n tar.close()
\n
\n plot.flush()
\n
\n plot.close()<\/p>\n

class MainPlotFrame(wx.Frame):
\n
\n def __init__(self, *args, **kwds):
\n
\n # begin wxGlade: MainPlotFrame.__init__
\n
\n kwds[“style”] = wx.DEFAULT_FRAME_STYLE
\n
\n wx.Frame.__init__(self, *args, **kwds)
\n
\n self.dig = wx.FileDialog(None, message=”Choose data files”, style = wx.FD_MULTIPLE)
\n
\n self.dia = MyDialog(self, -1, ‘Close after click!’)
\n
\n self.__set_properties()
\n
\n self.__do_layout()
\n
\n if self.dig.ShowModal() == wx.ID_OK:
\n
\n for fn in self.dig.GetPaths():
\n
\n self.dia.graphroutine(fn)
\n
\n self.dia.ShowModal()
\n
\n self.dia.converting(fn)
\n
\n self.Close()
\n
\n # end wxGlade<\/p>\n

def __set_properties(self):
\n
\n # begin wxGlade: MainPlotFrame.__set_properties
\n
\n self.SetTitle((“frame_1”))
\n
\n self.SetSize((displaysizeX, displaysizeY))
\n
\n # end wxGlade<\/p>\n

def __do_layout(self):
\n
\n # begin wxGlade: MainPlotFrame.__do_layout
\n
\n sizer_1 = wx.BoxSizer(wx.VERTICAL)
\n
\n self.SetSizer(sizer_1)
\n
\n self.Layout()
\n
\n # end wxGlade<\/p>\n

# end of class MainPlotFrame
\n
\n if __name__ == “__main__”:
\n
\n app = wx.App(0)
\n
\n #~ wx.InitAllImageHandlers()
\n
\n frame_main = MainPlotFrame(None, -1, “”)
\n
\n frame_main.Show()
\n
\n app.MainLoop()
\n <\/p>\n

\n\n<\/p>\n

\n
\n
\n\"\uc2e0\uace0\"
\n<\/a>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"

#!\/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 […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_crdt_document":"","_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[12],"tags":[],"class_list":["post-10795","post","type-post","status-publish","format-standard","hentry","category-12"],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/p8o6gA-2O7","jetpack-related-posts":[],"jetpack_likes_enabled":true,"_links":{"self":[{"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=\/wp\/v2\/posts\/10795","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=10795"}],"version-history":[{"count":0,"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=\/wp\/v2\/posts\/10795\/revisions"}],"wp:attachment":[{"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=10795"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=10795"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/melotopia.net\/b\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=10795"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}