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

Add cmap and norm for all non positive images in plt_slices

Browse Source
This commit is contained in:
Yin Li 2020-09-23 18:10:53 -04:00
parent e40ea52190
commit c102fc588e
1 changed files with 10 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
map2map/utils/figures.py
View File

@ -1,5 +1,4 @@
from math import log2, log10, ceil from math import log2, log10, ceil
import warnings
import torch import torch
import numpy as np import numpy as np
import matplotlib import matplotlib
@ -46,15 +45,14 @@ def plt_slices(*fields, size=64, title=None, cmap=None, norm=None):
) )
for f, (field, cmap_col, norm_col) in enumerate(zip(fields, cmap, norm)): for f, (field, cmap_col, norm_col) in enumerate(zip(fields, cmap, norm)):
all_non_neg = (field >= 0).all() all_non_neg = np.all(field >= 0)
all_non_pos = (field <= 0).all() all_non_pos = np.all(field <= 0)
if cmap_col is None: if cmap_col is None:
if all_non_neg: if all_non_neg:
cmap_col = 'viridis' cmap_col = 'inferno'
elif all_non_pos: elif all_non_pos:
warnings.warn('no implementation for all non-positive values') cmap_col = 'inferno_r'
cmap_col = None
else: else:
cmap_col = 'RdBu_r' cmap_col = 'RdBu_r'
@ -68,8 +66,12 @@ def plt_slices(*fields, size=64, title=None, cmap=None, norm=None):
else: else:
norm_col = LogNorm(vmin=quantize(l2), vmax=quantize(h2)) norm_col = LogNorm(vmin=quantize(l2), vmax=quantize(h2))
elif all_non_pos: elif all_non_pos:
warnings.warn('no implementation for all non-positive values yet') if l1 < 0.1 * l2:
norm_col = None norm_col = Normalize(vmin=-quantize(-l2), vmax=0)
else:
norm_col = SymLogNorm(linthresh=quantize(-h2),
vmin=-quantize(-l2),
vmax=-quantize(-h2))
else: else:
vlim = quantize(max(-l2, h2)) vlim = quantize(max(-l2, h2))
if w1 > 0.1 * w2 or l1 * h1 >= 0: if w1 > 0.1 * w2 or l1 * h1 >= 0: