當前位置: 首頁>>代碼示例>>Python>>正文


Python Map.meta['crval1']方法代碼示例

本文整理匯總了Python中sunpy.map.Map.meta['crval1']方法的典型用法代碼示例。如果您正苦於以下問題:Python Map.meta['crval1']方法的具體用法?Python Map.meta['crval1']怎麽用?Python Map.meta['crval1']使用的例子?那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在sunpy.map.Map的用法示例。


在下文中一共展示了Map.meta['crval1']方法的1個代碼示例,這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚,您的評價將有助於係統推薦出更棒的Python代碼示例。

示例1: enumerate

# 需要導入模塊: from sunpy.map import Map [as 別名]
# 或者: from sunpy.map.Map import meta['crval1'] [as 別名]
#rect = patches.Rectangle([25.0, 5.6], 1.0, 1.0, color='black', fill=True, clip_on=False)
# Run synthetic data through 1param tempmap method
for w, wid in enumerate(widths):#heights):
    print '\nWidth:', wid
    fig = plt.figure(figsize=(30, 12))
    for wl, wlength  in enumerate(['94', '131', '171', '193', '211', '335']):
        #emiss = Map(emission[wl, :, :, w], mapmeta)
        emiss = Map(emission[wl, :, w, :].copy(), mapmeta)
        emiss.cmap = sunpy.cm.get_cmap('sdoaia{}'.format(wlength))
        emiss.meta['naxis1'] = emiss.shape[1]
        emiss.meta['naxis2'] = emiss.shape[0]
        #emiss.meta['cdelt1'] = widths[1] - widths[0]
        emiss.meta['cdelt1'] = heights[1] - heights[0] #np.log10(heights[1]) - np.log10(heights[0])
        emiss.meta['cdelt2'] = temps[1] - temps[0]
        #emiss.meta['crval1'] = widths[0]
        emiss.meta['crval1'] = heights[0] #np.log10(heights[0])
        emiss.meta['crval2'] = temps[0]
        emiss.meta['crpix1'] = 0.5
        emiss.meta['crpix2'] = 0.5
        if wlength == '94': wlength = '094'
        fits_dir = path.join(CThome, 'data', 'synthetic', wlength)
        if not path.exists(fits_dir): makedirs(fits_dir)
        emiss.save(path.join(fits_dir, 'model.fits'), clobber=True)
        #print '----', emission[2, :, :, :].min(), emission[2, :, :, :].max()
        #print '------', emission[2, :, w, :].min(), emission[2, :, w, :].max()
        emiss.data /= emission[2, :, w, :]
        #print '--------', emiss.min(), emiss.max()
        ax = fig.add_subplot(1, 6, wl+1)
        emiss.plot(aspect='auto', vmin=emiss.min(), vmax=emiss.max())
        plt.title('{}'.format(wlength))
        plt.xlabel('Input EM')
開發者ID:Cadair,項目名稱:CoronaTemps,代碼行數:33,代碼來源:test_vs_model_DEMs.py


注:本文中的sunpy.map.Map.meta['crval1']方法示例由純淨天空整理自Github/MSDocs等開源代碼及文檔管理平台,相關代碼片段篩選自各路編程大神貢獻的開源項目,源碼版權歸原作者所有,傳播和使用請參考對應項目的License;未經允許,請勿轉載。