@@ -33,7 +33,6 @@ identical to the one presented here.
33
33
In order to query information for a specific Solar System body, a
34
34
``Horizons `` object has to be instantiated:
35
35
36
- .. code-block :: python
37
36
.. doctest-remote-data ::
38
37
39
38
>>> from astroquery.jplhorizons import Horizons
@@ -61,7 +60,7 @@ location is on Earth if it has not been specifically set. The following example
61
60
uses the coordinates of the `Statue of Liberty
62
61
<https://www.google.com/maps/place/Statue+of+Liberty+National+Monument/@40.6892534,-74.0466891,17z/data=!3m1!4b1!4m5!3m4!1s0x89c25090129c363d:0x40c6a5770d25022b!8m2!3d40.6892494!4d-74.0445004> `_
63
62
as the observer's location:
64
- .. code-block :: python
63
+
65
64
.. doctest-remote-data ::
66
65
67
66
>>> statue_of_liberty = {' lon' : - 74.0466891 ,
@@ -106,8 +105,6 @@ In the case of ambiguities in the name resolution, a list of matching objects
106
105
will be provided. In order to select an object from this list, provide the
107
106
respective id number or record number as ``id `` and use ``id_type=None ``:
108
107
109
-
110
- .. code-block :: python
111
108
.. doctest-remote-data ::
112
109
113
110
>>> from astroquery.jplhorizons import Horizons
@@ -145,7 +142,6 @@ for a given observer location (``location``) and epoch or range of epochs
145
142
(``epochs ``) in the form of an astropy table. The following example queries the
146
143
ephemerides of asteroid (1) Ceres for a range of dates as seen from Mauna Kea:
147
144
148
- .. code-block :: python
149
145
.. doctest-remote-data ::
150
146
151
147
>>> from astroquery.jplhorizons import Horizons
@@ -163,8 +159,13 @@ ephemerides of asteroid (1) Ceres for a range of dates as seen from Mauna Kea:
163
159
1 Ceres (A801 AA) 2010-Jan-31 00:00 2455227.5 ... 17.2067 247.2518 3.7289
164
160
1 Ceres (A801 AA) 2010-Feb-10 00:00 2455237.5 ... 18.5029 250.0576 3.4415
165
161
1 Ceres (A801 AA) 2010-Feb-20 00:00 2455247.5 ... 19.5814 252.7383 3.1451
166
- >>> # The following fields are available for each ephemerides query:
167
- >>> print (eph.columns)
162
+
163
+
164
+ The following fields are available for each ephemerides query:
165
+
166
+ .. code-block :: python
167
+
168
+ >> > print (eph.columns) # doctest: +REMOTE_DATA
168
169
< TableColumns names= (' targetname' ,' datetime_str' ,' datetime_jd' ,' H' ,' G' ,' solar_presence' ,' flags' ,' RA' ,' DEC' ,' RA_app' ,' DEC_app' ,' RA_rate' ,' DEC_rate' ,' AZ' ,' EL' ,' AZ_rate' ,' EL_rate' ,' sat_X' ,' sat_Y' ,' sat_PANG' ,' siderealtime' ,' airmass' ,' magextinct' ,' V' ,' surfbright' ,' illumination' ,' illum_defect' ,' sat_sep' ,' sat_vis' ,' ang_width' ,' PDObsLon' ,' PDObsLat' ,' PDSunLon' ,' PDSunLat' ,' SubSol_ang' ,' SubSol_dist' ,' NPole_ang' ,' NPole_dist' ,' EclLon' ,' EclLat' ,' r' ,' r_rate' ,' delta' ,' delta_rate' ,' lighttime' ,' vel_sun' ,' vel_obs' ,' elong' ,' elongFlag' ,' alpha' ,' lunar_elong' ,' lunar_illum' ,' sat_alpha' ,' sunTargetPA' ,' velocityPA' ,' OrbPlaneAng' ,' constellation' ,' TDB-UT' ,' ObsEclLon' ,' ObsEclLat' ,' NPole_RA' ,' NPole_DEC' ,' GlxLon' ,' GlxLat' ,' solartime' ,' earth_lighttime' ,' RA_3sigma' ,' DEC_3sigma' ,' SMAA_3sigma' ,' SMIA_3sigma' ,' Theta_3sigma' ,' Area_3sigma' ,' RSS_3sigma' ,' r_3sigma' ,' r_rate_3sigma' ,' SBand_3sigma' ,' XBand_3sigma' ,' DoppDelay_3sigma' ,' true_anom' ,' hour_angle' ,' alpha_true' ,' PABLon' ,' PABLat' )>
169
170
170
171
The values in these columns are the same as those defined in the Horizons
@@ -214,7 +215,6 @@ Horizons) and for a given epoch or a range of epochs (``epochs``) in the form of
214
215
an astropy table. The following example queries the osculating elements of
215
216
asteroid (433) Eros for a given date relative to the Sun:
216
217
217
- .. code-block :: python
218
218
.. doctest-remote-data ::
219
219
220
220
>>> from astroquery.jplhorizons import Horizons
@@ -226,8 +226,13 @@ asteroid (433) Eros for a given date relative to the Sun:
226
226
--- d ... AU d
227
227
------------------ ------------- ... ----------------- -----------------
228
228
433 Eros (A898 PA) 2458133.33546 ... 1.782442696867877 642.9387350660577
229
- >>> # The following fields are queried:
230
- >>> print (el.columns)
229
+
230
+
231
+ The following fields are queried:
232
+
233
+ .. code-block :: python
234
+
235
+ >> > print (el.columns) # doctest: +REMOTE_DATA
231
236
< TableColumns names= (' targetname' ,' datetime_jd' ,' datetime_str' ,' H' ,' G' ,' e' ,' q' ,' incl' ,' Omega' ,' w' ,' Tp_jd' ,' n' ,' M' ,' nu' ,' a' ,' Q' ,' P' )>
232
237
233
238
Optional parameters of :meth: `~astroquery.jplhorizons.HorizonsClass.elements `
@@ -253,7 +258,6 @@ the form of an astropy table. The following example queries the state
253
258
vector of asteroid 2012 TC4 as seen from Goldstone for a range of
254
259
epochs:
255
260
256
- .. code-block :: python
257
261
.. doctest-remote-data ::
258
262
259
263
>>> from astroquery.jplhorizons import Horizons
@@ -275,8 +279,13 @@ epochs:
275
279
(2012 TC4) 2458028.493055556 ... 0.03910796987365932 -0.004063569840032833
276
280
(2012 TC4) 2458028.5 ... 0.03907974874463358 -0.004064045454467904
277
281
Length = 145 rows
278
- >>> # The following fields are queried:
279
- >>> print (vec.columns)
282
+
283
+
284
+ The following fields are queried:
285
+
286
+ .. code-block :: python
287
+
288
+ >> > print (vec.columns) # doctest: +REMOTE_DATA
280
289
< TableColumns names= (' targetname' ,' datetime_jd' ,' datetime_str' ,' H' ,' G' ,' x' ,' y' ,' z' ,' vx' ,' vy' ,' vz' ,' lighttime' ,' range' ,' range_rate' )>
281
290
282
291
@@ -304,17 +313,21 @@ We provide some examples to illustrate how to use them based on the following
304
313
JPL Horizons ephemerides query of near-Earth asteroid (3552) Don Quixote since
305
314
its year of Discovery:
306
315
307
- .. code-block :: python
308
316
.. doctest-remote-data ::
309
317
310
318
>>> from astroquery.jplhorizons import Horizons
311
319
>>> obj = Horizons(id = ' 3552' , location = ' 568' ,
312
320
... epochs= {' start' :' 2010-01-01' , ' stop' :' 2019-12-31' ,
313
321
... ' step' :' 1y' })
314
322
>>> eph = obj.ephemerides()
315
- >>> # As we have seen before, we can display a truncated version of table
316
- >>> # ``eph`` by simply using
317
- >>> print (eph)
323
+
324
+
325
+ As we have seen before, we can display a truncated version of table
326
+ ``eph `` by simply using
327
+
328
+ .. code-block :: python
329
+
330
+ >> > print (eph) # doctest: +REMOTE_DATA
318
331
targetname datetime_str ... PABLon PABLat
319
332
-- - -- - ... deg deg
320
333
-------------------------- ---------------- - ... -------- --------
@@ -394,7 +407,6 @@ let's calculate the total rate of the object by summing 'RA_rate' and 'DEC_rate'
394
407
in quadrature:
395
408
396
409
397
- .. code-block :: python
398
410
.. doctest-remote-data ::
399
411
400
412
>>> import numpy as np
@@ -425,9 +437,8 @@ available, too, e.g., the ``RA_rate`` column is expressed in ``arcsec /
425
437
h `` - arcseconds per hour:
426
438
427
439
.. code-block :: python
428
- .. doctest-remote-data ::
429
440
430
- >>> print (eph[' RA_rate' ])
441
+ >> > print (eph[' RA_rate' ]) # doctest: +REMOTE_DATA
431
442
RA_rate
432
443
arcsec / h
433
444
----------
@@ -446,7 +457,6 @@ h`` - arcseconds per hour:
446
457
The unit of this column can be easily converted to any other unit describing the
447
458
same dimensions. For instance, we can turn ``RA_rate `` into ``arcsec / s ``:
448
459
449
- .. code-block :: python
450
460
.. doctest-remote-data ::
451
461
452
462
>>> eph[' RA_rate' ].convert_unit_to(' arcsec/s' )
@@ -568,7 +578,6 @@ For example, get the barycentric coordinates of Jupiter as an astropy
568
578
569
579
570
580
571
-
572
581
Acknowledgements
573
582
================
574
583
0 commit comments