Skip to content

SSDataLab/tess-ephem

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

81 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

tess-ephem

Where are Solar System objects located in TESS FFI data?

pypi pytest black flake8 mypy

tess-ephem is a user-friendly package which enables users to compute the positions of Solar System objects -- asteroids, comets, and planets -- in the data archive of NASA's TESS Space Telescope.

Installation

python -m pip install tess-ephem

Example use

tess-ephem allows you to search the entire archive of TESS FFI's for the presence of a known minor planet, and obtain the result as a Pandas DataFrame. The output pixel coordinates (column and row) follow the TESS convention, with (1,1) being the middle of the pixel in the lower left corner of the FFI. For example:

>>> from tess_ephem import ephem
>>> ephem("Sedna")
           sector  camera  ccd       column          row  pixels_per_hour        ra      dec    vmag  hmag  sun_distance  obs_distance  phase_angle
time
2458437.5       5       1    4  1541.303105  1102.826838         0.069068  57.12707  7.65227  20.803  1.49     84.948291     83.983859       0.1482
2458438.5       5       1    4  1543.100766  1102.946276         0.074363  57.11724  7.64961  20.803  1.49     84.947515     83.981039       0.1450
2458439.5       5       1    4  1544.975690  1102.998615         0.075353  57.10692  7.64724  20.802  1.49     84.946739     83.978850       0.1423
2458440.5       5       1    4  1546.862705  1103.025319         0.075362  57.09651  7.64500  20.801  1.49     84.945964     83.977315       0.1403
2458441.5       5       1    4  1548.747559  1103.038136         0.075132  57.08610  7.64284  20.801  1.49     84.945188     83.976314       0.1394
...           ...     ...  ...          ...          ...              ...       ...      ...     ...   ...           ...           ...          ...
2460254.5      71       2    4  1984.592691  1003.520110         0.072528  60.21660  8.36089  20.750  1.49     83.592715     82.661747       0.2395
2460255.5      71       2    4  1984.822681  1001.701217         0.072889  60.20673  8.35803  20.749  1.49     83.591992     82.656583       0.2304
2460256.5      71       2    4  1985.051503   999.875207         0.073002  60.19682  8.35517  20.747  1.49     83.591270     82.651896       0.2214
2460257.5      71       2    4  1985.272750   998.048780         0.072664  60.18690  8.35235  20.745  1.49     83.590547     82.647726       0.2127
2460258.5      71       2    4  1985.475079   996.244021         0.071257  60.17708  8.34965  20.744  1.49     83.589824     82.644125       0.2044

[78 rows x 13 columns]

You can also obtain the ephemeris for one or more specific times by passing the time parameter:

>>> ephem("Sedna", time="2018-11-21 17:35:00")
                         sector  camera  ccd       column          row  pixels_per_hour        ra      dec    vmag  hmag  sun_distance  obs_distance  phase_angle
time
2018-11-21 17:35:00.000       5       1    4  1553.855013  1103.032894         0.074053  57.05786  7.63721  20.802  1.49     84.943069     83.975873       0.1419

>>> from astropy.time import Time
>>> ephem("Sedna", time=Time([2458441.5,2460258.5], format='jd'))
           sector  camera  ccd       column          row  pixels_per_hour        ra      dec    vmag  hmag  sun_distance  obs_distance  phase_angle
time
2458441.5       5       1    4  1548.747559  1103.038136         0.075132  57.08610  7.64284  20.801  1.49     84.945188     83.976314       0.1394
2460258.5      71       2    4  1985.475079   996.244021         0.071257  60.17708  8.34965  20.744  1.49     83.589824     82.644125       0.2044

Orbital elements can be obtained by passing the orbital_elements=True parameter. The function returns the average orbital elements of the target during the queried time. Perihelion distance is in AU and orbital inclination is in degrees.

>>> df_ephem, orbital_elements = ephem("Sedna", time="2018-11-21 17:35:00", orbital_elements=True)
>>> orbital_elements
{'perihelion_distance': 76.13300496827914, 'eccentricity': 0.8413468330602263, 'orbital_inclination': 11.929585189509927}

You can alternatively obtain the ephemeris during a specific sector by passing the sector parameter:

>>> ephem("Sedna", sector=70)
           sector  camera  ccd       column          row  pixels_per_hour        ra      dec    vmag  hmag  sun_distance  obs_distance  phase_angle
