def test_getMapRatio(self, galaxy): maps = Maps(galaxy.plateifu) map_ratio = maps.getMapRatio('emline_gflux', 'nii_6585', 'ha_6564') map_arith = maps.emline_gflux_nii_6585 / maps.emline_gflux_ha_6564 assert map_ratio.value == pytest.approx(map_arith.value, nan_ok=True) assert map_ratio.ivar == pytest.approx(map_arith.ivar, nan_ok=True) assert map_ratio.mask == pytest.approx(map_arith.mask, nan_ok=True)
hbd = maps['emline_gsigma_hb_4862'] mapplot.plot(dapmap=hbd, fig=fig, ax=ax, cb_kws=cb_kws, log_cb=1) #cbrange=(0.1,100)) ax = fig.add_subplot(2,2,2) o1d = maps['emline_gsigma_oi_6302'] mapplot.plot(dapmap=o1d, fig=fig, ax=ax, cb_kws=cb_kws, log_cb=1) #cbrange=(0.1,100)) ax = fig.add_subplot(2,2,3) n2d = maps['emline_gsigma_nii_6549'] mapplot.plot(dapmap=n2d, fig=fig, ax=ax, cb_kws=cb_kws, log_cb=1) #cbrange=(0.1,100)) ax = fig.add_subplot(2,2,4) s2d = maps['emline_gsigma_sii_6718'] mapplot.plot(dapmap=s2d, fig=fig, ax=ax, cb_kws=cb_kws, log_cb=1) #cbrange=(0.1,100)) fig.tight_layout() pdf.savefig() plt.close() # page 5 # ha_6564:oiid_3728 normal plot fig = plt.figure() ax = fig.add_subplot(2, 2, 1) h_o_ratio = Maps.getMapRatio(maps,'emline_sflux','oiid_3728','ha_6564') mapplot.plot(dapmap=h_o_ratio, fig=fig, ax=ax, cb_kws=cb_kws) # ha_6564:oiid_3728 logarithmic plot ax = fig.add_subplot(2, 2, 2) h_o_ratio = Maps.getMapRatio(maps,'emline_sflux','oiid_3728','ha_6564') mapplot.plot(dapmap=h_o_ratio, fig=fig, ax=ax, cb_kws=cb_kws, log_cb=1, cbrange=(0.1,100)) fig.tight_layout() pdf.savefig() plt.close() os.system("open %s &" % filename)