Misc_Routines.py
81 lines
| 2.1 KiB
| text/x-python
|
PythonLexer
Fiorella Quino
|
r178 | """ | ||
| The module MISC_ROUTINES gathers classes and functions which are useful for daily processing. As an | ||||
| example we have conversion factor or universal constants. | ||||
| MODULES CALLED: | ||||
| NUMPY, SYS | ||||
| MODIFICATION HISTORY: | ||||
| Created by Ing. Freddy Galindo (frederickgalindo@gmail.com). ROJ, 21 October 2009. | ||||
| """ | ||||
| import numpy | ||||
| import sys | ||||
| class CoFactors(): | ||||
| """ | ||||
| CoFactor class used to call pre-defined conversion factor (e.g. degree to radian). The cu- | ||||
| The current available factor are: | ||||
| d2r = degree to radian. | ||||
| s2r = seconds to radian?, degree to arcsecond.? | ||||
| h2r = hour to radian. | ||||
| h2d = hour to degree | ||||
| """ | ||||
| d2r = numpy.pi/180. | ||||
| s2r = numpy.pi/(180.*3600.) | ||||
| h2r = numpy.pi/12. | ||||
| h2d = 15. | ||||
| class Redirect: | ||||
| def __init__(self,stdout): | ||||
| self.stdout = stdout | ||||
| def write(self,message): | ||||
| self.stdout.insertPlainText(message) | ||||
| class WidgetPrint: | ||||
| """ | ||||
| WidgetPrint class allows to define the standard output. | ||||
| """ | ||||
| def __init__(self,textid=None): | ||||
| self.__stdout = sys.stdout | ||||
| self.textid = textid | ||||
| self.wPrint() | ||||
| def wPrint(self): | ||||
| if self.textid == None: sys.stdout = self.__stdout | ||||
| if self.textid != None: sys.stdout = Redirect(self.textid) | ||||
| print ("") | ||||
| class Vector: | ||||
| """ | ||||
| direction = 0 Polar to rectangular; direction=1 rectangular to polar | ||||
| """ | ||||
| def __init__(self,vect,direction=0): | ||||
| nsize = numpy.size(vect) | ||||
| if nsize <= 3: | ||||
| vect = vect.reshape(1,nsize) | ||||
| self.vect = vect | ||||
| self.dirc = direction | ||||
| def Polar2Rect(self): | ||||
| if self.dirc == 0: | ||||
| jvect = self.vect*numpy.pi/180. | ||||
| mmx = numpy.cos(jvect[:,1])*numpy.sin(jvect[:,0]) | ||||
| mmy = numpy.cos(jvect[:,1])*numpy.cos(jvect[:,0]) | ||||
| mmz = numpy.sin(jvect[:,1]) | ||||
| mm = numpy.array([mmx,mmy,mmz]).transpose() | ||||
| elif self.dirc == 1: | ||||
| mm = [numpy.arctan2(self.vect[:,0],self.vect[:,1]),numpy.arcsin(self.vect[:,2])] | ||||
| mm = numpy.array(mm)*180./numpy.pi | ||||
| return mm | ||||