catalog = astropy_tabular_data(apath('fitted_line_parameters_Chi2Constraints.ipac'), format='ascii.ipac') catalog.label='FittedLineParameters' catalog.style.color = 'green' catalog.style.marker = 'o' cube = load_data(hpath('APEX_H2CO_303_202_bl.fits')) cube.label='H2CO 303/202' cube2 = load_data(molpath('APEX_SiO_54.fits')) cube2.label='SiO' cube3 = load_data(hpath('APEX_13CO_matched_H2CO.fits')) cube3.label='13CO' higaltem = load_data('/Users/adam/work/gc/gcmosaic_temp_conv36.fits') dc = DataCollection([cube, catalog, cube2, cube3, higaltem]) dc.merge(cube,cube2,cube3) dc.add_link(LinkSame(cube.id['Galactic Longitude'], catalog.id['GLON'])) dc.add_link(LinkSame(cube.id['Galactic Latitude'], catalog.id['GLAT'])) def ms_to_kms(x): return x/1e3 def kms_to_ms(x): return x*1e3 dc.add_link(LinkTwoWay(cube.id['Vrad'], catalog.id['center'], ms_to_kms, kms_to_ms)) subset_tem_lt_60 = (catalog.id['temperature_chi2'] < 60) & (catalog.id['temperature_chi2'] > 10) & (catalog.id['area'] < 0.015) subset_tem_gt_60 = (catalog.id['temperature_chi2'] > 60) & (catalog.id['area'] < 0.015) app = GlueApplication(dc)
catalog = astropy_tabular_data( apath('fitted_line_parameters_Chi2Constraints.ipac'), format='ascii.ipac') catalog.label = 'FittedLineParameters' catalog.style.color = 'green' catalog.style.marker = 'o' cube = load_data(hpath('APEX_H2CO_303_202_bl.fits')) cube.label = 'H2CO 303/202' cube2 = load_data(molpath('APEX_SiO_54.fits')) cube2.label = 'SiO' cube3 = load_data(hpath('APEX_13CO_matched_H2CO.fits')) cube3.label = '13CO' higaltem = load_data('/Users/adam/work/gc/gcmosaic_temp_conv36.fits') dc = DataCollection([cube, catalog, cube2, cube3, higaltem]) dc.merge(cube, cube2, cube3) dc.add_link(LinkSame(cube.id['Galactic Longitude'], catalog.id['GLON'])) dc.add_link(LinkSame(cube.id['Galactic Latitude'], catalog.id['GLAT'])) def ms_to_kms(x): return x / 1e3 def kms_to_ms(x): return x * 1e3 dc.add_link( LinkTwoWay(cube.id['Vrad'], catalog.id['center'], ms_to_kms, kms_to_ms))
dendcube.label='Dendrogram Cube' # this label is redundant, it will be deleted upon merge # cube contains real WCS information; dendcube does not h2cocube = load_data(hpath('APEX_H2CO_303_202_bl.fits')) h2cocube.label='H2CO 303202 Cube' catalog = astropy_tabular_data(hpath('PPV_H2CO_Temperature.ipac'), format='ipac') catalog.label='Fitted Catalog' h2cocube.add_component(dendcube.get_component('structure'), 'structure') dc = DataCollection(dendrogram) dc.append(h2cocube) dc.append(catalog) dc.append(dendcube) dc.merge(h2cocube,dendcube) dc.merge(dendro, catalog) app = GlueApplication(dc) cube_viewer = app.new_data_viewer(ImageWidget) cube_viewer.add_data(h2cocube) h2cocube.join_on_key(dendro, 'structure', dendro.pixel_component_ids[0]) scatter = app.new_data_viewer(ScatterWidget) scatter.add_data(dendro) scatter.yatt = dendro.id['temperature_chi2'] scatter.xatt = catalog.id['r321303'] dendview = app.new_data_viewer(DendroWidget)