|
|
"""
|
|
|
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
|
|
|
|
|
|
|
|
|
|
