Update visualization codes

exps/LFNA/vis-synthetic.py

@@ -24,6 +24,7 @@ if str(lib_dir) not in sys.path:
     sys.path.insert(0, str(lib_dir))
 
 
+from datasets.synthetic_core import get_synthetic_env
 from datasets.synthetic_example import create_example_v1
 from utils.temp_sync import optimize_fn, evaluate_fn
 
@@ -169,6 +170,52 @@ def compare_cl(save_dir):
     os.system("{:} -pix_fmt yuv420p {xdir}/vis.webm".format(base_cmd, xdir=save_dir))
 
 
+def visualize_env(save_dir):
+    save_dir = Path(str(save_dir))
+    save_dir.mkdir(parents=True, exist_ok=True)
+
+    dynamic_env = get_synthetic_env()
+    min_t, max_t = dynamic_env.min_timestamp, dynamic_env.max_timestamp
+    for idx, (timestamp, (allx, ally)) in enumerate(tqdm(dynamic_env, ncols=50)):
+        dpi, width, height = 30, 1800, 1400
+        figsize = width / float(dpi), height / float(dpi)
+        LabelSize, LegendFontsize, font_gap = 80, 80, 5
+        fig = plt.figure(figsize=figsize)
+
+        cur_ax = fig.add_subplot(1, 1, 1)
+        allx, ally = allx[:, 0].numpy(), ally[:, 0].numpy()
+        cur_ax.scatter(
+            allx,
+            ally,
+            color="k",
+            linestyle="-",
+            alpha=0.99,
+            s=10,
+            label="timestamp={:05d}".format(idx),
+        )
+        cur_ax.set_xlabel("X", fontsize=LabelSize)
+        cur_ax.set_ylabel("Y", rotation=0, fontsize=LabelSize)
+        for tick in cur_ax.xaxis.get_major_ticks():
+            tick.label.set_fontsize(LabelSize - font_gap)
+            tick.label.set_rotation(10)
+        for tick in cur_ax.yaxis.get_major_ticks():
+            tick.label.set_fontsize(LabelSize - font_gap)
+        cur_ax.set_xlim(-10, 10)
+        cur_ax.set_ylim(-60, 60)
+        cur_ax.legend(loc=1, fontsize=LegendFontsize)
+
+        save_path = save_dir / "{:05d}".format(idx)
+        fig.savefig(str(save_path) + ".pdf", dpi=dpi, bbox_inches="tight", format="pdf")
+        fig.savefig(str(save_path) + ".png", dpi=dpi, bbox_inches="tight", format="png")
+        plt.close("all")
+    save_dir = save_dir.resolve()
+    base_cmd = "ffmpeg -y -i {xdir}/%05d.png -vf scale=1800:1400 -pix_fmt yuv420p -vb 5000k".format(
+        xdir=save_dir
+    )
+    os.system("{:} {xdir}/env.mp4".format(base_cmd, xdir=save_dir))
+    os.system("{:} {xdir}/vis.webm".format(base_cmd, xdir=save_dir))
+
+
 if __name__ == "__main__":
 
     parser = argparse.ArgumentParser("Visualize synthetic data.")
@@ -180,4 +227,5 @@ if __name__ == "__main__":
     )
     args = parser.parse_args()
 
+    visualize_env(os.path.join(args.save_dir, "vis-env"))
     compare_cl(os.path.join(args.save_dir, "compare-cl"))

lib/datasets/synthetic_core.py

@@ -0,0 +1,29 @@
+import copy
+from .synthetic_env import SyntheticDEnv
+from .math_dynamic_funcs import DynamicQuadraticFunc
+from .math_adv_funcs import ConstantFunc, ComposedSinFunc
+
+
+def get_synthetic_env(total_timestamp=1000, num_per_task=1000, mode=None):
+    mean_generator = ComposedSinFunc()
+    std_generator = ComposedSinFunc(min_amplitude=0.5, max_amplitude=1.5)
+    dynamic_env = SyntheticDEnv(
+        [mean_generator],
+        [[std_generator]],
+        num_per_task=num_per_task,
+        timestamp_config=dict(
+            min_timestamp=-0.5, max_timestamp=1.5, num=total_timestamp, mode=mode
+        ),
+    )
+    function = DynamicQuadraticFunc()
+    function_param = dict()
+    function_param[0] = ComposedSinFunc(
+        num_sin_phase=4, phase_shift=1.0, max_amplitude=1.0
+    )
+    function_param[1] = ConstantFunc(constant=0.9)
+    function_param[2] = ComposedSinFunc(
+        num_sin_phase=5, phase_shift=0.4, max_amplitude=0.9
+    )
+    function.set(function_param)
+    dynamic_env.set_oracle_map(copy.deepcopy(function))
+    return dynamic_env

lib/datasets/synthetic_env.py

@@ -43,6 +43,14 @@ class SyntheticDEnv(data.Dataset):
 
         self._oracle_map = None
 
+    @property
+    def min_timestamp(self):
+        return self._timestamp_generator.min_timestamp
+
+    @property
+    def max_timestamp(self):
+        return self._timestamp_generator.max_timestamp
+
     def set_oracle_map(self, functor):
         self._oracle_map = functor
 
@@ -61,7 +69,7 @@ class SyntheticDEnv(data.Dataset):
         index, timestamp = self._timestamp_generator[index]
         mean_list = [functor(timestamp) for functor in self._mean_functors]
         cov_matrix = [
-            [cov_gen(timestamp) for cov_gen in cov_functor]
+            [abs(cov_gen(timestamp)) for cov_gen in cov_functor]
             for cov_functor in self._cov_functors
         ]
 

lib/datasets/synthetic_utils.py

@@ -53,6 +53,14 @@ class TimeStamp(UnifiedSplit, data.Dataset):
         self._total_num = num
         UnifiedSplit.__init__(self, self._total_num, mode)
 
+    @property
+    def min_timestamp(self):
+        return self._min_timestamp
+
+    @property
+    def max_timestamp(self):
+        return self._max_timestamp
+
     def __iter__(self):
         self._iter_num = 0
         return self