| import os |
| from typing import Optional, Sequence, Union, Dict |
| import math |
| import numpy as np |
| from matplotlib import pyplot as plt |
| from matplotlib.font_manager import FontProperties |
| from matplotlib.colors import ListedColormap |
| from matplotlib.patches import Patch |
| from .sevir_cmap import get_cmap, VIL_COLORS, VIL_LEVELS |
|
|
|
|
| HMF_COLORS = np.array([ |
| [82, 82, 82], |
| [252, 141, 89], |
| [255, 255, 191], |
| [145, 191, 219] |
| ]) / 255 |
|
|
| THRESHOLDS = (0, 16, 74, 133, 160, 181, 219, 255) |
|
|
|
|
| def plot_hit_miss_fa(ax, y_true, y_pred, thres): |
| mask = np.zeros_like(y_true) |
| mask[np.logical_and(y_true >= thres, y_pred >= thres)] = 4 |
| mask[np.logical_and(y_true >= thres, y_pred < thres)] = 3 |
| mask[np.logical_and(y_true < thres, y_pred >= thres)] = 2 |
| mask[np.logical_and(y_true < thres, y_pred < thres)] = 1 |
| cmap = ListedColormap(HMF_COLORS) |
| ax.imshow(mask, cmap=cmap) |
|
|
|
|
| def plot_hit_miss_fa_all_thresholds(ax, y_true, y_pred, **unused_kwargs): |
| fig = np.zeros(y_true.shape) |
| y_true_idx = np.searchsorted(THRESHOLDS, y_true) |
| y_pred_idx = np.searchsorted(THRESHOLDS, y_pred) |
| fig[y_true_idx == y_pred_idx] = 4 |
| fig[y_true_idx > y_pred_idx] = 3 |
| fig[y_true_idx < y_pred_idx] = 2 |
| |
| fig[np.logical_and(y_true < THRESHOLDS[1], y_pred < THRESHOLDS[1])] = 1 |
| cmap = ListedColormap(HMF_COLORS) |
| ax.imshow(fig, cmap=cmap) |
|
|
|
|
| def vis_sevir_seq( |
| save_path, |
| seq: Union[np.ndarray, Sequence[np.ndarray]], |
| label: Union[str, Sequence[str]] = "pred", |
| norm: Optional[Dict[str, float]] = None, |
| interval_real_time: float = 10.0, plot_stride=2, |
| label_rotation=0, |
| label_offset=(-0.06, 0.4), |
| label_avg_int=False, |
| fs=10, |
| max_cols=10, ): |
| """ |
| Parameters |
| ---------- |
| seq: Union[np.ndarray, Sequence[np.ndarray]] |
| shape = (T, H, W). Float value 0-1 after `norm`. |
| label: Union[str, Sequence[str]] |
| label for each sequence. |
| norm: Union[str, Dict[str, float]] |
| seq_show = seq * norm['scale'] + norm['shift'] |
| interval_real_time: float |
| The minutes of each plot interval |
| max_cols: int |
| The maximum number of columns in the figure. |
| """ |
|
|
| def cmap_dict(s): |
| return {'cmap': get_cmap(s, encoded=True)[0], |
| 'norm': get_cmap(s, encoded=True)[1], |
| 'vmin': get_cmap(s, encoded=True)[2], |
| 'vmax': get_cmap(s, encoded=True)[3]} |
|
|
| |
| |
| |
| |
|
|
| fontproperties = FontProperties() |
| fontproperties.set_family('serif') |
| |
| fontproperties.set_size(fs) |
| |
|
|
| if isinstance(seq, Sequence): |
| seq_list = [ele.astype(np.float32) for ele in seq] |
| assert isinstance(label, Sequence) and len(label) == len(seq) |
| label_list = label |
| elif isinstance(seq, np.ndarray): |
| seq_list = [seq.astype(np.float32), ] |
| assert isinstance(label, str) |
| label_list = [label, ] |
| else: |
| raise NotImplementedError |
| if label_avg_int: |
| label_list = [f"{ele1}\nAvgInt = {np.mean(ele2): .3f}" |
| for ele1, ele2 in zip(label_list, seq_list)] |
| |
| seq_list = [ele[::plot_stride, ...] for ele in seq_list] |
| seq_len_list = [len(ele) for ele in seq_list] |
|
|
| max_len = max(seq_len_list) |
|
|
| max_len = min(max_len, max_cols) |
| seq_list_wrap = [] |
| label_list_wrap = [] |
| seq_len_list_wrap = [] |
| for i, (seq, label, seq_len) in enumerate(zip(seq_list, label_list, seq_len_list)): |
| num_row = math.ceil(seq_len / max_len) |
| for j in range(num_row): |
| slice_end = min(seq_len, (j + 1) * max_len) |
| seq_list_wrap.append(seq[j * max_len: slice_end]) |
| if j == 0: |
| label_list_wrap.append(label) |
| else: |
| label_list_wrap.append("") |
| seq_len_list_wrap.append(min(seq_len - j * max_len, max_len)) |
|
|
| if norm is None: |
| norm = {'scale': 255, |
| 'shift': 0} |
| nrows = len(seq_list_wrap) |
| fig, ax = plt.subplots(nrows=nrows, |
| ncols=max_len, |
| figsize=(3 * max_len, 3 * nrows)) |
|
|
| for i, (seq, label, seq_len) in enumerate(zip(seq_list_wrap, label_list_wrap, seq_len_list_wrap)): |
| ax[i][0].set_ylabel(ylabel=label, fontproperties=fontproperties, rotation=label_rotation) |
| ax[i][0].yaxis.set_label_coords(label_offset[0], label_offset[1]) |
| for j in range(0, max_len): |
| if j < seq_len: |
| x = seq[j] * norm['scale'] + norm['shift'] |
| ax[i][j].imshow(x, **cmap_dict('vil')) |
| if i == len(seq_list) - 1 and i > 0: |
| ax[-1][j].set_title(f"Min {int(interval_real_time * (j + 1) * plot_stride)}", |
| y=-0.25, fontproperties=fontproperties) |
| else: |
| ax[i][j].axis('off') |
|
|
| for i in range(len(ax)): |
| for j in range(len(ax[i])): |
| ax[i][j].xaxis.set_ticks([]) |
| ax[i][j].yaxis.set_ticks([]) |
|
|
| |
| num_thresh_legend = len(VIL_LEVELS) - 1 |
| legend_elements = [Patch(facecolor=VIL_COLORS[i], |
| label=f'{int(VIL_LEVELS[i - 1])}-{int(VIL_LEVELS[i])}') |
| for i in range(1, num_thresh_legend + 1)] |
| ax[0][0].legend(handles=legend_elements, loc='center left', |
| bbox_to_anchor=(-1.2, -0.), |
| borderaxespad=0, frameon=False, fontsize='10') |
| plt.subplots_adjust(hspace=0.05, wspace=0.05) |
| plt.savefig(save_path) |
| plt.close(fig) |
|
|