[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Orekit Users] orekit and Jython

To: orekit-users@orekit.org
Subject: Re: [Orekit Users] orekit and Jython
From: beowulf zhang <beowulf.zhang@gmail.com>
Date: Tue, 24 May 2011 20:20:16 +0800
Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:date :message-id:subject:from:to:content-type:content-transfer-encoding; bh=c/TE8Q6tindeefmwV2Tpeuzq/AzO2xl6rN/k2AAAHAI=; b=AgAQwmQZxyFzUohyHzTEQMUwWIdqrjpPDHQ4SIM073Y87TjFwTZpFwfCM9nb6M/sjw hwAZF/FKfmPLWujRuF9TWhVWwrRII3MWU0Pondro8SEP12MFaTk6lqX6rdJ7+ZLqPbIO RhLB3bZmwtDVSZ17SIRvOBFvswu+WCi40EUF0=
Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=nNVg4ndIgO+grJ0c09I8Xk/HHQc8a76DQRkF+ARNJV4LsqfpDMLUbxGl062CyUg6RS C67NQqhcelSqzhqZbCJQ0APebltiAywdeBJr9O61W0yVxCnuj/mRVfjNq7+TbAFjbkZH 3g660PUd0gLO0wlks0Ht6hq1xM7GvVKjOab6I=
In-reply-to: <4DDB596A.8080102@c-s.fr>
References: <FECF5338-38C5-44B9-9206-182300A2CEAC@gmail.com> <BANLkTikMEk7jc-nvc-b0ebuvoaVBA0EKfw@mail.gmail.com> <20110519094630.38174ueiadnjewgm@messagerie.si.c-s.fr> <BAY164-w2791EF23B3D4036D9984DBD58E0@phx.gbl> <BANLkTikLw=vp_2AzTPV7ge4DByTmi8WY2Q@mail.gmail.com> <BAY164-w20B6D1F8956C019564FD5DD5710@phx.gbl> <BANLkTi=XMbPK=x1xpSnx2ReVxcNJmjrQZg@mail.gmail.com> <BAY164-w266308D39058BF7B7304E5D5710@phx.gbl> <4DDB596A.8080102@c-s.fr>

OK

2011/5/24, Luc Maisonobe <Luc.Maisonobe@c-s.fr>:
> Hi All,
>
> Petrus, Brian, Beowulf, as you have already done investigating about
> language interfaces for Orekit, would you like to set up some Wiki pages
> about it (mainly reporting the examples in this thread) ?
>
> If so, then I can set up the wiki permissions on the forge so you can
> create and modify pages by yourself.
>
> Does this sound reasonable to you ?
>
> Luc
>
> Le 20/05/2011 18:17, Brian H a écrit :
>> You should check "javaclasspath" to make sure that the commons-math
>> directorie appear at the top of the class path
>>
>>> Date: Sat, 21 May 2011 00:04:17 +0800
>>> From: beowulf.zhang@gmail.com
>>> To: orekit-users@orekit.org
>>> Subject: Re: [Orekit Users] orekit and Jython
>>>
>>> Thank you!
>>>
>>> I add line
>>>
>>> javaclasspath '$matlabroot/java/orekit/commons-math-2.1.jar'
>>> javaclasspath '$matlabroot/java/orekit/orekit-5.0.jar'
>>>
>>> in script head, so all is fine!
>>>
>>> 2011/5/21, Brian H <brianh4321@hotmail.com>:
>>> >
>>> > PS. Make sure you add the orekit jar file. This has to be once per
>> session,
>>> > with the command
>>> >
>>> >
>> javaclasspath(javaclasspath('-dynamic'),'/Users/bhodges/Documents/myMatlab/orekit/orekit-5.0.jar')
>>> >
>>> >
>>> >> Date: Fri, 20 May 2011 23:00:38 +0800
>>> >> From: beowulf.zhang@gmail.com
>>> >> To: orekit-users@orekit.org
>>> >> Subject: Re: [Orekit Users] orekit and Jython
>>> >>
>>> >> I'm working in Matlab 7.7, I want to test your code, I have seted the
>>> >> orekit's and common-math's classpath in classpath.txt, but I get error
>>> >> message
>>> >>
>>> >> ??? Undefined variable "TimeScalesFactory" or class
>>> >> "TimeScalesFactory.getUTC".
>>> >>
>>> >> Error in ==> orekit_test1 at 42
>>> >> utc = TimeScalesFactory.getUTC();
>>> >>
>>> >> whether you can help me or not?
>>> >> Thank you!
>>> >>
>>> >> 2011/5/19, Brian H <brianh4321@hotmail.com>:
>>> >> >
>>> >> > Concerning the matlab interface(?)
>>> >> >
>>> >> > While its somewhat crude, here is my early test script of the first
>>> >> > tutorial. If I remember correctly you have to set the orekit data
>> path
>>> >> > using
>>> >> > the System call below. There was also a problem with matlab's
>> imported
>>> >> > version of
>>> >> > the commons-math library being somewhat out of date. This was fixed
>>> >> > by downloading the current version and adding that directory at
>> the top
>>> >> > of
>>> >> > your classpath.txt located in the MATLAB directory. Anyway, here is
>>> >> > a
>>> >> > copy
>>> >> > of the original mfile. It was
>>> >> > my first attempt at the matlab-java interface, pardon it's
>> crudeness it
>>> >> > was
>>> >> > meant as
>>> >> > an exercise.
>>> >> >
>>> >> >
>>> >> > import java.lang.*
>>> >> > import java.lang.Enum
>>> >> > import java.text.*
>>> >> > import java.util.*
>>> >> >
>>> >> > import org.apache.commons.*
>>> >> > import org.apache.commons.math.*
>>> >> > import org.apache.commons.math.geometry.*
>>> >> > import org.apache.commons.math.util.*
>>> >> > import org.apache.commons.math.util.MathUtils.*
>>> >> > import org.orekit.bodies.*
>>> >> > import org.orekit.bodies.BodyShape.*
>>> >> > import org.orekit.bodies.GeodeticPoint.*
>>> >> > import org.orekit.bodies.OneAxisEllipsoid.*
>>> >> > import org.orekit.errors.OrekitException.*
>>> >> > import org.orekit.frames.*
>>> >> > import org.orekit.frames.Frame.*
>>> >> > import org.orekit.frames.FramesFactory.*
>>> >> > import org.orekit.frames.TopocentricFrame.*
>>> >> > import org.orekit.frames.LocalOrbitalFrame.*
>>> >> > import org.orekit.orbits.*
>>> >> > import org.orekit.orbits.CartesianOrbit.*
>>> >> > import org.orekit.orbits.Orbit.*
>>> >> > import org.orekit.propagation.*
>>> >> > import org.orekit.propagation.Propagator.*
>>> >> > import org.orekit.propagation.analytical.*
>>> >> > import org.orekit.propagation.analytical.KeplerianPropagator.*
>>> >> > import org.orekit.time.*
>>> >> > import org.orekit.time.AbsoluteDate.*
>>> >> > import org.orekit.time.TimeScale.*
>>> >> > import org.orekit.time.TimeScalesFactory.*
>>> >> > import org.orekit.utils.*
>>> >> > import org.orekit.utils.PVCoordinates.*
>>> >> > import org.orekit.data.*
>>> >> > import org.orekit.data.DataProvidersManager.*
>>> >> > import LOFTypeWrapper
>>> >> > System.setProperty('orekit.data.path','orekit-data')
>>> >> >
>>> >> >
>>> >> > symbols = DecimalFormatSymbols(Locale.US);
>>> >> > d3 = DecimalFormat('0.000', symbols);
>>> >> > utc = TimeScalesFactory.getUTC();
>>> >> > initialDate = AbsoluteDate(2008, 10, 01, 0, 0, 00.000, utc);
>>> >> > mu = 3.986004415e+14; % gravitation coefficient
>>> >> > inertialFrame = FramesFactory.getEME2000();
>>> >> > posisat = Vector3D(-6142438.668, 3492467.560, -25767.25680);
>>> >> > velosat = Vector3D(505.8479685, 942.7809215, 7435.922231);
>>> >> > pvsat = PVCoordinates(posisat, velosat);
>>> >> > initialOrbit = CartesianOrbit(pvsat, inertialFrame, initialDate,
>>> >> > mu);
>>> >> > kepler = KeplerianPropagator(initialOrbit);
>>> >> > lof = LocalOrbitalFrame(inertialFrame, LOFTypeWrapper.getType(1),
>>> >> > kepler,
>>> >> > 'LOF');
>>> >> > ae = 6378137.0; % equatorial radius in meter
>>> >> > f = 1.0 / 298.257223563; % flattening
>>> >> > ITRF2005 = FramesFactory.getITRF2005(true); % terrestrial frame at
>>> >> > an
>>> >> > arbitrary date
>>> >> > earth = OneAxisEllipsoid(ae, f, ITRF2005);
>>> >> >
>>> >> > % Station
>>> >> > longitude = Math.toRadians(45.);
>>> >> > latitude = Math.toRadians(25.);
>>> >> > altitude = 0.;
>>> >> > station = GeodeticPoint(latitude, longitude, altitude);
>>> >> > staF = TopocentricFrame(earth, station, 'station1');
>>> >> > System.out.println(' time doppler (m/s)');
>>> >> > finalDate = AbsoluteDate(initialDate, 6000, utc);
>>> >> > extrapDate = initialDate;
>>> >> > while (extrapDate.compareTo(finalDate) <= 0)
>>> >> >
>>> >> > % We can simply get the position and velocity of station in LOF
>>> >> > frame at
>>> >> > any time
>>> >> > pv = staF.getTransformTo(lof,
>>> >> > extrapDate).transformPVCoordinates(PVCoordinates.ZERO);
>>> >> >
>>> >> > % And then calculate the doppler signal
>>> >> > doppler = Vector3D.dotProduct(pv.getPosition(), pv.getVelocity()) /
>>> >> > pv.getPosition().getNorm();
>>> >> > abc= d3.format(doppler);
>>> >> > def= String(' ');
>>> >> > ghi= extrapDate.toString();
>>> >> > jkl= ghi.concat(def.concat(abc));
>>> >> > System.out.println(jkl);
>>> >> >
>>> >> > extrapDate = AbsoluteDate(extrapDate, 600, utc);
>>> >> >
>>> >> > end
>>> >> >
>>> >> >
>>> >> >
>>> >> >> Date: Thu, 19 May 2011 09:46:30 +0200
>>> >> >> From: luc.maisonobe@c-s.fr
>>> >> >> To: orekit-users@orekit.org
>>> >> >> Subject: Re: [Orekit Users] orekit and Jython
>>> >> >>
>>> >> >> beowulf zhang <beowulf.zhang@gmail.com> a écrit :
>>> >> >>
>>> >> >> > Very good! Thank you...
>>> >> >> >
>>> >> >> > I konw the m language of Matlab support java programming now,
>> but I
>>> >> >> > was failure to test the orekit examples. I want to know whether
>>> >> >> > anyone
>>> >> >> > are interesting testing the orekit examples in Matlab or not.
>>> >> >>
>>> >> >> If you want to test these examples, I would be interested to know
>>> >> >> if
>>> >> >> they work.
>>> >> >> We could set up some wiki pages (yes, we know, the wiki is quite
>> empty
>>> >> >> for now) showing how to interface with several scripting languages.
>>> >> >>
>>> >> >> Luc
>>> >> >>
>>> >> >> >
>>> >> >> >
>>> >> >> > 2011/5/17, Petrus Hyvönen <petrus.hyvonen@gmail.com>:
>>> >> >> >> Hi all,
>>> >> >> >>
>>> >> >> >> I was playing around with the idea of using Python for mission
>>> >> >> >> design,
>>> >> >> >> and
>>> >> >> >> found Jython and the orekit library. As new to both Python
>> and Java,
>>> >> >> >> I
>>> >> >> >> had
>>> >> >> >> some initial troubles of getting things to work, primary the
>> access
>>> >> >> >> to
>>> >> >> >> orekit-data, what was current dir and javapaths etc. I am now
>> using
>>> >> >> >> Eclipse
>>> >> >> >> with pydev extension, and can add the orekit jars and
>> orekit-data to
>>> >> >> >> the
>>> >> >> >> eclipse project.
>>> >> >> >>
>>> >> >> >> As an exercise I took the liberty to translate a few of the
>> orekit
>>> >> >> >> examples
>>> >> >> >> from java to jython, and copy them below as it could be of
>> interest
>>> >> >> >> for
>>> >> >> >> others. I have been exploring scipy before (based on regular
>> python)
>>> >> >> >> and it
>>> >> >> >> contains very good plotting and array routines, which is not
>>> >> >> >> available
>>> >> >> >> in
>>> >> >> >> jython. So next step would be to find a suitable java
>> plotting and
>>> >> >> >> arrays
>>> >> >> >> library to be able to represent the results.
>>> >> >> >>
>>> >> >> >> Regards
>>> >> >> >> /Petrus
>>> >> >> >>
>>> >> >> >> First, a jython translation of the SlaveMode.java:
>>> >> >> >> __________________________________________
>>> >> >> >>
>>> >> >> >> # -*- coding: utf-8 -*-
>>> >> >> >> '''
>>> >> >> >> /* Copyright 2002-2010 CS Communication & Systèmes
>>> >> >> >> * Licensed to CS Communication & Systèmes (CS) under one or more
>>> >> >> >> * contributor license agreements. See the NOTICE file
>>> >> >> >> distributed
>>> >> >> >> with
>>> >> >> >> * this work for additional information regarding copyright
>>> >> >> >> ownership.
>>> >> >> >> * CS licenses this file to You under the Apache License, Version
>>> >> >> >> 2.0
>>> >> >> >> * (the "License"); you may not use this file except in
>>> >> >> >> compliance
>>> >> >> >> with
>>> >> >> >> * the License. You may obtain a copy of the License at
>>> >> >> >> *
>>> >> >> >> * http://www.apache.org/licenses/LICENSE-2.0
>>> >> >> >> *
>>> >> >> >> * Unless required by applicable law or agreed to in writing,
>>> >> >> >> software
>>> >> >> >> * distributed under the License is distributed on an "AS IS"
>> BASIS,
>>> >> >> >> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express
>>> >> >> >> or
>>> >> >> >> implied.
>>> >> >> >> * See the License for the specific language governing
>>> >> >> >> permissions
>>> >> >> >> and
>>> >> >> >> * limitations under the License.
>>> >> >> >> */
>>> >> >> >>
>>> >> >> >> Translated from SlaveMode.java to jython by Petrus Hyvönen
>>> >> >> >> 2011-05-04
>>> >> >> >> '''
>>> >> >> >>
>>> >> >> >> # orekit.jar,orekit-data.zip and commons-maths in CLASSPATH
>> through
>>> >> >> >> eclipse
>>> >> >> >> project
>>> >> >> >> import java, os
>>> >> >> >>
>>> >> >> >> from org.orekit.errors import OrekitException
>>> >> >> >> from org.orekit.frames import Frame
>>> >> >> >> from org.orekit.frames import FramesFactory
>>> >> >> >> from org.orekit.orbits import KeplerianOrbit
>>> >> >> >> from org.orekit.orbits import Orbit
>>> >> >> >> from org.orekit.propagation import SpacecraftState
>>> >> >> >> from org.orekit.propagation.analytical import
>>> >> >> >> KeplerianPropagator
>>> >> >> >> from org.orekit.data import DataProvidersManager
>>> >> >> >> from org.orekit.data import ZipJarCrawler
>>> >> >> >> from org.orekit.time import AbsoluteDate
>>> >> >> >> from org.orekit.time import TimeScalesFactory
>>> >> >> >>
>>> >> >> >> from math import radians
>>> >> >> >>
>>> >> >> >> # Configure Orekit. The file orekit-data.zip must be in
>> current dir
>>> >> >> >> DM = DataProvidersManager.getInstance()
>>> >> >> >> crawler=ZipJarCrawler("orekit-data.zip")
>>> >> >> >> DM.clearProviders()
>>> >> >> >> DM.addProvider(crawler)
>>> >> >> >>
>>> >> >> >> #Initial orbit parameters
>>> >> >> >> a = 24396159 # semi major axis in meters
>>> >> >> >> e = 0.72831215 # eccentricity
>>> >> >> >> i = radians(7.0)# inclination
>>> >> >> >> omega = radians(180) # perigee argument
>>> >> >> >> raan = radians(261) #right ascension of ascending node
>>> >> >> >> lM = 0.0 # mean anomaly
>>> >> >> >>
>>> >> >> >> #Inertial frame
>>> >> >> >> inertialFrame = FramesFactory.getEME2000()
>>> >> >> >>
>>> >> >> >> #Initial date in UTC time scale
>>> >> >> >> utc = TimeScalesFactory.getUTC();
>>> >> >> >> initialDate = AbsoluteDate(2004, 01, 01, 23, 30, 00.000, utc)
>>> >> >> >>
>>> >> >> >> #gravitation coefficient
>>> >> >> >> mu = 3.986004415e+14
>>> >> >> >>
>>> >> >> >> #Orbit construction as Keplerian
>>> >> >> >> initialOrbit = KeplerianOrbit(a, e, i, omega, raan, lM,
>>> >> >> >> KeplerianOrbit.MEAN_ANOMALY,
>>> >> >> >> inertialFrame, initialDate, mu)
>>> >> >> >>
>>> >> >> >> #Simple extrapolation with Keplerian motion
>>> >> >> >> kepler = KeplerianPropagator(initialOrbit)
>>> >> >> >>
>>> >> >> >> #Set the propagator to slave mode (could be omitted as it is the
>>> >> >> >> default
>>> >> >> >> mode)
>>> >> >> >> kepler.setSlaveMode()
>>> >> >> >>
>>> >> >> >> #Overall duration in seconds for extrapolation
>>> >> >> >> duration = 90*60.0
>>> >> >> >>
>>> >> >> >> #Stop date
>>> >> >> >> finalDate = AbsoluteDate(initialDate, duration, utc)
>>> >> >> >>
>>> >> >> >> #Step duration in seconds
>>> >> >> >> stepT = 30.0
>>> >> >> >>
>>> >> >> >> #Extrapolation loop
>>> >> >> >> cpt = 1
>>> >> >> >> extrapDate = initialDate
>>> >> >> >> while extrapDate.compareTo(finalDate) <= 0:
>>> >> >> >> currentState = kepler.propagate(extrapDate)
>>> >> >> >> print "step %d: time %s %s" % (cpt, currentState.getDate(),
>>> >> >> >> currentState.getOrbit())
>>> >> >> >> extrapDate = AbsoluteDate(extrapDate, stepT, utc)
>>> >> >> >> cpt=cpt+1
>>> >> >> >>
>>> >> >> >>
>>> >> >> >> ____________________________
>>> >> >> >>
>>> >> >> >> An interesting example of jython inheritance of java class:
>>> >> >> >> _____________________________
>>> >> >> >>
>>> >> >> >>
>>> >> >> >> # -*- coding: utf-8 -*-
>>> >> >> >> '''
>>> >> >> >> /* Copyright 2002-2010 CS Communication & Syst?mes
>>> >> >> >> * Licensed to CS Communication & Syst?mes (CS) under one or more
>>> >> >> >> * contributor license agreements. See the NOTICE file
>>> >> >> >> distributed
>>> >> >> >> with
>>> >> >> >> * this work for additional information regarding copyright
>>> >> >> >> ownership.
>>> >> >> >> * CS licenses this file to You under the Apache License, Version
>>> >> >> >> 2.0
>>> >> >> >> * (the "License"); you may not use this file except in
>>> >> >> >> compliance
>>> >> >> >> with
>>> >> >> >> * the License. You may obtain a copy of the License at
>>> >> >> >> *
>>> >> >> >> * http://www.apache.org/licenses/LICENSE-2.0
>>> >> >> >> *
>>> >> >> >> * Unless required by applicable law or agreed to in writing,
>>> >> >> >> software
>>> >> >> >> * distributed under the License is distributed on an "AS IS"
>> BASIS,
>>> >> >> >> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express
>>> >> >> >> or
>>> >> >> >> implied.
>>> >> >> >> * See the License for the specific language governing
>>> >> >> >> permissions
>>> >> >> >> and
>>> >> >> >> * limitations under the License.
>>> >> >> >> */
>>> >> >> >>
>>> >> >> >> Translated from SlaveMode.java to jython by Petrus Hyv?nen
>>> >> >> >> 2011-05-04
>>> >> >> >> '''
>>> >> >> >> # orekit and common maths in CLASSPATH through eclipse project
>>> >> >> >> import java, os
>>> >> >> >>
>>> >> >> >> from org.orekit.data import DataProvidersManager
>>> >> >> >> from org.orekit.data import ZipJarCrawler
>>> >> >> >> from org.apache.commons.math.geometry import Vector3D
>>> >> >> >> from org.orekit.bodies import BodyShape
>>> >> >> >> from org.orekit.bodies import GeodeticPoint
>>> >> >> >> from org.orekit.bodies import OneAxisEllipsoid
>>> >> >> >> from org.orekit.errors import OrekitException;
>>> >> >> >> from org.orekit.frames import Frame
>>> >> >> >> from org.orekit.frames import FramesFactory
>>> >> >> >> from org.orekit.frames import TopocentricFrame
>>> >> >> >> from org.orekit.orbits import KeplerianOrbit
>>> >> >> >> from org.orekit.orbits import Orbit
>>> >> >> >> from org.orekit.propagation import Propagator
>>> >> >> >> from org.orekit.propagation import SpacecraftState
>>> >> >> >> from org.orekit.propagation.analytical import
>>> >> >> >> KeplerianPropagator
>>> >> >> >> from org.orekit.propagation.events import ElevationDetector
>>> >> >> >> from org.orekit.propagation.events import EventDetector
>>> >> >> >> from org.orekit.time import AbsoluteDate
>>> >> >> >> from org.orekit.time import TimeScalesFactory
>>> >> >> >> from org.orekit.utils import PVCoordinates
>>> >> >> >>
>>> >> >> >> from math import degrees, radians, pi
>>> >> >> >>
>>> >> >> >> # Configure Orekit
>>> >> >> >> DM = DataProvidersManager.getInstance()
>>> >> >> >> crawler=ZipJarCrawler("orekit-data.zip")
>>> >> >> >> DM.clearProviders()
>>> >> >> >> DM.addProvider(crawler)
>>> >> >> >>
>>> >> >> >> # Initial state definition: date, orbit
>>> >> >> >> initialDate = AbsoluteDate(2004, 01, 01, 23, 30, 00.000,
>>> >> >> >> TimeScalesFactory.getUTC())
>>> >> >> >> mu = 3.986004415e+14
>>> >> >> >> inertialFrame = FramesFactory.getEME2000() # inertial frame for
>>> >> >> >> orbit
>>> >> >> >> definition
>>> >> >> >> position = Vector3D(-6142438.668, 3492467.560, -25767.25680)
>>> >> >> >> velocity = Vector3D(505.8479685, 942.7809215, 7435.922231)
>>> >> >> >> pvCoordinates = PVCoordinates(position, velocity)
>>> >> >> >> initialOrbit = KeplerianOrbit(pvCoordinates, inertialFrame,
>>> >> >> >> initialDate, mu)
>>> >> >> >>
>>> >> >> >> # Propagator : consider a simple keplerian motion (could be more
>>> >> >> >> elaborate)
>>> >> >> >> kepler = KeplerianPropagator(initialOrbit)
>>> >> >> >>
>>> >> >> >> #Earth and frame
>>> >> >> >> ae = 6378137.0 # // equatorial radius in meter
>>> >> >> >> f = 1.0 / 298.257223563 #; // flattening
>>> >> >> >> ITRF2005 = FramesFactory.getITRF2005() #; // terrestrial
>> frame at an
>>> >> >> >> arbitrary date
>>> >> >> >> earth = OneAxisEllipsoid(ae, f, ITRF2005)
>>> >> >> >>
>>> >> >> >> # Station
>>> >> >> >> longitude = radians(45.0)
>>> >> >> >> latitude = radians(25.0)
>>> >> >> >> altitude = 0.0
>>> >> >> >> station1 = GeodeticPoint(latitude, longitude, altitude)
>>> >> >> >> sta1Frame = TopocentricFrame(earth, station1, "station1")
>>> >> >> >>
>>> >> >> >> # Event definition
>>> >> >> >> maxcheck = 1.0
>>> >> >> >> elevation = radians(5.0)
>>> >> >> >>
>>> >> >> >>
>>> >> >> >> class VisibilityDetector(ElevationDetector):
>>> >> >> >> # Class for handling the eventOccured java. Example of
>> subclassing
>>> >> >> >> # a java class in jython
>>> >> >> >> def __init__(self, maxCheck, elevation, topo):
>>> >> >> >> ElevationDetector.__init__(self,maxCheck, elevation, topo)
>>> >> >> >>
>>> >> >> >> def eventOccurred(self, s, increasing):
>>> >> >> >> if (increasing):
>>> >> >> >> print "Visibility on",
>>> >> >> >> self.topocentricFrame.getName(),"begins
>>> >> >> >> at" , s.getDate()
>>> >> >> >> else:
>>> >> >> >> print "Visibility on", self.topocentricFrame.getName(),
>>> >> >> >> "ends
>>> >> >> >> at" , s.getDate()
>>> >> >> >> return self.CONTINUE
>>> >> >> >>
>>> >> >> >> sta1Visi = VisibilityDetector(maxcheck, elevation, sta1Frame)
>>> >> >> >>
>>> >> >> >> #Add event to be detected
>>> >> >> >> kepler.addEventDetector(sta1Visi)
>>> >> >> >>
>>> >> >> >> #Propagate from the initial date to the first raising or for the
>>> >> >> >> fixed
>>> >> >> >> duration
>>> >> >> >> finalState = kepler.propagate(initialDate.shiftedBy(1500.0))
>>> >> >> >>
>>> >> >> >> print "Final state : " ,
>>> >> >> >> finalState.getDate().durationFrom(initialDate)
>>> >> >> >>
>>> >> >> >> _________________________________
>>> >> >> >>
>>> >> >> >>
>>> >> >> >>
>>> >> >> >
>>> >> >> >
>>> >> >>
>>> >> >>
>>> >> >>
>>> >> >> ----------------------------------------------------------------
>>> >> >> This message was sent using IMP, the Internet Messaging Program.
>>> >> >>
>>> >> >>
>>> >> >
>>> >
>
>

