From 4b95042ad2099b3100163b50090b68b2d6fe1002 Mon Sep 17 00:00:00 2001
From: Adam Tyson <code@adamltyson.com>
Date: Tue, 11 Jul 2023 10:43:57 +0100
Subject: [PATCH] Add save option on analysis (#114)

* Add save option on analysis

* Add (broken) save test

* Update tracks test

* update testing

* Add sleep to fix tests

* Update test data

* Add remaining save tests

* set layer editability later to prevent it being reset

* sleep for longer to see if tests pass on windows/macos

* rename tests

* increase sleep
---
 brainreg_segment/layout/gui_constants.py      |   2 +
 brainreg_segment/segment.py                   |  25 +++---
 .../segmentation_panels/regions.py            |  19 ++++-
 .../segmentation_panels/tracks.py             |  29 +++++--
 tests/conftest.py                             |  11 ++-
 .../standard_space/regions/summary.csv        |   2 +-
 .../standard_space/tracks/test_track.h5       | Bin 39240 -> 0 bytes
 .../standard_space/tracks/test_track.npy      | Bin 0 -> 24128 bytes
 ...ons.py => test_region_ui_functionality.py} |  76 +++++++++++++++---
 ...acks.py => test_track_ui_functionality.py} |  60 +++++++++++---
 10 files changed, 178 insertions(+), 46 deletions(-)
 delete mode 100644 tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.h5
 create mode 100644 tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.npy
 rename tests/tests/test_integration/test_gui/{test_regions.py => test_region_ui_functionality.py} (58%)
 rename tests/tests/test_integration/test_gui/{test_tracks.py => test_track_ui_functionality.py} (69%)

diff --git a/brainreg_segment/layout/gui_constants.py b/brainreg_segment/layout/gui_constants.py
index c4b22cf..f82ac8a 100644
--- a/brainreg_segment/layout/gui_constants.py
+++ b/brainreg_segment/layout/gui_constants.py
@@ -14,6 +14,8 @@
     "left"  # Alignment of text in pushbuttons in methods chooser panel
 )
 
+SAVE_DEFAULT = False
+
 POINT_SIZE = 100
 SPLINE_SIZE = 50
 NUM_COLORS = 10
diff --git a/brainreg_segment/segment.py b/brainreg_segment/segment.py
index ba786aa..49fc95b 100644
--- a/brainreg_segment/segment.py
+++ b/brainreg_segment/segment.py
@@ -451,10 +451,9 @@ def initialise_loaded_data(self):
             ]
             self.hemispheres_data = self.hemispheres_layer.data
 
-        self.prevent_layer_edit()
-
         self.initialise_segmentation_interface()
         self.status_label.setText("Ready")
+        self.prevent_layer_edit()
 
     def collate_widget_layers(self):
         """
@@ -486,6 +485,7 @@ def collate_widget_layers(self):
         ]
 
     def prevent_layer_edit(self):
+        print("Preventing layer edit")
         self.collate_widget_layers()
         for layer in self.non_editable_widget_layers:
             layer.editable = False
@@ -568,18 +568,21 @@ def save(self, override=True):
             else:
                 choice = True  # for debugging
             if choice:
-                print("Saving")
-                worker = save_all(
-                    self.paths.regions_directory,
-                    self.paths.tracks_directory,
-                    self.label_layers,
-                    self.track_layers,
-                    track_file_extension=TRACK_FILE_EXT,
-                )
-                worker.start()
+                self.run_save()
             else:
                 print('Not saving because user chose "Cancel" \n')
 
+    def run_save(self):
+        print("Saving")
+        worker = save_all(
+            self.paths.regions_directory,
+            self.paths.tracks_directory,
+            self.label_layers,
+            self.track_layers,
+            track_file_extension=TRACK_FILE_EXT,
+        )
+        worker.start()
+
     def export_to_brainrender(self, override=False):
         if not override:
             choice = display_warning(
diff --git a/brainreg_segment/segmentation_panels/regions.py b/brainreg_segment/segmentation_panels/regions.py
index 11aea94..d92af15 100644
--- a/brainreg_segment/segmentation_panels/regions.py
+++ b/brainreg_segment/segmentation_panels/regions.py
@@ -7,6 +7,7 @@
     COLUMN_WIDTH,
     IMAGE_FILE_EXT,
     NUM_COLORS,
+    SAVE_DEFAULT,
     SEGM_METHODS_PANEL_ALIGN,
     SUMMARIZE_VOLUMES_DEFAULT,
 )
@@ -34,6 +35,7 @@ def __init__(
         parent,
         calculate_volumes_default=CALCULATE_VOLUMES_DEFAULT,
         summarise_volumes_default=SUMMARIZE_VOLUMES_DEFAULT,
+        save_default=SAVE_DEFAULT,
         brush_size=BRUSH_SIZE,
         image_file_extension=IMAGE_FILE_EXT,
         num_colors=NUM_COLORS,
@@ -43,6 +45,7 @@ def __init__(
 
         self.calculate_volumes_default = calculate_volumes_default
         self.summarise_volumes_default = summarise_volumes_default
+        self.save_default = save_default
 
         # Brushes / ...
         self.brush_size_default = BRUSH_SIZE  # Keep track of default
@@ -60,7 +63,7 @@ def add_region_panel(self, row):
             "Add new region",
             region_layout,
             self.add_new_region,
-            row=2,
+            row=3,
             column=0,
             tooltip="Create a new empty segmentation layer "
             "to manually segment a new region.",
@@ -70,7 +73,7 @@ def add_region_panel(self, row):
             "Analyse regions",
             region_layout,
             self.run_region_analysis,
-            row=2,
+            row=3,
             column=1,
             tooltip="Analyse the spatial distribution of the "
             "segmented regions.",
@@ -79,7 +82,7 @@ def add_region_panel(self, row):
             "Add region from selected layer",
             region_layout,
             self.add_region_from_existing_layer,
-            row=3,
+            row=4,
             column=0,
             tooltip="Adds a region from a selected labels layer (e.g. "
             "from another plugin). Make sure this region "
@@ -105,6 +108,13 @@ def add_region_panel(self, row):
             tooltip="Summarise each segmented region "
             "(e.g. center, volume etc.).",
         )
+        self.save_checkbox = add_checkbox(
+            region_layout,
+            self.save_default,
+            "Save segmentation",
+            row=2,
+            tooltip="Save the segmentation layers during analysis.",
+        )
 
         region_layout.setColumnMinimumWidth(1, COLUMN_WIDTH)
         self.region_panel.setLayout(region_layout)
@@ -204,6 +214,9 @@ def run_region_analysis(self, override=False):
                 choice = True  # for debugging
 
             if choice:
+                if self.save_checkbox.isChecked():
+                    self.parent.run_save()
+
                 print("Running region analysis")
 
                 if check_segmentation_in_correct_space(
diff --git a/brainreg_segment/segmentation_panels/tracks.py b/brainreg_segment/segmentation_panels/tracks.py
index 3213e52..ba1b9ae 100644
--- a/brainreg_segment/segmentation_panels/tracks.py
+++ b/brainreg_segment/segmentation_panels/tracks.py
@@ -9,6 +9,7 @@
     COLUMN_WIDTH,
     FIT_DEGREE_DEFAULT,
     POINT_SIZE,
+    SAVE_DEFAULT,
     SEGM_METHODS_PANEL_ALIGN,
     SPLINE_POINTS_DEFAULT,
     SPLINE_SIZE,
@@ -47,6 +48,7 @@ def __init__(
         spline_smoothing_default=SPLINE_SMOOTHING_DEFAULT,
         fit_degree_default=FIT_DEGREE_DEFAULT,
         summarise_track_default=SUMMARISE_TRACK_DEFAULT,
+        save_default=SAVE_DEFAULT,
     ):
         super(TrackSeg, self).__init__()
         self.parent = parent
@@ -62,6 +64,7 @@ def __init__(
         self.spline_size = spline_size  # Initialise
         self.spline_smoothing_default = spline_smoothing_default
         self.fit_degree_default = fit_degree_default
+        self.save_default = save_default
 
         # File formats
         self.track_file_extension = track_file_extension
@@ -78,7 +81,7 @@ def add_track_panel(self, row):
             "Add track",
             track_layout,
             self.add_track,
-            row=5,
+            row=6,
             column=0,
             tooltip="Create a new empty segmentation layer "
             "to manually annotate a new track.",
@@ -88,7 +91,7 @@ def add_track_panel(self, row):
             "Trace tracks",
             track_layout,
             self.run_track_analysis,
-            row=5,
+            row=6,
             column=1,
             tooltip="Join up the points using a spline fit "
             "and save the distribution of the track in "
@@ -99,7 +102,7 @@ def add_track_panel(self, row):
             "Add track from selected layer",
             track_layout,
             self.add_track_from_existing_layer,
-            row=6,
+            row=7,
             column=0,
             tooltip="Adds a track from a selected points layer (e.g. "
             "from another plugin). Make sure this track "
@@ -111,7 +114,7 @@ def add_track_panel(self, row):
             "Add surface points",
             track_layout,
             self.add_surface_points,
-            row=6,
+            row=7,
             column=1,
             tooltip="Add an additional first point at the surface of the "
             "brain. Selecting this option will add an additional "
@@ -128,14 +131,21 @@ def add_track_panel(self, row):
             "each part of the interpolated track "
             "(determined by the number of spline points). ",
         )
+        self.save_checkbox = add_checkbox(
+            track_layout,
+            self.save_default,
+            "Save tracing",
+            row=1,
+            tooltip="Save the traced layers during analysis.",
+        )
 
         self.fit_degree = add_int_box(
             track_layout,
             self.fit_degree_default,
-            1,
+            2,
             5,
             "Fit degree",
-            row=1,
+            row=2,
             tooltip="Degree of polynomial to fit to the track.",
         )
 
@@ -146,7 +156,7 @@ def add_track_panel(self, row):
             1,
             "Spline smoothing",
             0.1,
-            row=2,
+            row=3,
             tooltip="How closely or not to fit the points "
             "(lower numbers fit more closely, for "
             "a less smooth interpolation).",
@@ -158,7 +168,7 @@ def add_track_panel(self, row):
             1,
             10000,
             "Spline points",
-            row=3,
+            row=4,
             tooltip="How many points are sampled from the "
             "interpolation (used for the summary).",
         )
@@ -294,6 +304,9 @@ def run_track_analysis(self, override=False):
                 choice = True  # for debugging
 
             if choice:
+                if self.save_checkbox.isChecked():
+                    self.parent.run_save()
+
                 print("Running track analysis")
                 self.splines, self.spline_names = track_analysis(
                     self.parent.viewer,
diff --git a/tests/conftest.py b/tests/conftest.py
index 136c86e..be98ad2 100644
--- a/tests/conftest.py
+++ b/tests/conftest.py
@@ -27,8 +27,15 @@ def segmentation_widget(make_napari_viewer):
     return widget
 
 
-@pytest.fixture
+@pytest.fixture(scope="function")
 def segmentation_widget_with_data_atlas_space(tmp_path, segmentation_widget):
+    """
+    Fixture to load a brainreg directory into the segmentation widget.
+    Data is copied to tmpdir so that when it's loaded, so all the paths
+    are set correctly.
+    The manual segmentation data is then deleted so that saving/export
+    can be properly tested
+    """
     tmp_input_dir = tmp_path / "brainreg_output"
     shutil.copytree(brainreg_dir, tmp_input_dir)
     segmentation_widget.standard_space = True
@@ -37,4 +44,6 @@ def segmentation_widget_with_data_atlas_space(tmp_path, segmentation_widget):
     )
     segmentation_widget.directory = Path(tmp_input_dir)
     segmentation_widget.load_brainreg_directory()
+    # delete manual segmentation data to ensure it's saved correctly in tests
+    shutil.rmtree(segmentation_widget.paths.main_directory)
     return segmentation_widget
diff --git a/tests/data/brainreg_output/manual_segmentation/standard_space/regions/summary.csv b/tests/data/brainreg_output/manual_segmentation/standard_space/regions/summary.csv
index 3eba61f..0ff5030 100644
--- a/tests/data/brainreg_output/manual_segmentation/standard_space/regions/summary.csv
+++ b/tests/data/brainreg_output/manual_segmentation/standard_space/regions/summary.csv
@@ -1,2 +1,2 @@
 region,volume_mm3,axis_0_min_um,axis_1_min_um,axis_2_min_um,axis_0_max_um,axis_1_max_um,axis_2_max_um,axis_0_center_um,axis_1_center_um,axis_2_center_um
-test_region,1.95,123,29,46,142,65,80,132,46,62
+test_region,1.95,6150.0,1450.0,2300.0,7100.0,3250.0,4000.0,6600.0,2329.8173076923076,3141.2339743589746
diff --git a/tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.h5 b/tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.h5
deleted file mode 100644
index 2d51cd231ff8a4c0f041f609efb574e189edf2a5..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 39240
zcmeFZ2{@JizW+^#N>Zpam?<fxL>lywA}T{96f%WKLP{x9Dk+p{Ev#j&WylZ-WoV*-
ziZYbugp{OER4VDY)o1PRy7qsc^PKBk&viZLdCp$f-o4x}Z=cWi{;uymtULW8*RD1g
zKX%$!zLA%R2%jLI_@6JifBeSRTKZqZ-R6FO8-Br$Jci63_4f;YzR`Rm-<yhF9P7&o
z{1vU}EjO-N!$&2#^|{1v^Hq-i`;`b+{G0!)W?=0aJ+qM=w#L2YuF~7M&bX7|vi<K@
z-1WGY`;Gr5C*pej=XEX0C-<k~$od{RYAlT`*Ha6aPl6iQNN;Lbe?}(CH=8FX!l%Dp
z&%y$8LKL~<G~8&qX=DdUUWBSE!e_N+qlJ;F$)CC-hYUBbp}v6H8p$vDLw@QA8_Ap5
zSn91>zsADBz;xpVJxhoG=;s%zYv9jK{8RIvk5B4U@aOdo8~)!NpM2x~=DBtLZvFeu
z|1C4{&&TKARmsQaj^}OeJUH(2@EcjVBOF<QPF}m5Jr_6xIC|}I^&{IJ{;oSmHvh5S
z<wiFC#mmCPbTvPf!sW_G9g2Tm+~Z_6Zm<2Par}!?|Em@Ey1Yc{?=*~@i$-4mbIx;r
z&v*RPh(^wHK91flj(!gQfj+K(#yN7{<vpKn_jTMk@*Pn=?$1r&{kZ^FUq6ps-ha>6
zBTIw3UKXfpELLC0`*UOd{G6+|^DY+;Z?`|qs6YO>-zH0gc_V9{oWF(t{+zGxE?;UT
zkK}p#`1$*KIQ#SV5&E-d-Yy<HM|y};3;OSxhLO#GjE{S{Bl{)xr_Z0`iaLyZJzV_*
zsr*PkVV<A!b#>e0<>*V*neew3?lAX#Un@D%g^#cAFaGhT09Ab?Dnb2!ghIGIYSgGd
z6-4mPPJn-{=-(g08@s<7e3Ni5M#8-orVeC&My{Lx@%2z|<3>Fj>SsuPB>f*>4;QLy
z{^|H%U9UJizQ5-WL8|iKSKk=!cm6f0Y5ns$hR4ajPVwPYpN}tF`aeqi_1x=a*_8i0
z_dor+`vUh<|NEal%72gV-!t&<8Tj`M{CfueJp=!qf&V9GVC4D3$cvx<&ut+e_b!l2
zav%5izSv;Bo|Wmw|MNcDd(Tdvk^5KjTw~my`(l4zM{hr$U4E`Ef7<Yk+^+uZ^!Lu~
z@7G+%d*3*6$11|->EZq7C;!!b_{d#78Q<T$_ur2Xf9`un?zy@3P1OJP=iWEF<9(ai
zUrp6n@aHRjK5}2K@n5<x7yZ+t0UQ1Ee7?hp<WhioE;SPU=g$fM$^H3^ze5CB9O@tZ
zPrW}U@5lf9>xK7ol!v#!tDCFue|x={;pYKsd=toe5&5&uz5hHA@%umCSMK^F9Q)_}
zKlYFHGW$<O{=Xl~_{ci>+vuOS1*qe7<VPpUt;Ekhz^6HmPxCjfv;TM<9b5LdVx`?*
z{{8dc>*()1aGw|d^K&P1K9l-$KEmf-f1Y#koF4v@&+UKD54ZpEe*fp^+yDCK_I3ZC
z`rQ6M>AdjI&vX7e*Ad=h=>HF2bI(uQ=TzM5KR$;2^Nem}Ke>-%N6xJyFYfaf?&Fug
ze<8`o{c8p4-1zr39w~;LfaE@xnTXqykR#8nMqc=Mi;w%*X7XRn$9?R=eGbEYTr&;7
zKOKn|PrjMBEraA<OyzJ}9;tv-L@FVbkt)bp$iHWAK0a04o`Y0F&PC2c;%x}ueBAzf
z*TKj4_iuIh_^3yie_o4_f1i%=@h!$}Eu=Pb2~r2S6p1H$zU8>R0;!8!i3Fq`aure^
zxf;0!X@E3Du0<Lljgjk+>yaCfCP-5x_x{Tqx2appKQ9ZUC2|we3b`4{y&l`(wk>iC
zax0R1(X+>G2jn)SBhm@!jQn%+Nqs%?OnT(C{V(|3knTti<PM}K(hK?L;V|{JH*W7j
z`XF~BeUW}hf8-uy0CF#KA2JZheUx_qw}X%ek%y4M$Pi>GG7K4xJdBJ$9zh;O9z#YV
zk0YayCy*zR(a0ENEb<icG%^kuk4!)&B9o9bBpsQIOhKk18Av8F4aq{XksM??G6R{3
z%tD?)W+QWuxyU@^S!6!40C^5sh&+!hLKY)SkQa~_k)_BoWI3_|c?o$Lc?EeDc@0^K
ztU_K#-ay_&-a=L*Ymm2*caV3Hwa7Z;J>-34J+cAWh-^YWKt4n^BOf6jBcC9jBA+3j
zBVQn2B3qEJkgdozWIOUTvIF@B`4-uUe245pzDIT=KOjFMdyu`zKIA9lXJkKe0Qm*^
z75NQ0i2RQHf&7UaLjFPyBYz|Re3+oV#*hDe_%Vb}fd6km5GjNljTA<TAjcrbBF7;`
zk>inK$O%Ys<V55oqy$nDDTSPjoPwN+ltxZNPDjo_&P2)}Ws!17d87hT5vhb!MyeoZ
zA!j31k#mr0$hpXQNOk0V<N~Ayav^dNQWLotsfE-=E<x%bmm-%Tmm^mob&)HPfYd{-
zLh2)_zoh^3T7xt|8Y0&sjgZF3b;$L|4M-EDDbfsSj@*c}Kw2U<A+3;`k=95Xq%Cp_
zax2mfX^(V3ZbLdEosiB*7o;n4JJJp5j`Tq8Kzbs*kUNpy$X!Su<Zh%d(hup6+=C22
z?nUlH1|s((4<Lh(2a$)6!N?F~C^8Hgjy#NvKpsIJMIJ*&B99}ZkSCBQk<rK)WGwO&
z@-#9I8IMdrCL)uNG$b9Fj7&kMA{j^~G7ZT>vXLBQIx+*9iOfQtL1rU!kh#b_<XL1s
zvH*DwS%^H3EJ79|OOO|k7m=mNGGsZj0(l8}8F>YH6?qL=iL63iN8Ui*MBYMHBWsYi
zk#~@Hk+sM=<UQnlWIeJ0*@$dHK0rQ1HX|P)A0wY2pCX?jpCex&Um{zOuaK?CHe@^U
zHL?Tw2Kg4*iF}9bLcT|KBR?QNB72a%$UfvJ<Y#0*asc@S`4#yMIf(p@{DJ(5976s=
z4kLdf|Ngckb^Rap_dGueDS#A23L!@$g^?o2G03sVaY#|*c%&F|0#Y0~5jhDdfs{l_
zAtxiJAg3awk<*aVku#7pkupeGq#RNnsen{ODj}7TD#%&L*+^C79Hbg@E^;1H9XTJl
z0I7joh+Kr!L@q{ZA+?c9kUGev$Ysdo$Q4Lk<VqwU^^mKO`pDJDHAn-bA#yF!2x*L5
zhg^@`fHXmxBF&KI$c;z~q$P3_(h9j5X^pf&+9J0gw<7J3_DBchHl!oc3F(Y<LAoNh
zBi)ehNDt%=q$kn~xfAJ)+=cW(?ne3|{gD31J;(s$UgSPxAaXzQ05S-95P1k0j0{1B
zBEyj3$iv78<Pqdi<S}F<@;EXIc>;M78I6oV#v)H4Pb1@y@yG;ZA~Fd{L(-AS$P{EM
zl7VC*(~vAA8_7YYBQubh$SmX;WHvGfnTyOro<-&(3y|lKg~;>BB4jbL1bG2@5m|~X
zLzW{eke85`kynsck=Kxw$SUM@<PGFa<Sk@1vIcn@c?WqHS&OVg-b3C;)*~B`jmReC
z1LQ+wGx8DgG4cuWDe@WeIr0VaC9(zi3fYQmL$)JdBRi09kZ+Nl$alyt<a=Z{@&ocC
zvIp6V>_dJ+en$2q2asQoUy<LCgUIj5AIP7`A>=RQF!DE&`xpNF0(kvL3LpiMLdelb
zVWbFh400@T98wfH9w~;LfD}hgL{36VASIDf$jQhl$f-za<TT`T<P79YqzqCP$^G|C
z^0=*lR75Hvm60mQS;*N)RpcC`8gedj9#S1SAGrXjfn11Ogw#YXMrt9okxP&|$fd|-
z$mPftNL}PgBp~&WtC0H0)yOqS1Ee8xEz$^Sj9iCYkKBMXL7F1XkmksZNDHJTaud=D
zxfyAVv_aY;w;;D7?U43J2jn)SBhm@!jC4V|BDW*mknTti<PM}K(hIp0>5bfl^g-@M
z`Xc?1{>VMZ0OVfeK4c(rKk@)F2zd~B2pNnFL53p3km1O~$Oz;S<Wb}?WF+!9G75PD
zc@i0oj6udCPa#hu<B;*l1Y{yI2}wiJk;%vuWGa$@WFpg$EF>GrL8c=!keSFV<QZf(
zG6$K9%tM|<<|7M`=a7ZS^T;A(F|q`C0eKNwiY!BxBP)=Xke88HkXMn{kd??P<aOi?
z<W1x)WHqt|c^i2Lc^6rWtV7;I-bdCW8<362CgcO;Lu51Z5%Mwe3GylO8S**u1@a}b
z1^Eiuiflu+BVQvskZ+K0k)6mtxLpnzYjcenCE}npJZSe_p>rTKayg~E^z8ERs;O~M
z(lYKul<+y==3^8GRf=ZDL0?<(r1fLY!Jo@9m6z4gSd^|12V*~66A~Rq>a=fH&yJrJ
z2fkYlxYml2I=#8&5l806!Rcs?g~!E+-y*=M-KP-;?XTx~traJAR(f*+J+<QCQTZi_
z(UVBu!mhX&`(<&^;;o{3LxSXM{g#d22ytNVSbjEGig?ag=g(`_P<48*F_%pteNE3k
z)Lm*E2Zj40et(@x@|{+;pVdvNI<}kY&rc)yaFMH`N|te8+rMATYX<R4uiZH*Z4(Eg
z9K9>bGQ=wf8t$23PmS~D!1ESa;%zkfZwNWZfpNvF9ZB+}Z%<EUz>r%URJLx|VX8>>
z_h*(_Y`<3=R3sjMHd>i@n(h7x-Mix;eremKD=MV^3t^Fg_5f;s=P<<lXOsG^XK$={
zejpBB+#2&hbq?_rhI>Ce2#Euw9lm2;sgXJdAAR^;7ZC>v0SZNN^N2Uv+_Sg(cpT`i
zb-iaepLq4#Yb+~c;^4jRku5(L5N|l^N_<6p9MH-dwiYZTeb4&J=v<&v`y~<ExJ{Gn
z*Y4~^dWFn5xLq)zVv-iAbH2LzY+gDw&UeX^Z)y`CQp-MYCOZz+%N*?V*CBqeBsxDm
zKMu~G<DaLzjQIJXn^&=l;-Kh<U+2T+=itxXB(<KGsaNYTO5@;(o$Qpux+Jevkr9`C
zIS!6!n_gZBWSm-xS4?SDaiIRz<Z-JW$y>TT-I-V&2XNcNJzAgaZ`R=(ALDD|VEOXy
zUCY*xIuN*~KCU4S4zXl<It@sErM~K<xaK$ru`qodx0d8fTbEos{fwHo;@hSxjY;1s
zQ{79aU&TSJ?4w2R*PVmWcwFy^7LbqYh=YJH*JI*0kUDRYH+IH##X(QpGiO~>QfGnR
zm+=WbabRZSUf5|y@*Ywb;}iSiKrW+nPwYnG&ji~apbb*%wn?I7i6xm&{=%}V<X_bB
zrGMMAbrZ>-&WLtnjEV=7O>&H*n~C4GCD)5391jZH-PWqx5TD=F+>kCB577_b?R;oU
zyobp1xHA*uf$w+k`2AZ+-=ZJk7xE^@!=V1{*|PQ|pL-^GapClMc)@?=>J107zmDU(
zBreFsgG`cigPSAqJ*qlp70U6j&`H&9oHOw|e_Fo3HYXm!s<gKkxRAbm!K2%6&8NmG
z5!-FHo%B6@VTyL0W;|3}Joxdm8>!=&ys@iECm!~VS??9^LEiV+xuJtkb>m^l4g>ea
zo+LkY$4-|v{dn-#E^BM{B7VAg+uW|T@sL{o;mICv()awUEb~tr;$iZs`6g06#1BXu
zYx=%19&Df0>Xqy!b-r|)oEv2w56{MoFWKlv@{O*yI>*_?!>B{LnSK6bf9b&&4oEr0
z!^VU>sptTb->KcTL)I-GlnT;>)b^43_EB>#sCrTFi^}#rHG#xO1nS#r?v4jtdt+DU
z1EkK6-@j~D2E>C)QBw185b+OYD3ur=h=-O?dNutJsiS&7*xf259wgOLmS~2Myu_sE
z0ZtL|5F+|c?S3f9XTSa3;(a_GVh`7{+``Fz4H=cO4#dPm<)_#rz6erBZR!0>$KvC`
zt0-KYeuVh-7sf3}q{oA}a_=~eW5n+%U-~PZ84q&VGRJR6lD=10D=QYI$HVAkyC{b!
zl3%tdfnAv$4`B0c!q*dIJY%{n_tode!`#){Q%^*bIxB|fooOkGht=io8S=5D&im8v
z=k=6Q@7I*aB^OVT{QOKy(c#PSaQo6dqjhm4KXZ}RDzWR-IvGuKeicvh^-FE)WooG7
z)G~CGUm_XLmWrdr^XuZ_gJ84DXd0>W{9IHIG{%E<KuZISPU<*H&kx!7C?3RwyZF>n
zNFCuX%iLX_$3u`5lwM9Hb>2WwRzPb!i15uH*uW%pCjBg58u^Br2W9@UmNb$-BY#aP
z^?f{8z2+C##U}Z(qD?O6da3gRJ*at@L;BjCHT-^kfI3byF3*q7Aa%}*I=p`Lof?nJ
z_(?Ofh_7yRRp=g$2TkpYpu980mws2i_*)<Wln#jb=;RQe66V2_9FqXm(VIS0=aRn1
zcKmpzDwY7p8=if#I7@uuIgaN_i3E7Y><VekC*yx%I@fH|)C5qszDM^uN9y|-IUe(z
znE(T}jmrlMiT6n!CmE)Y0I5QHJHm^I?}9I%X|oc5cF|2itc3LavGko%;oJn+S{Lg|
zyFl_4M+$OkG!nopQ&>}>l#FM&!lckwS_#mTHYGZ*jQHenFLS>wO90k1QGd+}QvW?X
zRT0-qfOmci-d?^$>Kjxx49qr20J*UZ9jmVpfBoo5$yMtTpl+F#Pu*4G_a!`y*kYCd
z#%=}2EGmf)?mV=1&!z;J(H*Y-tcv*V1BP2;Y!e_?d6d((8^k{~`&F6akO1y`{1xBc
zB;&Vt*%5HlB>}E4QStMvCi$2W>!4R239#_KlX+ha$sb*Ly6J~^0w~1qx$1w1<gEqT
z_D%Lr0Fkboiov_Y%M6R{S{RrBmTKz^g6fFxoA<WF<WK_mlxH~pzDMeJ#4a%Q3{QY3
zpZLdz)e|qFT&8~PSOR=rSg=;Gk<^KKXW*21G69x9ACNrKMEu^TSsm9;C%{xAd0*j&
z#J7B5UuaEA0KVskwjFCG{^hb!&BLh);M8l~D)N|IkD_u;t(wV6fTE5cpN>5t@9)p$
z(}tIwNr3Irg$W|hh&LW^72a|-0ouE7)E<3KeDpnC`~Bw=fLRe5Kl&xvueNIMQOOq*
zV04FO-{BV0cT>x?{_;x+aJqW-t5L0_&hLps8qX>dprSI{F{F*hJAAzS<5mK+F5uhu
ztDWqZd$nEs^x6an*t1_@e+TjA`cul5HzdH|2Pem`Z^(EaJbXRNzL`3HzX&e&?IiUD
z0#CjUewF|`dh!?_-x2>~Udwy-s|3(poSx$Tp47kpDQ(5gjs&<oTWkL7ZsO_JQ=h->
zN&u@_h1PZ-iT7bLZ;t6rfI~4`r#|i>_1(jF4bB}%fQ)NLyv_QEXY6_!VD_C_pOg1*
zx$}wCQE|*#=|4>6KZU;1>nHm))rz%|CXfjBg<U_&2Z%qp{?xh4B8kvaIFzsPmDEpM
zrEk?XJ`qlvkN%$hjpQx0LJUSrBtl1@`3u?aB;Tw&DMW2*BA9%fW1aYe<d;miG0J>q
zA}q;#;Uzjmd~Dy=`vD4xFm-Ryw6I@f{QWgK?Wwa8VSCj#$8W=={-HR(l{e-lg2AT)
z@>+a_|M(oLZ+^I=OCu2iXZ*JBqdfO>>z3b{jpEvgaAIPfM=n3{uTnUpmn=_&^TC;V
zE~7}D4yD!M4tj~8I$5%FrU3EbG)u#y28p0@c!up`9^W*;*m!PTBDC%6NRAaG`7Z7E
z8IR2p;Z*vm93vspx7|TNm(MB@{ECf^j2cby_jVM|RNayYRht)TT;lO}tZwRW+?EJ=
zD|hGm36p%6JYV4f*F^B0Hutld2=Rt@7dU0^NQAXX$3!}K{PzuVw)MLb!Nc|J$CNQ7
zFO^42`tF|y6Ev@|Eyfo9<8#o(SI*Ov_9uc*QLECpaU|c{wXWSHI1z5dEZ=*b^4!l&
z&uXUi1|Cj?H3z=Z0!7LG#_Uao%*aHTqojF!!FZCF{}%tKAvzHXb*8WD<nj5d9~KV9
zCBj(G#^($%($~;lp?)?!5q7x+Ew`LNyznX`jZMr%kZaM~IaZwX73!JvI5a&G)*mwT
zuB1G7oSU{8RTgAZ>+oCA%7BTaj#Sdc56|-xL3}W$O>Gj%Tg57E9aBvC1sg5fc>F?R
zJ2CCDMDRChOiGj>eQ#W!Bj9=^5q6H+aCV&}sUtmTwKVQ}BCP!R{uG}S>1**Wuehov
z5vIIrSbKr;+;tmqs^7b>E)lf$DK&abCUvfV`sOLyln70><J4rPkUDdP!*fj@CxXn%
zJgX)i-&k)y=g^BpuzRR#d1NZd*BNc<&TmVED|gN+E|n(rHT@QReEBvJg6`Yg`pDzu
zj4S4geMp1{r>u0?)5tip>n<1ReoBPF>$ig~r;|EugZu%zz9zz~iH~T)Gf3agI*UNo
z&qRn<%Q##?dG7mi&mib>BVQ5(R$1%soJoA1%bJw}qmw`)OuAWChV*491@bQ*mjtmh
zM3**Ep8Gti=y2O8H}NFU1}opgvc!vouhpkZCBdT9nTNFGNS%_I`>xkcOM)A(j@rGW
zJa;^G8<&O*%O-(hMTK~>d?5^v{D<Od%3+6`O^}9i5*&H<JZYT*@tKF27hUF1`9AUf
zA<A>dlW_2%I&FRuXyyhh6etpJI_6eKt!5H@DczvrphWt<Q(pSyw@wnI-DV9>P$vFR
z<1Z1-l}X?d#muYZ@d{&~1-P$Hg0OR5^LDF{`k(L7<QYauFt205Vfk6a^G_EX(`1qa
znNP$E8!6A-U)SC@bA>IFVEXcf=R;?c{aRuqa&oy%65PE~aa2Q<)H$uPSk>1)3A#la
z7PV2HTjyAW7Js&L5?p+ebRl*QndgNrk0e{%lfbn>dpxL-d}d{fv!r(tXrwnU?d9>m
zY0^E$eo0`$m!!*_OY#=oX2qfVlHf-od+LUHWIUx8Op416CV~6p>34qec!RGmKK6zs
z!NYl~X1VHQzb56~a+-B43HF`%ly5zsc%7!Ei4G@|AhM0#K6(Kef4#^3ABm@vK=*}V
zPYLC@`}=NJ*v$H*BzUYQ-{_=4#=~hZ3lwHh>m+zxl=woD7pq)00XRu;^vH)9S1He3
zC&r_jU+z1T1mpWuk9aO3-o4o5&H1xQu-Q|!W~wHs6LQN}x%+$)q@EFPy~X46myTm9
zl_r6>rBI{qVzOVOCbpZ{UrvI4*?GybTEx=~s<dfUNwBvfe_0*nx#KzZ)YrPHItd;%
zoxQM6oA|dim4%{pNw99JnTYZd;_pq^1BQ*%_?KQ;)JS>m{Ft7z5e|DqjlWZR$w3|B
zqcWBVU3s1aw)Ig`s!NGK_3^pRU@MhpFR6P-dG2_Ej*ZV;^fn1XMGdxwEF<~ten|$q
zx~cO+<<nbg%Zax;t1pw+mjp(h(qkU;_;8IEYMoz_!2HS_m9P~g|DI26m*USP@OE7^
zbDl2oxl@01I`Glp;Ut@$$CT%e=S6@-VyX}g%KfY8VJk`gl(b&p)3G$LTY6e{9*}(Z
z@cL7e#c9wt=~M7y%5&#uP>H|ALW%}O`2v?h^@x8MIm_YHG#WfE?`=?9Md~D~9ha<^
zr9sl$XI0IV=bo3$Mf!h<D$`)W>-Q&v^+_EinYq)~&7nb){hGz9t4SR>&)=Sr^J#Ed
zM7*qt$2)w~>Z{SDfl92;<e)WV{r}*ko*BKA2KpH%%#;mCeg5M=POVu<1FM0^?e{6q
zou5S3v~!17)8NGXXzRU(B!9c3diaJB4HSj0smQG*`6A0ZAp)i}*rc?o@iyhT^W)Xi
zt*LKG1DoY5talrc{S_%3oEC0F1Gk~(OVY-qPG(a1!W;H9xFE4`sEWsLN*_ER;6j5C
zjlI%4){(qMWAImf4;uUytd*Fwo_GscY0hD98n~~p=q{%`cRY8aq7U8jqrs2J#Vp4S
z#9x%mjT#+DgSYggg=0)eoq1)8?-(AU>dZ*v6jGjBN8kUaR%AG}ekynOZ8jzONrtx`
z-aSTx@mh;#{H8oNU!ZQs5sjunrrHFREHmOYRZpgw#L=MGOj&G$Iq@Z$vURaE8u(f=
z9`;k7yH3`bzMS=dL4#>;Qv;JXlD^`5PZdd})1bY?=oeU!yi~4`mvuG`M*G}x=%hS%
zJO|ziSf=FDptM?#ebSP6^Ep;KT8e0(<Nl;-(I!%-(672czKjN!1d6MkQl7g`E>z0M
zyI!He-nHLY!B!;iIYG55`#KHIzPan5yqV+&?aoDisG&jRaf{)*l;@6TNA7{h1@~xR
z^7dkok2Ud6KW{$o-$VnA>}k)XY)Bo)#bTq&AJf1iI79jp<+<}C>EnL%*9#gfn!bO@
zHd~S(&t9;2RXYu?&Z*KC+Cu6(Rxg-vw37y=c4<>`dAw@NLaDkBG??tP^ohw<Qs>}X
zKEsKhX|VoXX<$F)x#RbF>5^gfjRxvzD~6Knh?lLJJuh{L2A>vxa#(IpJTqo;TN^(e
z?r#~)Xr(-N{6dC<*JcURfv<93?O_Mfcjud=J3B?`P&SeFVD>gr=e_&2U&RyYu#_+8
zN*#|+s$a78$7DKOSZW;U<4E#Hrrhs<8Fbk0Gh18IiTDS)MoCBI=@4PIyrz`L)AarQ
z?yJzj;CP<S7H2Y^w~Eq%Qgi9>zAhztn8!bQ$)CMdgAN%r3#-{KBp<qK@OQcv9nPvI
zwyt(1`#bcgI;49U9aeC<Uv^TSyRMujY0p`tM~9QG)>n^jC*!P<l^EJ@K!+!B0%7yq
zh}Zw6^!55WI#_y6RBhn#J=t1Q$C=Y%+2rhEUw6{CRc?}#g%usH>wlD&@*sT|Ifp(-
z*+K_}^wIvMl;^&`Uv`>!b!?+U%4_Ldn;j&-dZqr{dE4o*wPgLxpOok3b(&Qr0zB!!
zxq0(C!;|>Vd3AEveCUvU<4opCFH*<TaKKU|fDUQfrtfT{Jhx7>tk*g719Z4;yhb{F
zC&}ww_%%Nzgbr33>vEO6N&cpg-=mHQI)vTWC|6B+?z#$?kjPPwqC-?+=mED~q>gXH
zg7|&0bm~6j*Ci1jQm58neBSj0I{euDp)H5UC;NPPE1E*hzmMNX<K3jrz`A0?%`7_X
zVwbgaQ=Z%R^CR(l>`bZ-t?j}IUy?sJ;Yq-UTsq{Bmh_wFN4)Mp(`xNPI^3e~nS7t}
z+<CYVR;v|ufe!ol=d!&0iPtEPT7R#C4*Vatici=>`c@^pjF?(Uhxl>oZs#b^U025@
zN3}cOqQef+TN$PSB!61R-s#+3YCSLw@AOih+xORyoY2pDIyCHl{xD`Q@h|NecMY29
zP^+0$Gk+iP#&(RN)6eM8n|P01PkC<NZC~GCZh1uqZ^OOLy8=mGeYRQGoHukh_T{$d
zg#E-f^(AWWeNP8Jvv-UF9`9%pdFDnh9o&U~PBu9}#{Wt8(-QG7bYR|5+4qt1+`f!T
zTKCo;bXae2s_H}#slT_;=vvNiIuy*^)~|Mu_>^7ui@ypc!;o6}_*%+y>&VgjZtIUt
zhMoOY;yVtJ{M7aPhhrxs!-|EhA(3E`&pfi&vPCi(Qj-E3&rqH_pAQT^+)<mB40TO@
zF$N(df9<%tL!fLj2*#CYz2Wg8X8w{jO3C1V)UhT!l+;m8eELpuPBJ(-u=EtdNSzlp
zPA{D1C&TwUiJ6r=K7;<Gr%*E)f&zB*+lP}nX%?g9e(NN|Oi4YpUzF$W@7stAyEd#$
z2FZ=EA>}a1N2wRSOj(@_>90=O=tPiw_PR=s_eRO^#naU6G3B}YYw`S?gpO%4m}TW^
z>^Va6XP&Qpe$+A<+E?8AHR&kv<;I&TAK4_so3@)Jg_P$WM=2_4*OeWTAz829%H$ZS
z@3XV7&EF*%=8BH(=;rZjW=hE1^hk!bzT~ZuktDCB*%%<PD;YYOXDU^WlfGLE2fCg7
zli}d4P2x8x&mDhFjN!qez+`wTQM=ADig;Q>+?-Lt$)Kq-&u^H=xB8^_n;%Yw>JIji
zloKRh>J3lWk;!o9mdTMNCyCd(t@-3rbTarx+4we7p1Z%h4jA|A#U}&bn@(fjXp%pC
zWS#11dNOcAWyL4N5FgRWzrUSH&5yXs^*kPbJXQP6!i;2?rYG;XHkRyfUi!ry;W^2$
z^}5!_4$5=;4*ys?_F+LXF!WDs4?RWdr?JXzD3v6`M1RA3GfxxGI9;E(r#u;YM@7q5
z@c1e9o5QQGQtLJ?!r3a0+=nY!PLG^&GZ~r(E#f{?o;%MEWjoHe-$@3Ot1{=K<B4Cp
zAm~;3{bZ0z8+T<+0`W27`SZs;Oona7)GZv3zcpnQZR=C&INjBr=9EbGcT~54bU{lp
z?7lK(?=a=LeQkbTKL7i5GUUD%T9KSY{96NUC-bgkSbM7J(_$L&S8A8fN$*L9!+X*b
z8hQM;qe)}F^e02wO6KgHbkesVa3Y`Kcj|p-zRML!CiC3O^dC<jPKGtBvS(yao;#kl
zFUK{y1yW$y(@}?3rVxKJNXBc$m=x%rBvb#4$4hB^zZ@f$0!*f?ct9%Y8#(cfV!K2N
z>`|MiH<3a53d}gm(v(ht5_d1Bvy|udO);6e=%`Ez#0m-c8#0M^JJbK{sbUJOIzKs}
zoyUhNiO}Y%rhvfod+rC*NZ+9wvO9v+Q()}pb?c|HNZ-%lg62(&QXuh7nOrgDxqS~h
zg{)WBNrC1Jowpm<#GkXX+qzFT1w>x;r@rIyY(LkKTKyDo(;vSkoJ0D0bDYX#j8Z@$
zG_Pf5I_azOPDaqzBn1u`4tSJNp4+!;m5cK&ixkLSGrn_U2JyEGhnlC_q(IA@8CD-D
z&pj@l&oA2IWuF4}Sq7y?Gf7{#YBF%mIR!Sl6^&KMB6S39Q`00osQH;0VRo6vui<yK
zbMsDt)PRm-R%b|lT9NO(OMX;+#jnMEl;`$M+7UBhLLfE%rdN%}vq|4)vm*qY52e7O
zMK$e8IV9ive&(3c@D#W@(D?Q$<+=O&y+BcJTx1HYt5$BY&Lw_~-?p_5(J3%(^Sj%h
zDbMY@WHB?eI4%WTtX7^mkw^Nb^=H)z(^J6K&n!^oEXkXlU9GZ}nF9KDdMhg_&+Y4B
zr+2C_Jq6a@(Ee_dPyDc;o2XDu3glVZoav`Lx33~QILW3U1=1^PjZPMjzTf6=T3Ap_
zjpt%XtI9c&KYs4~i&5n%p#Db3tCI5EzTam$ByPT%0`rGLKG+lzZxcAn;p~kR$jv?H
z*iU(G-)q+%YVqGmfrn|`wI|M#zH56I%(S|n0&T}_)l`Z|9ih>qC*?gzfl1R9Lay=n
zNuy+>_@1Uf(r>rB*2N?r*0xGzQ%efucbbp=OnGkKT^!vtxvx_oCQJyTN=SZS)$`rI
z-=)CL;A~H&3&bl_wdYy(q`>i$+fH8L@e<Qz`*ZqJVDeUrOsk6|-(q5G_<Jw~Op>#T
zdwG1St$nuTa0<A83obiWO7_=v`4_2NfmC2B2$jf}kvhho-^KhENrgGJTe2#6yvKfl
z$(zJdVHGPU#-g0$A2o?(=Srl4oVldW2g-BzH~ag;HGI;k@asgCVMGP#JM=-ecauyi
zj1zb~S>_U{)0$hJo~M`!(N6ncUgYt98eY!)s;O|%wlT%@GRfaL#IL?tJryn|^KX1d
zd2ZkDvHAk}i&CMbaRFcG71CFB`%wERom5zV?r@g$RZ?gD(3d)E-Bd7YS6EZT<5SFk
zR2Jx`g1|J57sl5}Ucu$TRY4=FuTZn+Ysz!`e$!Z2ZEKPWS(cgu2P#Ql?Q*rpg%+tW
zc)@><WEH9N>6`AS(Ke~@`N^f;vpjzF2NOv<`&78?HQ#y7b&_8duMEY`sZb>MrT!)5
zxqXFtPWX@UNQDpS3$+4nkbK5Blk*PVsc<H9e4^M*;<e8${(8YL6_lFudb4=^gGuA`
zMFUeIscoVD%3CBKD@RXqI+O~__TP?u%;S$e%IPl)PXz;$(@%C)lYHYs;ms2wQ=wP%
z;xv&O;zJf4sB(=?1wU7t%`6^Y@Th*_<+xNBE6j;ra+}Osz_j|@N%T~B$GTA2KzZ)|
zz8zYn>cLEf4ejM0-S3dRT+o;NYw4-rzW<>V|6S5oGqG^V<eXI4KUGbK&f}LG>E7`w
zNQF|**PAtJNu7;x?hZGKQ^8tM(f1DJxqTyizyFXfPX$`fgkz3%BtO_uk>YbT6|#fs
z>E9{OJuj8L3(%{%kqY;Mj-EMnkNEZj-@0esNrlBsyYs5|NuB7_(JB7-Q(=1S-qI?{
zbNk+1>upl|AQdY4f0o<Sll;{~gF*^VQ(+xveEBCHzkB(vs(me~ASmBl64^ld3WlVH
z)W1#znR$Ej<r_(zlZM9nDqX2?e~)x}Ipw+IpMGl9n1em3V3pgJVBSRXk8Dkw9`vU|
zTi}<7F3NNBH#tk0bG}o@Y21Fl&<Df|yWRH+9Zm(IGl}-n4@sSp%!ezU2rxj{d6(gN
z%5(R(W$ToQ^T#mYd%zrZqh^xNI=<~~gct+Hoc=hejq=>SQ|Il!@<M_E2c28K?0ZDK
z#VYm0MbZp7lC1Dp{4w!fw+(!b%P?Te+gBH}DbMX2aOLE>Hbn+3|4RLG&=ZnBEqia)
z5>*CR2X6I!!s8o1+!u^dXF!jPn~Bd;k~dPcYk#wd0h)=Ssv^(G_@`9qRW8?IfOB-C
z0E_b6@tfQ<V8-b(;Nr!gm)g&X7gW;<daqCAuV!DY=kc5RpW5gdF@R-xAja(l8ULo~
zLs~Qw2IS73=Kh=V-1^hm+a-D}7*N)g3rR1DzaI9wXN?U5_zf+j=eH1VcXj!lRC@+o
z>s|7>n#YT5T+ZotX29II&|3RfB>#BX*+63t1_ZTk%>2sZTLZ03+1}Lr9BB%UZY6o$
zu|kT2ehesCtg=O=jd<bi&;gS`28>Ph)4s~%&8H>Z$UMY=SstyDR_!FeI8iKiD4YQX
z>reN8<ned57TH)tGJr8q@F3zf$*Xv+m&=W2fa}Dq7iM-4?{q!uHD5df9+_XGUEuKv
zc2Vh@>D2zZ4j<g`hU6bEyXjECWI*x{UDr1}e)*Q^Q-v}Zu(c}7<ltM9e^qW<w<U)G
zdG1@5N_7(dXw&+jq5^83)52!w^Z2FnDrSu-VZc49RnlwTk-Utq_kD+Q229$2P4p#?
zXB=7RbMYzzst?Es?CB!;2Ej(r@i!SDd-wQ|=zHStUX;jizQcgcl1kq*czj0k&BYb>
z8E|IIr0>hSN#4I8ze@Ze1J>H?A9~2+t@>z&ZciEDT@)s;^8?9mu{W!^(n8h0X(lH0
zk$Bcu21s@=V5#odnG7DUQ@XBrM;8N3>mSZv+(YuZOhwtM9tPB@JTkb)<EwLe!>0^T
z>#Bc|lWQ-@cRdbl+xeXVAKplW{o?VJDkrpW4m03qUP@L%AIUq!_eD$-WP<INzWRBe
z$Z^^#v%Ah`3=@{`))~G<dG2wv`R3!XHDXM7c|Tdh{xk8dVLj9{6(;yzS-RsZkB_zv
z_V$xzg4PG!?C5@y->+{OcSnW^PGj$Ns|*l-Sx@_%tP&IE)r?(yjmLYLf3Di2%7li-
znTI!jA$hisX>F}K6T+UVJ?-W3K8&fg@|sKtIcckT>?_HCZk|%LSBDA1sr3nR--ur^
zbxQs{T_$Lk7XB*Z@jIPbqZL;(!7)nF%Y2aJ<BL{o4>V#z>Fo_|T|B;a^!J7JCQR5Z
zv&k;(JIT+ndoiG7$%I?W99pOSApUmo`Mmu$Oz5jWwWFBFUlP&Z++fdyS;r?0t@}yx
zGPbtglwFvRE-jqU!Q+Lps-g~fFo8KQVe!Esl8=19O0v<L31)lRUrPNV{+e=jtcpL?
zH%&6SfX83obcru0kO^g5%?%8PNj`I*flK2dCcLt}I=Yp|hn<kPq;i-EETho+eZNWG
zH|B$2P$Ux;J-No3$anr9>nAceMXxcM3HL6f?a$@$t+wxdRN|SCcgny{pP%HTWz$a|
zpff>m`m?n!czkYcOjZMv2^083R_qx?@*C{-pH<Faf^fgxGBE+-Pc+TV*q_4$1r~kP
z86F?3y8C2(0TYr-tt~*1<X3K+=&n@4gq(MkJD>9SSN@L|1(q`*y#3oLUm=n|<a)C2
z-c=^JSSnVI8%_M&1-}>yH<=Le(S0z3#|u8_TekNO6Q<ksEY%ez`I+k%-mSe)9Tyv)
z9(~N?Ptc>*$~|O4Y>??I9}$w@QC@w~|0xsNj|b?D9YcI$PQK)w7A8c^^)E>0@ekd{
zI?e20LbbTJ%8Iche_|kuxx0%AwQnPH9`ShbbsO&9>S02Y-A8SoaU{QQM@@(H023xF
zUu_&Cdj22x4=>GLecbt-dVjm+_oh>x`+U`M$o@^$FcaoG>dsy<p492zAzmvbm<D?_
zb>2MU@y*Fsm>y%&;DzLp0v|Dw7eBsW+ZC}iNIJhQeC!0`%d{3vm?)VBZRJl~GI;!Z
zW=x)obQ<)%xV2GNoa9R@CoC?LNdvyS8#X-Q@w&+e^TsKq!P*6hmcA27K1MZN%t4jP
zU%TlpI*IsdFJt>6^)zr){vDac<G(wuNEXse1I{GRGLRtomBI$st##61c(MAxGai55
zR{c?)ZW`q2dau|cN%Bs%mCeKYY0&Z9EPjF%@!o>fmF7li(0=*nU=EM}U^ONw!z2x!
z#KwEBo=oyRrk2*<EYd(wSnop%kB|H4BWP@s2CVxk2Lh*%eD^%nlvMjPc$MfmRdOov
z5)%w(^g5@(53$8J3wXR-O#D$jk2LDvRb?MCk|y~%XP>`_^G<`p+7E^uJl@FV?ew>P
zX>iF~MJ{+6$)CP(X!Vl7G&p?dfbg{G#Lr}UTOB`?1_rqkhf8?8;EWLS7vX8Jge@j$
zHiP7C=N(+E9+?ISlW5Z2Jbp>T{81s%X%L_PK=;T@k}n_SRN5Gq243==e)2NJ_g;T#
zsYp))+nIx9mwEhGv7OKSnQ0JhHdoAAmgL(yo*GuCr-AW~9KU`Ze|J?+>Xe)`xFMa;
z6)i{dH7fp3-3!uSftZxf9C_kDM=$<WUYrJZznTbF^LYEfA>nc5X;2&He$G*W<dr&3
zf8TO74cr}%y8q(w_V=zgp1qL<t$K|M>53$;u*fR**X=ZLe&8v#Sc&)ui7AUs?x(?`
zTR}e?cs$3e{u=W^8iZD>{_s*J`4h{H7W6(zqkcb%Hg=2(@eZLTM^?5_`I4TwnLK{J
zPEJ+K>oiFJZSJTyi{z!86x(0DOM_Wc-ko{L<7fRi(mB6}Iv#T$@$H{Y@&|u28iV`O
zV41tV*JM@V)qNJS>Ic(cIiKbGA|7vHEM_4)oCYEC>IY2ckUW3&FaJ`07U&=NCi8*E
zXG=RRbP!}g#pUMaV`?NnPnO;`T9^fL2|a1bbBU)n9}O!Q!-D)`w}aPt{EuZ4QdXiY
zNUGWH?J$qzZOrB9zr<J|9Wv2xh{xZQj2FwC$O7*ZC!^@<B)=xYV$XU>7TihUE7O`!
z{Cl?kz5dB8*lrg)_<+Y>@-Uj5BF%!J5T^~h7m$4Qb0f1=Ggx3yl6p~0gZQc6R1b8?
zuwd|vi{@D#Ur}*CF;<=h^<!JBjTVx8x$?|(Z6y|rI-%*^$>Zg;f26gju;BKyqjMq_
zk^Im{<C77pEZCa$`=^2?@zzEIJLb-1!B<9AdnJ$8UHg2=19cYM_33Z1Urh3syE?uG
zYOr9+{ffRJ9{*zQ*9<vL7EE+qF*QYt<nzP_R@G>+V4#DwNk^M_&G;J)J9VhOLt}Fu
z^Z4+KMkbQWSzx}EJt<%b$#VicFIDKWV0lsGX-OU8*$>7~b<|@)Z+_S8Vjh32PRd1i
zH4Czb-!_{sC3&Z9?rHf3)OxTuL+|79E27`uu`pu6&n;&DG0R9^aZ~8)@9S7Fwnlin
z`f}psOD29uGokYO{oeO@{HolL*K5pKAR4wL#%l%1$7vPc>9(NO)q-Eu<8+CiX&ucx
zWyJ!A;q8<2czn<JZ;skFERae`IcU6+<i$5PNxa;`g4zp}LtQ-n)Iy<>FnbnU{Pi;~
z5=ef3n*XZVjx3m@;HR&uNBomip=<Y?srpUp#@*rZrK${7zwIpeemSq(eHF>euaOIz
z>OsAar^-4+^odvM@w<J+lRAF+XkW8={9VO^Lr&f-m|nU|-e@(+%f6A55!p?R)790r
zi^tE`SUoS_j|Ik;nAaoMko;)jI@OH<)cD=>bmkZkf6kR5@ogV<98LARbC<_!t!KPX
zIluy$LpFY%h9sX=URXf=V_|q`I89r0E%D>>FFU*mVL_VGbIE)jzgs!$S7bN~N^ew&
zni!FMhy2E4>PJ|x`C-3IFOT2F&ivhUj0G2jlGeo<ll+YI2<N>~)OrwZJgu>g_#x%T
z#nYo%aQbM?KogIjtQy>REtUl#TCCmv>q$N!&q>-Pjs+8Lgp8TIfq2!k+to!9SrGUl
z`f@3cUpB{Q{#hCeK5O(Hw=p64i%Kpt%u-lT=6l!wCy&ocp7o`l!Ghl#^!GDNNnScC
ztvrdv0s+5}l$B=0Px(^lu`HbhL$PO`wDI_Y#~}17lLfMyTQm-vll*%P%@bkSEXegY
zKBKac_(!9BhL!WE^BZlV&K(|q(CF!wJNYcQ?`zZPX+iScf}1kD3R&>vykU`;CGqn#
zZ$F+;%!0W>qms|_`1Oybd@Z`bf?v}Mn3kJJe)Id;zc-b!;ArOV%5RkC-Y0L~)%tm`
zf;z97`~FV0B6({U#m1B?EXcaoyhV32@ou*t(RC|XkaxGLxsAuiCcj$WcAZ+!-?r?G
zuqOEzPZi#W-(ta)QFj$oZHTvth}f-in+3~KHNMvI_=%|>+Hc=w!4KcEP9Iy6pEz@c
zuID`#EYUjrO>zqvzpL=p@bL{Sn3GVUTuyoJbAze2>P6?8STJT~sIT2t;ulNO@0vGL
z>+|)mm;83b-_{JR>wnAwp9NO#Sv<a_qP8^Q84KJWotwDMp5#L_?M`UFWP!%v1<!hU
z{MG*H#!p|dfWPZkRlEbqr)6~a?Qf@!>yj0BbhZ(%QF=aP#v2yU<E8pqczp02i|<!D
zspD5CPwTKF$tV6Qv2}P)T{kx>F;tz151xNHgYP2?bYrp<>v_D<wGEHiy)1ae*nHRD
znata8p2_D`pQ-DYy4|U1E~HNOK>C;011yl}eX*~S$8Vx7Y&-mox(>z;g=}{vdGW+o
zCCWcoU|rdkBf6dVw;y7Bs)txGZhldJ5s!bBGF#H^Hw$)0#o61qk^DGLbeixeHdGwW
zd;goq`wwSL$`)jU!d+QfmOIHS7Fz9GE6j%Wv6J04co0APV%WKNW7x3hYw_j*9)Dtj
z_|s#eY%pFXyCY=>$#1#zx@-0XHaJHfWUlfgo?X`7c6%ZlF2)QDbn$r4(Vr{bCD{=8
zD`)2^FOt8W{`H8+6gIpIy&<BrllYRD=0(}kY|yUTUfs&$jqmu?7|virW>k6BF>jJT
zZF*6!Q-%$$yV`OV>>{2nAI3Q%&xRru!G<S1-c9#%m$DKYZfrA@4)r1L@7PSEDYsPE
z&{)?Hp}L#YIbk_L%|(?Blc!ZnHSzeL>pm+9%w@yfZ|m>u_a%A9<BtV7>TFPyief4H
z5r0=K<1T2hVXJ&rN*#~ad$K&>)gm?&)ZVDr<4^Jr-W(A=q{Rl)8)JrL_7ESX<#uR>
z4mF>p4A<K{KH$#vXBErXuzhppr`-W0&tM#yVWZ22bsG&h)AtfDS`@7R4XF2VbJpHl
zJibXN&on`w4R58ockkLq^52r)z#;=Sh^|<AN;;7ENH@tTjceI3A*j6h29H1OH>1gW
z9UD}_qE~wFC-W?KseQMo2^%7ZTkcOiK<X@AvHok085{gr4dFL<{HMVc4u%$Nm@y^N
z-aClYG5wf&wqp~U`u$sV2kC>v+e@3Ygj%!V@X9w4H!06Ozo`Xj4ajVv@}>`)b{!(V
zN^)7xB|A2p{3K^EJ(&1;?NDdEjSWj5_Pwg+@qw@BMt^Z)!}lx?mTw5j7ftC{blQ~-
zzo(u*Diccl>x1jAsJpY_*wYUQcX|92skL)zcd(%@X;p1N7|CyWadp4jPBsV#Jej2!
zPJBgEe1U)u8{F<T<<(Q3dp%n9t>9LgFB{Yt)|-P4ll+ghpRX<3!-kkSF>_QSi0`Y`
zXFT4^hHCj#<wrbT!{?Zt-+t=-U6{Ex{0OOEqo2t);UF7mV;kZY93}oL1Ru!_W&_K1
zm2eB?x$`W){b$dbFzR@0luwU3M*MKmuuj`yHb~0EyXiy{@1WxC7Ic&iAA7G_cJg@B
zYwCNa9B0F<uQfXoj+6ZJAZ3rD6KweL>S2z46!GKjwO4G2p~kay+@#Nx=g!Z7qfvk7
zDK>1Fekz@Hg5+sOAu24M4J&J->`hM+FMOY0Xl4={a;MAd5A*nSn!ejg>D2LlDZo5G
zn$+K29>6e9Wka)VQRLPbQs?m2flJ*?HV9{b?;aOR^06}-E=91}5WmrP-zCa(=Z!ts
zlPa6RhSD*s=6jwZ{(A3Xhw?0HJVDXo(@zusi|_a-iySs&*dClw%j2EY8jgI(W5d&p
znidD*NZ+43KXykHu;HlW(-O6KQm3J&T0{0c8@fzR>Am3bpS1O?%Zk~caX#$(i3E~&
zyQ1e{evu6+@0(t%NF?=PbH2g*GB()fp6TwPJa@m22zrQxUt&XXnu-c5iR2A$+FYD|
zl?@L9_na`Nk^J!l7u}LdHhk+Wml2?o{I<Q`S?h1G!HO>ZvV_Msxo_?1xW$IqD=KQX
zCzCqw-6bX;yiLtpSXIl^6yi$*_RW;4WrLXXCWTrauYbx%F#jGKPM=FYeK3{elS^Cg
zu5Mt1=a<s?3mC-TewXj_qKO(mz3F>9kAI?g<%?f48wA&je@bAIe6w4>&A2Djb@o!C
zlu;V-b;XMroM&uEPK<IJ;_<gfdtB3b$%Z{AB|FctNdA$7Z)L+PHheYN9Pi8~{`VCv
zmRma;Y>x--n!+J|qg;ypa0eS6_ptWW^7tnahkwR*vLRco|6Fi7$umDb4p4i~2BE$L
zDa{Nr|Akij8>&9ALAkpi=N*rqQfVe)-OC2OI%h9NCaKdt8)SMuu_12NGrNshq|Wr%
zd6Nzgu)+M7_95XjB!7S4*K6r-)cr#q=iX%=KiVoW?%a3kd^@yqy>B+jHxxaWTQkH4
z3)bQvvvP>n9$OssY?uwdR_8o^NqO#hz+emi!(IFwn44VL5tmE+*v!1|f`S~7HD${i
z=aD+mZxe^;qd72EVs|v(S(1M#ll?+{3<p#aHmQ|Tp1aNkHr%A&7{`ILhy4e4<`X~O
zr(uz`7zdhT`@bj^kUBy~&*y#+=fGq4LB$u8=hhjFT|XQm!GYdO*<o?#i0?abL|tk!
z2R7==nY6Bu_#c6*m*q_5fKDu@Y1DaAUtzPS(u(OEc(F_SW(AM;Q*7;On8|^?N5wjQ
zi%6Zpmq(*q<T$YM$cANV#l&BYFdQ?ezya^+)2d$c`1{g=&QZ!7>h}=n`%+6tehHmJ
zlbOwd)?1$)H(elJP(UR2{2UI%=@=gryGX_fR?<nU=W*aPU-jJ@%5&%8>#0i)kLPn>
z!rdt2&{E<TU3eg{V<86&*HjHHEhC=)nnC!GCI@`wPqYnCp8LM<EcpI1TAKqk<0pMD
zC?`HVtWiRKDF=ee_vw37kUAzyi<OF(bD-(Ln@Yt?#6Nu)F1}_Z2L{Il?{4Ao8s@Vf
zKhmSt)x;f}=$A=-$AXA`?yEVFak|dS@(S^Lx<dxO8*m_XWnJlntHcKi|2B>?;=tOP
zXsucvAHe_Sz>M`AXn&&9b>te!?|(AoXub&tvK0I4S5*>!Ix5IZ*PH{@R;}-cDbJnH
zYUx3>dJ7JO%FkMQsfy%P!gk$qu;PH|qK{YhUMKmA8#PP%tU171Y8{|?gZP-2Zz4jr
za3G-br1K}rbLaW|{erXN_8hqRt~jFLCh>C<h4We4sQLfc`@-`UsWWC#^N|J49Jr9_
z=Al|mymiKqPNgd~{%fGt$>Zg_oT^RTIk2dkrI1-e>ihA1RC~391B)N$n>ybnb*#FA
zcX;pQfY73oMe=uuf5PDl`?ZS$&+}c@wNjotexbRRz9)P+@Zqb?w6wd#cZNGGo4$wI
z-z%qO?Q2P$<2k}FbM|teuuyZ$%sP_)k)pAA$$kzr1iRjQ!Q=I9ujf<;aX|c;aX{)l
zlHdKB`Oq?$1G(N~t+w7LK3ePS<MvPvJo<3QYeqfs*&FMNd=7KquE4p97d-y1)<@Uj
zBOJKy^lBxef#iR+z3+>P<iJ2#$dFwl@n#nJ7Sbm;K)X=ZC(}gs>t3XEbXGJ6x`gv*
zwDS1S!;-8;r>OI!LwE@L0jU$PDm=0>jsreB9A`K`B)<H*vGImP4orRO)u+@<>RSq$
zwLhhCz{Iw8=q=^B^S^2IKJ)D<)N#$==;b^*|BwG}Lru+@q|XcvtTvx|&GRwIml)nG
zImF_?9_5XD)SnQq8c~!zHk|{L@9c5;M0xJ|Ul6?6FENt?j@{WON}iJZQ%#zrLN*7E
zFPqr2_Zi8nW!jv`&*i{<_lL)pJtz6a^&!ue<a5ATVXED4%5&@U6@4DQagGBKQ$~5-
zctP^?_3__Lia4<STWQhJm&A`_8sB|d!htor9W>XskUGwpKYUzDIbf|R(<lCl<j*=r
ze(Wiyt{1)1Z<{I4?c3M*d}ZKe4*b}-Z*od2@hiSOauc}5fuEudq1)PsPrK{jda{ZG
z68HL5RN9G`GvBgg%1!FJuxhLD2g-Bbm%5_`t*mMexXie!Q1qJQ3wG+-&$&ab=dmjT
z_jM5ekiW06xRwLO6@uedz9DrC^|PNWzt4dqFEyVCy(M`mjrGm98mRMz{Fl}`%5(d!
zy<3@O@_+;BTAs@Boy30_-@X2EGY7WTUe4U|j`;CE&(%0Q;ec|?1SiEV;xEPyOn&!_
z1IsibZM%7V|2@W1pO@5id7EiW@q3bQIv=}m@D&H<+&jt_)J^=BWktim?Huq|wk}xx
zfs8-dqmwE8h6CMG-_pi^B%ZNMN+PC{1GlX<S3aUVcfWiJ<+e<I&jJ1DMkQ7c@iS|B
z0vI2t>!A1kYL8yx3l6z#SL&sX->kivi~30Yg^#wY=6#~BBX4`Eeo>xVzo_onHH`rd
zRJ{<Iar+bT%N;XS%f51ev!uQF^k?FqozB~?`<*(Da&jW=`iZA+YYMpele+G=gfeCi
z5P$g2u#NFB2WpSR_V)AmtU4+22EKH-ZvHd)+7~j;aUMgoO#<l<)6;1Z^_BQ*jHH2Q
zLh0Zl8nkuuH&Xw;N1ehpk#sm{`!ZW`kmOg@1uMQ8n+{WsX=wIQp1Z$lHxs{lj86y0
z)y(0`--(Yk`<l`-Aste5jYN+BApW3FxP;%NboilcZNB*@@%)O4cHbn^p)9BGsnQVf
zN^TYV4opc0vA`vq&pdu$nfmtM(&?~Be^YMdFOm;pIV*?HNQZ5uUZ10fiLZ%Ud2zIC
zI@~{;vES}D@$HglrXH71r+#l-L4WT5)wJElG<8OC0B5j4rCM7kYAa9$Y|A^obQ77`
zS>v=b9SG_s3KLujNQMT*76MZu?^D!HZ~+SfXSE<WXD75orr8F{G$~3<N*qGUV7yf$
z3bLhy?0%W2Uw`@M-22>faz5S1#rXDI)__yWWkCqB6Kti!$=kcy;xC91`Ugv^AG~9H
zy9bJFe5rRZf@a(-s&%D0aTk~$sX`HQP4}k`yOE1~rHWM%2ul|R-bBAkb-3mA+S(|D
zBPy?Lp6^lKuWzK)cOx?Q5o-FMGXDSm`9k;fQoW3MZgN#qt2^a;Ut}h%*^F>bDZN$T
zL4K+ycc^g-!d~Bek60FYv1oEj;8p~;0;5i}gnZP|kp9C*2#<ub$7h_pHL*msPL42V
zNUQwalk$ImkiY)wHiVY`u*qsRd0ln7MZ5z+yz^*T#!_-yc4Um<Qw04l2^lfV=v<u-
zO&W<3A$f4&42MJcjZI@khFFA(y%%5EoIG*1?w}+d;md&-UHfv%NB`R6Y1qxoE#J4I
zf=fQ=Aukbsj_}-+*qO{DkH0@}(eFdx9~Xvi@FKsm_;5Vv0K$XWqEUBma%&s*VKTyY
z`|aUBo!sT|SET_d2zOZ%{x|t_E@|cryYUdy-+#lKN`0tKl4gC1C=H=9>*vEsE6DZp
ztfIPfgx&o56B|~NbElGWR%ar-(Rtrvts-y!%xhE4*UUV3L?0M-^1QX2-&Pz&=;ppJ
z>kv@>yh`R@mCbzT{SkXAh2(QSSqfebLc3gHR;?yCZK{QHY6Ny(OX`Qd<Rf<*-Pw5v
z+BVHDz90F?j*Qx}e1zh&1I2bH*R6cwa65@0?44Qc_NTmjq&K9nkl9z9$vbr-@<N|5
z$--%-kIP59)oaL8e;RqIDM6UdW$|{cC4X`AaoOuK1cf6Xg9FHAf}!an-y$6PL{;Dx
zNdEDd-FNCN!ZG*p%g>zL&*v{;+Ia-y^I%tV5aq1~{n)S;A@N$mjqlcx=jjCB?5{#-
znl8mGF*)Z|pzX;egq{@pkU~P<XlUWZ)ga89n^S|N<Z|W~-yfF|a?lp)5={P2#BHm*
z9^ph^*VaKNf3(QW>S{oUvLF23976Lp5Tdro5aC(#%gxoHR7YU!bNr}BFx{#OI2J~J
z<)Z7Hq?vi2JK1~U!pYT^*PRX5n7yaIZV8JZZ%r-VDY%Ys`b?JAGm>1P)L+%!K#;dt
XbrVj$uQJlI#K`oeoHuhX$_W1f$1_U;

diff --git a/tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.npy b/tests/data/brainreg_output/manual_segmentation/standard_space/tracks/test_track.npy
new file mode 100644
index 0000000000000000000000000000000000000000..887f32587569a5afb98820229879850c7d174e19
GIT binary patch
literal 24128
zcmeIYdo)$w_wb)6T~QRepeRvExkQSIa!X1fQk1Buq)17jB#Dr7I2Y#}=MpKQ=)yT$
zq!5ZoN*6^a6cvgRzxf^G+5em|8t>2Nd;9(IJpJ<;bFDe&nrp7P=HB*+;^4etlbght
z;4z_^-adi5k7%w^(_FpBNOOsr=AJ`Gf{yI+_dImO+vh*VosJy!!QsGty8?VrT(n$Q
zS9gh;p4JkzQ)>Uof0HD~;J-(OFz)}5l*j&MoSY~1#r4Q>aXnK08%aCTA4yW5oF_?+
z|Ba-cxSzl2|HpR!BM;=fxShB^ahbFm|9||5+mkX$ay?T18^!g=b;S8MJu*&leew0k
zd2;-3{>1f2eerQp7WYHylO**>nba2_|C>zellG+E-^Tw*kDM1L>5nA2j<_BvlX|2~
z>XQ_oC-we`q&=x8E|cRV$$66G_}@s@3mKm{Nk4zn6ZcE%i~AARBjvy8k$y;$dViD2
zadE$-OwN-e$N$EERZrZWjFTKEWzruhlO**>itCZ%q&~^N>HTe-^h1*LL-KES<UC1I
zkMu*594GZgM)6-oA&?~XNRoQuB<D$z^WuKUaq)G@agwATlB6G!;`XH8-|We8lB6G!
zq`f#veQ}a{;^Tjl#r4Hy(vBP_;}Rby=l|9HllG)P(*AEG^+}R?f0M~^lH%)-GC5C@
z94GyY%cLI3e^-yR6X!pzBd$-%<oe?LSN;7>kMu{bLsEQPT>hIsaeHw+aX+NKxE?8!
z^WyxQA5xFB|C^rp_}}Dzcbr`RpX|u^NPTgN+l$YO%m1psf6^1b*W&At^1m7%srOIo
z{%`A%@sR71`x2+P9x4BC>xr*N%72T8oEPW6yI*m8QYP2^@AgOf|Id<)=bz&Huj-NZ
zWSsvbi|hSUe55|9N6O+9UsrtopX|l;|8J51DR2L`<l%pMpJe_0zkFk+@w?IgbiV(e
z|Nocg|IRx5uio4Ls{OxO&;MKN{HyVZ>;3-~-{00H$4QFE`=6En$sajS^50$OpX|iv
z#pS=P_fLA_^M8}Y$N#q8fA=`KE=khAxJ-_Vlhh+g&Wp?cG%h~>H+xb~od2Xx&XW|k
zCuMT}U!}MmDU%fULyn7+)cf1G_&h29XZi20L;ig)z8>j^9RDZzpY%w3af;iK^Wx*A
z{O^+VNAhp>|1|$kcH;Kp^55+Lv*UmBFFr0V|97wV@A~`iis#>5_kT+NKf9mHHrH(H
z<SC%Ms4()r5(N^}2N&kerGTpOay`L(3bgLg$)svhp!Zr~)VU=TaF=9PA74%ZGk<=_
zenSe}=%gfgm{6cP&}8anpa6f;?(gf?QDEcs2Gdn`6!_?}?YquK3Peu0H$}sV0;`HE
zFDkoHU_prWQCSZPxamtg7`K}OJv}AC!+R<4J7ZH~-vJ8Hx1EsfK1_i((i8eygDLR-
z`hzvELn$D%KKA)p1O-kteja@sMFHEbJ;%x~P(bI=I;WBZ3Rr6vQH4|rTtA{>e49am
z!<l;KH#iiySzdlUOF)5$YF_A-ObQI%iq95gQ()QYtrpxo3Yh4W8?*0Hz)*OZ!X%LD
zMH}h&Dd4oX>0nAZ1v1V)_(H9w0EDcor98#`TWU?GJg30j_eOV;UQyuwMpIGJI|>|)
zmz+mwqd+@z^G8Z21#|-krY84Lz<0(SYDymku7rDp(g!IZaiicq^A81*F30a>kEOyP
z_iG2Z(p1Q<fBsZ3g$lkCUQ@12qe9-V$lF;nsc`KI)8NKzD(srtDSumy3O!oJn}iFf
zaPYTNSFtt~zJ`x)En7x~q}z%sAL~<LKXY?uy)hLMmz*AKGNpo&-$CydODfEEY+2H2
zLxr2~G915cphCa=`TC*FR2b-X$Q|oSg*u<IkCS&$ff0T?bcPrDU)kxWx|a$bvD&w_
z4p1Q~#B!VdVJiF@{o`g9OoeB2G;Z02QbDUGe4k4M738PC4)%(oLT1Nkv;PGu4BHlP
zLlda5cvX$)Je3Oew<qhSFsN{}aMd3HhYI)BE>O=CP{G7GnOBlYg^8U`$7{2xux4cG
zmF7GubR~5y?YT>ZW$8{+MhGg*)ibn^DWyW~Dz{pd3M#k@qw|+Nro!vo*dD0E{3tJr
z*xX13=>d~{-mj_94kt5$-&0}w@4{6v?NqojuULcCMTLSqXYbrzD(tv!JyhC{`IzGQ
zp>c=`<#j&uyGN*ScWA+#QORUTJ>$z!n3xRb{eHdDl1YY?T>gGD`DA#ramglUrDO=S
z-F1Hd+++w~PF6TGKN&h<;5%J28Gfz$q;X?OG6>J*T&-B146El)KmE=y8JLsbUi)d1
z3|-K)P|iFVtV*8tFS1I8(&&o{7VDGYcyd#e`=(?F{}^U-)HxZtL#?+bY)giwO@B(R
zdM1O1w_or>?_@~0<$Cg+Z!$!iCDs4(PX<@Xmg6&nl0jue=Fp0eWca$Y<Ca5MGPEV?
z+U<`_2KfaAI_J-$-6t!r^oz-GtSO_dG$|RH2Y9#N(~{w>$=K&3tYnD2ddXrgKN)_n
zo;6}}B^m7ceWcy5Cqw)bQ;!ojlcA;De=PG(GB`IE_Z5ngVRYK?@>eBD&o9KUvSip-
zzJofuIvE@-6$@5BO@?O;A2dClCxc|j#Sh`HFfU7+yLj)Cq2^oKnul%4z$j(C?&!q)
zaX1er_M&~{p}{5n$-q9;6ufCD85Zu!Fbf(<1|=8nX1ZhwT(rNGOH4?C6bnn2mZ>S=
zWEf#JUOok$FPI*&SSbaJWFL%moSOoPU$@r;&rgBAimR<G%@hcu1e!fqk^+{FXWKfL
zr-1utdY#<L6nK})A8)ie1==*$MtYj3KqTGCI@&4)>igN7bJwRp=gO|jFE*t>tEZ&t
z7?%`y;+H>H%RL1uT^?F)-ih<`#@q<?NdfCBU9WV%6d2l|wypM13e3o&(}sdl;Lf|<
zvoubmKyT;zR)_EucvUg0H{@&zBphQxdQ1w;`$=tl8lM8$Ya4D1Q_zpc{Mto~6p*&C
zUGB_DfpI<4Uz`@Cfb&k<lAKKROS|*=b#@B4Bx-D*n4bbG>~_ekEJ%R|rP`7{MJb>>
z=(37ZngUxsW#!+mNP(zJeg0n`r$FhqpZitoQ(%d7<TZyEDNtAIp&j-n1s>ZvcW1Yx
zfPSFvm$w}$pk^mrCi6K3G9MFprr%PaWl$pc(2o@O_;^bQ_c!j}>LgJ&h6Y*Q>&zs_
z(}1lJB4IF@26<;&#_pA)!GuT~OU4Wu+?sp5w0af|9&cV1F``NXKQ*_L%NNif{VeB>
zw>AxuPc`V!m(gJ3#Im+(eHwgPu<GTgF%81%{!Gv}qrv=1uY&ijrNQ?n^m%Ms8c58R
znpp2Z1E1dZCDKkbFtVL+akU!_k`@|h9r2*S`6xq)%-u9t;oYdvybsqYY2K~iPlHjq
za*tgg4YqV|%0GRa26t~e<QJTxfz^*UUwY5bV0u=Wm*#mI7}wWN^SnrdTHl&qDM=X5
zp)+%9X*A%s7KBK%X&_T#A`5&P*o}Yw_V^VV?8|p(zj>Vo5oJdfbl;>w_NvJojk`1$
zJFV5hgP;N3>i$aleT=K=h--a04LT}HZcKShgD0{_!MYCjo$DrjrjZ7*myO01y{5sX
zFRzUU-_zja`Kjr89W>YzNV7WBO~d|Vsgm`T29e9(EdDq^gRU`)52^pA!8N5{A3ev=
z;Z3l73Tqr4E-9ObG)<yI_lW(anR0X()R1W2GJ_5cF`Ay2X3-&K@dbriRXR9M==(Ee
z0Ub0{mdswSO$W98qx)l)(ZTbJ;nxa%I{f0Zu8dzr2Z;;6F0C^|JMOI9vuo+_sI7eD
zfh`^KoGK$EH`3uiPq)6M6CEbp7@QU9Mu&7-;fe<ybZ{07hDv(VVR>En4@+M<Y*AI_
zpY^9hRcw4%Ss)$mC|rvjA3}$v1qPMYVRWzz`fV5!Nr&myWerv5>5!wdgFhvn4*FUb
z(;O&txTUZ9IFU|=^&8)6K4a5C?%1Wg83H=YxHj&9YbG5+J0zV_vvFVAE`H5<bhr^%
zo-?n24!28G=J^!SVZW@_{mfF#$Bx|i&lPm=yA~Rw`-Bcn-?!cjs;9#)1DWxKFX*5U
zurK=08#+wo=^9wH;&{2PZ1hJuXf>~x@%S?xPIZj2p7xy%0dKu8yZoes>itDaS;KVT
zB+P7Skzl|$>%n48X$Bl@r&S)D!ho)Q#eedrF@RQUwEfpi2Aor@X@fZocxGx#i=M{-
zx@Ayc%|Zrrs4IueSi*qDw-TA#mown>p}|psAp?dR%Oko?7~nO1rM8|q11ePIhfi2B
z!19mAkJ9xF=*cuxn(WAc_iEF<wzx1rZ};hEsqPFIIOuSoZ6^cPnprMcx`zSA&$Z-(
z_cP$*(#LAW0Swr0-S0FZm;ujK59B(AGQj-ypJk~L40t#$u(2(Q0l~5<e4SVZ9C@ls
zJ)X#bYOASPrO6Cn1`KpeNo7F)`h4rHTn2n^lzYNUXF#s(%izz~7+~68zjoyf2HZYV
zZFuH31C|%Y+CLF8;C}MkD5VkxoSzxny1R@4NyZ*txm67KGpsi5cP#_nda$dko-@F1
z2Rkq66$5tu=pvfmF+gY0rcUj42GoB|S$V9B0lOSxt~}^vz&F!xMsfoTV3*LlxBp_m
zu{zu0tD_7^JyetbLy8G%jIWO@Cow_E{=`V4EcV!Q22RZiOn7GXxpMI=CKO-T=NY2P
zgosT+3Kbem&`N3gq@c}&H*Q{Uyp}QHGUH{>4Sgo0Z5gjNYRrU=sN07&m@y$%=f>OA
zwM@udU*g+k%Y-ejbLEXUGU3W=o7d4!OenN-E@^aQLh8b_(gmJO2srq)CCHlz>*mU<
zJoIIPp?+Pk{2?Z6ztP|A6~u&!1nX0IAxxk@rD%^0XTtjcL7(GUCfxOhm%JF%f29BN
zTRamEh1&I*Q<-q~+<L7f1{0z_N}OorFyR7gW&4WDOvn=4@r%65gppr1lNxd`Zf?N?
zjay7$C)K7LEo6dX?bgWhdrZ)9njNF~kO>W{ZCCbHGQmdo<h#NeCL}~=FPq%Jgfc}7
z`t~L!jO~_G&S_?X+i&9L=m*@t)k-hNP9_vpuhtgyFyUu3ZQ?*56LLbPjj<kLLN6z1
zDr1BREdv*pcT1+i#gi%r)=W%=8Ow)=1esKbnbxk}DxV6HbFOpsl~W<rK}j!KB^CBw
z>3h|ro(iWmCeW8?rNZz7RljhZRQMh#<yfzm3a7jx>=ziPLhg<o+mD;30x!@fqS`VQ
z{P|u5Dz>Sx@RN$v5r<UhwD8_twj~uF-57qZ<dzDrm*u%1@JI!?=g?ogI~An1vC`yy
zQ(>Ki&yKzRsW2<=$Wl>YD$x8AWMxBA;nkuj$=zY8(Azb4;@!wp7|Kysn;er069aeJ
zc*duKkpbs)J|z{>`ZB7e8L6=8`r3utIjP{g!~D_>K`M+$?VU32YAUGn!f9?fsqk}|
z^NO7ORJf3P^Uc^o^f%WtW$V3Ea0}7&y#62+ii;ZzB`Q<FrdM~aOHC?_8$W(}*0WTY
zK2~MMn5I-X(_*pExj7XM@lDrV`;ZF#7OxMFeoBS;_gk}^dQw4twrbzizEo&-*k(OC
zm<p_I9+^%fso=AEV8%5`7TA6NlrTDh1->UFXFAKUps`*i^O`&hGKU)0jZtR7@MqOt
zXB8H_zg8&7QfEQn@&mgiv{+CXCZV%chXq54mXg_eEKuDu+&b2n1@^-Ok6le!km~rW
zB*&5k^Sz(lld@&O&*kgOw>hw2W#Qr%H@2|g+fUPP<K0-W_OpY+4i6SsP!~Y{ZWg}3
z#05?CWx?g24mUmhQ9sw<$L&BCJexk*a!Lpb64e+quP_#zZ@k)95XpkN71CQ}V^|QO
z7kb|(o&`CLwJQh;3xa3WUYpKfLC3I#mM@0|Y9|M>iv=vuPq|?<^C}Bg+nQD$$YH@|
z%0ABr`7DS!GW1KikOhNn!qk9!EGYOCY+mtz1+Ae!yXRK2Af7coHK>LKRW|+()z4UP
zFKlq!{3aG0UVZ5P@n#lCMX(}jKd|7UwVma{P8Q6*WHIqn59V{*#`<S{EO^LY#nB#O
zf%M**1E)t=a4WrG%}Yr(G;dd&wrnCBGA``wh>~H0($W(}Z{*pq#X=`VPnit|9$5uk
zP+>!G(Z%&G>TF0-eWGck#fG}iHBt#WY_L_|(b}%Z2F3g4C993u;PB9jLosE8)MCSv
zU6yQc?rU;0w?)n#Hl#bSVK#4%eD4-El%Dy}v(AkTJMOKnWO=Y*!sgZdzTIqiVRAjh
z&X*1ELtHlU{MleSX`K3CARDB+Pxm{7u)%Q_?Lm4N8w^rp5{Dz%P`N$NZA%Or79Oxy
zyB5y|uhNVUW2kIM-VrO<%3wo-$wSW^4jWdx&sG|DnGNp>w>@^h%7#W~`;&P&Y*@8a
zs5$W#8&pjFYdj0tz&^Jk;LbfZJe8`OGW8)Fs_w{N-d)Lt%gl!c!WuRN<z(NNYhZ)R
zSGx7yCN{WxY^o?~W&`&J8x-2uV0?FdzF#NC%~4mq-@}ITtG$tm{cPxb5z=yShz%tR
z<5oT#VS^|2Yt$?$4!l=6{5W7D2ev+ZF{wfZ<ug4nSAhfZuHpWH${aXh=|rhi;Q(R2
zGFMfD0|y+x-#@Cw0bany>S`Sh1hAE==jn6c`*X$n$Ba3!MoBUIi75vTdcBWVU&{e%
z-kLojwj8J)KeVFOfdkq*-t=oYao|z@&8!n{91u;g-1^jm11j6L{aoP9fpHo4V?%v8
z5YuI$Q0LEq;su$B3xhaN>MR;_GK2%^$F02U!Z>gyPF}R|EC-Gxe3m*HgZkmA=5_HH
z|F*7xg;b1RRge_Q;J~%&gp8*g4s6(Q;`)Ni9Ehu5a{0tn4uoqRzF3>X0kfUc_G#SW
zz`LLqdLe}z2=lr4^~pUBEYbbLp8t>oQqR7c9INEOT>BMO)ir3(h_z8`;DBFY`JJF9
z4%A)ER;X;|K*PRCUUS=UzCVL=xRV2N>o-)D^>E<0Uqzd8KL_4!diUkv5C@!wJ=*V&
zaG=jmu6l-48a!tfaeODHf!qb%okcQfP++JhJ53=Cx;P10-pXlkw?x*UKqU=KnWwTQ
zYotMfmO#c+D-FsI*m>mXq`}axHOz7PX<#O8RqAS-1~<3qG-jEmf!DT@=OdPB;4f8P
z;%J)&pIs)>E<2<_&qgQLpIg#^@+CmZ&MghPmujW5Jkmg3)_P9w?le$Kphlbfrh)eL
z*Ka8PX<+N!F}ou$4Zo9wtu+csgSi}kmkVKOAT=k#@l9kJY+HKDKqn>*jGitV8xfxd
zg=4+$)=|=+xAd)(IwK8!$Q*nX#7TpewkGRxK^jz9^sp4Kra{HRpr(B}Y4APX;7?(G
z8hC^ZOHVFL1C6$%A@_SIKPj!d{vZwJ>$xQUd6WiG@`?ryHEF;<P+QD-mIiUFZFPHJ
zrh(^ahjV7lX|T-r>ivWdX)wb}z4hIvG|>EY?&GqaG&nWNt_$x=1D${s+}go3*kUKM
zMRg<%B+CB;-Id@%rn1)xPbn@eRb{k{m*xWFdGwi_iCp+)B0s}r3K!fQ)fj(dxbRSc
zDw95q3+tk`9JN>ALKn~SN#9H^h<t5lrYdux{I%^Si#c3S`l%J#slo-JuqGvO9v2E1
z%oVKE;DXx9Uun$?xiI|2?qZY{7hG-o{gy1@f~n1Gqh~r?AUfNAge>R6n<YOktLbwg
zORnFd!jKF4)CW%w8gn7?j;({jYA*1D_ltz4T;M&Ms<g|T3+Ert@RnZ71+Sg^(z30%
zuqM8%a*Hh&)SXX%7+TMTdAFv0PIKUb#kGhJ>m0d2G0d;*-ok~suJN2pE?nrD`g7Mx
zH!jF+t(Sl6&V@@W#@#yO!G#lPK^BX4ap5IvT=A3LT)3;n)(YH%{?@6TQS#-2Mo-|Q
zqWxT`Qa?59<<AAxb~Tj=hq<t9`P!w~fm|3b{a9;rFc)%t*z!M*bHUJ_-IW^3g^a=*
zIS|H$$=QTwdjuB_FUa^46Ul|Pd7IDcoa4eFUixVLc`nQmMD0Eni}_mcB7b%~7iMaO
ze=SbrLU5LsvNwebTG#jJOi1B^iFUxU>vS&M(eR$LDHZ)QH3$0GTu_cpE2MF`pg3^D
z*F?aDoEI?eT{@oIa{ah7nOtZb8!)mUiwmc0o7^k2xzH`;eEGl)tONZ=FJ$w%VBe_t
zBkwjBwstKVbuK_2djEY;i09_5UYA<Lg}sj&8Kxy%NMydVZz;udYntB`S%&oz6?J&w
zBQ8v1eQteJg`B&_bpI1B_(*S$oca{&W1CK1Zao(a6c|;G&$;kO|McU&7hGr-R^Ls2
z#f8)9JK|Qp<-)zb*>+9uxR7?a`|F8TE-b%$Gh$9V7s9o-3>AIkLdu_8+dR8aUPfFV
z^Mwnx#T#Dmdbu!?<?!9&JJxl-VBkYP7o6!UTF(B&^G$hoYr!urbbL+-C?DqHcbr8E
zdq?pc`0;7dV|frXk}*9~iU;a9Tn^eu^I*}PGr6B8@*qxDuIc;~9=MA>bS{$R0k5F7
z<<T@A>>vNVWS;^LN(Dd8O;F^)t%T=#naVt{s|>8Lp2LHrjd#pHs_=j@?+pLkJRa;d
zB{~;q@L*E9tzy|i^y?(M*jtMSzt?})kX*upT^C+R@pX7mC3Cq7mh-^;<?7&fdOUd6
zeoi{fkO#4bd&A}!^B}mg^p((r2W)oO95+)QOvwwk{0Tg$AD6W;*^&o8Y3-nA#e<l=
z3W{|$7)RiodjIu2P!qKtp5nlRRcq{jT;0TjZ-Z+*t+()C;}_QTHfJ7qC_6NtcICn1
zlZO2&?mW1!VA3Pn!Goo(@Yr=H4?;dHi68LdL7rlpUXl+FKAc)#q_d9)iZg7qtNk!P
zZ|<GgdyoghdTNej01qsF=9Hx!;lVGP@5LrZdGNK|lKtWs50(d<-w}9%2Njm-V`NY9
z;2wmZyB5xad3%2MtUJSlk-QP(ma{zYT<E{|WHiQCtaD890uNrFT;Q7*$AhUJE7xpD
z;6eW`+rE#NctDSa*fUffNY_Y=n@i)tdi^6i?=pCBxG<RQ$l}4-?fs%|4i6INJQYRp
z@cagQSgM!t+<FZ>3p04|XUv7MTdwlJYw^=_pR@2h4u0v5%HhG&it^>EH!)8u&0P!f
z@%)O+Jst1x0JdgZbrta7oRqIjq=*OB58dv}zQ=<L$=;^7N_epIn16=-13Zs?+k4u|
zc=&g{{LE92cyKrP_}m%QJg~eJAeH@u2h6+gtJXf{!J|*v0dMMg@KU{KAn-X4p6~5*
zoBWaokMj-K{8u<X-nZEJEe{$!kCZ%p$Ai%#L+;*I9&A=i^%!YGyQs6jsULa3`SK!o
zaTgD6xSV);|1+LT$)*Xey*y}L1S&n>crZJ0>GZSxJgDja^FjG155~KsP;!TGpLws<
z)(!Ju<)r-BS0g;Ia+hd0B*Diz%^H%D;zN8#$}nR*AKt2DzR{V;2gwZ&84o7&VQRor
zJy#h%jElZ`?Xw&o1{3W^BINn-HSV0wjG27+6|&aksuCaOZ{4r4W;P!J)!%nMor`*J
zqhr0*_;4@Edg9=GKGY~n?T%f*2O|cbuCj;^Qj!y{-PGm-xG2-sF6G0IOCrw~%lM#o
z`<bNQ3O+~_TSX4*^Wk95(A)Txe26?#CqHi$9~@TYYvixyLsLVf+&VKpEO%V=qS2fW
zyK<tA?_0}<qumkxLsonUksh^+wdF&wM0@BQdp?|asTiH@z=x#RlLt&4`B1JrxVUx;
z9~PcDROacz2jPPXqpz-fNO;>G72(c@n>9Jt<vjSXBt<%#yOR$S^`4*8-OUHr%frSc
zK75$h=~cdQAMSU|m&Nb=_~6nR?&p7y4=?#+&ipxq`C8(1B<=_wK1O<)%s$G83s<Gz
zUOk5SN?pFy=mZ}ux0PNkKgox~A2<z8;h2BJ>n~bQ^C5F{ZC=1xK2#Zg@fkVChd(V{
zUt?qVuxU&77UeiTJc&_`&xq&4>9Y!4y-T=1i^!M~3Lgqf?d&$B@WDxHQ)?5Q56<I{
zJMKy4Lq)|E+IKb|ZdpIP6~^U*R#aZ*Bmp0E!?y;eq+=fS=?e2R`LH<MEiU^S9}H?E
zUKwTcVNBllkq5baVAxZK9P;=unPXSgbc+wW(|-kc-{r&gT`^yJ3b79AUYmswd{{Hk
zxL2~64{24NJ{RxvLCw+KNbw;b9!KZA=a%y!=b)*Fb|oJgB);CrujYfbW#-G(HGKRY
zVE?@ADIaXAN-`Xt@!|8-Zu=L{vF?87R(QVT!=1$bnV(+qq3;QMRlr+5)Za{8G5C%T
z?k0I7;jMgN`hMg{xATE*GDAM$BOd~8sJYMV;=|mEo?!N8TsO>nk47&aRy6L>%KFBK
zyvMJKm-q8w^{#ZSf**W%m6Wx|bO`skvp)FYZ$6x@9C5Q7;X}sb8FEj@2w<}BFx^>F
z0L5%t|EqBVsPTO~f9C`NSUe6_Z=WQ9#Sc?{_)Zl-;-(*|J+cBg6%Z*OI9&h|>gqdw
zDhNPhukd)Nq5%3$boPuY3*aPo_kzed0;r&v-5IYcfL4VoN*Cq{;KlrK`>7fNxJ_Jl
zOk5~{$QFe~@>&Aen8>P4UMzrxj~3Y|=?FmXJm(TqR{(QsTeIfq31BPhRSL&I0FQe7
zHmey4AiV5&GjEju&Z>5<QePv0w%^>nm(2u_<rd|=z(N2WD?^uOtQ7z|d2h`^YXP(z
zUS@K|Rse<7caLe>3qX(G9(lz<0B-$kPfbSw7~e0ObY+VGe{Z#jTIeDG$BLgd>8=7u
z8>#-N;f`@%7t{)N2w+Uz?sM~Z3Sg1nx5+#&0c=0Ee5b090MdU>h-B}@xVyFo&Gr*O
z=7#C|i~|A)nI}^^^N;`%ybV;ShXpV&<Dkv7AOYOfXF_~1=B1c7esYKa!tU&2Mu!T3
zrP%OWG7R(H`baw>LI4pnmS_&25x}k+tzVBt383d@U36cx01T9yyZtW+z`XtQoUS;`
z(~J+Zd=fCuE9vd;F5$W?<x|_J0+9Kx-qMgF0N$0W(>5~DFMgG(OclVOj_8F2M*vSp
z<PH*C0ld@>t}ze@Ksm4WcTPI)n=*DlGgAP;HGPHLYnazv=k1iU1z=;Jo^dHx0InC8
z)lAI8{oXULm^H2lGUR?4uaYW)39=bK+$4+O_e8ZnA!Cccdb7m_rbH3kHTu{g98(0c
zGM$yLM+umrKL5w?2mzNK9oVmeblfP=f<HJufvf8^OhBgPnuxQ%38>rbtaRlU0b9<M
z4pt5k@Jq*Y-j_kNV~<OcL%yD$Xs7>^0IqD4o$C(*1h<4qCkF_)<0U(v*H1vQ-_60Y
zJ_4rC4V~Hjoq(Pyt57*)@57#jhTjNiao?c4{VM^n?^ka=(@VfD0j(~xhXC~-D`~Y~
z2w*K5<p21LcCmT=bCItC%;VN|6QI2)-1lG?0awOa#U*zVpx_wLU-XH9HQq0I9Ulqk
zjB};RBHz!ht~Tl*pmo<I8?SZ(?&p4;5Z6ZF^T1W+&Iep4W47!2RsxK^Ds)Uju8UGB
zFlZrQ;+^R=JKq!F`m9<r?j3=@XKIw)ZN`1h{E`0QEdh&~XIIH0-JTy>zWNRB$Nur_
zeXj{9OrBhs{EC2xF%yPLn+T9<PTl|YCGJ0*r>25Dtly_-|AK&z4*MNXG!kH@+V(i}
zIRVB$R$h3~fceufj~S0#`@5;q;28mJRd3e&)DvLUrY@UYN5BetjC|Qs0*<cy<T_AG
z;CHN9?-n5CLx<8`Y6uuD*_j>vgn+FE+5-iTF>i}>PjyujaQ3tPx;e;}r`t9-R1r|K
zP(L}Ml7MMu`E&9g5fJQ`_o2OlfRdrfT}sG;HOyu9<piijue=prMnKiVd&lx05^((B
zgcBbhU_P?G3TGq#tZ-blsg!`7Tcifh-Y1}IM_XTE2>~90>ZQHK1cV1s(>0KD<}15z
zyGOvlu%1I=5dlrk#bM<H0T)gse;yVQ(C{#kriTnU85`~=B*0?8^D?iHfZ^~@<C+Qx
z@H;pvkV87lgooMO#r)=*Cr8}Dd}q)5RB)Ss(hWnw-)<4Gr;Tr}jr1=awBDUh0Ndz!
z3@eWS5pBA3<4po6y0OKw$j&Xw`L;I*IJcSB5}Av6J@wA+UJllu#^avBYy#>TE1&2g
zU6&5G`(MX;(Cs(6mPNpg??Ywn*9e$W{xNPI^6M|!#T{1(nAO@=!N?@Q{c)yA;}rsq
zOuF}JI&$%sk(!Me1lX6Y>xoMzAnn5etE$U*{&a^I(#T$eGVVG70VgMC@gw;JxcI$#
zbDu|mjbQ1<F-QgJ@F7z!0nr_quTP~B;FT0LKyU~Ma__YoW@8?bN@`7z8_z#W2xZ}Y
zxHvXiluE#cKF!;|nFPpAQkk_Hd2@d8?UM{VuZ!)`MRWr0vi2vA(6HXp7T1~~6F-hx
zo=(B~FDvaSO~!hOaj2C-&W$?NZbc=adV==y7zzQqe01`vlL+8<d+(D$I)A>h#qkmW
z<^^?u)I<VwxWuCu3An$0i_No<I?`rxp7Gd6G-JmKE@I!fBBAv$j)0XfPDU<9#!J7L
z77&Z~D5mS#tqTNrXJ4uuioyLo%=m7K+@LD8@$7j5gfeHFE2A+FAI9)xkwHC9DbD9`
z-x1sHrA86(qO)76^(^*>&=2X0kgCzMeGWwuaA>=~%bha>{%(IK<j-lmZ<$YO)*_ut
z;+^9naQ*gYQyapu-uE5*IU7l1Zkg^KM!>gf3y<rkFfSWaU-zHHz7p`@{2HX_?8Lq2
zLJ4SmU469n1Oc1M+ioc#e^7Nayh3n)^vDNU$Fc9syO%U@4D0E1!6h@~qX`>IVuG<=
zq+L`Rjv{Nbvs93AH+A;;24P+P$=h-(5c`wV4FAz1Sf7f^?%E-%qF?Kz9L7AR_cycz
zU>_@Sy1Eo8SlpNua)<y+!!2*h{PDb&XRMlzT<U+d+~Xht6V|MZ$U1=ctVe6_;C=!Q
zG`NSZLq?lkA(H$E_;RX4`@Js#0U!FxmLe;QuU-t@N5GTat*0yZ;(pbR^A(YA*RAXI
z*+amI1seNq`JgP@J4F)t*g4{@qc;JQZStD9yRj}-dP;ot!hRjJ%gGEGwuw?7zYF`*
z?Ap-gop>)adpGJJ+aEk}2=gSs&|=;38V`Kl^eua=g0$1zXLo1^0jbL#eJ|clK<#$z
zGC5@4=)U{A+zA-(>+vyf8v*prHbW`olKRG6XEy@Q2g+=|;!40Vl|6HZwqhSpFPOF-
zsnyV|!*;>ABW?xtI%8i6k^f+h?2~N_p*j(8mD*_ZaSPU)Yp?1m<kzfUIu|z+F!g=C
z&wED#V%BYbsE5>Ey2UPf69KnB=Kg%O5$($?9_t{lDs@zxaUdW=H>j&|1Iiv1dW(_a
zeT9N>d(5xd!c`6H@t#pTek?-fT0Cz$Wk)~~^KomvEyn*PW7;C5b*$;0Q#N>>We4Ba
zTN5DnMJ8eq(qZ%P+Au5Jhn3v2hIIt^DlahAMn0W0<4D9>tOuU3s?idk6F+;cmm*&s
zO6ZKVAfU)SF!z-?0oNu6GnONRo*uq(9`M|LLDM@k0<vBg>8wOvc~E*K-jo1_r1sj5
zH3Y2x{6x|e8Bue-BY8Ce62J97e=)&)mtI$}LXHov^5m?-KDTGzlOM)7{$b~ijYw;=
zc-c%N%u8j}hcQU^hObR-EAiguzWSDDi1#bYYw;9hhw|=oJ_ZD+tiLn1SRdab5)8RY
zNL%B39!K;DsGw<@K3;+SXV~6!Au{If1kZ@&1USty;J(tu^F4Z9!T_14Y#W)d4C`ZS
zyGEA|-q#5Cj<v`qZBZ3zOYyn!-1O0}CD`xyJ3l!iyVzx#Hy7hQRJ1ua1zEqdRC1p-
z_Qfh0?gK6C`<JENRgq&LT0e9V_T8SnOI~OaP(B)Ip@&>wb~`v>Apy*U7mvCZ5D=RC
zz{mzEJnB?_S%ZL?dfta6klT%(7<;H=yp6u9qWJ_o*r70A2^pw*#r)_zd|x;)e5ziJ
zfa}zr#%0K|&aL(rRS7U9rZ;z~5cs>zP^=BoXHJK6##{m}ug^D;LVB2(tnr$I&!N?>
z+e&8RJxgm#Qb8Ks7IvJRh38yz%lWl3-nZ4+?~IX0H-??1EB!}Zk6rsy5zj4CUd08u
zC?Rh8?U~r0Jgny_Ajj`ju?w7m?{!1narFuWnC!XGxg2?Z<eBFsd3=wUF6jL_9s5Rq
z9eoorrEit{jcFL~t^z|jq;9;%;s7}U465CXp335TgzpLuU1U<hc~+tfz9&pt+4prS
z_T5}XFGnPOP^$Ch6nuUeHbqZI-YWDr4V+BC=4WYB8Yba;QiIiZeI#XyX>amG%-<7F
z$)6K&AN$VfxFYqJ#sw5g<NLsY^_8=bamP7EC&v@O;d<VCJB|QpL*H#?$g?rYQ+QGY
zeCWQ{Fc$eayExBVl7Jm*M>5LCV!xu^$X$fI;pSX(K?3Wsbk*1|V+fe7VzPcCa?d^8
zt2akQ(9kz9TLGym>mvvr5rOh9SF@%+B8Xgdt=9yp`C0!KXIKP}%j@GLkgH~Ugn0iJ
z!OB%PLn?lWAU-JLqBe5OY^B>5heY5J_P+P)pa}N#Hmz|&cHC>`-2Ewn%WEp9&qnSh
z0vTaHM9>>su%dN91k?T;ZeE9s+j)X@wO<5_rLP~Dip(=_aX8#3g6^ID7SF$nfP1^x
z#u!OSoV%0#O#}%RpUz7lUFQ5M_4z7-0{^8lRlOoOd11Ss4sut?t`|u?80RiU=bv9h
zAS&AW!5vw)@XzUz&mx%mrA2oE(refI-k5F?#5UC5>**4KJ&nTNf(-aI99z&S!k@Q0
zW9B0Nw5_0>`6PlLQ*Y&e{3rsyS=R6Dk&E5SmGe48FsL?Vw<3~3d>5W<7eRBF)4bL;
z5u7QQoMMf16@FLB{vd++&p%w9j^w))+8l2cK@NSxz}pt|-+Cg~66qZ?``Fd@IR9p%
zy)5#~QZ=2U??mwAvC-Vu%_6XvtfOp>ylwny-j%nQ-%Q1oQ<1EzC9Z*QM351ed#>rV
z2#WakO3aW2s~?TecqIbotPIyF$Q@?-caAinzU_L|rk5h9pH)sVL)J-#%V)g6d^&I9
zOhGQ$sh}U&h;?x^y`ky32zuJyh6B=f4Bg;Lg9x@L<qS+k9zOM}IOrLkOT<>{>v|Cc
zA4xl6f$S1o-*&YQ&(X)tUJgm05Mh1nDdzj_ul3EfB52gV?6wYh&(kH~dW{I)=X58?
zBP9b(3r;){;rsKV&X&g_IKqyaZ;P}`*YUnlErPW^BiV|``8u_uVO1iywJFrCqf!KS
z=e`wgK>lg3QNR5N^?qqn=ODLbRZl)!A%Y2Zu@}0_MIfWOR&EQ@i1L&rEE9plyGgq9
zkqX<NRK`3Mf!nt4zTX~*;ElGEoEy?t<q1(ziggpU^3Wor@a93y#QRvceYWjCOGI#y
znUU#<{4u!Wa7D2QlsxUsmLc=!T1}_j!+l3B?iwi)!MTqO<$I8})H#_oga|U<yt!zI
z+&As(F^&l9BT{$Ic%)wBaaVt#2&#AS9U2Qo(0=H=jVTiU1B6vZ0iOGg?^ZI%;b!gi
z!FNT_T7Q0X^BobaxX`iP3VF?0#V_YJp4+)eCuSn;4n)O;-4elXS6^O7z6gd`i3JWw
z>6qtrck)E=^WgL^D#%FIyh+hFMR3CB?V{crB1l(Ov)zh(IPsxhaW2MtPH<5ZIlYB`
zGa(22fp1vj&ukH}59<tjBJFG==RdlRviVk9UF6QJvjOxhtfv?4v@uAVna;)gu8Cl@
z;Cjc?t0FL$tC(txY&a06&&w3y_m%8z6Oq$HmPQ=Dg7>hh>2gyB=5=@6BXeX^ifZ@O
zbj*ML!70;`yC*v?IdNG8sXYR_R)Gk#b6~F>vVkdaCXbK(XWV_}EaYWVx0`2pSpT%k
zWnElccju`Wn~{25MIFL4tm`P~n~z*|_TBIWjtE)|j!O2iMR2wI^>}xr+mA2fAF%Kq
zWkg9XMkZ}IHJp@+^Cr39elbN*et)3V3wcv_@zW}X2=tq@3RfVfPkGE{(y^{)eLgi7
zInK%3&5wrnbCRrXeG29^zFcAzQhm8k4KEq{s6|=wL}Zox{jG<oBH#?HP-vov;F-}t
z6(DtWT0=9FL=d*HVWAu{;z)Y+u}dNtyIa@qU7`poZ8nazMy3V5q2wfBovqAMnTbqj
zxyC;gF9MT?T-ElAB9PPdPqs%k9sMGCD-QcX$DrSAq*uLP)tOik+_*4ze%A%m+u42C
z5osgHDJ+N)!Mt?yAywqgJ7K!f=dmAYf@e=O=Fx6vlrz%VzmQLy!@6u+a917qZ6v((
zLKN296#m0+XR#l8?kRFb?y<RYy*N??TNTZc79vm0o#21*4ECd}H)j2(MG!ijChd-F
z=T9rUAAxzuUU+;FG9dDXT0*!8RF1~n8weA@(Py!a+mWk!FEy2(65)5MU>7Z<4uATi
zgp=5RUYYCs2*vmmn$K@X)^YB?C_N#9nCkLbT1buV3de*Hyx$IKCkKw>eX-I^bw|pn
zugboE44-F)d3l=1f?D6Ki@_+Ho6!1>Vt!wZ33fwze4$zw2VwrJ>LxEhI^VxkdLd8*
zb2M`Ie?5Y6T1rW~AQv6lPZb@;ybRuxpNHJJyy@(@0L-T^^<(!T5sd#<W#WjeSKfO1
zj=u;t>B|?)My|Wm^e*BcKF>~{9c({<b?@9Gu^#!z<J<4t{UX@^Wp9T9^1$XL_e1<d
zKy6lvYxc$ZYL&LOL>|d~_U+0(5!^|3iJO94)8RMs@Lmx-{nnN6e2)kWC*B-1MjD)=
zZQ}S~pW+;+jYX>6(|NMT8}IQF>bI)hc+RULrIsPXjT^^NyhQl-#omD*yHH*~?b$Zu
zjxjl$dpoh7J}J7+M;171jYfH5zoZeAP7kc}(X1I8kP*#wXK(Dl=hE2;E%Hb!OKt6A
z+p+KaP22g(9iJ!m^fPOazDA5o+-)NG7~T^niA)MS612w+pX1)6)|IXzVDb-1FF~FS
zXXnLl#k{QFZTrmy`@3su^%msv0afEW&LS8ItiGg#G+J={`3WcN18!ra-fj`WBrnss
zW=O+?Df4)nMbNo?ovb9Xms$GJ+ff8lW_w;M--LRS11^h@xdM+bF&l9o8-fpZJ7C}6
zJoAM;(p5rPJ7)vNd);-9EOOlC+}Hqn>{mvEjJowG6H_j&K=z-i2%y;EK5EAr_uGnq
zC&B1+LUwz5pSf*=`L{EZoQbp>D|H~)8vD`PIHJ)C=a=5<H$+NXIw{cBVO;*!vwyC|
z{EhCNxE0xMlh}O66785w=1k<BwY$~?TcADtUiou#%*)3(Ljz>%>5Mok;Pbqo{ZXG8
zz8AgL=-!O9U2WQT(-iyI;wde1$o0SW-8!@e;~msFShHFL{>P0bFGWUFddA0^;C)^l
z{-bl12%^;wTiPIHYGjXQ7-N69Xb>xfl*l8Zyo~Vq8yD_fvQh+6hg>yPkS&?>YC;Y1
zIS_l>waEa_XKzorA#%mV3?-^Q-VcX@4d3*zZ~yw^<A9WGbKI4+0_$vByx9cgNY2GE
zKFjerzC!!*eO<g~2Bz97$U74*C!JU(!vCjh@U~G0>uq<<L_OsAi`GgBOYxpmhEC{S
zg6|oA>t0zQf0l$@<S)kk%^5p?L|X)jwO46wNV?Xho?BWXs12=MAcq{CIo5RlBFy9Z
zQp0jh5$xUHB|8tfq}-`8bRpJR_G_=k1tKU7aC@VRoadEk5UYXpQ^)sjR~Nw-IrV5W
zWalBtaK?OmZknmM^v%QfnWr13Y(ORi9pGiCVf}|X%@~7hHh90+O%?maKqoy<1@GHL
z`^?G6I}hI{d(Xvt7<^_|(HuPIfXS0*BBRH;B>B(A_l&UBKg(xf-`8Upt0LWltlW+%
z<MS|4#idpW-)owFFVjS>Nh|0GQxw56y=5mFXJTJ$H||=BG@Ouc7&Su#Z@Yp#-zeby
zXm<0~LvGezz3zfMKF?PFkZ75X>)CiQS0b62+hpRW;rmB+;f^*rj4RXa{3_(EE`|6+
zS<Gu`nSX~2)=8w5_G;v-Q-@0~O%=g0iFwOEPQmwk&t)f9BilX&^Da%sy6LIoc1*(l
zRx;~^3DV}q;-!fbMGz?SuCQ$aK35mZsu?4t-BLU+N@IO}bU5BJUIfP%%KI50_4gI&
z$Be^z8`gaJMhf4%4VF6TBA<I*$&8f5_tNJ9&l|_$^Pze7D=lRIlWliTO5pp}z_J6i
zW3UgW7Y)usrs_0r4;mH1=a&8S@)04(%b8~>A)6)|CHwsm!g#gCF2t}9b_)8xOhx9m
z`EU39Erk2}LA)Ekgy7QFOh_ON*T&>*9ung3ej9=XgF+bEI$ojwrw|TgUyNRhwEa=^
zl==hpBM<ks4G3X-`k4><NX_i(>u39g5YoMHeM6rR6mR<9Q%BZ1oEjhaT?jrima3P2
z69V7aRbd)(!_KBA&#yvQdhSACPOlIqDp|D*_Xy$Io{stJkuOWv88N=#I;TwMef%th
z?;Y3P>LdL(`5ZjcErjGrHcxB2(0<R-@pF;G`MolJokDoB*0<}<Cn4<d4WdgTV~=zz
zZTcvLNjEe1vpR&3ar9i`r*<J|P4kS_N4hpl*%Z+xgx|7$UDY3i@LIvqMG;x7$D(+(
z3gKYU0)BRj5OT*pIz9Ma2zE|8R#+mXB}?)zz7s;WXVtiu%|ehAxh$HGZ16Fgec-JS
zT5fnW-+m*6AD`n7kG>YdE5B<EHb|$?Sm~r!nBPNFCcSRLIB3#u)REKHZasD2r4Vw;
zpMSaaLI@Ryb}aqVC<KkfM<(l#F~`3tU3`w`6_WR`p+N{qZSAIWkPj#BiQfH82x_aQ
za<0{5yydsge5=DejT##mA-@>}-41&yghHXIWLd2c;+q3iCL@i%e;MysBZRY7(L~A<
zA?PfbW%lN=5GM3T#Lq)E^dw)~Q;p|)P>+&TB?RvhJ)5tULMSjgT&a&#(JY@C@<<5V
z-b~Rgsu04etEw|3km1cQ%hr|)q5tu8i}PhzZ)I0vsvipBg3=L|405%i{Rzhhxc+oG
zb!w>))^DSyG~P#hC~i<fK2h*|yRAeB)r&UYU=`!JJU+6a`JNCKY)-wUid2}G(Y~_?
z{aGA*&m)AeX=C-}R*?{NBbY1HkzW_wj@d25y16uPRZu7dZFxh=hXSnY;Z5t*k)~|5
z60f^Lm|`6+$-g6neRkHeEw`~=2RDCEMYeEHhI`x+LSi!K11nz$pB7SNU*`!SFI`_!
z8CmGH<DSb+ArLtqfqDb$>s9uJx?CZ28D2=0L8d8e4P2iigy4rg)6QoLVFhjRv4_`%
zz<&FYKL#24=?`^vmJpZ|Y&IOfCInGDw>tl-5HhP<WP3BQ9^^WdwUOhqqdUE?2w_|w
zH!Lkf2wzrAX=+NxIw%qh$|G4H%-?RfjQbrZITtO!`?T(Ne+gd*C+1gc{^DVtn<uI2
zA$L?#U+w1#;j+Hvw##W)-vM9QZ#Y6Is|ik5KrT+Z5n#_2LgR|#qme8j#J#%WC`uK=
zjE~BDzcYnkbT!jV3t4nn=d~vT>*`XJaS9#t=XQSEV;a_l@-)-2NHg_cA6BLa!J}fX
z{lR3+i;GJzpNe%lFU;!&MF<8VhvcRpO>?9oAV~;^rzzbGx`g}aixgxg3L!Fm{?<1M
zLg+kJQ!R^36YY|)jK{ul|IW{&7cn2NuH<CK3BkricH!$-><8y6c1%S+A+9-_T@V7-
zhBEbV4CXu3J&J!G?`is5;#oAFn__>y6jDb$#LwUy=1<A`y-$=7&f6HxNI8pjDyccK
zJQCwf{URDTBgF3odzLOh#_;#=a5^o7Bl$1ZM@9%C>%-Z>o8dzE6S-mc`!KAx_Q`QF
z$jPqEBdbpdVe8UoGyP8DJ-?lJoE9nscZ+P=!xKUfs7nTY3lYMxYX_#PAo-dPciSDu
zeRh8c4nBr`+9tz_AB=rFuD$-rQLNL+Cua`@2|*~oR!swG{!_8VArRvmU%w;dh!6sk
z>{A4XvF~*6y6_|b``Yf$#y=0?er?Myt07nSR#w^i<GSf}Sw{}y{VKX^oq9kB1Ls^g
z5B3Y;(3bPL-F`wi8K@a5kL*|}BVpo;=kaZ(xz|37>+5Wjn7w%Z(LZ}~_uze(I=j8m
zM+h8YYUHpt_KzVY@A*iF>5YG^cMIW)bM~f#UU=>bA={I73Be>ye%{@k*w3xEa$b94
zJ#ANM81=w<(6ua8M?RJ3?z7$@1of@56@J_Cd9r^0+l%g~Zyuj@V;i1R%;x3K+=LJo
Zl;l6)DujDm7VTC<798u5G2SYK{|~G4o$vqv

literal 0
HcmV?d00001

diff --git a/tests/tests/test_integration/test_gui/test_regions.py b/tests/tests/test_integration/test_gui/test_region_ui_functionality.py
similarity index 58%
rename from tests/tests/test_integration/test_gui/test_regions.py
rename to tests/tests/test_integration/test_gui/test_region_ui_functionality.py
index 18e1334..abf05b6 100644
--- a/tests/tests/test_integration/test_gui/test_regions.py
+++ b/tests/tests/test_integration/test_gui/test_region_ui_functionality.py
@@ -1,8 +1,11 @@
 from filecmp import cmp
 from pathlib import Path
+from time import sleep
 
+import numpy as np
 import pandas as pd
 import pytest
+from tifffile import imread
 
 brainreg_dir = Path.cwd() / "tests" / "data" / "brainreg_output"
 validate_regions_dir = (
@@ -11,8 +14,8 @@
 
 
 @pytest.fixture
-def test_regions_dir(tmpdir):
-    tmp_input_dir = tmpdir / "brainreg_output"
+def test_regions_dir(tmp_path):
+    tmp_input_dir = tmp_path / "brainreg_output"
     test_regions_dir = (
         tmp_input_dir / "manual_segmentation" / "standard_space" / "regions"
     )
@@ -59,23 +62,57 @@ def test_add_existing_region(
     assert len(segmentation_widget_with_data_atlas_space.label_layers) == 2
 
 
-def test_region_analysis(
+def test_region_analysis_without_save(
     segmentation_widget_with_data_atlas_space, test_regions_dir
 ):
+    segmentation_widget_with_data_atlas_space.region_seg.calculate_volumes_checkbox.setChecked(
+        True
+    )
+    segmentation_widget_with_data_atlas_space.region_seg.summarise_volumes_checkbox.setChecked(
+        True
+    )
     segmentation_widget_with_data_atlas_space.region_seg.run_region_analysis(
         override=True
     )
-    region_csv_validate = pd.read_csv(validate_regions_dir / "test_region.csv")
-    region_csv_test = pd.read_csv(test_regions_dir / "test_region.csv")
-    pd.testing.assert_frame_equal(region_csv_test, region_csv_validate)
 
-    summary_csv_validate = pd.read_csv(validate_regions_dir / "summary.csv")
-    summary_csv_test = pd.read_csv(test_regions_dir / "summary.csv")
-    pd.testing.assert_frame_equal(summary_csv_test, summary_csv_validate)
+    # ensure data is saved before it is loaded again
+    sleep(8)
+    check_analysis(test_regions_dir, validate_regions_dir)
+
+    # check saving didn't happen (default)
+    test_saved_region = Path(test_regions_dir / "test_region.tiff")
+    assert test_saved_region.exists() is False
+
+
+def test_region_analysis_with_save(
+    segmentation_widget_with_data_atlas_space, test_regions_dir
+):
+    segmentation_widget_with_data_atlas_space.region_seg.calculate_volumes_checkbox.setChecked(
+        True
+    )
+    segmentation_widget_with_data_atlas_space.region_seg.summarise_volumes_checkbox.setChecked(
+        True
+    )
+    segmentation_widget_with_data_atlas_space.region_seg.save_checkbox.setChecked(
+        True
+    )
+    segmentation_widget_with_data_atlas_space.region_seg.run_region_analysis(
+        override=True
+    )
+
+    # ensure data is saved before it is loaded again
+    sleep(8)
+    check_analysis(test_regions_dir, validate_regions_dir)
+    check_saving(test_regions_dir, validate_regions_dir)
 
 
-def test_region_save(segmentation_widget_with_data_atlas_space):
+def test_region_save(
+    segmentation_widget_with_data_atlas_space, test_regions_dir
+):
     segmentation_widget_with_data_atlas_space.save(override=True)
+    # ensure data is saved before it is loaded again
+    sleep(8)
+    check_saving(test_regions_dir, validate_regions_dir)
 
 
 def test_region_export(
@@ -84,7 +121,26 @@ def test_region_export(
     segmentation_widget_with_data_atlas_space.export_to_brainrender(
         override=True
     )
+
+    # ensure data is saved before it is loaded again
+    sleep(8)
     cmp(
         validate_regions_dir / "test_region.obj",
         test_regions_dir / "test_region.obj",
     )
+
+
+def check_analysis(test_regions_dir, validate_regions_dir):
+    region_csv_validate = pd.read_csv(validate_regions_dir / "test_region.csv")
+    region_csv_test = pd.read_csv(test_regions_dir / "test_region.csv")
+    pd.testing.assert_frame_equal(region_csv_test, region_csv_validate)
+
+    summary_csv_validate = pd.read_csv(validate_regions_dir / "summary.csv")
+    summary_csv_test = pd.read_csv(test_regions_dir / "summary.csv")
+    pd.testing.assert_frame_equal(summary_csv_test, summary_csv_validate)
+
+
+def check_saving(test_regions_dir, validate_regions_dir):
+    image_validate = imread(validate_regions_dir / "test_region.tiff")
+    image_test = imread(test_regions_dir / "test_region.tiff")
+    np.testing.assert_array_equal(image_test, image_validate)
diff --git a/tests/tests/test_integration/test_gui/test_tracks.py b/tests/tests/test_integration/test_gui/test_track_ui_functionality.py
similarity index 69%
rename from tests/tests/test_integration/test_gui/test_tracks.py
rename to tests/tests/test_integration/test_gui/test_track_ui_functionality.py
index b9f25e1..0534f6c 100644
--- a/tests/tests/test_integration/test_gui/test_tracks.py
+++ b/tests/tests/test_integration/test_gui/test_track_ui_functionality.py
@@ -1,4 +1,5 @@
 from pathlib import Path
+from time import sleep
 
 import numpy as np
 import pandas as pd
@@ -11,8 +12,8 @@
 
 
 @pytest.fixture
-def test_tracks_dir(tmpdir):
-    tmp_input_dir = tmpdir / "brainreg_output"
+def test_tracks_dir(tmp_path):
+    tmp_input_dir = tmp_path / "brainreg_output"
     test_tracks_dir = (
         tmp_input_dir / "manual_segmentation" / "standard_space" / "tracks"
     )
@@ -72,32 +73,67 @@ def test_add_surface_point(
     )
 
 
-def test_track_analysis(
+def test_track_analysis_without_save(
     segmentation_widget_with_data_atlas_space, test_tracks_dir
 ):
     segmentation_widget_with_data_atlas_space.track_seg.run_track_analysis(
         override=True
     )
-    regions_validate = pd.read_csv(validate_tracks_dir / "test_track.csv")
-    regions_test = pd.read_csv(test_tracks_dir / "test_track.csv")
-    pd.testing.assert_frame_equal(regions_validate, regions_test)
+    # check saving didn't happen (default)
+    test_saved_track = Path(test_tracks_dir / "test_track.points")
+    assert test_saved_track.exists() is False
+
+    check_analysis(test_tracks_dir, validate_tracks_dir)
+
+
+def test_track_analysis_with_save(
+    segmentation_widget_with_data_atlas_space, test_tracks_dir, rtol=1e-10
+):
+    segmentation_widget_with_data_atlas_space.track_seg.save_checkbox.setChecked(
+        True
+    )
+    segmentation_widget_with_data_atlas_space.track_seg.run_track_analysis(
+        override=True
+    )
+
+    check_analysis(test_tracks_dir, validate_tracks_dir)
+    check_saving(test_tracks_dir, validate_tracks_dir, rtol)
 
 
 def test_track_save(
     segmentation_widget_with_data_atlas_space, test_tracks_dir, rtol=1e-10
 ):
     segmentation_widget_with_data_atlas_space.save(override=True)
-    points_validate = pd.read_hdf(validate_tracks_dir / "test_track.points")
-    points_test = pd.read_hdf(test_tracks_dir / "test_track.points")
-    np.testing.assert_allclose(points_validate, points_test, rtol=rtol)
+
+    # ensure data is saved before it is loaded again
+    sleep(8)
+    check_saving(test_tracks_dir, validate_tracks_dir, rtol)
 
 
 def test_track_export(
     segmentation_widget_with_data_atlas_space, test_tracks_dir
 ):
+    segmentation_widget_with_data_atlas_space.track_seg.run_track_analysis(
+        override=True
+    )
     segmentation_widget_with_data_atlas_space.export_to_brainrender(
         override=True
     )
-    spline_validate = pd.read_hdf(validate_tracks_dir / "test_track.h5")
-    spline_test = pd.read_hdf(test_tracks_dir / "test_track.h5")
-    pd.testing.assert_frame_equal(spline_test, spline_validate)
+
+    # ensure data is saved before it is loaded again
+    sleep(8)
+    spline_validate = np.load(validate_tracks_dir / "test_track.npy")
+    spline_test = np.load(test_tracks_dir / "test_track.npy")
+    np.testing.assert_equal(spline_validate, spline_test)
+
+
+def check_analysis(test_tracks_dir, validate_tracks_dir):
+    regions_validate = pd.read_csv(validate_tracks_dir / "test_track.csv")
+    regions_test = pd.read_csv(test_tracks_dir / "test_track.csv")
+    pd.testing.assert_frame_equal(regions_validate, regions_test)
+
+
+def check_saving(test_tracks_dir, validate_tracks_dir, rtol):
+    points_validate = pd.read_hdf(validate_tracks_dir / "test_track.points")
+    points_test = pd.read_hdf(test_tracks_dir / "test_track.points")
+    np.testing.assert_allclose(points_validate, points_test, rtol=rtol)