def test_smap(self):
target = 'SMAP'
et0 = spice.utc2et( '2016-06-01T12:00:00' )
dpr = spice.dpr()
a,b,c = [spice.gdpool('BODY399_RADII',i,1)[1] for i in range(3)]
ods = []
for deltatime in [0.1 * i for i in range(1080)]:
et = et0 + deltatime
stSmap,lsSmap = spice.spkezr( target, et, 'IAU_EARTH', 'NONE', 'EARTH' )
posn, veloc = stSmap[:3], stSmap[3:]
stSun,lsSun = spice.spkezr( 'SUN', et0, 'IAU_EARTH', 'LT', 'EARTH' )
mtx = spice.pxform( 'SMAP_REFLECTOR', 'IAU_EARTH', et)
boreEbf = spice.mxv( mtx, [1.0,0,0] )
point = spice.surfpt( posn, boreEbf, a, b, c)
rsurfpt,lon,lat = spice.reclat( point )
utc = spice.et2utc( et, 'ISOC', 3 )
ods += [ OD(deltatime=deltatime,posn=posn,veloc=veloc,boreEbf=boreEbf
,utc=utc,point=point,rsurfpt=rsurfpt
,rsmap=spice.vnorm(posn),lat=lat,lon=lon
,raynge=spice.vnorm(spice.vsub(point,posn))
,sunsep=spice.vsep( spice.ucrss(posn,veloc), stSun[:3] )
)
]
try:
### Moved matplotlib import to here so test runs to here at least
from matplotlib import pyplot as plt
plt.figure(1)
keys = 'lat lon raynge'.split()
secs = [od['deltatime'] for od in ods]
for idx in range(len(keys)):
scal = 1.0 if keys[idx] in 'rsurfpt rsmap raynge sunsep rp ecc t0 mu a P eccmean amean Pmean'.split() else dpr
ordinate = [od[keys[idx]]*scal for od in ods]
plt.subplot( 221+idx )
plt.plot( secs, ordinate )
plt.plot( secs, ordinate, '.')
plt.title( keys[idx] )
plt.ylabel( '%s%s' % (keys[idx],'' if scal==1.0 else ', deg') )
if idx>1: plt.xlabel( 'T-T0, s' )
abscissa = [od['lon']*dpr for od in ods]
ordinate = [od['lat']*dpr for od in ods]
plt.subplot( 221+idx+1 )
plt.title( 'lon vs. lat' )
plt.plot( abscissa, ordinate )
plt.xlabel( 'lon, deg' )
plt.ylabel( 'lat, deg' )
plt.show()
except:
print( "Bypassed, or failed, matplotlib tests" )
def test_bodvXds(self):
periState = spice.conics( self.elts, self.t0)
self.assertAlmostEqual( spice.vnorm(periState[:3]), self.rpKM, places=11)
helioIson = MPC_HELIOCENTRIC_ISON()
for jdtt,pos,vel in helioIson.jpv:
et = spice.str2et( jdtt )
conicState = spice.conics( self.elts, et )
print( conicState )
print( pos+vel )
self.assertAlmostEqual( spice.vdist(conicState[:3], pos), 0.0, places=3)
def test_smap(self):
target = 'SMAP'
et0 = spice.utc2et('2016-06-01T12:00:00')
dpr = spice.dpr()
a, b, c = [spice.gdpool('BODY399_RADII', i, 1)[1] for i in range(3)]
ods = []
for deltatime in [0.1 * i for i in range(1080)]:
et = et0 + deltatime
stSmap, lsSmap = spice.spkezr(target, et, 'IAU_EARTH', 'NONE',
'EARTH')
posn, veloc = stSmap[:3], stSmap[3:]
stSun, lsSun = spice.spkezr('SUN', et0, 'IAU_EARTH', 'LT', 'EARTH')
mtx = spice.pxform('SMAP_REFLECTOR', 'IAU_EARTH', et)
boreEbf = spice.mxv(mtx, [1.0, 0, 0])
point = spice.surfpt(posn, boreEbf, a, b, c)
rsurfpt, lon, lat = spice.reclat(point)
utc = spice.et2utc(et, 'ISOC', 3)
ods += [
OD(deltatime=deltatime,
posn=posn,
veloc=veloc,
boreEbf=boreEbf,
utc=utc,
point=point,
rsurfpt=rsurfpt,
rsmap=spice.vnorm(posn),
lat=lat,
lon=lon,
raynge=spice.vnorm(spice.vsub(point, posn)),
sunsep=spice.vsep(spice.ucrss(posn, veloc), stSun[:3]))
]
try:
### Moved matplotlib import to here so test runs to here at least
from matplotlib import pyplot as plt
plt.figure(1)
keys = 'lat lon raynge'.split()
secs = [od['deltatime'] for od in ods]
for idx in range(len(keys)):
scal = 1.0 if keys[
idx] in 'rsurfpt rsmap raynge sunsep rp ecc t0 mu a P eccmean amean Pmean'.split(
) else dpr
ordinate = [od[keys[idx]] * scal for od in ods]
plt.subplot(221 + idx)
plt.plot(secs, ordinate)
plt.plot(secs, ordinate, '.')
plt.title(keys[idx])
plt.ylabel('%s%s' %
(keys[idx], '' if scal == 1.0 else ', deg'))
if idx > 1: plt.xlabel('T-T0, s')
abscissa = [od['lon'] * dpr for od in ods]
ordinate = [od['lat'] * dpr for od in ods]
plt.subplot(221 + idx + 1)
plt.title('lon vs. lat')
plt.plot(abscissa, ordinate)
plt.xlabel('lon, deg')
plt.ylabel('lat, deg')
plt.show()
except:
print("Bypassed, or failed, matplotlib tests")
