def setUp(self): import pkg_resources import tempfile file_name = 'data/sample_data_NGC7023.tbl' file_path = pkg_resources.resource_filename('pypahdb', file_name) self.observation = Observation(file_path) self.decomposer = Decomposer(self.observation.spectrum) self.tmpdir = tempfile.gettempdir()
def test_read_fits(self): """Can we create an instance of Observation from a FITS file?""" file_name = 'resources/sample_data_NGC7023.fits' file_path = pkg_resources.resource_filename('pypahdb', file_name) assert isinstance(Observation(file_path), Observation)
def test_read_ascii(self): """Can we create an instance of Observation from an ASCII file?""" file_name = 'data/sample_data_NGC7023.tbl' file_path = pkg_resources.resource_filename('pypahdb', file_name) self.assertIsInstance(Observation(file_path), Observation)
#!/usr/bin/env python3 """ example.py Example of using pypahdb to decompose an astronomical PAH spectrum. """ import pkg_resources from pypahdb.decomposer import Decomposer from pypahdb.observation import Observation if __name__ == '__main__': # Sample data (in FITS format). file_path = 'resources/sample_data_NGC7023.fits' data_file = pkg_resources.resource_filename('pypahdb', file_path) # Construct an Observation object. obs = Observation(data_file) # Pass the Observation's spectrum to Decomposer, which performs the fit. pahdb_fit = Decomposer(obs.spectrum) # Write the results to file. # doplots is disabled due to significant CPU overhead pahdb_fit.save_pdf('NGC7023_pypahdb_fits_example.pdf', header=obs.header, doplots=False) pahdb_fit.save_fits('NGC7023_pypahdb_fits_example.fits', header=obs.header)