Follow-Ups:
- Re: [Orekit Users] orekit and Jython
  - From: Luc Maisonobe <Luc.Maisonobe@c-s.fr>

References:
- [Orekit Users] orekit and Jython
  - From: "Petrus Hyvönen" <petrus.hyvonen@gmail.com>
- Re: [Orekit Users] orekit and Jython
  - From: beowulf zhang <beowulf.zhang@gmail.com>
- Re: [Orekit Users] orekit and Jython
  - From: MAISONOBE Luc <luc.maisonobe@c-s.fr>
- RE: [Orekit Users] orekit and Jython
  - From: Brian H <brianh4321@hotmail.com>
- Re: [Orekit Users] orekit and Jython
  - From: beowulf zhang <beowulf.zhang@gmail.com>
- RE: [Orekit Users] orekit and Jython
  - From: Brian H <brianh4321@hotmail.com>
- Re: [Orekit Users] orekit and Jython
  - From: beowulf zhang <beowulf.zhang@gmail.com>
- RE: [Orekit Users] orekit and Jython
  - From: Brian H <brianh4321@hotmail.com>
- Re: [Orekit Users] orekit and Jython
  - From: Luc Maisonobe <Luc.Maisonobe@c-s.fr>

Prev by Date: Re: [Orekit Users] orekit and Jython
Next by Date: Re: [Orekit Users] orekit and Jython
Previous by thread: Re: [Orekit Users] orekit and Jython
Next by thread: Re: [Orekit Users] orekit and Jython
Index(es):
- Date
- Thread