本文整理汇总了Python中sunpy.map.Map.superpixel方法的典型用法代码示例。如果您正苦于以下问题:Python Map.superpixel方法的具体用法?Python Map.superpixel怎么用?Python Map.superpixel使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类sunpy.map.Map
的用法示例。
在下文中一共展示了Map.superpixel方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1:
# 需要导入模块: from sunpy.map import Map [as 别名]
# 或者: from sunpy.map.Map import superpixel [as 别名]
ax = fig.add_subplot(111, axisbg='gray')
plt.imshow(np.log10(model[:, side/2, :]), cmap='afmhot', interpolation='none',
vmin=0.0, vmax=1.0,
extent=[zrange[0], zrange[1], xrange[1], xrange[0]])
ax.add_artist(patches.Circle((0, 0), 1.0, color='black', fill=False))
plt.colorbar()
plt.savefig('rad-dist-model')
plt.close()
# Normalise each LOS
model /= np.nansum(model, axis=2)[:, :, None]
model.flush()
print 'Distribution model calculated'
# Distribute EM over LOS according to model
EMmap = EMmap.superpixel((4096./side, 4096./side), method='average')
EMlog = EMlog.superpixel((4096./side, 4096./side), method='average')
print np.nanmin(EMlog.data), np.nanmax(EMlog.data)
dz = ((z[0, 0, 1] - z[0, 0, 0]) * u.solRad).to(u.cm).value
model = (model * EMmap.data[:, :, None]) / dz
model = np.log10(np.sqrt(model))
np.save('density', model[:].copy())
print 'Density cube calculated'# and saved'
#for j in np.linspace(0, model.shape[1]-1, 7):
fig = plt.figure(figsize=(32, 24))
ax = fig.add_subplot(111, axisbg='gray')
print side/2, x[side/2]
slice = model[:, side/2, :]
#slice = 10.0**model[:, side/2, :]
plt.imshow(slice, cmap=emcm, interpolation='none',