def read_cube(filename, **kwargs):
cube_data = CubeData.read(filename)
data = Data()
data.add_component(Component(cube_data), label="cube")
print("Loaded successfully")
return data
from cube_tools.core.data_objects import CubeData, SpectrumData
import numpy as np
if __name__ == '__main__':
cube_data = CubeData.read("/Users/nearl/Desktop/cube_tools_demo_data/manga-7443-12703-LOGCUBE.fits",
format='fits')
spec_data = SpectrumData.read(
"/Users/nearl/Desktop/cube_tools_demo_data/manga-7443-12703-LOGCUBE.fits",
format='fits')
print('disp', cube_data[:, 0, 0].dispersion[0:10])
print('flux', cube_data[:, 0, 0].get_flux())
# import matplotlib.pyplot as plt
#
# plt.plot(cube_data[:, 50, 50].dispersion,
# cube_data[:, 50, 50].flux)
# plt.show()
