Misc_Routines.py
60 lines
| 1.6 KiB
| text/x-python
|
PythonLexer
| r1738 | """ | ||
| 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 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 | |||
