Python normalize示例

示例#1

文件： plot_radial_distributions.py 项目： astrofrog/gammapy

"""Plot radial surface density distributions of Galactic sources."""
import matplotlib.pyplot as plt
import numpy as np
from gammapy.astro.population import radial_distributions
from gammapy.utils.distributions import normalize

# TODO: use Quantity
max_radius = 20  # kpc
r = np.linspace(0, max_radius, 100)

for key in radial_distributions:
    model = radial_distributions[key]()
    if model.evolved:
        linestyle = '-'
    else:
        linestyle = '--'
    label = model.__class__.__name__
    plt.plot(r, normalize(model, 0, max_radius)(r), linestyle=linestyle, label=label)
plt.xlim(0, max_radius)
plt.ylim(0, 0.28)
plt.xlabel('Galactocentric Distance [kpc]')
plt.ylabel('Normalized Surface Density [kpc^-2]')
plt.legend(prop={'size': 10})
plt.show()

示例#2

"""Plot velocity distributions of Galactic sources."""
import numpy as np
import matplotlib.pyplot as plt
from gammapy.astro.population import velocity_distributions
from gammapy.utils.distributions import normalize

# TODO: use Quantity
v_min, v_max = 10, 3000  # km / s
v = np.linspace(v_min, v_max, 200)

for key in velocity_distributions:
    model = velocity_distributions[key]()
    label = model.__class__.__name__
    plt.plot(v, normalize(model, v_min, v_max)(v), linestyle='-', label=label)

plt.xlim(v_min, v_max)
plt.ylim(0, 0.004)
plt.xlabel('Velocity [km/s]')
plt.ylabel('Probability Density [(km / s)^-1]')
plt.semilogx()
plt.legend(prop={'size': 10})
plt.show()

示例#3

"""Plot radial surface density distributions of Galactic sources."""
import matplotlib.pyplot as plt
import numpy as np
import astropy.units as u
from gammapy.astro.population import radial_distributions
from gammapy.utils.distributions import normalize

radius = np.linspace(0, 20, 100) * u.kpc

for key in radial_distributions:
    model = radial_distributions[key]()
    if model.evolved:
        linestyle = "-"
    else:
        linestyle = "--"
    label = model.__class__.__name__
    x = radius.value
    y = normalize(model, 0, radius[-1].value)(radius.value)
    plt.plot(x, y, linestyle=linestyle, label=label)

plt.xlim(0, radius[-1].value)
plt.ylim(0, 0.26)
plt.xlabel("Galactocentric Distance [kpc]")
plt.ylabel("Normalized Surface Density [kpc^-2]")
plt.legend(prop={"size": 10})
plt.show()

示例#4

文件： plot_radial_distributions.py 项目： vorugantia/gammapy

"""Plot radial surface density distributions of Galactic sources."""
import matplotlib.pyplot as plt
import numpy as np
from gammapy.astro.population import radial_distributions
from gammapy.utils.distributions import normalize

# TODO: use Quantity
max_radius = 20  # kpc
r = np.linspace(0, max_radius, 100)

for key in radial_distributions:
    model = radial_distributions[key]()
    if model.evolved:
        linestyle = '-'
    else:
        linestyle = '--'
    label = model.__class__.__name__
    plt.plot(r,
             normalize(model, 0, max_radius)(r),
             linestyle=linestyle,
             label=label)
plt.xlim(0, max_radius)
plt.ylim(0, 0.28)
plt.xlabel('Galactocentric Distance [kpc]')
plt.ylabel('Normalized Surface Density [kpc^-2]')
plt.legend(prop={'size': 10})
plt.show()

示例#5

文件： plot_velocity_distributions.py 项目： contrera/gammapy-docs

"""Plot velocity distributions of Galactic sources."""
import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u
from gammapy.astro.population import velocity_distributions
from gammapy.utils.distributions import normalize

velocity = np.linspace(10, 3000, 200) * u.km / u.s

for key in velocity_distributions:
    model = velocity_distributions[key]()
    label = model.__class__.__name__
    x = velocity.value
    y = normalize(model, velocity[0].value, velocity[-1].value)(velocity.value)
    plt.plot(x, y, linestyle='-', label=label)

plt.xlim(velocity[0].value, velocity[-1].value)
plt.ylim(0, 0.005)
plt.xlabel('Velocity [km/s]')
plt.ylabel('Probability Density [(km / s)^-1]')
plt.semilogx()
plt.legend(prop={'size': 10})
plt.show()

示例#6

文件： plot_velocity_distributions.py 项目： klepser/gammapy

"""Plot velocity distributions of Galactic sources."""
import numpy as np
import matplotlib.pyplot as plt
from gammapy.astro.population import velocity_distributions
from gammapy.utils.distributions import normalize

v_min, v_max = 10, 3000  # km / s
v = np.linspace(v_min, v_max, 200)
colors = ['b', 'k', 'g']

for color, key in zip(colors, velocity_distributions):
    model = velocity_distributions[key]()
    label = model.__class__.__name__
    plt.plot(v,
             normalize(model, v_min, v_max)(v),
             color=color,
             linestyle='-',
             label=label)

plt.xlim(v_min, v_max)
plt.ylim(0, 0.004)
plt.xlabel('Velocity [km/s]')
plt.ylabel('Probability Density [(km / s)^-1]')
plt.semilogx()
plt.legend(prop={'size': 10})
plt.show()

示例#7

"""Plot radial surface density distributions of Galactic sources."""
import matplotlib.pyplot as plt
import numpy as np
from gammapy.astro.population import radial_distributions
from gammapy.utils.distributions import normalize

max_radius = 20  # kpc
r = np.linspace(0, max_radius, 100)

for key in radial_distributions:
    model = radial_distributions[key]()
    if model.evolved:
        linestyle = '-'
    else:
        linestyle = '--'
    label = model.__class__.__name__
    plt.plot(r, normalize(model, 0, max_radius)(r), linestyle=linestyle, label=label)
plt.xlim(0, max_radius)
plt.ylim(0, 0.28)
plt.xlabel('Galactocentric Distance [kpc]')
plt.ylabel('Normalized Surface Density [kpc^-2]')
plt.legend(prop={'size': 10})
plt.show()