time
2460208.5      70       4    2  1965.969985  1827.030603         0.025550  60.53693  8.49483  20.815  1.49     83.625984     83.171959       0.6161
2460209.5      70       4    2  1966.279804  1826.423362         0.029493  60.53386  8.49248  20.814  1.49     83.625260     83.156129       0.6112
2460210.5      70       4    2  1966.614906  1825.724165         0.032766  60.53030  8.48989  20.813  1.49     83.624537     83.140186       0.6058
2460211.5      70       4    2  1966.971035  1824.956892         0.035144  60.52638  8.48711  20.812  1.49     83.623813     83.124437       0.6001
2460212.5      70       4    2  1967.339271  1824.141075         0.037045  60.52220  8.48421  20.811  1.49     83.623090     83.108950       0.5940
2460213.5      70       4    2  1967.718547  1823.280470         0.038666  60.51778  8.48120  20.810  1.49     83.622366     83.093746       0.5878
2460214.5      70       4    2  1968.103146  1822.383332         0.040092  60.51316  8.47812  20.809  1.49     83.621642     83.078838       0.5813
2460215.5      70       4    2  1968.491641  1821.451723         0.041354  60.50835  8.47498  20.808  1.49     83.620919     83.064240       0.5747
2460216.5      70       4    2  1968.884752  1820.489031         0.042441  60.50337  8.47178  20.806  1.49     83.620195     83.049971       0.5680
2460217.5      70       4    2  1969.276403  1819.501819         0.043287  60.49825  8.46856  20.805  1.49     83.619472     83.036058       0.5611
2460218.5      70       4    2  1969.665886  1818.495541         0.043699  60.49302  8.46533  20.804  1.49     83.618748     83.022538       0.5541
2460219.5      70       4    2  1970.045345  1817.485982         0.043112  60.48776  8.46215  20.803  1.49     83.618025     83.009457       0.5471
2460220.5      70       4    2  1970.388690  1816.508633         0.040265  60.48264  8.45920  20.802  1.49     83.617302     82.996750       0.5403
2460221.5      70       4    2  1970.660834  1815.574303         0.039919  60.47768  8.45668  20.800  1.49     83.616578     82.983656       0.5337
2460222.5      70       4    2  1970.920335  1814.571748         0.043872  60.47232  8.45415  20.799  1.49     83.615855     82.969867       0.5266
2460223.5      70       4    2  1971.201557  1813.486767         0.046801  60.46652  8.45141  20.798  1.49     83.615131     82.956089       0.5192
2460224.5      70       4    2  1971.503820  1812.342507         0.048934  60.46041  8.44849  20.796  1.49     83.614408     82.942565       0.5115
2460225.5      70       4    2  1971.815329  1811.153797         0.050643  60.45406  8.44547  20.795  1.49     83.613685     82.929353       0.5035
2460226.5      70       4    2  1972.135377  1809.926088         0.052095  60.44750  8.44236  20.794  1.49     83.612961     82.916471       0.4954
2460227.5      70       4    2  1972.459314  1808.663881         0.053363  60.44075  8.43919  20.792  1.49     83.612238     82.903932       0.4871
2460228.5      70       4    2  1972.786794  1807.374325         0.054473  60.43385  8.43597  20.791  1.49     83.611515     82.891750       0.4787
2460229.5      70       4    2  1973.114235  1806.058161         0.055414  60.42680  8.43272  20.790  1.49     83.610791     82.879944       0.4701
2460230.5      70       4    2  1973.442721  1804.720474         0.056123  60.41963  8.42944  20.788  1.49     83.610068     82.868542       0.4615
2460231.5      70       4    2  1973.766531  1803.369529         0.056414  60.41238  8.42617  20.787  1.49     83.609345     82.857585       0.4528
2460232.5      70       4    2  1974.077449  1802.019428         0.055711  60.40512  8.42297  20.785  1.49     83.608622     82.847125       0.4442

When passing the sector parameter, the time_step is by default 1 day. This can be changed as follows:

>>> ephem("Sedna", sector=70, time_step=0.1)
           sector  camera  ccd       column          row  pixels_per_hour         ra       dec       vmag  hmag  sun_distance  obs_distance  phase_angle
time
2460207.6      70       4    2  1965.673822  1827.524121         0.026894  60.539373  8.496977  20.815973  1.49     83.626635     83.185495     0.620233
2460207.7      70       4    2  1965.710011  1827.468236         0.026487  60.539100  8.496720  20.816002  1.49     83.626563     83.184025     0.619771
2460207.8      70       4    2  1965.745339  1827.412939         0.026102  60.538829  8.496468  20.816045  1.49     83.626491     83.182549     0.619312
2460207.9      70       4    2  1965.779762  1827.358109         0.025764  60.538559  8.496221  20.816058  1.49     83.626418     83.181067     0.618856
2460208.0      70       4    2  1965.813240  1827.303629         0.025498  60.538290  8.495980  20.816000  1.49     83.626346     83.179575     0.618400
...           ...     ...  ...          ...          ...              ...        ...       ...        ...   ...           ...           ...          ...
2460233.0      70       4    2  1974.223608  1801.353369         0.054524  60.401530  8.421430  20.785000  1.49     83.608260     82.842086     0.439900
2460233.1      70       4    2  1974.250923  1801.222026         0.054170  60.400820  8.421134  20.784886  1.49     83.608188     82.841090     0.439050
2460233.2      70       4    2  1974.277356  1801.091524         0.053786  60.400114  8.420844  20.784707  1.49     83.608115     82.840098     0.438204
2460233.3      70       4    2  1974.302887  1800.961890         0.053386  60.399412  8.420560  20.784486  1.49     83.608043     82.839105     0.437364
2460233.4      70       4    2  1974.327497  1800.833148         0.052985  60.398714  8.420282  20.784243  1.49     83.607971     82.838112     0.436529

[259 rows x 13 columns]

About

Where are Solar System objects located in TESS data?

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Contributors 4

  •  
  •  
  •  
  •