Load spectral data

In [1]:

from bh_molecule.instruments.vis133m import Vis133M
from bh_molecule import use_dark
use_dark()

from bh_molecule.instruments.vis133m import Vis133M from bh_molecule import use_dark use_dark()

Prep data¶

In [2]:

FITS = './data/example.fits'
wavecal_csv = './data/133mVis_wavcal.csv'
s26 = Vis133M(FITS,wavecal_csv)
s26.set_scale(0.0001837)
s26.set_dark(frame=49, channel=43)
s26.set_baseline_zero(True)
s26.set_time_linspace(0, 10)

FITS = './data/example.fits' wavecal_csv = './data/133mVis_wavcal.csv' s26 = Vis133M(FITS,wavecal_csv) s26.set_scale(0.0001837) s26.set_dark(frame=49, channel=43) s26.set_baseline_zero(True) s26.set_time_linspace(0, 10)

Pixel range C-F map¶

• C: channel
• F: frame (time)

In [3]:

s26.plot_pixel_range(456,642,log_scale=True,vmin=5e-3,vmax=10,cmap='gist_gray')

s26.plot_pixel_range(456,642,log_scale=True,vmin=5e-3,vmax=10,cmap='gist_gray')

Out[3]:

<Axes: xlabel='time (s)', ylabel='channel'>

Band map C-F¶

In [4]:

s26.plot_band_map((433.9, 435),log_scale=True,vmin=5e-3,vmax=10,cmap='gist_gray')

s26.plot_band_map((433.9, 435),log_scale=True,vmin=5e-3,vmax=10,cmap='gist_gray')

Out[4]:

<Axes: xlabel='time (s)', ylabel='channel'>

Channel stack¶

In [5]:

s26.plot_channel_stack(30,log_scale=True,vmin=0.002,cmap='gist_gray')

s26.plot_channel_stack(30,log_scale=True,vmin=0.002,cmap='gist_gray')

Out[5]:

<Axes: xlabel='time (s)', ylabel='wavelength (nm)'>

In [6]:

# Only for mkdocs, you don't need this in your own code
from IPython.display import HTML
import plotly.io as pio

def show_plotly(fig, include_js="cdn"):
    """Return an HTML chunk MkDocs/nbconvert will render."""
    return HTML(pio.to_html(fig, include_plotlyjs=include_js, full_html=False))

# Only for mkdocs, you don't need this in your own code from IPython.display import HTML import plotly.io as pio def show_plotly(fig, include_js="cdn"): """Return an HTML chunk MkDocs/nbconvert will render.""" return HTML(pio.to_html(fig, include_plotlyjs=include_js, full_html=False))

Interactive spectrum¶

In [7]:

fig = s26.plot_spectrum_plotly(31, 30,theme='dark')
show_plotly(fig) # only for mkdocs, you don't need this in your own code
# fig.show() # for local use

fig = s26.plot_spectrum_plotly(31, 30,theme='dark') show_plotly(fig) # only for mkdocs, you don't need this in your own code # fig.show() # for local use

Out[7]: