Statistics on the Asteroids table

The table structure is the same as that of the MPC's MPCORB with the addition of the Tisserand (TJ) parameter

https://minorplanetcenter.net/iau/info/MPOrbitFormat.html

import SciServer 
import SciServer.CasJobs as CasJobs
import matplotlib.pyplot as plt
import numpy as np

Notes on the taxonomy

                    1  Atira
                    2  Aten
                    3  Apollo
                    4  Amor
                    5  Object with q < 1.665 AU
                    6  Hungaria
                    7  Unused or internal MPC use only
                    8  Hilda
                    9  Jupiter Trojan
                   10  Distant object

Atens, number / TJ

https://minorplanetcenter.net/iau/lists/t_atens.html

points = 100
Tjmin = 5.5
Tjmax = 9
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    Tj = Tjmin + (Tjmax-Tjmin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='2' and Tisserand >="+str(Tj)+" and Tisserand < "+str(Tj+(Tjmax-Tjmin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = Tj
    ypoints[i] = np.log10(result['n']+1)

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
    
plt.scatter(xpoints, ypoints)
plt.xlabel("Tj")
plt.ylabel("Log10 (n)")
plt.title('Aten number/Tj',fontsize=16)

Text(0.5, 1.0, 'Aten number/Tj')

Atens, number / inclination

points = 100
imin = 0
imax = 180
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    inc = imin + (imax-imin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='2' and i >="+str(inc)+" and i< "+str(inc+(imax-imin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = inc
    ypoints[i] = np.log10(result['n']+1)
    
pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(xpoints, ypoints)
plt.xlabel("i")
plt.ylabel("Log10 (n)")
plt.title('Aten number/inclination',fontsize=16)

Text(0.5, 1.0, 'Aten number/inclination')

Atens, inclination / eccentricity

query = "select e,i from Asteroids where substring(Taxonomy,5,1)='2'" 
result = SciServer.CasJobs.executeQuery(query, context="MyDB")

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(result['e'], result['i'])
plt.xlabel("e")
plt.ylabel("i")
plt.title('Aten inclination/eccentricity',fontsize=16)

Text(0.5, 1.0, 'Aten inclination/eccentricity')

Apollo, number / Tj

https://minorplanetcenter.net/iau/lists/Apollos.html

points =100
Tjmin = 1.6
Tjmax = 6.5
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    Tj = Tjmin + (Tjmax-Tjmin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='3' and Tisserand >="+str(Tj)+" and Tisserand < "+str(Tj+(Tjmax-Tjmin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = Tj
    ypoints[i] = np.log10(result['n']+1)

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(xpoints, ypoints)
plt.xlabel("Tj")
plt.ylabel("Log10 (n)")
plt.title('Apollo number/Tj',fontsize=16)

Text(0.5, 1.0, 'Apollo number/Tj')

Apollo, number / inclination

points = 100
imin = 0
imax = 180
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    inc = imin + (imax-imin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='3' and i >="+str(inc)+" and i< "+str(inc+(imax-imin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = inc
    ypoints[i] = np.log10(result['n']+1)

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(xpoints, ypoints)
plt.xlabel("i")
plt.ylabel("Log10 (n)")
plt.title('Apollo number/inclination',fontsize=16)

Text(0.5, 1.0, 'Apollo number/inclination')

Apollo, inclination / eccentricity

query = "select e,i from Asteroids where substring(Taxonomy,5,1)='3'" 
result = SciServer.CasJobs.executeQuery(query, context="MyDB")


pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(result['e'], result['i'])
plt.xlabel("e")
plt.ylabel("i")
plt.title('Apollo inclination/eccentricy',fontsize=16)

Text(0.5, 1.0, 'Apollo inclination/eccentricy')

Amor, number / TJ

https://minorplanetcenter.net/iau/lists/Amors.html

points = 100
Tjmin = 2
Tjmax = 6
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    Tj = Tjmin + (Tjmax-Tjmin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='4' and Tisserand >="+str(Tj)+" and Tisserand < "+str(Tj+(Tjmax-Tjmin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = Tj
    ypoints[i] = np.log10(result['n']+1)

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(xpoints, ypoints)
plt.xlabel("i")
plt.ylabel("Log10 (n)")
plt.title('Amor number/Tj',fontsize=16)

Text(0.5, 1.0, 'Amor number/Tj')

Amor, number / inclination

points = 100
imin = 0
imax = 180
xpoints = np.zeros(points)
ypoints = np.zeros(points)

for i in range(0, points):
    
    inc = imin + (imax-imin)*i/points
    
    query = "select count(*) as n from Asteroids where substring(Taxonomy,5,1)='4' and i >="+str(inc)+" and i< "+str(inc+(imax-imin)/points)
    result = SciServer.CasJobs.executeQuery(query, context="MyDB")
    
    xpoints[i] = inc
    ypoints[i] = np.log10(result['n']+1)

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))
plt.scatter(xpoints, ypoints)
plt.xlabel("i")
plt.ylabel("Log10 (n)")
plt.title('Amor number/inclination',fontsize=16)

Text(0.5, 1.0, 'Amor number/inclination')

Amor, inclination / eccentricity

query = "select e,i from Asteroids where substring(Taxonomy,5,1)='4'" 
result = SciServer.CasJobs.executeQuery(query, context="MyDB")

pixelsx = 1920
pixelsy = 1080
dpi = 90
fig = plt.figure(figsize=(pixelsx/dpi, pixelsy/dpi))

plt.scatter(result['e'], result['i'])
plt.xlabel("e")
plt.ylabel("i")
plt.title('Amor inclination/eccentricy',fontsize=16)

Text(0.5, 1.0, 'Amor inclination/eccentricy')