예제 #1
0
파일: tycho.py 프로젝트: lucabaldini/ximpol
def plot_doc():
    """
    """
    img = xFITSImage(os.path.join(XIMPOL_DATA, 'tycho_cmap.fits'))
    fig = img.plot(show=False)
    xFITSImage.add_label(fig, 'XIPE 300 ks')
    plt.show()
예제 #2
0
파일: tycho.py 프로젝트: pabell/ximpol
def plot_doc():
    """
    """
    img = xFITSImage(os.path.join(XIMPOL_DATA, 'tycho_cmap.fits'))
    fig = img.plot(show=False)
    xFITSImage.add_label(fig, 'XIPE 300 ks')
    plt.show()
예제 #3
0
def plot_doc():
    """
    """
    img = xFITSImage(os.path.join(XIMPOL_DATA, 'crab_complex_cmap.fits'))
    fig = img.plot(show=False)
    xFITSImage.add_label(fig, 'XIPE %d ks' %DURATION/1000.)
    plt.show()
예제 #4
0
파일: casa.py 프로젝트: lucabaldini/ximpol
def plot_doc():
    """
    """
    img = xFITSImage(os.path.join(XIMPOL_DATA, 'casa_cmap.fits'))
    fig = img.plot(show=False)
    #Option to draw the psf circle on the count map
    RAD_PSF = 11/60.
    fig.show_circles(350.769, 58.873, RAD_PSF/60., lw=2, color='white')
    fig.add_label(0.73,0.90, 'PSF', relative=True, size='x-large',
                                    color='white', horizontalalignment='left') 
    xFITSImage.add_label(fig, 'XIPE 250 ks')
    plt.show()
예제 #5
0
파일: casa.py 프로젝트: pabell/ximpol
def display():
    """Display the source model.
    """
    from ximpol.utils.matplotlib_ import pyplot as plt
    from ximpol.srcmodel.img import xFITSImage

    print(ROI_MODEL)
    fig = plt.figure('Energy spectrum')
    total_spectral_model.plot(logy=True, show=False, label='Total')
    nonthermal_spectral_model.plot(logy=True, show=False, label='Non-thermal')
    thermal_spectral_model.plot(logy=True, show=False, label='Thermal')
    plt.legend(bbox_to_anchor=(0.95, 0.95))
    fig = thermal_component.image.plot(show=False)
    xFITSImage.add_label(fig, 'Chandra 1.5-3.0 keV')
    fig = nonthermal_component.image.plot(show=False)
    xFITSImage.add_label(fig, 'Chandra 4.0-6.0 keV')
    img = xFITSImage(he_img_file_path)
    fig = img.plot(show=False)
    polarization_map.build_grid_sample(ROI_MODEL.ra, ROI_MODEL.dec)
    polarization_map.overlay_arrows(fig)
    plt.show()
예제 #6
0
파일: casa.py 프로젝트: pabell/ximpol
def display():
    """Display the source model.
    """
    from ximpol.utils.matplotlib_ import pyplot as plt
    from ximpol.srcmodel.img import xFITSImage

    print(ROI_MODEL)
    fig = plt.figure('Energy spectrum')
    total_spectral_model.plot(logy=True, show=False, label='Total')
    nonthermal_spectral_model.plot(logy=True, show=False, label='Non-thermal')
    thermal_spectral_model.plot(logy=True, show=False, label='Thermal')
    plt.legend(bbox_to_anchor=(0.95, 0.95))
    fig = thermal_component.image.plot(show=False)
    xFITSImage.add_label(fig, 'Chandra 1.5-3.0 keV')
    fig = nonthermal_component.image.plot(show=False)
    xFITSImage.add_label(fig, 'Chandra 4.0-6.0 keV')
    img = xFITSImage(he_img_file_path)
    fig = img.plot(show=False)
    polarization_map.build_grid_sample(ROI_MODEL.ra, ROI_MODEL.dec)
    polarization_map.overlay_arrows(fig)
    plt.show()