def import_iris_obs(): caption = "Select a directory containing files from one IRIS OBS, and stack all raster scans." data_path = Path(pick_directory(caption)) rasters = list(data_path.glob("*raster*")) sji = list(data_path.glob("*SJI*")) sji_data = [] for s in sji: sji_data.append(load_sji_fits(s)) raster_data = read_iris_spectrograph_level2_fits(rasters, spectral_windows=['Mg II k 2796'], memmap=False, uncertainty=False) raster_data = {window: stack_spectrogram_sequence(seq) for window, seq in raster_data.data.items()} result = [] for window, window_data in raster_data.items(): w_data = Data(label=f"{window.replace(' ', '_')}") w_data.coords = WCSCoordinates(wcs=window_data.wcs) w_data.add_component(Component(window_data.data), f"{window}") result.append(w_data) return result + sji_data
def load_rasters(self, windows): raster_data = read_iris_spectrograph_level2_fits(self.raster_files, spectral_windows=windows, memmap=False, uncertainty=False) raster_data = raster_data.data if self.stack.checkState() > 0: raster_data = {window: stack_spectrogram_sequence(seq) for window, seq in raster_data.items()} self.load_stacked_sequence(raster_data) else: self.load_sequence(raster_data)
def import_iris_obs(): caption = "Select a directory containing files from one IRIS OBS." data_path = Path(pick_directory(caption)) rasters = list(data_path.glob("*raster*")) sji = list(data_path.glob("*SJI*")) sji_data = [] for s in sji: sji_data.append(load_sji_fits(s)) raster_data = _parse_iris_raster( read_iris_spectrograph_level2_fits(rasters, uncertainty=False, memmap=True), 'iris') return raster_data + sji_data
def iris_l2_test_raster(): return read_iris_spectrograph_level2_fits( os.path.join(testpath, 'iris_l2_20170222_153635_3690215148_raster_t000_r00000.fits'))
def read_iris_raster(raster_file): raster_data = _parse_iris_raster( read_iris_spectrograph_level2_fits(raster_file), 'iris') return raster_data
def read_iris_raster(raster_file): raster_data = _parse_iris_raster( read_iris_spectrograph_level2_fits(raster_file, uncertainty=False, memmap=True), 'iris') return raster_data
from glue.app.qt.application import GlueApplication from glue.core import Data, DataCollection from glue.core.data_factories import load_data from glue.plugins.wcs_autolinking.wcs_autolinking import wcs_autolink from glue_solar.instruments.iris import _parse_iris_raster from irispy.spectrograph import read_iris_spectrograph_level2_fits from glue.viewers.image.qt import ImageViewer # create some data data_path = Path("/home/stuart/sunpy/data/iris_glue/") rasters = list(data_path.glob("*raster*")) sji = list(data_path.glob("*SJI*")) raster_data = _parse_iris_raster(read_iris_spectrograph_level2_fits(rasters), 'iris') sji_data = [] for s in sji: sji_data += load_data(s) dc = DataCollection(raster_data + sji_data) ga = GlueApplication(dc) im1 = ga.new_data_viewer(ImageViewer) im1.add_data(raster_data[0]) im2 = ga.new_data_viewer(ImageViewer) im2.add_data(sji_data[0])