Aquila-Consortium/map2map
32
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add plt_slices(size=None)

Browse Source
This commit is contained in:
Yin Li 2020-09-09 14:04:27 -04:00
parent 6c62cb09db
commit c4ab7e065b
1 changed files with 15 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
map2map/data/figures.py
View File

@ -79,19 +79,22 @@ def plt_slices(*fields, size=64, title=None, cmap=None, norm=None):
norm_col = SymLogNorm(linthresh=linthresh, vmin=-vlim, vmax=vlim) norm_col = SymLogNorm(linthresh=linthresh, vmin=-vlim, vmax=vlim)
for c in range(field.shape[0]): for c in range(field.shape[0]):
s0 = (c,) + tuple(d // 2 for d in field.shape[1:-2]) s = (c,) + tuple(d // 2 for d in field.shape[1:-2])
s1 = ( if size is None:
slice( s += (slice(None),) * 2
(field.shape[-2] - size) // 2, else:
(field.shape[-2] + size) // 2, s += (
), slice(
slice( (field.shape[-2] - size) // 2,
(field.shape[-1] - size) // 2, (field.shape[-2] + size) // 2,
(field.shape[-1] + size) // 2, ),
), slice(
) (field.shape[-1] - size) // 2,
(field.shape[-1] + size) // 2,
),
)
axes[c, f].pcolormesh(field[s0 + s1], cmap=cmap_col, norm=norm_col) axes[c, f].pcolormesh(field[s], cmap=cmap_col, norm=norm_col)
axes[c, f].set_aspect('equal') axes[c, f].set_aspect('equal')